用窗函数法设计 FIR 数字滤波器
一、实验目的
1.掌握用窗函数法设计 FIR 数字滤波器的原理和方法。
2.熟悉线性相位 FIR 数字滤波器特征。
3.了解各种窗函数对滤波特性的影响。
二、实验仪器
微型计算机 matlab 软件
三、实验原理和方法
如果所希望的滤波器的理想频率响应函数为
H
d
(e
h
d
(n)
=
j
)
,则其对应的单位脉冲响应为
1
2
H
d
(e
j
)e
j
d
(2-1)
窗函数设计法的基本原理是用有限长单位脉冲响应序列
h(n)
逼近
h
d
(n)
。由于
h
d
(n)
往往是无限长序列,且是非因果的,所以用窗函数
(n)
将
h
d
(n)
截断,并进行加权处理,
得到:
h(n)
=
h
d
(n)
(n)
(2-2)
j
其频率响应函数
H
d
(e )
为:
h(n)
就作为实际设计的FIR 数字滤波器的单位脉冲响应序列,
H
d
(e
j
)
=
h(n)e
j
(2-3)
n0
N 1
式中,N 为所选窗函数
(n)
的长度。
由第七章可知,用窗函数法设计的滤波器性能取决于窗函数
(n)
的类型及窗口长度 N
的取值。设计过程中,要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和
窗口长度 N。各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见第七章。
这样选定窗函数类型和长度 N 后,求出单位脉冲响应
h(n)
=
h
d
(n)
·
(n)
,并按式
(2-3)求出
H (e
果
H (e
j
j
)
。
H (e
j
)
是否满足要求,要进行验算。一般在
h(n)
尾部加零使长度
j
满足于 2 的整数次幂,以便用 FFT 计算
H (e )
。如果要观察细节,补零点数增多即可。如
)
不满足要求,则要重新选择窗函数类型和长度 N,再次验算,直至满足要求。
如果要求线性相位特性,则
h(n)
还必须满足
h(n) h(N 1 n)
(2-4)
根据上式中的正负号和长度 N 的奇偶性又将线性相位 FIR 滤波器分成四类。要根据设计的滤
波特性正确选择其中一类。例如,要设计线性低通特征,可选择
h(n) h(N 1 n)
一类,
而不能选
h(n) h(N 1 n)
一类。
四、实验内容
1.复习用窗函数法设计 FIR 数字滤波器一节内容,阅读本实验原理掌握设计步骤。
2.编写程序
① 编写能产生矩型窗、哈明窗、汉宁窗、莱克曼窗的窗函数子程序。