LXXTLJ连续系统频域分析.doc
上传人:sy****28 上传时间:2024-09-10 格式:DOC 页数:2 大小:43KB 金币:16 举报 版权申诉
预览加载中,请您耐心等待几秒...

LXXTLJ连续系统频域分析.doc

LXXTLJ连续系统频域分析.doc

预览

在线预览结束,喜欢就下载吧,查找使用更方便

16 金币

下载此文档

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

连续系统频域分析模块该模块实现功能为对输入的连续系统函数进行频域响应分析,如图所示:输入分子、分母:num=input('num[n1,n2,..]=');den=input('den[d1,d2,..]=');handles.num=num;handles.den=den;guidata(hObject,handles);绘制零极点:A=handles.den;B=handles.num;axes(handles.axes1);[z,p]=tf2zp(B,A)zplane(z,p)guidata(hObject,handles);频响特性b=handles.num;a=handles.den;sys=tf(b,a);[p,z]=pzmap(sys);p=p';z=z';f1=0;f2=10;k=0.1;p=p';z=z';f=f1:k:f2;w=f*(2*pi);y=i*w;n=length(p);m=length(z);ifn==0%nozeropointyq=ones(m,1)*y;vq=yq-z*ones(1,length(w));bj=abs(vq);ai=1;elseifm==0%noppointyp=ones(n,1)*y;vp=yp-p*ones(1,length(w));ai=abs(vp);bj=1;elseyp=ones(n,1)*y;yq=ones(m,1)*y;vp=yp-p*ones(1,length(w));vq=yq-z*ones(1,length(w));ai=abs(vp);bj=abs(vq);endHw=prod(bj,1)./prod(ai,1);axes(handles.axes2);plot(f,Hw);title('fpline')xlabel('w')ylabel('H(jw)')guidata(hObject,handles);