【正文】 }。iN。采樣時(shí)間為 define f1=500 。 X[j]=X[i]。 while(k=j){j=k。 /* FFT 的點(diǎn)數(shù) */ unsigned int N2=N1。 =。 X[i].re=X[i].re*scale+。 =(X[id].re*[id].im*)*scale。jLE1。 LE1=LE1。 scale = 。 /* FFT 的點(diǎn)數(shù) */ unsigned int L。 /* 復(fù)變量臨時(shí)存儲(chǔ)器 */ plex U。 i++) { = X[i].re*X[i].re。 } /* 啟動(dòng) FFT */ bit_rev(X,M)。 i++) { /* 構(gòu)造輸入信號(hào)樣本 */ X[i].re =xin[i]。 } /* */ generator(xin,N)。 L++) { LE=1L。 /* 說(shuō)明功率譜信號(hào)數(shù)組 ,為實(shí)數(shù) */float re1[N],im1[N]。 /* 說(shuō)明輸入信號(hào)數(shù)組 ,為復(fù)數(shù) */plex W[M]。（ 1）主程序： include include /*包含浮點(diǎn)復(fù)數(shù)結(jié)構(gòu)體定義頭文件*/extern void bit_rev(plex *,unsigned int)。l 按照自然排序的時(shí)域信號(hào)數(shù)據(jù)是 x(0)、 x(1)、 x(2)、x(3)、 x(4)、 x(5)、 x(6)、 x(7)，其序號(hào)寫成二進(jìn)制數(shù)分別為 000b、 001b、 010b、 011b、 100b、101b、 110b、 111b，將這些二進(jìn)制數(shù)前后倒轉(zhuǎn)，即得到進(jìn)行 FFT前數(shù)據(jù)所對(duì)應(yīng)的實(shí)際二進(jìn)制數(shù)地址：000b、 100b、 010b、 110b、 001b、 101b、 011b、 111b，對(duì)應(yīng)的十進(jìn)制數(shù)是： 0、 7?；?2的 DIT又被稱為庫(kù)利一圖基算法。通過(guò) FFT算法， DFT的計(jì)算量大大減少，運(yùn)算時(shí)間縮短 1~2個(gè)數(shù)量級(jí)。按 Q15格式存放結(jié)果sample_loop popboth XAR7 。wk(n2)=wk(n) || mpym *AR7+,T3,AC0 。啟動(dòng) IIR濾波器環(huán) mov *AR0+ 14,AC0 。初始化系數(shù)基地址 amov 0,AR3 。設(shè)置循環(huán)緩沖區(qū)大小 sfts T0,1 mov mmap(T0),BK03 。 樣本數(shù) 1 mov T0,BRC0 。 保存 T3 pshboth XAR7 。 : : 。 w1j a2k。參數(shù) 5:AR3 – 延遲線指針。參數(shù) 1:T0 輸入信號(hào)緩沖區(qū)的樣本數(shù)。void iir(int *,unsigned int,int *,int *,unsigned int,int *)。 IIR濾波器的匯編語(yǔ)言實(shí)現(xiàn) 例 712，采用指針編寫的 IIR濾波器匯編語(yǔ)言程序。 w[1][k]=w[0][k]。kN_IIR。l 二階 IIR濾波器傳遞函數(shù)為：圖 75 二階 IIR濾波器的直接型 Ⅰ 的實(shí)現(xiàn)圖 76 的信號(hào)流圖圖 77 二階 IIR濾波器的直接型 Ⅱ 的實(shí)現(xiàn)由于直接型 Ⅱ 對(duì)于給定的傳遞函數(shù)具有最小可能的延遲數(shù)、加法器數(shù)和乘法器數(shù)，所以被稱為標(biāo)準(zhǔn)型。 恢復(fù) ST1, ST2和 ST3popm ST2_55 popm ST1_55 ret.end IIR濾波器l IIR濾波器的優(yōu)點(diǎn)：結(jié)構(gòu)簡(jiǎn)單，運(yùn)算量小，可以用較少的階數(shù)獲得很高的選擇性。設(shè)置外部循環(huán)次數(shù)為 L1 mpym *AR0+,*AR2+,AC0 。AR0和 AR2為循環(huán)緩沖區(qū)指針mov T1,AR0 。置位 SMULmov mmap(AR0),BSA01 。 參數(shù) 2: T0 FIR 濾波器的階數(shù) L； 參數(shù) 3: T1 輸入信號(hào)當(dāng)前值在循環(huán)緩沖區(qū)的序數(shù)。 該程序用于實(shí)現(xiàn) FIR濾波器 ,可被 C語(yǔ)言程序調(diào)用。 out[index] = fir(in,coeff,L,index)。 signal += a2*cos((float)w2*n)。iL。 float signal。/* Lowpass FIR filter coefficients */int coeff[L]={ 26,13,14,36,31,8,58,71,15,83,139,76,90,231,194,50,331,383,78,405, 654,347,403,1024,863,228,1577,1972,453,2910,6836,9470,9470,6836, 2910,453,1972,1577,228,863, 1024,403,347,654,405,78,383,331,50,194,231,90,76,139,83,15,71, 58,8,31,36,14,13,26}。l 在循環(huán)緩沖器中，采取數(shù)據(jù)保持固定、反時(shí)針?lè)较蛞苿?dòng)地址的方式，代替保持緩沖器地址固定且正方向移動(dòng)數(shù)據(jù)。}u 直接型 FIR濾波器的實(shí)現(xiàn)涉及到兩個(gè)基本操作，一個(gè)是輸入信號(hào)向量與濾波器系數(shù)向量的內(nèi)積計(jì)算，另一個(gè)是輸入信號(hào)向量的更新處理。for(i=0。i0。/********************************************************* 該程序用于實(shí)現(xiàn) FIR濾波器* L——濾波器的階數(shù)* b[i]——濾波器的系數(shù) ,i=0,1,…,L1* x[i]——輸入信號(hào)向量 ,i=0,1,…,L1。l 本節(jié)主要討論 FIR濾波器的 DSP實(shí)現(xiàn)方法，有關(guān)IIR濾波器的實(shí)現(xiàn)將在下一節(jié)中介紹。 AC0=X0*Y1 MACM *AR0,uns(*AR1),AC0 。限制條件： W1被指定為偶地址。 SXMD = 1 (允許符號(hào)擴(kuò)展 ) W1。 Y1 Y0 Q31操作數(shù) X0。u 為了解決冗余符號(hào)位的問(wèn)題，需要在程序中設(shè)定狀態(tài)寄存器 ST1中的 FRCT（小數(shù)方式）為 1，這樣當(dāng)乘法器將結(jié)果傳送至累加器時(shí)就會(huì)自動(dòng)左移 1位。u 可以用于遞歸運(yùn)算。l 為了克服這些缺陷，在實(shí)際應(yīng)用中更多采用的是小數(shù)運(yùn)算。 否則 ,對(duì)商取反 .NEG AC1,AC1MOV AC1,dbl(*AR2)skip: B L1 小數(shù)乘法l 在定點(diǎn) DSP的某些應(yīng)用中，整數(shù)運(yùn)算很難滿足要求。 執(zhí)行 subc 16次SUBC *AR3,AC1,AC1MOV AC1,*AR2 。 最后一次執(zhí)行 subc|| MOV 11,AR4 。 QuotH = abs(NumH) 。 AC1 = abs(Den)MOV AC1,*AR3 。*******************************************************.sect .text.align 4.global start.sym start,start,36,2,0start:MOV 0100h,AR0 MOV 0102h,AR1 MOV 0104h,AR2 MOV 0106h,AR3 MOV 0108h,AR4L1:BSET SXMD 。第一次除法運(yùn)算之前，把被除數(shù)的高半部分存入 AC1。 注：。 AR1 除數(shù) (Den)。******************************************************。 對(duì)商求反XCCPART label,AC00 。AR2 |除數(shù) |MOV HI(AC1),*AR3 。 暫時(shí)保存 |除數(shù) |MOV *AR0,AC1 。 置位 SXMD (打開(kāi)符號(hào)擴(kuò)展 )MPYM *AR0,*AR1,AC0 。打開(kāi)符號(hào)擴(kuò)展，被除數(shù)、除數(shù)可為負(fù)數(shù)。 AR3 余數(shù)。 指針?lè)峙洹?置位 SXMD (打開(kāi)符號(hào)擴(kuò)展 ) B L1例 76，帶符號(hào) 16位除 16位整數(shù)除法。 把被除數(shù)低位裝入 AC0_LRPT (161) 。 AR1指向除數(shù)SUBC *AR1,AC0,AC0 。************************************************.sect .text.align 4.global start.sym start,start,36,2,0start: MOV 0100h,AR0 MOV 0102h,AR1 MOV 0104h,AR2 MOV 0106h,AR3L1: BCLR SXMD 。第一次除法之前，把被除數(shù)高位存入 AC0。 AR1除數(shù)。***********************************************。 保存商MOV HI(AC0),*AR3 。 清零 SXMD (關(guān)閉符號(hào)擴(kuò)展 )MOV *AR0,AC0 。 無(wú)符號(hào)除法，被除數(shù)、除數(shù)均為 16位。 AR2商。***********************************************。下面舉例說(shuō)明如何在 C55x DSP中實(shí)現(xiàn)除法運(yùn)算。保存 W3 B L1 除法運(yùn)算l C55x沒(méi)有提供硬件除法器，也沒(méi)有提供專門的除法指令，要實(shí)現(xiàn)除法運(yùn)算需借助于條件減法指令 SUBC和重復(fù)指令 RPT。AC0=X0*Y016+X1*Y0+X0*Y1MOV AC0,*AR2+ 。AR1指向 Y0MPYM uns(*AR0),uns(*AR1),AC0 。限制條件：延遲鏈和輸入序列必須指定為長(zhǎng)字類型。 SXMD = 1 (允許符號(hào)擴(kuò)展 ) W1。 Y1 Y0 32位操作數(shù) X0。*************************************************************