《自-基于MATLAB 一种IIR数字带通滤波器的设计与仿真》由会员分享，可在线阅读，更多相关《自-基于MATLAB 一种IIR数字带通滤波器的设计与仿真（9页珍藏版）》请在金锄头文库上搜索。
1、基于MATAB一种IIR数字带通滤波器的设计与仿真摘要:根据一实例提出基于ALA的两种设计IIR数字带通滤波器的方法,给出了TB完整程序与Sool工具界面设计的详细步骤，并用Siulink进行仿真。关键词：ALA；数字滤波器;IIR，SPTool，设计；仿真；l T Applicaon of MTLA in heDeign and imultin ofne kind IIR igial Dadps Filter Abstact：Tispaerpresents tw meth ason MTAB to desn IRdigial bnpss filtr effectvely acodin o o
2、 xal。he detaled program and proces o designin an inerface by us he PToI isgvn，and t ftr simulatibsed o ALA lnguag nd Simuinki so itrdcd in tis paper。 wods:MALB；digitl filtr；R;Sol;desgn;sui; 引言随着信息时代和数字世界的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。目前数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理中起着重要
3、的作用并已获得广泛应用的是数字滤波器（DF，igital Fler),根据其单位冲激响应函数的时域特性可分为两类:无限冲激响应IR(Ininite mulse Rspons)滤波器和有限冲激响应IR(Fnte Implse Repnse)滤波器。与FIR滤波器相比,IIR的实现采用的是递归结构，极点须在单位圆内，在相同设计指标下,实现IR滤波器的阶次较低,即所用的存储单元少，从而经济效率高。 MATLAB是英文MTrix LABorator(矩阵实验室)的缩写。它是美国的athWors公司推出的一套用于科学计算和图形处理可视化、高性能语言与软件环境。它的信号处理工具箱包含了各种经典的和现代的数
4、字信号处理技术，是一个非常优秀的算法研究与辅助设计的工具。在设计数字滤波器时，通常采用MAT来进行辅助设计和仿真。本文以MATLAB.1为设计平台，利用MLAB完全工具函数El与SPTol工具进行IIR数字带通滤波器的设计,并加以仿真。2 数字滤波器的基本概念滤波器的种类很多,从功能上可分为低通、高通、带通和带阻滤波器，每一种又有模拟滤波器和数字滤波器两种形式。如果滤波器的输人和输出都是离散时间信号，则该滤波器的冲击响应也必然是离散的,这种滤波器称之为数字滤波器。数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。数字滤波器也是具有一定传输选择特性的
5、数字信号处理装置,其输入、输出均为数字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。它的基木工作原理是利用离散系统特性对系统输人信号进行加工和变换，改变输入序列的频谱或信号波形，让有用频率的信号分量通过，抑制无用的信号分量输出。数字滤波器除了具有数字信号处理的固有优点外,还有滤波精度高（与系统字长有关)、稳定性好(仅运行在与I两个电平状态)、灵活性强等优点。 3 数字滤波器设计的基本步骤数字滤波器的设计要经过如下三个步骤：（1）确定指标： 在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标的形式一般在
6、频域中给出幅度和相位响应。（2）模型逼近：一旦确定了技术指标，就可利用已学习过的基本原理和关系式，提出一个滤波器模型来逼近给定的指标体系。这是滤波器设计所要研究的主要问题。(）实现性能分析和计算机仿真 以上两步的结果得到的滤波器，通常是以差分方程、系统函数或脉冲响应描述的。根据这个描述就可以分析其频率特性和相位特性、以验证设计结果是否满足指标要求;或者利用计算机仿真实现设计的滤波器，再分析滤波结果来判断。4 设计实例4.1 MATLAB程序设计法MATLAB信号处理工具箱提供了各种滤波器设计函数及滤波器实现函数，lip是完全设计函数之一,其调用格式为b，=ellip(，Rp,R,n,ptons
7、)，为滤波器阶数，Rp表示通带波纹最大衰减，Rs表示阻带波纹最小衰减，n为滤波器的归一化截止频率（Nyist频率为1Hz)；函数默认为低通或带通滤波器：低通滤波器时,n为截止频率；带通滤波器时,截止频率Wn=W1,W2。a分别为滤波器传递函数分子和分母系数向量;opions为滤波器类型参数:ig为高通滤波器,截止频率为Wn;sop为带阻滤波器，截止频率n1，W2。例:产生有三个正弦成分（5Hz，1 Hz和0z)的信号,设计一滤波器来去除z和30的正弦信号,保留15 的信号。根据题目要求结合所学知识首先选用elip函数直接设计此滤波器(IR数字带通滤波器），其完整MATLAB程序如下：图二：II
8、R带通滤波器同幅频响应 图四：滤波前、后的幅频图图一：含有3个正弦分量的信号 图三：滤波后的信号波形1) 产生含有3个正弦分量的信号（如图一所示)Fs=100;t=(1：10)F;s1=sin(2*i*t5);2sin(2*p*t1);s3sin（*pi*t*30);s=s+;plo(t,s)；ile(含有3个正弦分量的信号);t(gcf，coor,wite);xlabel(Te(second); ylael(Tie wveform）;2) 产生一个4阶IR带通滤波器,通带为10Hz到20Hz,并得出其幅频响应(如图二所示）b，a=ei(4，1,，0,20*2s);H,=freqz(,a，12
9、）;plot(*Fs/(2pi)，s(H);titl(IR带通滤波器幅频响应)；st(gcf,colr,w);label(Frequeny（Hz)； ylae(Mag.oequeny rspne);gd；3) 对信号进行滤波（滤波后信号如图三所示）4)sflter(b,a,s);t(t,sf);tle(滤波后的信号波形);et(f,oor,whe）；xll(Time(econs); ylabe（ime werm); axis(01 -1 1);5) 绘出信号滤波前、后的幅频图（如图四所示)S=t（s，512);SFfft(s，52); =(:55）/256*（F/2);plot(w,ab（S(
10、1:26),SF(1：);title(滤波前、后的幅频图)；et(gcf,color,whit）;xlabel（rqen(Hz);ylabel(Mag.offequency resonse);grid;eged（滤波前的幅频，滤波后的幅频)4.2利用STl工具进行设计Pl是信号处理工具箱中一个具有交互式图形用户界面的信号处理工具,专门用于完成常用的数字信号处理任务。通过这个工具,只需要鼠标简单的操纵鼠标（点击或拖动），就可以完成载入、观察、分析、实现和设计数字滤波器并进行谱分析等甚至十分复杂的数字信号处理任务,而不需要用户对数字滤波器的设计原理非常熟悉。利用PTol工具进行本例的设计有如下几个
11、步骤:1) 创建并导入信号源;在MATAB命令窗口输入命令：Fs=100;t(1:100)/Fs;s= sin(2pi*5) i(2*pi*t15) n（2*i*t*30)；此时，变量Fs、t、将显示在MATLA的workspace列表中。在MTLA命令窗口中输入spool,即可弹出So的主界面（如图五所示),点击菜单则打 图五：SPTooI的主界面 图六：滤波器设计界面开了“Impot to SPTool”对话框，将信号s和采样频率F导入并命名为ig_s。 2) 单击Fits栏中的New按钮打开滤波器设计窗口（如图六所示),在界面的最上面指定滤波器的名字ilt1，采样频率Fs和设计方法Ell
12、t I;左侧的peiicatio面板用来指定设计指标;order为4，tp为adpas,aand中的FP1为1，F2为20,RP为.1，spbnd为40；完成以上输入之后单击ly按钮，则完成了滤波器的设计；中间的rquency espo面板则根据设计指标自动绘出滤波器的幅频响应,而右侧的easureents面板显示了滤波器设计完成后的实测参数。此时在Filters栏中，选中filt1按下View按钮即弹出Vool界面（如图七所示)，我们还可以观察到滤波器的相位响应、幅频相位响应、脉冲响应、阶跃响应等等（只须点击界面上的相应快捷按钮)。3) 将滤波器应用到s信号序列。分别在igls、Filter
13、栏中选择s、ilt1单击Filtrs栏列表下的ply按钮,在弹出的ApplyFite对话框中将输出信号命名为Sig_s（滤波后信号）。同时按下Alt与Shift键选中Sig与Sg_sf点击Sls栏下的Vie可观察到它们的时域波形,如图八所示。4) 进行频谱分析。在Sgnals中选择Si_sf，单击Specr栏下的Crate按钮,在弹出的Specra ieer界面中选择mtho为F,nfft=512，单击pply按钮生成s的频谱spectsf,同理得到s的频谱pct_,同时按下Alt与Shf键选中ect_s与sect_f点击Spectr栏下的Viw可观察到它们的频谱,如图九所示。 图七：FVTool界面 图九：滤波前、后的信号的频谱 图八：滤波前、后信号波形由图八可以看出,所设计的II带通滤波器il使输入信号s中频率为15h的正弦波信号通过,而将频率为5hz和0h的正弦波信号大大衰减,从而可知所设计的滤波器满足题目要求。.3 两种方法的比较 将Ellip函数设计的滤波器从工作空间导入到Po工具的Filters栏中，通过滤波器分析工具FVToo可以看出两种方法设计出来的滤波器性能完全相同,相比于利用滤波器设
