实验6_离散时间系统的z域分析报告.docxVIP

实验6离散时间系统的z域分析 、实验目的 掌握z变换及其反变换的定义，并掌握 MATLAB实现方法。 学习和掌握离散时间系统系统函数的定义及 z域分析方法。 掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 、实验原理 Z变换 序列x(n)的z变换定义为 X(z)x( n)z X(z) x( n)z Z反变换定乂为 x(n)? x(n) ?X(z)zn 1dz ■ 在MATLAB 中，可以采用符号数学工具箱的 ztrans函数和iztrans 函数计算z变换 和z反变换： Z=ztrans(F) 求符号表达式F的z变换。 F=ilaplace(Z) 求符号表达式Z的z反变换。 离散时间系统的系统函数 离散时间系统的系统函数 H(z)定义为单位抽样响应 h(n)的z变换 H(z) h(n)z z变换之比得到 z变换之比得到 此外，连续时间系统的系统函数还可以由系统输入和输出信号的 H(z) 丫(z)/X(z) 由上式描述的离散时间系统的系统函数可以表示为 H(z)bo dz1 …bMZ M H(z) 1 N a° a〔z …a” z 离散时间系统的零极点分析 离散时间系统的零点和极点分别指使系统函数分子多项式和分母多项式为零的点。在 MATLAB中可以通过函数 roots来求系统函数分子多项式和分母多项式的根，从而得到系 统的零极点。 此外，还可以利用 MATLAB的zplane函数来求解和绘制离散系统的零极点分布图， zplane函数调用格式为： zpla ne(b,a) b,a 为系统函数的分子、分母多项式的系数向量(行向量) 。 zpla ne(z,p)乙p 为零极点序列(列向量)。 系统函数是描述系统的重要物理量，研究系统函数的零极点分布不仅可以了解系统单位抽样响应的变化，还可以了解系统的频率特性响应以及判断系统的稳定性： 系统函数的极点位置决定了系统单位抽样响应 h(n)的波形，系统函数零点位置只影 响冲激响应的幅度和相位，不影响波形。 系统的频率响应取决于系统的零极点， 根据系统的零极点分布情况，可以通过向量分 析系统的频率响应。 因果的离散时间系统稳定的充要条件是 H(z)的全部极点都位于单位圆内。 、实验内容 (1)已知因果离散时间系统的系统函数分别为: ① H(z)3 2z 0.5z 0.005z 0.3② H(z)3z4 3z ① H(z) 3 2 z 0.5z 0.005z 0.3 ② H(z) 3z4 3z3 z2 3z 1 864202 4 0n m 山 Q J-Q.d上 864202 4 0n m 山 Q J-Q.d 上 mLL_rra匚 一-JEE- 恳 * ■■■O 0.5198 - 0.5346i -0.5396 + 0.0000i □ Real Pg* I- 0.0000 + 0.0000i [] 实验结果分析：由零极点分布可得 冲激响应： h(n)=((-1.5272 - 2.2795*i)*(0.5198 + 0.5346i)An+(-1.5272 0.5346*i)A n+(-0.2790)*(-0.5396)A n)*heaviside( n) 频率响应： + 2.2795*i)*(0.5198 试采用MATLAB 画出其零极点分布图，求解系统的冲激响应 h(n)和频率响应 H( e )， 并判断系统是否稳定。 MATLAB 代码如下： b=[1 2 1]; a=[1 -0.5 -0.005 0.3]; zplan e(b,a); b1=[1 2 1]; a1=[1 -0.5 -0.005 0.3 0]; [r,p,k]=residue(b1,a1) r = -1.5272 - 2.2795i -1.5272 + 2.2795i -0.2790 + 0.0000i 3.3333 + 0.0000i 0.5198 + 0.5346i H(ejw)(ejw)2 2ejw 1(ejw)3 0.5(ejw)2 0.005ejw 0.3由于该系统所有极点位于②MATLAB H(ejw) (ejw)2 2ejw 1 (ejw)3 0.5(ejw)2 0.005ejw 0.3 由于该系统所有极点位于 ②MATLAB 代码如下： Z平面单位圆内，故系统是稳定的。 b=[1 -1 0 2]; a=[3 3 -1 3 -1]; zplan e(b,a); b1=[1 -1 0 2]; a仁[3 3 -1 3 -1 0]; [r,p,k]=residue(b1,a1) -0.1375 + 0.0000i 0.2628 + 0.3222i 0.2628 - 0.3222i 1.6119 + 0.0000i -2.0000 + 0.0000i P = -1.6462 + 0.0000i