基于dsp的fft實現設計報告-資料下載頁
2025-06-20 12:31本頁面
【正文】 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 */七、實驗過程1．編譯并下載程序。2．打開觀察窗口：*選擇菜單ViewGraphTime/Frequency…進行如下圖所示設置。圖1圖2圖33．清除顯示：在以上打開的窗口中單擊鼠標右鍵，選擇彈出式菜單中“Clear Display”功能。4．設置斷點： 中有注釋“break point”的語句上設置軟件斷點。圖45．運行并觀察結果。⑴ 選擇“Debug”菜單的“Animate”項，或按Alt+F5 鍵運行程序。⑵ 觀察“Test Wave”窗口中時域圖形；圖5⑶ 在“Test Wave”窗口中點擊右鍵，選擇屬性，更改圖形顯示為FFT。觀察頻域圖形。圖6⑷ 觀察“FFT”窗口中的由CCS 計算出的正弦波的FFT。圖7(5)改變輸入函數 INPUT[i]=(sin(PI*2*i/SAMPLENUMBER*3)+sin(PI*2*i/SAMPLENUMBER*4+sin(PI*2*i/SAMPLENUMBER*8)))*1024。圖8(5)改變輸入函數 INPUT[i]=(sin(PI*2*i/SAMPLENUMBER*3)+sin(PI*2*i/SAMPLENUMBER*10)+sin(PI*2*i/SAMPLENUMBER*20)))*1024。圖9八、實驗結果通過觀察頻域和時域圖，程序計算出了測試波形的功率譜，與CCS 計算的FFT 結果相近。九、問題與思考 (1)觀察圖6和圖7，可以看到二者波形相似，但橫縱坐標均不相同，縱坐標大約是二倍的關系，橫坐標大約為142倍。 (2)觀察圖8，因為兩個頻率比較相近，因此出現了前兩個頻譜交疊的現象。十、心得體會通過這次DSP課程設計，熟悉了FFT的算法原理和FFT子程序的算法流程和應用，掌握了DSP中FFT的設計和編程思想，以及用FFT對連續(xù)信號和時域信號進行頻譜分析的方法，和使用CCS的波形觀察器觀察波形和頻譜情況，收獲頗豐。
環(huán)評公示相關推薦
鄂ICP備17016276號-1