【毕业设计】基于matlab的自适应滤波器设计管理资料

1、基于MATLAB的自适应滤波器设计(吉首大学物理科学与信息工程学院，湖南吉首416000)张海摘要关键词：MATLAB；自适应；滤波器BasedonMATLABauto-adaptedfilterdesign(CollegeofPhysicsScienceandInformationEngineering,JishouUniversity,Jishou,Hunan416000)ZhanghaiAbstractKeyword：Matlab；Auto-adapted；Filter第一章绪论引言滤波器根据其逼近函数的形式不同,可设计出多种滤波器.常用的有巴特沃思滤波器、切比雪夫总型滤波器、切比雪夫E

2、型滤波器、椭圆滤波器、巴塞尔滤波器。对于这些滤波器的设计，都是先给定其副频特性的模平方uH(jX)u2,再求出系统函数H(s)。设计滤波器时，需由经典式求出滤波器的系统函数H(s),求出极点Sk(k=1,2,2N)，给定N,Xc,E,即可求得2N个极点分布。然后利用归一化函数,得出归一化的电路组件值,即可得到满足要求的滤波器。此种设计中,需要进行烦琐、冗长的数字计算,这对于电路设计者来说,不仅费时费力,准确性不易把握,而且不符合当今高速发展的时代要求。自适应滤波器是近30年来发展起来的关于信号处理方法和技术的滤波器,其设计方法对滤波器的性能影响很大。自适应滤波器能够得到比较好的滤波性能,当输入

3、信号的统计特性未知,或者输入信号的统计特性变化时,自适应滤波器能够自动地迭代调节自身的滤波器参数,以满足某种准则的要求,从而实现最优滤波。自适应滤波器一般包括滤波器结构和自适应算法两个部分,这两部分不同的变化与结合,可以导出许多种不同形式的自适应滤波器。MATLAB简介Matlab是由美国MathWorks公司推出的软件产品。它是一完整的并可扩展的计算机环境,是一种进行科学和工程计算的交互式程序语言。它的基本数据单元是不需要指定维数的矩阵,它可直接用于表达数学的算式和技术概念,而普通的高级语言只能对一个个具体的数据单元进行操作。在Matlab内部还配备了涉及到自动控制、信号处理、计算机仿真等种

4、类繁多的工具箱，所以Matlab的应用非常广泛,它可涉足于数值分析、控制、信号分析、通信等多种领域。自适应滤波器的应用适应滤波器在信号检测、信号恢复、数字通信等许躲领域中被广泛应用，因而一直是学术界一个重要研究课题。近年来，微电子技术和超大规模集成(VLSI)电路技术的飞速发展又促进了自适应滤波技术的进步。自适应滤波技术正是由于具有对干扰频率不敏感且其权值调整是基于对系统参数的优化等特点，而越来越多地受到人们的关注。传统的自适应滤波器主要在时域中实现，采用抽头延迟线(TappedDelayLine)结构及Widrow-Hoff自适应LMS算法。这种方法算法简单，稳健性也比较好,因而被广泛的应用

5、。但是滤波器的阶数可能会很高，步长系数可能会很小，收敛性能不理想，对输入信号的自相关矩阵有很强的依赖性，因而不具有高自适应率。当输入信号的自相关矩阵的特征值分布发散度很大时，算法的收敛速度很慢，跟踪性能不好。第二章自适应滤波原理自适应滤波器所谓自适应滤波，就是利用前一时刻已获得的滤波器参数等结果，自动地调节现时刻的滤波器参数，已适应信号和噪声未知的或随时间变化的统计特性，从而实现最优滤波。自适应滤波器由参数可调的数字滤波器(或称为自适应处理器)和自适应算法两部分组成，如图1所示：x(n)参数可调y(n)数字滤波器1/11L_自适应算法小汽)4图1自适应滤波器输入信号x(n)通过参数可调的数字滤

