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

通信系统综合设计与实践题目基于MATLAB的低通滤波器设计院（系）名称信院通信系专业名称通信工程学生姓名学生学号指导教师2013 年 5 月25 日摘要 (2)1．巴特沃斯低通数字滤波器简介 (3)1.1 选择巴特沃斯低通滤波器及双线性变换法的原因 (4)1.2 巴特沃斯低通滤波器的基本原理 (4)1.2.1 巴特沃斯低通滤波器的基本原理 (4)1.2.2 双线性变换法的原理......... . (5)1.3 数字滤波器设计流程图......... .. (7)1.4 数字滤波器的设计步骤.......... . (7)2. 巴特沃斯低通数字滤波器技术指标的设置 (8)3. .................................................................................................................. 用matlab 实现巴特沃斯低通数字滤波器的仿真并分析 .. (9)3.1 用matlab 实现巴特沃斯低通数字滤波器的仿真 (9)3.2 波形图分析......... ........ (10)4. .................................................................................................................. 用Simulink 实现巴特沃斯低通数字滤波器的仿真并分析 (11)4.1 Simulink 简介........ ....... .. (11)4.2 用Simulink 实现巴特沃斯低通数字滤波器的仿真 (11)4.3 波形图分析........ ........ . (11)4.3.1 Simulink 波形图分析........... ............ .. 114.3.2 与matlab 波形的比较........ . (14)5. 总结与体会 ................. . (15)6. 附录 (16)摘要低通滤波器是让规定频率以下的信号分量通过，而对该频率以上的信号分量抑制的电容、电感与电阻等器件的组合装置。

三峡大学课程设计报告专业班级 20091421课程数字信号处理课程设计学号 2009142116 学生姓名姜祥奔指导教师王露2012年 5 月平时成绩（20%）报告成绩（40%）答辩成绩（40%）总成绩数字信号处理课程设计实验一：用双线性变换法和脉冲响应不变法设计IIR 数字滤波器采用双线性变换法设计一个巴特沃斯数字低通滤波器，要求：通带截止频率100Hz ，阻带截止频率200Hz ，通带衰减指标Rp 小于2dB ，阻带衰减Rs 大于15dB ，滤波器采样频率Fs=500Hz 。

绘制频率响应曲线。

理论部分：（原理及设计过程）第一步：论ω和f 的关系及数字域性能的公式表示。

