function?Fa=DFPei（fa）

%?Compute?discrete?fractional?Fourier?transform
%?of?order?a?of?vector?f?according?to?Pei/Yeh/Tseng

N=length（f）;?f=f（:）;
shft?=?rem（（0:N-1）+fix（N/2）N）+1;

global?hn_saved?p_saved
if?（nargin==2）?p?=?2;?end;
p?=?min（max（2p）N-1）;

if?（length（hn_saved）?~=?N?|?p_saved?~=?p）
????hn?=?make_hn（Np）;
????hn_saved?=?hn;?p_saved?=?p;
else
????hn?=?hn_saved;
end;
Fa（shft1）=hn*（exp（-j*pi/2*a*[0:N-1]）.‘.*（hn‘*f（shft）））;

function?hn?=?make_hn（Np）
even?=?rem（N2）==0;
shft?=?rem（（0:N-1）+fix（N/2）N）+1;
%?Gauss-Hermite?samples
u?=?（-N/2:（N/2）-1）‘/sqrt（N/2/pi）;
ex?=?exp（-u.^2/2）;?
hn（:1）?=?ex;?r?=?norm（hn（:1））;?hn（:1）?=?hn（:1）/r;
hn（:2）=2*u.*ex;?s?=?norm（hn（:2））;?hn（:2）?=?hn（:2）/s;
r?=?s/r;
for?k?=?3:N+1
?????hn（:k）=2*r*u.*hn（:k-1）-2*（k-

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-07-25?16:30??FRFT\
?????文件????????2646??2008-07-29?18:29??FRFT\cdpei.m
?????文件??????502153??2008-07-29?18:29??FRFT\Computation?of?the?Fractional?Fourier?Transform.pdf
?????文件????????1522??2008-07-29?18:30??FRFT\Disfrct.m
?????文件????????2130??2008-07-29?18:29??FRFT\Disfrft.m
?????文件????????1506??2008-07-29?18:30??FRFT\Disfrst.m
?????文件????????1514??2008-07-29?18:29??FRFT\frft.m
?????文件????????1200??2008-07-29?18:30??FRFT\rescale.m
?????文件?????????789??2008-07-29?18:30??FRFT\testfrct.m
?????文件?????????364??2008-07-29?18:30??FRFT\testfrft.m
?????文件??????545336??2010-12-06?21:37??FRFT\分數階傅里葉變換.pdf
?????文件??????301172??2011-05-26?00:37??FRFT\分數階傅里葉變換域上帶通信號的采樣定理.pdf
?????文件??????440745??2010-12-06?21:31??FRFT\分數階傅里葉變換的快速計算新方法.pdf
?????文件??????208288??2010-12-13?22:58??FRFT\基于分數階Fourier變換的雷達目標識別.pdf
?????文件??????249617??2010-12-06?21:38??FRFT\基于分數階傅里葉變換的正交頻分復用系統同步分析.pdf
?????文件??????348998??2010-12-06?21:30??FRFT\基于分數階傅里葉變換的線性調頻信號二維波達方向估計.pdf
?????文件?????3411870??2010-12-06?21:31??FRFT\基于變標處理和分數階傅里葉變換的運動目標檢測算法.pdf
?????文件?????3136742??2010-12-06?21:06??FRFT\小波分析與分數傅里葉變換及應用.pdf