實驗一用FFT進行信號譜分析

1、實驗一用FFT進行信號譜分析一. 實驗?zāi)康模篈. 進一步加深DFT算法原理和基本性質(zhì)的理解（因為FFT只是DFT的一種快速算法，所以FFT的運算結(jié)果必然使DFT得基本性質(zhì)）。B. 熟悉FFT算法原理和FFT子程序的應(yīng)用。C. 學(xué)習(xí)用FFT對連續(xù)信號和時域離散信號進行譜分析的算法，了解可能出現(xiàn)的分析誤差及其原因，以便在實際中正確應(yīng)用FFT。二. 實驗原理：DFT是離散傅里葉級數(shù)（DFS）在時域和頻域內(nèi)的兩個主值序列構(gòu)成的變換對：Xk=DFTxn=n=0N-1xn 0kN-1xn=IDFTXk=1Nk=0N-1XkWn-kn 0kN-1用DFT來分析有限長序列的頻譜特性。FFT是DFT的快速算法，

2、常用的兩種快速算法有按時間抽取的FFT算法和按頻率抽取的FFT算法，其中按時間抽取的FFT算法如下式：Xk=X1k+WNkX2k k=0,1,2N2-1Xk+N2=X1k-WNkX2k k=0,1,2N2-1按頻率抽取的FFT算法如下：Xk=n=0N2-1xn+-1kxn+N2WNk k=0,1,2N-1三. 實驗數(shù)據(jù)1、x1(n)及其8點和16點DFT2、x2(n)及其8點和16點DFT3、x3(n)及其8點和16點DFT4、x4(n)的8點和16點波形及其DFT5、x5(n)的8點和16點波形及其DFT6、x6(n)的16點、32點和64點采樣序列波形及其DFT選7時，計算并圖示x7n=x

3、4n+x5n*R8(n)和x7n=x4n+x5n*R16(n)及其DFT。程序自動計算并繪圖驗證DFT的共軛對稱性。當(dāng)N=16時， x4n=x4N-n，x5n=x5N-n。即x4n為x7n的共軛對稱分量，而x5n是x7n的共軛反對稱分量。根據(jù)DFT的共軛對稱性，應(yīng)有以下結(jié)果：x7k=DFTx7n16點=ReX7k+jImX7kx7n=x4n+x5n的8點和16點波形及其DFTx4k=DFTx4n16點=ReX7kx5k=DFTx5n16點=jImX7k繪出ReX7k和jImX7k的模。它們正是圖中16點的|X4k|和|X5k|。選8時，計算并圖示x8n=x4n+jx5n*R8(n)和x8n=x

4、4n+jx5n*R16(n)及其DFT。程序自動計算并繪圖驗證DFT的共軛對稱性的第二種形式：如果xn=xrn+jxin， Xk=DFTxn=Xepk+Xopk，則Xepk=DFTxrn，Xopk=DFTjxin。其中Xepk=12Xk+X*N-k， Xopk=12Xk-X*N-k。x8n的8點和16點DFT程序計算結(jié)果如下：X8ek=12|X8k+X8*N-k|及x4n=x8rn=IDFTX8ek，正好與圖中x4n的16點|x4n|及x4n相同。jX8ok=12|X8k-X8*N-k|及x5n=x8in=IDFTX8ekj，正好與圖中16點的|x5n|及x5n相同。四. 思考題（1）當(dāng)N=8

5、時，x1n和x2n的幅頻特性會相同嗎？為什么？N=16 呢？（2）對于周期序列，如果周期不知道，如何用FFT進行譜分析？ 答：（1）當(dāng)N=8時，x1n和x2n的幅頻特性會相同. 當(dāng)N=16時，x1n和x2n的幅頻特性會不相同。（2）、如果xn的周期預(yù)先不知道，可截取M點進行DFT，即xMn=xnRMn XMn=DFTxMn 0kM-1 再將截取長度擴大1倍，截取x2Mn=R2Mn X2Mn=DFTx2Mn 0k2M-1 比較xMk和x2Mk，如果兩者的主譜差別滿足分析誤差要求，則以xMk或x2Mk近似表示 xn的頻譜，否則，繼續(xù)將截取長度加倍，直至前后兩次分析所得主譜頻率差別滿足誤差要求。設(shè)最

6、后截取長度為iM，則xiMk0表示=2iMk0點的譜線強度。通過實驗，我知道了用FFT對信號作頻譜分析是學(xué)習(xí)數(shù)字信號處理的重要內(nèi)容。經(jīng)常需要進行譜分析的信號是模擬信號和時域離散信號。對信號進行譜分析的重要問題是頻譜分辨率D和分析誤差。頻譜分辨率直接和FFT的變換區(qū)間N有關(guān)，因為FFT能夠?qū)崿F(xiàn)的頻率分辨率是2ND。可以根據(jù)此式選擇FFT的變換區(qū)間N。誤差主要來自于用FFT作頻譜分析時，得到的是離散譜，而信號（周期信號除外）是連續(xù)譜，只有當(dāng)N較大時，離散譜的包絡(luò)才能逼近于連續(xù)譜，因此N要適當(dāng)選擇大一些。周期信號的頻譜是離散譜，只有用整數(shù)倍周期的長度作FFT，得到的離散譜才能代表周期信號的頻譜。如果不知道信號周期，可以盡量