发布时间 : 星期一 文章云南师范大学-dsp（数字信号处理）试验报告10更新完毕开始阅读1a390f5502020740be1e9b6a
m=[1 1 1 1 0.589 0 0 0]; b2=fir2(15,f,m); [H2,w]=freqz(b2,1,128); legend('Ideal','fir2 Designed'); subplot(2,2,3); plot(f,m,w/pi,abs(H2)); title('Comparison of Frequency Response Magnitudes'); subplot(2,2,4); H_db=20*log10(abs(H2)); plot(w,H_db); grid on 运行波形如下显示： Comparison of Frequency Response Magnitudes1.550Ideal10 0.5-50000.51-100 01234Comparison of Frequency Response Magnitudes1.550100.5-50000.51-10001234 8
（3） 利用频率抽样法设计一带通FIR数字滤波器，并绘制其衰耗特性，已知阶数M=15，给定指标为 0 m=0,1,7 Ad[m]= 0.456 m=2,6 1 m=3,4,5 改变Ad[2]或Ad[6]的值，观察该FIR低通数字滤波器的衰耗特性的变化。 【解】 f=[0 1/pi 1/pi 1/pi 1/pi 1/pi 1/pi 1]; m=[0 0 0.456 1 1 1 0.456 0]; b=fir2(15,f,m); [H,w]=freqz(b,1,128); legend('Ideal','fir2 Designed'); subplot(2,2,1); plot(f,m,w/pi,abs(H)); title('Comparison of Frequency Response Magnitudes'); subplot(2,2,2); H_db=20*log10(abs(H)); plot(w,H_db); grid on m=[0 0 0.223 1 1 1 0.456 0]; b2=fir2(15,f,m); [H2,w]=freqz(b2,1,128); legend('Ideal','fir2 Designed'); subplot(2,2,3); plot(f,m,w/pi,abs(H2)); title('Comparison of Frequency Response Magnitudes'); subplot(2,2,4); H_db=20*log10(abs(H2)); plot(w,H_db); grid on 运行波形如下显示 Comparison of Frequency Response Magnitudes10Ideal-200.5-40-60000.51-80 01234 Comparison of Frequency Response Magnitudes10-200.5-40-60000.51-8001234 9
2．对实验现象、实验结果的分析及其结论 通过该试验加深理解FIR数字滤波器的时域特性和频域特性，掌握FIR数字滤波器的设计原理和设计方法，以及FIR数字滤波器的应用。 在FIR窗口法设计中，为何采用不同特性的窗函数？选用窗函数的依据是什么？ 答：FIR滤波器的通带和阻带波动一致，均由窗函数的旁瓣面积确定。因此，采用窗函数设计FIR滤波器，其通带衰减Ap和阻带衰减As不能分别单独控制，而是依赖于所选择的窗函数。在工程中应用一些旁瓣幅度较小的窗函数，来提高FIR滤波器的阻带衰减。 经验：学习使用matlab快捷键很重要，能大大提高编程效率。比如方向键的上键能快速输入上次的内容；ctrl+r键全部注释， ctrl+t键取消全部注释（选中的部分）；另外，每次编写之前都写一句：clc,clear,close all,清除工作空间及面板，这样不会受上一次的操作影响等，使用matlab的help能查找到不熟悉的函数使用方法。 另外，截图时采用Figure中的Edit->Copy Figure以减少图片所占容量，便于编辑。还可以通过File->Save As…保存不同格式的figure图片，便于其他地方其他格式时使用。 教师评语及评分： 签名： 年 月 日
10