实验六连续时间系统的零极点分析.docxVIP

实验六连续时间系统的零极点分析 实验目的： 1、 学会用Matlab求解系统函数的零极点； 2、 学会用Matlab分析系统函数的极点分布与系统稳定性的关系。 实验原理： 1、 系统零极点绘制 系统函数H(s)通常是一个有理分式，其分子和分母均为多项式。利用Matlab 中的roots函数，可以求出分子和分母多项式的根，即可计算出 H(s)的零极点。 例如：多项式N(s) s4 2s2 4s 5的根可以由下列语句求出： N = [1 0 2 4 5] ; r=roots(N); 求出零极点后以零极点的实部和虚部作图， 即可得出零极点的分布图。例如: 执行zs=roots(b) ； ps=roots(a) ； (b，a分别为分子分母多项式系数向量)，再 执行 plot(real(zs),imag(zs), o ,real(ps),imag(ps), x , markersize ,12); 就能够画出系统的零极点分布图。 绘制系统零极点的分布图再 Matlab中还有一种更加简便的方法，即利用函 数pzmap，调用形式为： pzmap(sys) 它表示画出由sys所描述的系统的零极点分布图。利用 sys = tf(b,a)来构建系统 模型，这在实验2中已经介绍过，b,a分别为系统函数H(s)的分子分母多项式系 数向量。 2、 系统函数的零极点与系统的稳定性 根据信号与线性系统中的知识我们知道：当系统函数的极点全部位于 s平面 的左平面时，系统是稳定的。在绘制好系统零极点分布图后，就可以根据这个知 识点判断系统的稳定性。 注意：在绘制系统零极点分布图时，可以适当变换坐标的显示范围，来达到 增强零极点分布图可读性的效果。 实验内容: 、用两种方法绘制如下系统函数的零极点分布图，并且判断系统是否稳定。 “ s 1 1、H(s)— s 2s 2 方法一，程序代码:零极图:b=[1,-1];a=[1 2 2];zs=roots(b);ps=roots(a);plot(real(zs),imag(zs).ro,real(ps),imag(ps). 方法一，程序代码: 零极图: b=[1,-1]; a=[1 2 2]; zs=roots(b); ps=roots(a); plot(real(zs),imag(zs). ro,real(ps),imag(ps). 零极点图 gx, ?? markersize,12); axis([-3,3,-2,2]); grid on; title(零极点图) xlabel(实部) ylabel(虚部) 方法二，程序代码: 零极图: 零极图:零极点图 零极图: 零极点图 2、 H(s)s 2、 H(s) ~3 2 s3 2s2 2s 1 方法一，程序代码： b=[1 3 2]; a=[1 2 2 1]; zs=roots(b); ps=roots(a); plot(real(zs),imag(zs),o,real(ps), imag(ps), gx,markersize,12); axis([-2 0 -2 2]); grid on; title(零极点图) xlabel(实部) ylabel(虚部) 方法二，程序代码： 零极图： figure; b=[1 3 2]; a=[1 2 2 1]; sys=tf(b,a); pzmap(sys) axis([-2 0 -2 2]); Prie-lETa I4s: 1 ￥乍 -i i K 卜 ? ° ■ - — * ] ] IB .1 d .H .1 町| dQQ 41 3、H(s) s 2s 2 方法一，程序代码： b=[1 2 2]; zs=roots(b); plot(real(zs),imag(zs),o,marker size,12); axis([-2 0 -2 2]); grid on; title(零极点图) xlabel(实部) ylabel(虚部) 零极图：21 510.5部 0 零极图： 2 1 5 1 0.5 部 0 虚 -0.5 -1 -2 -1.8 -1.6 -1.4 -1.2 - 1 -0.8 -0.6 -0.4 -0.2 零极点图 实部 零极图: figure; b=[1 2 2]; a=[1]； sys=tf(b,a); pzmap(sys); axis([-2 0 -2 2]); -■D-L电、已知系统函数为H(s)s3 2s -■ D - L电 、已知系统函数为H(s) 系统的稳定性，并画出该系统的冲激响应和幅频响应 程序代码： 输出曲线图： b=[1]； a=[1 2 2 1]; sys=tf(b,a); subplot(3,1,1); pzmap(sys); t=0:0.01:2*pi; y=impu