4-信号与系统实验.pdfVIP

第四次实验 实验名称：z 变换及离散系统的z 域分析 实验目的： （1）学会运用Matlab 求离散信号的Z 变换和反变换； （2 ）会会运用Matlab 分析离散信号系统函数的零极点； （3 ）学会运用Matlab 分析系统函数零极点分布与其时域特性的关系； （4 ）运用Matlab 进行离散时间系统的频率特性分析； 实验原理： (1)离散信号的Z 变换和反变换 ∞ X (z ) Z [x (n)] ∑x (n)z −n , n −∞ x (n) Z −1[X (z )]; 在Matlab 中分别使用：Z=ztrans(x); X=iztrans(Z); (2)系统函数的零极点分布 Y z B z ( ) ( ) H (z ) ，B(z)=0 时的根即为零点，A(z)=0 时的根即为极点； X z A z ( ) ( ) Matlab 中的命令为：[Z,P,K]=tf2zp(B,A); (3)系统函数零极点分布与其时域特性的关系: 与拉氏变换在连续系统中的作用类似，在离散系统中，z 变换建立了时域函数h(n)与Z 域函 数H(z)之间的关系，因此，Z 变换H(z)可以在形式上反应h(n) 的特点。在此，可以通过Matlab 作图更加直观地反映。 (4)离散时间系统的频率特性分析: 在稳定的因果离散系统中，令z= ejw ,则有系统的频率响应为： H (ejw ) | H (ejw ) | ej ϕ( w) ； 其中，| H (ejw ) | 表示幅频特性，ej ϕ( w) 表示相频特性。 在Matlab 中可以使用：[H,w]=freqz(B,A,N) [H,w]=freqz(B,A,N,’whole’) 实验内容： 1. 使用计算机仿真实验样例的例题 例4-1,试用ztrans 函数求下列函数的z 变换。 (1)x(n)=ancos(Пn)u(n); x=sym(a^n*cos(pi*n)); Z=ztrans(x); simplify(Z) ans = z/(a + z) n-1 n-1 (2)x(n)=[2 -(-2) ]u(n) x=sym(2^(n-1)-(-2)^(n-1)); Z=ztrans(x); simplify(Z) ans = 4/(z^2 - 4) + 1 例4-2 ，试用iztrans 函数求下列函数的z 反变换。 2 (1) X(z)=(8z-19)/(z -5z+6); Z=sym((8*z-19)/(z^2-5*z+6)); x=iztrans(Z); x x = (3*2^n)/2 + (5*3^n)/3 - (19*kroneckerDelta(n, 0))/6 2 3 (2)X(z)=z(2z -11z+12)/(z-1)(z-2) Z=sym(z*(2*z^2-11*z+12)/(z-1)/(z-2)^3); x=iztrans(Z) x = (5*2^n)/2 - (2^n*binomial(n - 1, 2))/2 - 2^n*(n - 1) - 3 例4-3 ，试用Matlab 命令对函数X(z)=18/(18+3z-1-4z-2-z-3)进行部分分式展开，并求出其z 反变换。 B=[18]; A=[18,3,-4,-1]; [R,P,K]=residuez(B,A) R = 0.3600 0.2400 0.4000 P = 0.5000 -0.3333 -0.3333 K = [] 例4-4 ，已知一离散因果LTI 系统的系统函数为H(z)=(z+0.32)/(z2+z+0.16),试用matlab 命令 求该系统的零极点。