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

《数字信号处理》实验报告五 系别： 计系 班级： 通信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; %避免被零除 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。程序代码和仿真结果如下: clf; Fs=16000; [n,fo,mo,w1]=firpmord([1500 2500],[1 0],[0.06 0.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。程序代码和仿真结果如下: clf; Fs=16000; [n,fo,mo,w1]=firpmord([3500 4500],[1 0],[0.01 0.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，采样频率为16000HZ。程序代码和仿真结果如下: clf; Fs=16000; [n,fo,mo,w1]=firpmord([1500 2500 3500 4500], [1 0], [0.06 0.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.1，阻带纹波的最大允许值为0.01，采样频率为16000HZ。程序代码和仿真结果如下: clf; Fs=16000; [n,fo,mo,w1]=firpmord([1500 2500 3500 4500], [1 0], [0.06 0.1], 16000); b=firpm(n,fo,mo,w1); [h,w]=freqz(b); plot(w*Fs /(2*pi),abs(h)); xlabel(频率/Hz); ylabel(幅度); 三、调试和运行程序过程中产生的问题及采取的措施： MATLAB程序运行出现错误后，切换到MATLAB命令行中，观察命令行中的错误信息，确定错误原因，在错误信息中，查找源程序中错误代码的行数，点击行数链接，可以直接定位错误语句，错误语句定位后，自动跳转到代码编辑器中，在代码编辑器中通过光标标识错误代码的位置