【正文】 ** DATE: 02/00 ** ************************************************************************** .title 54xDSKPLUS Vector Table Initialization .ref _c_int00, _DMAC2ISR .sect .vecsRESET: BD _c_int00 。 fir(p_inp,coeffs,p_out,amp。 init_dma3()。 fir(p_inp,coeffs,p_out1,amp。 int *p_out1。 dma_init(DMA_CH3,dmsefc,dmmcr,dmctr,SPACE_DATA,src_addr,SPACE_DATA,dst_addr)。 dmctr=0xFF。void init_dma3(void) { while(DMPRECamp。extern int *delayptr1。extern int delaybuff[16]。extern unsigned int dst_addr。extern unsigned int dmctr。 }}（）includeincludeincludeincludeextern void delay(s16 period)。iperiod。 for(。 DMA_ENABLE(DMA_CH2)。 dmidx1=0。 dmpre=((HIGH_PRIORITY10)|(INTSEL_016))。 dmgsa=src_addr。 dst_addr=(unsigned int)amp。 dmmcr=((AUTOINIT_ENABLE15)|(DINM_ENABLE14)|(IMOD_HALFBLOCK13)|(CTMOD_DEC12)| (INDEXMODE_NOMOD8)|(SPACE_DATA6)|(INDEXMODE_INC2)|(SPACE_DATA))。 codec_sample_rate(hHandset,SR_16000)。/*設(shè)置codec參數(shù)*/ codec_dac_mode(hHandset,CODEC_DAC_15BIT)。 delay(1000)。 }/******************************************************************************/ /*****************************************************************************/ if(brd_init_bios()) return。 q=y[iii+1]。 mm=1。 /************************計(jì)算濾波器的FFT******************************************/ cbrev(coeffs,y,SS/2)。 TIMER_HALT(0)。 BSCR=0x8806。double p,q,temp1。int *delayptr1= amp。int flag=0。/*低通濾波器*/DATA coeffs[16]={ 0, 113, 419, 878, 1031, 0, 3731, 19656, 19656, 3731, 0, 1031, 878, 419, 113, 0}。/*存放輸出數(shù)據(jù)*/pragma DATA_SECTION(out_buffer,outt_buffer)。unsigned int dmsefc,dmmcr,dmctr,src_addr,dst_addr。const real64_T B[16] = {1737, 22, 494, 1762, 4117, 2011, 5080, 11684, 11684, 5080, 2011, 4117, 1762, 494, 22, 1737}。在FDA中選擇帶通FIR濾波器，階數(shù)選擇15，抽樣頻率選擇16kHz， kHz， kHz。 x[k]=x[k1]。其實(shí)現(xiàn)程序update(DATA x[],DATA dk)如下：void update(DATA x[],DATA dk) //dk為當(dāng)前采樣值{ s16 j,k。x輸入信號(hào)數(shù)組hFIR濾波器的系數(shù)矢量，也是要逼近的對(duì)象r輸出數(shù)組d延遲緩沖區(qū)，用來存放上一時(shí)刻的輸出數(shù)據(jù)nx向量x中的個(gè)數(shù)nh系數(shù)的個(gè)數(shù)注意：h、d的地址要求同dlms()函數(shù)中的h、d參數(shù)。其使用方法如下：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í)無數(shù)據(jù)溢出。 if (eflag1 != PASS) { exit(1)。x[3*NX/4], h, amp。x[2*NX/4], h, amp。x[NX/4], h, amp。db[0]。 // clear output buffer (optional) for (i=0。 /* test */ eflag1 = test (r, rtest, NX, MAXERROR)。 iNH。void main(void){ /* 1. Test for singlebuffer */ /* clear */ for (i=0。short eflag1= PASS。 // for h if (eflagr != PASS) { exit(1)。 // clear delay buffer (a must) /* pute */ dlms(x,h,r,amp。i++) r[i] =0。iNH。Dsplib庫中的dlsm函數(shù)和fir函數(shù)的子程序：，路徑：C:\ti\c5400\dsplib\EXAMPLES\DLMS//*****************************************************************************// Filename: // Version: // Description: test for dlms routine//*****************************************************************************include include include include short i。i++) { w[i]=w[i]+uek*x[i]。 } ek=dkyk。 yk=0。遠(yuǎn)端信號(hào)通過自適應(yīng)濾波器產(chǎn)生回聲預(yù)測(cè)信號(hào)y(n)。這個(gè)實(shí)驗(yàn)采用的LMS算法，即最小均方誤差算法設(shè)計(jì)的自適應(yīng)濾波器進(jìn)行未知系統(tǒng)識(shí)別，以將回聲信號(hào)濾除。 *(volatile u16*)DXR1_ADDR(HANDSET_CODEC) = dk+bf[i+1]/2。回波產(chǎn)生程序如下：while (!MCBSP_RRDY(HANDSET_CODEC)) {}。這次實(shí)驗(yàn)我們的回聲信號(hào)定義為語音信號(hào)的延遲再進(jìn)行適當(dāng)?shù)乃p。利用DSP產(chǎn)生回波在打電話的過程中，常常能聽見自己說話的聲音在話筒里重復(fù)，實(shí)際上就是自己的聲音泄露到了接收系統(tǒng)中。相應(yīng)的h0(k), h1 (k), ...， hm1（k)為抽頭權(quán)向量h(k)的元素。首先有一個(gè)橫向?yàn)V波器(圍繞它構(gòu)造LMS算法)，該部件的作用在于完成濾波過程。LMS算法的一個(gè)顯著特點(diǎn)就是它的簡(jiǎn)單性。濾波器系數(shù)可通過MATLAB設(shè)計(jì)得到。下面為FIR濾波器的算法實(shí)現(xiàn)方法。要使用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）等。A/D轉(zhuǎn)換器能夠?qū)⒔邮盏恼Z音信號(hào)（模擬信號(hào)）轉(zhuǎn)化為數(shù)字信號(hào)，供CPU進(jìn)行處理。輸入信號(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)行初始化配置。當(dāng)一組數(shù)據(jù)處理完成后，將數(shù)據(jù)存放在數(shù)據(jù)存儲(chǔ)區(qū)OUTBUFFER中，選擇OUTBUFFER為DMA通道3的首地址，并選擇源工作地址工作在不調(diào)整方式。本實(shí)驗(yàn)利用DMA通道2與MCBSP1通道結(jié)合來讀取A/D轉(zhuǎn)換器的數(shù)據(jù)，利用DMA通道3與MCBSP1通道結(jié)合來將處理后的數(shù)據(jù)發(fā)送至D/A轉(zhuǎn)換器。DINM位和IMOD位設(shè)置DMA中斷產(chǎn)生方式。目的地址寄存器DMDST規(guī)定DMA要傳送數(shù)據(jù)目的地址和首地址。濾波后的數(shù)據(jù)利用DMA方式送到D/A轉(zhuǎn)換器轉(zhuǎn)換為模擬信號(hào)。二、設(shè)計(jì)要求基本部分：（1） 對(duì)DMA進(jìn)行初始化；（2） 對(duì)A/D、D/A進(jìn)行初始化；（3） 編寫DMA中斷服務(wù)程序和濾波算法程序，實(shí)現(xiàn)語音信號(hào)的實(shí)時(shí)濾波；發(fā)揮部分：（1） 使用DSP產(chǎn)生帶回波的語音信號(hào)；（2） 利用自適應(yīng)濾波實(shí)現(xiàn)語音信號(hào)的回波對(duì)消。語音通信的過程主要分為三個(gè)部分，即發(fā)送部分——將音頻信號(hào)轉(zhuǎn)化為電信號(hào)，經(jīng)發(fā)送端設(shè)備變換為適合傳輸?shù)男问桨l(fā)送出去；傳輸部分——通過傳輸信道將電信號(hào)進(jìn)行傳輸；接收部分——經(jīng)接收端設(shè)備恢復(fù)出原來的語音信號(hào)，經(jīng)耳機(jī)或者喇叭轉(zhuǎn)換為接收者可以聽到的聲音信號(hào)。但是由于環(huán)境的原因，我們采集到語音信號(hào)經(jīng)常含有不同程度的噪聲。（3） 比較采用不同窗和階數(shù)時(shí)濾波器的濾波效果；（4） 測(cè)試所設(shè)計(jì)濾波器的幅頻特性和相頻特性，并與MATLAB下的設(shè)計(jì)結(jié)果進(jìn)行比較。設(shè)計(jì)流程框圖如下：上電復(fù)位A/D和D/A初始化DMA通道初始化DMA方式接收A/D數(shù)據(jù)