数字信号处理实验四 - 图文

实验四 线性时不变离散时间系统的频域分析

实验室名称： 信息学院2204 实验时间：2015年11月5日

姓 名：蒋逸恒 学号：20131120038 专业： 通信工程 指导教师：陶大鹏

成绩 教师签名： 年 月 日 一、实验目的 1、对线性时不变离散频域系统进行matlab分析，总结离散时不变系统的特点。 2、掌握matlab对离散是不变系统的应用，熟悉数字滤波器的设计过程，掌握其原理。 二、实验内容 Q4.1修改程序P3.1，取三个不同的M值，当w属于[0 2π]时计算并画出 式（2.13）所示滑动平均滤波器的幅度和相位谱。证明由幅度和相位谱表现出的对称类型。它表示了哪种类型的滤波器？你现在能解释习题Q2.1的结果吗？ Q4.2使用修改后的程序P3.1，计算并画出当w属于[0 π]时的传输函数 H（z)=0.15（1-z-2）/（1-0.5z-1+0.7z-2)（4.36) 的因果线性时不变离散时间系统的频率响应。它表示那种类型的滤波器？ Q4.3对下面的传输函数重做习题 H（z)=0.15（1-z-2）/（0.7-0.5z-1+0.7z-2)(4.37) 式（4.36）和式（4.37）给出的两个滤波器之间的区别是什么？你将选择哪一个滤波器来滤波，为什么？ Q4.6 使用zplane分别生成式（4.36）和式（4.37）确定的两个滤波器的极点零图。讨论你的结果。 Q4.8 修改程序P4.1计算并画出式（4.39）所示长度为20，截止角频率为wt=0.45的有限冲激响应低通滤波器的冲激响应。 Q4.9修改程序P4.1计算并画出式（4.39）所示长度为15，截止角频率为wt=0.65的有限冲激响应低通滤波器的冲激响应。 Q4.10编写一个MATLAB程序，计算并画出式（4.39）所示有限冲激应低通滤波器的振幅响应。使用这个程序，选取几个不同的N值，画出振幅响应并讨论你的结果。 Q4.11 运行程序P4.2，计算并画出一个长度为2的滑动平均滤波器的增益响应。从图中验证3dB截止在π/2处。 Q4.23用MATLAB产生如下两个因果系统传输函数的极零点图： H1（z)=1/（1-1.848z-1+0.85z-2) H2（z)=1/（1-1.851z-1+0.85z-2) 研究生成的极零点图，你可以判断它们的稳定性吗？ 三、实验器材及软件

1. 微型计算机1台 2. MATLAB 7.0软件 四、实验原理 4.1；4.2；4.3 传递函数是指零初始条件下线性系统响应（即输出）量的拉普拉斯变换（或z变换）与激励（即输入）量的拉普拉斯变换之比。传输函数决定了一个系统的性能和用处，他的频率响应特性将直接影响该系统对输入的信号的作用。传递函数的傅立叶逆变换实际上就是该系统的冲激响应，所以冲激响应又称为系统的系统响应。 4.6 传递函数的零极点图反映了该系统的稳定性，通过理论学习，我们知道当使用离散傅立叶变换时（z域变换），必须保证该系统函数的零极点在单位圆之内才会证明该系统是一个稳定的系统，否则系统是一个不稳定的系统。 4.7；4.8；4.10 滤波器是对输入信号进行处理的一种器件，分为数字滤波和模拟滤波，数字滤波是使用软件的方式对离散化的数据进行滤波，而模拟滤波是使用硬件电路实现对连续信号的滤波。滤波器的系统函数可表现出滤波器的类型，包括低通、高通、带通、带阻四种，而所有的滤波器都可以通过修改低通滤波器来获得相应的幅频特性，所以系统滤波器是所有滤波器的基础。理想的滤波器幅频特性的波形是矩形窗，在时域的波形将是一个无限长的波形，这是不可能实现的，所以我们不可能做到理想的情况，只可以尽可能的接近优良状态。 4.23 系统函数的零极点图将直接反映系统的稳定性，只有稳定的系统才能被应用起来，在S域（复频域）的傅立叶变换中，我们要求极点必须在虚轴的左侧才是稳定的系统，而在Z域（离散域）中，我们要求极点必须在单位圆之内才能保证该系统函数反映的系统是一个稳定的系统。 五、实验步骤 1、进行本实验，首先必须熟悉matlab的运用，所以第一步是学会使用matlab。 2、学习相关基础知识，根据《数字信号处理》课程的学习理解实验内容和目的。 3、在充分熟悉基础知识的情况下进行实验，利用matlab完成各种简单的波形产生和观察，理解各种波形产生的原理和方法。 4、从产生的图形中学习新的知识，掌握实验的目的，充分学习数字信号处理的运用。 5、最后思考各种波形的联系和建立完整的知识体系。 六、实验记录（数据、图表、波形、程序等） 4.1 2|H(ej?)|幅度谱振幅M = input('滤波器所需的长度= '); M=2 1.5num = ones(1,M);w = 0:8*pi/511:2*pi; h = freqz(num, 1, w);subplot(2,1,1) plot(w/pi,abs(h));grid title('|H(e^{j\\omega})|幅度谱') xlabel('\\omega /\\pi'); 以弧度为单位的相位10.5000.20.40.60.811.21.41.61.82? /?相位谱arg[H(ej?)]210-1-200.20.40.60.811.21.41.61.82ylabel('振幅'); subplot(2,1,2) plot(w/pi,angle(h));grid title('相位谱arg[H(e^{j\\omega})]') xlabel('\\omega /\\pi'); ylabel('以弧度为单位的相位');

