巴特沃斯低通滤波器

《数字信号处理》课程设计报告

设计课题滤波器设计与实现

专业班级

姓名

学号

报告日期 2012年12月

《数字信号处理》

课程设计任务书

题目滤波器设计与实现

学生姓名学号专业班级

设计内容与要求一、设计内容：

见所选题目。

二、设计要求

1 设计报告一律按照规定的格式，使用A4纸，格式、封面统一给出模版。

2 报告内容

（1）设计题目及要求

（2）设计原理 (包括滤波器工作原理、涉及到的matlab函数的说明) （3）设计内容（设计思路，设计流程、仿真结果）

（4）设计总结（收获和体会）

（5）参考文献

（6）程序清单

起止时间2012年 12 月 3日至 2011年 12月11 日指导教师签名2011年 12月 2日系（教研室）主任签名年月日

学生签名

年 月 日

《数字信号处理》课程设计报告

一、设计题目及要求

设计题目 基于MATLAB 的巴特沃斯低通滤波器的设计

设计要求

1. 通过实验加深对巴特沃斯低通滤波器基本原理的理解。

2．学习编写巴特沃斯低通滤波器的MATLAB 仿真程序

3. 滤波器的性能指标如下：设计一个模拟低通巴特沃斯滤波器，

技术指标：通带截止频率10000/rad s ，通带最大衰减3dB ；阻带起始频率

30000/rad s

，阻带最小衰减40dB ，画出其幅度谱和相位谱。

二、设计原理

1. 巴特沃斯低通滤波器简介:

巴特沃斯滤波器是电子滤波器的一种，特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。这种滤波器最先由英国工程师斯替芬·巴特沃斯（Stephen Butterworth ）在1930年发表在英国《无线电工程》期刊的一篇论文中提出的，可以构成低通、高通、带通和带阻四种组态，是目前最为流行的一类数字滤波器 ,经过离散化可以作为数字巴特沃思滤波器 ,较模拟滤波器具有精度高、稳定、灵活、不要求阻抗匹配等众多优点 ,因而在自动控制、语音、图像、通信、雷达等众多领域得到了广泛的应用，是一种具有最大平坦幅度响应的低通滤波器。

2.巴特沃斯低通滤波器的设计原理:

巴特沃斯滤波器以巴特沃斯函数来近似滤波器的系统函数。巴特沃斯滤波器是根据幅频特性在通频带内具有最平坦特性定义的滤波器。

巴特沃思滤波器的低通模平方函数表示

()

ΩΩ+

=Ωc N

/22

a 11

)

(j H N=1,2，…… （2-6）

下面归纳了巴特沃斯滤波器的主要特征 a 对所有的N ，()

1a j H 20

=Ω=Ω。

b 对所有的N ，()

707.0a

j 2c

=ΩΩH =

Ω即()dB

3a lg 20j H c

=Ω

=ΩΩ

c()Ωj

H a2是Ω的单调下降函数。

d()Ωj

H a2随着阶次N的增大而更接近于理想低通滤波器。

如下图2所示，可以看出滤波器的幅频特性随着滤波器阶次N的增加而变得越来越好，在截止频率Ωc处的函数值始终为1/2的情况下，通带内有更多的频带区的值接近于1；在阻带内更迅速的趋近于零。

图2 巴特沃思低通滤波平方幅频特性函数

3、系统函数

设巴特沃斯的系统函数为H a（s），则：

4、函数说明

1buttord函数

(1)[N,wc]=buttord(wp，ws，αp，αs)

用于计算巴特沃斯数字滤波器的阶数N和3dB截止频率wc。其中，调用参数wp，ws分别为数字滤波器的通带、阻带截止频率的归一化值，要求：0≤wp≤1，0≤ws≤1(1表示数字频率pi)。当ws≤wp时，为高通滤波器；当wp和ws为二元矢量时，为带通或带阻滤波器，这时wc也是二元向量。

αp，αs分别为通带最大衰减和组带最小衰减（dB）。

N，wc为butter函数的调用参数。

(2)[N,Ωc]=buttord(Ωp，Ωs，αp，αs，‘s’)

用于计算巴特沃斯模拟滤波器的阶数N和3dB截止频率Ωc。其中，Ωp，Ωs，Ωc均为实际模拟角频率。

说明：buttord函数使用阻带指标计算3dB截止频率，这样阻带会刚好满足要求，而通带会有富余。

2 butter函数

(1)[b，a]=butter（N，wc，‘ftype’）

计算N阶巴特沃斯数字滤波器系统函数分子、分母多项式的系数向量b、a。其中，调用参数N和wc分别为巴特沃斯数字滤波器的阶数和3dB截止频率的归一化值（关于pi归一化），一般是调用buttord(1)格式计算N和wc。系数b、a是按照z-1的升幂排列。

(2)[B，A]=butter（N，Ωc，‘ftype’，‘s’）

计算巴特沃斯模拟滤波器系统函数的分子、分母多项式系数向量ba、aa。其中，调用参数N和Ωc分别为巴特沃斯模拟滤波器的阶数和3dB截止频率(实际角频率)，可调用buttord(2)格式计算N和Ωc。系数B、A按s的正降幂排列。 tfype为滤波器的类型：

○1ftype=high时，高通；Ωc只有1个值。

○2ftype=stop时，带阻阻；此时Ωc=[Ωcl,Ωcu]，分别为带阻滤波器的通带3dB下截止频率和上截止频率。

○3ftype缺省时，若Ωc只有1个值，则默认为低通；若Ωc有2个值，则默认为带通；其通带频率区间Ωcl < Ω < Ωcu。

所设计的带通和带阻滤波器系统函数是2N阶。因为带通滤波器相当于N阶低通滤波器与N阶高通滤波器级联。

三、设计内容

1.设计思路

巴特沃斯低通滤波器的设计步骤大致为：