基于dsp芯片的自適應(yīng)濾波器實現(xiàn)課程設(shè)計-在線瀏覽
2024-08-05 01:26本頁面
【正文】 (25)而由式（24）給出的正交條件則變?yōu)椋篍{[d(n)y(n)]x(n)}=0把式（25）代入上式后，有 E[d(n)x(n)]= E[x(n)xT(n)]h (26)如果令：r=E[d(n)x(n)],Фxx= E[x(n)xT(n)],那么最佳權(quán)向量 h*=Фxx 1r （27）3 詳細(xì)設(shè)計過程 LMS自適應(yīng)算法自適應(yīng)算法是根據(jù)某種判斷來設(shè)計的。LMS 算法是以最小均方誤差為判據(jù)的最典型的算法，也是應(yīng)用最廣泛的一種算法。誤差序列的均方值又叫“均方誤差”(Mean Sqluare Error，MSE)。由此產(chǎn)生的算法稱為LMS。FIR濾波器的最主要特點是沒有反饋回路，因此是無條件穩(wěn)定系統(tǒng)，其單位脈沖響應(yīng)h(n)是一個有限長序列。不斷地輸入樣本x(n)，經(jīng)延時(z1)做乘法累加，再輸出濾波結(jié)果Y(n)對式(1)進(jìn)行z變換．整理后可得FIR濾波器的傳遞函數(shù)為：。 其中，X(n )=[x(n)，x(n1)，…，x(nN+1)] 為自適應(yīng)濾波器的輸入矢量；W( n)=[w0(1 )，w1 (n )， …，wn1(n )]T 為系統(tǒng)的權(quán)系數(shù)矢量；為時間序列；N 為濾波器的階數(shù)；T為矩轉(zhuǎn)置符。 自適應(yīng)濾波器DSP設(shè)計實現(xiàn)實驗中采用的自適應(yīng)濾波器采用16階FIR濾波器，采用相同的信號作為參考d(n)和輸入信號x(n)，并采用上一時刻的誤差值來修正本時刻的濾波器系數(shù)，對濾波器輸出除128進(jìn)行幅度限制。 (2) 根據(jù)輸入的采樣值計算濾波器的輸出并求出誤差； 這一步是最重要的濾波計算，我們用FIRLMS(int *nx,float *nh,int nError,int nCoeffNumber) 這個函數(shù)來實現(xiàn)。 float fWork。 for ( i=0。i++ ) { fWork=nx[i]*nError*fU。 r+=(nx[ii]*nh[i])。 return r。nError是上一次的誤差值, fWork是當(dāng)前的濾波器權(quán)系數(shù),輸出值r。輸出值y（n）等于輸入值x(n1)*w(i)的積分。189。 (3) 根據(jù)LMS算法的迭代公式更新濾波器參數(shù)，有新的采樣輸入后轉(zhuǎn)到第二步循環(huán)執(zhí)行。把采樣點放在數(shù)字中。通過for循環(huán)語句進(jìn)行一個一個樣點值濾波計算。for ( i=COEFFNUMBER+1。i++ ) { nLastOutput=FIRLMS(xx+i,h,nLastOutputxx[i1],COEFFNUMBER)。 wc[i]=rr[i]xx[i]。for循環(huán)一個一個得到濾波過的信號值第INPUTNUMBER=1024個為止。通過本課程設(shè)計加深了對DSP原理的理解，初步掌握了DSP 芯片的開發(fā)應(yīng)用，為接下來的深入學(xué)習(xí)打下了堅實的基礎(chǔ)。 參考文獻(xiàn) [1] 吳湘淇,肖熙,[M].北京:電子工業(yè)出版社,2002[2] 萬建偉,王 [M].長沙:國防科技大學(xué)出版社,2008[3] John G. Proakis,Dimitris [M].方艷梅,:電子工業(yè)出版社,2007[4] Sanjit K [M].:電子工業(yè)出版社,2006 附件源程序清單：/******************************************************************************Main()主函數(shù)。 no Onus ///////////////////////////////////////////////////includedefine PI define COEFFNUMBER 16define INPUTNUMBER 1024int FIRLMS(int *nx,float *nh,int nError,int nCoeffNumber)。int xx[INP
外語相關(guān)推薦
鄂ICP備17016276號-1