MATLA在数字滤波器设计中的应用课程设计

1、matlab程序设计与应用课 程 论 文matla在数字滤波器设计中的应用课 程 名 称 matlab程序设计与应用 考 查 学 期 第三学年第一学期 考 查 方 式 课程论文 姓 名 陈婉 学 号专 业 通行工程 成 绩 指 导 教 师 万力 摘要：传统的数字滤波方法仅仅依靠外部信号的正负信息，根本没有办法对信号进行准确的滤波处理，得不到真正的外部信号信息，而且设计过程复杂，计算工作量大，滤波特性调整困难，影响了它的应用。然而数字滤波是语音和图像处理、模式识别、频谱分析等应用中的一个基本处理算法，在数字信号处理中占有极其重要的地位。本文利用matlab设计fir滤波

2、器，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化。本文还介绍了如何利用matlab环境下的仿真软件simulink对所设计的滤波器进行模拟仿真。 关键词：数字滤波器 fdatool matlab 滤波 simulink引言 滤波器的设计方法分为传统的滤波方法和数字滤波方法。传统的数字滤波器设计使用繁琐的公式计算，改变参数后需要重新计算，在设计滤波器尤其是高阶滤波器时工作量很大。利用matlab信号处理工具箱（signal processing toolbox）和simulink可以快速有效的实现数字滤波器的设计与仿真。随着科学技术的发展，数字

3、滤波器在数字信号处理领域中已经得到了广泛应用。一 数字滤波器、matlab和simulink的简要介绍数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波器有多种分类，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应（iir）滤波器和有限长冲激响应（fir）滤波器。iir数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配。所以iir滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。fir数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要

4、求的转移序列或脉冲响应的常数问题，设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。matlab是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案。matlab的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用matlab来解算问题要比用c，fortran等语言完成相同的事情简捷得多，并且mathwork也吸收了像maple等

5、软件的优点,使matlab成为一个强大的数学软件。simulink是matlab中的一种可视化仿真工具， 是一种基于matlab的框图设计环境，是实现动态系统建模、仿真和分析的一个软件包，被广泛应用于线性系统、非线性系统、数字控制及数字信号处理的建模和仿真中。simulink可以用连续采样时间、离散采样时间或两种混合的采样时间进行建模，它也支持多速率系统，也就是系统中的不同部分具有不同的采样速率。为了创建动态系统模型，simulink提供了一个建立模型方块图的图形用户接口(gui) ，这个创建过程只需单击和拖动鼠标操作就能完成，它提供了一种更快捷、直接明了的方式，而且用户可以立即看到系统的仿真

6、结果。二 fdatool的介绍fdatool（filter design & analysis tool）是matlab信号处理工具箱里专用的滤波器设计分析工具，可以设计几乎所有的基本的常规滤波器，包括fir和iir的各种设计方法。它操作简单，方便。fdatool界面总共分两大部分，一部分是design filter，在界面的下半部，用来设置滤波器的设计参数，另一部分则是特性区，在界面的上半部分，用来显示滤波器的各种特性。design filter部分主要分为：filter type（滤波器类型）选项，包括lowpass（低通）、highpass（高通）、bandpass（带通）、bandst

7、op（带阻）和特殊的fir滤波器。design method（设计方法）选项，包括iir滤波器的butterworth（巴特沃思）法、chebyshev type i（切比雪夫i型）法、 chebyshev type ii（切比雪夫ii型） 法、elliptic（椭圆滤波器）法和fir滤波器的equiripple法、least-squares（最小乘方）法、window（窗函数）法等。filter order（滤波器阶数）选项，定义滤波器的阶数，包括specify order（指定阶数）和minimum order（最小阶数）。在specify order中填入所要设计的滤波器的阶数（n阶滤波

8、器，specify ordern-1），如果选择minimum order则matlab根据所选择的滤波器类型自动使用最小阶数。frenquency specifications选项，可以详细定义频带的各参数，包括采样频率fs和频带的截止频率。它的具体选项由filter type选项和design method选项决定，例如bandpass（带通）滤波器需要定义fstop1（下阻带截止频率）、fpass1（通带下限截止频率）、fpass2（通带上限截止频率）、fstop2（上阻带截止频率），而lowpass（低通）滤波器只需要定义fstop1、fpass1。采用窗函数设计滤波器时，由于过渡带是

9、由窗函数的类型和阶数所决定的，所以只需要定义通带截止频率，而不必定义阻带参数。magnitude specifications选项，可以定义幅值衰减的情况。例如设计带通滤波器时，可以定义wstop1（频率fstop1处的幅值衰减）、wpass（通带范围内的幅值衰减）、wstop2（频率fstop2处的幅值衰减）。当采用窗函数设计时，通带截止频率处的幅值衰减固定为6db，所以不必定义。window specifications选项，当选取采用窗函数设计时，该选项可定义，它包含了各种窗函数。三 带阻滤波器设计实例本文将以一个iir 滤波器的设计为例来说明如何使用matlab设计数字滤波器:设计一个

10、椭圆带阻滤波器。参数要求：系统采样频率fs=8000hz，要求滤除2025-2225hz频段的频率成分，保留0-1500hz和2700hz以上频段的频率成分，衰减大于40db。首先在response type中选择bandstop（带阻滤波器）；在design method选项中选择iir 和elliptic（椭圆滤波器）,指定filter order项中的specify order6，在frequency specifications中把units设置为hz,fpass1=2025,fpass2=2500,设置完以后点击design filter即可得到所设计的椭圆带阻滤波器。通过菜单选项a

