DSPFFT試驗(yàn)報(bào)告

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

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

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