资源描述：
《【精品】系统的复频域分析.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实验六系统的复频域分析§6.1MATLAB函数lsim(用于系统函数)目的用lsim仿真由系统函数表征的因果LTI系统的输出-基本题1.定义系数向量al和bl用以描述由下面系统函数表征的因果LTI系统:2.C-2定义系数向量a2和b2用以描述由下而系统函数表征的因果LTI系统:3.3H.(s)=-5+0.3定义系数向量a3和b3用以描述由下面系统函数表征的因果LTI系统:2$比(巧=二一35+0.84•稠lsim和前面部分定义的向量求这些因果LTI系统对由匸[0:0」:0.5],x=cos(t)给出的输入的输出。以上四题解:t=[0:0x=cos(t);bl=[l-2];al=[l2];b2=
2、3;a2=[l0.3]b3=2a3=[l0.8]yl=lsim(bl,al,x,t);subplot(2,2,l)plot(t,yl);y2=lsim(b2,a2,x,t);subplot(2,2,2)plot(t,y2);y3=lsim(b3,a3,x,t);subplot(2,2,3)plot(t,y3);>Figure1-nxFileEditViewInsertToolsDesktopWindowHelp3easi墜1a
3、Q头◎®題/“Ia
4、日io口nE9§6.2作连续时间的零极点图目的这一节要学习如何在一个零极点图上展现有理系统函数的零极点。基本题1.下列每个系统函数都对应于稳定的LT
5、I系统。用roots求每个系统函数的零极点，如上所示的利用plot画出零极点图并作适当标注。(i)H(s)=/+5+2s+3/Til2s?+5s+1200H(s)=s+2s+10(iii)H(s)=2&+5s+12(52+25+10)(5+2)若干不同的信号能有相同的拉普拉斯变换有理表达式，但有不同的收敛域。例如，具有单位冲激响bv.为hcit)=(r)=-e^a,u(-t)的因果和反因果LTI系统就有相同的分子和分母多项式的有理系统函数为H<($)=—-—Re(s)>-as+aH”(s)=—Re(s)<-as+a然而，它们有不同的系统函数，因为它们有不同的收敛域。解：代码如下bl=[l5]
6、;al=[l23];zsl=roots(bl)psl=roots(al)subplot(3,2,l)plot(real(zsl),imag(zsl),'o');gridaxis([-1010-101()])subplot(3,2,2)plot(real(psl),imag(pslgridaxis([-1010-1010])b2=[2512];al=ll210J;zs2=roots(b2)ps2=roots(a2)subplot(3,2,3)plot(real(zs2),imag(zs2)/o,);gridaxis(卜1010-1010J)subplot(3,2,4)plot(real(ps2)
7、,imag(ps2)/xf);gridaxis(l-1010-1010J)b3=[l5];a3=[l23];zs3=roots(b3)ps3=roots(a3)subplot(3,2,5)plot(real(zs3),imag(zs3)/or);gridaxis([-1010-1()1()])subplot(3,2,6)plot(real(ps3),imag(ps3)/x,);gridFileEditViewInsertToolsDesktopWindowHelp3nad殳％報⑥心尿▼□目■旦10axis([-1010-1010])^^jFigure1n叵X-1O1O100-505100•50
8、510-5-1005101.对1中每个有理表达式，确定它们的收敛域。解：代码如下i.Re(s)>-l;ii.Re(s)为任意值i.Re(s)>-22.对输入和输出满足下面微分方程：如_3曲)=単+2如+5无⑴dtdrdt的因果LTI系统，求系统函数的零点和极点，并完成一幅适当标注的零极点图。解：代码如下b=[l-3]；a=[125];zs=roots(b)ps=roots(a)subplot(2,1,1)plot(real(zs),imag(zs),'o');gridaxis([-1010-1010])subplot(2,1,2)plot(real(ps),imag(ps)；x');grida
9、xis([-1010-1010])中等题下面的练习需耍用到函数plotpz,其具体实现如下。函数plotpz是用于对分子和分母多项式的系数存入向量a和b屮的LTI系统画岀零极点图用的。这个函数除了画图以外，还产生零极点的值。function[ps,zs]=plotpz(b,a,ROC)ps=roots(a);zs=roots(b);ps=ps(:);zs=zs(:);Maxl=max(abs(im
