【正文】 temp1=sqrt(p*p+q*q)。 rfft(y,SS,0)。 XPC=0。(delaybuff[0])。/*高通濾波器*//*這里COEFF的數(shù)據(jù)類型曾經(jīng)是 int*//*帶通濾波器*//*帶阻濾波器*//*全通濾波器*//*存放延遲數(shù)據(jù)*/pragma DATA_SECTION(delaybuff,delayb)。unsigned int dmpre,dmsrcp,dmdstp,dmidx0,dmidx1,dmfri0,dmfri1,dmgsa,dmgda,dmgcr,dmgfr。點擊“Desgin filter” 完成濾波器設(shè)計，濾波器的波形如下圖所示。 for(j=1。x輸入信號數(shù)組h自適應(yīng)濾波器系數(shù)矢量r輸出數(shù)組d延遲緩沖區(qū)，用來存放上一時刻的輸出數(shù)據(jù)des理想輸出信號數(shù)組step為了適應(yīng)DSP定點數(shù)的運算一般取STEP為327nx向量x中的個數(shù)nh系數(shù)的個數(shù)注意：h和d的存儲地址要求起點必須位于Kbit邊界處，即起始地址值的最低有效位必須是0，必須在CMD文件中開辟對齊1K邊界的空間，并且只能在頭文件中進行段的配置。r[3*NX/4], amp。r[NX/4], amp。 iNH。 i++) db[i] = 0。short eflag2= PASS。dp,des,STEP, NH, NX)。i++) h[i] =0。 } return(ek)。 for(i=0。該自適應(yīng)濾波器是FIR橫向濾波器，可以根據(jù)輸出自動修改濾波器的權(quán)系數(shù)，從而逼近未知系統(tǒng)回聲通道。 //codec句柄如未準備好則等待dk = *(volatile u16*)DRR1_ADDR(HANDSET_CODEC)。在一些特定的系統(tǒng)中出現(xiàn)這種現(xiàn)象將會對信號的正常收發(fā)產(chǎn)生不利影響。其次，有一個對橫向濾波器抽頭權(quán)值進行自適應(yīng)控制過程的算法。 LMS算法LMS(leastmeansquare)算法，即最小均方算法，是由Widrow和Hof在1960年創(chuàng)建的，直到現(xiàn)在它還是應(yīng)用最廣泛的自適應(yīng)濾波算法。FIR濾波器算法數(shù)字濾波器原理一般具有如下差分方程:y(n)=∑akx(nk)+ ∑bky(nk) （1）式中x(n)為輸人序列，y(n)為輸出序列，ak和bk為濾波器系數(shù)，,N是濾波器階數(shù)。詳細的DMA初始化程序請見后面的程序清單。選擇MCBSP1通道的接收寄存器DRR11（41H）為DMA傳送數(shù)據(jù)的首地址，并選擇DMA源地址工作在不調(diào)整方式，選擇MCBSP1接收事件為DMA同步事件，以實現(xiàn)DMA和MCBSP的結(jié)合。單元計數(shù)寄存器DMCTR規(guī)定DMA傳送數(shù)據(jù)的個數(shù)為DMCTR寄存器的值加1。（3） 比較采用不同窗和階數(shù)時濾波器的濾波效果；（4） 測試所設(shè)計濾波器的幅頻特性和相頻特性，并與MATLAB下的設(shè)計結(jié)果進行比較。語音通信的過程主要分為三個部分，即發(fā)送部分——將音頻信號轉(zhuǎn)化為電信號，經(jīng)發(fā)送端設(shè)備變換為適合傳輸?shù)男问桨l(fā)送出去；傳輸部分——通過傳輸信道將電信號進行傳輸；接收部分——經(jīng)接收端設(shè)備恢復(fù)出原來的語音信號，經(jīng)耳機或者喇叭轉(zhuǎn)換為接收者可以聽到的聲音信號。濾波后的數(shù)據(jù)利用DMA方式送到D/A轉(zhuǎn)換器轉(zhuǎn)換為模擬信號。DINM位和IMOD位設(shè)置DMA中斷產(chǎn)生方式。當(dāng)一組數(shù)據(jù)處理完成后，將數(shù)據(jù)存放在數(shù)據(jù)存儲區(qū)OUTBUFFER中，選擇OUTBUFFER為DMA通道3的首地址，并選擇源工作地址工作在不調(diào)整方式。A/D轉(zhuǎn)換器能夠?qū)⒔邮盏恼Z音信號（模擬信號）轉(zhuǎn)化為數(shù)字信號，供CPU進行處理。下面為FIR濾波器的算法實現(xiàn)方法。LMS算法的一個顯著特點就是它的簡單性。相應(yīng)的h0(k), h1 (k), ...， hm1（k)為抽頭權(quán)向量h(k)的元素。這次實驗我們的回聲信號定義為語音信號的延遲再進行適當(dāng)?shù)乃p。 *(volatile u16*)DXR1_ADDR(HANDSET_CODEC) = dk+bf[i+1]/2。遠端信號通過自適應(yīng)濾波器產(chǎn)生回聲預(yù)測信號y(n)。 } ek=dkyk。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++) r[i] =0。 // for h if (eflagr != PASS) { exit(1)。void main(void){ /* 1. Test for singlebuffer */ /* clear */ for (i=0。 /* test */ eflag1 = test (r, rtest, NX, MAXERROR)。db[0]。x[2*NX/4], h, amp。 if (eflag1 != PASS) { exit(1)。x輸入信號數(shù)組hFIR濾波器的系數(shù)矢量，也是要逼近的對象r輸出數(shù)組d延遲緩沖區(qū)，用來存放上一時刻的輸出數(shù)據(jù)nx向量x中的個數(shù)nh系數(shù)的個數(shù)注意：h、d的地址要求同dlms()函數(shù)中的h、d參數(shù)。 x[k]=x[k1]。const real64_T B[16] = {1737, 22, 494, 1762, 4117, 2011, 5080, 11684, 11684, 5080, 2011, 4117, 1762, 494, 22, 1737}。/*存放輸出數(shù)據(jù)*/pragma DATA_SECTION(out_buffer,outt_buffer)。int flag=0。double p,q,temp1。 TIMER_HALT(0)。 mm=1。 }/******************************************************************************/ /*****************************************************************************/ if(brd_init_bios()) return。/*設(shè)置codec參數(shù)*/ codec_dac_mode(hHandset,CODEC_DAC_15BIT)。 dmmcr=((AUTOINIT_ENABLE15)|(DINM_ENABLE14)|(IMOD_HALFBLOCK13)|(CTMOD_DEC12)| (INDEXMODE_NOMOD8)|(SPACE_DATA6)|(INDEXMODE_INC2)|(SPACE_DATA))。 dmgsa=src_addr。 dmidx1=0。 for(。 }}（）includeincludeincludeincludeextern void delay(s16 period)。extern unsigned int dst_addr。extern int *delayptr1。 dmctr=0xFF。 int *p_out1。 init_dma3()。} }（）MEMORY{ PAGE 0: VEC1: origin = 0080h, length = 0078h /* Internal Program RAM */ PRAM: origin = 2000h, length = 2000h /* Internal Program RAM */ VEC2: origin = 0xff80, length = 0x78 PAGE 1: SCRATCH: origin = 0060h, length = 0020h /* Scratch Pad Data RAM */ DMARAM: origin = 0100h, length = 0400h /* DMA buffer */ STACK: origin = 0600h, length = 0400h /* Stack Memory Space */ EXRAM: origin = 1100h, length = 0e00h /* External Data RAM */}/*****************************************************************************//* DSP Memory Allocation *//*****************************************************************************/SECTIONS{ .cinit PRAM PAGE 0 .text PRAM PAGE 0 .vectors VEC2 PAGE 0 .vecs VEC1 PAGE 0 .stack STACK PAGE 1 .trap SCRATCH PAGE 1 .bss EXRAM PAGE 1 .cio EXRAM PAGE 1 .fifo EXRAM PAGE 1 .sintab EXRAM PAGE 1 audio_buffer DMARAM PAGE 1 outt_buffer DMARAM PAGE 1 .const EXRAM PAGE 1 coefficients EXRAM ,align (512)