双线变换法数字低通滤波器

双线变换法数字低通滤波器

燕山大学

课程设计说明书

题目：双线性变换法设计数字低通滤波器

学院（系）：电气工程学院

年级专业： 10级检测2班

学号： 100103020122

学生姓名：刘培露

指导教师：王娜

教师职称：讲师

课程名称：数字信号处理课程设计

基层教学单位：仪器科学与工程系指导教师：王娜

21

说明：1、此表一式四份，系、指导教师、学生各一份，报送院教务科一份。

2、学生那份任务书要求装订到课程设计报告前面。

电气工程学院教务科

目录

21

一.课题描述 (2)

二.设计原理 (2)

2.1 IIR数字滤波器设计原理 (2)

2.2巴特沃斯低通滤波器的原理 (3)

2.3双线性变换法 (4)

2.4用双线法设计巴特沃斯数字低通滤波器的步

21

骤 (9)

三. MATLAB程序 (10)

四.程序中命令介绍 (12)

五. 图像结果……………………………………………………………………………. ..14

六. 结果分析 (16)

七.心得体会: (16)

21

一.课题描述

数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。可以设计系统的频率响应，让它满足一定的要求，从而对通过该系统的信号的某些特定的频率成分进行过滤，这就是滤波器的基本原理。如果系统是一个连续系统，则滤波器称为模拟滤波器。如果系统是一个离散系统，则滤波器称为数字滤波器。

数字滤波实质上是一种运算过程，实现对信号的运算处理。输入数字信号（数字序列）通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为是一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性:

()()()ω

ωj

ω

j e

j

e

X

Y=

e

H

其中()ωj e Y、()ωj e X分别是数字滤波器的输出序列和输入序列的频域特性（或称为频谱特性）,()ωj e H是数字滤波器的单位取样响应的频谱，又称为数字滤波器的频域响应。输入序列的频谱()ωj e X经过滤波后()ωj e X()ωj e H,因此，只要按照输入信号频谱的特点和处理信号的目的，适当选择()ωj e H，使得滤波后的()ωj e X()ωj e H满足设计的要求，这就是数字滤波器的滤波原理。

二.设计原理

21

2.1 IIR数字滤波器设计原理

IIR数字滤波器的设计一般是利用目前已经很成熟的模拟滤波器的设计方法来进行设计，通常采用模拟滤波器原型有butterworth函数、chebyshev 函数、bessel函数、椭圆滤波器函数等。

IIR数字滤波器的设计步骤：

① 按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标；

② 根据模拟滤波器技术指标设计为响应的模拟低通滤波器；

③ 跟据脉冲响应不变法和双线性不变法把模拟滤波器转换为数字滤波器；

④ 如果要设计的滤波器是高通、带通或带阻滤波器，则首先把它们的技术指标转化为模拟低通滤波器的技术指标，设计为数字低通滤波器，最后通过频率转换的方法来得到所要的滤波器。

在MATLAB中，经典法设计IIR数字滤波器主要采用以下步骤：

图2.1 IIR数字滤波器设计步骤

21

21

2.2巴特沃斯低通滤波器的原理

巴特沃斯滤波器的特点是同频带内的频率响应曲线最为平坦，没有起伏，而在组频带则逐渐下降为零。在振幅的对数对角频率的波特图上，从某一边界见频率开始，振幅随着角频率的增加而逐渐减少，趋向于负无穷大。

一阶巴特沃斯滤波器的衰减率为每倍频20分贝，二阶巴特沃斯滤波器的衰减率为每倍频12分贝，三阶的衰减率为每分贝18分贝，如此类推，巴特沃斯滤波器的振幅对角频率单调下降，并且滤波器的结束越高，在组频带振幅衰减速度越快，其他滤波器高阶的振幅对角频率图和低阶数的振幅对角频率有不同的形状。

N c

s s H s H )(

11)()(2

2Ω-+=

-

上述函数的特点是等距离分布在半径为Ω的圆上。

因此，极点用下式表示为

N

k j j c k e

e s )12(2+∏

Ω=

1,2,1,0-=N k

)

(s H a 的表示式：

∏-=-Ω=

1

)

()(N k k

n c

a s

s s H

为了使设计公式和图表统一，将频率归一化。巴特沃斯滤波器采用3dB 截止频率c Ω归一化，归一化后的系统函数为

∏-=Ω-Ω=Ω1

0)(1

)(

N k c

k c c

a s s s G

令c c s j p ΩΩ=Ω=+=λλη,,λ称为归一化频率，p 称为归一化复变

21

量，这样巴特沃斯滤波器的归一化低通原型系统函数为

∏-=-=

10

)

(1

N k k

a p

p G

式中，c k s p Ω=，为归一化极点，用下式表示:

)

21221(N

k j k

e

p

++=π 1,2,1,0-=N k

2.3双线性变换法

1.双线性变换法的基本原理

由于从s 平面到z 平面的映射sT

z e =具有多值性，使得设计出来的数字滤波器不可避免的出现频谱混迭现象。为了克服脉冲响应不变法可能产生的频谱混叠效应的缺点，我们使用一种新的变换——双线性变换。双线性变换法可认为是基于对微分方程的积分，利用对积分的数值逼近的思想 。 仿真滤波器的传递函数()H s 为

1

(),M

k

k k N

k

k k c s

H s M N d s

===

>∑∑

将展开为部份分式的形式，并假设无重复几点，则

1()N

k

k pk

A H s s s ==-∑