数字信号处理第四章分析.pdf

第四章 线性时不变离散时间系统的频域分析 一、传输函数和频率响应 例 4.1 传输函数分析 Q4.1 clear; M = input(Enter the filter length M: ); w = 0:2*pi/1023:2*pi; num = (1/M)*ones(1,M); den = [1]; h = freqz(num, den, w); subplot(2,1,1) plot(w/pi,abs(h));grid title(Magnitude Spectrum |H(e^{j\omega})|) xlabel(\omega /\pi); ylabel(Amplitude); subplot(2,1,2) plot(w/pi,angle(h));grid title(Phase Spectrum arg[H(e^{j\omega})]) xlabel(\omega /\pi); ylabel(Phase in radians); M=2 M=10 M=15 幅度谱为偶对称，相位谱为奇对称，这是一个低通滤波器。 M 越大，通带越窄且过渡带越 陡峭。 Q4.2 使用修改后的程序 P3.1 ，计算并画出当 w=[0,pi] 时传输函数 的 因果线性时不变离散时间系统的频率响应。它表示哪种类型的滤波器？ w = 0:pi/511:pi; num = [0.15 0 -0.15]; den = [1 -0.5 0.7]; 如下图 1 这是一个带通滤波器。 图 1 图 2 Q4.3 对下面的传输函数重做习题 Q4.2 ： ，式（ 4.36 ）和式（ 4.37 ） 给出的两个滤波器之间的区别是什么？你将选择哪一个滤波器来滤波，为什么？ w = 0:pi/511:pi; num = [0.15 0 -0.15]; den = [0.7 -0.5 1]; 如上图 2 也是一个带通滤波器， 这两个滤波器的幅度谱是一样的， 相位谱不太一样， 我会选 择第一个带通滤波器，因为它的相位谱更加平滑，相位失真小。 Q4.4 使用 MATLAB计算并画出当 w=[0,pi] 时因果线性时不变离散时间系统的群延迟。系统 的传输函数为 。 clf; w = 0:pi/511:pi; num = [1 -1.2 1]; den = [1 -1.3 1.04 -0.222]; h= grpdelay(num,den,w); plot(w/pi,h); xlabel(w/pi); ylabel( 群延迟 ); Q4.5 使用 Q3.50 中编写的程序，分别计算并画出式（ 4.36 ）和式（ 4.37 ）确定的两个滤波 器的冲激响应中的前一百个样本。讨论你的结果。 clf; num = [0.15 0 -0.15]; den = [0.7 -0.5 1]; L = input( 输入样本数 L: ); [g t] = impz(num,den,L); stem(t,g); title([ 前 ,num2str(L), 脉冲响应的样本 ]); xlabel( 时间序号 n); ylabel(h[n]); （4.36 ）式 （4.37 ）式 由图可知：这些情节由 impz 给生成的因果的脉冲响应实