資源簡介
分析二維數據在頻率-空間域、頻率-波數域、拉東域的特征,適用于地震、雷達、超聲數據的分析
代碼片段和文件信息
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
- 上一篇:數值分析 5套卷
- 下一篇:電力電子課程設計 直流斬波電路的設計與仿真
評論
共有 條評論