? /? M=4 M=8 |H(ej?)|幅度谱|H(ej?)|幅度谱8 46 34 22 1000.20.40.60.811.2 000.20.40.60.811.21.41.61.82? /?? /? 相位谱arg[H(ej?)]相位谱arg[H(ej?)]44 22 00 -2-2 -4-400.20.40.60.811.21.41.61.8200.20.40.60.811.2? /?? /? 4.2 num = [0.15 0 -0.15]; |H(ej?)|幅度谱den = [1 -0.5 0.7] 1w = 0:8*pi/511:pi; h = freqz(num, den, w); 0.5subplot(2,1,1) plot(w/pi,abs(h));grid 0振幅振幅1.41.61.82以弧度为单位的相位以弧度为单位的相位1.41.61.82振幅title('|H(e^{j\\omega})|幅度谱') xlabel('\\omega /\\pi'); ylabel('振幅'); subplot(2,1,2) plot(w/pi,unwrap(angle(h)));grid title('相位谱arg[H(e^{j\\omega})]') xlabel('\\omega /\\pi'); ylabel('以弧度为单位的相位'); 以弧度为单位的相位00.10.20.30.40.5? /?0.60.70.80.91相位谱arg[H(ej?)]210-1-200.10.20.30.40.50.60.70.80.91? /? 4.3 num = [0.15 0 -0.15]; den = [0.7 -0.5 1] w = 0:8*pi/511:pi; 振幅1|H(ej?)|幅度谱h = freqz(num, den, w); subplot(2,1,1) plot(w/pi,abs(h));grid title('|H(e^{j\\omega})|幅度谱') xlabel('\\omega /\\pi'); ylabel('振幅'); subplot(2,1,2) plot(w/pi,unwrap(angle(h)));grid title('相位谱arg[H(e^{j\\omega})]') xlabel('\\omega /\\pi'); ylabel('以弧度为单位的相位'); 0.5000.10.20.30.40.50.60.70.80.91? /?相位谱arg[H(ej?)]6420以弧度为单位的相位00.10.20.30.40.50.60.70.80.91? /?

4.6 num = [0.15 0 -0.15]; Imaginary Part10.50-0.5-1-3-2-1式4.36den = [1 -0.5 0.7]; subplot(2,1,1); zplane(num,den); title('式4.36'); den1 = [0.7 -0.5 1]; subplot(2,1,2); zplane(num,den1); title('式4.37'); 0Real Part式4.371231Imaginary Part0.50-0.5-1-3-2-10Real Part123 4.7 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('时间序号n'); ylabel('振幅');grid; N = 130.60.50.40.3振幅0.20.10-0.1-0.2 4.8 clf; fc = 0.45/(2*pi); n = [-9.5:1:9.5]; y = 2*fc*sinc(2*fc*n); k = n+9.5; stem(k,y);title('N = 19'); axis([0 19 -0.2 0.6]); xlabel('时间序号n'); ylabel('振幅');grid; 02468时间序号n1012N = 190.60.50.40.3振幅0.20.10-0.1-0.20246810时间序号n12141618

4.10 wc = 0.25*2*pi; N=9; n = [-N/2:1:N/2]; y1=fft(y); k = 2*(n+N/2)/N; plot(k,abs(y1)); title('N = 9'); xlabel('w/pi'); ylabel('振幅');grid; 振幅N = 91.41.21y = sin(wc.*n)./(pi.*n); 0.80.60.40.2000.20.40.60.81w/pi1.21.41.61.82 1.4 1.2 1 0.8 0.6 0.4 0.2 00 4.11 振幅N = 991.4N = 9991.210.8振幅0.60.40.2000.20.40.60.81w/pi1.21.41.61.820.20.40.60.81w/pi1.21.41.61.82M = 20-5-10-15单位为dB的增益X: 0.502Y: -3.037M = 2; num = ones(1,M)/M; [g,w] = gain(num,1); plot(w/pi,g);grid axis([0 1 -50 0.5]) xlabel('\\omega /\\pi'); ylabel('单位为dB的增益'); title(['M = ', num2str(M)]) -20-25-30-35-40-45-5000.10.20.30.40.50.60.70.80.91Imaginary Part 4.23 den = [1 -1.848 0.85]; den1 = [1 -1.851 0.85]; subplot(2,1,1); zplane(1,den); title('H1(z)'); subplot(2,1,2); zplane(1,den1); title('H2(z)'); ? /?H1(z)10.50-0.5-1-3-2-10Real PartH2(z)12321Imaginary Part0.50-0.5-1-3-2-10Real Part1232

本文来源：https://www.bwwdw.com/article/fzyt.html