实验5 连续时间系统的复频域分析.docx
实验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中时间变量为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函数可以将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+%夕-16"1+.+Z.v.V-ICInS+NTS+.+%3.连续时间系统的零极点分析系统的零点是指式(7)的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统函数的值无穷大。通常将系统函数的零极点绘在S平面上，零点用°表示，极点用X表示，这样得到的图形称为零极点的分布图。由零极点的定义可知，零点和极点分别指式(7)的分子多项式和分母多项式的根。利用MATLAB求多项式的根可以通过函数roots来实现，该函数的调用格式为：r=roots(c)C为多项式的系数向量，返回值r为多项式的根向量。分别对式(7)的分子多项式和分母多项式求根即可得到零极点。此外，在MATLAB中还提供了更简便的方法来求取零极点和绘制系统函数的零极点分布图，即利用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为传递函数的分子多项式和分母多项式的系数向量,返回值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 .已知因果连续时间系统的系统函数分别如下:(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);plot(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因果稳定系统的充要条件是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);plot(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*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种不同情况下，系统函数的零极点分布图，并绘制相应冲激响应的时域波形，观察并分析系统函数极点位置对冲激响应时域特性的影响。(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(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)subplot(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-2-1.4-1.2-1-0.806-0.4-0.20Real Axifeeconds'1)-410.5A00.50123456(6) P=I+4j,p=l-4j时:b=l;a=l-217;sys=tf(b,a);subplot(211);PZm叩(SyS);%绘制零极点分布图t=0:0.01:15;h=(sin(4*t).*exp(t)4subplot(212);plot(t,h);xlabel('t');title('h(t)')；Pole-ro Mapspuo。8 S) S×< AJeu'&euq由以上六幅图可知，系统函数无零点时，若所有极点均位于S平面的左半平面，那么冲激响应的时域波形是收敛的；若所有极点位于S平面的右半平面，那么冲激响应的时域波形是发散的；若所有极点位于虚轴上，那么冲激响应的时域波形既不发散也不收敛。4.已知3个连续时间系统函数(DHs) =?+5 + 17 H(S) =5 + 852+5 + 17(3)” (S) =s-SS2 +5 + 17上述三个系统具有相同的极点，只是零点不同，试用MATLAB分别绘制系统的零极点分布图及相应冲激响应的时域波形，观察并分析系统函数零点位置对冲激响应时域特性的影响。(l)"(s) =$2 +s + 17b=l;a=l217;sys=tf(b,a);subplot(211);PZmaP(Sys);%绘制零极点分布图Ht=ilaplace(l(s2+2*s+17)subplot(212);symst;ezplot(t,Ht)xlabel(,t')jtitle('h(t)1);Pole Zero Map-1-0.8-0.6-0.4-0.20Real Axieconds'1)23456 H(S) =5 + 852÷5 + 17b=l 8;a=l 2 17; sys=tf(bza);subplot(211);PZmaP(Sys);%绘制零极点分布图 Ht=ilaplace(s+8)(s2+2*s+17) subplot(212);syms t;ezplot(t,Ht)iXabel('t')jtitle(,h(t),);Pole-Zero Map-6-5-4-3Real Axig (seconds'1)-1 O426c-8(3)/(5)=-/°-r+5+17b=l-8;a=l217;sys=tf(b,a);subplot(211);PZm叩(SyS);%绘制零极点分布图Ht=ilaplace(s-8)(s2+2*s+17)subplot(212);symst;ezplot(t,Ht)label('t')jtitle(,h(t),);Pole-Zero Map12345678(SPuooeS) s×< Juew 一-2-4-2由以上三幅图可以看出，当系统极点相同时，零点的不同并不会改变冲激响应在时域的收敛与否；零点不同仅仅影响冲激函数在时域的振幅与相位。四、心得体会本实验主要研究连续时间系统的复频域分析问题。其中分析系统函数的零极点分布可以得到很多结论。系统的极点位置影响冲激函数的波形，而系统的零点只影响冲激响应的幅度与相位；可以根据系统的零极点分布判断系统的稳定性，因果稳定系统的充要条件是H(三)的全部极点位于s平面的左半面。这些结论在上述实验中得到了很好的验证。总之，本实验原理简单易懂，操作起来也不更杂；我在进行实验的过程中体会到了验证原理和巩固已有知识的乐趣。
