matlab中fir1中的边界频率,实验四 参考 基于MATLAB的FIR数字滤波器设计

实验4 基于MATLAB 的FIR 数字滤波器设计

一、 实验目的

(1) 加深对数字滤波器的常用指标理解。 (2) 学习数字滤波器的设计方法。 二、 实验原理 低通滤波器

1()1()P P P S

S H H δδδπ?-≤Ω≤+Ω≤Ω??

Ω≤Ω≤Ω≤??

低通滤波器的常用指标：

1、 通带边缘频率P Ω，

2、 阻带边缘频率S

Ω ，

3、 通带起伏P δ，

4、

通

带

峰

值起伏

]

)[1(log

2010

dB p p

δα

--=，

5、

阻带起伏

s

δ，最小阻带衰减

]

)[(log 2010dB s S δα-=。

三、 要求 (1)在MA TLAB 中，熟悉函数fir1、kaiserord 、remezord 、remez 的使用；

B = fir1(n,Wn,'high','noscale')设计滤波器；

[n,Wn,beta,ftype] = kaiserord(f,a,dev)估计滤波器阶数；

[n,fo,ao,w] = remezord (f,a,dev ,fs)计算等波纹滤波器阶数n 和加权函数w(ω)； B=remez(n,f,a)进行等波纹滤波器的设计

(2)阅读附录中的实例，学习FIR 滤波器的设计方法及其在MA TLAB 中的实现；

(3)给出FIR 数字滤波器的冲激响应，绘出它们的幅度和相位频响曲线，讨论它们各自的实现形式和特点。 四、 实验内容

利用MA TLAB 编程，分别用窗函数法和等波纹滤波器法设计两种FIR 数字滤波器，指标要求如下：

通带边缘频率：π

45.01=ΩP ，

π

65.02=ΩP ，通带峰值起伏：

]

[1dB p

≤α

。

阻带边缘频率：π3.01=ΩS ，

π

75.02=ΩS ，最小阻带衰减：

]

[40dB S ≥α。

附录：

例1 用凯塞窗设计一FIR 低通滤波器，通带边界频率π

3.0=Ωp ，阻带边界频率

π

5.0=Ωs ，阻带衰减 不小于50dB 。

S Ω-P

Ω-P S

ΩPassband Stopband

Transition

band

Fig 1 Typical magnitude

specification for a digital LPF