基于MATLAB的巴特沃斯低通滤波器的设计.docVIP

姓名： 班级： 学号： 时间：2011年6月 设计题目 基于MATLAB的巴特沃斯低通滤波器的设计 设计要求 通过实验加深对巴特沃斯低通滤波器基本原理的理解。 2．学习编写巴特沃斯低通滤波器的MATLAB仿真程序 3. 滤波器的性能指标如下：通带截止频率fp=5kHz，通带最大衰减=2dB，阻带截止频率fs=12kHz，阻带最小衰减=30dB 设计过程 设计原理 1.1 巴特沃斯低通滤波器简介: 巴特沃斯滤波器是 \o 电子滤波器 电子滤波器的一种，特点是 \o 通频带（尚未撰写） 通频带内的 \o 频率响应 频率响应曲线最大限度平坦，没有起伏，而在 \o 阻频带（尚未撰写） 阻频带则逐渐下降为零。这种滤波器最先由 \o 英国 英国 \o 工程师 工程师斯替芬·巴特沃斯（Stephen Butterworth）在 \o 1930年 1930年发表在英国《无线电工程》期刊的一篇论文中提出的，可以构成低通、高通、带通和带阻四种组态，是目前最为流行的一类数字滤波器 ,经过离散化可以作为数字巴特沃思滤波器 ,较模拟滤波器具有精度高、稳定、灵活、不要求阻抗匹配等众多优点 ,因而在自动控制、语音、图像、通信、雷达等众多领域得到了广泛的应用，是一种具有最大平坦幅度响应的低通滤波器。 1.2巴特沃斯低通滤波器的设计原理: 巴特沃斯低通滤波器的幅度平方函数用下式表示： 其中 N为滤波器的阶数。当=0时，=1；=时，=1/,是3dB截止频率。=时，逐渐增大，幅度下降非常迅速。、N同幅度特性关系如图1.1所示。N决定了幅度下降速度，N越大，通带就越平坦，过渡带也随之变窄，阻带幅度同过渡带下降的速度越迅速，总体010.707N=2N=4N=8频响特性同理想低通滤波器的实际误差越小。 0 1 0.707 N=2 N=4 N=8 图1.1 、N同幅度特性关系 用s代替，把幅度平方函数变成s的函数： s=，此公式说明了幅度平方函数有2N个极点，极点可以用下面的公式来表达： k=0,1,2，···，2N-1。2N个极点等间隔分布在半径为的圆上，间隔是/N rad。如图1.2所示： 0 0 1 图1.2 三阶巴特沃斯滤波器极点分布 为形成稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成，而右半平面的N个极点构成。的表示式为 例如N=3,通过下式可以计算出6个极点 ，，，，， 当N=3时，6个极点中位于左半平面的三个分别为: ，， 取s平面左半平面的极点组成： 将对3dB截止频率Ωc 归一化后的表示为： 令，p称为归一化拉氏复变量。 ，称为归一化频率。 经过归一化后巴特沃斯滤波器的传输函数为： (1.1) 式中，为归一化极点，为位于左半平面的极点用下式表示： 将极点表示式代入(1.1)式，得到的的分母是p的N阶多项式，用下式表示： 下面来确定N： 由技术指标， 和确定。 在定义 (1.2) (1.3) (1.4) 中，将Ω=和Ω=分别代入(1.4)式中，得到和，再将和代入(1.2)和(1.3)式中，得到： 整理得： (1.5)