数字信号实验第四章答案.docx

1、数字信号实验第四章答案数字信号处理实验报告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 lengthM = input(Enter the fil

2、ter length M: );% Compute the frequency samples of the DTFTw = 0:2*pi/1023:2*pi;num = (1/M)*ones(1,M);den = 1;% Compute and plot the DTFTh = freqz(num, den, w);subplot(2,1,1)plot(w/pi,abs(h);gridtitle(Magnitude Spectrum |H(ejomega)|)xlabel(omega /pi);ylabel(Amplitude);subplot(2,1,2)plot(w/pi,angle(h

3、);gridtitle(Phase Spectrum argH(ejomega)xlabel(omega /pi);ylabel(Phase in radians); 所得结果如图示：M=2M=7 幅度和相位谱表现出对称性的类型是由于冲激响应是实数，因此频率响应是周期且对称的，幅度谱是周期甚至对称的，相位响应是周期奇对称。 采用移动平均滤波表示过滤器的类型 低通滤波器 Q2.1的结果可以解释为 它是一个低通滤波器，输入是一个两个正弦分量的总和，一个高频和低频。结果依赖于过滤器的长度，但总的结果是更高频率的正弦输入分量衰减超过较低的频率正弦输入分量。Q4.2 因果LTI离散时间系统的频率响应曲线

4、得到使用修改后的程序如下: 通过此传递函数表示的过滤器类型 带通滤波器Q4.3 对问题的因果LTI离散时间系统的频率响应图q4.3得到使用修改后的程序如下 通过此传递函数表示的过滤器类型- 带通滤波器 4.2和4.3滤波器的不同在于 幅度谱两个是一样的，但是第二个相位谱不是连续的，原因是它们的极点不是一样的，4.36的极点在单位圆内，所以是稳定的，4.37的极点在单位圆外，所以不稳定。 我会选择4.36滤波器的原因是 因为4.37没有4.36稳定Q4.4 过滤器的群延迟的问题q4.4指定和使用功能如下所示: 从图像可以观察到: 这是一个窄阻带的带阻滤波器，在大多数的带通滤波器中，群延迟是恒定的

5、。Q4.5 4.2和4.3两个滤波器的前100个样本冲击响应的图像如图示： 可以得到观察: 4.36给出的滤波器是稳定的，即hn是可求和的，并且冲激响应呈指数型衰减，而4.37给出的是不稳定的，所以hn随n呈指数型增长。Q4.6 使用zplane生成式4.36和4.37确定的两个滤波器的零极点图:从图像可以观察到: 上图的极点都在单位圆内，所以是因果稳定的，下图的极点在单位圆外，不是因果稳定。4.2 传输函数的类型Project 4.2 滤波器 P4_1的代码如下：% Program P4_1% Impulse Response of Truncated Ideal Lowpass Filte

6、rclf;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的有

7、限冲击响应低通滤波器的冲击响应:% Program Q4_8% Impulse Response of Truncated Ideal Lowpass Filterclf; 9 fc = 0.45;n = -9.5:1:9.5;y = 2*fc*sinc(2*fc*n);k = n+9.5;stem(k,y);title(N = 20);axis(0 19 -0.2 0.7);xlabel(Time index n);ylabel(Amplitude);grid; 得到的结果是:Q4.9 必要的修改程序p4_1计算并画出长度与15和0.65的截止频率的FIR低通滤波器的脉冲响应:% Progr

8、am Q4_9% Impulse Response of Truncated Ideal Lowpass Filterclf;fc = 0.65;n = -7.0:1:7.0;y = 2*fc*sinc(2*fc*n);k = n+7.0;stem(k,y);title(N = 14);axis(0 14 -0.4 1.4);xlabel(Time index n);ylabel(Amplitude);grid; 得到的图像是：Q4.10 MATLAB计算程序和绘图的FIR低通滤波器的的幅度响应如下： % Program Q4_10% Compute and plot the amplitud

9、e response% of Truncated Ideal Lowpass Filterclear;% Get N from the user command lineN = input(Enter the filter time shift N: );% compute the magnitude spectrumNo2 = N/2;fc = 0.25;n = -No2:1:No2;y = 2*fc*sinc(2*fc*n);w = 0:pi/511:pi;h = freqz(y, 1, w);plot(w/pi,abs(h);grid;title(strcat(|H(ejomega)|,

10、 N=,num2str(N);xlabel(omega /pi);ylabel(Amplitude); 低通滤波器的幅度相应（若干个n值）:从图像可以得到观察 随着滤波器长度的增加，从通过到不通过变得更加陡峭，我们也可以看到吉布斯现象：当滤波器增加时，幅度相应更加趋向一个理想的低通特征。然而随着w增长，峰值是增加而不是降低。 P4_2的代码:% Program P4_2% Gain Response of a Moving Average Lowpass Filterclf;M = 2;num = ones(1,M)/M;g,w = gain(num,1);plot(w/pi,g);grida

11、xis(0 1 -50 0.5)xlabel(omega /pi);ylabel(Gain in dB);title(M = , num2str(M)回答: Q4.11 长度为2的滑动平均滤波器的增益相应的图像： 从图中可以看出，3-dB截止频率是- pi/2Q4.12 必要的修改程序p4_2计算并画出一个级联的K长度为2的滑动平均滤波器的增益响应如下： % Program Q4_12% Gain Response of a cascade connection of K % two-point Moving Average Lowpass Filtersclear;K = input(Ent

12、er the number of sections K: );Hz = 1;% find the numerator for H(z) = cascade of K sectionsfori=1:K;Hz = conv(Hz,1 1);end;Hz = (0.5)K * Hz;% Convert numerator to dBg,w = gain(Hz,1);% make a horizontal line on the plot at -3 dBThreedB = -3*ones(1,length(g);% make a vertical line on the plot at the% t

13、heoretical 3dB frequencyt1 = 2*acos(0.5)(1/(2*K)*ones(1,512)/pi;t2 = -50:50.5/511:0.5;plot(w/pi,g,w/pi,ThreedB,t1,t2);grid;axis(0 1 -50 0.5)xlabel(omega /pi);ylabel(Gain in dB);title(K = ,num2str(K),; Theoretical omega_c = ,num2str(t1(1); 使用修改后的程序级联部分的增益响应曲线得到如下所示： 从图中可以看出，级联的3-dB截止频率是 -0.30015piQ4.

14、13 必要的修改程序p4_2计算高通滤波器的增益响应（4.42）如下:% Program Q4_13% Gain Response of Highpass Filter (4.42)clear;M = input(Enter the filter length M: );n = 0:M-1;num = (-1).n .* ones(1,M)/M;g,w = gain(num,1);plot(w/pi,g);grid;axis(0 1 -50 0.5);xlabel(omega /pi);ylabel(Gain in dB);title(M = , num2str(M); 通过修改程序，增益响应

15、M=5的曲线为： 我们可以得到3-dB 的截止频率在 大约在0.8196piQ4.14 从式子. (4.16) 3-dB截止频率c 在0.45 我们可以得到 = 0.078702 取代 在式子 (4.15) 和(4.17) 我们得到的一阶IIR低通和高通滤波器的传递函数，分别给出了 HLP(z) = HHP(z) = 我们得到的增益响应如图示:从这些图中我们看到，所设计的滤波器满足规格。幅度响应HLP（Z）HHP（Z）获得使用MATLAB如下： 这两个滤波器是 全通滤波器 HLP(z) 和 HHP(z)一个图的平方幅度响应的总和如图示：这两个滤波器是 - 功率互补Q4.15 从式（4.24），

16、我们得到K = 10, B = 1.8660660 在式子 (4.23)替代B 和 c = 0.3 我们可以得到 = -0.3107925 使用这种 公式值（4.22）我们到达10 IIR低通滤波器级联的传递函数 Substituting c = 0.3 in Eq. (4.16) we obtain = 0.3249197 Using this value of in Eq. (4.15) we arrive at the transfer function of a first-order IIR lowpass filter 增益响应 and 如图示: 我们可以得到观察是 - 第一阶II

17、R低通部分串联实现比单一的第一阶滤波器设计的相同的3 dB截止频率更清晰的过渡带。Q4.16 在式子(4.19)取代 o = 0.61 我们可以得到 -0.3387 在式子(4.20)中 = 0.15 我们可以得到 解得 = 1.6319 和 = 0.6128. . 替代值 和值 在式子 (4.18) 我们得到的的IIR带通的传递函数的传递函数 HBP,1(z) = 取代值 和第二个值 在式子 (4.18) 我们得到的的IIR带通的传递函数的传递函数 HBP,2(z) = 使用 zplane 我们可以发现极点位置 HBP,1(z) 和 HBP,2(z) 我们可以得到稳定的传输函数HBP(z)

18、-稳定的传递函数的增益响应的HBP（Z）的图像如图所示： 使用等值的 和 在式子. (4.21) 接下来我们获得一个稳定的IIR带阻滤波器的传递函数 HBS(z) =传递函数的增益响应HBs（Z）的图像如下所示：从这些图中我们看到，所设计的滤波器不满足规格 幅度相应之和HBP(z) + HBS(z) 的图像如图示:我们从图可以得到两个滤波器是 全通 平方幅度相应之和 HBP(z) and HBS(z) 的图像如图示我们可以得到这两个滤波器是 - 功率互补Q4.17 一个来自原型FIR低通滤波器的梳状滤波器的传递函数方程（4.38）是由 G(z) = H0(zL) = 0.5+0.5z-L为下列

19、值L以上的梳状滤波器的幅度响应图如下所示：从这些图中我们看到，梳状滤波器有缺口 k = (2k+1)pi/L =和 k =2kpi/L, 其中k = 0, 1, . . ., L-1.Q4.18 一个来自原型FIR高通滤波器的梳状滤波器的传递函数式（4.41）M2由下式给出 G(z) = H1(zL) = 0.5+0.5z-L为下列值L以上的梳状滤波器的幅度响应图如下所示我们可以得到梳状滤波器有L个缺口 k = 2kpi/ L 并且 L 的峰值在 k = (2k+1)pi/L.Q4.19 P4_3的代码是% Program P4_3% Zero Locations of Linear Phas

20、e FIR Filtersclf;b = 1 -8.5 30.5 -63;num1 = b 81 fliplr(b);num2 = b 81 81 fliplr(b);num3 = b 0 -fliplr(b);num4 = b 81 -81 -fliplr(b);n1 = 0:length(num1)-1;n2 = 0:length(num2)-1;subplot(2,2,1); stem(n1,num1);xlabel(Time index n);ylabel(Amplitude); grid;title(Type 1 FIR Filter);subplot(2,2,2); stem(n2

21、,num2);xlabel(Time index n);ylabel(Amplitude); grid;title(Type 2 FIR Filter);subplot(2,2,3); stem(n1,num3);xlabel(Time index n);ylabel(Amplitude); grid;title(Type 3 FIR Filter);subplot(2,2,4); stem(n2,num4);xlabel(Time index n);ylabel(Amplitude); grid;title(Type 4 FIR Filter);pausesubplot(2,2,1); zp

22、lane(num1,1);title(Type 1 FIR Filter);subplot(2,2,2); zplane(num2,1);title(Type 2 FIR Filter);subplot(2,2,3); zplane(num3,1);title(Type 3 FIR Filter);subplot(2,2,4); zplane(num4,1);title(Type 4 FIR Filter);disp(Zeros of Type 1 FIR Filter are);disp(roots(num1);disp(Zeros of Type 2 FIR Filter are);dis

23、p(roots(num2);disp(Zeros of Type 3 FIR Filter are);disp(roots(num3);disp(Zeros of Type 4 FIR Filter are);disp(roots(num4); 通过运行程序的p4_3如下产生四的FIR滤波器的脉冲响应图： 我们可以从中得到观察: 滤波器 #1 长度 9并且对称冲击响应并且类型1是线性相位滤波器 。 滤波器 #2 长度 10并且对称冲击响应并且类型 2 是线性相位FIR 滤波器. 滤波器 #3 长度 9 并且 不对称冲击响应并且类型3线性相位FIR 滤波器. 滤波器#4长度10并且不对称冲击响应

24、并且类型 4 线性相位 FIR滤波器 从P4_3的滤波器的零点可以看出: 滤波器#1有零点 z = 2.9744, 2.0888, 0.9790 + 1.4110i, 0.9790 - 1.4110i, 0.3319 + 0.4784i, 0.3319 - 0.4784i, 0.4787, 0.3362 滤波器 #2 有零点 z = 3.7585 + 1.5147i, 3.7585 - 1.5147i, 0.6733 + 2.6623i, 0.6733 - 2.6623i, -1.0000, 0.0893 + 0.3530i, 0.0893 - 0.3530i, 0.2289 + 0.0922

25、i, 0.2289 - 0.0922i 滤波器#3有零点 z = 4.7627, 1.6279 + 3.0565i, 1 .6279 - 3.0565i, -1.0000, 1.0000, 0.1357 + 0.2549i, 0.1357 - 0.2549i, 0.2100 滤波器 #4 有零点 z = 3.4139, 1.6541 + 1.5813i, 1.6541 - 1.5813i, -0.0733 + 0.9973i, -0.0733 - 0.9973i, 1.0000, 0.3159 + 0.3020i, 0.3159 - 0.3020i, 0.2929 对每个这些滤波器的相位响应曲

26、线，利用MATLAB如下所示： 我们可以得到每个滤波器都有线性相位 滤波器# 1群延迟是- 4 滤波器# 2群延迟是 4.5 滤波器# 3群延迟是- 4 滤波器#4群延迟是 4.5Q4.20 通过运行程序的p4_3如下产生四的FIR滤波器的脉冲响应图： 从图像可以得到观察 滤波器#1 是长度 9并且对称冲击响应的类型1线性相位滤波器 滤波器#2 是长度 10并且对称冲击响应的类型2线性相位滤波器 滤波器#3是长度 9并且不对称冲击响应的类型3线性相位滤波器 滤波器#4 是长度10并且不对称冲击响应的类型4线性相位滤波器 通过观察P4_3滤波器的零点得到: 滤波器 #1有零点z = 2.3273

27、 + 2.0140i 2.3273 - 2.0140i -1.2659 + 2.0135i -1.2659 - 2.0135i-0.2238 + 0.3559i -0.2238 - 0.3559i 0.2457 + 0.2126i 0.2457 - 0.2126i 滤波器#2 的零点z = 2.5270 + 2.0392i 2.5270 - 2.0392i -1.0101 + 2.1930i -1.0101 - 2.1930i -1.0000 -0.1733 + 0.3762i -0.1733 - 0.3762i 0.2397 + 0.1934i 0.2397 - 0.1934i 滤波器#3的

28、零点z = -1.0000 0.2602 + 1.2263i 0.2602 - 1.2263i 1.0000 0.6576 + 0.7534i 0.6576 - 0.7534i 0.1655 + 0.7803i 0.1655 - 0.7803i 滤波器#4的零点z = 2.0841 + 2.0565i 2.0841 - 2.0565i -1.5032 + 1.9960i -1.5032 - 1.9960i 1.0000 -0.2408 + 0.3197i -0.2408 - 0.3197i 0.2431 + 0.2399i 0.2431 - 0.2399i 对每个这些滤波器的相位响应曲线，利用

29、MATLAB如下所示： 从这些图像可以得出结论这些滤波器是 线性 相位 滤波器1的群延迟是 - 4 滤波器2的群延迟是 4.5 滤波器3的群延迟是 - 4 滤波器4的群延迟是 4.5Answers:Q4.21 H1(z) 的幅度响应如图示: 有图想知道幅度响应的最大值 = 0.5049pi 并且value = 0.3289 使用zplane H1(z)的极点是在单位圆的外面因此是不稳定的Q4.22 幅度相应 G1(z) 的图像是: 我们可以得到幅度响应的最大值 = 1.7854 G1(z) 的极点在单位圆的里面因此传输函数是稳定的 因为G1(z)幅度相应的最大值= 1.7854 , 我们分割G1(z) 1/1.7854得到一个有界实函数 G2(z) = 4.3 稳定性测试 P4_4的代码：% Program P4_4% Stability Testclf;den = input(Denominator coefficients = );ki = poly2rc(den);disp(Stability test parameters are);disp(ki); 回答:Q4.23 H1(z) 和H2(z) 的零极点图: