1、电信类课程试验报告学 院:基础信息工程 系 别:电子信息工程 课程名称:数字信号处理姓 名: 学 号: 日 期:实验七 实验名称:用窗口法设计 FIR 数字滤波器一、实验目的了解一个实际滤波器设计过程,加深掌握用窗口法设计 FIR 数字滤波器的原理和窗函数对数字滤波器性能的影响。二、主要函数简介在 MATLAB 的数字信号处理工具箱中提供了函数 fir1.fir1 是采用经典窗口法设计线性相位 FIR 数字滤波器,且具有标准低通,带通,高通和带阻等类型,调用格式为 B=fir1(n, ,ftype,windows) nw输入参数:n 为 FIR 数字滤波器的阶数,对于高通,带阻滤波器 n 取偶
2、数. 为归一化的滤波器截止频率;,取n值范围 01;对于带通,带通滤波器 , = , ,且 。 ftype为滤波器类型,缺省时为地通或带通滤波nw1212器,为high时是高通滤波器,为stop时为带阻滤波器.window 为窗函数,其长度为 n+1,缺省时,自动取海明(hamming)窗.输出参数:B 为 FIR 数字滤波器系数向量,长度为 n+1.三、实验程序实验内容:1编写用改进的升余弦窗设计 FIR 线性相应低通数字滤波器的程序,已知 0.5,N=21.这里 为低通滤波器cc的截止频率,N 为滤波器的长度 ,滤波器的阶数为 N-1.2调试运行程序,要求在屏幕上显示出单位脉冲响应 h(n
3、)的数值.画出其幅度响应 |H |及 20log |H | je10je的曲线.3画出窗口函数 及其频谱 W 的幅度|W |和 20log |H |的曲线.)(njeje10je实验程序:首先新建一个 M 文件,然后把下面一段代码复制进去 M 文件,保存为 ideal_lp.m。在该目录下可以调用 ideal_lp 函数。function hd=ideal_lp(wc,M);alpha = (M-1)/2;n = 0:1:(M-1);m = n - alpha +eps; hd = sin(wc*m)./(pi*m);运行程序:clear all;wc=0.5*pi;N=21;hd =ideal_lp(wc,N);n=0:1:N-1;stem(n,hd);figure;freqz(hd,1,512);w_box=hamming(N);h_box=fir1(N-1,wc/pi,w_box);figure;freqz(h_box,1,512);figure;stem(n,w_box);实验图形:四、实验小结五、教师评语教师签字: 年 月 日