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

信号与线性系统实验报告姓名：xx学号：xxxxxxx2015年12月信号与线性系统实验六第1页共8页实验六连续时间系统的零极点分析实验目的：1、学会用Matlab求解系统函数的零极点；2、学会用Matlab分析系统函数的极点分布与系统稳定性的关系。实验原理：1、系统零极点绘制系统函数H(s)通常是一个有理分式，其分子和分母均为多项式。利用Matlab中的roots函数，可以求出分子和分母多项式的根，即可计算出H(s)的零极点。例如：多项式542)(24ssssN的根可以由下列语句求出：N＝[10245]；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平面的左平面时，系统是稳定的。在绘制好系统零极点分布图后，就可以根据这个知识点判断系统的稳定性。注意：在绘制系统零极点分布图时，可以适当变换坐标的显示范围，来达到增强零极点分布图可读性的效果。实验内容：一、用两种方法绘制如下系统函数的零极点分布图，并且判断系统是否稳定。信号与线性系统实验报告姓名：xx学号：xxxxxxx2015年12月信号与线性系统实验六第2页共8页1、221)(2ssssH方法一，程序代码：零极图：b=[1,-1];a=[122];zs=roots(b);ps=roots(a);plot(real(zs),imag(zs),'ro',real(ps),imag(ps),'gx',…'markersize',12);axis([-3,3,-2,2]);gridon;title('零极点图')xlabel('实部')ylabel('虚部')-3-2-10123-2-1.5-1-0.500.511.52零极点图实部虚部方法二，程序代码：零极图：figure;b=[1,-1];a=[122];sys=tf(b,a);pzmap(sys);axis([-3,3,-2,2]);2、12223)(232ssssssH信号与线性系统实验报告姓名：xx学号：xxxxxxx2015年12月信号与线性系统实验六第3页共8页方法一，程序代码：零极图：b=[132];a=[1221];zs=roots(b);ps=roots(a);plot(real(zs),imag(zs),'o',real(ps),imag(ps),'gx','markersize',12);axis([-20-22]);gridon;title('零极点图')xlabel('实部')ylabel('虚部')-2-1.8-1.6-1.4-1.2-1-0.8-0.6-0.4-0.20-2-1.5-1-0.500.511.52零极点图实部虚部方法二，程序代码：零极图：figure;b=[132];a=[1221];sys=tf(b,a);pzmap(sys)axis([-20-22]);3、22)(2sssH方法一，程序代码：零极图：信号与线性系统实验报告姓名：xx学号：xxxxxxx2015年12月信号与线性系统实验六第4页共8页b=[122];zs=roots(b);plot(real(zs),imag(zs),'o','markersize',12);axis([-20-22]);gridon;title('零极点图')xlabel('实部')ylabel('虚部')-2-1.8-1.6-1.4-1.2-1-0.8-0.6-0.4-0.20-2-1.5-1-0.500.511.52零极点图实部虚部方法二，程序代码：零极图：figure;b=[122];a=[1];sys=tf(b,a);pzmap(sys);axis([-20-22]);二、已知系统函数为321()221Hssss，画出该系统的零极点分布图，分析系统的稳定性，并画出该系统的冲激响应和幅频响应。程序代码：输出曲线图：信号与线性系统实验报告姓名：xx学号：xxxxxxx2015年12月信号与线性系统实验六第5页共8页b=[1];a=[1221];sys=tf(b,a);subplot(3,1,1);pzmap(sys);t=0:0.01:2*pi;y=impulse(b,a,t);subplot(3,1,2);plot(t,y);xlabel(‘t’)ylabel(‘h(t)’)title('冲激响应')[H,w]=freqs(b,a);subplot(3,1,3);plot(w,abs(H));xlabel(‘w’)ylabel(‘|H(jw)|’)title('幅频特性')-1-0.9-0.8-0.7-0.6-0.5-0.4-0.3-0.2-0.10-101Pole-ZeroMapRealAxisImaginaryAxis01234567-0.500.5th(t)冲激响应01234567891000.51w|H(jw)|幅频特性三、已知系统函数为21()21Hssas，画出2,1,0,1,2a时该系统的零极点分布图。如果系统稳定时，画出响应系统的幅频响应，并分析系统极点位置对系信号与线性系统实验报告姓名：xx学号：xxxxxxx2015年12月信号与线性系统实验六第6页共8页统的幅频特性有何影响。零极图，程序代码：零极图：fora=-2:2num=[1];den=[12*a1];sys=tf(num,den);subplot(3,2,a+3);pzmap(sys);end幅频响应，程序代码：响应图：figure;fora=0:2num=[1];den=[12*a1];[H,w]=freqs(num,den);subplot(3,1,a+1);plot(w,abs(H));xlabel(‘w’)ylabel(‘|H(jw)|’)title(‘幅频特性')end0123456789100100200w|H(jw)|幅频特性01234567891000.51w|H(jw)|幅频特性010203040506070809010000.51w|H(jw)|幅频特性信号与线性系统实验报告姓名：xx学号：xxxxxxx2015年12月信号与线性系统实验六第7页共8页-3-2-10123-2-1.5-1-0.500.511.52零极点图实部虚部-2-1.8-1.6-1.4-1.2-1-0.8-0.6-0.4-0.20-2-1.5-1-0.500.511.52零极点图实部虚部-2-1.8-1.6-1.4-1.2-1-0.8-0.6-0.4-0.20-2-1.5-1-0.500.511.52零极点图实部虚部信号与线性系统实验报告姓名：xx学号：xxxxxxx2015年12月信号与线性系统实验六第8页共8页-1-0.9-0.8-0.7-0.6-0.5-0.4-0.3-0.2-0.10-101Pole-ZeroMapRealAxisImaginaryAxis01234567-0.500.5th(t)冲激响应01234567891000.51w|H(jw)|幅频特性0123456789100100200w|H(jw)|幅频特性01234567891000.51w|H(jw)|幅频特性010203040506070809010000.51w|H(jw)|幅频特性