实验5FIR滤波器的设计.doc

《实验5FIR滤波器的设计.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、《数字信号处理》实验报告五系别：计系班级：通信13002学号：姓名王泓媛日期：2015.4.23指导教师：黄泽鑫一、实验内容：FIR低通滤波器的设计，FIR带阻滤波器的设计。学会使用系统提供的库函数firpmord、firpm、freqs，并理解其输入和输出的各个参数的含义。理解滤波器各项技术指标的含义。二、调试程序，并完成下面习题1.矩形窗低通FIR滤波器的设计，阶数为13。N=________13_____;a=(N-1)/2;Wc=pi/2;n=[0:1:(N-1)];m=n-a+eps;%避

2、免被零除hd=sin(Wc*m)./(pi*m);[_______H1______,W]=freqz(___hd,1__________);plot(W/pi,20*log10(abs(H1)/abs(H1)));xlabel('频率');ylabel('幅频响应');title('FIR加矩形窗时的频谱图');2.FIR低通数字滤波器:截止频率为1500HZ，阻带的起始频率为2500HZ,通带纹波的最大允许值为0.06，阻带纹波的最大允许值为0.1，采样频率为16000HZ。程序代码和仿真结果如下

3、:clf;Fs=16000;[n,fo,mo,w1]=firpmord([15002500],[10],[0.060.1],16000);b=firpm(n,fo,mo,w1);[h,w]=freqz(b);plot(w*Fs/(2*pi),abs(h));xlabel('频率/Hz');ylabel('幅度');3.FIR高通数字滤波器:导通频率为4500HZ，阻带的起始频率为3500HZ,通带纹波的最大允许值为0.1，阻带纹波的最大允许值为0.01，采样频率为16000HZ。程序代码和仿真结果如

4、下:clf;Fs=16000;[n,fo,mo,w1]=firpmord([35004500],[10],[0.010.1],16000);b=firpm(n,fo,mo,w1);[h,w]=freqz(b);plot(w*Fs/(2*pi),abs(h));xlabel('频率/Hz');ylabel('幅度');4.FIR带通数字滤波器:通带的带宽为2500HZ~3500HZ，阻带的起始频率分别为1500HZ和4500HZ，通带纹波的最大允许值为0.1，阻带纹波的最大允许值为0.01，采样频率为

5、16000HZ。程序代码和仿真结果如下:clf;Fs=16000;[n,fo,mo,w1]=firpmord([1500250035004500],[10],[0.060.1],16000);b=firpm(n,fo,mo,w1);[h,w]=freqz(b);plot(w*Fs/(2*pi),abs(h));xlabel('频率/Hz');ylabel('幅度');5.FIR带阻数字滤波器:阻带频率为2500HZ~3500HZ，通带起始频率分别为1500HZ和4500HZ，通带纹波的最大允许值为0

6、.1，阻带纹波的最大允许值为0.01，采样频率为16000HZ。程序代码和仿真结果如下:clf;Fs=16000;[n,fo,mo,w1]=firpmord([1500250035004500],[10],[0.060.1],16000);b=firpm(n,fo,mo,w1);[h,w]=freqz(b);plot(w*Fs/(2*pi),abs(h));xlabel('频率/Hz');ylabel('幅度');三、调试和运行程序过程中产生的问题及采取的措施：MATLAB程序运行出现错误后，切换到

7、MATLAB命令行中，观察命令行中的错误信息，确定错误原因，在错误信息中，查找源程序中错误代码的行数，点击行数链接，可以直接定位错误语句，错误语句定位后，自动跳转到代码编辑器中，在代码编辑器中通过光标标识错误代码的位置