6、波器后产生输出信号(或响应)y(n),将其与参数信号(或称期望信号)d(n)进行比较，形成误差信号e(n)。e(n)(有时还要利用x(n),通过某种自适应算法对滤波器进行参数调整，最终使e(n)的均方误差最小。因此，自适应滤波器实际上是一种能够自动调整本身参数的特殊维纳滤波器，在设计时不需要事先知道关于输入信号和噪声的统计特性知识，它能够在自己工作过程中逐渐了解或估计所需的统计特性，并一次为根据自动调整自己的参数，以达到最佳滤波效果。一旦输入信号统计特性发生变化，它又能跟踪这种变化，自动调整参数，使滤波器性能达到最佳。自适应算法自适应算法主要是根据滤波器输入的统计特性进行处理，他可能还与滤波器

7、的输入及其它数据有关，据此，存在开环算法和闭环算法。开环算法的控制输出仅取决与滤波器的输入和其它输入数据；闭环的控制输出则是滤波器输出及其它输入信号的函数。闭环控制利用输出反馈，它不但能在滤波器输入信号变化时保持最佳输出，且还能在某种程度上补偿滤波原件参数的变化和误差及运算误差。缺点：存在稳定性问题及收敛度不高；开环算法的优点是调整速度快，一般无稳定性问题，但通常它要求的计算量大且不能补偿组件参数误差及运算误差，所以多用闭环算法。图2开环算法图3闭环算法自适应滤波去噪原理一种自适应去噪滤波器原理如图1所示，信号X(k)受到宽带噪声n(k)的污1染，与噪声相关的宽带信号n(k)是可以测量的。其中

8、n(k)是和n1(k)彼此相关22的噪声信号，而与信号x(k)不相关。如果将n(k)作为自适应滤波器的输入，而将2x(k)作为期望输出信号，通过对输出误差e(k)的控制，可以调整自适应滤波器的权系数w(k),权系数的更新使得它的输出y(k)趋于等于n1(k),则当滤波器稳定以后，自适应滤波器输出误差e(k)就是滤除了n1(k)的期望输出信号x(k)。+e(k)Y(k)X(k)+N(k)%（k）自适应滤波图4自适应去噪滤波原理图第三章自适应滤波器的LMS算法LMS算法最小均方算法(LMS-LeastMeanSquare)是一种很有用且很简单的估计梯度的方法，其突出的特点是计算量小，易于实现，且不

9、要求脱线计算。最陡下降法每次迭代都需要知道性能曲面上某点的梯度值，而实际上梯度值只能根据观测数据进行估计，LMS算法是一种很有用且很简单的估计梯度的方法。其关键技术在于按照e(n)及各x(n)值，通过某种算法，确定Ee2(n)=min时的各h*(n)值，从而自动调节各h(n)值至h*(n)值。LMS算法最核心的思想是用平方误差代替均方误差，即V(n)V(n)=Le2(n)/dwde2(n)/dwde2(n)/dw-T01LV=2e(n)=-2e(n)x(n)5w将上式带如最陡下降法迭代计算权矢量的公式，得w(n+1)=w(n)+卩(-V(n)LMS算法的基本关系式w(n+1)=w(n)+2卩e

10、(n)x(n)其中卩是收敛因子，决定了收敛速度和稳定性，卩满足：0卩1九max，九max是R的最大特征值，R=EX(n)Xt(n)5XXXX许多学者对LMS算法进行了研究，提出了传统LMS算法的许多有效的改进措施：如采用变步长LMS算法、变换域LMS算法、QR分解LMS算法等，有效地克服了其性能局限性。为了比较直观地观察和分析各种LMS算法的收敛性能，借助MATLAB强大的工程计算和绘图功能，,用计算机仿真,对输入信号做相应的处理,并分析仿真结果。仿真结果中收敛曲线均是采用蒙特卡罗仿真，独立运行100次求其统计平均得到的。收敛曲线的横轴均为滤波器迭代次数。时域LMS算法传统的LMS算法具有计算

11、量小，结构简单，易于实现等诸多优点，尤其是这种算法是最先由统计分析法导出的一种实用算法，它是一类自适应滤波器的基础。所以，详细分析时域LMS算法中个参数对算法的影响具有重要意义。下面就针对时域LMS算法各参数做一下讨论。步长U0仿真抽样点数为N=512,滤波器阶数K=&单频信号为s=a*esp(i*pi*t),a=l,加入均值为零的高斯白噪声后信噪比为snr=10dB。由仿真结果可看出：当U二，LMS算法不能收敛，U二，算法收敛较慢，不能跟踪输入信号的变化，而当U=，既能快速达到收敛，有较小的起始阶段误差，并达到较小的稳态误差。步长U的大小决定着算法的收敛速度和达到稳态的失调量的大小，对于常数

12、U来说，算法的失调与自适应收敛过程是一对矛盾，要想得到较快的收敛速度可选用大的U值，这将导致较大的失调量(如果要满足失调量的要求，则收敛过程将受到制约。这里，仿真结果与理论值相一致。下面介绍的变步长LMS算法能有效解决这一矛盾。3.2.2阶数K。仿真抽样点数为M=512,收敛步长U二二a*esp(i*pi*t),a=l,加入均值为零的高斯白噪声，信噪比为snr=10Db，做出k=6，&10,12时的收敛曲线。该曲线表明：对于不同的滤波器阶数可得到不同的滤波效果，当K=8时稳态误差最小，信号输出波形最好。这是因为LMS滤波器阶数K与稳态误差及输入信号特性有关，对于具体的输入信号，有一个最佳(或准

13、最佳)的加权数目K使稳态误差最小，再增加权数目时，稳态误差有变大的可能.3.2.3信噪比snr。仿真抽样点数为N=512，步长U二二a*esp(i*pi*t),a=1,阶数k=8。由计算机仿真图看出，当信噪比snr升高时,LMS算法性能将急剧恶化。这可由频域LMS算法来克服时域LMS算法的性能局限。自适应滤波由于具有对干扰频率不敏感且其权值调整是基于对系统参数的优化等特点，广泛地应用于信号检测、信号恢复、数字通信等领域。传统自适应滤波器主要在时域实现，该算法简单，稳健性能较好，因而被广泛应用。第四章基于LMS算法的自适应滤波器的仿真及性能分析原理图原始语音采用文件“audio.wav”，噪声采

14、用一组频率为50Hz,500Hz和1000Hz的正弦信号模拟,其结构图如图4所示。FromWaveFileAudio.Wav8000Hz/1Ch/16b原始语音信号噪声信号原始语音信号噪声信号混合语音信号去噪后语音信号FromWaveFilenoise.Wav8000Hz/lCh/16b噪声信号InOutRLSErrTapeRLSAdaptiveFilter去噪后语音信导Results混合后语音信号FFTFFT1FFTFFT去噪后语音的频谱混合后语音的频谱原始语音信号Results1图4语音去噪结构图自适应滤波源程序程序清单：Functiony,e=1ms(x,d,u,N,e)%y,e=1ms

15、(x,d,u,N,e)%u=2*收敛因子%w=估计的FIR滤波器%y=输出数组y(n)%x=输入数组x(n)%d=期望数组d(n)，长度与x相同%u=步长%N=FIR滤波器的长度%M=length(x)；Y=zeros(1,M)；W=zeros(1,N)ezeros(1,M)；forn=N:M；x1=s(n:-1:n-N+1);y(n)=h*x1;e(n)=d(n)-y(n);w(n)=w+u*e(n)*x1;end源程序：Functionlmsmain()%主函数%XS是周期性信号%XN是干扰噪声%T=0:;Xs=10*sin(*t);Figure(1);Subplot(211);Plot(

16、t,xs);grid;Ylabei(幅值)；Title(i输入周期性信号);Xn=randn(1,1000)；Subplot(212);Plot(t,xn);grid;Ylabel(幅值)；Xlabel(时间)；Title(it随即噪声信号);d=xs+xn；x=xs+xn;u=2e-3;N=10;W=0000000000;ye=lms(x,d,u,N,e);Figure(2);%subplot(313)Plot(t,d,t,y,t,e);grid;Xlabel(时间)；Ylabel(幅值)；波形仿真部分源程序如下:%设置滤波器参数M=32；lam=1；delta=0.1；w0=zeros(M

17、,1)；P0=(1/delta)3eye(M,M)；Zi=zeros(M-1,1)；%RLS自适应滤波S=initrls(w0,P0,lam,Zi)；y,e=adaptrls(noise,signal,S)；其仿真过程中的图形如下:T,|ijfIi.MPPI_*4.JB图5原始语音信号的波形图6引入噪音的波形图7混合后的语音信号图8混合后信号FFT的图形图9滤波后信号FFT的图形图10滤波前后信号的比较由仿真过程可以看出，通过图7和图8可以清晰的看到在50,500,1000处的尖峰被滤除了。再由图10比较可以发现除在开始时信号有很小一段延时外,在其它时间都能很好与原始语音信号重合。性能分析通过

18、反复试验,使用Matlab能很好的解决各种噪声混合后的语音信号去噪问题。当在现场使用时，只要保证两个麦克风同步，所拾取的语音从同一时间端点开始进行处理，就能得到清晰的语音信号。第五章结论自适应滤波由于具有对干扰频率不敏感且其权值调整是基于对系统参数的优化等特点，广泛地应用于信号检测、信号恢复、数字通信等领域。传统自适应滤波器主要在时域实现，该算法简单，稳健性能较好，因而被广泛应用。随着人们对现代通信过程中的信号质量要求越来越高,自适应滤波器的算法越来越复杂。本文通过自适应去噪滤波器的设计,为自适应滤波器的硬件实现提供了一种切实可行的实验依据。变换域自适应滤波获得了比时域自适应滤波更好的收敛性能

19、。目前已有多种形式的正交策略，如KLT、离散Fourier变换、离散余弦变换(DCT)等。近年来又兴起了被称作“数学显微镜”的小波变换(DWT),可利用它的多分辨和时频局部化特性在不同的分辨水平上重现信号。将小波变换引入自适应滤波结构，具有小波分析优于Fourier分析的所有优点，在时变信号和快速变化信号的自适应滤波方面有着广阔的应用前景，是未来自适应滤波发展的新方向。利用Matlab来设计要比传统设计步骤简便、直观的多.在实际应用中，应根据需要来选择是采用分步设计,还是采用一次性设计,不论是哪一种,与传统设计相比，计算量都要减少得多。Matlab作为一种交互式的程序语言，在滤波的应用中有着自己独特的优点.运用给定的信号处理工具箱函数,编制相应的程序，可实现各种形式的模拟、数字滤波器，这对于信号的