模拟频率与数字频率之间为线性关系，T f T fT T s,101,24-===Ω=πω为抽样周期，故 ππω4.0500110021=⨯⨯==c c kHz f 对应于 修正后)2t a n (21c T ω=Ω ππω8.0500120021=⨯⨯==st st kHz f 对应于 修正后)2tan(21st T ω=Ω按衰减的定义2)()(log 204.0010≤πj j e H e H 15)()(log 208.0010≥πj j e H e H 设0=ω处频率响应幅度归一化为1,即1)(0=j e H ，则上两式变成2)(log 204.010-≥πj e H (1)15)(log 208.010-≤πj e H (2)这就是数字滤波器的性能指标的表达式。

2 下面把数字低通滤波器的性能要求转变为“样本”模拟低通滤波器的性能要求。

由T Ω=ω，按修正式)()(Tj H e H a j ωω≈，设没有混叠效应（即混叠效应设计完成后再进行校验） 则有πωωω≤Ω==),()()(j H TjH eH a a j (3)利用（3）式，由（1）、（2）式可写出模拟低通滤波器的指标为2)102(log 20)4.0(log 2021010-≥⨯=ππj H TjH a a …….（4） 15)104(log 20)8.0(log 2021010-≤⨯=ππj H TjH a a …（5） 3 计算“样本”模拟低通滤波器所需的阶数N 及3dB 截止频率C Ω。

用MATLAB 设计巴特沃斯低通滤波器1 巴特沃斯低通滤波器的特性一个理想低通滤波器的幅频特性如图3-80的阴影部分所示。

为了实现这个理想低通特性，需要在从0~ωC 的整个频带内增强增益，在ω>ωC 增益要降到0。

实际上，理想滤波器是不可能实现的。

图3-78是实际滤波器的幅频特性。

但是实际滤波器的特性愈接近理想特性愈好，巴特沃斯（Butterworth ）滤波器就是解决这个问题的方法之一。

巴特沃斯滤波器以巴特沃斯函数来近似滤波器的系统函数，巴特沃斯的低通模平方函数为：221|()|1,2,,1(/)NC H j N j j ωωω==+ (3-138)式中以C ω是滤波器的电压-3dB 点或半功率点。

不同阶次的巴特沃斯滤波器特性如图3-79(a)所示。

4阶巴特沃斯滤波器的极点分布如图3-79(b)所示。

巴特沃斯滤波器幅频响应有以下特点：• 最大平坦性：在0=ω附近一段范围内是非常平直的，它以原点的最大平坦性来逼近理想低通滤波器。

• 通带、阻带下降的单调性。

这种滤波器具有良好的相频特性。

• 3dB 的不变性：随着N 的增加，频带边缘下降越陡峭，越接近理想特性。

但不管N 是多少，幅频特性都通过-3dB 点。

• 极点配置在半径为ωC 的圆上，并且均匀分布。

左半平面上的N 个极点是)(s H 的极点，右半平面上的N 个极点是)(s H -的极点。

2 巴特沃斯低通滤波器的实现为使巴特沃斯滤波器实用，我们必须能够实现它。

一个较好的方法是将巴特沃斯滤波器函数化成若干二阶节级联，其中每一节实现一对共轭复极点。

通过将极点以共轭复数的形式配对，对所有的每一个二阶节都具有实系数。

1图3-78 低通滤波器的幅频特性图3-80所示运算放大器电路为实现一对共轭极点提供了很好的方法。

电路的系统函数为202202121121122121)(1)11(1)(ωωω++=+++=s Qs C C R R s C R C R s C C R R s H (3-139)式中，ω0是S 平面原点与极点之间的距离，Q 被称为电路的“品质因数”，它提供了对响应峰值尖锐程度的一种度量。

利用MATLAB设计巴特沃斯低通数字滤波器引言数字滤波器是数字信号处理中的重要组成部分，可以用于去除信号中的噪音和不需要的频率成分。

巴特沃斯滤波器是一种常见的数字滤波器，被广泛应用于信号处理领域。

本文将介绍如何利用MATLAB设计巴特沃斯低通数字滤波器，并给出详细的步骤和示例代码。

设计步骤利用MATLAB设计巴特沃斯低通数字滤波器主要包括以下步骤：1.设计滤波器的参数2.计算滤波器的传递函数3.绘制滤波器的幅频响应曲线4.通过频域图像观察滤波器的性能下面将分别介绍每个步骤的详细操作。

设计滤波器的参数巴特沃斯低通数字滤波器的参数包括截止频率和阶数。

截止频率决定了滤波器的通频带，阶数决定了滤波器的陡峭程度。

通过MATLAB的butter()函数可以方便地设计巴特沃斯低通数字滤波器。

该函数的参数为滤波器的阶数和截止频率。

示例代码如下：order = 4; % 阶数cutoff_freq = 0.4; % 截止频率[b, a] = butter(order, cutoff_freq);计算滤波器的传递函数通过设计参数计算得到滤波器的传递函数。

传递函数是一个复数，包括了滤波器的频率响应信息。

使用MATLAB的freqz()函数可以计算滤波器的传递函数。

该函数的参数为滤波器的系数b和a，以及频率取样点的数量。

示例代码如下：freq_points = 512; % 频率取样点数量[h, w] = freqz(b, a, freq_points);绘制滤波器的幅频响应曲线经过计算得到的传递函数能够提供滤波器的幅频响应信息。

通过绘制幅频响应曲线，可以直观地观察滤波器的频率特性。

使用MATLAB的plot()函数可以绘制滤波器的幅频响应曲线。

该函数的参数为频率点和传递函数的幅值。

示例代码如下：magnitude = abs(h); % 幅值plot(w/pi, magnitude);xlabel('归一化频率');ylabel('幅值');title('巴特沃斯低通数字滤波器幅频响应');通过频域图像观察滤波器的性能通过绘制滤波器的频域图像，可以直观地观察滤波器对不同频率的信号的响应情况。

巴特沃斯低通、切比雪夫低通、高通IIR滤波器设计05941401 1120191454 焦奥一、设计思路IIR滤波器可以分为低通、高通、带通、带阻等不同类型的滤波器，而以系统函数类型又有巴特沃斯、切比雪夫等滤波器。

其中巴特沃斯较为简单，切比雪夫较为复杂；低阶比高阶简单，但却有着不够良好的滤波特性。

在满足特定的指标最低要求下，低阶、巴特沃斯滤波器能更大程度地节省运算量以及复杂程度。

滤波器在不同域内分为数字域和模拟域。

其中数字域运用最广泛。

在设计过程中，一般是导出模拟域的滤波器，之后通过频率转换变为数字域滤波器，实现模拟域到数字域的传递。

在针对高通、带通、带阻的滤波器上，可以又低通到他们的变换公式来进行较为方便的转换。

综上，IIR滤波器的设计思路是，先得到一个满足指标的尽可能简单的低通模拟滤波器，之后用频域变换转换到数字域。

转换方法有双线性变换法、冲激响应不变法等。

虽然方法不同，但具体过程有很多相似之处。

首先将数字滤波器的指标转换为模拟滤波器的指标，之后根据指标设计模拟滤波器，再通过变换，将模拟滤波器变换为数字滤波器，是设计IIR滤波器的最基本框架。

以下先讨论较为简单的巴特沃斯低通滤波器。

二、巴特沃斯低通滤波假设需要一个指标为0~4hz内衰减小于3db、大于60hz时衰减不小于30db的滤波器。

其中抽样频率为400hz。

以双线性变换方法来设计。

首先将滤波器转换到模拟指标。

T =1f f ⁄=1400Ωf ′=2ff f =8ff f =Ωf ′f =0.02fΩf ′=2ff f =120ff f =Ωf ′f =0.3f根据双线性变换Ω=2f tan ⁡(f 2) 得到Ωf =25.14Ωf =407.62这就得到了模拟域的指标。

由巴特沃斯的方程Α2(Ω)=|f f (f Ω)|2=11+(ΩΩf )2f20ff |f f (f Ω)|=−10ff [1+(ΩΩf)2f] {20ff |f f (f Ωf )|≥−320ff |f f (f Ωf )|≤−30ff得到{ −10ff [1+(Ωf Ωf)2f ]≥−3−10ff [1+(Ωf Ωf )2f]≤−30当N取大于最小值的整数时，解出N=2，因此为二阶巴特沃斯低通滤波器。

基于MATLAB的低通巴特沃斯滤波器仿真作者：刘凌云赵鹏宇弓美桃来源：《数字技术与应用》2013年第02期摘要：本文介绍了低通巴特沃斯滤波器的原理，分析了低通巴特沃斯滤波器的特性，然后用MATLAB对低通巴特沃斯滤波器进行了仿真。

关键词：低通巴特沃斯滤波器 MATLAB 仿真中图分类号：TN911.7 文献标识码：A 文章编号：1007-9416（2013）02-0124-011 巴特沃斯滤波器原理巴特沃斯滤波器是电子滤波器的一种，这种滤波器由英国工程师斯替芬·巴特沃斯（Stephen Butterworth）在1930年提出，它的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。

可以构成低通、高通、带通和带阻四种组态，是目前最为流行的一类数字滤波器，经过离散化可以作为数字巴特沃思滤波器，较模拟滤波器具有精度高、稳定、灵活、不要求阻抗匹配等众多优点，因而在语音、图像、通信、雷达等众多领域得到了广泛应用，是一种具有最大平坦幅度响应的低通滤波器。

2 MATLAB简介MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。

它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案。

MATLAB主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。

在信号处理中低通巴特沃斯滤波器常用的函数包括：2.1 buttord函数[N，wn]=buttord（wp，ws，rp，rs，'s'）用于计算巴特沃斯数字滤波器的阶数N和3dB截止频率wn。

其中，调用参数wp，ws分别为数字滤波器的通带、阻带截止频率的归一化值，要求：0≤wp≤1，0≤ws≤1（1表示数字频率pi）。

基于MATLAB设计巴特沃斯低通滤波器
李钟慎
【期刊名称】《信息技术》
【年(卷),期】2003(027)003
【摘要】首先分析了巴特沃斯低通滤波器的特性,然后用MATLAB的信号处理工具箱提供的函数设计了巴特沃斯低通滤波器,使得巴特沃斯滤波器的设计变得更加简单、快捷、直观.
【总页数】3页(P49-50，52)
【作者】李钟慎
【作者单位】华侨大学机电及自动化学院,泉州,362011
【正文语种】中文
【中图分类】TP311
【相关文献】
1.基于Matlab的巴特沃斯数字低通滤波器的设计 [J], 张廷尉;陈红;王磊
2.基于Pspice的八阶巴特沃斯低通滤波器设计与优化 [J], 叶榆;贺国权
3.基于EWB的巴特沃斯有源低通滤波器的设计与仿真 [J], 张白莉;郭红英
4.基于WEBENCH的巴特沃斯低通滤波器的设计 [J], 徐志国;杨娟
5.基于CFA的巴特沃斯低通滤波器设计和仿真 [J], 刘宛;董金宁;池思慧;赵柏树因版权原因，仅展示原文概要，查看原文内容请购买。

各专业全套优秀毕业设计图纸吉首大学信息科学与工程学院课程设计报告书课程：数字信号处理教程_______________________课题：基于MATLAB的巴特沃斯低通滤波器以切贝雪夫低通滤波器的设计姓名： _________________________学号：—专业：通信工程___________________________________年级：2012 级 __________________________________________指导教师： _______________________________________基地指导教师： _________________________________________2014年12 月一、项目介绍与设计目的1. 通过实验加深对巴特沃斯低通滤波器以及切贝雪夫低通滤波器基本原理的理解。

2. 学习编写巴特沃斯低通滤波器和切比雪夫的MATLA仿真程序。

3. 滤波器的性能指标如下：通带截止频率fp=3kHz，通带最大衰减Rp=2dB阻带截止频率fst=6kHz，阻带最小衰减As=30db二、设计方案1. 项目环境要求MATLAB软件2. 设计内容一.理论设计：模拟巴特沃思低通滤波器的设计1.有技术指标可求的设计参数Q p=6000 n rad/s; Q st=12000n rad/s;Rp=2db; As=30 db2■求N;2 1g(Q p/QJ=5.369 取N=63■确定参数=1971.114. 求系统函数Han( s) =1/(1+3.8637033S+7.4641016S A2+9.141620S A3+7.4641016S A4+3.8637033S A5+S A6) 5. 去归一化(I严叶_］严xlabel( 'f(kHz)' );ylabel( 'dB' );axis([-1,12,-55,1])set(gca, 'xtickmode' , 'manual' ，‘xtick' ,[0,1,2,3,4,5,6,7,8,9,10]); Set(gca, 'ytickmode' , 'maunal' , 'ytick' ，[-50,-40,-30,-2,-10,0]);grid;输出结果N=6b=5.8650e+25a=1 7.5158e+04 2.9000e+09 7.0010e+13 1.1267e+18 1.1496e+22 5.8650e+25dbHx=2.0000 33.7692在matlab 中显示如下：a =L De+025 *0.0000 0.0000 0. 0000 0. 0000 0. 0000 0. 0011 b =5.86505.8650e+025dbHx =2.0000 33.7962程序运行结果如下故系统函数为：Ha(s)=5.8650*10A25/(s A6+7.6158*10A4s A5+2.9000*10A9s A4+7.0010*10A13s A3+1.1267*10*18 s A2+1.1496*10A32s+5.8650*10A25)(2)切贝雪夫型低通滤波器程序：clc , clear all ;OmegaR=2*pi*3000;OmegaS=2*pi*6000;Rp=2;As=30; g=sqrt((10A(As/10)-1”(10A(Rp/10)-1)); OmegaR=OmegaS/OmegaPN=ceil(log10(g+sqrt(g*g-1))/log10(OmegaR+sqrt(OmegaR*OmegaR-1)))OmegaC=OmegaS;[zO,pO,kO]=cheb2ap(N,As);aO=real(poly(pO))aNn=aO(N+1);p=pO*OmegaC;a=real(poly(p))aNu=a(N+1);bO=real(poly(zO));M=le ngth(bO);bNn=bO(M);z=zO*OmegaC;b=real(poly (z));bNu=b(M);k=k0*(aNu*b Nn )/(a Nn *bNu);b=k*bwO=[OmegaP,OmegaS];[H,w]=freqs(b,a);Hx=freqs(b,a,wO);dbHx=-20*log10(abs(Hx)/max(abs(H))) plot(w/(2*pi)/1000,20*log10(abs(H)));xlabel( 'f(kHz)' );ylabel( 'dB' );axis([-1,12,-55,1])set(gca, 'xtickmode' , 'manual' , 'xtick' ,[0,1,2,3,4,5,6,7,8,9,10]);Set(gca, 'ytickmode' , 'maunal' , 'ytick' ，[-50,-40,-30,-20,-10,0]);grid;输出结果:N=4b=0.0316 0 3.5954e+08 0 5.1099e+17A=1 6.6533e+04 2.2247e+09 4.3659e+13 5.1099e+17程序运行结果如下图此低通滤波器的系统函数为：Ha(s)=0.0316s A4+3.5954*10A8S A2+5.1099*10A17/(S A4+6.6533*10A4S A3+2.2247*10A9S A2+4.3659*10A13S+5.1099*10A17)三、总结和分析通过本次butterworth低通滤波器的设计，使我对低通数字滤波器的工作原理和特性有了深刻的认识和了解，实验主要用到了matlab软件。

基于MATLAB设计巴特沃斯低通滤波器课程设计专业：XXXXXX姓名：XXX学号： XX指导老师：XXX2011年11 月26日通信系统仿真课程设计任务书院（系）：电气信息工程学院目录1 绪论 (1)1.1 引言 (1)1.2 数字滤波器的设计原理 (1)1.3 数字滤波器的应用 (2)1.4 MATLAB的介绍 (3)1.5 本文的工作及安排 (3)2 滤波器分类及比较 (4)2.1 滤波器的设计原理 (4)2.2 滤波器分类 (4)2.3 两种类型模拟滤波器的比较 (6)3 巴特沃斯低通滤波器 (7)3.1 巴特沃斯低通滤波器简介 (7)3.2 巴特沃斯低通滤波器的设计原理 (7)4 MATLAB仿真及分析 (11)4.1 MATLAB工具箱函数 (11)4.2 巴特沃斯低通滤波器的MATLAB仿真 (11)另附程序调试运行截图： (13)5.1 总结 (13)5.2 展望 (13)1 绪论1.1 引言凡是有能力进行信号处理的装置都可以称为滤波器。

滤波器在如今的电信设备和各类控制系统里面应用范围最广、技术最为复杂，滤波器的好坏直接决定着产品的优劣。

自60年代起由于计算机技术、集成工艺和材料工业的发展，滤波器发展上了一个新台阶，并且朝着低功耗、高精度、小体积、多功能、稳定可靠和价廉方向努力，其中小体积、多功能、高精度、稳定可靠成为70年代以后的主攻方向。

使以数字滤波器为主的各种滤波器得到了飞速的发展，到70年代后期，数字滤波器的单片集成已被研制出来并得到应用。

80年代，致力于各类新型滤波器的研究，努力提高性能并逐渐扩大应用范围。

90年代至现在主要致力于把各类滤波器应用于各类产品的开发和研制。

当然，对数字滤波器本身的研究仍在不断进行。

[1]滤波器主要分成经典滤波器和数字滤波器两类。

从滤波特性上来看，经典滤波器大致分为低通、高通、带通和带阻等。

本文主要对低通数字滤波器做主要研究。

1.2 数字滤波器的设计原理所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相对比例，或者滤除某些频率成分的数字器件或程序。

基于MATLAB 设计巴特沃斯低通滤波器摘 要： 首先分析了巴特沃斯低通滤波器的特性。

然后用MATLAB 的信号处理工具箱提供的函 数设计了巴特沃斯低通滤波器，使得巴特沃斯滤波器的设计变得更加简单、快捷、直观。

巴特沃斯(Butterworth)滤波器是一种具有最大平坦幅度响应的低通滤波器，它在通信领域里已有广应用，在电测中也具有广泛的用途，可以作检测信号的滤波器。

MATLAB 语言是一种面向科学与工程计算的语言。

它编程效率高，测试程序手段丰富，扩展能力强，内涵丰富。

它的信号处理工具箱(Signal Processing Toolbox)提供了设计巴特沃斯滤波器的函数，本文充分利用这些函数，进行了巴特沃斯滤波器的程序设计，并将其作为函数文件保存，可方便地进行调用。

1. 巴特沃斯低通滤波器的特性巴特沃斯低通滤波器的平方幅度响应为： n c j H 22)(11)(ωωω+=其中，n 为滤波器的阶数，ωc 为低通滤波器的截止频率。

该滤波器具有一些特殊的性质：① 对所有的n ，都有当 ω=0时，|H(j0)|2 =1；② 对所有的n ，都有当ω=ωc 时，|H(j ωc )|2 =0.5 ，即在ωc 处有3dB 的衰减；③|H(j ω)|2 是ω的单调递减函数，即不会出现幅度响应的起伏；④ 当n →+∞时，巴特沃斯滤波器趋向于理想的低通滤波器；⑤ 在ω=0处平方幅度响应的各级导数均存在且等于0，因此|H(j ω)|2 在该点上取得最大值，且具有最大平坦特性。

图l 展示了2阶、4阶、8阶巴特沃斯低通滤波器的幅频特性。

可见阶数n 越高，其幅频特性越好，低频检测信号保真度越高。

巴特沃斯与贝塞尔(Besse1)、切比雪夫(Cheby.shev)滤波器的特性差异如图2所示。

从图2可以看出，巴特沃斯滤波器在线性相位、衰减斜率和加载特性三个方面具有特性均衡的优点，因此在实际使用中，巴特沃斯滤波器已被列为首选。

2 .巴特沃斯低通滤波器的MATLAB实现MATLAB的信号处理工具箱提供了有关巴特沃斯滤波器的函数buttap、buttord、butter。

fp=0.2*pi;fs=0.3*pi;Rp=1;Rs=15;Fs=1000;Wp=2*pi*fp;Ws=2*pi*fs;Nn=256;n=(0:100-1);[N,Wn]=buttord(Wp,Ws,Rp,Rs,'s'); [b,a]=butter(N,Wn,'s');w=linspace(1,300,100)*2*pi;H=freqs(b,a,w);figure(1);plot(w/(2*pi),20*log10(abs(H)));title('巴特沃斯模拟滤低通波器幅频特性'); xlabel('频率/Hz');ylabel('幅度/db');[bz,az]=impinvar(b,a,Fs);figure(2);zplane(bz,az);title('冲击响应不变法零极点分布图') figure(3)Subplot(2,1,1);freqz(bz,az,Nn,Fs);grid on;title('巴特沃斯数字滤波器幅频特性'); xlabel('幅度/Hz');ylabel('频率/db'); Subplot(2,1,2);title('巴特沃斯数字滤波器相频特性'); xlabel('幅度/Hz');ylabel('频率/db');grid on;ws1=0.2*pi;wp1 = 0.35*pi;wp2 = 0.65*pi;wc = [wp1/pi,wp2/pi];B=wp1-ws1;N = ceil(8/0.15);n=0:N-1;%%rectangle窗window=rectwin(N);[h1,w]=freqz(window,1);subplot(121);stem(window);xlabel('n');title('rectangle窗函数');hn = fir1(N-1,wc,Boxcar(N));[h2,w]=freqz(hn,1,512);figure(1)subplot(122);plot(w/pi,20*log(abs(h2)/abs(h2(1)))); grid;xlabel('w/pi');ylabel('幅度(dB)');title('rectangle滤波器的幅度特性'); ;%%blackman窗window=hamming(N);[h1,w]=freqz(window,1);figure(2)subplot(121);stem(window);xlabel('n');title('hamming窗函数');hn = fir1(N-1,wc,hamming(N));[h2,w]=freqz(hn,1,512);subplot(122);plot(w/pi,20*log(abs(h2)/abs(h2(1)))); grid;xlabel('w/pi');ylabel('幅度(dB)');title('hamming滤波器幅度特性'); %%hanning窗window=hanning(N);[h1,w]=freqz(window,1);figure(3)subplot(121);stem(window);xlabel('n');title('blackman窗函数');hn = fir1(N-1,wc,hanning(N));[h2,w]=freqz(hn,1,512);subplot(122);plot(w/pi,20*log(abs(h2)/abs(h2(1))));grid;xlabel('w/pi');ylabel('幅度(dB)');title('hanning滤波器的幅度特性'); %%blackman窗window=blackman(N);[h1,w]=freqz(window,1);figure(4)subplot(121);stem(window);xlabel('n');title('blackman窗函数');hn = fir1(N-1,wc,blackman(N));[h2,w]=freqz(hn,1,512);subplot(122);plot(w/pi,20*log(abs(h2)/abs(h2(1)))); grid;xlabel('w/pi');ylabel('幅度(dB)');title('blackman滤波器的幅度特性'); 第四题第五题。

分类号编号烟台大学毕业论文基于MATLAB设计巴特沃斯低通滤波器The Design of Butterworth Low-passing Filter Based on MA TLAB申请学位：院系：专业：姓名：学号：指导老师：2011年05 月26日烟台大学基于MA TLAB设计巴特沃斯低通滤波器姓名：导师：2011年05月26日烟台大学烟台大学毕业论文任务书院（系）：光电信息科学技术学院[摘要]滤波器设计是数字信号处理的重要内容。

在MATLAB软件中有丰富的滤波器设计的相关命令，掌握相关的方法后可以提高我们的工作效率。

首先对巴特沃斯低通滤波器的特性进行研究，然后用MATLAB信号处理工具箱提供的函数设计出巴特沃斯低通滤波器模型，并对具体实例进行分析，使得巴特沃斯滤波器的设计更加快捷、直观、简单。

[关键词]巴特沃斯低通滤波器； MATLAB仿真；[Abstract]First，analyse the characteristics of Butterworth low-pass filter, second use MATLAB signal processing toolbox design the mode of Butterworth low - pass filter ，to study it though an explme. The method makes the design of Butterw orth filter quicklier ,more intuitively,and simp -lier.[Keywords] Butterworth low-pass filter; MATLAB simulation;目录1 绪论 (1)1.1 引言 (1)1.2 数字滤波器的设计原理 (1)1.3数字滤波器的应用 (2)1.4MATLAB的介绍 (3)1.5本文的工作及安排 (3)2 滤波器分类及比较 (5)2．1滤波器的设计原理 (5)2.2 滤波器分类 (5)2.3四种类型模拟滤波器的比较 (9)3巴特沃斯低通滤波器 (11)3.1巴特沃斯低通滤波器的设计原理 (11)4 MATLAB仿真及分析 (15)4.1 MATLAB工具箱函数 (15)4.2 巴特沃斯低通滤波器的MATLAB仿真 (15)5 结论与展望 (19)5.1 总结 (19)5.2 展望 (19)致谢 (20)参考文献 (21)1 绪论1.1 引言凡是有能力进行信号处理的装置都可以称为滤波器。

基于m a t l a b的巴特沃斯低通滤波器的实现基于matlab 的巴特沃斯低通滤波器的实现一、课程设计的目的运用MATLAB实现巴特沃斯低通滤波器的设计以及相应结果的显示，另外还对多种低通滤波窗口进行了比较。

二、课程设计的基本要求1）熟悉和掌握MATLAB 的基本应用技巧。

2）学习和熟悉MATLAB相关函数的调用和应用。

3）学会运用MATLAB实现低通滤波器的设计并进行结果显示。

三、双线性变换实现巴特沃斯低通滤波器的技术指标：1.采样频率10Hz。

2.通带截止频率fp=0.2*pi Hz。

3.阻带截止频率fs=0.3*pi Hz。

4.通带衰减小于1dB，阻带衰减大于20dB四、使用双线性变换法由模拟滤波器原型设计数字滤波器程序代码：T=0.1; FS=1/T;fp=0.2*pi;fs=0.3*pi;wp=fp/FS*2*pi;ws=fs/FS*2*pi;Rp = 1; % 通带衰减As = 15; % 阻带衰减OmegaP = (2/T)*tan(wp/2); % 频率预计OmegaS = (2/T)*tan(ws/2); % 频率预计%设计巴特沃斯低通滤波器原型N = ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(OmegaP/OmegaS))); OmegaC = OmegaP/((10^(Rp/10)-1)^(1/(2*N)));[z,p,k] = buttap(N); %获取零极点参数p = p * OmegaC ;k = k*OmegaC^N;B = real(poly(z));b0 = k;cs = k*B; ds = real(poly(p));[b,a] = bilinear(cs,ds,FS);% 双线性变换figure(1);% 绘制结果freqz(b,a,512,FS);%进行滤波验证figure(2); % 绘制结果f1=50; f2=250;n=0:63;x=sin(2*pi*f1*n)+sin(2*pi*f2*n);subplot(2,2,1);stem(x,'.'); title ('输入信号');y=filter(b,a,x);subplot(2,2,2);stem(y,'.') ;wd1=boxcar(N)';b1=hd.*wd1;wd2=hanning(N)';b2=hd.*wd2;wd3=blackman(N)';b3=hd.*wd3;wd4=hamming(N)';b4=hd.*wd4;[h1,w]=freqz(b1,1);[h2,w]=freqz(b2,1);[h3,w]=freqz(b3,1);[h4,w]=freqz(b4,1);plot(w,20*log10(abs(h1)),w,20*log10(abs(h2)),':',w,20*log10(abs(h3)),'-.',w,20*log10(abs(h4)),'--')legend('矩形窗','汉宁窗','布莱克曼窗','汉明窗');结果显示;00.51 1.52 2.53 3.5-160-140-120-100-80-60-40-20020六、总结及心得:基于matlab 的巴特沃斯低通滤波器的实现姓名：学号： S20100 专业：日期： 2010年06月08日。

第一章数字滤波器及MATLAB语言概述数字滤波是语音和图像处理、模式识别、谱分析等应用中的一个基本处理算法，在数字信号处理中占有极其重要的地位。

研究基于Matlab 环境下的IIR数字滤波器的设计与实现，给出了相应的Matlab函数命令，并将滤波器应用于图像噪声的去除，取得了不同的效果，就其结果做出了进一步的解释和说明。

数字滤波器是具有一定传输选择性的数字信号处理装置，其输入、输出均为数字信号，实质上是一个由有限精度算法实现的线性不变离散系统。

它的基本工作原理是利用离散系统特性对系统输入信号进行加工和变换，改变输入序列的频谱或信号波形，让有用频率的信号分量通过，抑制无用的信号输出分量。

MATLAB是美国MathWorks公司推出的一套用于工程计算的可视化高性能语言与软件环境。

MATLAB为数字滤波的研究和应用提供了一个直观、高效、便捷的利器。

它以矩阵运算为基础，把计算、可视化、程序设计融合到了一个交互式的工作环境中。

MATLAB推出的工具箱使各个领域的研究人员可以直观方便地进行科学研究、工程应用，其中的信号处理(signalproeessing)、图像处理(imageproeessing)、小波(wavelet)等工具箱为数字滤波研究的蓬勃发展提供了有力的工具。

数字滤波器与模拟滤波器相比，具有精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配以及能实现模拟滤波器无法进行的特殊滤波等优点。

本文主要介绍有限冲激响应数字滤波器（FIR）和无限冲激响应数字滤波器（IIR）的设计原理、方法、步骤以及在MATLAB中的实现，并以实例形式列出设计程序和仿真结果。

第二章 方案提取和技术要求第一节 方案提取IIR 数字滤波器是一种离散时间系统，其系统函数为假设M ≤N ，当M ＞N 时,系统函数可以看作一个IIR 的子系统和一个(M-N)的FIR 子系统的级联。

IIR 数字滤波器的设计实际上是求解滤波器的系数和 ，它是数学上的一种逼近问题，即在规定意义上（通常采用最小均方误差准则）去逼近系统的特性。