第十八章 离散系统的Matlab Z域分析仿真,一、离散系统的系统函数,离散序列的Z变换定义 离散系统函数的定义及其零极点表示 设某离散系统函数为 零极点的求法：roots( )命令,表示系统函数分子、分母多项式系数时注意：,系统函数按Z的降幂排列时 系数向量应由最高次项系数开始，直到常数项，缺项补零。 系统函数按Z1的升幂排列时 分子、分母多项式应保证维数相同，缺项补零。,例如：,1） 则：A=1,2,-3,7,5;B=3,-5,10,0 2） A=2,-5,2;B=1,-0.5, 0,零极点z平面作图实用函数：ljdt.m,例：已知系统函数为： 绘制其零极点图。 由系统函数得： A=1 3 7 5;B=3 5 10 0 可调用ljdt(A,B)即可。,二、离散系统性能分析,零极点的分布对于系统性能分析有重要意义，可分析以下方面的系统特性： 系统的稳定性； 系统单位样值响应的时域特性； 系统的频率特性；,1 零极点分布于系统的稳定性,判定方法： 时域：绝对可和条件。 Z域：系统函数H（z）的所有极点在z平面的单位圆内。 实现： 只要调用ljdt（）直观判断极点位置就可以确定系统稳定性，并可使用abs（）求出各极点至原点的长度。,例：已知系统函数为: 判定系统的稳定性。 程序如ex1201.m所示。,2 由系统函数分析系统时域特性,理论分析 实现方法： 由系统函数H(z)的表达式可以得到分子、分母系数向量A和B，调用impz(B,A)即可得到单位样值时域响应。 例：已知系统函数为：,试分析 时的单位响应； 程序如ex1202.m所示。,3 离散系统频率响应特性,频率响应的概念 Matlab实现离散系统频率特性分析 1）直接法 提供了专门由系统函数求解频率响应的函数freqz（），调用格式如下： H,w=freqz(b,a,n) H,w=freqz(b,a,n,whole),例：已知系统函数 讨论其频率响应特性 程序如figure2.m所示。,2）几何矢量法,其理论依据与上次课介绍的s域的几何矢量法完全类似。这里给出通用函数dplxy（） 例：系统函数为 分析该系统的功能。 程序如ex1203.m所示。,