ricker.m

来自「这是matlab在地球物理数据处理方面的源码」· M 代码 · 共 21 行

function w = ricker(f,dt)
% RICKER(f,dt)
% Returns a Ricker wavelet of
% central freq. f.   
% dt is the sampling interval


% Mauricio D. Sacchi
% Jan/1997.
% [e-mail]

% f << 1/(2dt)....

nw=4./f/dt;
nw=2*floor(nw/2)+1;
nc=floor(nw/2);
i=1:nw;
alpha=(nc-i+1).*f*dt;
beta=alpha.^2;
w=(1.-beta.*2).*exp(-beta);