【正文】 { s16 j,k。其使用方法如下：short oflag = dlms (DATA *x, DATA *h, DATA *r, DATA **d, DATA *des, DATA step, ushort nh, ushort nx)其C子函數(shù)定義的各個變量的意義如下表所示： oflag溢出錯誤標(biāo)志，為1時在計算過程中發(fā)生了數(shù)據(jù)溢出；為0時無數(shù)據(jù)溢出。x[3*NX/4], h, amp。x[NX/4], h, amp。 // clear output buffer (optional) for (i=0。 iNH。short eflag1= PASS。 // clear delay buffer (a must) /* pute */ dlms(x,h,r,amp。iNH。i++) { w[i]=w[i]+uek*x[i]。 yk=0。這個實驗采用的LMS算法，即最小均方誤差算法設(shè)計的自適應(yīng)濾波器進(jìn)行未知系統(tǒng)識別，以將回聲信號濾除。回波產(chǎn)生程序如下：while (!MCBSP_RRDY(HANDSET_CODEC)) {}。利用DSP產(chǎn)生回波在打電話的過程中，常常能聽見自己說話的聲音在話筒里重復(fù)，實際上就是自己的聲音泄露到了接收系統(tǒng)中。首先有一個橫向濾波器(圍繞它構(gòu)造LMS算法)，該部件的作用在于完成濾波過程。濾波器系數(shù)可通過MATLAB設(shè)計得到。要使用A/D和D/A轉(zhuǎn)換器，必須首先對A/D和D/A轉(zhuǎn)換器進(jìn)行初始化設(shè)置，即設(shè)置A/D轉(zhuǎn)換器的工作模式（15+1bit模式，16bit模式）、輸入增益（0dB,6dB,12dB）以及抽樣頻率（8000Hz,16000Hz）等。輸入信號經(jīng)A/D轉(zhuǎn)換后，利用DMA通道2進(jìn)行傳輸進(jìn)入DSP處理，而濾波后的信號輸出利用了DMA通道3進(jìn)行傳輸進(jìn)入D/A轉(zhuǎn)換器等輸出，因此需要對DMA通道2和通道3進(jìn)行初始化配置。本實驗利用DMA通道2與MCBSP1通道結(jié)合來讀取A/D轉(zhuǎn)換器的數(shù)據(jù)，利用DMA通道3與MCBSP1通道結(jié)合來將處理后的數(shù)據(jù)發(fā)送至D/A轉(zhuǎn)換器。目的地址寄存器DMDST規(guī)定DMA要傳送數(shù)據(jù)目的地址和首地址。二、設(shè)計要求基本部分：（1） 對DMA進(jìn)行初始化；（2） 對A/D、D/A進(jìn)行初始化；（3） 編寫DMA中斷服務(wù)程序和濾波算法程序，實現(xiàn)語音信號的實時濾波；發(fā)揮部分：（1） 使用DSP產(chǎn)生帶回波的語音信號；（2） 利用自適應(yīng)濾波實現(xiàn)語音信號的回波對消。但是由于環(huán)境的原因，我們采集到語音信號經(jīng)常含有不同程度的噪聲。設(shè)計流程框圖如下：上電復(fù)位A/D和D/A初始化DMA通道初始化DMA方式接收A/D數(shù)據(jù)等待DMA接收中斷信號濾波信號經(jīng)D/A輸出主程序中斷服務(wù)程序模式選擇模式0普通濾波模式1自適應(yīng)濾波加有噪聲色語音信號由DSK板的MIC端口輸入，經(jīng)A/D轉(zhuǎn)換器轉(zhuǎn)換為數(shù)字信號后到達(dá)DSP的DMA通道2，DSP接收到信號后調(diào)用中斷服務(wù)程序進(jìn)行濾波，可選擇濾波模式進(jìn)行不同種類的濾波，濾波后的信號在CCS上顯示，或通過DMA通道3將信號輸出到D/A轉(zhuǎn)換器，轉(zhuǎn)換為模擬信號后從SPEAKER端口輸出。CTMOD位設(shè)置發(fā)送計數(shù)模式控制，CTMOD位為0時，DMA工作在多幀模式，CTMOD位為1時，DMA工作在自動緩沖ABU模式。選擇MCBSP1通道的發(fā)送寄存器DXR11（43H）為DMA通道3傳送數(shù)據(jù)的目的地址，并選擇DMA源地址工作在不調(diào)整方式，選擇MCBSP1接收事件為DMA同步事件，以實現(xiàn)DMA和MCBSP的結(jié)合。D/A轉(zhuǎn)換器是將數(shù)字信號轉(zhuǎn)換為模擬信號的電路。FIR濾波為有限沖擊響應(yīng)濾波，其濾波結(jié)構(gòu)是一個分節(jié)的延時線，每節(jié)的輸出加權(quán)累加,得到濾波器的輸出。此外，它不需要計算有關(guān)的相關(guān)函數(shù)，也不需要矩陣求逆運(yùn)算。通過LMS算法計算這個向量所得的值表示一個估計，當(dāng)?shù)螖?shù)趨于無窮時，該估計的期望值可能接近維納最優(yōu)解,在濾波過程中，期望響應(yīng)d(k)與抽頭輸入向量z(k)一道參與處理。即將麥克輸入的語音信號疊加上述麥克輸入語音信號的延時衰減后的信號，用此混合信號模擬真實的回波信號。 //延遲+衰減利用自適應(yīng)濾波實現(xiàn)回波對消聲學(xué)回聲消除的功能原理框圖如下圖所示：其基本原理可概括為，自適應(yīng)地合成回聲，并從有回聲干擾地信號中減去該合成回聲。在k時刻，它們之差為剩余回聲信號：用它來控制LMS自適應(yīng)濾波器的系數(shù)ω(n)其中μ為自適應(yīng)步長因子，一般取0μ1，可視為常數(shù)，k時刻濾波器的輸出為：以上三個方程是LMS最小均方誤差算法的核心方程，也是C語言編程的依據(jù)。 uek=2*mu*ek。short eflagr= PASS。 // clear output buffer (optional) for (i=0。 } if (eflagh != PASS) { exit(1)。 iNX。 /* 2. Tests for dualbuffer */ /* clear */ for (i=0。 /* pute */ if (NX=4) { fir(x, h, r, amp。r[2*NX/4], amp。 } if (eflag2 != PASS) { exit(2)。DATA 。 //從低位壓向高位 } x[0]=dk。七、程序設(shè)計實驗中用到的源程序有以下幾個：（）/*****************************************************************************//* */ /* */ /* *//* *//*****************************************************************************/include include include include include include include void delay(s16 period)。int out_buffer[0x200]。int temp。u16 f[SS/2+1]。 brd_set_wait_states(7,7,9)。 for(iii=2。 while (t) { brd_led_toggle(BRD_LED0)。 codec_adc_mode(hHandset,CODEC_ADC_15BIT)。 dmmcr=0xe045。 dmgda=dst_addr。 dmfri0=0。)。extern unsigned int channel。extern int inp_buffer[0x200]。extern int L=0。 src_addr=(unsigned int)amp。if(choose==1) { p_inp=inp_buffer+frame*0x100。 frame^=1。 RESET vector NOP NOPNMI: BD NMI 。但是該實驗的DSK板資源有限，當(dāng)階數(shù)過高時，超過64，效果反而大打折扣。首先感覺這次實驗的難度很大。但當(dāng)我開始運(yùn)行程序的時候才發(fā)現(xiàn)存在很多問題，比如在加載很多程序的時候，沒有考慮好它們和我們的程序的對應(yīng)性；；頭文件路徑不對，堆?？臻g未設(shè)置，變量未定義或者重復(fù)定義等等的各種問題。十一、參考文獻(xiàn)【1】高海林 錢滿義 DSP技術(shù)及其應(yīng)用 北京交通大學(xué)出版社，2005【2】譚浩強(qiáng) C程序設(shè)計（第二版）清華大學(xué)出版社，2003【3】冬雷 DSP原理及開放技術(shù) 北京：清華大學(xué)出版社 北京交通大學(xué)出版社 【4】王樹恩 宋彥 一種快速的回聲消除算法及DSP實現(xiàn) 合肥：中國科學(xué)技術(shù)大學(xué)電子工程與信息科學(xué)系 【5】亓淑敏 梁文家 張美娟 基于DSP的自適應(yīng)噪聲消除系統(tǒng) 西安：長安大學(xué)信息工程學(xué)院 【6】Thomas 著，—，2004【7】TMS320C54x DSP CPU and Peripherals. Texas Instrument Inc，2001.【8】TMS320C54x DSP Applications and Guide. Texas Instrument . 41 / 41。后來通過查閱書籍、與同學(xué)討論、向老師請教等方式，經(jīng)過多次的修改，終于可以拿出能夠正常運(yùn)行并得到正確結(jié)果的程序。在經(jīng)過老師三次的課堂講解之后，我只是基本上學(xué)會了CCS的操作，對DSP處理的程序根本不知其所以然。從blackman窗到 hamming窗和hann窗，阻帶衰減變小，通帶平坦度下降，主瓣長度變窄，但過渡帶變窄。打開CCS界面，新建任務(wù)fir，、。 p_out=out_buffer+frame*0x100。 p_out=out_buffer+frame*0x100。 dst_addr=DXR1_ADDR(HANDSET_CODEC)。extern int out1_buffer[0x200]。extern int coeffs[16]。extern unsigned int dmmcr。 for(i=0。 dma_global_init(dmpre,dmsrcp,dmdstp,dmidx0,dmidx1,dmfri0,dmfri1,dmgsa,dmgda,dmgcr,dmgfr)。 dmgfr=1。 src_addr=DRR1_ADDR(HANDSET_CODEC)。 codec_aout_gain(hHandset,CODEC_AOUT_MINUS_6dB)。 brd_led_toggle(BRD_LED1)。iii=iii+2) { p=y[iii]。 IMR=0。/*****************************************************************************//* MAIN *//*****************************************************************************/void main(){ s16 t=2。int choose=1。/*存放濾波器系數(shù)*/pragma DATA_SECTION(coeffs,coefficients)。define SS 16 /******************************