【文章內(nèi)容簡介】 if(i==5000) i=0。 if(j==5000) j=0。 if(m==5000) m=0。 bf1[i]=out_h[0]。i++。//放入緩沖區(qū) out_delay=*bf1[i+1]+*bf1[j++]+bf1[m++]。 yk=dk+out_delay。 *(volatile u16*)DXR1_ADDR(HANDSET_CODEC)=yk。 } else if(mode==2) { fir(x,coff_h,out_h,amp。dp_h,LENGTH_H,LENGTH_IN)。 if(i==5000) i=0。 if(j==5000) j=0。 if(m==5000) m=0。 bf1[i]=out_h[0]。i++。 out_delay=*bf1[i+1]+*bf1[j++]+bf1[m++]。 yk=dk+out_delay。 update(y,yk)。 dlms(y,coff_w,out_w,amp。dp_w,out_h,STEP,LENGTH_W,LENGTH_IN)。 //ek=out_h[0]out_w[0]。 DSP 課程設(shè)計(jì) 10 *(volatile u16*)DXR1_ADDR(HANDSET_CODEC)=out_w[i]。 } } } void delay(void) { long int j。 for(j=0。 j100000。 j++) asm(_nop)。 } void update(DATA x[],DATA dk) { s16 j,k。 for(j=1。jLENGTH_IN。++j) { k=LENGTH_INj。 x[k]=x[k1]。 } x[0]=dk。 } void initarray(DATA x[]) { s16 i。 for(i=0。iLENGTH_IN。++i) { x[i]=0。 } } DSP 課程設(shè)計(jì) 11 的程序流程圖如下圖所示 : 首 先 ，在輸入信號之前，應(yīng)先進(jìn)行一些必要的準(zhǔn)備工作，即 初始化 FIR延遲輸出數(shù)組 dbuffer_h[]、 FIR輸出數(shù)組、自適應(yīng)濾波器系數(shù)矢量 coff_w[]和自適應(yīng)濾波器的輸出數(shù)組 out_w[]，再進(jìn)行 DSK板的初始化。然后采集輸入的語音信號，根據(jù)模式的不同進(jìn)行不同的處理以輸出不同的信號來進(jìn)行效果的比較。 mode=0時(shí)為全通，即不做任何處理直接輸出； mode=1時(shí)信號通過 FIR濾波器后產(chǎn)生回聲，為使現(xiàn)象更明顯可以附加一個(gè)延遲； mode=2時(shí)用自適應(yīng)濾波器逼近 FIR濾 波器以使之輸出近似的信號，用信號通過兩個(gè)系統(tǒng)產(chǎn)生的信號之差 (即剩余回聲 )來調(diào)整自適應(yīng)濾波器的系數(shù)，使剩余回聲逐漸收斂于 0。這時(shí)DSP 課程設(shè)計(jì) 12 只聽到自己說話的聲音，因此這個(gè)系統(tǒng)可以看做是近端說話模式，即將遠(yuǎn)端信號設(shè)為 0。其系統(tǒng)框圖如下 圖所示 ： 六、 程序設(shè)計(jì)及調(diào)試 由于實(shí)驗(yàn)條件限制（實(shí)驗(yàn)室的 DSP 板只有 1 個(gè)輸入端），不能像實(shí)際線路中那樣有近端回波、遠(yuǎn)端回波， 所以 該實(shí)驗(yàn)中 ，我們 簡化了這個(gè)回聲產(chǎn)生。 實(shí)驗(yàn)中，我們 將麥克輸入的語音信號疊加上述麥克輸入語音信號的延時(shí)衰減后的信號，用此混合信號模擬真實(shí)的回波信號。 回波產(chǎn)生程序如下： while (!MCBSP_RRDY(HANDSET_CODEC)) {}。 //codec 句柄如未準(zhǔn)備好則等待 dk = *(volatile u16*)DRR1_ADDR(HANDSET_CODEC)。 //從 ADC 讀數(shù) bf[i]=dk。 //放入緩沖區(qū) if(i==10000) i=0。 *(volatile u16*)DXR1_ADDR(HANDSET_CODEC) = dk+bf[i+1]/2。 //延遲 +衰減 調(diào)試中 可以 模擬產(chǎn)生較大延時(shí)的回波。 當(dāng) 對著麥克說話時(shí) ， 不僅能耳塞中不僅能聽到自己及時(shí)的聲音也能聽到自己的回聲，而且延時(shí)特別大，感覺特別明顯。 但是 我們編寫的 回波抵消算法（ LMS 算法）程序 的效果卻并不DSP 課程設(shè)計(jì) 13 盡如人意。 噪聲特別大， 甚至 根本沒法聽清自己的聲音 。 分析原因可能是采樣輸入數(shù)據(jù)和 lms 算法效率不高，產(chǎn)生很大延時(shí)，也產(chǎn)生了很大的噪聲， 從而嚴(yán)重影響實(shí)驗(yàn)的結(jié)果 。 算法的開銷與濾波器的長度有關(guān)，因此應(yīng)該適當(dāng)減小濾波器的長度的值。也有可是 FIR 濾波器的階數(shù)不夠。 既然， 因?yàn)樽约壕幍?LMS 算法不能很好的抵消回波 ， 所以我們對于我們所編寫的程序進(jìn)行了一定地改進(jìn)。這一次， 我們采用了C5000 的 dsplib 庫中的函數(shù) short dlms(PARAM)取代自定的 lms()函數(shù)，同時(shí)用 dsplib 庫里的函數(shù) short fir(PARAM)產(chǎn)生回聲通道的輸出信號 (也就是用預(yù)定階數(shù)的 FIR 濾波器模擬回聲通道 ),該輸出信號作為自適應(yīng)濾波器的理想輸出信號，對之進(jìn)行逼近進(jìn)而得到預(yù)定結(jié)果。 這兩個(gè)函數(shù)是用 C5000 匯編語言編寫的，具有很高的執(zhí)行效率，使用這個(gè)函數(shù)可以大大 提高運(yùn)算速度，以解決用 C 語言編寫的程序效率不能滿足要求的問題。其使用方法如下： short oflag = dlms (DATA *x, DATA *h, DATA *r, DATA **d, DATA *des, DATA step, ushort nh, ushort nx) 參數(shù)如下表所示： DSP 課程設(shè)計(jì) 14 注意： h 和 d 的存儲地址要求起點(diǎn)必須位于 Kbit 邊界處，即起始地址值的最低有效位必