数字信号实验第四章答案解析.doc

数字信号处理实验报告4 线性时不变离散时间系统频域分析 实验目的 通过使用matlab做实验来加强对传输函数的类型和频率响应 和稳定性测试来强化理解概念。 4.1 传输函数分析 回答: Q4.1 修改程序P3_1去不同的M值，当0w2pi时计算并画出式（2.13）所示滑动平均滤波器的幅度和相位谱，代码如下： % Program Q4_1 % Frequency response of the causal M-point averager of Eq. (2.13) clear; % User specifies filter length M = input(Enter the filter length M: ); % Compute the frequency samples of the DTFT w = 0:2*pi/1023:2*pi; num = (1/M)*ones(1,M); den = [1]; % Compute and plot the DTFT 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=7 幅度和相位谱表现出对称性的类型是由于–冲激响应是实数，因此频率响应是周期且对称的，幅度谱是周期甚至对称的，相位响应是周期奇对称。 采用移动平均滤波表示过滤器的类型– 低通滤波器 Q2.1的结果可以解释为 – 它是一个低通滤波器，输入是一个两个正弦分量的总和，一个高频和低频。结果依赖于过滤器的长度，但总的结果是更高频率的正弦输入分量衰减超过较低的频率正弦输入分量。 Q4.2 因果LTI离散时间系统的频率响应曲线得到使用修改后的程序如下: – 带通滤波器 Q4.3 对问题的因果LTI离散时间系统的频率响应图q4.3得到使用修改后的程序如下 通过此传递函数表示的过滤器类型- 带通滤波器 4.2和4.3滤波器的不同在于 – 幅度谱两个是一样的，但是第二个相位谱不是连续的，原因是它们的极点不是一样的，4.36的极点在单位圆内，所以是稳定的，4.37的极点在单位圆外，所以不稳定。 我会选择4.36滤波器的原因是 – 因为4.37没有4.36稳定 Q4.4 过滤器的群延迟的问题q4.4指定和使用功能如下所示: : 这是一个窄阻带的带阻滤波器，在大多数的带通滤波器中，群延迟是恒定的。 Q4.5 4.2和4.3两个滤波器的前100个样本冲击响应的图像如图示： 可以得到观察: 4.36给出的滤波器是稳定的，即h[n]是可求和的，并且冲激响应呈指数型衰减，而4.37给出的是不稳定的，所以h[n]随n呈指数型增长。 Q4.6 使用zplane生成式4.36和4.37确定的两个滤波器的零极点图: 从图像可以观察到: 上图的极点都在单位圆内，所以是因果稳定的，下图的极点在单位圆外，不是因果稳定。 4.2 传输函数的类型 Project 4.2 滤波器 P4_1的代码如下： % Program P4_1 % Impulse Response of Truncated Ideal Lowpass Filter clf; fc = 0.25; n = [-6.5:1:6.5]; y = 2*fc*sinc(2*fc*n);k = n+6.5; stem(k,y);title(N = 13);axis([0 13 -0.2 0.6]); xlabel(Time index n);ylabel(Amplitude);grid; 回答: Q4.7 接近于理想低通滤波器的冲击响应的图像： FIR低通滤波器的长度是 - 14 P4_1决定低通滤波器长度的语句是 – n = [-6.5:1:6.5]; 控制截止频率的参数是- fc = 0.25; Q4.8 修改程序P4.1，计算并画出式（4.39）所示长度为20，截止频率为wc=0.45的有限冲击响应低通滤波器的冲击响应: % Program Q4_8 % Impulse Response of Truncated Ideal Lowpass Filter clf; 9 fc = 0.45; n =