半带FIR低通滤波器

数字信号处理课程设计

设计题目半带FIR低通滤波器设计

题目编号 0205 学院名称电气学院

指导教师

班级

学号

学生姓名

目录

一．半带滤波器的概述 (1)

二．半带FIR滤波器的性质及设计方法 (1)

2.1半带FIR滤波器的性质 (1)

2.2半带滤波器的设计方法 (2)

2.3滤波器指标的确定 (3)

2.4手工计算 (4)

三．在MATLAB中仿真 (4)

3.1滤波器的系数 (4)

3.2半带FIR滤波器的频率特性曲线 (4)

3.3基于MATLAB的幅频响应曲线 (5)

四．滤波器的结构不同对性能指标的影响 (6)

4.1利用直接型结构构建滤波器 (6)

4.2利用级联型结构构建数字滤波器 (7)

五．参数字长对性能指标的影响 (8)

5.1参数字长取2位对性能指标的影响 (8)

5.2参数字长取8位对性能指标的影响 (8)

5.3参数字长取12位对性能指标的影响 (9)

5.4参数字长取14位对性能指标的影响 (9)

5.5结论 (9)

六．心得体会 (9)

七．参考文献 (10)

八．附录 (10)

1

半带FIR 低通滤波器的设计

一．半带滤波器的概述

在数字滤波器的设计过程中，为了能够有效地进行抽取滤波，往往采用多级抽取的方法，这就需要采用内插和抽取的原理。如果对滤波器进行M 倍抽取，则输出的采样频率是输入采样频率的I ／M ；如果对滤波器进行M 倍内插，则输出的采样频率是输入采样频率的M 倍。

半带滤波器是一种基于抽取和内插原理的滤波器，它是一种特殊的低通FIR(有限冲激响应)数字滤波器——一种抽取因子为2的抽取滤波器，这种滤波器由于通带和阻带相对于二分之一Nyquist 频率对称，因而有近一半的滤波器系数为O 。由于系数为0的部分在运算的过程中不需要消耗运算量，所以运算量减少了一半，使实际滤波过程中的运算量大幅度减少。因此计算的效率高，实时性比较强，有利于滤波运算的实时实现，多速率信号处理中有着特别重要的地位，并且得到了广泛的应用。

二．半带FIR 滤波器的性质及设计方法

2.1半带FIR 滤波器的性质

图1为抽取因子为2的半带滤波器的框图，其中:)(z H h 为抽取滤波器，s f 为输入采样2/1s s f f =为输出采样率。

图1 抽取因子为2的半带滤波器框图

现在考虑频率特性如图2所示的一种特殊的FIR 滤波器，即为半带FIR 滤波器，它具有如下的特性：

（1）通带纹波p δ和阻带纹波s δ相等，即s p δδ=；

（2）通带边频p F 和阻带边频s F 相对于4/s f 对称，即2/s s p f F F =+，用数字频率表示为：πωω=+s p 。

f s

x(n)

H h (z)

f s

2

2/1s s f f =

2 2.2半带滤波器的设计方法

设计滤波器的过程中，根据滤波器的冲击响应)(n h 和节数N ，有以下4种类型的滤波器：

（1）类型I :)(n h 偶对称，N 为偶数，4种滤波器都可以设计。

（2）类型II ：)(n h 偶对称，N 为奇数，可以用来设计低通和带通滤波器，不可以用来设计高通和带阻滤波器。

（3）类型III ：)(n h 奇对称，N 为偶数，只能用来设计带通滤波器，其他类型的滤波器都不可以设计。

（4）类型IV :)(n h 奇对称，N 为奇数，可以用来设计高通和带通滤波器，不可以用来设计低通和带阻滤波器。

根据图2可以看出半带FIR 滤波器是一种特殊的低通滤波器，并且其相位要求是线性的，这里采用类型II 进行设计，即：

)1()(n N h n h --=

式中N 为滤波器的节数，且为奇数。

对于这种N 为奇数的偶对称FIR 滤波器，其频率响应H （e jw )可以表示为： H h (e jw )=H h (ω)e -jw(N-1)/2 式中：s f f /2πω=，为相对数字频率；

∑-==

2

/)1(0

)cos()()(N n h n n a H ωω (1)

由线性相位FIR 滤波器性质可知其系数首先要满足：

)2

1

()0(-=N h a )2

1

(2)(n N h n a --= n=1,2,...,21-N

其次，由半带滤波器的定义πωω=+s p ，s p σσ=和图2可以得到： )()(ωπω-+h h H H =1 (2) 将（1）式带到（2）式中得：

∑-=-+2

/)1(0

)]cos())[cos((N n n n n n a ωπω

=

∑-=2

/)1(0

))[cos((N n n n a ω+（-1）n

)]

cos(n ω

3

=

∑-2

/)1()cos()(2N n n n a 为偶数项

ω

=1

故有：5.0)2

1

(

)0(=-=N h a 0]22/)1[()2(=--=n N h n a n=1,2,...,4

1

-N 上式表明：半带滤波器除了2

1

-=

N n 点外，所有a(n)的偶次系数均为0，滤波器系数集的形式为：)(n h ={x,0,x,0,...,x,0.5,x,0,...,0,x}

设计半带FIR 数字滤波器有窗函数法和频率采样法2种方法，这里采用窗函数法来设计半带FIR 滤波器。

图2给出了半带滤波器的设计参数)(p p F ω和)(s s F ω以及s p δδ=。由图2可知：在用窗函数法设计半带滤波器时，其理想幅频特性可以表示为：

H h (e jw )=｛

可求得：

⎰-=2/2/21)(πππe n h d -jaw e jwn ωd =

)

()(2sin a n a n -⎥⎦⎤⎢⎣⎡-ππ,其中，21-=N a 。 其次是选择窗函数，为了满足半带滤波器对p F ,s F 及s p δδ=的要求，这里利用Kaiser 窗进行设计。

2.3滤波器指标的确定

（1）计算过度带：p s ωωω-=∆ （2）确定滤波器阶数N ：1285.28

lg 20+∆--≥

ω

δN ， N 取奇数

（3）计算凯塞窗：)

())112(

1()(020ββωI N n

I n k ---= ，10-≤≤N n 式中：

e -jaw 2

πω≤

π

ωπ

≤≤2