1、利用residue函数可将X(s)展成(4)式形式,调用格式为:其中b、a为分子和分母多项式系数向量,r、p、k分别为上述展开式中的部分分式系数、极点和直项多项式系数。2。连续时间系统的系统函数连续时间系统的系统函数是指系统单位冲激响应的拉氏变换 (5)连续时间系统的系统函数还可以由系统输入与输出信号的拉氏变换之比得到。 (6)单位冲激响应反映了系统的固有性质,而从复频域反映了系统的固有性质.由(6)描述的连续时间系统,其系统函数为的有理函数 (7)3.连续时间系统的零极点分析系统的零点指使式(7)的分子多项式为零的点,极点指使分母多项式为零的点,零点使系统的值为零,极点使系统的值为无穷大。通
2、常将系统函数的零极点绘在s平面上,零点用表示,极点用表示,这样得到的图形为零极点分布图。可以通过利用MATLAB中的求多项式根的roots函数来实现对(7)分子分母根的求解,调用格式如下:r=roots(c),c为多项式的系数向量,返回值r为多项式的根向量。求取零极点以及绘制系统函数的零极点分布图可以采用pzmap函数,调用格式如下:pzmap(sys)绘出由系统模型sys描述的系统的零极点分布图。p,z=pzmap(sys)这种调用方式返回极点与零点,不绘出零极点分布图.还有两个专用函数tf2zp和zp2tf可实现系统的传递函数模型和零极点增益模型的转换。调用格式如下:z,p,k=tf2zp
3、(b,a)b,a=tf2zp(z,p,k)研究系统函数的零极点分布不仅可以了解系统冲激响应的形式,还可以了解系统的频率特性以及判断系统的稳定性。1)零极点分布与冲激响应的关系系统的极点位置决定着系统冲激响应h(t)的波形,冲激响应的幅值是由系统函数的零点和极点共同确定的,系统的零点位置只影响冲激响应的幅度和相位,不影响波形。2)零极点分布与系统频率响应的关系系统函数的零极点分布不仅决定了系统函数H(s),也决定了系统的频率响应,根据系统的零极点分布情况,可以由几何矢量法分析系统的频率响应。3)零极点分布与系统稳定性的关系稳定性是系统的固有性质,与激励信号无关,由于系统函数包含了系统的所固有的性
4、质,因而可以根据系统函数的零极点分布判断系统的稳定性。因果系统稳定的充要条件是的全部极点位于s的左半平面。三实验内容(1)已知系统的冲激响应,输入信号,试采用复频域的方法求解系统的响应,编写MATLAB程序实现。代码:DFTfifth_2_1。msyms th=heaviside(t)heaviside(t2);x=heaviside(t);H=laplace(h);X=laplace(x);Y=HX;y=ilaplace(Y) DFTfifth_2_1 y =t heaviside(t 2)(t - 2)所以系统的响应为y(t)=t(t-2)u(t2)(2)已知因果连续时间系统的系统函数分别
5、如下:1)2)试采用MATLAB绘出其零极点分布图,求解系统的冲激响应h(t)和频率响应,并判断系统是否稳定。 b=1; a=1 2 2 1; sys=tf(b,a); p,z=pzmap(sys)p = 1.0000 0。5000 + 0。8660i5000 - 0.8660iz = Empty matrix: 0by1 pzmap(sys)所有极点都位于s平面的左半平面,所以系统是稳定的。 syms s Hs=1/(s3+2*s2+2s+1); h=ilaplace(Hs)h =exp(-t) exp(-t/2)*(cos(3(1/2)t)/2) - (3(1/2)*sin(3(1/2)t
6、)/2)/3)所以系统的冲激响应为绘制时域和频域的曲线:b=1;a=1 2 2 1;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);abs(H)subplot(313);plot(w,angle(H);xlabel(wphase);title(phase(H); b=1 0 1; a=1 2 -3 3 3 2; sys=tf(b,a)sys = s2 + 1 - s5 + 2 s4 -
7、3 s3 + 3 s2 + 3 s + 2Continuoustime transfer function。 p,z=pzmap(sys) -3.1704 0.9669 + 0。9540i 0。9669 - 0.9540i -0。3817 + 0.4430i3817 - 0.4430i 0 + 1。0000i 0 - 1。 pzmap(sys)由于s平面有半平面有极点,所以是不稳定系统.绘制冲激响应和频域响应的图形方法同上一题图形如下:(3)已知连续时间系统函数的极点位置分别如下所示(设系统无零点):分别绘制以下六种不同情况下,系统函数的零极点分布图,并绘制相应冲激响应的时域波形,观察并分析系
8、统函数极点位置对冲激响应时域特性的影响.1)p=0 b=1; a=1 0; sys=tf(b,a) 1 sContinuous-time transfer function。 syms t h=heaviside(t); ezplot(h,5 5) title(h(t)2)p=2 b=1; a=1 2; sys=tf(b,a) - s + 2Continuoustime transfer function. pzmap(b,a) syms t h=exp(2t)heaviside(t); ezplot(h)3)p=2 b=1; a=1 -2; - s 2 pzmap(b,a) h=exp(2*
9、t)heaviside(t);4) b=1; a=1 0 4; - s2 + 4 h=(1/2)*sin(2t)heaviside(t); ezplot(h)5) b=1; a=1 2 17; sys=tf(b,a) - s2 + 2 s + 17 pzmap(b,a) syms t h=(1/2)exp(t)sin(2t)*heaviside(t); ezplot(h)6) b=1; a=1 2 17; sys=tf(b,a) - s2 2 s + 17 pzmap(b,a) h=exp(t)sin(2*t)heaviside(t)/2;极点在左半平面时呈衰减趋势,在左半平面坐标轴上时呈指数
10、衰减,在非坐标轴位置上时成衰减振荡;在右半平面时成增加趋势,在右半平面坐标轴上时呈增加趋势,在非坐标轴上时呈增幅振荡;在纵轴上时,在非原点时呈等幅振荡,在原点时为单位阶跃响应。(4)已知三个连续时间系统的系统函数,极点相同,零点不同,试用MATLAB分别绘制系统的零极点分布图及相应冲激响应的时域波形,观察并分析系统函数零点位置对冲激响应时域特性的影响。 b=1; a=1 2 17; -Continuous-time transfer function. pzmap(b,a) h=(1/4)*exp(t)sin(4t)*heaviside(t); axis tight b=1 8; s + 8
11、- syms s H=(s+8)/(s2+2s+17); h=ilaplace(H)exp(t)(cos(4t) + (7*sin(4t)/4) clearh=exp(-t)(cos(4*t)+(7sin(4t)/4)heaviside(t); axis tight3) b=1 8; a=1 2 17; s - 8 - H=(s-8)/(s2+2s+17); h=ilaplace(H)exp(-t)*(cos(4*t) (9sin(4*t)/4)h=exp(-t)*(cos(4t)-(9sin(4t)/4)heaviside(t);极点位置不变时,响应是衰减还是增幅不会发生改变。但随着零点位置的变化,响应的振幅和相位会发生变化,零点由原点向正半轴或负半轴移动时,响应时域波形振幅增大.四、实验收获与体会本次实验用MATLAB实现了和两个变换,掌握了连续时间系统函数的复频域分析法,加深了对系统零极点分布于系统特性关系的理解。