【正文】 fir(p_inp,coeffs,p_out,amp。對(duì)于相位響應(yīng)，當(dāng)階數(shù)上升時(shí)，對(duì)于FIR濾波器，一直是比較理想的線性響應(yīng)。這次實(shí)驗(yàn)將理論上學(xué)習(xí)的DSP知識(shí)真正應(yīng)用到了實(shí)踐中，雖然掌握了理論的相關(guān)知識(shí)，但在與實(shí)踐相聯(lián)系的過程中還是遇到了不少問題。當(dāng)然自己寫程序有點(diǎn)超乎我的能力之外，所以程序思路大體上都是參考的資料。由頻域波形可看出，高頻的噪音在先經(jīng)過帶通濾波器后都先被濾掉了，輸出信號(hào)只有語音頻率帶有幅值； 不同低通濾波器的對(duì)比(1)FIR16階hann窗MATLAB設(shè)計(jì)結(jié)果：h[k]參數(shù)0, 41, 298, 761, 967, 0, 3696, 19616, 19616, 3696, 0, 967, 761, 98, 41, 0CCS中幅頻特性：(2) FIR32階hann窗MATLAB設(shè)計(jì)結(jié)果：h[k]參數(shù)0, 7, 18, 0, 83, 223, 341, 306, 0, 590, 1299, 1777, 1540, 0, 3993, 19789, 19789, 3993, 0, 1540, 1777, 1299, 590, 0, 306, 341, 223, 83, 0, 18, 7, 0CCS中幅頻特性(3)FIR64階hann窗MATLAB設(shè)計(jì)結(jié)果：h[k]參數(shù)0, 1, 3, 5, 0, 14, 34, 47, 39, 0, 65 ,132, 161, 121, 0, 172, 328, 385, 278, 0，376, 709, 827, 599, 0, 848, 1671, 2095, 1699, 0, 4064, 19828, 19828, 4064, 0, 1699, 2095, 1671, 848, 0, 599, 827, 709, 376, 0, 278, 385, 328, 172, 0, 121, 161, 132, 65, 0, 39, 47, 34, 14, 0, 5, 3, 1, 0CCS中幅頻特性：(4)FIR 16階blackman窗MATLAB設(shè)計(jì)結(jié)果：h[k]參數(shù)0, 16, 139, 443, 690, 0, 3471, 19484, 19484, 3471, 0, 690, 443, 139, 16, 0CCS中幅頻特性（5）FIR32階blackman窗MATLAB設(shè)計(jì)結(jié)果：h[k]參數(shù)0, 3, 7, 0, 38, 114, 194, 193, 0, 448, 1066, 1556, 1421, 0, 3935, 19756, 19756, 3935, 0, 1421, 1556, 1066, 448, 0, 193, 194, 114, 38, 0，7, 3, 0CCS中幅頻特性（6）FIR16階hamming窗MATLAB設(shè)計(jì)結(jié)果：h[k]參數(shù)0, 113, 419, 878, 1031, 0, 3731, 19656, 19656, 3731, 0, 1031, 878, 419, 113, 0CCS中幅頻特性（7）FIR32階hamming窗MATLAB設(shè)計(jì)結(jié)果：h[k]參數(shù)51, 61, 53, 0, 119, 280, 397, 339, 0, 617, 1338, 1809, 1555, 0, 3996, 19766, 19766, 3996, 0, 1555, 1809, 1338, 617, 0, 339, 397, 280, 119, 0, 53, 61, 51CCS中幅頻特性（8）FIR 64階hamming窗帶通濾波器MATLAB設(shè)計(jì)結(jié)果：h[k]參數(shù)22, 46, 1, 35, 26, 85, 31, 38, 154, 69, 135, 5, 326, 196, 261, 58, 307, 798, 127, 0, 170, 1421, 734, 189, 1161, 1218, 2895, 69, 2772, 2419, 10979, 9792, 9792, 10979, 2419, 2772, 69, 2895, 1218, 1161, 189, 734, 1421, 170, 0, 127, 798, 307, 58, 261, 196, 326, 5, 135, 69, 154, 38, 31, 85, 26, 35, 1, 46, 22CCS中幅頻特性（9）濾波器的相頻特性（FIR 16位blackman窗）MATLAB設(shè)計(jì)結(jié)果：設(shè)計(jì)相頻特性CCS中相頻特性結(jié)論：通過對(duì)各種窗函數(shù)的研究可以發(fā)現(xiàn)：在窗函數(shù)類型相同時(shí)，階數(shù)越高，通帶與阻帶越平坦，通帶與阻帶的過渡帶越窄。 init_dma3()。 dmctr=0xFF。extern unsigned int dst_addr。 for(。 dmgsa=src_addr。/*設(shè)置codec參數(shù)*/ codec_dac_mode(hHandset,CODEC_DAC_15BIT)。 mm=1。double p,q,temp1。/*存放輸出數(shù)據(jù)*/pragma DATA_SECTION(out_buffer,outt_buffer)。 x[k]=x[k1]。 if (eflag1 != PASS) { exit(1)。db[0]。void main(void){ /* 1. Test for singlebuffer */ /* clear */ for (i=0。i++) r[i] =0。 } ek=dkyk。 *(volatile u16*)DXR1_ADDR(HANDSET_CODEC) = dk+bf[i+1]/2。相應(yīng)的h0(k), h1 (k), ...， hm1（k)為抽頭權(quán)向量h(k)的元素。下面為FIR濾波器的算法實(shí)現(xiàn)方法。當(dāng)一組數(shù)據(jù)處理完成后，將數(shù)據(jù)存放在數(shù)據(jù)存儲(chǔ)區(qū)OUTBUFFER中，選擇OUTBUFFER為DMA通道3的首地址，并選擇源工作地址工作在不調(diào)整方式。濾波后的數(shù)據(jù)利用DMA方式送到D/A轉(zhuǎn)換器轉(zhuǎn)換為模擬信號(hào)。（3） 比較采用不同窗和階數(shù)時(shí)濾波器的濾波效果；（4） 測(cè)試所設(shè)計(jì)濾波器的幅頻特性和相頻特性，并與MATLAB下的設(shè)計(jì)結(jié)果進(jìn)行比較。選擇MCBSP1通道的接收寄存器DRR11（41H）為DMA傳送數(shù)據(jù)的首地址，并選擇DMA源地址工作在不調(diào)整方式，選擇MCBSP1接收事件為DMA同步事件，以實(shí)現(xiàn)DMA和MCBSP的結(jié)合。FIR濾波器算法數(shù)字濾波器原理一般具有如下差分方程:y(n)=∑akx(nk)+ ∑bky(nk) （1）式中x(n)為輸人序列，y(n)為輸出序列，ak和bk為濾波器系數(shù)，,N是濾波器階數(shù)。其次，有一個(gè)對(duì)橫向?yàn)V波器抽頭權(quán)值進(jìn)行自適應(yīng)控制過程的算法。 //codec句柄如未準(zhǔn)備好則等待dk = *(volatile u16*)DRR1_ADDR(HANDSET_CODEC)。 for(i=0。i++) h[i] =0。short eflag2= PASS。 iNH。r[3*NX/4], amp。 for(j=1。unsigned int dmpre,dmsrcp,dmdstp,dmidx0,dmidx1,dmfri0,dmfri1,dmgsa,dmgda,dmgcr,dmgfr。(delaybuff[0])。 rfft(y,SS,0)。 brd_led_toggle(BRD_LED2)。inp_buffer。 temp=*(volatile u16*)DRR1_ADDR(HANDSET_CODEC)。extern unsigned int src_page。0x0008){}。delayptr1,16,256)。將DSK板與計(jì)算機(jī)連接好，接通電源，在CCS中選擇File中的Load Program。于是我決定從老師上課用的例題程序開始，一點(diǎn)點(diǎn)的研讀程序，分析它的功能。我學(xué)的是自動(dòng)化專業(yè)，也許以后不會(huì)有很多機(jī)會(huì)接觸這樣的信號(hào)處理的任務(wù)，但是這次實(shí)驗(yàn)將給我以后完成其他實(shí)驗(yàn)時(shí)提供思路。其他特性變化不明顯。 frame^=1。}interrupt void DMAC2ISR(void){ int *p_inp,*p_out。extern int currbuff。j++)。 dmidx0=0。 /*初始化DMA2通道*/ dmsefc=((DSYNC_REVT112))。 mm++。 brd_set_cpu_freq(100)。int frame=0。/* * Filter Design and Analysis Tool Generated Filter Coefficients C Source * Generated by MATLAB Signal Processing Toolbox *//* General type conversion for MATLAB generated Ccode */include /* * Expected path to * C:\MATLAB6p5\extern\include\ */const int BL = 16。同理，：oflag = short fir (DATA *x, DATA *h, DATA *r, DATA **dbuffer, ushort nh, ushort nx)oflag溢出錯(cuò)誤標(biāo)志，為1時(shí)在計(jì)算過程中發(fā)生了數(shù)據(jù)溢出；為0時(shí)無數(shù)據(jù)溢出。 fir(amp。dbptr, NH, NX)。 // for r eflagh = test(h, htest, NH, MAXERROR)。同時(shí)替換原有的LMS算法的C語言實(shí)現(xiàn)程序，更換為CCS的diplib庫中的函數(shù)short dlms(PA