1、1.1 市场与市场营销市场与市场营销1.2 我国汽车市场的发展与现状我国汽车市场的发展与现状复习思考题复习思考题实验9 离散系统的零极点分析一、实验目的一、实验目的(1)了解离散系统的零极点与系统因果性和稳定性的关系。(2)观察离散系统零极点对系统冲激响应的影响。(3)熟悉MATLAB中进行离散系统零极点分析的常用子函数。二、实验涉及的二、实验涉及的MATLAB子函数子函数1.zplane功能:功能:显示离散系统的零极点分布图。调用格式:调用格式:zplane(z,p);绘制由列向量z确定的零点、列向量p确定的极点构成的零极点分布图。zplane(b,a);绘制由行向量b和a构成的系统函数确定
2、的零极点分布图。hz,hp,htzplane(z,p);执行后可得到3个句柄向量:hz为零点线句柄,hp为极点线句柄,ht为坐标轴、单位圆及文本对象的句柄。2.roots功能:功能:求多项式的根。调用格式:调用格式:rroots(a);由多项式的分子或分母系数向量求根向量。其中,多项式的分子或分母系数按降幂排列,得到的根向量为列向量。三、实验原理三、实验原理1.离散系统的因果性和稳定性离散系统的因果性和稳定性1)因果系统由理论分析可知,一个离散系统的因果性在时域中必须满足的充分必要条件是:h(n)0 n0即系统的冲激响应必须是右序列。在变换域,极点只能在z平面上一个有界的以原点为中心的圆内。如
3、果系统函数是一个多项式,则分母上z的最高次数应大于分子上z的最高次数。2)稳定系统在时域中,离散系统稳定的充分必要条件是:它的冲激响应绝对可加,即在变换域,则要求所有极点必须在z平面上以原点为中心的单位圆内。0nh(n)3)因果稳定系统综合系统的因果性和稳定性两方面的要求可知,一个因果稳定系统的充分必要条件是:系统函数的全部极点必须在z平面上以原点为中心的单位圆内。2.系统极点的位置对系统响应的影响系统极点的位置对系统响应的影响系统极点的位置对系统响应有着非常明显的影响。下面举例说明系统的极点分别是实数和复数时的情况,使用MATLAB提供的zplane子函数制作零极点分布图进行分析。例例9-1
4、 研究z右半平面的实数极点对系统响应的影响。已知系统的零-极点增益模型分别为:求这些系统的零极点分布图以及系统的冲激响应,判断系统的稳定性。1.5zz(z)H,1zz(z)H,0.85zz(z)H321解解 根据公式写出zpk形式的列向量,求系统的零极点分布图以及系统的冲激响应。程序如下:%在右半平面的实数极点的影响z10;p10.85;k1;b1,a1zp2tf(z1,p1,k);subplot(3,2,1),zplane(z1,p1);ylabel(极点在单位圆内);subplot(3,2,2),impz(b1,a1,20);z20;p21;b2,a2zp2tf(z2,p2,k);subp
5、lot(3,2,3),zplane(z2,p2);ylabel(极点在单位圆上);subplot(3,2,4),impz(b2,a2,20);z30;p31.5;b3,a3zp2tf(z3,p3,k);subplot(3,2,5),zplane(z3,p3);ylabel(极点在单位圆外);subplot(3,2,6),impz(b3,a3,20);由图9-1可见,这3个系统的极点均为实数且处于z平面的右半平面。由图可知,当极点处于单位圆内,系统的冲激响应曲线随着频率的增大而收敛;当极点处于单位圆上,系统的冲激响应曲线为等幅振荡;当极点处于单位圆外,系统的冲激响应曲线随着频率的增大而发散。图9
6、-1 处于z右半平面的实数极点对系统响应的影响例例9-2 研究z左半平面的实数极点对系统响应的影响。已知系统的零-极点增益模型分别为求这些系统的零极点分布图以及系统的冲激响应,判断系统的稳定性。解解 根据公式写出zpk形式的列向量,求系统的零极点分布图以及系统的冲激响应。程序如下:1.5zz(z)H,1zz(z)H,0.85zz(z)H321%在左半平面的实数极点的影响z10;p10.85;k1;b1,a1zp2tf(z1,p1,k);subplot(3,2,1),zplane(z1,p1);ylabel(极点在单位圆内);subplot(3,2,2),impz(b1,a1,20);z20;p
7、21;b2,a2zp2tf(z2,p2,k);subplot(3,2,3),zplane(z2,p2);ylabel(极点在单位圆上);subplot(3,2,4),impz(b2,a2,20);z30;p31.5;b3,a3zp2tf(z3,p3,k);subplot(3,2,5),zplane(b3,a3);ylabel(极点在单位圆外);subplot(3,2,6),impz(z3,p3,20);由图9-2可见,这3个系统的极点均为实数且处于z平面的左半平面。由图可知,当极点处于单位圆内,系统的冲激响应曲线随着频率的增大而收敛;当极点处于单位圆上,系统的冲激响应曲线为等幅振荡;当极点处于
8、单位圆外,系统的冲激响应曲线随着频率的增大而发散。图9-2 处于z左半平面的实数极点对系统响应的影响例例9-3 研究z右半平面的复数极点对系统响应的影响。已知系统的零-极点增益模型分别为求这些系统的零极点分布图以及系统的冲激响应,判断系统的稳定性。j)1j)(z1(z0.3)z(z(z)Hj)0.80.6j)(z0.80.6(z0.3)z(z(z)Hj)0.70.5j)(z0.70.5(z0.3)z(z(z)H321解解 根据公式写出zpk形式的列向量,求系统的零极点分布图以及系统的冲激响应。程序如下:%复数极点的影响z10.3,0;p10.50.7j,0.50.7j;k1;b1,a1zp2t
9、f(z1,p1,k);subplot(3,2,1),zplane(b1,a1);ylabel(极点在单位圆内);subplot(3,2,2),impz(b1,a1,20);z20.3,0;p20.60.8j,0.60.8j;b2,a2zp2tf(z2,p2,k);subplot(3,2,3),zplane(b2,a2);ylabel(极点在单位圆上);subplot(3,2,4),impz(b2,a2,20);z30.3,0;p31j,1j;b3,a3zp2tf(z3,p3,k);subplot(3,2,5),zplane(b3,a3);ylabel(极点在单位圆外);subplot(3,2,
10、6),impz(b3,a3,20);由图9-3可见,这3个系统的极点均为复数且处于z平面的右半平面。由图可知,当极点处于单位圆内,系统的冲激响应曲线随着频率的增大而收敛;当极点处于单位圆上,系统的冲激响应曲线为等幅振荡;当极点处于单位圆外,系统的冲激响应曲线随着频率的增大而发散。由系统的极点分别为实数和复数的情况,我们可以得到结论:系统只有在极点处于单位圆内才是稳定的。3.系统的因果稳定性实例分析系统的因果稳定性实例分析在MATLAB中提供了roots子函数,用于求多项式的根。配合使用zplane子函数制作零极点分布图,可以帮助我们进行系统因果稳定性的分析。例例9-4 已知离散时间系统函数为求
11、该系统的零极点及零极点分布图,并判断系统的因果稳定性。1z2.5z1zH(z)2图9-3 处于z右半平面的复数极点对系统响应的影响解解 该题给出的公式是按z的降幂排列。MATLAB程序如下:b0,1,1;a1,2.5,1;rzroots(b)%求系统的零点rproots(a)%求系统的极点subplot(1,2,1),zplane(b,a);%求系统的零极点分布图title(系统的零极点分布图);subplot(1,2,2),impz(b,a,20);title(系统的冲激响应);xlabel(n);ylabel(h(n);程序运行结果如下:rz 1rp 2.0000 0.5000零极点分布图
12、如图9-4所示,系统的冲激响应如图9-5所示。由运行结果和图9-4可见,该系统有一个极点rp12,在单位圆外;由图9-5可见,该系统的冲激响应曲线随着n增大而发散。因此,该系统不是因果稳定系统。图9-4 例9-4的零极点分布图图9-5 例9-4系统的冲激响应图例例9-5 已知离散时间系统函数为求该系统的零极点及零极点分布图,并判断系统的因果稳定性。解解 MATLAB程序如下:b0.2,0.1,0.3,0.1,0.2;a1,1.1,1.5,0.7,0.3;rzroots(b)rproots(a)43214321z0.3z0.7z1.5z1.11z0.2z0.1z0.3z0.10.2H(z)sub
13、plot(1,2,1),zplane(b,a);title(系统的零极点分布图);subplot(1,2,2),impz(b,a,20);title(系统的冲激响应);xlabel(n);ylabel(h(n);程序运行结果如下:rz 0.50000.8660i 0.50000.8660i 0.25000.9682i 0.25000.9682irp 0.23670.8915i 0.23670.8915i 0.31330.5045i 0.31330.5045i零极点分布图如图9-6所示,系统的冲激响应如图9-7所示。图9-6 例9-5的零极点分布图图9-7 例9-5系统的冲激响应图由运行结果和图
14、9-6零极点分布图可见,该系统的所有极点均在单位圆内;由图9-7可见,该系统的冲激响应曲线随着n增大而收敛。因此,该系统是一个因果稳定系统。四、实验任务四、实验任务(1)阅读并输入实验原理中介绍的例题程序,理解每一条语句的含义,观察程序输出结果,理解零极点对系统特性的影响。(2)已知系统的零-极点增益模型分别为j)0.80.6j)(z0.80.6(z0.3)(z(z)Hj)0.70.5j)(z0.70.5(z0.3)(z(z)H21求这些系统的零极点分布图以及系统的冲激响应,并判断系统的因果稳定性。j)1j)(z1(z0.3)(z(z)H3(3)已知离散时间系统函数分别为求该系统的零极点及零极点分布图,并判断系统的因果稳定性。4)2)(z(z3)1)(z(z5(z)H13213212z0.4z0.35z0.41z4z1.6z1.64(z)H1113z0.511z0.511z12(z)H五、实验预习五、实验预习(1)认真阅读实验原理,明确本次实验任务,读懂各函数和例题程序,了解实验方法。(2)根据实验任务预先编写实验程序。(3)预习思考题:因果稳定的离散系统必须满足的充分必要条件是什么?MATLAB提供了哪些进行零极点求解的子函数?如何使用?六、实验报告六、实验报告(1)列写调试通过的实验程序及运行结果。(2)思考题:回答实验预习思考题。系统函数零极点的位置与系统冲激响应有何关系?