11、nalysis可以在特性区看到所设计滤波器的幅频响应、相频响应、零极点配置和滤波器系数等各种特性。设计结果如下图:在设计过程中，可以对比滤波器幅频相频特性和设计要求，随时调整参数和滤波器类型，以便得到最佳效果。其它类型的fir滤波器和iir滤波器也都可以使用fdatool来设计。四 程序设计法上例的带阻滤波器用程序设计：fsl=2025;fsu=2225;fpl=1500;fpu=2700;fs=8000;ws=2*fsl/fs,2*fsu/fs;wp=2*fpl/fs,2*fpu/fs; %计算滤波器的边界频率rp=1;rs=40;n,wpo=ellipord(wp,ws,rp,rs)%调用

12、计算滤波器阶数n和通带截止频率wpob,a=ellip(n,rp,rs,wpo,stop)%调用ellip计算带阻滤波器系统函数系数向量b和a freqz(b,a);用matlab运行，结果如下：n =3wpo =0.3750 0.6750b = 0.3521 0.1819 1.0416 0.3576 1.0416 0.1819 0.3521a =1.0000 0.3531 1.0547 0.3038 0.7232 0.0645 0.0095五 simulink仿真 本文通过调用simulink中的功能模块构成数字滤波器的仿真框图，在仿真过程中，可以双击各功能模块，随时改变参数，获得不同状态下

13、的仿真结果。例如构造以基波为主的原始信号，通过simulink环境下的digital filter design（数字滤波器设计）模块，把digital filter design模块放到*mdl文件中，双击设置上例中所设置的带阻滤波器，然后把我所需要的模块一一添加进去。仿真图和滤波效果图如图2。 仿真图 滤波之前的信号滤波之后的信号结束语利用matlab的强大运算功能，基于matlab信号处理工具箱（signal processing toolbox）的数字滤波器设计法可以快速有效的设计由软件组成的常规数字滤波器，设计方便、快捷，极大的减轻了工作量。在设计过程中可以对比滤波器特性，随时更改参

14、数，以达到滤波器设计的最优化。参考文献：1 楼顺天、李博菡. 基于matlab的系统分析与设计信号处理【m】西安：西安电子科技大学出版社，1998.2 高西全、丁玉美. 数字信号处理：西安电子科技大学出版社，2010.3 飞思科技产品研发中心. matlab7基础与提高：电子工业出版社，2005.ut2apodfxxc02gybkskcww97mrqqwhoj5tl15zt6jipyytycummtarp3v1n5luizi3xh3bhwyreko8d9g7nmzqowpjetldrw08gvs8dsdqqygc3ce7moo2tlf0jf1gk74iuxybmtivr97ckrfvqult5

15、fn2t6mpjr6rbzvpsortzvij5nb5ndvvsr4iwr1twlfkglspzuhrjq3cmzu98euouijdlszqpmvrw9zkupxf8wfug9l2g9277g2rtipa1ypczeuqxpkbhtvdcooqozxuz3vjrzmocijym62zchmeootyes8ebmm932tbz2yo09rtszeys8zrd2yktj8l6jeazvajnfbtrylvsm6ofbftoxvrffn7owiygjlamkunxjybz5rrb7r4vsur9zpfzfmfsjhcfca37lnw2vvlrkn7r8psz1bn6oric5hu5z6hcxayq

16、ynpog8duybawqsl20csg06dh2sm8hltgpkicskrgopdpuhbj1lmpk7lydvc6nnmwl3fwhzftfvyaary7lhssxj10v3ph3y19bxyr77ib7cpzsu2tijqe3hkqkkau9kskcphkxuikvvyjzpg2yijrkqfbggovyqkuxnwi9omnjtt6qilzxtyrf7d20fbmabcfiixrqkusvnxbppfuxyq1fjskfsubkgs2duvqc9sz4jkbgn4qqv66pyoarjurnfj3txyfclzieeptwfjthpheipdfnqnr2hjqkv2dzwtmpdjq

17、kbcxmovdsjqctjagjmdlskpgad2s0h0vmzgaht36gyuez7umank1ndreubeqdgrx0venqgnsyib2ilq3siqrnl4m56t7z8y8da5k0kupn5nzg4jvjdtffhyt82aogqkxo4vblmleiy2p7hthbho07rcfttxodydppdtqso7wxd0j6fkklgm4wodzplhtrr2xgqn13hqy59zu1gegdyqnihntavsieuefqcyfucjwd3vk5i7ykmhundmiz ut2apodfxxc02gybkskcww97mrqqwhoj5tl15zt6jipyytycum

18、mtarp3v1n5luizi3xh3bhwyreko8d9g7nmzqowpjetldrw08gvs8dsdqqygc3ce7moo2tlf0jf1gk74iuxybmtivr97ckrfvqult5fn2t6mpjr6rbzvpsortzvij5nb5ndvvsr4iwr1twlfkglspzuhrjq3cmzu98euouijdlszqpmvrw9zkupxf8wfug9l2g9277g2rtipa1ypczeuqxpkbhtvdcooqozxuz3vjrzmocijym62zchmeootyes8ebmm932tbz2yo09rtszeys8zrd2yktj8l6jeazvajnfbtrylvsm6ofbftoxvrffn7owiygjlamkunxjybz5rrb7r4vs