北京理工大学信号与系统实验报告5-连续时间系统的复频域分析

实验5连续时间系统的复频域分析（综合型实验）一、实验目的1）掌握拉普拉斯变换及其反变换的定义并掌握MATLAB实现方法。2）学习和掌握连续时间系统函数的定义及复频域分析方法。3）掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。二、实验原理与方法1.拉普拉斯变换连续时间信号x(t)的拉普拉斯变换定义为(s)(t)estXxdt（1）拉普拉斯反变换为1(t)(s)e2jstjxXdsj（2）MATLAB中相应函数如下：(F)Llaplace符号表达式F拉氏变换，F中时间变量为t，返回变量为s的结果表达式。(F,t)Llaplace用t替换结果中的变量s。()FilaplaceL以s为变量的符号表达式L的拉氏反变换，返回时间变量为t的结果表达式。(,)FilaplaceLx用x替换结果中的变量t。拉氏变换还可采用部分分式法，当(s)X为有理分式时，它可以表示为两个多项式之比：110110...(s)(s)(s)...MMMMNNNNbsbsbNXDasasa（3）上式可以采用部分分式法展成以下形式1212(s)...NNrrrXspspsp（4）再通过查找常用拉氏变换对易得反变换。利用residue函数可将X(s)展成（4）式形式，调用格式为：[r,p,k]residue(b,a)其中b、a为分子和分母多项式系数向量，r、p、k分别为上述展开式中的部分分式系数、极点和直项多项式系数。2.连续时间系统的系统函数连续时间系统的系统函数是指系统单位冲激响应的拉氏变换(s)(t)estHhdt（5）连续时间系统的系统函数还可以由系统输入与输出信号的拉氏变换之比得到。(s)(s)/X(s)HY（6）单位冲激响应(t)h反映了系统的固有性质，而(s)H从复频域反映了系统的固有性质。由（6）描述的连续时间系统，其系统函数为s的有理函数110110...(s)...MMMMNNNNbsbsbHasasa（7）3.连续时间系统的零极点分析系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s平面上，零点用表示，极点用表示，这样得到的图形为零极点分布图。可以通过利用MATLAB中的求多项式根的roots函数来实现对（7）分子分母根的求解，调用格式如下：r=roots(c)，c为多项式的系数向量，返回值r为多项式的根向量。求取零极点以及绘制系统函数的零极点分布图可以采用pzmap函数，调用格式如下：pzmap(sys)绘出由系统模型sys描述的系统的零极点分布图。[p,z]=pzmap(sys)这种调用方式返回极点与零点，不绘出零极点分布图。还有两个专用函数tf2zp和zp2tf可实现系统的传递函数模型和零极点增益模型的转换。调用格式如下：[z,p,k]=tf2zp(b,a)[b,a]=tf2zp(z,p,k)研究系统函数的零极点分布不仅可以了解系统冲激响应的形式，还可以了解系统的频率特性以及判断系统的稳定性。1）零极点分布与冲激响应的关系系统的极点位置决定着系统冲激响应h(t)的波形，冲激响应的幅值是由系统函数的零点和极点共同确定的，系统的零点位置只影响冲激响应的幅度和相位，不影响波形。2）零极点分布与系统频率响应的关系系统函数的零极点分布不仅决定了系统函数H(s)，也决定了系统的频率响应()H，根据系统的零极点分布情况，可以由几何矢量法分析系统的频率响应。3）零极点分布与系统稳定性的关系稳定性是系统的固有性质，与激励信号无关，由于系统函数(s)H包含了系统的所固有的性质，因而可以根据系统函数的零极点分布判断系统的稳定性。因果系统稳定的充要条件是(s)H的全部极点位于s的左半平面。三．实验内容（1）已知系统的冲激响应(t)u(t)u(t2)h，输入信号(t)u(t)x，试采用复频域的方法求解系统的响应，编写MATLAB程序实现。代码：%DFTfifth_2_1.msymsth=heaviside(t)-heaviside(t-2);x=heaviside(t);H=laplace(h);X=laplace(x);Y=H*X;y=ilaplace(Y)DFTfifth_2_1y=t-heaviside(t-2)*(t-2)所以系统的响应为y(t)=t-(t-2)*u(t-2)（2）已知因果连续时间系统的系统函数分别如下：1）321(s)221Hsss2）54321(s)23332Hsssss试采用MATLAB绘出其零极点分布图，求解系统的冲激响应h(t)和频率响应()H，并判断系统是否稳定。1）b=[1];a=[1221];sys=tf(b,a);[p,z]=pzmap(sys)p=-1.0000-0.5000+0.8660i-0.5000-0.8660iz=Emptymatrix:0-by-1pzmap(sys)-1.4-1.2-1-0.8-0.6-0.4-0.20-1-0.8-0.6-0.4-0.200.20.40.60.81Pole-ZeroMapRealAxis(seconds-1)ImaginaryAxis(seconds-1)所有极点都位于s平面的左半平面，所以系统是稳定的。symssHs=1/(s^3+2*s^2+2*s+1);h=ilaplace(Hs)h=exp(-t)-exp(-t/2)*(cos((3^(1/2)*t)/2)-(3^(1/2)*sin((3^(1/2)*t)/2))/3)所以系统的冲激响应为2333(t)[e(costsint)]u(t)232tthe绘制时域和频域的曲线：b=[1];a=[1221];sys=tf(b,a);subplot(311);impulse(sys);xlabel('t');title('h(t)');subplot(312);[H,w]=freqs(b,a);plot(w,abs(H));xlabel('w');ylabel('Magnitude');title('abs(H)');subplot(313);plot(w,angle(H));xlabel('w');ylabel('phase');title('phase(H)');2）b=[101];a=[12-3332];sys=tf(b,a)sys=s^2+1-------------------------------------s^5+2s^4-3s^3+3s^2+3s+2Continuous-timetransferfunction.02468101214-0.500.5h(t)t(seconds)Amplitude01234567891000.51wMagnitudeabs(H)012345678910-4-2024wphasephase(H)[p,z]=pzmap(sys)p=-3.17040.9669+0.9540i0.9669-0.9540i-0.3817+0.4430i-0.3817-0.4430iz=0+1.0000i0-1.0000ipzmap(sys)由于s平面有半平面有极点，所以是不稳定系统。绘制冲激响应和频域响应的图形方法同上一题图形如下：（3）已知连续时间系统函数的极点位置分别如下所示（设系统无零点）：-3.5-3-2.5-2-1.5-1-0.500.51-1-0.8-0.6-0.4-0.200.20.40.60.81Pole-ZeroMapRealAxis(seconds-1)ImaginaryAxis(seconds-1)01234567891000.20.40.60.8wMagnitudeabs(H)012345678910-2-1012wphasephase(H)010203040506070-3-2-101x1028h(t)t(seconds)Amplitude分别绘制以下六种不同情况下，系统函数的零极点分布图，并绘制相应冲激响应的时域波形，观察并分析系统函数极点位置对冲激响应时域特性的影响。1）p=0b=[1];a=[10];sys=tf(b,a)sys=1-sContinuous-timetransferfunction.pzmap(sys)1(s)(t)u(t)Hhssymsth=heaviside(t);ezplot(h,[-55])title('h(t)')2）p=-2b=[1];a=[12];sys=tf(b,a)sys=1-----s+2Continuous-timetransferfunction.pzmap(b,a)-1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81Pole-ZeroMapRealAxis(seconds-1)ImaginaryAxis(seconds-1)-5-4-3-2-101234500.20.40.60.81th(t)Pole-ZeroMapRealAxis(seconds-1)ImaginaryAxis(seconds-1)-2-1.8-1.6-1.4-1.2-1-0.8-0.6-0.4-0.20-1-0.8-0.6-0.4-0.200.20.40.60.8121(s)(t)e*(t)2tHhussymsth=exp(-2*t)*heaviside(t);ezplot(h)3)p=2b=[1];a=[1-2];sys=tf(b,a)sys=1-----s-2Continuous-timetransferfunction.pzmap(b,a)21(s)(t)*(t)2tHheussymsth=exp(2*t)*heaviside(t);ezplot(h)00.511.522.500.10.20.30.40.50.60.70.80.91texp(-2t)heaviside(t)00.20.40.60.811.21.41.61.82-1-0.8-0.6-0.4-0.200.20.40.60.81Pole-ZeroMapRealAxis(seconds-1)ImaginaryAxis(seconds-1)01234560246810x104texp(2t)heaviside(t)4)122,2pjpjb=[1];a=[104];sys=tf(b,a)sys=1-------s^2+4Continuous-timetransferfunction.pzmap(b,a)211(s)(t)sin(2t)*(t)42Hhussymsth=(1/2)*sin(2*t)*heaviside(t);ezplot(h)5)1214,14pjpjb=[1];a=[1217];sys=tf(b,a)sys=1--------------s^2+2s+17Continuous-timetransferfunction.pzmap(b,a)-1-0.8-0.6-0.4-0.200.20.40.60.81-2.5-2-1.5-1-0.500.511.522.5Pole-ZeroMapRealAxis(seconds-1)ImaginaryAxis(seconds-1)0123456-0.6-0.4-0.200.20.40.6t(sin(2t)heaviside(t))/2-1.4-1.2-1-0.8-0.6-0.4-0.20-4-3-2-101234Pole-ZeroMapRealAxis(seconds-1)ImaginaryAxis(seconds-1)211(s)(t)sin(2)(t)21