ch1_6系统函数教案解析.pptVIP

FIR一阶低通数字滤波器产生的梳状数字滤波器，M=5 试由HLP1(z)获得一梳状滤波器。 解： 例：选取二阶阻带滤波器参数a=0.8, b= cos(0.4p)， M=4时，梳状数字滤波器的频幅响应 设计一梳状滤波器。 解： 信号卷积与解卷 N=500; k=0:N-1; x=cos(0.02*pi*k)+cos(0.5*pi*k); h=[1/4 1/4 1/4 1/4]; y=conv(x,h); % Display input signal x); w=linspace(-1.0*pi,pi,500); subplot(3,2,1); plot(x); grid on; title(Input signal x=cos(0.02*pi*k)+cos(0.5*pi*k)); % Display input signal frequency spectrum); X=freqz(x,[1],w); subplot(3,2,2); plot(w/pi,abs(X)); grid on; title(Input signal frequency spectrum X(ejW)); 信号卷积与解卷 % Display output signal); subplot(3,2,3); plot(y); title(output signal y=x*h); Y=freqz(y,[1],w); subplot(3,2,4); plot(w/pi,abs(Y)); title(Output signal frequency spectrum Y(ejW)); % Deconvolution and display the result); h0=deconv(y,x) x0=deconv(y,h); subplot(3,2,5); plot(x0-x); grid on; title(Deconvolution x0=deconv(y,h), Error x0-x); % Display system frequency response); H=freqz(h,[1],w); subplot(3,2,6); plot(w/pi,abs(H)); grid on; title(System amplitude frequency response |H(ejW)|); legend(h=[1/4,1/4,1/4,1/4],Location,NorthWest); 信号卷积与解卷 数字信号处理(Digital Signal Processing) 信号与系统系列课程组 国家电工电子教学基地 离散信号与系统分析基础 离散信号与系统的时域分析 离散信号的频域分析 离散系统的频域分析 双边z变换 系统函数 全通滤波器与最小相位系统 信号的抽样与重建 系统函数 离散系统z域分析 离散系统函数H(z) 差分方程和系统函数H(z) 系统函数H(z)的表示方式 利用MATLAB求解不同表示形式的系统函数 简单数字滤波器 系统函数 系统函数H(z) 对于离散LTI系统: H(z)=Z {h[k]}=Y(z)/X(z) H(z)称为离散LTI系统的系统函数 当H(z) 的ROC包含单位圆时 当h[k]是实序列时，由H(ejW)的对称性质可得 系统函数 差分方程和系统函数 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) 将零点、极点和增益常数表示转换为2阶因子表示。 利用MATLAB求解不同表示形式的系统函数 系统函数 利用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