您的当前位置：首页实验七用窗口法设计FIR数字滤波器

实验七用窗口法设计FIR数字滤波器

来源：欧得旅游网

电信类课程试验报告

学院：基础信息工程

系别：电子信息工程

课程名称：数字信号处理

姓名：

学号：

日期：

实验七

实验名称：用窗口法设计FIR 数字滤波器

一、实验目
了解一个实际滤波器设计过程，加深掌握用窗口法设计FIR 数字滤波器原理与窗函数对数字滤波器性能影响。

二、主要函数简介

在MATLAB 数字信号处理工具箱中提供了函数fir1.fir1 是采用经典窗口法设计线性相位FIR 数字滤波器,且

具有标准低通,带通,高通与带阻等类型,调用格式为B=fir1(n, w n ,’ftype,windows)

0～1;对于带通,带通滤波器,
输入参数:n 为FIR 数字滤波器阶数,对于高通,带阻滤波器n 取偶数.

w =[ w , w ] ,且w < w 。’ftype’为滤波器类型,缺省时为地通或带通滤波器, 为归一化滤波器截止频率；,取值范围

输出参数:B 为FIR 数字滤波器系数向量,长度为n+1. 为’high’时是高通滤波器,为’stop’时为带阻滤波器.window 为窗函数,其长度为n+1,缺省时,自动取海明(hamming)窗.

三、实验程序

实验内容：

1．编写用改进升余弦窗设计FIR 线性相应低通数字滤波器程序,已知

c



0.5,N=21.这里

c

为低通滤波器截止

频率,N 为滤波器长度,滤波器阶数为N-1.

2．调试运行程序,要求在屏幕上显示出单位脉冲响应h(n)数值.画出其幅度响应 |H

e

j



|及 20log

e

j



曲线.

3．画出窗口函数

(n

)

及其频谱W

e

j

幅度|W

e

j

|与20log

e

j

|曲线.

实验程序：

首先新建一个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;

第 1 页

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);
实验图形：

第2页

第3页

四、实验小结

五、教师评语
教师签字：年月日

希望以上资料对你有所帮助，附励志名言十条：
1、上帝说：你要什么便取什么，但是要付出相当的代价。

2、目标的坚定是性格中最必要的力量源泉之一，也是成功的利器之一。没有它，天才会在矛盾无定的迷径中徒劳无功。

3、当你无法从一楼蹦到三楼时，不要忘记走楼梯。要记住伟大的成功往往不是一蹴而就的，必须学会分解你的目标，逐步实施。

4、无论你觉得自己多么的了不起，也永远有人比你更强。

5、没有播种，何来收获；没有辛苦，何来成功；没有磨难，何来荣耀；没有挫折，何来辉煌。

6、相信就是强大，怀疑只会抑制能力，而信仰就是力量。

7、一个人除非自己有信心，否则带给别人信心。

8、除非想成为一流，否则就是二流。 9、在你内心深处，还有无穷的潜力，有一天当你回首看时，你就会知道这绝对是真的。

10、成功，往往住在失败的隔壁！

第4页