实验5 连续时间系统的复频域分析.docx

《实验5 连续时间系统的复频域分析.docx》由会员分享，可在线阅读，更多相关《实验5 连续时间系统的复频域分析.docx（11页珍藏版）》请在优知文库上搜索。

1、实验5连续时间系统的复频域分析一、实验目的1、掌握拉普拉斯变换及其反变换的定义，并掌握MATLAB实现方法。2、学习和掌握连续时间系统系统函数的定义及复频域分析方法。3、掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。二、实验原理与方法1、拉普拉斯变换连续时间信号XQ)的拉普拉斯变换定义为拉普拉斯变换定义为X(三)=x(t)e-stdt(1)J-XJ拉普拉斯反变换定义为x(t)-rX(s)es,ds(2)2用J”在MATLAB中，可以采用符号数学工具箱的Iaplace函数和iIaplace函数进行拉氏变换和反拉氏变换。1.=IaPIaCC(F)符号表达式F的拉氏变换，F中时间变量

2、为t,返回变量为S的结果表达式。1.=Iaplace(F,t)用t替换结果中的变量s。F=ilaplace(L)以S为变量的符号表达式L的拉氏反变换,返回时间变量为t的结果表达式。F=iIaplace(L,x)用X替换结果中的变量t。除了上述iIaplace函数，还可以采用部分分式法，求解拉普拉斯逆变换，具体原理如下：当X(三)为有理分式时，它可以表示为两个多项式之比：X(三)=祟=+。D(三)aNs+即_科+劭式(3)可以用部分分式法展成一下形式X(三)=/一+/一+.+(4)Pis-p2s-pN通过查常用拉普拉斯变换对，可以由式(1-2)求得拉普拉斯逆变换。利用MATLAB的residue

3、函数可以将I(三)展成式(1-2)所示的部分分式展开式，该函数的调用格式为：r,p,k=residuc(b,a)其中b、a为分子和分母多项式系数向量，r、p、k分别为上述展开式中的部分分式系数、极点和直项多项式系数。2、连续时间系统的系统函数连续时间系统的系统函数是系统单位冲激响应的拉氏变换HG)=h(t)e-sldt(5)J-OO此外，连续时间系统的系统函数还可以由系统输入和系统输出信号的拉氏变换之比得到H(三)=Y(三)ZX(三)(6)单位冲激响应反映了系统的固有性质，而($)从复频域反映了系统的固有性质。由式(6)描述的连续时间系统，其系统函数为S的有理函数H(S) =CWSAf+%夕-

4、161+.+Z.v.V-ICInS+NTS+.+%3.连续时间系统的零极点分析系统的零点是指式(7)的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统函数的值无穷大。通常将系统函数的零极点绘在S平面上，零点用表示，极点用X表示，这样得到的图形称为零极点的分布图。由零极点的定义可知，零点和极点分别指式(7)的分子多项式和分母多项式的根。利用MATLAB求多项式的根可以通过函数roots来实现，该函数的调用格式为：r=roots(c)C为多项式的系数向量，返回值r为多项式的根向量。分别对式(7)的分子多项式和分母多项式求根即可得到零极点。此外，在MATLAB中还提供了

5、更简便的方法来求取零极点和绘制系统函数的零极点分布图，即利用PZmaP函数，该函数的调用格式为：PZmaP(SyS)绘出由系统模型sys描述的系统的零极点分布图。p,z=pzmap(sys)这种调用方法返回极点和零点，而不绘出零极点分布图。其中SyS为系统传函模型，由t命令sys=tf(b,a)实现，b、a为传递函数的分子多项式和分母多项式的系数向量。MATLAB还为用户提供了两个专用函数tf2zp和zp2tf来实现系统传递函数模型和零极点增益模型的转换，其调用格式为：z,p,k=tf2zp(b,a)b,a=zp2tf(z,p,k)其中b、a为传递函数的分子多项式和分母多项式的系数向量,返回值

