下载

1下载券

加入VIP
  • 专属下载特权
  • 现金文档折扣购买
  • VIP免费专区
  • 千万文档免费下载

上传资料

关闭

关闭

关闭

封号提示

内容

首页 基于matlab的FIR数字滤波器设计~

基于matlab的FIR数字滤波器设计~.doc

基于matlab的FIR数字滤波器设计~

tyme
2013-04-04 0人阅读 举报 0 0 暂无简介

简介:本文档为《基于matlab的FIR数字滤波器设计~doc》,可适用于IT/计算机领域

基于MATLAB的FIR数字滤波器设计通信专业课程设计一太原科技大学课程设计(论文)设计(论文)题目:基于MATLAB的FIR数字滤波器设计计姓名xxxx学号班级通信学院电子信息工程指导教师xxxx年月日太原科技大学课程设计(论文)任务书学院(直属系):电子信息工程学院时间:年月日学生姓名Xxxx指导教师Xxxx设计(论文)题目基于MATLAB的FIR数字滤波器设计主要研究内容针对目前常用的窗函数法、频率采样法和最优化设计法等FIR滤波器的设计方法进行比较分析研究它们各自的优缺点及适用对象。研究方法窗函数法、频率采样法和最优化设计法主要技术指标(或研究目标)passrad=*piw=boxcar()w=hamming()n=::教研室意见教研室主任(专业负责人)签字:年月日目录摘要第章绪论数字滤波器的研究背景和意义数字滤波器的应用及现状数字滤波器的设计方法第章FIR滤波器FIR滤波器的基本概念FIR滤波器的特点FIR滤波器的种类第章FIR数字滤波器的设计窗函数法频率采样法最优化设计等波纹切比雪夫逼近准则仿真函数第章结果分析参考文献附录基于MATLAB的FIR数字滤波器设计摘要论文针对窗函数法、频率采样法和最优化设计法三种设计方法采用MATLAB进行FIR滤波器的设计与仿真并比较了三种方法的特点。结果表明在同样的设计指标下利用等波纹切比雪夫逼近法则的设计可以获得最佳的频率特性和衰耗特性具有通带和阻带平坦过渡带窄等优点。关键词:FIR滤波器MATLAB窗函数法频率采样法等波纹切比雪夫逼近第章绪论数字滤波器的研究背景和意义当今数字信号处理技术正飞速发展它不但自成一门学科更是以不同形式影响和渗透到其他学科它与国民经济息息相关与国防建设紧密相连它影响或改变着我们的生产、生活方式因此受到人们普遍的关注。数字化、智能化和网络化是当代信息技术发展的大趋势而数字化是智能化和网络化的基础实际生活中遇到的信号多种多样例如广播信号、电视信号、雷达信号、通信信号、导航信号等等。上述这些信号大部分是模拟信号也有小部分是数字信号。模拟信号是自变量的连续函数自变量可以是一维的也可以是二维或多维的。大多数情况下一维模拟信号的自变量是时间经过时间上的离散化(采样)和幅度上的离散化(量化)这类模拟信号便成为一维数字信号。因此数字信号实际上是数字序列表示的信号语音信号经采样和量化后得到的数字信号是一个一维离散时间序列而图像信号的经采样和量化后得到的是数字信号是一个二维离散空间序列。数字滤波技术是数字信号分析、处理技术的重要分支。无论是信号的获取、传输还是信号的处理和交换都离不开滤波技术它对信号安全可靠和有效灵活地传输时至关重要的。在所有的电子系统中使用最多技术最复杂的要算数字滤波器了。数字滤波器的优劣直接决定产品的优劣。数字滤波器的应用及现状数字滤波器精确度高、使用灵活、可靠性高具有模拟设备所没有的许多优点已广泛应用于各个学科技术领域例如数字电视、语音、通信、雷达、声纳、遥感、图像、生物医学以及许多工程应用领域。随着信息时代数字时代的到来数字滤波技术已经成为一门极其重要的学科和技术领域。以往的滤波器大多采用模拟电路技术但是模拟电路技术存在很多难以解决的问题例如模拟电路元件对温度的敏感性等等。而采用数字技术则避免很多类似的难题当然数字滤波器在其他方面也有很多突出的优点所以采用数字滤波器对信号进行处理是目前的发展方向。数字滤波器的设计方法数字滤波器按照单位取样响应的时域特性可以分为无限脉冲响应(IIR)系统和有限脉冲响应(FIR)系统。FIR数字滤波器的优点在于它可以做成具有严格线性相位,而同时可以具有任意的幅度特性它的传递函数没有极点这保证了设计出的FIR数字滤波器一定是平稳的。所谓数字滤波器设计简单地说就是要找到一组能满足特定滤波要求的系数向量a和b。而滤波器设计完成后还需要进一步考虑如何将其实现即选择什么样的滤波器结构来完成滤波运算。FIR数字滤波器的设计方法很多其中较为常用的是窗函数设计法、频率采样设计法和最优化设计法。本文讨论利用窗函数法、频率采样法和等波纹切比雪夫逼近法(调用remez函数)来分别实现各种FIR滤波器的设计。窗函数法设计的基本思想是把给定的频率响应通过IDTFT(InverseDiscreteTimeFourierTransform)求得脉冲响应然后利用加窗函数对它进行截断和平滑以实现一个物理可实现且具有线性相位的FIR数字滤波器的设计目的。其核心是从给定的频率特性通过加窗确定有限长单位取样响应频率采样法设计的基本思想是把给出的理想频率响应进行取样通过IDFT从频谱样点直接求得有限脉冲响应最优化设计方法是指采用最优化准则来设计的方法在FIRDF的最优化设计中,最优化准则有均方误差最小化准则和等波纹切比雪夫逼近(也称最大误差最小化)准则两种。第章FIR滤波器FIR滤波器的基本概念FIR滤波器:有限长单位冲激响应滤波器是数字信号处理系统中最基本的元件它可以在保证任意幅频特性的同时具有严格的线性相频特性同时其单位抽样响应是有限长的因而滤波器是稳定的系统。因此FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。滤波器设计是根据给定滤波器的频率特性求得满足该特性的传输函数。FIR滤波器的特点有限长单位冲激响应(FIR)滤波器有以下特点:  ()系统的单位冲激响应在有限个n值处不为零  ()系统函数在处收敛极点全部在z=处(因果系统)  ()结构上主要是非递归结构没有输出到输入的反馈但有些结构中(例如频率抽样结构)也包含有反馈的递归部分。  设FIR滤波器的单位冲激响应为一个N点序列则滤波器的系统函数为()就是说它有(N)阶极点在z=处有(N)个零点位于有限z平面的任何位置。优点:()很容易获得严格的线性相位避免被处理的信号产生相位失真这一特点在宽频带信号处理、阵列信号处理、数据传输等系统中非常重要()可得到多带幅频特性()极点全部在原点(永远稳定)无稳定性问题()任何一个非因果的有限长序列总可以通过一定的延时转变为因果序列所以因果性总是满足()无反馈运算运算误差小。缺点:()因为无极点要获得好的过渡带特性需以较高的阶数为代价()无法利用模拟滤波器的设计结果一般无解析设计公式要借助计算机辅助设计程序完成。FIR滤波器的种类()数字集成电路FIR滤波器一种是使用单片通用数字滤波器集成电路这种电路使用简单但是由于字长和阶数的规格较少不易完全满足实际需要。虽然可采用多片扩展来满足要求但会增加体积和功耗因而在实际应用中受到限制。设计数字滤波器的任务就是寻求一个因果稳定的线性时不变系统使其系统函数H(z)具有指定的频率特性。()DSP芯片FIR滤波器另一种是使用DSP芯片。DSP芯片有专用的数字信号处理函数可调用实现FIR滤波器相对简单但是由于程序顺序执行速度受到限制。而且就是同一公司的不同系统的DSP芯片其编程指令也会有所不同开发周期较长。图DSPFIR滤波器()可编程FIR滤波器还有一种是使用可编程逻辑器件FPGA/CPLD。FPGA有着规整的内部逻辑块整列和丰富的连线资源特别适合用于细粒度和高并行度结构的FIR滤波器的实现相对于串行运算主导的通用DSP芯片来说并行性和可扩展性都更好。第章FIR数字滤波器设计FIR滤波器设计的任务是选择有限长度的使传输函数满足一定的幅度特性和线性相位要求。由于FIR滤波器很容易实现严格的线性相位所以FIR数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。设计过程一般包括以下三个基本问题:()根据实际要求确定数字滤波器性能指标()用一个因果稳定的系统函数去逼近这个理想性能指标()用一个有限精度的运算去实现这个传输函数。窗函数法设计FIR数字滤波器的最简单的方法是窗函数法通常也称之为傅立叶级数法。FIR数字滤波器的设计首先给出要求的理想滤波器的频率响应,设计一个FIR数字滤波器频率响应去逼近理想的滤波响应。然而窗函数法设计FIR数字滤波器是在时域进行的因而必须由理想的频率响应推导出对应的单位取样响应再设计一个FIR数字滤波器的单位取样响应去逼近。设计过程如下:()加窗的作用是通过把理想滤波器的无限长脉冲响应乘以窗函数来产生一个被截断的脉冲响应即并且对频率响应进行平滑。MATLAB工具箱提供的窗函数有:矩形窗(Rectangularwindow)、三角窗(Triangularwindow)、布拉克曼窗(Blackmanwindow)、汉宁窗(Hanningwindow)、海明窗(Hammingwindow)、凯塞窗(Kaiserwindow)、切比雪夫窗(Chebyshevwindow)。窗函数主要用来减少序列因截断而产生的Gibbs效应。但当这个窗函数为矩形时得到的FIR滤波器幅频响应会有明显的Gibbs效应并且任意增加窗函数的长度(即FIR滤波器的抽头数)Gibbs效应也不能得到改善。为了克服这种现象,窗函数应该使设计的滤波器:()频率特性的主瓣宽度应尽量窄,且尽可能将能量集中在主瓣内()窗函数频率特性的旁瓣ω趋于π的过程中其能量迅速减小为零。频率采样法频率采样法是从频域出发,根据频域采样定理对给定的理想滤波器的频率响应加以等间隔的抽样得到:k=,,…,N()再利用可求得FIR滤波器的系统函数及频率响应。而在各采样点间的频率响应则是其的加权内插函数延伸叠加的结果。但对于一个无限长的序列用频率采样法必然有一定的逼近误差,误差的大小取决于理想频响曲线的形状,理想频响特性变换越平缓,则内插函数值越接近理想值,误差越小。为了提高逼近的质量可以通过在频率相应的过渡带内插入比较连续的采样点扩展过渡带使其比较连续从而使得通带和阻带之间变换比较缓慢以达到减少逼近误差的目的。选取w∈,π内N个采样点的约束条件为:()()增大阻带衰减三种方法:)加宽过渡带宽以牺牲过渡带换取阻带衰减的增加。)过渡带的优化设计利用线性最优化的方法确定过渡带采样点的值得到要求的滤波器的最佳逼近(而不是盲目地设定一个过渡带值)。)增大N。如果要进一步增加阻带衰减但又不增加过渡带宽可增加采样点数N。代价是滤波器阶数增加运算量增加。直接从频域进行设计物理概念清楚直观方便适合于窄带滤波器设计这时频率响应只有少数几个非零值但是截止频率难以控制。典型应用:用一串窄带滤波器组成多卜勒雷达接收机覆盖不同的频段多卜勒频偏可反映被测目标的运动速度最优化设计最优化设计方法是指采用最优化准则来设计的方法。在FIRDF的最优化设计中,最优化准则有均方误差最小化准则和等波纹切比雪夫逼近(也称最大误差最小化)准则两种。实际设计中,只有采用窗函数法中的矩形窗才能满足前一种最优化准则,但由于吉布斯(Gibbs)效应的存在使其根本不能满足设计的要求。为了满足设计的要求,可以采用其它的窗函数来消除吉布斯效应,但此时的设计已经不能满足该最优化准则了。因此,要完成FIRDF的最优化设计,只能采用后一种优化准则来实现。 等波纹切比雪夫逼近准则在滤波器的设计中,通常情况下通带和阻带的误差要求是不一样的。等波纹切比雪夫逼近准则就是通过对通带和阻带使用不同的加权函数,实现在不同频段(通常指的是通带和阻带)的加权误差最大值相同,从而实现其最大误差在满足性能指标的条件下达到最小值。尽管窗函数法与频率采样法在FIR数滤波器的设计中有着广泛的应用,但两者不是最优化的设计。通常线性相位滤波在不同的频带内逼近的最大容许误差要求不同。等波纹切比雪夫逼近准则就是通过通带和阻带使用不同的加权函数,实现在不同频段(通常指的是通带和阻带)的加权误差最大值相同,从而实现其最大误差在满足性能指标的条件下达到最小值即使得和之间的最大绝对误差最小。等波纹切比雪夫逼近是采用加权逼近误差它可以表示为:()其中为逼近误差加权函数在误差要求高的频段上可以取较大的加权值否则应当取较小的加权值。尽管按照FIR数字滤波器单位取样响应h(n)的对称性和N的奇、偶性FIR数字滤波器可以分为种类型但滤波器的频率响应可以写成统一的形式:()其中k∈{,},H(ω)为幅度函数且是一个纯实数表达式也可以写成统一的形式:()其中为ω的固定函数为M个余弦函数的线性组合。仿真函数利用数字信号处理工具箱中的remezord和remez函数可以实现FIRDF的最优化设计。在此先介绍这两个函数:()n,fo,ao,weights=remezordf,a,dev功能:利用remezord函数可以通过估算得到滤波器的近似阶数n,归一化频率带边界fo,频带内幅值ao及各个频带内的加权系数weights。输入参数f为频带边缘频率,a为各个频带所期望的幅度值,dev是各个频带允许的最大波动。()h=remez(n,fo,ao,weights‘,ftype’)功能:利用remez函数可以得到最优化设计的FIRDF的系数,输入参数n是滤波器的阶数,fo,ao,weights参数含义说明同()。ftype是所设计的滤波器类型,它除了可以设计普通的滤波器外,它还可以设计数字希尔钞特变换器以及数字微分器。实际设计中,由于remezord函数可跑高估或低估滤波器的阶数n,因此在得到滤波器的系数后,必须检查其阻带最小衰减是否满足设计要求。如果此时的技术指标不能满足设计要求,则必须提高滤波器的阶数到n,n等。故等波纹切比雪夫逼近法设计FIR数字滤波器的步骤是:①给出所需的频率响应加权函数和滤波器的单位取样响应的长度N。②由①中给定的参数来形成所需的、和的表达式。③根据Remez算法,求解逼近问题。④利用傅立叶逆变换计算出单位取样响应。第章结果分析窗函数法仿真结果图窗函数设计的FIR低通滤波器频率响应结果分析:可以看出采用特殊的窗函数如Hamming窗可以减小Gibbs效应但同时也会使滤波器的过度带变宽。波动幅度取决于窗函数幅度频谱旁瓣的相对幅度,而波纹的多少取决于窗函数旁瓣的多少。以上两点是就是窗函数直接截断Hd(n)引起的截断效应在频域的反映,截断效应直接影响滤波器的性能,因为通带内的波动会影响滤波器痛带中的平稳性,阻带内的波动则影响阻带最小衰减,因此,减少截断效应也是FIR数字滤波器设计的关键之一。频率采样法图FIR的单位取样响应图FIR的低通衰减幅频特性结果分析:图为在间断点处增加一个过渡点后的情况。从图,中可以看出滤波器的带外衰减指标有了明显的改善!但这同时增加了滤波器的过渡带宽。所以在带外衰减和过渡带宽这两个指标之间需要有一个折衷。因频率取样点都局限在πN的整数倍点上所以在指定通带和阻带截止频率时这种方法受到限制比较死板。充分加大N可以接近任何给定的频率但计算量和复杂性增加。频率采样法偏离设计指标明显阻带衰减最小只有适当选取过渡带样点值才会取得较好的衰耗特性最优化设计图 滤波器输出的幅频及相频响应特性结果分析:在设计中,如果该滤波器的特性不满足要求,那么,原有参数必须作适当调整。这在程序中很容易实现,只需对参数进行重新设定,就可以得到新条件下滤波器的特性。采用最优化设计方法时大大减小了滤波器的阶数,从而减小了滤波器的体积,并最终降低了滤波器的成本。这样使得设计出来的滤波器更为简单经济。因而在实际的滤波器设计中,这种最优化方法是完全可行的。在实际应用中,如果需要对某一信号源进行特定的滤波,并要检验滤波效果,应用传统方法实施起来比较繁琐。在Matlab环境下,可先用软件模拟产生信号源,再设计滤波器对其进行滤波。同样是设计一个FIR低通数字滤波器综合分析可以看出:()窗函数法在阶数较低时阻带特性不满足设计要求只有当滤波器阶数较高时使用海明窗和凯塞窗基本可以达到阻带衰耗要求()频率采样法偏离设计指标最明显阻带衰减最小而且设计比采用窗函数法复杂。只有适当选取过渡带样点值才会取得较好的衰耗特性()利用等波纹切比雪夫逼近法则的设计可以获得最佳的频率特性和衰耗特性具有通带和阻带平坦过渡带窄等优点。综上所述FIR滤波器很容易实现具有严格线性相位的系统,使信号经过处理后不产生相位失真,舍入误差小,而且稳定因此越来越受到广泛的重视。MATLAB软件的诞生,使数字信号处理系统的分析与设计得简单,它已经成为电子工程师必备的一个工具软件。参考文献张亚妮基于MATLAB的数字滤波器设计J辽宁工程技术大学学报,,():饶志强,叶念渝FIR和IIR数字滤波器的探讨与实现J计算机与数字工程,,():孙强运用MATLAB实现数字滤波器的设计J电脑学习,,():附录窗函数实现程序:passrad=*piw=boxcar()w=hamming()n=::hd=sin(passrad*(n))(pi*(n))hd()=passradpih=hd*rot(w)h=hd*rot(w)mag,rad=freqz(h)mag,rad=freqz(h)subplot(,,)plot(rad,*log(abs(mag)))gridontitle('designedbyRectangularwindow')subplot(,,)plot(rad,*log(abs(mag)))gridontitle('designedbyHammingwindow')h,w=freqz(h,,,)subplot(,,)plot(w,unwrap(angle(h)))gridonh,w=freqz(h,,,)subplot(,,)plot(w,unwrap(angle(h)))gridon频率采样法实现程序:f=a=dev=给出滤波器的参数nfaw=remezord(f,a,dev)N=nalpha=(N)k=:Nwp=*piws=*pi计算理想低通滤波器的截止频率wc=(wpws)m=fix(wc*N(*pi))在两边过渡带取值为的采样点T=Hrs=ones(,m),T,zeros(,N*m),T,ones(,m)k=:floor(alpha)k=floor(alpha):Nphai=alpha*(*pi)N*k,alpha*(*pi)N*(Nk)H=Hrs*exp(j*phai)计算单位冲激响应h=ifft(H,N)h,w=freqz(h,,,)hr=abs(h)h=*log(hr)画出FIRDF的单位取样响应figure()k=:Nstem(k,h,'k')axis(,N,*min(real(h)),*max(real(h)))xlabel('n')ylabel('h(n)')gridon画出FIRDF的低通衰减幅频特性figure()plot(w,h)xlabel('Normalized')Frequency(×πradsample)')ylabel('Magnitude(dB)')gridon最优化设计实现程序fs=    采样频率rp=通带波纹rs=阻带波纹f=截止频率a=期望的幅度dev=(^(rp))(^(rp))^(rs)n,fo,ao,w=remezord(f,a,dev,fs)b=remez(n,fo,ao,w)调用最优设计法中remez函数reqz(b,,,fs)�EMBEDEquationDSMT����EMBEDEquationDSMT����EMBEDEquationDSMT����EMBEDEquationDSMT����EMBEDEquationDSMT����EMBEDEquationDSMT����EMBEDEquationDSMT���……y(n)(a)�EMBEDEquationDSMT���………………y(n)(b)………unknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknownunknown

用户评价(0)

关闭

新课改视野下建构高中语文教学实验成果报告(32KB)

抱歉,积分不足下载失败,请稍后再试!

提示

试读已结束,如需要继续阅读或者下载,敬请购买!

文档小程序码

使用微信“扫一扫”扫码寻找文档

1

打开微信

2

扫描小程序码

3

发布寻找信息

4

等待寻找结果

我知道了
评分:

/20

基于matlab的FIR数字滤波器设计~

VIP

在线
客服

免费
邮箱

爱问共享资料服务号

扫描关注领取更多福利