function?imag_fk_spectrum（ddtdxffkk）
%?ff?頻率范圍?[f_min?f_max]
%?kk??波數(shù)范圍

N=size（d1）;?%t
M=size（d2）;?%x???%道數(shù)
NFFT=2^nextpow2（N）;?%?Next?power?of?2?from?length?of?N
MFFT=2^nextpow2（M）;
tmax=（NFFT-1）*dt;????
t=0:dt:tmax;
ntpad=NFFT;??
nxpad=MFFT;
d?=?[d;zeros（ntpad-NM）];
d?=?[d?zeros（NFFTnxpad-M）];
[specfxfx]=?fftrl（dt0ntpad）;??%t-x域變換到f-x域
specfk=ifft（specfx[]2）;??%%?等價(jià)于spec1?=?ifft（specfx1.‘MFFT）.‘;??%將f-x域變換到f-k域

[rp?cp]=size（specfk）;
kx=fftshift（1/2/cp/dx*（-cp:2:cp-2））;


specfk_shift=fftshift（specfk2）;

rev_kx=fftshift（kx）;

f_ind=near（fxff（1）ff（2））;
k_ind=near（rev_kxkk（1）kk（2））;

specfk_cut=specfk_shift（f_indk_ind）;

imagesc（rev_kx（k_ind）fx（f_ind）abs（specfk_cut））;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件?????????763??2017-10-13?10:43??imag_fk_spectrum.m
?????文件?????????371??2017-10-13?10:43??imag_spectrum.m
?????文件????????2234??2017-11-02?20:43??smooth_spectrum.m