matlab设计长度为50的滤波器,实验5 基于Matlab的数字滤波器设计

实验五 基于Matlab 的数字滤波器设计

实验目的：加深对数字滤波器的常用指标和设计过程的理解。

实验原理：低通滤波器的常用指标如下所述，其典型规格如下：

??

????????≤≤≤≤+≤≤-πδδδw w e G w w e G s s jw p p jw p ,)(,1)(1 通带边缘频率：p ?

阻带边缘频率：s ? 通带起伏：p δ 通带峰值起伏：))(1(log 2010dB p p δα--=

阻带起伏： s δ 最小阻带衰减：))((log 2010dB s s δα-=

数字滤波器有IIR 和FIR 两种类型，他们的特点和设计方法不同。

一、窗函数法设FIR 滤波器

在Matlab 中产生窗函数十分简单：

(１) 矩形窗

调用格式：w=boxcar(n) ，根据长度n 产生一个矩形窗w 。

(２) 三角窗

调用格式：w=triang(n) ，根据长度n 产生一个三角窗w 。

(3)汉宁窗

调用格式：w=hanning(n) ，根据长度n 产生一个汉宁窗w 。

(4)汉明窗

调用格式：w=hamming(n) ，根据长度n 产生一个汉明窗w 。

(5)布莱克曼窗

调用格式：w=Blackman(n) ，根据长度n 产生一个布莱克曼窗w 。

基于窗函数的FIR 滤波器设计利用Matlab 提供的函数firl 来实现。

调用格式：firl(n,wn,’ftype ’,Window)，n 为阶数,wn 是截止频率(如果输入是形如[w1 w2]的矢量时，本函数将设计带通滤波器，其通带为w11、设计一个长度为8的线性相位FIR 滤波器。

Window=boxcar(8); b=fir1(7,0.4,Window); freqz(b,1); Window=blackman(80;

b=fir1(7,0.4,Window); freqz(b,1);

2、设计线性相位带通滤波器，其长度N=15，上下边带截止频率分别为w1=0.3π,w2=0.5π Window=blackman(16); B=firl(15,[0.3 0.5],Window); Freqz(b,1);

3、设计指标为，dB A dB R s a p p 50,3.0,25.0,2.0====π?π?

的底通滤波器。

见课本P197。

二、IIR 滤波器的实现

(1)脉冲响应不变法设计数字ButterWorth 滤波器

调用格式：[bz,az]=impinvar(b,a,s F )，再给定模拟滤波器参数b,a 和取样频率s F 的前提下，计算数字滤波