基于dsp的fft實(shí)現(xiàn)設(shè)計(jì)報(bào)告-資料下載頁(yè)
2025-06-20 12:31本頁(yè)面
【正文】 b*2。 i。 } /* b= 2^(L1) */ for ( j=0。j=b1。j++ ) /* for (2) */ {p=1。 i=7L。 while ( i0 ) /* p=pow(2,7L)*j。 */ {p=p*2。 i。 } p=p*j。 for ( k=j。k128。k=k+2*b ) /* for (3) */ { TR=dataR[k]。 TI=dataI[k]。 temp=dataR[k+b]。 dataR[k]=dataR[k]+dataR[k+b]*cos_tab[p]+dataI[k+b]*sin_tab[p]。 dataI[k]=dataI[k]dataR[k+b]*sin_tab[p]+dataI[k+b]*cos_tab[p]。 dataR[k+b]=TRdataR[k+b]*cos_tab[p]dataI[k+b]*sin_tab[p]。 dataI[k+b]=TI+temp*sin_tab[p]dataI[k+b]*cos_tab[p]。 } /* END for (3) */ } /* END for (2) */ } /* END for (1) */ for ( i=0。iSAMPLENUMBER/2。i++ ) { w[i]=sqrt(dataR[i]*dataR[i]+dataI[i]*dataI[i])。 }} /* END FFT */七、實(shí)驗(yàn)過(guò)程1．編譯并下載程序。2．打開觀察窗口：*選擇菜單ViewGraphTime/Frequency…進(jìn)行如下圖所示設(shè)置。圖1圖2圖33．清除顯示：在以上打開的窗口中單擊鼠標(biāo)右鍵，選擇彈出式菜單中“Clear Display”功能。4．設(shè)置斷點(diǎn)： 中有注釋“break point”的語(yǔ)句上設(shè)置軟件斷點(diǎn)。圖45．運(yùn)行并觀察結(jié)果。⑴ 選擇“Debug”菜單的“Animate”項(xiàng)，或按Alt+F5 鍵運(yùn)行程序。⑵ 觀察“Test Wave”窗口中時(shí)域圖形；圖5⑶ 在“Test Wave”窗口中點(diǎn)擊右鍵，選擇屬性，更改圖形顯示為FFT。觀察頻域圖形。圖6⑷ 觀察“FFT”窗口中的由CCS 計(jì)算出的正弦波的FFT。圖7(5)改變輸入函數(shù) INPUT[i]=(sin(PI*2*i/SAMPLENUMBER*3)+sin(PI*2*i/SAMPLENUMBER*4+sin(PI*2*i/SAMPLENUMBER*8)))*1024。圖8(5)改變輸入函數(shù) INPUT[i]=(sin(PI*2*i/SAMPLENUMBER*3)+sin(PI*2*i/SAMPLENUMBER*10)+sin(PI*2*i/SAMPLENUMBER*20)))*1024。圖9八、實(shí)驗(yàn)結(jié)果通過(guò)觀察頻域和時(shí)域圖，程序計(jì)算出了測(cè)試波形的功率譜，與CCS 計(jì)算的FFT 結(jié)果相近。九、問(wèn)題與思考 (1)觀察圖6和圖7，可以看到二者波形相似，但橫縱坐標(biāo)均不相同，縱坐標(biāo)大約是二倍的關(guān)系，橫坐標(biāo)大約為142倍。 (2)觀察圖8，因?yàn)閮蓚€(gè)頻率比較相近，因此出現(xiàn)了前兩個(gè)頻譜交疊的現(xiàn)象。十、心得體會(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