【正文】 1,:)*x(in+1:i)39。 w(i,:)=w(i1,:)+e(1,i)*K(i,:)。 p=pK(i,:)39。*x(in+1:i)*p。 y(1,i)=w(i,:)*x(in+1:i)39。 end figure(1) subplot(2,1,1)。plot(0:Nn,xr(1,n:N),39。r39。)ylabel(39。幅值39。)。xlabel(39。時間39。)。title(39。原始輸入信號xr時域圖39。) grid onsubplot(2,1,2)。plot(abs(fft(xr)))。grid。ylabel(39。幅值39。)。xlabel(39。頻率39。)。title(39。原始輸入信號xr頻域圖39。)。figure(2) subplot(2,1,1)。plot(0:Nn,x(1,n:N))title(39。濾波器輸入信號x時域圖39。) ylabel(39。幅值39。)。xlabel(39。時間39。)。gridsubplot(2,1,2)。plot(abs(fft(x)))。grid。ylabel(39。幅值39。)。xlabel(39。頻率39。)。title(39。濾波器輸入信號x頻域圖39。)。figure(3) plot(0:Nn,e(1,n:N))。hold on plot(0:Nn,y(1,n:N),39。r39。)。hold on grid on。title(39。誤差信號e(t)和輸出信號y(t)39。) 。ylabel(39。幅值39。)。xlabel(39。時間39。)。gtext(39。e(k)39。),gtext(39。y(k)39。) figure(4) subplot(2,1,1)。plot(abs(fft(e)))。grid。ylabel(39。幅值39。)。xlabel(39。頻率39。)。title(39。誤差信號e頻域圖39。)。subplot(2,1,2)。plot(abs(fft(y)))。grid。ylabel(39。幅值39。)。xlabel(39。頻率39。)。title(39。濾波器輸出信號y頻域圖39。)。 原始信號時域和頻域圖 濾波器輸誤差信號和輸出信號頻域圖 濾波器輸誤差信號和輸出信號頻域圖 濾波器輸誤差信號和輸出信號頻域圖如圖 從時域觀察，原始信號峰值為10，噪聲峰值為2。濾波之后誤差信號值接近于0。從頻域上看，接近原始信號頻譜。濾波效果較好。 輸入信號為非周期信號本實驗原始輸入信號為純背景音頻信號，格式為wav,將其放在matlab文件，存為“”。所加噪聲方差為2。實驗程序如下：clear。clc。randn(39。seed39。, 0) 。rand(39。seed39。, 0) 。 N=1000。n=200。 k=12。 b=2*randn(1,N)。 %純背景音樂的音頻信號F=500。c=16。[Y,F,c]=wavread(39。39。) 。Y1=Y(((F*10+1):F*12),:)。wavwrite(Y1,F,16,39。39。) 。wavplay(Y1,F)。Y2=(zeros(N,1))。for j=1:N Y2(j,1)=100*Y1(j,1)。endxr=Y239。%噪聲疊加x=(zeros(1,N))。for j=1:N x(1,j)=xr(1,j)+b(j)。end%RLS算法Delta=。p=(1/Delta)*eye(n)。K=zeros(N,n)。 w=zeros(N,n)。 e=zeros(1,N)。 y=zeros(1,N)。 tr=zeros(1,N)。 for i=n+1:N K(i,:)=(p*x(in+1:i)39。./(1+x(in+1:i)*p*x(in+1:i)39。))39。 e(1,i)=xr(i)w(i1,:)*x(in+1:i)39。 w(i,:)=w(i1,:)+e(1,i)*K(i,:)。 p=pK(i,:)39。*x(in+1:i)*p。 y(1,i)=w(i,:)*x(in+1:i)39。 end figure(1) subplot(2,1,1)。plot(0:Nn,xr(1,n:N),39。r39。)ylabel(39。幅值39。)。xlabel(39。時間39。)。title(39。原始輸入信號xr時域圖39。) grid onsubplot(2,1,2)。plot(abs(fft(xr)))。grid。ylabel(39。幅值39。)。xlabel(39。頻率39。)。title(39。原始輸入信號xr頻域圖39。)。figure(2) subplot(2,1,1)。plot(0:Nn,x(1,n:N))title(39。濾波器輸入信號x時域圖39。) ylabel(39。幅值39。)。xlabel(39。時間39。)。gridsubplot(2,1,2)。plot(abs(fft(x)))。grid。ylabel(39。幅值39。)。xlabel(39。頻率39。)。title(39。濾波器輸入信號x頻域圖39。)。figure(3) plot(0:Nn,e(1,n:N))。hold on plot(0:Nn,y(1,n:N),39。r39。)。hold on grid on。title(39。誤差信號e(t)和輸出信號y(t)39。) 。ylabel(39。幅值39。)。xlabel(39。時間39。)。gtext(39。e(k)39。),gtext(39。y(k)39。) figure(4) subplot(2,1,1)。plot(abs(fft(e)))。grid。ylabel(39。幅值39。)。xlabel(39。頻率39。)。title(39。誤差信號e頻域圖39。)。subplot(2,1,2)。plot(abs(fft(y)))。grid。ylabel(39。幅值39。)。xlabel(39。頻率39。)。title(39。濾波器輸出信號y頻域圖39。)。： 濾波器輸誤差信號和輸出信號頻域圖 濾波器輸誤差信號和輸出信號頻域圖 濾波器輸誤差信號和輸出信號頻域圖 濾波器輸誤差信號和輸出信號頻域圖如圖，從時域分析，背景音樂信號加入相同均值為2噪聲，誤差信號幅值較三角波實驗大。如圖，從頻域上看，濾波之后，在100~200一級800~900之間含有少量噪聲頻率。濾波效果不如輸入為周期信號。第五章 結(jié)論與展望 結(jié)論⑴通過對比實驗，分析收斂性和濾波效果，確定了濾波器正則化系數(shù)的選取，信噪比高正則化系數(shù)選取小正常數(shù)，反之取大正常數(shù)。⑵對不同類型濾波器輸入信號獨立實驗，周期型號和非周期信號的濾波器效果都較好，但周期信號的濾波效果更好。但從程序的運行時間上看，其復(fù)雜度高，計算量比較大。 對進一步研究的展望自適應(yīng)信號處理技術(shù)廣泛應(yīng)用于系統(tǒng)辨識、回波消除、自適應(yīng)譜線增強、自適應(yīng)信道均衡、語音線性預(yù)測、預(yù)測解卷、信號檢測、自適應(yīng)噪聲消除、自適應(yīng)天線陣等諸多信號處理領(lǐng)域中。自從Widrow提出著名的LMS算法以來，人們已經(jīng)提出了許多自適應(yīng)濾波算法，對自適應(yīng)濾波算法的研究己經(jīng)成為當(dāng)今自適應(yīng)信號處理中最為活躍的研究課題之一。對自適應(yīng)算法的優(yōu)化改進將大大的提高自適應(yīng)濾波器的性能。也是今后學(xué)者們致力的方向。參 考 文 獻[1] 于洋. 一種改進的LMS自適應(yīng)濾波器設(shè)計[D].[本科學(xué)位論文].北京：北京石油化工學(xué)院，2009[2] 丁玉美. 數(shù)字信號處理[M].西安：西安電子科技大學(xué)出版社，2002[3] 王姣. 自適應(yīng)濾波器的MATLAB仿真[D]. [本科學(xué)位論文].北京：北京石油化工學(xué)院，2002[4] 薛小梅. 一種基于LMS算法的自適應(yīng)濾波器設(shè)計[D]. [本科學(xué)位論文].北京：北京石油化工學(xué)院，2006[5] 吳燕. 一種自適應(yīng)濾波器的設(shè)計與仿真[J].科技廣場，2010，第7期：5961[6] 程玉柱. 自適應(yīng)濾波器的算法比較研究[J].大慶師范學(xué)報，2008，第2期：2326[7] 張立萍. RLS自適應(yīng)濾波器的Matlabe設(shè)計與仿真[J].赤峰學(xué)院學(xué)報，2011，第4期：2526[8] 蒙俊甫. LMS自適應(yīng)濾波器的仿真與實現(xiàn)[J].設(shè)計參考，2009，第9期：6773[9] 程玉柱. 自適應(yīng)濾波器的算法比較研究[J].大慶師范學(xué)報，2008，第2期：2326[10] 吳正茂. 自適應(yīng)濾波器及其應(yīng)用研究[J].南昌水專學(xué)報，2004，第2期：3640[11] 向大威、顧亞平. 自適應(yīng)濾波器的幾種全更新算法[J]中國科學(xué)院東海研究站[12] Simon Haykin. Adaptive Filter Thearo[M].北京：電子工業(yè)出版社，2006[13] 李勇、徐震. MATLAB輔助現(xiàn)代工程數(shù)字信號處理[M].西安電子科技大學(xué)出版社，[14] 王世一. 數(shù)字信號處理[M].北京理工大學(xué)出版社，2012[15] 張立萍. RLS自適應(yīng)濾波器的Matlabe設(shè)計與仿真[J]. 赤峰學(xué)院學(xué)報，2011，第4期：252646