实验Z变换、零极点分析

《实验Z变换、零极点分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验二Z变换、离散系统零极点分布和频率分析实验目的1.学会运用MATLAB求离散时间信号的z变换和z反变换；2．学会运用MATLAB分析离散时间系统的系统函数的零极点；3.学会运用MATLAB分析系统函数的零极点分布与其时域特性的关系；4.学会运用MATLAB进行离散时间系统的频率特性分析。实验内容一、实验原理及实例分析（一）离散时间信号的Z变换1．利用MATLAB实现域的部分分式展开式MATLAB的信号处理工具箱提供了一个对进行部分分式展开的函数residuez()，其调用形式为：[r,p,k]=residuez(num,den)式中，num和den分别为的

2、分子多项式和分母多项式的系数向量，为部分分式的系数向量，p为极点向量，k为多项式的系数向量。【实例1】利用MATLAB计算的部分分式展开式。解：利用MATLAB计算部分分式展开式程序为%部分分式展开式的实现程序num=[18];den=[183-4-1];[r,p,k]=residuez(num,den)2．Z变换和Z反变换MATLAB的符号数学工具箱提供了计算Z变换的函数ztrans()和Z反变换的函数iztrans（），其调用形式为上面两式中，右端的和分别为时域表示式和域表示式的符号表示，可应用函数sym来实现，其调用格式为式中，A为待分析的表示式的字符串

3、，S为符号化的数字或变量。【实例2】求（1）指数序列的Z变换;(2)的Z反变换。解（1）Z变换的MATLAB程序%Z变换的程序实现f=sym('a^n');F=ztrans(f)程序运行结果为：z/a/(z/a-1)可以用simplify()化简得到 :-z/(-z+a)（2）Z反变换的MATLAB程序%Z反变换实现程序F=sym('a*z/(z-a)^2');f=iztrans(F)程序运行结果为f=a^n*n（二）系统函数的零极点分析1.系统函数的零极点分布离散时间系统的系统函数定义为系统零状态响应的z变换与激励的z变换之比，即（3-1）如果系统函数的有理

4、函数表示式为：（3-2）那么，在MATLAB中系统函数的零极点就可通过函数roots得到，也可借助函数tf2zp得到，tf2zp的语句格式为：[Z,P,K]=tf2zp(B,A)其中，B与A分别表示的分子与分母多项式的系数向量。它的作用是将的有理分式表示式转换为零极点增益形式，即：（3-3）【实例3】已知一离散因果LTI系统的系统函数为试用MATLAB命令求该系统的零极点。解：用tf2zp函数求系统的零极点，MATLAB源程序为B=[1,0.32];A=[1,1,0.16];[R,P,K]=tf2zp(B,A)R=-0.3200P=-0.8000-0.2000

5、K=1因此，零点为，极点为与。若要获得系统函数的零极点分布图，可直接应用zplane函数，其语句格式为：zplane(B,A)其中，B与A分别表示的分子和分母多项式的系数向量。它的作用是在Z平面上画出单位圆、零点与极点。【实例4】已知一离散因果LTI系统的系统函数为，试用MATLAB命令绘出该系统的零极点分布图。解：用zplane函数求系统的零极点，MATLAB源程序为B=[1,0,-0.36];A=[1,-1.52,0.68];zplane(B,A),gridonlegend('零点','极点')title('零极点分布图')程序运行结果如图3-1所示。可见

6、，该因果系统的极点全部在单位圆内，故系统是稳定的。图3-1零极点分布图2、系统函数的零极点分布与其时域特性的关系与拉氏变换在连续系统中的作用类似，在离散系统中，z变换建立了时域函数与z域函数之间的对应关系。因此，z变换的函数从形式可以反映的部分内在性质。我们仍旧通过讨论的一阶极点情况，来说明系统函数的零极点分布与系统时域特性的关系。【实例5】试用MATLAB命令画出现下列系统函数的零极点分布图、以及对应的时域单位取样响应的波形，并分析系统函数的极点对时域波形的影响。（1）解：MATLAB源程序为b3=[1,0];a3=[1,-1.2,0.72];subplot

7、(1,2,1)zplane(b3,a3)title('极点在单位圆内的共轭复数')subplot(1,2,2)impz(b3,a3,30);gridon;figure程序运行结果分别如图3-2的（a）所示。(a)图3-2系统函数的零极点分布与其时域特性的关系当极点位于单位圆内时，为衰减序列；当极点位于单位圆上时，为等幅序列；当极点位于单位圆外时，为增幅序列。若有一阶实数极点，则为指数序列；若有一阶共轭极点，则为指数振荡序列；若的极点位于虚轴左边，则序列按一正一负的规律交替变化。（三）离散时间LTI系统的频率特性分析对于因果稳定的离散时间系统，如果激励序列为正弦

8、序列，则系统的稳态响应为。其中，通常是