【導(dǎo)讀】畢業(yè)論文（設(shè)計(jì)）。課題基于LMS的自適應(yīng)濾波器的設(shè)計(jì)。院部電氣工程學(xué)院。專業(yè)班級(jí)09電子信息工程班

　　

【正文】 應(yīng)用 [M]．北京：北京國(guó)防工業(yè)出版社， 20xx [11]鄒國(guó)良．自適應(yīng)濾波理論及應(yīng)用 [M]．河北：河北大學(xué)出版社， 1997 [12]高西全 丁玉美．?dāng)?shù)字信號(hào)處理 [M]. 陜西：西安電子科 技大學(xué)出版社， 20xx [13]邱天爽等 通信中的自適應(yīng)信號(hào)處理 [M] 北京：電子工業(yè)出版社， 20xx [14]Simon Haykin． Adaptive filter theory [M]． Electronic industry press， 20xx [15]zhang jianhua. Digital signal processing[M]. Electronic industry press， 20xx 熊偉：基于 LMS 的自適應(yīng)濾波器的設(shè)計(jì) 16 致 謝 本科學(xué)習(xí)生涯即將畫上句號(hào)，四年短暫的時(shí)光還歷歷在目，無(wú)論在 生活上還是在學(xué)習(xí)上，我都有許多收獲與提高。再次我衷心感謝那些幫助我的人。 首先，我要感謝我的導(dǎo)師馮老師對(duì)的我的悉心指導(dǎo)，馮老師是一個(gè)學(xué)識(shí)淵博、治學(xué)嚴(yán)謹(jǐn)?shù)膬?yōu)秀老師，對(duì)我的影響很大，我受益匪淺。她不僅在學(xué)術(shù)上耐心指導(dǎo)和教育我，而且在生活中給我不少幫助。借此機(jī)會(huì)，謹(jǐn)向馮老師表達(dá)我真摯的敬意和衷心的感謝。 其次感謝銅陵學(xué)院領(lǐng)導(dǎo)和老師給我的幫助及創(chuàng)造的良好的學(xué)習(xí)氛圍，讓我在大學(xué)四年的時(shí)間中，獲益匪淺。 同時(shí)，感謝與我在同一個(gè)教室的同班同學(xué)，感激他們?cè)谏詈蛯W(xué)習(xí)中給我提供的幫助，讓我覺(jué)得大家是一個(gè)整體，有家一樣的感覺(jué)。 感謝各位教授和老師在百忙之中評(píng)閱我的論文。 最后，向所有給我?guī)椭椭С值娜吮硎局孕牡母兄x，謝謝大家 ！ 銅陵學(xué)院畢業(yè)論文 (設(shè)計(jì) ) 17 附 錄 基本的 LMS 算法的實(shí)現(xiàn)程序 %基本的 LMS算法 %u（ K） =常數(shù) ,u是一個(gè)控制因子 %設(shè)輸入隨機(jī)信號(hào) v（ n），樣本間相互獨(dú)立，用高斯白噪聲代替 %且功率為 ，即 E((r(n))^2)=，則 v(n)=*randn(1,1024)。 %設(shè) W(0)=[0,0]， u（ K） = clear。 close all。 tic。 %計(jì)時(shí) M=1024 。 % M是濾波器階數(shù) ,可以改變，直到計(jì)算到達(dá)到 e(n)預(yù)期小的值為止，或 w(n)達(dá)到穩(wěn)定為止 u=。 %控制因子 vn=*randn(1,M)。%輸入隨機(jī)信號(hào) t=1:1024。 dn=sin(4*pi*t/100)。%標(biāo)準(zhǔn)信號(hào) w1=zeros(1,M)。 w2=zeros(1,M)。 for n=2:M xn=sin(4*pi*n/100)+vn。%輸入信號(hào) yn(n)=w1(n)*xn(n)+w2(n)*xn(n1)。 e(n)=xn(n)yn(n)。 w1(n+1)=w1(n)+2*u*e(n)*xn(n)。 w2(n+1)=w2(n)+2*u*e(n)*xn(n1)。 end figure(1) % LMS算法的權(quán)值隨時(shí)間 n變換的軌跡： W1(n)～ n； W2(n)～ n subplot(2,1,1)。 plot(w1,39。r39。)。 hold on plot(w2,39。k39。)。title(39。W1(n)～ n(紅線 ),W2(n)～ n（黑線） 39。)。 hold off % 誤差 e(n)隨時(shí)間變化的關(guān)系曲線： e(n)～ n subplot(2,1,2) 熊偉：基于 LMS 的自適應(yīng)濾波器的設(shè)計(jì) 18 plot(e)。title(39。e(n)～ n39。)。 figure(2) subplot(3,1,1)。 plot(dn)。title(39。dn標(biāo)準(zhǔn)信號(hào) 39。)。 axis([0 1200 2 2])。 subplot(3,1,2)。 plot(xn)。title(39。xn輸入信號(hào) 39。)。 subplot(3,1,3)。 plot(yn)。title(39。yn輸出信號(hào) 39。)。 axis([0 1200 2 2])。 toc 歸一化 LMS 算法實(shí)現(xiàn)程序 %歸一化的 LMS %u(k)=a/(b+x(k)x(k)39。) 0 a2,b=0 %設(shè)輸入隨機(jī)信號(hào) v（ n），樣本間相互獨(dú)立，用高斯白噪聲代替 %設(shè)功率為 ，即 E((r(n))^2)=，則 v(n)=*randn(1,1024)。 %設(shè) W(0)=[0,0]， a=1,b=0 clear。 close all。 tic。 %計(jì)時(shí) M=1024 。 % M 是濾波器階數(shù) ,可以改變，直到計(jì)算到達(dá)到 e(n)預(yù)期小的值為止，或 w(n)達(dá)到穩(wěn)定為止 u=。 %控制因子 vn=*randn(1,M)。%輸入隨機(jī)信號(hào) t=1:1024。 dn=sin(4*pi*t/100)。%標(biāo)準(zhǔn)信號(hào) w1=zeros(1,M)。 w2=zeros(1,M)。 a=1。b=0。 for n=2:M xn=sin(4*pi*n/100)+vn。%輸入信號(hào) yn(n)=w1(n)*xn(n)+w2(n)*xn(n1)。 e(n)=xn(n)yn(n)。 u(n)=1/(xn*xn39。)。 w1(n+1)=w1(n)+2*u(n)*e(n)*xn(n)。 銅陵學(xué)院畢業(yè)論文 (設(shè)計(jì) ) 19 w2(n+1)=w2(n)+2*u(n)*e(n)*xn(n1)。 end figure(1) % LMS算法的權(quán)值隨時(shí)間 n變換的軌跡： W1(n)～ n； W2(n)～ n subplot(2,1,1)。 plot(w1,39。r39。)。 hold on plot(w2,39。k39。)。title(39。W1(n)～ n(紅線 ),W2(n)～ n（黑線） 39。)。 hold off % 誤差 e(n)隨時(shí)間變化的關(guān)系曲線： e(n)～ n subplot(2,1,2) plot(e)。title(39。e(n)～ n39。)。 figure(2) subplot(3,1,1)。 plot(dn)。title(39。dn標(biāo)準(zhǔn)信號(hào) 39。)。 axis([0 1200 2 2])。 subplot(3,1,2)。 plot(xn)。title(39。xn輸入信號(hào) 39。)。 subplot(3,1,3)。 plot(yn)。title(39。yn輸出信號(hào) 39。)。 axis([0 1200 2 2])。 toc 功率歸一化 LMS 算法的實(shí)現(xiàn)程序 %功率歸一化 LMS算法 %u(k)=a/(gx(k)*gx(k)) gx(k)表示 x(k)的方差 %其中 gx(k)*gx(k)=l*gx(k1)*gx(k1)+e(k)*e(k)可以遞推算出， 0l=1,0a2/M, M是濾波器階數(shù) %設(shè)輸入隨機(jī)信號(hào) v（ n），樣本間相互獨(dú)立，用高斯白噪聲代替 %設(shè)功率為 ，即 E((r(n))^2)=，則 v(n)=*randn(1,1024)。 %設(shè) W(0)=[0,0]， a=1/M,l= clear。 close all。 tic。 %計(jì)時(shí) M=1024 。 % M是濾波器階數(shù) ,可以改變，直到計(jì)算到達(dá)到 e(n)預(yù)期小的值為止，或 w(n)達(dá)到熊偉：基于 LMS 的自適應(yīng)濾波器的設(shè)計(jì) 20 穩(wěn)定為止 u=。 %控制因子 a=1/M。 l=1。 vn=*randn(1,M)。%輸入隨機(jī)信號(hào) t=1:1024。 dn=sin(4*pi*t/100)。%標(biāo)準(zhǔn)信號(hào) w1=zeros(1,M)。 w2=zeros(1,M)。 gx2=zeros(1,M)。%gx2(k)=gx(k)*gx(k) a=1。b=0。 for n=2:M xn=sin(4*pi*n/100)+vn。%輸入信號(hào) yn(n)=w1(n)*xn(n)+w2(n)*xn(n1)。 e(n)=xn(n)yn(n)。 gx2(n)=l*gx2(n1)+e(n)*e(n)。 u(n)=a/(gx2(n))。 w1(n+1)=w1(n)+2*u(n)*e(n)*xn(n)。 w2(n+1)=w2(n)+2*u(n)*e(n)*xn(n1)。 end figure(1) % LMS算法的權(quán)值隨時(shí)間 n變換的 軌跡： W1(n)～ n； W2(n)～ n subplot(2,1,1)。 plot(w1,39。r39。)。 hold on plot(w2,39。k39。)。title(39。W1(n)～ n(紅線 ),W2(n)～ n（黑線） 39。)。 hold off % 誤差 e(n)隨時(shí)間變化的關(guān)系曲線： e(n)～ n subplot(2,1,2) plot(e)。title(39。e(n)～ n39。)。 figure(2) subplot(3,1,1)。 plot(dn)。title(39。dn標(biāo)準(zhǔn)信號(hào) 39。)。 axis([0 1200 2 2])。 subplot(3,1,2)。 plot(xn)。title(39。xn輸入信號(hào) 39。)。 銅陵學(xué)院畢業(yè)論文 (設(shè)計(jì) ) 21 subplot(3,1,3)。 plot(yn)。title(39。yn輸出信號(hào) 39。)。 axis([0 1200 2 2])。 Toc