《课程设计（论文）基于MATLAB的IIR滤波器设计.doc》由会员分享，可在线阅读，更多相关《课程设计（论文）基于MATLAB的IIR滤波器设计.doc（22页珍藏版）》请在三一办公上搜索。
1、物理与电子信息工程系 学生学年论文(设计)题 目 基于MATLAB的IIR滤波器设计 学生姓名 学 号 专 业 电子信息工程年 级 10级 指导教师 二一 二 年 月 日物理与电子信息工程系本科生学年论文指导记录、指导教师审阅意见表学生姓名学号院、系物理与电子信息工程学院专 业电信班级10级指导教师职 称论文题目基于MATLAB的IIR数字滤波器的设计指导记录：指导教师签名： 学生签名： 年 月 日指导记录：指导教师签名： 学生签名： 年 月 日指导记录：指导教师签名： 学生签名： 年 月 日指导记录：指导教师签名： 学生签名： 年 月 日指导教师审阅意见：（对学年论文完成情况及质量、工作能力
2、及态度、思想表现、论文学术水平等进行总体评价）成 绩 ： 指导教师签名： 年 月 日基于MATLAB的IIR数字滤波器的设计摘要：利用MATLAB 设计滤波器,可以按照设计要求非常方便地调整设计参数,极大地减轻了设计的工作量,有利于滤波器设计的最优化。Matlab因其强大的数据处理功能被广泛应用于工程计算,其丰富的工具箱为工程计算提供了便利,利用Matlab信号处理工具箱可以快速有效地设计各种数字滤波器,设计简单方便。本文介绍了在MATLAB R2007a 环境下滤波器设计的方法和步骤。关键词：滤波器, matlab, fdatool, IIR目录摘要1第一章、引言11.1目的11.2意义1第
3、二章、IIR滤波器的设计理12.1数字滤波器 1第三章、滤波器的技术标2第四章、IIR型滤波器的设计34.1.巴特沃斯低通滤波器设计34.2.模拟高通滤波器的设计 64.3.模拟带通滤波器的设计 84.4.模拟带阻滤波器的设计10第五章、模拟与数字滤波器的转换方法（双线性变换法）11第六章、结论17参文考献 18第一章、引言1.1目的利用MATLAB 设计IIR数字滤波器1.2意义 利用Matlab信号处理工具箱可以快速有效地设计各种数字滤波器,设计简单方便。第二章、IIR滤波器的设计原理2.1数字滤波器滤波器是指用来对输入信号进行滤波的硬件和软件。所谓数字滤波器，是指输入、输出均为数字信号，
4、通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。数字滤波器和模拟滤波器相比，因为信号的形式和实现滤波的方法不同，数字滤波器具有比模拟滤波器精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配等优点。 一般用两种方法来实现数字滤波器:一是采用通用计算机，把滤波器完成的运算编成程序通过计算机来执行，也就是采用计算机软件来实现;二是设计专用的数字处理硬件。 MATLAB的信号处理工具箱是专门应用于信号处理领域的专用工具箱，它的两个基本组成就是滤波器的设计与实现部分以及谱分析部分。工具箱提供了丰富而简便的设计，使原来繁琐的程序设计简化成函数的调用。只要以正确的指标参数调用相
5、应的滤波器设计程序或工具箱函数，便可以得到正确的设计结果，使用非常方便。2.2.数字滤波器的分类 数字滤波器从功能上分类:可分为低通滤波器、高通滤波器、带通滤波器、带阻滤波器。 从滤波器的网络结构或者从单位脉冲响应分类:可分为IIR滤波器(即无限长单位冲激响应滤波器)和FIR滤波器(即有限长单位冲激响应滤波器)。它们的函数分别为:第一个公式中的H (z)称为N阶IIR滤波器函数，第二个公式中的H (z)称为(N-1)阶FIR滤波器函数。第三章、滤波器的技术指标数字滤波器的传输函数为：式中，|称为幅频特性函数，Q()称为相频特性函数。幅频特性表示信号通过该滤波器后各频率成分的衰减情况，而相频特性
6、反映各频率通过滤波器后在时间上的延时情况。一般IIR数字滤波器，通常只用幅频响应函数，相频特性一般不作要求。 1 1- 0.707 0 图3.1低通滤波器的技术要求 图3.1表示可实现的低通滤波器的幅频特性，和分别称为通带截止频率和阻带截止频率。通带频率范围为，其不一定水平，在通带中要求(1-)1；阻带频率范围为，其不一定衰减到零，在阻带中要求。从到称为过渡带，一般为单调下降的。通带内和阻带内允许的衰减一般用dB数表示，通带内允许的最大衰减用表示，阻带内允许的最小衰减用表示，和用dB分别定义为 如将归一化为1，式和式表示为当幅度下降到时，此时=3dB，称为3dB通带截止频率，它是滤波器设计的重
7、要参数之一。第四章、IIR型滤波器的设计IIR型滤波器设计的最通用的方法是借助于模拟滤波器的设计方法来完成的。IIR型滤波器的设计步骤如下：（1）按一定规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标。（2）根据转换后的技术指标设计模拟低通滤波器。（3）在按一定规则将转换成。（4）将高通、带通、或带阻数字滤波器的技术指标转化为低通模拟滤波器的技术指标，然后按上述步骤（2）设计出低通，再将转换为所需的。4.1.巴特沃斯低通滤波器设计Butterworth（巴特沃斯）低通滤波器的设计方法可按以下三个步骤来进行：（1）.确定阶数N将实际频率归一化，得到归一化幅频平方特性式中，N为滤波器
8、的阶数；为3dB截止频率。当= 时,=当 时，=1,即 在处无衰减当0时，随着迅速衰减，N越大，衰减的速度越快，即过渡带越窄。将=和=分别代入式，再将所得的幅频平方函数代入式可得：1+ 1+根据式和式可得令，可求得： N=当N有小数时，取大于N的最小整数（2）求得归一化传输函数将幅度平方函数写成s的函数：可知幅度平方函数有2N个极点： 2N个极点等间隔分布在半径为的圆上，称该圆为巴特沃斯圆。为了形成稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成，其表达式为对3dB截止频率归一化，归一化后的为 式中， ，令 , , p称为归一化拉氏复变量，得巴特沃斯传输函数为式中,为归一化极点表示为
9、因为N已求得所以也可以通过查表来确定（3）求出,将去归一化得到实际的滤波器传输函数将 带入中得到用MATLAB实现:Wp=30;Ws=50;Ap=1;As=50;Fs=200;n,Wn=buttord(Wp,Ws,Ap,As,s);%选择Butterworth滤波器参数b,a=butter(n,Wn,s);%采用双线性变化法进行离散化处理bn1,an1=bilinear(b,a,2000);%双线性变化H1,W=freqs(bn1,an1);plot(W,abs(H1);grid;xlabel(频率);ylabel(幅度响应);4.2.模拟高通滤波器的设计设低通滤波器 和高通滤波器的幅频特性如
10、图4.2.1、4.2.2所示。图中、分别称为低通滤波器的归一化通带截止频率和归一化阻带截止频率，、分别称为高通滤波器的归一化通带下限频率和归一化阻带上限频率。 0 图4.2.1 低通滤波器的幅频特性 0 图4.2.2 高通滤波器的幅频特性由于和都是频率的偶函数，可以把右边曲线和曲线对应起来，低通从经过和到0时，高通的则从0经过和到因此和之间的关系为： 从低通滤波器转换到模拟高通滤波器的设计步骤如下。（1）确定高通滤波器的技术指标：通带下限频率，阻带上限频率。通带最大衰减，阻带最小衰减 （2）确定相应低通滤波器的技术指标：按照式，将高通滤波器的边界频率转换成低通滤波器的边界频率，各项设计指标为低
11、通滤波器的通带截止频率 低通滤波器的阻带截止频率 通带最大衰减仍为，阻带最小衰减仍为（3）设计归一化低通滤波器（4）求模拟高通的。将按照式转换成归一化高通，为去归一，将 代入中，得上式，就是由归一化低通直接转换成模拟高通的转换公式。用MATLAB实现:%模拟高通演示%直接采用MATLAB设计的高通模拟滤波器Wp=2*pi*200;Ws=2*pi*100;Ap=2;As=15;%HPF性能指标n,Wn=buttord(Wp,Ws,Ap,As,s);%获得滤波器参数b,a=butter(n,Wn,high,s);%设计Butterworth HPEfigure(1);freqs(b,a,20000
12、);%画出设计出的HPE频率响应曲线Wc=2*pi*200;b1=1 0 0 0;a1=1 2*Wc 2*(Wc.2) Wc.3figure(2);freqs(b1,a1,20000);%画出理论求出的HPE频率响应特性曲线 4.3.模拟带通滤波器的设计低通滤波器与带通滤波器的幅频特性如图所示 0 图4.3.1 带通滤波器的幅频特性 0 图4.3.2 低通滤波器的幅频特性从低通滤波器转换到模拟带通滤波器的设计步骤如下。（1）确定带通滤波器的技术指标：带通上限频率，带通下限频率；下阻带上限频率，上阻带下限频率；通带中心频率， ；通带宽带。与以上边界频率对应的归一化边界频率如下： 将带通和低通的幅
13、度特性对应起来，可以得到和的对应关系， 还需要确定的技术指标有：通带最大衰减阻带最小衰减（2）确定归一化模拟低通技术要求： 与的绝对值可能不相等，一般取绝对值较小的，这样保证在较大的处能满足要求。通带最小衰减仍为，阻带最小衰减亦为 （3）设计归一化模拟低通（4）直接将转换成带通滤波器由于 代入式得 将 代入上式，得 去归一化，将 代入上式，得因此4.4.模拟带阻滤波器的设计低通滤波器与带阻滤波器的幅频特性如图所示。 0 图4.4.1低通滤波器的幅频特性 0 图4.4.2带阻滤波器的幅频特性从低通滤波器转换到模拟带阻滤波器的设计步骤如下。（1）确定模拟带阻滤波器的设计要求：上通带截止频率 ，下通
14、带截止频率，阻带下限频率，阻带上限频率，阻带中心频率，阻带宽度。以上边界频率对应的归一化边界频率如下 将阻带和低通的幅度特性对应起来，可以得到和的对应关系还需要确定的技术指标有：通带最大衰减,阻带最小衰减（2）确定归一化模拟低通技术要求： 取与的绝对值较小的，这样保证在较大的处能满足要求。通带最大衰减仍为，阻带最小衰减亦为 （3）设计归一化模拟低通（4）直接将转换成阻带滤波器可以得到因此第五章、 模拟与数字滤波器的转换方法（双线性变换法）双线性变换是采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到之间，再用 转换到Z平面上。设, ,经过非线性频率压缩后用， 表示，用正切变换实现频率压缩：
15、 式中，T为采样间隔,当从 经过0变化到 时,则由 经过0变化到 实现了s平面上整个虚轴的 之间的变换。这样便有再通过转换到z平面上。得到用双线性变换法设计一个巴特沃斯低通数字滤波器，频率f在500Hz处衰减3dB，在750Hz处衰减至少15dB，采样频率为2000Hz,确定系统函数H(z)。解：（1）求数字指标和： （2）求、 ：（3）确定滤波器阶数：（4）确定系统函数：利用MATLAB完成实现%数字低通双线性变化Wp=2*pi*500;Ws=2*pi*150;Ap=3;As=15;Fs=2000;n,Wn=buttord(Wp,Ws,Ap,As,s);%选择Butterworth滤波器参数
16、b,a=butter(n,Wn,s);%采用双线性变化法进行离散化处理bn1,an1=bilinear(b,a,2000);%双线性变化H1,W=freqs(bn1,an1);plot(W,abs(H1);grid;xlabel(频率);ylabel(幅度响应);假设原始信号为： ，由于某些原因，信号被一个另外频率的信号 及白噪声污染，这样实际的信号变为 。不妨取采样频率为 ，通过下面的MATLAB演示，可以清楚地看到这种滤波选择性滤波器的功能。%频率选择性滤波器功能实现Fs=1000;t=0:1/Fs:2;%取2秒长度的信号x0=sin(2*pi*80*t)+2*sin(2*pi*140*t
17、);%原始信号xN=sin(2*pi*300*t);x=x0+xN+randn(size(t);%受污损信号x(t)figure(1);subplot(211);plot(t,x0);xlabel(t);ylabel(x0(t);title(原始信号x0(t);subplot(212);plot(t,x);xlabel(t);ylabel(x(t);title(受污损信号x(t);%设计一个特定的频率选择性滤波器n=100;%滤波器改为取100f=0 0.13 0.15 0.17 0.19 0.25 0.27 0.29 0.31 1;m=0 0 1 1 0 0 1 1 0 0 ;b=firls
18、(n,f,m);H,W=freqz(b,1,512,2);figure(2);plot(W,20*log(abs(H);gridxlabel(归一化频率);ylabel(滤波器的对数幅频特性);%完成滤波功能x1=filter(b,1,x);figure(3);subplot(211);plot(t,x1);xlabel(t);ylabel(x1(t);title(滤波输出信号x1(t);subplot(212);plot(t,x1-x0);xlabel(t);ylabel(x1(t)-x0(t);title(误差信号xe(t);第六章、结论采用MATLAB设计滤波器，使原来非常繁琐复杂的程序
19、设计变成了简单的函数调用，为滤波器的设和实现开辟了广阔的天地，尤其是Matlab工具箱使各个领域的研究人员可以直观方便地进行科学研究与工程应用。其中的信号处理工具箱、图像处理工具箱、小波工具箱等更是为数字滤波研究的蓬勃发展提供了可能。MATLAB 信号处理工具箱为滤波器设计及分析提供了非常优秀的辅助设计工具, 在设计数字滤波器时, 善于应用MATLAB进行辅助设计, 能够大大提高设计效率。参文考献1倪养华,王重玮编著， 数字信号处理与实现，上上海交通大学出版社，1998。2 施阳，MATLAB语言工具箱，西安：西北工业大学出版社,1999。3薛年喜，MATLAB在数字信号处理中的应用，清华大学
20、出版社，2003.11.1。4张葛祥、李 娜，MATLAB仿真技术与应用，北京：清华大学出版社，2003。5 陈桂明，应用MATLAB语言处理数字信号与数字图像，北京：科学出版社，2001。Using MATLAB design IIR牛强 10131280129赤峰学院 物理与电子信息工程学院 赤峰 024000ABSTRACTBy using MATLAB , we can design filters and modify the filtersparameters conveniently according to our demands. This relieves greatly
21、design work loads and makes for optimization of filter designing. Matlab can be widely used in engineering calculations because of its powerful functions of data processing. Its rich toolbox makes the calculations easy. With Matlab signal processing toolbox, various digital filters can be designed effectively in simple way. This article introduce the methods and processes in the circumstance of MATLAB R2007a.Keywords：filter , matlab, fdatool, IIR
