【正文】 //向下權(quán)值依次比上級減半 } n=n+k。 //最高位權(quán)值 while(n=k) { n=nk。xi[n]=T。 T=xi[m]。xr[m]=xr[n]。m=N12。 //最高位加1相當于十進制加上最高位的權(quán)N1/2 n=N2。 //m為正序列；n為到序列；k為各個權(quán)值；//N2為最高位的權(quán)值 float T。i++){ INPUT[i]=sin(PI*2*i/N*3)*1024。for(i=0。 cos_tab[i]=cos(PI*2*i/N)。iN。}void InitForFFT() //FFT初始化函數(shù)，建立正余弦函數(shù)表{int i。i++)DATA[i]=w[i]。for(i=0。}Mum=(int)(+log(N)/log(2))。 fWaveI[i]=。iN。 //FFT初始化函數(shù)MakeWave()。 //計算功率普}}main(){int i。mN/2。 Xi[k]=Xi[k]+Y。 Xi[k+B]=Xi[k]Y。 Y=Xi[k+B]*cos_tab[S]Xr[k+B]*sin_tab[S]。k=N1。j++) //每級需要進行B種蝶形運算 { S=j*(int)pow(2,Mumm)+。 //B=2^(m1) for(j=0。m=Mum。 finv(N,Xr,Xi)。 //S為旋轉(zhuǎn)因子的冪數(shù)，B為蝶形圖運算輸入數(shù)據(jù)的距離， //也即各級旋轉(zhuǎn)因子的個數(shù) int m,j,k。 //正余弦函數(shù)int Mum。float fWaveR[N],fWaveI[N],w[N]。 //波形發(fā)生函數(shù)void finv(int N1,float *xr,float *xi)。而這次，呂老師不但要帶著我們搞課設(shè)，還要指導(dǎo)那些參加飛思卡爾比賽的同學(xué)，呂老師這一周過的真的很累，老師，您辛苦了在這里我也同樣要感謝雷老師在最后一天對我們的答辯，因為呂老師臨時有事，故雷老師給我答辯，這本不是雷老師的事，雷老師本可對我們寬松些，這樣自己也會輕松許多，然而雷老師卻非常的認真，總是問的詳詳細細的，好多同學(xué)都被問的無言以對，而且因為大家都急著答辯，害的雷老師很晚才去吃午飯，炎熱的夏天，悶熱的教室，豆大的汗珠順著雷老師的臉頰流下來，在次，我怎么忘記對雷老師的感謝同時，也感謝和自己做相同課設(shè)題目的黃玉亮、石菊等人，我每當越到不懂的地方都是先和大家一切探討，所以，這次課設(shè)的成功同樣離不開我的同學(xué)們，謝謝你們在此，也感謝圖書館為我們提供便利的條件查詢資料，謝謝參考文獻[1] 俞一彪，孫兵. ：東南大學(xué)出版社，2005[2] 李行一. ：重慶大學(xué)出版社，2002[3] 、結(jié)構(gòu)及應(yīng)用基礎(chǔ)——TMS320F28x[M].北京：機械工業(yè)出版社，2007.[4]姜沫岐，許涵，俞鵬，段國強 . DSP原理與應(yīng)用—從入門到提高. 北京: 機械工業(yè)出 版社， 2007[5] 俞卞章. 數(shù)字信號處理. 西安: 西北工業(yè)大學(xué)出版社， 2002[6] 附錄A1 程序清單include //數(shù)字函數(shù)頭文件define PI define N 128 //采樣次數(shù)Nvoid InitForFFT()。一周的課設(shè)結(jié)束了，這次的課設(shè)意義又非比尋常，課設(shè)一結(jié)束，我們大三的生活也即將結(jié)束，我們也將進入找工作的階段，我們原本迷茫、原本惆悵，可一想到呂老師那充滿鼓舞的微笑，我們也充滿了動力，一切困難我們都會努力的解決致謝春蠶到死絲方盡，蠟炬成灰淚始干，我用這兩句來說我們的呂老師實在是在恰當不過了，在這一周的課設(shè)中，呂老師每天都要到實驗室來看看同學(xué)們的設(shè)計過程，如果哪位同學(xué)有什么不能解決的問題，呂老師都會不厭其煩的引導(dǎo)，給出思路。原本我對老師有種懼怕或者說抵觸，可是呂老師讓我的這種感覺不見了，甚至很想和呂老師多多交流、溝通。比較仿真圖，可以對設(shè)計進行更深入的學(xué)習(xí)，驗證了程序設(shè)計的正確性。在設(shè)計部分，成功的在CCS環(huán)境下檢測程序，做出軟件仿真波形圖：正弦波輸入波形、PC機FFT計算處理后的仿真波形圖。通過這次的課設(shè)，我更加熟悉了CCS開發(fā)的環(huán)境，它提供了環(huán)境配置、源文件編輯、程序調(diào)試、跟蹤和分析等工具，幫助用戶在一個軟件環(huán)境下完成編輯、編譯鏈接、調(diào)試和數(shù)據(jù)分析等工作。功率信號的功率譜反應(yīng)了信號功率隨頻率分布的特點，功率譜是信號先自相關(guān)再作FFT變換小結(jié) 單片機、DSP、智能儀器、ARM是我們專業(yè)的幾門主要課程，都涉及到編程的知識，雖然我們以前也做過單片機的不少實驗和課設(shè)，但是以前都是用匯編編的，但這次的DSP課設(shè)卻用C語言編程，而且DSP中有很多指令單片機中都沒有，這些都給我?guī)砹瞬簧俚穆闊?。利用此系?shù)表可以完成8～1024點的FFT運算。正弦系數(shù)表和余弦系數(shù)表由文件給出， 匯編命令將正弦系數(shù)和余弦系數(shù)與程序代碼匯編在一起。而旋轉(zhuǎn)因子又是一復(fù)數(shù)運算，因此在程序中設(shè)置了正弦表和余弦表，通過查這兩個表來計算旋轉(zhuǎn)因子。在第三級到最后一級蝶形運算中，組內(nèi)蝶形運算次數(shù)逐漸增多，而且多是乘法運算，單獨處理是不可取的，因此這段運算采用相同的程序塊來實現(xiàn)。 FFT計算當2=256時，128點復(fù)數(shù)FFT運算過程分成七級、3個階段實現(xiàn)，這三個階段是：第一級蝶形運算，第二級蝶形運算，第三級～第級蝶形運算，以下分別分析這三個階段蝶形運算的工作過程。 位倒序子程序在編寫位倒序子程序時，使用提供的位倒序?qū)ぶ贩绞剑梢苑奖愕膶崿F(xiàn)位倒序操作。正弦波形：FFT變換之后波形： 方波： FFT變換之后的波形： 三角波： FFT變換之后的波形： 7 工作過程分析 程序的初始化128點實數(shù)FFT的程序主要由4部分組成，分別是位反轉(zhuǎn)子程序，F(xiàn)FT核心計算子程序，奇偶分離及結(jié)果產(chǎn)生子程序，還有功率譜計算子程序。并在原來程序的基礎(chǔ)上加入了指示燈的信號程序：當程序進入到FFT變換時，黃色指示燈亮；執(zhí)行完或執(zhí)行之前，紅色指示燈亮。都采用16bit signed integer。其中輸入在single time下完成。(如上圖)面的例子展示了一個單時刻(時間/頻率)圖表屬性對話框點擊，一旦屬性被設(shè)定了，選擇OK按鈕，開啟圖窗口，就會畫出指定的數(shù)據(jù)點根據(jù)芯片和所編寫的程序，我做了兩個個仿真圖：方波輸入波形、PC機FFT計算處理后的仿真波形圖。選擇View Graph，可以得到這些圖，然后選擇需要的圖。點擊Open加載。加載過程是將上述構(gòu)建成功，生成的可執(zhí)行文件加載到目標板，目標板可以是軟件仿真環(huán)境，也可以是硬件目標板。按照以下步驟創(chuàng)建一個新工程并運行程序：（1） 建一個文本文檔，在里面寫下程序，（2） , 選擇ProjectNew（3） 在Project Name框里輸入工程名：(例如:1)（4） 在Location框里輸入或者瀏覽