【正文】 { s16 j,k。其使用方法如下：short oflag = dlms (DATA *x, DATA *h, DATA *r, DATA **d, DATA *des, DATA step, ushort nh, ushort nx)其C子函數(shù)定義的各個(gè)變量的意義如下表所示： oflag溢出錯(cuò)誤標(biāo)志，為1時(shí)在計(jì)算過程中發(fā)生了數(shù)據(jù)溢出；為0時(shí)無(wú)數(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。這個(gè)實(shí)驗(yàn)采用的LMS算法，即最小均方誤差算法設(shè)計(jì)的自適應(yīng)濾波器進(jìn)行未知系統(tǒng)識(shí)別，以將回聲信號(hào)濾除?；夭óa(chǎn)生程序如下：while (!MCBSP_RRDY(HANDSET_CODEC)) {}。利用DSP產(chǎn)生回波在打電話的過程中，常常能聽見自己說(shuō)話的聲音在話筒里重復(fù)，實(shí)際上就是自己的聲音泄露到了接收系統(tǒng)中。首先有一個(gè)橫向?yàn)V波器(圍繞它構(gòu)造LMS算法)，該部件的作用在于完成濾波過程。濾波器系數(shù)可通過MATLAB設(shè)計(jì)得到。要使用A/D和D/A轉(zhuǎn)換器，必須首先對(duì)A/D和D/A轉(zhuǎn)換器進(jìn)行初始化設(shè)置，即設(shè)置A/D轉(zhuǎn)換器的工作模式（15+1bit模式，16bit模式）、輸入增益（0dB,6dB,12dB）以及抽樣頻率（8000Hz,16000Hz）等。輸入信號(hào)經(jīng)A/D轉(zhuǎn)換后，利用DMA通道2進(jìn)行傳輸進(jìn)入DSP處理，而濾波后的信號(hào)輸出利用了DMA通道3進(jìn)行傳輸進(jìn)入D/A轉(zhuǎn)換器等輸出，因此需要對(duì)DMA通道2和通道3進(jìn)行初始化配置。本實(shí)驗(yàn)利用DMA通道2與MCBSP1通道結(jié)合來(lái)讀取A/D轉(zhuǎn)換器的數(shù)據(jù)，利用DMA通道3與MCBSP1通道結(jié)合來(lái)將處理后的數(shù)據(jù)發(fā)送至D/A轉(zhuǎn)換器。目的地址寄存器DMDST規(guī)定DMA要傳送數(shù)據(jù)目的地址和首地址。二、設(shè)計(jì)要求基本部分：（1） 對(duì)DMA進(jìn)行初始化；（2） 對(duì)A/D、D/A進(jìn)行初始化；（3） 編寫DMA中斷服務(wù)程序和濾波算法程序，實(shí)現(xiàn)語(yǔ)音信號(hào)的實(shí)時(shí)濾波；發(fā)揮部分：（1） 使用DSP產(chǎn)生帶回波的語(yǔ)音信號(hào)；（2） 利用自適應(yīng)濾波實(shí)現(xiàn)語(yǔ)音信號(hào)的回波對(duì)消。但是由于環(huán)境的原因，我們采集到語(yǔ)音信號(hào)經(jīng)常含有不同程度的噪聲。設(shè)計(jì)流程框圖如下：上電復(fù)位A/D和D/A初始化DMA通道初始化DMA方式接收A/D數(shù)據(jù)等待DMA接收中斷信號(hào)濾波信號(hào)經(jīng)D/A輸出主程序中斷服務(wù)程序模式選擇模式0普通濾波模式1自適應(yīng)濾波加有噪聲色語(yǔ)音信號(hào)由DSK板的MIC端口輸入，經(jīng)A/D轉(zhuǎn)換器轉(zhuǎn)換為數(shù)字信號(hào)后到達(dá)DSP的DMA通道2，DSP接收到信號(hào)后調(diào)用中斷服務(wù)程序進(jìn)行濾波，可選擇濾波模式進(jìn)行不同種類的濾波，濾波后的信號(hào)在CCS上顯示，或通過DMA通道3將信號(hào)輸出到D/A轉(zhuǎn)換器，轉(zhuǎn)換為模擬信號(hào)后從SPEAKER端口輸出。CTMOD位設(shè)置發(fā)送計(jì)數(shù)模式控制，CTMOD位為0時(shí)，DMA工作在多幀模式，CTMOD位為1時(shí)，DMA工作在自動(dòng)緩沖ABU模式。選擇MCBSP1通道的發(fā)送寄存器DXR11（43H）為DMA通道3傳送數(shù)據(jù)的目的地址，并選擇DMA源地址工作在不調(diào)整方式，選擇MCBSP1接收事件為DMA同步事件，以實(shí)現(xiàn)DMA和MCBSP的結(jié)合。D/A轉(zhuǎn)換器是將數(shù)字信號(hào)轉(zhuǎn)換為模擬信號(hào)的電路。FIR濾波為有限沖擊響應(yīng)濾波，其濾波結(jié)構(gòu)是一個(gè)分節(jié)的延時(shí)線，每節(jié)的輸出加權(quán)累加,得到濾波器的輸出。此外，它不需要計(jì)算有關(guān)的相關(guān)函數(shù)，也不需要矩陣求逆運(yùn)算。通過LMS算法計(jì)算這個(gè)向量所得的值表示一個(gè)估計(jì)，當(dāng)?shù)螖?shù)趨于無(wú)窮時(shí)，該估計(jì)的期望值可能接近維納最優(yōu)解,在濾波過程中，期望響應(yīng)d(k)與抽頭輸入向量z(k)一道參與處理。即將麥克輸入的語(yǔ)音信號(hào)疊加上述麥克輸入語(yǔ)音信號(hào)的延時(shí)衰減后的信號(hào)，用此混合信號(hào)模擬真實(shí)的回波信號(hào)。 //延遲+衰減利用自適應(yīng)濾波實(shí)現(xiàn)回波對(duì)消聲學(xué)回聲消除的功能原理框圖如下圖所示：其基本原理可概括為，自適應(yīng)地合成回聲，并從有回聲干擾地信號(hào)中減去該合成回聲。在k時(shí)刻，它們之差為剩余回聲信號(hào)：用它來(lái)控制LMS自適應(yīng)濾波器的系數(shù)ω(n)其中μ為自適應(yīng)步長(zhǎng)因子，一般取0μ1，可視為常數(shù)，k時(shí)刻濾波器的輸出為：以上三個(gè)方程是LMS最小均方誤差算法的核心方程，也是C語(yǔ)言編程的依據(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è)計(jì)實(shí)驗(yàn)中用到的源程序有以下幾個(gè)：（）/*****************************************************************************//* */ /* */ /* *//* *//*****************************************************************************/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 。但是該實(shí)驗(yàn)的DSK板資源有限，當(dāng)階數(shù)過高時(shí)，超過64，效果反而大打折扣。首先感覺這次實(shí)驗(yàn)的難度很大。但當(dāng)我開始運(yùn)行程序的時(shí)候才發(fā)現(xiàn)存在很多問題，比如在加載很多程序的時(shí)候，沒有考慮好它們和我們的程序的對(duì)應(yīng)性；；頭文件路徑不對(duì)，堆棧空間未設(shè)置，變量未定義或者重復(fù)定義等等的各種問題。十一、參考文獻(xiàn)【1】高海林 錢滿義 DSP技術(shù)及其應(yīng)用 北京交通大學(xué)出版社，2005【2】譚浩強(qiáng) C程序設(shè)計(jì)（第二版）清華大學(xué)出版社，2003【3】冬雷 DSP原理及開放技術(shù) 北京：清華大學(xué)出版社 北京交通大學(xué)出版社 【4】王樹恩 宋彥 一種快速的回聲消除算法及DSP實(shí)現(xiàn) 合肥：中國(guó)科學(xué)技術(shù)大學(xué)電子工程與信息科學(xué)系 【5】亓淑敏 梁文家 張美娟 基于DSP的自適應(yīng)噪聲消除系統(tǒng) 西安：長(zhǎ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。后來(lái)通過查閱書籍、與同學(xué)討論、向老師請(qǐng)教等方式，經(jīng)過多次的修改，終于可以拿出能夠正常運(yùn)行并得到正確結(jié)果的程序。在經(jīng)過老師三次的課堂講解之后，我只是基本上學(xué)會(huì)了CCS的操作，對(duì)DSP處理的程序根本不知其所以然。從blackman窗到 hamming窗和hann窗，阻帶衰減變小，通帶平坦度下降，主瓣長(zhǎng)度變窄，但過渡帶變窄。打開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 /******************************