实验十一z变换及离散时间系统z域分析分析解析.docxVIP

PAGE PAGE 63 南昌大学实验报告 学生姓名： 周倩文 学 号： 6301712010 班级： 通信 121 班 实验类型： ■验证□综合□设计□创新 实验日期： 5 月 30 号 实验成绩： z 变换及离散时间系统的 Z 域分析 一、目的 掌握利用 MATLAB 绘制系统零极点图的方法 掌握离散时间系统的零极点分析方法 掌握用 MATALB实现离散系统频率特性分析的方法 掌握逆 Z 变换概念及 MATLAB实现方法 二、离散系统零极点 线性时不变离散系统可用线性常系数差分方程描述，即 ?N a y(n ? i) ? ?M b x(n ? j) （8-1） i j i?0 j ?0 其中 y(k ) 为系统的输出序列， x(k ) 为输入序列。 将式（8-1）两边进行 Z 变换的 Y (z) ?M b z? j j  B(z) H (z) ? ? j ?0 ? （8-2） 将式（8-2）因式分解后有： X (z) ?N i?0 a z?i i A(z) ?M (z ? q ) j H (z) ? C j ?1 （8-3） ?N i?1 (z ? p ) i 其中C 为常数， q 极点。 ( j ? 1,2, j , M ) 为 H (z) 的 M 个零点， p (i ? 1,2, , N ) 为 H (z) 的 N 个 i 系统函数 H (z) 的零极点分布完全决定了系统的特性，若某系统函数的零极点已知， 则系统函数便可确定下来。 因此，系统函数的零极点分布对离散系统特性的分析具有非常重要意义。通过对系统函数零极点的分析，可以分析离散系统以下几个方面的特性： 系统单位样值响应h(n) 的时域特性； 离散系统的稳定性； 离散系统的频率特性； 三、离散系统零极点图及零极点分析 零极点图的绘制 设离散系统的系统函数为 H (z) ? B(z) A(z) 则系统的零极点可用 MATLAB的多项式求根函数 roots()来实现，调用格式为： p=roots(A) 其中 A 为待根求多项式的系数构成的行矩阵，返回向量 p 则是包含多项式所有根的列向量。多项式根的 MATLAB命令举例如下： A=[1 3/4 1/8]; P=roots(A) 运行结果为： P = -0.5000 -0.2500 需注意的是，在求系统函数零极点时，系统函数可能有两种形式：一种是分子、分母多项式均按 z 的降幂次序排列；另一种是分子、分母多项式均按 z?1 的升幂次序排列。这两种方式在构造多项式系数向量时稍有不同。 H (z) 按 z 的降幂次序排列：系数向量一定要由多项式最高次幂开始，一直到常数 项，缺项要用 0 补齐；如 H (z) ? z3 ? 2z z4 ? 3z3 ? 2z2 ? 2z ?1 其分子、分母多项式系数向量分别为 A=[1 0 2 0]、B=[1 3 2 2 1]。 H (z) 按 z?1 的升幂次序排列：分子和分母多项式系数向量的维数一定要相同，不足的要用 0 补齐，否则 z ? 0 的零点或极点就可能被漏掉。如 1? 2z?1 H (z) ? 1 1 1? z?1 ? z?2 2 4 其分子、分母多项式系数向量分别为 A=[1 2 0]、B=[1 1/2 1/4]。 用 roots()求得 H (z) 的零极点后，就可以用plot()函数绘制出系统的零极点图。下面是求系统零极点，并绘制其零极点图的MATLAB实用函数 ljdt()，同时还绘制出了单位圆。 例 1：绘制如下系统函数的零极点 （1） H (z) ? 3z3 ? 5z2 ?10z z3 ? 3z2 ? 7z ? 5 （2） H (z) ? 1? 0.5z?1 1? 3 4 z?1 ? 1 8 z?2 解：MATLAB命令如下 （1） 绘制的零极点图如图 8-1（a）所示。 （2） 绘制的零极点图如图 8-1（b）所示。 （a） (b) 离散系统零极点分析 离散系统零极点分布与系统稳定性 《信号与系统》课程已讲到离散系统稳定的条件为： 时域条件：离散系统稳定的充要条件为 对可和； ?? h(n) ? ? ，即系统单位样值响应绝 n??? Z 域条件：离散系统稳定的充要条件为系统函数 H (z) 的所有极点均位于 Z 平面的单位圆内。 对于三阶以下的低阶系统，可以利用求根公式求出系统函数的极点，从而判断系统 的稳定性，但对于高阶系统，手工求解则显得十分困难，这时可以利用 MATLAB 来实现。实现方法是调用前述的函数 ljdt()绘出系统的零极点图，然后根据极点的位置判断系统的稳定性。 例 2：系统函数如例 1 所示，判断两个系统的稳定性。 解：由例 1 绘出的零极点图可以看出两个系统的稳定性分别