source2.m

时间域有限差分法（TD-FDM）源代码

%time derivative of gaussian pulse

function src = source2(step,dt)

sigma = 1/(16*pi);
mean = 4*sigma; %offset

t = step * dt;
a = -(t-mean)^2 / (2*sigma^2);

src = -(1/(sqrt(2*pi))) * ((t-mean)/sigma^3) * exp(a);

end