巴特沃斯滤波器的设计与仿真

  • 格式:doc
  • 大小:159.00 KB
  • 文档页数:14

下载文档原格式

  / 27
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

信号与系统课程设计论文

摘要

传统的数字滤波器的设计过程复杂,计算工作量大,滤波特性调整困难,影响了它的应用。本文介绍了一种利用matlab提供的巴特沃斯滤波器设计函数“buttord”可求出所需的滤波阶数和3dB 截止频率的方法。利用matlab设计滤波器设计函数,可以随时对比设计要求和滤波器特性调整参数,直观简便,极大的减轻了工作量,有利于滤波器设计的最优化。

关键词:巴特沃斯滤波器 Matlab 截止频率

I

信号与系统课程设计论文

Abstract

Design for traditional digital filteris very complicated ,count also very complicated.it is very hard to readjust filtering character, influence apply,the main body of the book introduce use matlab provide butwosto design function“buttord” filteris 3dB end e matlab to design could compete reqire parameter burden,contribute to optimize.

Keywords: filteris matlab end frequency

II

信号与系统课程设计论文

目录

摘要 (Ⅰ)

Abstrct

第1章绪论 (1)

1.1 课题背景 (1)

第2章巴特沃斯滤波器的设计 (2)

2.1 巴特沃斯滤波器阶数的选择 (2)

2.2 巴特沃斯滤波器系数计算 (2)

2.2.1 巴特沃斯低通滤波器系数计算 (2)

2.2.2巴特沃斯高通滤波器系数计算 (3)

2.2.3巴特沃斯带通滤波器系数计算 (3)

2.2.4巴特沃斯带阻滤波器系数计算 (4)

第3章巴特沃斯滤波器设计仿真 (5)

3.1巴特沃斯滤波器设计仿真 (5)

3.1.1 巴特沃斯低通滤波器实例仿真 (5)

3.1.2巴特沃斯高通滤波器实例仿真 (6)

3.1.3巴特沃斯带通滤波器实例仿真 (7)

3.1.4巴特沃斯带阻滤波器实例仿真 (8)

结论 (9)

参考文献 (10)

III

第1章绪论

1.1课题背景

巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦,没有起伏,而在阻频带则逐渐下降为零。在振幅的对数对角频率的波得图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。巴特沃斯滤波器特别适用于低频应用,其对于维护增益的平坦性来说非常重要。这种滤波器最先由英国工程师斯替芬·巴特沃斯(Stephen Butterworth)在1930年发表在英国《无线电工程》期刊的一篇论文中提出的。在电力系统微机保护和二次控制中,很多信号的处理与分析都是基于对正弦基波和某些整次谐波的分析,而系统电压电流信号(尤其是故障瞬变过程)中混有各种复杂成分,所以滤波器一直是电力系统二次装置的关键部件。目前微机保护和二次信号处理软件主要采用数字滤波器。传统的数字滤波器设计使用繁琐的公式计算,改变参数后需要重新计算,在设计滤波器尤其是高阶滤波器时工作量很大。利用matlab提供的巴特沃斯滤波器设计函数“buttord”可以快速有效的实现数字滤波器的设计与仿真。巴特沃斯滤波器是电子滤波器的一种。【1】

第2章巴特沃斯滤波器的设计

2.1巴特沃斯滤波器阶数的选择

ω,sω,Rp,Rs之后,利用MATLAB提供的巴特在已知设计参数p

沃斯滤波器设计函数“buttord”即可求出所需要的滤波器阶数和3dB 截止频率。“buttord”命令的格式是:[n,Wn]= buttord(Wp,Ws,Rp,Rs)其中,Wp是通带截止频率,如Wp=500/4000或Wp=[650,1200]/4000等等。Ws是阻带起始频率,Wp也可以是标量或2个元素的向量。Rp是通带内波动(dB)。Rs是阻带内最小衰减(dB)。返回值n是巴特沃斯滤波器最低阶数,不同的n 有不同的频率响应。n越大频率响应越理想,但n越大。根据要求的指标是多少再选择相应的阶数。这里我选用2阶就够了。Wn是巴特沃斯滤波器截止频率,Wn可以是标量或2个元素的向量。【2】

2.2巴特沃斯滤波器系数计算

[n,Wn]= buttord(Wp,Ws,Rp,Rs)

其中Wp和Ws分别是通带和阻带的拐角频率(截止频率),其

取值范围为0至1之间。当其值为1时代表采样频率的一半。Rp 和Rs分别是通带和阻带区的波纹系数。不同类型(高通、低通、带通和带阻)滤波器对应的Wp和Ws值遵循一定的规则。【3】

2.2.1巴特沃斯低通滤波器系数计算

低通滤波器遵循的规律:Wp和Ws为一元矢量且Wp<Ws。稍后我会在仿真部分阐述设计一个巴特沃斯低通滤波器,通带截止频率为2KHz,阻带截止频率为3KHz,通带波纹系数为1,阻带衰减系数为20,采样频率为10KHz。[b,a]= buttord(n,Wn) % n为低通滤波器阶数

% Wn为低通滤波器截止频率

% b为H(z)的分子多项式系数

% a为H(z)的分母多项式系数

2.2.2 巴特沃斯高通滤波器系数计算

高通滤波器是容许高频信号通过、但减弱(或减少)频率低于截止频率信号通过的滤波器。对于不同滤波器而言,每个频率的信号的减弱程度不同。它有时被称为低频剪切滤波器;在音频应用中也使用低音消除滤波器或者噪声滤波器。高通滤波器与低通滤波器特性恰恰相反。一个滤波器滤除一个复杂信号中不想要的低频成份同时让高频信号通过是很有用的。当然,'低'和'高'频率的含义是相对于滤波器设计者所选择的截止频率而言的。利用巴特沃斯(Butterworth)高通滤波器对图像进行锐化处理。

高通滤波器遵循规律:Wp和Ws为一元矢量且Wp>Ws。

[b,a]= buttord(n,Wn,‘high’)

% n为高通滤波器阶数

% Wn为高通滤波器截止频率

% b为H(z)的分子多项式系数

% a为H(z)的分母多项式系数

2.2.3 巴特沃斯带通滤波器系数计算

带通滤波器是指能通过某一频率范围内的频率分量、但将其他范围的频率分量衰减到极低水平的滤波器,与带阻滤波器的概念相对。这些滤波器也可以用低通滤波器同高通滤波器组合来产生。实际上,并不存在理想的带通滤波器。滤波器并不能够将期望频率范围外的所有频率完全衰减掉,尤其是在所要的通带外还有一个被衰减但是没有被隔离的范围。这通常称为滤波器的滚降现象,并且使用每十倍频的衰减幅度dB来表示。通常,滤波器的设计尽量保证滚降范围越窄越好,这样滤波器的性能就与设计更加接近。然而,随着滚降范围越来越小,通带就变得不再平坦—开始出现“波纹”。

带通滤波器遵循规律:Wp和Ws为二元矢量且Wp<Ws,如