实验四z变换(通信11用).pptVIP

2.实验原理与方法 函数RESIDUEZ 功能：部分分式的展开成如下的形式： 调用格式：[R,P,K] = RESIDUEZ(B,A) ％说明：由矢量a和b构成的数字滤波器的差分方程为 其中： B=[b0,b1,…,bM,bM-1],A=[a0,a1,…,aN,aN-1]，要求A的首项不为0 ％R-以列矢量形式出现的对应极点的留数R=[r(1) r(2) r(3) r(4)] P-以列矢量形式出现的极点矢量P=[P(1) P(2) P(3) P(4)] K-以行矢量出现的直接项多项式K=[k(1) k(2)] 例 b = [1,0.4*sqrt(2)]; a=[1,-0.8*sqrt(2),0.64]; [R,p,k] = residuez(b,a) Mp=abs(p) % pole magnitudes Ap=angle(p)/pi % pole angles in pi units  运行结果： R = 0.5000 - 1.0000i 0.5000 + 1.0000i p = 0.5657 + 0.5657i 0.5657 - 0.5657i k = [] Mp = 0.8000 0.8000 Ap = -0.2500 0.2500 最小相位系统 定义----零极点都在z平面单位圆内的因果系统称为 最小相位系统。记为Hmin(z)。 M1-6 实验四信号与系统的Z域分析 1.实验目的 （1）实现信号的Z变换 (2)了解Z变换的应用： 信号、系统的z域表示 部分分式法的应用 (3) 零点、极点和增益常数表示 (4) 二阶因子表示 [z,p,k]=tf2zp(b,a) 将z的有理函数表示转换为零点、极点和增益常数表示。 [b,a]=zp2tf(z,p,k) 将零点、极点和增益常数表示转换为有理函数表示。 sos=zp2sos(z,p,k) (second-order section) 将零点、极点和增益常数表示转换为2阶因子表示。 利用MATLAB求解不同表示形式的系统函数 例： 试求下面系统函数的零极点形式二阶因子形式。 %Determination of the factored form and %the second order section form of a % rational z-transform b =[1 0 0.04 0]; a =[1 -0.8 0.16 -0.128]; [z,p,k]=tf2zp(b,a); disp(Zeros are at); disp(z); disp(Poles are at); disp(p); disp(Gain constant);disp(k); sos=zp2sos(z,p,k); disp(Second-order sections); disp(real(sos)); 程序的运行结果为 Zeros are at 0 0 + 0.2000i 0 - 0.2000i Poles are at 0.8000 0.0000 + 0.4000i 0.0000 - 0.4000i Gain constant 1 例： 试求下面系统函数的零极点形式二阶因子形式。 程序的运行结果为 Second-order sections 0.2000 0 0 1.0000 -0.8000 0 5.0000 0 0.2000 1.0000 0.0000 0.1600 例： 试求下面系统函数的零极点形式二阶因子形式。 系统函数与系统稳定性 LTI系统稳定的充要条件： H(z)的收敛域ROC包含单位圆 系统函数H(z)的极点位于Z平面单位圆内时,系统稳定 对于因果LTI系统，系统稳定的充要条件可以通过系统函数的极点位置来判断，即 利用MATLAB画出z平面的零极点分布 zplane(b,a) b，a ：正幂表示的H(z)分子和分母多项式系数矩阵 %zeros and poles of the transfer function b =[1 2 1 0]; a =[1 -0.5 -0.005 0.3]; zplane(b,a); 例： 已知一离散因果的LTI系统的系统函数为 画出z平面的零极点分布图。 -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 Real part Imagi