【正文】 這些問題有時會導(dǎo)致編譯時報錯，有時不報錯但會影響程序運行的結(jié)果。在階數(shù)相同，窗函數(shù)不同時，濾波器特性會有不同的效果。 }if (choose==0) { p_inp=inp_buffer+frame*0x100。out_buffer+((unsigned int)frame*0x100)。extern int out_buffer[0x200]。}void delay(s16 period){ int i,j。 dmgcr=0xFF。 codec_ain_gain(hHandset,CODEC_AIN_6dB)。iiiSS。int iii,mm。int out1_buffer[0x200]。 // 當(dāng)前采樣值存入最低位}六、用MATLAB設(shè)計FIR濾波器用MATLAB設(shè)計FIR流程如下：打開MATLAB，在命令窗中輸入fdatool,打開濾波器設(shè)計與分析工具。 } return。dbptr, NH, NX/4)。 i++) r[i] = 0。 iNH。 for(i=0。一般回聲消除算法通過自適應(yīng)濾波來完成，其基本原理如下圖所示：其中，遠(yuǎn)端信號x（k)通過回聲信道h產(chǎn)生回聲y(k)，近端信號d(k)是由回聲y(k)和近端聲音信號(可包含噪聲信號)得到。在這種情況下，給定一個輸入，橫向濾波器產(chǎn)生一個輸入，橫向濾波器產(chǎn)生一個輸出y(k)作為期望響應(yīng)d(k)的估計。由前面的分析知，F(xiàn)IR濾波器數(shù)學(xué)上可表示為:y(n)= ∑h(k)x(nk) （5）式（5）中 x[n]為最近(t=nT)的輸人信號，x[nk]是延時了k個取樣周期的輸人信號，h[k]是第k個延時的加權(quán)值(即濾波器系數(shù))，可由MATLAB設(shè)計出來，y[n]是時刻t=nT時濾波器的輸出信號，N是濾波器的階數(shù)也稱濾波器的抽頭數(shù)，為實整數(shù)，且須滿足N*t1/fs，其中fs為采樣頻率。在CCS集成開發(fā)環(huán)境中。四、算法原理直接存儲器訪問DMA直接存儲器訪問（Direct Memory Access）是C54x DSP非常重要的片上外設(shè)，DMA控制器可以在不影響CPU的情況下完成數(shù)據(jù)的傳輸，因此數(shù)據(jù)傳輸速度快，在要求信號實時采集和處理的系統(tǒng)中常采用DMA方式進行信號采集和傳輸。當(dāng)語音信號受到背景噪聲干擾而使語音通信質(zhì)量變得不可接受時，要對語音信號中的噪聲進行濾除，也就是本次實驗要完成的語音噪聲濾波。此外，DMA的6個通道還受通道優(yōu)先級和使能控制寄存器DMPREC控制。D/A轉(zhuǎn)換器能夠?qū)?shù)字信號轉(zhuǎn)換為語音信號送SPEAKER端口輸出。這兩個過程一起工作組成一個反饋環(huán)，如下圖所示。聽到的回波是混合信號，此時不僅能聽到自己的聲音，也能聽到幾個衰減后疊加的回聲，該回聲為近端輸入的延遲和衰減。 double uek,yk。void main(void){ /* clear */ for (i=0。}，路徑：C:\ti\c5400\dsplib\EXAMPLES\FIR//*****************************************************************************// Filename: // Version: // Description: test for fir routine//*****************************************************************************include include include include short i。 i++) r[i] = 0。 fir(amp。采用數(shù)據(jù)壓入方式，將最新的聲音采樣值存入數(shù)組x[0],然后將存儲的數(shù)據(jù)從低位壓入高位。define SS 16 /*****************************************************************************//* 全局變量 *//*****************************************************************************/HANDLE hHandset。int choose=1。 IMR=0。 brd_led_toggle(BRD_LED1)。 src_addr=DRR1_ADDR(HANDSET_CODEC)。 dma_global_init(dmpre,dmsrcp,dmdstp,dmidx0,dmidx1,dmfri0,dmfri1,dmgsa,dmgda,dmgcr,dmgfr)。extern unsigned int dmmcr。extern int out1_buffer[0x200]。 p_out=out_buffer+frame*0x100。打開CCS界面，新建任務(wù)fir，、。在經(jīng)過老師三次的課堂講解之后，我只是基本上學(xué)會了CCS的操作，對DSP處理的程序根本不知其所以然。十一、參考文獻【1】高海林 錢滿義 DSP技術(shù)及其應(yīng)用 北京交通大學(xué)出版社，2005【2】譚浩強 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。首先感覺這次實驗的難度很大。 RESET vector NOP NOPNMI: BD NMI 。if(choose==1) { p_inp=inp_buffer+frame*0x100。extern int L=0。extern unsigned int channel。 dmfri0=0。 dmmcr=0xe045。 while (t) { brd_led_toggle(BRD_LED0)。 brd_set_wait_states(7,7,9)。int temp。七、程序設(shè)計實驗中用到的源程序有以下幾個：（）/*****************************************************************************//* */ /* */ /* *//* *//*****************************************************************************/include include include include include include include void delay(s16 period)。DATA 。r[2*NX/4], amp。 /* 2. Tests for dualbuffer */ /* clear */ for (i=0。 } if (eflagh != PASS) { exit(1)。short eflagr= PASS。在k時刻，它們之差為剩余回聲信號：用它來控制LMS自適應(yīng)濾波器的系數(shù)ω(n)其中μ為自適應(yīng)步長因子，一般取0μ1，可視為常數(shù)，k時刻濾波器的輸出為：以上三個方程是LMS最小均方誤差算法的核心方程，也是C語言編程的依據(jù)。即將麥克輸入的語音信號疊加上述麥克輸入語音信號的延時衰減后的信號，用此混合信號模擬真實的回波信號。此外，它不需要計算有關(guān)的相關(guān)函數(shù)，也不需要矩陣求逆運算。D/A轉(zhuǎn)換器是將數(shù)字信號轉(zhuǎn)換為模擬信號的電路。CTMOD位設(shè)置發(fā)送計數(shù)模式控制，CTMOD位為0時，DMA工作在多幀模式，CTMOD位為1時，DMA工作在自動緩沖ABU模式。但是由于環(huán)境的原因，我們采集到語音信號經(jīng)常含有不同程度的噪聲。目的地址寄存器DMDST規(guī)定DMA要傳送數(shù)據(jù)目的地址和首地址。輸入信號經(jīng)A/D轉(zhuǎn)換后，利用DMA通道2進行傳輸進入DSP處理，而濾波后的信號輸出利用了DMA通道3進行傳輸進入D/A轉(zhuǎn)換器等輸出，因此需要對DMA通道2和通道3進行初始化配置。濾波器系數(shù)可通過MATLAB設(shè)計得到。利用DSP產(chǎn)生回波在打電話的過程中，常常能聽見自己說話的聲音在話筒里重復(fù)，實際上就是自己的聲音泄露到了接收系統(tǒng)中。這個實驗采用的LMS算法，即最小均方誤差算法設(shè)計的自適應(yīng)濾波器進行未知系統(tǒng)識別，以將回聲信號濾除。i++) { w[i]=w[i]+uek*x[i]。 // clear delay buffer (a must) /* pute */ dlms(x,h,r,amp。 iNH。x[NX/4], h, amp。其使用方法如下：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ù)溢出。在FDA中選擇帶通FIR濾波器，階數(shù)選擇15，抽樣頻率選擇16kHz， kHz， kHz。/*低通濾波器*/DATA coeffs[16]={ 0, 113, 419, 878, 1031, 0, 3731, 19656, 19656, 3731, 0, 1031, 878, 419, 113, 0}。 BSCR=0x8806。 q=y[iii+1]。 codec_sample_rate(hHandset,SR_16000)。 dmpre=((HIGH_PRIORITY10)|(INTSEL_016))。iperiod。extern int delaybuff[16]。 dma_init(DMA_CH3,dmsefc,dmmcr,dmctr,SPACE_DATA,src_addr,SPACE_DATA,dst_addr)。