基于dsp的fft實(shí)現(xiàn)設(shè)計(jì)報(bào)告-閱讀頁(yè)
2025-07-05 12:31本頁(yè)面
【正文】 aveI[SAMPLENUMBER],w[SAMPLENUMBER]。void InitForFFT(){int i。iSAMPLENUMBER。 cos_tab[i]=cos(PI*2*i/SAMPLENUMBER)。 for ( i=0。i++ ) { INPUT[i]=sin(PI*2*i/SAMPLENUMBER*3)*1024。 InitForFFT()。 for ( i=0。i++ ) { fWaveR[i]=INPUT[i]。 w[i]=。 for ( i=0。i++ ) { DATA[i]=w[i]。 // break point}void FFT(float dataR[SAMPLENUMBER],float dataI[SAMPLENUMBER]){ int x0,x1,x2,x3,x4,x5,x6,xx。 float TR,TI,temp。iSAMPLENUMBER。 x0=iamp。 x1=(i/2)amp。 x2=(i/4)amp。 x3=(i/8)amp。x4=(i/16)amp。 x5=(i/32)amp。 x6=(i/64)amp。 xx=x0*64+x1*32+x2*16+x3*8+x4*4+x5*2+x6。 } for ( i=0。i++ ) { dataR[i]=dataI[i]。 } for ( L=1。L++ ) { /* for(1) */ b=1。 while ( i0 ) { b=b*2。 } /* b= 2^(L1) */ for ( j=0。j++ ) /* for (2) */ {p=1。 while ( i0 ) /* p=pow(2,7L)*j。 i。 for ( k=j。k=k+2*b ) /* for (3) */ { TR=dataR[k]。 temp=dataR[k+b]。 dataI[k]=dataI[k]dataR[k+b]*sin_tab[p]+dataI[k+b]*cos_tab[p]。 dataI[k+b]=TI+temp*sin_tab[p]dataI[k+b]*cos_tab[p]。iSAMPLENUMBER/2。 }} /* END FFT */七、實(shí)驗(yàn)過(guò)程1．編譯并下載程序。圖1圖2圖33．清除顯示：在以上打開(kāi)的窗口中單擊鼠標(biāo)右鍵，選擇彈出式菜單中“Clear Display”功能。圖45．運(yùn)行并觀察結(jié)果。⑵ 觀察“Test Wave”窗口中時(shí)域圖形；圖5⑶ 在“Test Wave”窗口中點(diǎn)擊右鍵，選擇屬性，更改圖形顯示為FFT。圖6⑷ 觀察“FFT”窗口中的由CCS 計(jì)算出的正弦波的FFT。圖8(5)改變輸入函數(shù) INPUT[i]=(sin(PI*2*i/SAMPLENUMBER*3)+sin(PI*2*i/SAMPLENUMBER*10)+sin(PI*2*i/SAMPLENUMBER*20)))*1024。九、問(wèn)題與思考 (1)觀察圖6和圖7，可以看到二者波形相似，但橫縱坐標(biāo)均不相同，縱坐標(biāo)大約是二倍的關(guān)系，橫坐標(biāo)大約為142倍。十、心得體會(huì)通過(guò)這次DSP課程設(shè)計(jì)，熟悉了FFT的算法原理和FFT子程序的算法流程和應(yīng)用，掌握了DSP中FFT的設(shè)計(jì)和編程思想，以及用FFT對(duì)連續(xù)信號(hào)和時(shí)域信號(hào)進(jìn)行頻譜分析的方法，和使用CCS的波形觀察器觀察波形和頻譜情況，收獲頗豐
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1