ch1_5离散系统复频域分析.pptVIP

数字信号处理(Digital Signal Processing) * * 信号与系统系列课程组 国家电工电子教学基地 离散信号与系统分析基础 离散信号与系统的时域分析 离散信号的频域分析 离散系统的频域分析 双边z变换与反变换 离散系统的系统函数 全通滤波器与最小相位系统 信号的抽样与重建 离散系统z域分析 离散系统函数H(z) 差分方程和系统函数H(z) 系统函数H(z)的表示方式 利用MATLAB求解不同表示形式的系统函数 系统函数与系统稳定性 简单数字滤波器 系统函数H(z) 对于离散LTI系统: y[k]=x [k]*h[k] Y(z)=H(z)X(z) H(z)称为离散LTI系统的系统函数 当H(z) 的ROC包含单位圆时 H(z)=Y(z)/X(z)=Z{h[k]} 差分方程和系统函数 N=0, a0?0 时，系统称FIR (Finite Impulse Response); N0，{ak ;k=1,2,...,N}中至少有一项非零时，系统被称 为IIR (Infinite Impulse Response)系统。 系统函数H(z)的表示方式 (1) z-1的有理函数表示 (2) z的有理函数表示 系统函数H(z)的表示方式 (3) 零点、极点和增益常数表示 (4) 二阶因子表示 [z,p,k]=tf2zp(b,a) 将z的有理函数表示转换为零点、极点和增益常数表示。 [b,a]=zp2tf(z,p,k) 将零点、极点和增益常数表示转换为有理函数表示。 sos=zp2sos(z,p,k) (second-order section) 将零点、极点和增益常数表示转换为二阶因子表示。 利用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 例： 试求下面系统函数的零极点形式二阶因子形式。 利用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 Imaginary part 2 例： 已知一离散因果的LTI系统的系统函数为 画出z平面的零极点分布图。 系统函数与系统稳定性 LTI系统稳定的充要条件： H(z)的收敛域ROC包含单位圆 因果系统H(z)的极点位于z平面单位圆内时,系统稳定 对于因果LTI系统，系统稳定的充要条件可以通过系统函数的极点位置来判断，即 例： 已知一离散LTI系统的系统函数为 (1) |z|3 系统不稳定、因果 (2) 2|z