DSPFFT試驗報告

1、DSP技術課程設計(FFT)作者姓名：鐘欣專業(yè)名稱：電子科學與技術指導教師：張雙目錄一、實驗結(jié)果3正弦信號輸入的結(jié)果3余弦信號輸入的結(jié)果4心得5參考文獻6附件1正弦信號程序7正弦信號文件7文件11附件2余弦信號程序12余弦信號文件12實驗結(jié)果正弦信號輸入的結(jié)果余弦信號輸入的結(jié)果心得通過本次的DSP設計課程，我學會了在CCS軟件平臺 上運用C語言編程來實現(xiàn)傅里葉運算。在課程中老師講解了 的基本原理，以及如何通過運用C語言編程來實現(xiàn)。實驗過程中，通過對傅里葉變換的深入了解，以及老師 對程序的講解，我成功的實現(xiàn)了正弦和余弦的傅里葉變換運 算，完成了實現(xiàn)要求。這次課程讓我對傅里葉變換有了更加 深入的了

2、解。在使用CCS軟件時，要注意加上文件，避免岀錯；還有就是文件要放在工程文件夾里，不然要出錯，還有就是 iiit和float之間的區(qū)分。我經(jīng)過此次的學習后，無論是傅里葉變換的原理，還是 DSP編程方面，都有了很大的提升。這對以后的課程學習提 供了堅實的基礎。我也從中獲益良多。參考文獻1 楊明廣.C程序設計教材.四川科學技術出版社,2012.2 鄒彥DSP原理及應用（修訂版）.電子工業(yè)岀版社,2002.3 吳京.信號分析與處理（修訂版）.電子工業(yè)出版社,2014.附件1正弦信號程序正弦信號文件* * * * * * * * * * * * * * * * * * */include o #mcl

3、ude o # mclude,HI#defiiie PI#define SAMPLENUMBER 128 void FFT();void IiiitFoiFFTO；void MakeWave();float INPUTSAMPLENUMBER,DATASAMPLENUMBER;float fVaveRSAMPLENUMBER,fWaveISAMPLENUMBER,wSAMPLENUMBER;float sin_tabSANIPLENUlVIBER,cos_tabSAMPLENUMBER;mam()int 1；IiutForFFTO；MakeWave();foi(i=0;iSAMPLENUMBE

4、R;i+ )fVaveRi=INPUTi;fVaveIi=;wi=;FFT(fWaveR,fVaveI);for(i=0;iSAMPLENUMBER;i+)DATAi=wi;wlule(l);void FFT(float dataRSAMPLENUMBER,float dataISAMPLENUMBER)int xO,x 1 ,x2,x3,x4,x5,x6,xx;mt i,j,k,b,p,L；float TR,TI,temp;for( i=O;iSAMPLENUMBER;i-H-)x0=xl=x2=x3=x4=x5=x6=0;x0=i&0x01; xl=(i/2)&0x01 ; x2=(i/4

5、)&0x01 ; x3=(i/8)&0x01; x4=(i/16)&0x01; x5=(i/32)&0x01;x6=(i/64)&0x01;xx=x0*64+xl *32+x2* 16+x3*8+x4*4+x5*2+x6; dataIxx=dataRi;for( i=O;i0 )b=b*2; i-;for(j=0;j0)p=p*2;i-;ip=p*i;for(k=j ;kv 12 8 ;k=k+2 *b)fITR=dataRk ;TI=dataIk ;temp=dataRk+b; dataRk=dataRk+dataRk+b*cos_tabp+dataIk+b*sm_tabp; datalk=

6、datalk -dataRk+b *sin_tabp+datalk+b *cos_tabp; dataRk+b=TR-dataRk+b *cos_tabp -datalk+b *sin_tab p;dataIk+b=TI+temp*sin_tabp-dataIk+b *cos_tabp;iI JfoT(i=0 ；1VS AMPLENUMBER/2 ;i+)wi=sqrt(dataRi *dataRi+dataIi *datali);void IiMtFoiFFT()lilt 1;foi(i=0;iSAMPLENUMBER;i+) sm_tabi=sin(PI*2*i/SAMPLENUMBER)

7、;cos_tabi=cos(PI*2*i/SAMPLENUMBER); void MakeWave()int 1；foi(i=0;iPROGlPAGE0.cinit: PROG1 PAGE 0.switch: PROGlPAGE0.vectois:VECTPAGE0.const: DARAM1PAGE 1.bss :DARAM1 PAGE 1.stack : DARAM2PAGE 1.system: DARAM2PAGE 1.data :DARAM2PAGE 1/* * * * * * * * * * *End of File* * * * * * * * */附件2余弦信號程序余弦信號文件/

8、* * * * * * * * * */#mclude o #mclude o # mcludeHU #define PI #ckflnc SAMPLENUMBER 128void FFT();void IiiitFoiFFTO；void MakeWave();float INPUTSAMPLENUMBER,DATASAMPLENUMBER; float fWavcRS AMPLENUMBER,fWaveISAMPLENUMBER,wSAMPLENUMBER; float sin_tabSAlVIPLENUlVIBER,cos_tabSAMPLENUMBER; inam()int 1；Iiut

9、FoiFFTO；MakeWave();for(i=0;iSAMPLENUMBER;i+ )fWaveRi=INPUTi;fVaveIi=;wi=;FFT(fWaveR,fVaveI); for(i=0;iSAMPLENUMBER;i+) DATAi=wi;wlule(l);void FFT(float dataRSAMPLENUMBER,float dataISAMPLENUMBER) mt x(xl.x2、x3,x4,x5、x6,xx;uit ij,k,bp，L;float TR.TI.temp;for( i=O;iSAMPLENUMBER;i+ )x0=xl=x2=x3=x4=x5=x6=

10、0;x0=i&0x01; xl=(i/2)&0x01 ; x2=(i/4)&0x01 ; x3=(i/8)&0x01; x4=(i/16)&0x01; x5=(i/32)&0x01;x6=(i/64)&0x01;xx=x0*64+xl *32+x2* 16+x3*8+x4*4+x5*2+x6; dataIxx=dataRi;for( i=O;i0 )b=b*2; i-;for(j=0;j0)p=p*2;i-;p=p*j；for(k=j ;kv 12 8 ;k=k+2 *b)TR=dataRk ;TI=dataIk ;temp=dataRk+b;dataRk=dataRk+dataRk+b *c

11、os_tabp+datalk+b *sin_tabp;datalk=datalk -dataRk+b *sin_tabp+datalk+b *cos_tabp;dataRk+b=TR-dataRk+b*cos_tabp-dataIk+b*sin_tabp;dataIk+b=TI+tcmp*sin_tabp -datalk+b *cos_tabp; for(i=0;iSAMPLENUMBER/2;i-l-+)wi=sqrt(dataRi *dataRi+dataIi *datali);void IiiitFoiFFT()lilt 1;foi(i=0;KSAMPLENUMBER;i+) sm_tabi=sin