matlab的窗长和形状对滤波器性能,实验六 用窗函数设计FIR滤波器(附思考题程序)...

实验六 用窗函数设计FIR 滤波器

1.实验目的

(1) 熟悉FIR 滤波器设计的方法和原理

(2) 掌握用窗函数法设计FIR 滤波器的方法和原理，熟悉滤波器的特性 (3) 了解各种窗函数滤波器特性的影响

2.实验原理

FIR 滤波器的设计方法主要有三种：窗函数法、频率取样法、切比雪夫等波纹逼近法。FIR 滤波器的设计是要寻求一系统函数)(z H ，使其频率响应)(ωj e H 逼近滤波器要求的理想频率响应()j d H e ω

，其对应的单位脉冲响应)(n h d 。 (1)用窗函数设计FIR 滤波器的基本方法

在时域用一个窗函数截取理想的)(n h d 得到)(n h ，以有限长序列)(n h 近似逼近理想的

)(n h d ；在频域用理想的)(ωj d e H 在单位圆上等角度取样得到h(k)，根据h(k)得到H(z)将

逼近理想的Hd(z)。

设理想滤波器)(ωj d e H 的单位脉冲响应为)(n h d 。以低通线性相位FIR 数字滤波器为例。

?∑-

-∞

-∞

==

=π

π

ωωω

ω

ω

π

d e e H n h e n h e H jn j d d jn n d

j d )(21)()()(

)(n h d 一般是无限长的、非因果的，不能直接作为FIR 滤波器的单位脉冲响应。要想得

到一个因果的有限长的滤波器h(n)，最直接的方法是截断)()()(n w n h n h d =，即截取为有限长因果序列，并用合适的窗函数进行加权作为FIR 滤波器的单位脉冲响应。按照线性相位滤波器的要求，h(n)必须是偶对称的。对称中心必须等于滤波器的延时常数，即

??

?-==2

/)1()

()()(N a n w n h n h d 用矩形窗设计的FIR 低通滤波器，所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象，且最大波纹大约为幅度的9%，(现象称为吉布斯(Gibbs )效应)。

(2)典型的窗函数

(a )矩形窗(Rectangle Window)

)()(n R n w N =