fk.m

来自「著名的seismiclab的代码 是地震学研究人员必备的工具」· M 代码 · 共 21 行

function [S,f,k] = fk(d,dt,dx);

[nt,nx]=size(d);

nk = 2*(2^nextpow2(nx));
nf = 2*(2^nextpow2(nt));
S = fftshift(fft2(d,nf,nk));

S = abs(S(nf/2+1:nf,:));

f = (0:1:nf/2-1)/nf/dt;
k = (-nk/2+1:1:nk/2)/nk/dx;

Smax = max(max(S));

Smin = min(min(S));



S = (S-Smin)/(Smax-Smin):