【正文】 hpibuff0 HPRAM0 PAGE 1 hpibuff1 HPRAM1 PAGE 1 hpibuff2 HPRAM2 PAGE 1 dma_buff DMARAM PAGE 1 } 七、 結(jié)果分析 經(jīng)過上述的調(diào)試之后，下面我們對于得到的最終程序，進(jìn)行一下系統(tǒng)分析： （ 1） 在回波抵消過程中 ， 通過 dlms設(shè)置斷點觀察 自適應(yīng)濾波器的 系DSP 課程設(shè)計 20 數(shù)矢量 coff_w趨向于 Fir濾波器的系數(shù)矢量 coff_h的過程 斷點設(shè)置如下圖所示： coff_h觀察結(jié)果如下圖所示： coff_w觀察結(jié)果如下圖所示： (2)當(dāng)輸入語音信號變化時， 觀察輸入變量 dk的變化 如下圖所示： DSP 課程設(shè)計 21 (3)在回波產(chǎn)生過程中， Mode=1 時，輸入 語音 信號 變化時 coff_w 不變 ，過程圖如下所示： (4)在 回波產(chǎn)生過程中 ， 當(dāng) Mode=2時， out_w和 out_h的對比圖如下圖示： 綜上所述， 當(dāng) mode=1 時，如果不加附加延遲產(chǎn)生的聲音信號有輕微不易察覺的延遲和比 較明顯的雜音，這是聲音經(jīng)過 10 階 FIR 濾DSP 課程設(shè)計 22 波器濾波即回聲通道的結(jié)果，當(dāng)加入附加延時后可以模擬比較逼真的回聲效果，當(dāng)然這并不是真正的回聲，也不是 LMS 算法要逼近的部分，F(xiàn)IR 的輸出 out_h[]才是真正的回聲。 //定義 FIR濾波器的系數(shù)即回聲通道的權(quán)系數(shù) pragma DATA_SECTION(dbuffer_h,.dbufferh) //將數(shù)組 dbuffer_h指定到內(nèi)存段 .dbufferh DATA dbuffer_h[LENGTH_IN]。 pragma DATA_SECTION(coff_w,.coffw) //將數(shù)組 coff_w指定到內(nèi)存段 .coffw DATA coff_w[LENGTH_W]。 //FIR濾波器即回聲消除通道輸出 DATA out_w[LENGTH_IN]。 (1)添加頭文件 /**/ define LENGTH_W 10 //LMS濾波器長度 , define LENGTH_H 10 //自適應(yīng)濾波器長度 define STEP 100 //自適應(yīng)步長 define LENGTH_IN 5 //輸入緩沖數(shù)組 x[]長度 DSP 課程設(shè)計 17 DATA x[LENGTH_IN]。 在編完主程序后，因引入了新的數(shù)據(jù)和變量 ，必須在工程文件中添加一個 頭文件，將這些變量放入指定 的 內(nèi)存中。而使用宏 progma 來指定自定義段時應(yīng)該在頭 文件中進(jìn)行，否則編譯會提示警告： , line 46: warning: pragma can only be applied to a file level symbol, not 39。 用 C 語言頭文件輸出數(shù)據(jù)，數(shù)據(jù)類型選 16 位有符號整型（以對應(yīng) DSP 程序中的 DATA 類型），生成的系數(shù)如下 (非實驗采用的數(shù)據(jù) )： const int16_T B[10] = {1218, 1611, 2317, 2900, 3231, 3231, 2900, 2317, 1611,1218}。 由于聲音為低頻信號，用 MATLAB的 FDATOOL工具產(chǎn)生 1個低通 FIR濾波器。 //從低位壓向高位 } x[0]=dk。++j) { k=LENGTH_INj。 for(j=1。采用數(shù)據(jù)壓入方式，將最新的聲音采樣值存入數(shù)組x[0],然后將存儲的數(shù)據(jù)從低位壓入高位。 DATA 為 中定義的 short 整型。 oflag = short fir (DATA *x, DATA *h, DATA *r, DATA **dbuffer, ushort nh, ushort nx) oflag 溢出錯誤標(biāo)志， =1 時在計算過程中發(fā)生了數(shù)據(jù)溢出； =0 時無數(shù)據(jù)溢出。其使用方法如下： short oflag = dlms (DATA *x, DATA *h, DATA *r, DATA **d, DATA *des, DATA step, ushort nh, ushort nx) 參數(shù)如下表所示： DSP 課程設(shè)計 14 注意： h 和 d 的存儲地址要求起點必須位于 Kbit 邊界處，即起始地址值的最低有效位必須是 0，必須在 CMD 文件中開辟對齊 1K邊界的空間，并且只能在頭文件中進(jìn)行段的配置。這一次， 我們采用了C5000 的 dsplib 庫中的函數(shù) short dlms(PARAM)取代自定的 lms()函數(shù)，同時用 dsplib 庫里的函數(shù) short fir(PARAM)產(chǎn)生回聲通道的輸出信號 (也就是用預(yù)定階數(shù)的 FIR 濾波器模擬回聲通道 ),該輸出信號作為自適應(yīng)濾波器的理想輸出信號，對之進(jìn)行逼近進(jìn)而得到預(yù)定結(jié)果。也有可是 FIR 濾波器的階數(shù)不夠。 分析原因可能是采樣輸入數(shù)據(jù)和 lms 算法效率不高，產(chǎn)生很大延時，也產(chǎn)生了很大的噪聲， 從而嚴(yán)重影響實驗的結(jié)果 。 但是 我們編寫的 回波抵消算法（ LMS 算法）程序 的效果卻并不DSP 課程設(shè)計 13 盡如人意。 //延遲 +衰減 調(diào)試中 可以 模擬產(chǎn)生較大延時的回波。 //放入緩沖區(qū) if(i==10000) i=0。 //codec 句柄如未準(zhǔn)備好則等待 dk = *(volatile u16*)DRR1_ADDR(HANDSET_CODEC)。 實驗中，我們 將麥克輸入的語音信號疊加上述麥克輸入語音信號的延時衰減后的信號，用此混合信號模擬真實的回波信號。這時DSP 課程設(shè)計 12 只聽到自己說話的聲音，因此這個系統(tǒng)可以看做是近端說話模式，即將遠(yuǎn)端信號設(shè)為 0。然后采集輸入的語音信號，根據(jù)模式的不同進(jìn)行不同的處理以輸出不同的信號來進(jìn)行效果的比較。++i) { x[i]=0。 for(i=0。 } x[0]=dk。++j) { k=LENGTH_INj。 for(j=1。 j++) asm(_nop)。 for(j=0。 DSP 課程設(shè)計 10 *(volatile u16*)DXR1_ADDR(HANDSET_CODEC)=out_w[i]。dp_w,out_h,STEP,LENGTH_W,LENGTH_IN)。 update(y,yk)。 out_delay=*bf1[i+1]+*bf1[j++]+bf1[m++]。 bf1[i]=out_h[0]。 if(j==5000) j=0。dp_h,LENGTH_H,LENGTH_IN)。 *(volatile u16*)DXR1_ADDR(HANDSET_CODEC)=yk。//放入緩沖區(qū) out_delay=*bf1[i+1]+*bf1[j++]+bf1[m++]。 bf1[i]=out_h[0]。 if(j==5000) j=0。dp_h,LENGTH_H,LENGTH_IN)。 if(mode==0) { *(volatile u16*)DXR1_ADDR(HANDSET_CODEC)=dk。 // *(volatile u16*)DXR1_ADDR(HANDSET_CODEC)=dk。 while(1) //讀數(shù)循環(huán) { while (!MCBSP_RRDY(HANDSET_CODEC)) {}。 j=3000。 initarray(x)。 /* 12dB gain on analog output from DAC */ codec_sample_rate(hHandset,SR_8000)。 /* ADC in 15bit mode */ codec_ain_gain(hHandset, CODEC_AIN_0dB)。 /* Acquire handle to codec */ /* Set codec para