6、z为零点列向量,P为极点列向量，k为系统函数零极点形式的增益。三、实验内容1 .已知系统的冲激响应z(f)=(r)-。-2),输入信号XQ)=试采用复频域的方法求解系统的响应，编写MATLAB程序实现。fl=sym(heaviside(t)-heaviside(t-2);Fl=laplace(fl)f2=sym(,heaviside(t)jF2=laplace(f2)Y=Fl*F2jy=ilaplace(Y)symst;ezplot(t,y)ixlabel(t)ititle(,y(t),);结果：Fl=Vs-exp(-2*s)sF2=lsy=t-heaviside(t-2)*(t-2)2 .已

7、知因果连续时间系统的系统函数分别如下:(I)H(S) =s + 2s + 2s +1 H(S)s5 + 2s4-3s3+3s2+3s + 2试采用MATLAB画出其零极点分布图，求解系统的冲激响应力和频率响应HGy),并判断系统是否稳定。HG)=SymSsw;b=la=lz2z2,l;sys=tf(bza);p,z=pzmap(sys)%求取零极点subplot(211);PZmaP(Sys);%绘制零极点分布图Hwl=SUbs(1/(s3+2*s2+2*s+1)zszj*w)Htl=Haplace(1/(s3+2*s2+2*s+l)Hwlw=freqs(bza);subplot(223);p

8、lot(wzabs(Hwl);xlabel(,w,);title(,IH(w),);gridon;subplot(224);plot(w,angle(Hwl);xlabel(,w,);title(,phi(w),);结果：P=-1.0000+0.0000i-0.5000+0.8660i-0.5000-0.8660iz=Emptymatrix:O-by-1Hwl=l(w3*(-i)-2*w2+w*(2*i)+1)Htl=exp(-t)-exp(-t2)*(cos(3(l2)*t)2)-(3A(l/2)*sin(3A(l/2)*t)/2)/3)-0.5Pole-Zero Map0.5-1(w)10

9、因果稳定系统的充要条件是H(三)的全部极点位于S平面的左半平面，由系统的零极点分布图可知，系统的全部极点位于S平面的左半面，则系统稳定。 H(S) =?+254-3?+3?+35 + 2SymSsw;b=l0l;a=l2-3332;sys=tf(b,a);p,z=pzmap(sys)%求取零极点subplot(211);PZmaP(Sys);%绘制零极点分布图Hwl=SUbs(s2+1)(s5+2*s4-3*s3+3*s2+3*s+2)zs,j*w)Htl=Haplace(s2+1)(s5+2*s4-3*s3+3*s2+3*s+2)Hwlw=freqs(bza);subplot(223);pl

10、ot(w,abs(Hwl);xlabel(,w,);title(,IH(w),);gridon;subplot(224);plot(wzangle(Hwl);xlabel(,w,);title(,phi(w),);结果：P =-3.17040.96690.9669-0.3817-0.3817 z =0.00000.0000+ 0.0000i+ 0.9540i-0.9540i+ 0.4430i-0.4430i+ 1.0000i-1.0000iHwl=-(wa2-1)/(w5*i+2*wa4+wa3*(3*i)-3*w2+w*(3*i)+2)Htl=symsum(r32*exp(r3*t)/(5*

11、r34+8*r33-9*r32+6*r3+3)zr3inRootOf(s35+2*s34-3*s33+3*s32+3*s3+2,s3)+symsum(exp(r3*t)(5*r34+8*r33-9*r32+6*r3+3),r3inRootOf(s35+2*s34-3*s33+3*s32+3*s3+2,s3)Pole-Zero Map因果稳定系统的充要条件是H(三)的全部极点位于s平面的左半平面，由系统的零极点分布图可知，系统的部分极点位于S平面的左半面，则系统不稳定。3 .已知连续时间系统函数的极点位置分别如下所示(设系统无零点)，试用MATLAB绘制6种不同情况下，系统函数的零极点分布图，并

12、绘制相应冲激响应的时域波形，观察并分析系统函数极点位置对冲激响应时域特性的影响。(I)P=O时:b=l;a=l0;sys=tf(b,a);subplot(211);PZm叩(SyS);%绘制零极点分布图Ht=ilaplace(ls)subplot(212);symst;ezplot(t,Ht)label(t)jtitle(,h(t),);spu8s) s JBU)eE 一Pole-ZeroM叩10.50-0.5-1-1-0.8-0.6-0.4-0.200.20.40.60.810123456RealAxifeeconds1)21.5A10.50(2) p=2时：b=l;a=l2;sys=tf(

13、bza);subplot(211);PZmaP(Sys);%绘制零极点分布图Ht=ilaplace(1(s+2)subplot(212);symst;0.5Pole-Zero Map 1-2-1.8 -1.6 -1.4 -1.2-1-0.8 -0.6 -0.4 0.20Real AX枭伊CondS-I) 1.5-0.5 -1ezplot(t,Ht)label(t)jtitle(,h(t),);1A0.500.50123456(3) p=2时Sb=l;a=l-2;sys=tf(bza);subplot(211);PZmaP(Sys);%绘制零极点分布图t=0:0.01:3;h=exp(2*t)s

14、ubplot(212);plot(t,h);xlabel(t)；title(h(t)；Pole-ZeroMap0.5、0-0.5LfflIiflflIf-00.20.40.60.811.21.41.61.82RealAX枭俨CondS-B600400(4) p=2j,p=2j时：b=l;a=l04;sys=tf(bza);subplot(211);PZmaP(Sys);%绘制零极点分布图Ht=ilaplace(l(s2+4)subplot(212);symst;ezplot(t,Ht)jxlabel(,t)ititle(,h(t),);Pole-ro Map(5) p=-l+4j,p=-l-4jW:b=l;a=l217;sys=tf(bza);subplot(211);PZmaP(Sys);%绘制零极点分布图Ht=ilaplace(l/(sA2+2*s+17)subplot(212);symst;ezplot(t,Ht);XIabeIPole-roMap420-