一种分数阶傅里叶变换快速算法的研究.docVIP

一种分数阶傅里叶变换快速算法的研究 　　摘 要：介绍了分数阶傅里叶变换的定义，接着提出了一种分数阶傅里叶变换的快速算法，其中分数阶傅里叶变换快速算法分三步进行：线性调频信号乘法，线性调频信号卷积，另一个线性调频信号乘法,从而利用FFT来计算FRFT。这种算法思想直观，结果与连续FRFT的输出接近。最后用具体的信号作了计算机仿真，并给出Matlab仿真结果图。 　　?? 　　关键词：分数阶傅里叶变换；FFT；时频分析；卷积?? 　　中图分类号：TN911??7 文献标识码：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 分数阶傅里叶变换的定义?? 　　传统的傅里叶变换(FFT)对平稳信号的处理效果很好，但当信号频率随时间变化时，FFT就显得有些力不从心了。分数阶傅里叶变换(Fractional Fourier Transform,FRFT)可以很好地弥补FFT的不足，特别是处理线性调频信号(LFM)时，能够得到令人满意的结果。?? 　　FRFT也称为角度傅里叶变换(AFT) 或者旋转傅里叶变换(RFT),其定义式为: 　　?И? 　　Xp(u)=R??α［x(t)］=∫??∞??－∞??Kp(t,u)x(t)??d??t 　　(1) 　　?И? 　　式中变换核取作： 　　?お? 　　Kp(t,u)=1－??jcot?? α2??π?? 　　??e????j??(12u??2??cot?? α－ut??csc?? α+12t??2??cot?? α)??, α≠n??π???? 　　δ(t－u), α=2n??π???? 　　δ(t+u), α=(2n+1)??π?? 　　其中??n为整数，即n∈Z。α=p??π??/2称为分数阶??Fourier??变换的阶数，并有R??α=R????pπ2??=F??p。Kp(t,u)称为??FRFT??的核函数。Xp(u)称为x(t)的p阶??Fourier??变换。??FRFT??是一种线性算子,记为Fp,他满足以下性质： ?И? 　　(1) FRFT变换为线性算子；?? 　　(2) ??F??0［x(t)］=F??4［x(t)］=x(t)??(恒等变换)；?? 　　(3) ??F??1［x(t)］=F??5［x(t)］=X(ω)??(标准Fourier变换)；?? 　　(4) 广义Fourier变换算子为加性算子，即有??F????p+q??=F??pF??q。?И?