matlab 滤波窗函数,FIR滤波器窗函数设计法详细步骤以及Matlab代码

采用窗函数法设计理想低通，高通滤波器，参考北京交通大学陈后金主编的【数字信号处理】5.2节 窗函数法设计线性相位FIR数字滤波器P164，和P188。

设计步骤如下：

1) 确定滤波器类型，不同的FIR类型可设计不同类型的滤波器，I型可设计LP(低通滤波器)，HP(高通滤波器)，BP(带通滤波器)，BS(带阻滤波器)。

Fir I型

Fir II型

Fir III型

Fir IV型

LP，HP，BP，BS

LP，BP

BP

HP，BP，BS

2) 确定设计的滤波器的参数

Eg：若要设计一个低通滤波器，fp=20，fs=30；Ap=1，As=40，则3db截频Wc = 2*pi*(fs-fp)/Fs；Fs为采样频率。

当选定某一窗函数时，衰耗Ap和As就已经确定，凯撒窗除外。Ap和As的计算方法可参看另外一篇博客： https://www.cnblogs.com/xhslovecx/p/10118570.html

3) 确定窗函数

窗的类型

主瓣宽度

近似过渡带宽度

δp，δs

Ap(dB)

As(dB)

矩形窗

4pi/N

1.8pi/N

0.09

0.82

21

Hann

8pi/N

6.2pi/N

0.0064

0.056

44

Hamming

8pi/N

7pi/N

0.0022

0.019

53

Blackman

12pi/N

11.4pi/N

0.0002

0.0017

74

Kaiser

可调窗，需要确定 β值

50

21<=A<=50,   β=0.5842(A-21)^0.4 + 0.07886(A-21)；

A<21,             β = 0；

4) 确定滤波器的阶数M，首先确定滤波器的长度N。对于除凯撒窗以外的窗函数&#x