快速傅立葉變換FFT實驗

1、快速傅立葉變換（FFT）算法實驗一 實驗?zāi)康? 加深對 DFT 算法原理和基本性質(zhì)的理解；2 熟悉 FFT 算法原理和 FFT 子程序的應(yīng)用；3 學(xué)習(xí)用 FFT 對連續(xù)信號和時域信號進行譜分析的方法， 了解可能出現(xiàn)的分析誤差及 其原因，以便在實際中正確應(yīng)用 FFT 。二 實驗設(shè)備計算機，CCS 3.1版軟件，E300實驗箱，DSP仿真器，導(dǎo)線三 基本原理1 離散傅立葉變換 DFT 的定義：將時域的采樣變換成頻域的周期性離散函數(shù)，頻域 的采樣也可以變換成時域的周期性離散函數(shù)，這樣的變換稱為離散傅立葉變換，簡 稱 DFT 。2. FFT是DFT的一種快速算法， 將DFT的N2步運算減少為（N/2

2、） log?”步，極大的 提高了運算的速度。3. 旋轉(zhuǎn)因子的變化規(guī)律。4. 蝶形運算規(guī)律。5. 基 2FFT 算法。四. 實驗步驟1. E300底板的開關(guān)SW4的第1位置ON，其余置OFF。其余開關(guān)不用具體設(shè)置。2. E300板子上的SW7開關(guān)的第1位置OFF，其余位置ON3. 閱讀本實驗所提供的樣例子程序；4. 運行 CCS 軟件，對樣例程序進行跟蹤，分析結(jié)果；記錄必要的參數(shù)。5. 填寫實驗報告。6. 提供樣例程序?qū)嶒灢僮髡f明A .實驗前準(zhǔn)備用導(dǎo)線連接“ Signal expansion Unit ”中 2 號孔接口“ SIN ”和“ A/D 單元”的 2 號孔接 口“ AD_IN0 ”。（

3、實驗采用的是外部的 AD 模塊）B.實驗1. 正確完成計算機、 DSP 仿真器和實驗箱的連接后，系統(tǒng)上電。2. 啟動 CCS3.1, Project/Open 打開 “ algorithm01_fft 子目錄下 “ fft.pjt 工程文件；雙擊“fft.pjt ” 及Source”可查看各源程序；加載 Debugfft.out ” ；3.單擊DebugGo main ”進入到主程序，在主程序flag=O ；處設(shè)置斷點;* /C5416 bSI via FDSSin PP E*ulat o/DSP_C5Jii - C54F - Cade Ca&posec Sludia - F:. J旦其4.

4、單擊“ Debug Run ”運行程序，或按 F5運行程序；程序?qū)⑦\行至斷點處停止;5. 用View / Graph / Time/Frequency 打開一個圖形觀察窗口；設(shè)置該觀察圖形窗口變量及參數(shù)；采用雙蹤觀察在啟始地址分別為px和pz，長度為128,數(shù)值類型為16位整型，px：存放經(jīng)A/D轉(zhuǎn)換后的輸入信號；pz :對該信號進行FFT變換的結(jié)果。6. 單擊“ Debug Animate 運行程序，或按F10運行；觀察窗口并觀察輸入信號波形及其FFT變換結(jié)果；注意：下圖觀察窗口中px波形在上面，pz波形在下面。7. 單擊“ Debug Halt暫停程序運行，關(guān)閉窗口，本實驗結(jié)束。實驗結(jié)果：

5、在CCS3.1環(huán)境，同步觀察輸入信號波形及其FFT變換結(jié)果；語音信號的FFT分析、實驗?zāi)康?. 加深對DFT算法原理和基本性質(zhì)的理解；2. 熟悉 FFT 算法原理和 FFT 子程序的應(yīng)用；3. 學(xué)習(xí)用 FFT 對連續(xù)信號和時域信號進行譜分析的方法， 了解可能出現(xiàn)的分析誤差及 其原因，以便在實際中正確應(yīng)用 FFT 。二、實驗設(shè)備計算機，CCS3.1版本軟件，DSP仿真器，E300實驗箱，54xpCPU板，音頻線，耳機三、實驗步驟1 、復(fù)習(xí) DFT 的定義、性質(zhì)和用 DFT 作譜分析的有關(guān)內(nèi)容；2、復(fù)習(xí) FFT 算法原理與編程思想，并對照 DIT-FFT 運算流程圖和程序框圖，了解本實 驗提供的

6、FFT 子程序；3、閱讀本實驗所提供的樣例子程序；4、運行 CCS 軟件，對樣例程序進行跟蹤，分析結(jié)果；記錄必要的參數(shù)。5、填寫實驗報告。1 提供樣例程序?qū)嶒灢僮髡f明四、實驗操作步驟 :1. E300板上的開關(guān)SW4的第1位置ON其余OFF.開關(guān)SW6全部撥到ON其余開關(guān)全部 置 OFF。2. 利用自備的音頻信號源，或把計算機當(dāng)成音頻，從E300板子上的音頻接口“ MIC IN ”輸入音頻信號,同時將耳機插到“ EARPHONE_OU孔內(nèi).3. 運行 Code Composer Studio (CCS3.1)， 進入后點擊“ DebugConnect ”；4. 用“ ProjectOpen ”

7、 打開系統(tǒng)項目文件19_audiofft fft.pjt ，編譯并裝載“audio_fft.out”，并播放音源信號，輸入到 “ MIC IN ”中.5. 單擊“ DebugGo main ”進入到主程序中，并在k+處設(shè)置斷點；單擊“ DebugRun運行程序，程序?qū)⑦\行至斷點處停止；6. 用View / Graph / Time/Frequency打開一個圖形觀察窗口；設(shè)置該觀察圖形窗口變量及參數(shù)；采用雙蹤觀察在啟始地址分別為 px和pz,長度為128的單元中數(shù)值的變化，數(shù)值類型為16位有符號整型變量，這兩段存儲單元中分別存放的是經(jīng)A/D轉(zhuǎn)換后的語音信號和對該信號進行 FFT變換的結(jié)果；（由

8、于音源不同，故波形也不同，下面僅供參考）7. 單擊“ Debug Animate ”運行程序，或按 F10運行；調(diào)整觀察窗口并觀察輸入信號波形 及其FFT變換結(jié)果；8. 單擊“ Debug Halt ”暫停程序運行，關(guān)閉窗口，本實驗結(jié)束。五、實驗報告要求1. 簡述實驗原理及目的；2. 結(jié)合實驗中所給定典型序列幅頻特性曲線，與理論結(jié)果比較，并分析說明誤差產(chǎn)生的原因以及用FFT作譜分析時有關(guān)參數(shù)的選擇方法。3總結(jié)實驗所得主要結(jié)論。六、程序參數(shù)說明void kfft(pr,pi,n,k,fr,fi,l,il):基2快速傅立葉變換子程序，n為變換點數(shù)，應(yīng)滿足 2的整數(shù)次幕，k為幕次(正整數(shù))；數(shù)組x :輸入信號數(shù)組，A/D轉(zhuǎn)換數(shù)據(jù)