毕业论文范本
本科毕业论文(设计)
题 目 基于MATLAB 窄带带通滤波的设计 院(系) 电子工程与电气自动化学院 专 业 电子科学与技术 学生姓名 宋 俊 涛 学 号 10029034 指导教师 李素平 职称 讲师 论文字数 8210
完成日期:2014年5月10日
2
巢湖学院本科毕业论文(设计) 诚信承诺书
本人郑重声明:所呈交的本科毕业论文(设计) ,是本人在导师的指导下,独立进行研究工作所取得的成果。除文中已经注明引用的内容外,本论文不含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体,均已在文中以明确方式标明。本人完全意识到本声明的法律结果由本人承担。
本人签名:
日期: 年 月 日
巢湖学院本科毕业论文 (设计) 使用授权说明
本人完全了解巢湖学院有关收集、保留和使用毕业论文 (设计) 的规定,即:本科生在校期间进行毕业论文(设计) 工作的知识产权单位属巢湖学院。学校根据需要,有权保留并向国家有关部门或机构送交论文的复印件和电子版,允许毕业论文 (设计) 被查阅和借阅;学校可以将毕业论文(设计) 的全部或部分内容编入有关数据库进行检索,可以采用影印、缩印或扫描等复制手段保存、汇编毕业,并且本人电子文档和纸质论文的内容相一致。
保密的毕业论文(设计) 在解密后遵守此规定。
本人签名:
导师签名: 日期: 年 月 日 日期: 年 月 日
基于MATLAB 窄带带通滤波器的设计
摘 要
随着现代网络、通信及无线电技术的快速发展,频率资源变得日益紧张,
滤波器作为分离有用信号和干扰信号的重要部件,其性能的优劣直接影响信息系统的质量。论文讨论窄带带通滤波器的设计问题,这对实际应用具有重要的指导意义,文中讨论了FIR 滤波器和IIR 滤波器设计中的原理方法和它们之间的优缺点。为了深入讨论这些问题,论文利用MATLAB 语言进行窄带带通滤波器的设计和并利用Simulink 模块进行带通滤波器的仿真。
设计中选用FIR 滤波器和IIR 滤波器进行设计,设计的结果突出了它们之
间的特点,对设计的方法和算法有启示意义,不同性能指标应选用不同设计方法以实现实际应用中的经济效益。
关键词:FIR ;IIR ;MATLAB ;Simulink
Narrow Band Flter Design Based on MATLAB
Abstract
Frequency resources get more and more nervously with the modern
communication network and radio technology fast development. Filter is an important section to take part of the useful signal and disturbed signal. And the filter ’s function influences the signal system directly. The thesis discusses the design question of narrow band bandpass filter. It ’s important to guide reality utilization. The paper discussed FIR filter’s and IIR filter’s theory of the designing way and the advantage and shortcoming of them. Making use of MATLAB and simulink model toolbox of MATLAB to design the narrow band bandpass filter .
In the article, FIR filter and IIR filter are designed. The characteristics of them
are discussed in the conclusion. We should use different design ways for different aims in the reality utilization.
Key word: FIR ; IIR ; MATLAB; simulink
研究背景
目前网络信息,移动通信,卫星导航等产业迅速繁荣,对国民经济乃至世
界各国经济的发展起着巨大的推动作用,深刻影响着人们的生活方式。然而这些产业无一不是占用着一定的频率资源,因此频率资源变得日益紧张,滤波器作为分离有用信号和干扰信号的重要部件,其性能的优劣直接影响通信系统的质量。滤波器在电子系统的各个方面,在航天,导航等重要领域都有应用。我们知道每一通信线路都占用一段非常短的频率段,在通信过程中,接收端之间的通道的抗噪滤波性能的质量是保证通信畅通的重要部分,因此研究带通滤波器更具有研究意义。
数字信号处理技术是当今飞速发展的一门学科,她更是以不同的形式影响
着和渗透到其他学科,它与国民经济息息相关,与国防建设紧密相连,影响和改变着人们的生产,生活方式。数字滤波技术是数字信号分析、处理技术的重要分支。无论是信号的获取、传输,还是信号的处理和交换都离不开滤波技术,它对信号安全可靠和有效灵活地传输是至关重要的。在所有的电子系统中,使用最多技术最复杂的要算数字滤波器了。数字滤波器的优劣直接决定产品的优劣。MATLAB 是数字信号处理应用和计算的重要工具,利用MATLAB 可对滤波器的性能和特点进行研究设计。
1. 滤波器
1.1滤波器的概念
滤波器是一种用来消除干扰杂讯的器件,将输入或输出经过过滤而得到纯净的交流电。您可以通过基本的滤波器积木块——二阶通用滤波器传递函数,推导出最通用的滤波器类型:低通、带通、高通、帯阻和椭圆型滤波器。
1.2滤波器的发展过程 滤波器就是能进行信号的处理和提取的装置。在近代电信装备和各类控制
系统中,滤波器应用极为广泛;在所有的电子部件中,使用最多,技术最复杂
要算滤波器了。滤波器的优劣直接决定产品的优劣,所以,对滤波器的研究和生产历来为各国所重视。
1917年美国和德国科学家分别发明了LC 滤波器,次年导致了美国第一个
多路复用系统的出现。50年代无源滤波器日趋成熟。自60年代起由于计算机技术、集成工艺和材料工业的发展,滤波器发展上了一个新台阶,并且朝着低功耗、高精度、小体积、多功能、稳定可靠和价廉方向努力,其中小体积、多功能、高精度、稳定可靠成为70年代以后的主攻方向,导致RC 有源滤波器 、数字滤波器、开关电容滤波器和电荷转移器等各种滤波器的飞速发展。到70年代后期,上述几种滤波器的单片集成被研制出来并得到应用。80年代致力于各类新型滤波器性能提高的研究并逐渐扩大应用范围。90年代至今在主要致力于把各类滤波器应用于各类产品的开发和研制。当然,对滤波器本身的研究仍在不断进行。
我国广泛使用滤波器是50年代后的事,当时主要用于话路滤波和报路滤
波。经过半个世纪的发展,我国滤波器在研制、生产应用等方面已有一定进步,但由于缺少专门研制机构,集成工艺和材料工业跟不上来,使许多新型滤波器的研制应用与国际水平有一段距离。
2滤波器原理
一个给定的输入输出关系,可以用多种不同的数字网络来实现。在不考虑
量化影响时,这些不同的实现方法是等效的;但在考虑量化影响时,这些不同的实现方法性能上就有差异。因此,运算结构是很重要的,同一系统函数H (z ) ,运算结构的不同,将会影响系统的精度、误差、稳定性、经济性以及运算速度等许多重要性能。IIR(无限冲激响应) 滤波器与FIR(有限冲激响应) 滤波器在结构上有自己不同的特点,在设计时需综合考虑。数字滤波器又可以理解为是一个计算程序或算法,将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列,并在转化过程中,使信号按预定的形式变化。数字滤波器有多种分类,根据数字滤波器冲激响应的时域特征,可将数字滤波器分为两种,即无限长冲激响应(IIR )滤波器和有限长冲激响应(FIR )滤波器。理想带通滤波器幅频特性为:
低通H (e
-2-
高通02H (ej ω
基于MATLAB 窄带带通滤波的设计
-2π-π
带通0πH (ej ω) 2π
ωH (ej ω) 带阻
图2-1 理想带通滤波器幅频特性
-π0 -2π2ω
实际带通滤波器的幅频特性为
图2-2 实际带通滤波器幅频特性
如果δ1和δ2分别为通带和阻带的容限值。在实际设计中,给出的具体指
标一般为通带允许的最大衰减αp 和阻带最小衰减αs 。
(2-1)
(2-2)
在设计中我们需要根据所要设计的各项性能指标,根据已有的知识计算出
各种变量。对于设计一种滤波器要根据所要设计的性能指标来进行选择设计方法,IIR 和FIR 滤波器设计方法各有其优缺点,本论文设计主要是进行两种滤波器的特点的研究及验证,查资料可知两种滤波器的特点如表2-1所示:
表2-1两种滤波器特点比较分析
基于MATLAB 窄带带通滤波的设计
2.1 IIR滤波器
2.1.1 IIR滤波器原理
IIR 数字滤波器具有无限宽的冲激响应,与模拟滤波器相匹配。所以IIR 滤
波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。设计IIR 滤波器有比较成熟的规律和程序化方法,因为与模拟滤波器匹配,设计时又可以采用模拟滤波器设计的思路。IIR 滤波器可以用系统函数表示:
H (z ) =∑b z k
k =0
N
k =1M -k =1-∑a k z -k Y (z ) (2-3) X (z )
由这样的系统函数可以得到表示系统输入与输出关系的常系数线形差分程为
y (n ) =∑a k y (n -k ) +∑b k x (n -k ) (2-4) k =0k =0N M
有差分方程可以看到数字滤波器的功能就是把输入序列x (n ) 通过一定的运
算变换成输出序列y (n ) 。不同的运算处理方法决定了滤波器实现结构的不同。无限冲激响应滤波器的单位抽样响应h (n ) 是无限长的,其差分方程是递归式的,即结构上存在着输出信号到输入信号的反馈,其系统函数在z 平面的有限区间(0
2.1.2 IIR滤波器特点
IIR 滤波器设计各种的方法中,都有其公共的特点:第一,IIR 数字滤波器
的系统函数可以写成封闭的函数形式。第二,IIR 数字滤波器采用递归型结构,即结构上带有反馈环路,IIR 滤波器运算结构通常由延时,乘以系数和相加等基本运算组成,可以组合成直接型,级联型,并联型等结构形式,都具有反馈回路,由于运算中的舍入处理,使误差不断积累,有时会产生微弱的寄生震荡。
基于MATLAB 窄带带通滤波的设计
第三,IIR 数字滤波器在计算上可以借助成熟的模拟滤波器的成果,如巴特沃斯设计法,切比雪夫设计法,椭圆法,贝塞尔法等,有现成的设计数据或图表可查,其设计工作量比较小,对计算工具的要求不高。在设计一个IIR 数字滤波器时,我们根据指标先写出模拟滤波器的公式,然后通过一定的变换,将模拟滤波器的公式转换成滤波器的公式。第四,IIR 数字滤波器的相位特性不好控制,对相位要求较高时,需加相位校准网络。
在于FIR 滤波器相比,最突出的特点是IIR 滤波器的相位不具有线性,若需要线性相位,则要用全通网络进行校正。其脉冲响应序列是无限的,其在物理实现上不可能无限脉冲,因此现实设计的IIR 滤波器会产生混叠效应和泄露效应,在设计同样的滤波器设计指标上,IIR 滤波器所需要的阶数会比FIR 滤波器所需要的阶数低。
2.2 FIR滤波器
2.2.1 FIR滤波器原理
FIR 滤波器,即是有限长单位冲激响应滤波器,是数字信号处理系统中最基本的元件,可以保证任意幅频特性的同时具有严格的线性相频特性其单位冲激响应是有限长的,因此可以构成稳定的滤波器系统。其设计是根据给定滤波器的频率特性,求的满足该特性的传输函数
窗函数设计法:
一般是先给定所要求的理想滤波器频率响应H d (e j ω) ,由H d (e j ω) 导出h d (n ) ,我们知道理想滤波器的冲击响应h d (n ) 是无限长的非因果序列,而我们要设计的是h d (n ) 是有限长的FIR 滤波器,所以要用有限长序列h d (n ) 来逼近无限长序列h d (n ) ,设: h d (n ) =1
2π⎰π-πH d (e j ω) e j ωd ω (2-5)
常用的方法是用有限长度的窗函数w (n ) 来截取
即: h (n ) =ω(n ) h d (n ) (2-6)
基于MATLAB 窄带带通滤波的设计
这里窗函数就是矩形序列R N (n ), 加窗以后对理想低通滤波器的频率响应将产生什么样的影响呢? 根据在时域是相乘关系, 在频域则是卷积关系: H (e j ω) =1
2π⎰π+-πH d (e j ω) W R [e j (ω-θ) ]d θ (2-7)
其中, 为矩形窗谱, H (e j ω) 是FIR 滤波器频率响应.
的幅度函数H (ω) 的起伏现象, 可知, 加窗处理后, 通过频域卷积过程看
对理想矩形的频率响应产生以下几点影响:
(1)使理想频率特性不连续点处边沿加宽,形成一个过渡带,其宽度等于窗的频率响应的主瓣宽度。
(2)在截止频率的两边的地方即过渡带的两边,出现最大的肩峰值,肩峰的两侧形成起伏振荡,其振荡幅度取决于旁瓣的相对幅度,而振荡的多少,则取决于旁瓣的多少。
(3)改变N ,只能改变窗谱的主瓣宽度,改变ω的坐标比例以及改变的绝对值大小,但不能改变主瓣与旁瓣的相对比例(此比例由窗函数的形状决定) 。
(4)对窗函数的要求
a 、窗谱主瓣尽可能窄,以获取较陡的过渡带;
b 、尽量减小窗谱的最大旁瓣的相对幅度;即能量集中于主瓣,使肩峰和波纹减小,增大阻带的衰减。
窗函数设计法是从时域出发,把理想的h d (n ) 用一定形状的窗函数截取成有限长的h (n ) ,来近似理想的h d (n ) ,这样得到的频率响应
的理想的频率响应
频率采样法:
频率抽样法则是从频域出发,把给定的理想频率响应
抽样得到
H (k ) ,即 H d (k ) =H d (e jw ) |ω=2πk N 逼近于所要求。 加以等间隔,然后以此
作为实际FIR 滤波器的频率特性的抽样值(2-8)
基于MATLAB 窄带带通滤波的设计
知道H (k ) 后,由DFT 定义可唯一确定有限长序列 h(n) ,利用这N 个频域抽样值H (k ) 同样利用频率内插公式可得FIR 滤波器的系统函数H (z ) ,及频率响应
,频率抽样法内插公式:
频率抽样法小结:
优点:可以在频域直接设计,并且适合于最优化设计。
缺点:抽样频率只能等于 2π/N 的整数倍,或等于2π/N 的整数倍加上π/N 。因而不能确保截止频率
的自由取值,要想实现自由地选择截止频率,必须增1-z -N H (z ) =N H (k ) ∑-k -11-W z k =0N N -1(2-9) 加抽样点数N ,但这又使计算量增大。
为了提高逼近质量,减少通带边缘由于抽样点的陡然变化而引起的起伏振荡。有目的地在理想频率响应的不连续点的边缘,加上一些过渡的抽样点,增加过渡带,减少起伏振荡。
2.2.2 FIR滤波器特点
FIR 滤波器单位冲激响应h (n ) 的特点:
其单位冲激响应h (n ) 是有限长(
N -1
n =0), 系统函数为: H (z ) =∑h (n ) z -n (2-10)
因此,它在有限Z 平面有(N -1) 个零点,而它的(N -1) 个极点均位于原点z =0 处。
FIR 滤波器线性相位的特点:
如果FIR 滤波器的单位抽样响应h (n ) 为实数,而且满足以下任一条件:
偶对称h (n ) =h (N -1-n )
(2-11) 奇对称h (n ) =-h (N -1-n ) (2-12) 其对称中心在n =(N -1)/2处,则滤波器具有准确的线性相位。
与IIR 滤波器相比,FIR 滤波器有它突出的优点,但同样也有它设计中的缺点,其优点在于:第一,很容易获得严格的线性相位,避免被处理的信号产生
基于MATLAB 窄带带通滤波的设计
相位失真,在宽频带信号处理、阵列信号处、数据传输等应用系统中非常重要。第二,可容易获得多带幅频特性。第三,极点全部在原点,无稳定性问题。第四,任何一个非因果的有限长序列,总可以通过一定的延时,转变为因果序列,所以因果性总是满足。第五,无反馈运算,运算误差小。它的缺点在于:第一,因为没有极点,要获得好的过渡带特性,需以较高的阶数为代价。第二,无法利用模拟滤波器的设计结果,一般无解析公式,要借助计算机辅助设计程序完成。
3设计和仿真
3.1 IIR数字滤波器的MATLAB 设计
IIR 数字滤波器的设计方法有:巴特沃斯设计法,切比雪夫设计法,椭圆法,贝塞尔法等,但根据自己所了解的滤波器设计不同方法之间的不同特点,我这里选择椭圆法来进行IIR 滤波器的设计。在设计IIR 滤波器中有脉冲冲激响应法和双线性不变法两种设计方法,在设计中采用脉冲冲激响应法进行设计。下面是我采用的椭圆法设计的MATLAB 源程序:设计指标为wp1=5000Hz,wp2=6000Hz,ws1=4800Hz,ws2=6200Hz,带通通带两边第一下降沿的第一旁瓣的下降幅度为60dB, 因此选用椭圆法进行设计。
wp1=5000*2*pi;wp2=6000*2*pi;ws1=4800*2*pi;ws2=6200*2*pi;
rp=0.5;rs=60;Fs=20000;
wp=[wp1 wp2];ws=[ws1 ws2];Nn=128;
[N Wc]=ellipord(wp,ws,rp,rs,'s');
[z,p,k]=ellipap(N,rp,rs);
[b,a]=zp2tf(z,p,k);
wo=sqrt(wp1*wp2);Bw=wp2-wp1;
[b1,a1]=lp2bp(b,a,wo,Bw);
[bz az]=impinvar(b1,a1,Fs);
figure;
freqz(bz,az);
grid on;
figure;
zplane(bz,az);
xlabel('实部');ylabel('虚部');grid on;
基于MATLAB 窄带带通滤波的设计
figure;
grpdelay(bz,az,Nn);
xlabel('归一化频率');ylabel('群延迟/采样数');
[H W]=freqz(bz,az);
figure;
plot(W*Fs/(2*pi),abs(H));
xlabel('频率/Hz');ylabel('幅值');grid on
运行源程序得到幅度相位响应图,零极点图,群延迟图和幅值与频率的关系图分别如下:
图3-1 IIR滤波器幅度相位响应关系图
由图3-1可知,采用设椭圆法设计IIR 滤波器反应了IIR 滤波器设计中的一些特点,这些并不是设计方法的问题,而是IIR 滤波器的特性,方法的优化只能让性能接近预期值。看图可知,IIR 滤波器相频特性并不具有线性关系,幅度和相位在频率边界有较大的起伏波动。
巢湖学院2014届本科毕业论文(设计)
图3-2 IIR滤波器零极点图
从图3-2可知,设计的IIR 滤波器极点很接近单位圆,零点则有的在圆心,有的则在单位圆附近。因此设计中的指标的改变,例如带通通带两边第一下降沿的第一旁瓣的下降幅度增大时,可能极点会移到单位圆外,就会出现稳定性问题。
图3-3 IIR滤波器群延迟图
由图3-3可知,在设计所得的IIR 滤波器群延迟图中,群延迟在带通边界
基于MATLAB 窄带带通滤波的设计
频率的范围内出现很大的尖峰和波动,这与相频响应图相对应,反应IIR 滤波器泄露效应。
图3-4 IIR滤波器幅值与频率的关系图
由图3-4可知,利用椭圆法设计IIR 滤波器的幅频响应图,图显示了设计的IIR 滤波器达到了设定的设计指标,通带波动幅度rp 和阻带最小下降rs 都达到了指标。
3. 2 IIR数字滤波器的simulink 仿真
利用MATLAB 中simulink 模块进行IIR 滤波器设计及滤波性能的仿真,仿真设计中用到了正弦波信号源SineWave, 信号混合器add ,数字滤波器Fdatool 设计工具箱(DigitalFilter Design),显示器scope 等。其中三个SineWave 信号源的 频率设置分别为1kHz,5.5kHz,7kHz, 其中一个设置如图:
基于MATLAB 窄带带通滤波的设计
图3-5 SineWave信号源频率设置图
数字滤波器fdatool 工具箱设置如下:
图3-6 IIR数字滤波器fdatool 设置图
基于MATLAB 窄带带通滤波的设计
从利用fdatool 工具箱中设定指标参数进而设计的带通滤波器中发现,采用IIR 滤波器设计方法中的椭圆法设计带通滤波器,所需的最小阶数为16,可见设计IIR 滤波器所需的阶数是较小的。
整个设计仿真模型的搭建如图:
图3-7 IIR滤波器仿真模型的搭建图
设计仿真运行的合成波和经过滤波的波形如图:
图3-8 过滤波与混合波对比图(IIR )
从上图中可以看出,三种信号波形成混合波和经过滤波器滤波的波形相比,通过使用fdatool 工具箱设定设计参数指标来进行设计的IIR 滤波器进行滤波,7kHz 和1kHz 的信号波形已被滤除,仿真设计实现了预定的功能。
基于MATLAB 窄带带通滤波的设计
3.3 FIR数字滤波器的MATLAB 设计
相比与IIR 滤波器设计FIR 数字滤波器设计的方法,显得复杂,无序没有特定程序化的方法,有窗函数法,频率采样法,最优化设计法等设计方法,为了使设计的滤波器性能最优,以及设计的程序的方便性,采用最优化的FIR 滤波器设计方法。下面是我采用的最优化设计方案的MATLAB 源程序:设计指标为wp1=5000Hz,wp2=6000Hz,ws1=4800Hz,ws2=6200Hz,带通通带两边第一下降沿的第一旁瓣幅度为60dB ,因此选用FIR 滤波器设计法最优化的kaiser 窗函数法。
Fs=20000;wp1=5000;wp2=6000;ws1=4800;ws2=6200;
f=[4800,5000,6000,6200];rp=0.5;rs=60;
A=[0,1,0];
a1=(10^(rp/20)-1)/(10^(rp/20));
b1=10^(-rs/20);
dev=[b1,a1,b1];
[n,wn,beta,type]=kaiserord(f,A,dev,Fs);
n=n+rem(n,2);
b=fir1(n,wn,type,kaiser(n+1,beta),'noscale');
figure;
freqz(b,1,1024);
figure;
[h,w]=freqz(b);
plot(w*Fs/(2*pi),abs(h));
title('带通滤波器的幅频响应');
xlabel('频率');ylabel('幅值');grid on;
figure;
impz(b);grid on;
figure;zplane(b,1);
xlabel('实部');
ylabel('虚部');
title('带通滤波器的零极点图');
figure;grpdelay(b,1,512);
xlabel('归一化频率');ylabel('群延迟/采样数');
title('带通滤波器的群延迟');
根据设计的指标设计的源程序,运行源程序得到对应的幅度相频响应,零极点图,群延迟图分别如下:
基于MATLAB 窄带带通滤波的设计
图3-9 FIR幅度和相位响应图
有图3-9可知,利用FIR 滤波器设计方法中的最优化kaiser 窗函数法进行设计的滤波器,上图的设计结果反映了相比IIR 滤波器,FIR 滤波器具有线性相位特性,避免了相位失真的现象。在滤波器频率边界下降旁瓣较小,也可以看出,设计的FIR 滤波器所需的滤波器阶数较大。
图3-10 FIR滤波器零极点图
基于MATLAB 窄带带通滤波的设计
从图3-10可知,设计的FIR 滤波器的极点全在原点,且极点的阶数为364,零点在单位圆附近,因此从设计的结果可以得出结论,FIR 滤波器不存在稳定性的问题,无论设计的指标为多少,FIR 滤波器都是稳定的。 o
图3-11 FIR滤波器的群延迟图
由图3-11可知,设计的FIR 滤波器群延迟中,群延迟与频率成以水平直线关系,这一关系与相频响应一一对应,充分的反映了FIR 滤波器的严格的线性相位关系。
图3-12 FIR滤波器幅值与频率关系图
基于MATLAB 窄带带通滤波的设计
由图3-12可知,利用最优化的kaiser 窗函数法设计的FIR 滤波器幅度与频率的关系图与设计理论图相比,设计达到了设计的预期指标,从图中知,设计的通带特性成理想的特性,阻带内也达到接近理想的状态。
3.4 FIR数字滤波器的simulink 仿真
同样,在MATLAB 中利用simulink 模块设计FIR 数字滤波器和进行滤波器性能的仿真,在仿真中也用的了正弦波信号源SineWave ,信号混合器add ,数字滤波器设计工具fdatool (DigitalFilter Design), 显示器scope 等模型,三个信号源SineWave 的频率设置也是1kHz,5.5kHz,7kHz, 不同的是数字滤波器fdatool 工具箱的设置,选择采用kaiser 窗,最小阶数,采样频率20kHz.ws1=4.8kHz,wp1=5kHz,wp2=6kHz,ws2=6.2kHz,最小下降边沿幅度为60dB. 其数字滤波器fdatool 参数设置如图:
图3-13 FIR数字滤波器fdatool 设置图
从利用fdatool 工具箱中设定指标参数进而设计的带通滤波器中发现,利用kaiser 窗函数法设计的最小阶数为363,而在设计中所选用的方法是FIR 滤波器设计法中所需阶数较小的之一。
整个FIR 滤波器设计和仿真模型的搭建如图:
基于MATLAB 窄带带通滤波的设计
图3-14 FIR滤波器仿真模型的搭建图
模型仿真运行的混合波和经过滤波器后的波形如下:
图3-15 过滤波与混合波对比图(FIR )
从仿真设计模型生成的混合波和经过滤波器的波形对比可以看出,经过滤波器过滤的信号已经把7KHz 和1kHz 的信号波过滤掉了,
因此仿真设计成功的实现
基于MATLAB 窄带带通滤波的设计
了预定的功能。
4. 设计和仿真结果分析
上面进行的设计主要研究了IIR 滤波器和FIR 滤波器两种滤波器的设计方法的不同优缺点,并利用了两种滤波器中设计比较好的方法进行设计,从而发现它们的优缺点的不同,为人们的使用提供指导。
在前面我分别设计研究了IIR 滤波器和FIR 滤波器,发现IIR 滤波器设计出的幅度相频响应,不具有线性相位,相位随着频率的改变而不成线性变化,极点是有可能随着参数的改变而移出单位圆的,群延迟是在过渡带附近有很大的波动,不过其所需的滤波器阶数较小。相比IIR 滤波器,FIR 滤波器设计所得到的各种响应图,从中发现FIR 滤波器具有严格的线性,相位随着频率的改变成线性变化,极点均在原点,不存在稳定问题,群延迟一直保持不变。
从上面的设计所得到结果,可以发现,IIR 滤波器相比FIR 滤波器,在设计相同的性能指标,IIR 滤波器所需的阶数较小,有程序化的设计方法和公式,但是IIR 滤波器不具有线性相位,可能会出现稳定性问题,相位延迟出现较大变化。而FIR 滤波器具有线性相位,不存在稳定性问题,相位延迟保持不变,这对应用中所需要相位要求严格的可以使用FIR 滤波器,但是FIR 滤波器也有缺点,在设计相同性能指标的情况下,FIR 滤波器所需的阶数较大。如设计中的IIR 滤波器所需阶数仅需16阶,而FIR 滤波器所需的阶数则需363阶。
基于MATLAB 窄带带通滤波的设计
参考文献
[1] 王艳芬,王刚,张晓光,刘卫东. 数字信号处理原理及实现[M].北京: 清华大学出版社, 2008:131-210.
[2] 唐向宏,岳恒立,郑雪峰. MATLAB及在电子信息类课程中的应用[M]. 北京: 电子工业出版社, 2011:15 -178.
[3] 万永革]\.数字信号处理的MATLAB 实现[M]. 北京:科学出版社,2007,4:110-234.
[4] 张明照,刘政波,刘斌等. 应用MATLAB 实现信号分析和处理[M]. 北京:科学出版社,2006:176-235.
[5] 吴大正. 信号与线性系统分析[M].北京: 高等教育出版社, 2005:115-312.
[6] 王世一. 数字信号处理[M].北京: 北京理工大学出版社,2011:123-268.
[7] 郝 娟, 徐沛文. MATLAB 在FIR 数字带通滤波器设计中的应用[J].应用技术与研究.
[8] 张学敏. FIR 低通和带通滤波器的关系分析与仿真[J].现代电子技术, 2008,19(282).
[9] 王艳文, 杨 楠. 基于matlab 窗函数的FIR 带通滤波器设计[J]. 信息产业 (141).
[10] 凌生强, 廖柏林等. 基于matlab 的FIR 带通滤波器设计及DSP 实现[J].现代电子技术,2012,35(5).