分析不同数据在不同频率域、频率波数域的特征

分析二维数据在频率-空间域、频率-波数域、拉东域的特征，适用于地震、雷达、超声数据的分析

function imag_fk_spectrum（ddtdxffkk）
% ff 频率范围 [f_min f_max]
% kk  波数范围

N=size（d1）; %t
M=size（d2）; %x   %道数
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）;  %% 等价于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））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         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