一種分?jǐn)?shù)階傅里葉變換快速算法的研究
2008-04-12 00:00:00黃瓊玲劉振興尉宇
現(xiàn)代電子技術(shù) 2008年9期
摘 要：介紹了分?jǐn)?shù)階傅里葉變換的定義，接著提出了一種分?jǐn)?shù)階傅里葉變換的快速算法，其中分?jǐn)?shù)階傅里葉變換快速算法分三步進(jìn)行：線性調(diào)頻信號乘法，線性調(diào)頻信號卷積，另一個線性調(diào)頻信號乘法，從而利用FFT來計算FRFT。這種算法思想直觀，結(jié)果與連續(xù)FRFT的輸出接近。最后用具體的信號作了計算機仿真，并給出Matlab仿真結(jié)果圖。
關(guān)鍵詞：分?jǐn)?shù)階傅里葉變換；FFT；時頻分析；卷積
中圖分類號：TN9117 文獻(xiàn)標(biāo)識碼：A
文章編號：1004-373X(2008)09-156-02
Research on Fast Algorithm for Fractional Fourier Transform
HUANG Qiongling，LIU Zhenxing，WEI Yu
(Department of Information，Wuhan University of Science and Technology，Wuhan，430081，China)
Abstract：The definition of the Fractional Fourier Transform (FRFT) is presented in the paper.A new algorithm for efficient and accurate computation of FRFT is given.The new algorithm of FRFT includes three steps:The multiplication of linear frequency modulation signal;the convolve of linear frequency modulation signal;another multiplication of linear frequency modulation signal;so as tomake use of FFT to compute FRFT.This kind of calculate waykeeps a view and the output is close to the continuous FRFT.Finally，a few simulation results for some typical signals are provided to compare with previous ones by other methods in the end.
Keywords：fractional fourier transform;FFT;time-frequency analysis;convolve
1 分?jǐn)?shù)階傅里葉變換的定義
傳統(tǒng)的傅里葉變換(FFT)對平穩(wěn)信號的處理效果很好，但當(dāng)信號頻率隨時間變化時，F(xiàn)FT就顯得有些力不從心了。分?jǐn)?shù)階傅里葉變換(Fractional Fourier Transform，F(xiàn)RFT)可以很好地彌補FFT的不足，特別是處理線性調(diào)頻信號(LFM)時，能夠得到令人滿意的結(jié)果。
FRFT也稱為角度傅里葉變換(AFT) 或者旋轉(zhuǎn)傅里葉變換(RFT)，其定義式為:
Xp(u)=Rα［x(t)］=∫∞－∞Kp(t，u)x(t)dt
(1)
式中變換核取作：
Kp(t，u)=1－jcot α2π
ej(12u2cot α－utcsc α+12t2cot α)， α≠nπ
δ(t－u)， α=2nπ
δ(t+u)， α=(2n+1)π
其中n為整數(shù)，即n∈Z。α=pπ/2稱為分?jǐn)?shù)階Fourier變換的階數(shù)，并有Rα=Rpπ2=Fp。Kp(t，u)稱為FRFT的核函數(shù)。Xp(u)稱為x(t)的p階Fourier變換。FRFT是一種線性算子，記為Fp，他滿足以下性質(zhì)：
(1) FRFT變換為線性算子；
(2) F0［x(t)］=F4［x(t)］=x(t)(恒等變換)；
(3) F1［x(t)］=F5［x(t)］=X(ω)(標(biāo)準(zhǔn)Fourier變換)；
(4) 廣義Fourier變換算子為加性算子，即有Fp+q=FpFq。
2 采用分解方法計算FRFT的步驟
分?jǐn)?shù)階Fourier變換可以具體分解為以下三個主要的計算步驟：線性調(diào)頻信號乘法；線性調(diào)頻信號卷積；另一個線性調(diào)頻信號乘法。假定p∈［－1，1］，則我們可以將經(jīng)過量綱歸一化的信號f(x)的分?jǐn)?shù)階Fourier變換式(2)分解為以下三步運算：
fp(x)=e－jπx2tan(α/2)g′(x)
(3)
和:
即是說，分?jǐn)?shù)階Fourier變換的數(shù)值計算的順序如下：先計算式(式(5))，再計算式(4)，最后計算式(3)。下面是每一步計算的有關(guān)細(xì)節(jié)。
第一步：將函數(shù)f(x)與線性調(diào)頻函數(shù)相乘(式(5))。注意，g(x)的頻率帶寬與時間帶寬乘積可以是f(x)的相應(yīng)帶寬乘積的兩倍，所以要求g(x)的采樣間隔為1/(2Δx)。如果f(x)樣本值的采樣間隔是1/Δx，那么就需要對這些樣本值進(jìn)行插值，然后再與線性調(diào)頻函數(shù)的離散采樣值相乘，以得到所希望的g(x)的采樣。
第二步：將g(x)與一線性調(diào)頻函數(shù)作卷積式(式(4))。注意，由于g(x)是帶限信號，所以線性調(diào)頻函數(shù)也可以用其帶限形式代替而不會有任何影響。也就是說，我們可以取:
g′(x)=Aα∫∞－∞ejπβ(x－x′)2g(x′)dx′
=Aα∫∞－∞h(x－x′)g(x′)dx′
(7)
是函數(shù)ejπβx2的Fourier變換。于是，式(7)的離散形式為:
g′(m2Δx)=∑Nn=－Nhm－n2Δxgn2Δx
(10)
這一離散卷積可以利用快速Fourier變換計算。
第三步：計算式(3)得到f(x)的分?jǐn)?shù)階Fourier變換fp(x)的采樣值fpm2Δx。由于假定f(x)的所有變換都是帶限的，他們位于區(qū)間－12Δx，12Δx，所以需要用因子2對fpm2Δx進(jìn)行二抽一采樣，以得到離散采樣fpm2Δx。因此，對于不是π/2整數(shù)倍的角度，分?jǐn)?shù)階Fourier變換的計算對應(yīng)以下步驟：
(1) 原信號與一線性調(diào)頻函數(shù)相乘；
(2) Fourier變換(其變元乘以尺度系數(shù)csc α)；
(3) 再與一線性調(diào)頻函數(shù)相乘；
(4) 乘以一復(fù)幅值因子。
3 對信號的FRFT處理及仿真圖
首先要給出一個輸入信號x，然后根據(jù)分解方法編出FRFT快速算法，根據(jù)不同的p值輸入信號x會生成不同的曲線，其中p∈(0，4)。找出每個p值時FRFT的輸出最大值點，組成一個一維數(shù)組m1，再從m1中找出一個最大值，該值所對應(yīng)的p值就是FRFT變換的最佳角度α=pπ/2。
由圖1(d)中可以看出當(dāng)p=1時，α=π/2就是普通的傅里葉變換，這也驗證了分?jǐn)?shù)階傅里葉變換的正確性。
圖2 chirp信號隨著p變化的FRFT變換仿真結(jié)果
由圖2可以看出，在圖2(e)中p=15的時候，即α=3π/4時形成了一個沖擊信號，說明了在此角度上信號的能量最好地集聚在一點上，由此可以識別出信號的調(diào)頻系數(shù)，檢測出信號的參數(shù)，這就是FRFT處理LFM信號的顯著作用。
4 結(jié) 語
分?jǐn)?shù)階傅里葉變換是近二十年來發(fā)展起來的一種全新的信號時頻分析工具，在很多方面得到了十分廣泛的應(yīng)用。而其快速算法的研究則對擴展其應(yīng)用領(lǐng)域有著十分重要的意義。本文提出了一種有效并能準(zhǔn)確計算FRFT 的新算法。該算法具有易實現(xiàn)、易理解、精度較高等優(yōu)點，相信FRFT將會受到更廣泛的重視，在信號處理領(lǐng)域會有良好的應(yīng)用前景。
參 考 文 獻(xiàn)
［1］平先軍，陶然，周思永，等.一種新的分?jǐn)?shù)階傅里葉變換快速算法\\[J\\].電子學(xué)報，2001，29(3):406-408.
［2］Lufs B Almeida.The Fractonal Fourier Transform and Time-Frequency Representations ［J］.IEEE Transactions on Signal Processing，1994，42(11).
［3］尉宇.線性調(diào)頻和非線性調(diào)頻信號的檢測與參數(shù)估計\\[D\\].武漢：華中科技大學(xué)，2005.
作者簡介 黃瓊玲 1982年出生，碩士研究生。主要研究方向為數(shù)字信號處理。
注：本文中所涉及到的圖表、注解、公式等內(nèi)容請以PDF格式閱讀原文。
