source1.m

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

%gaussian pulse source
function src = source1(step,dt)

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

t = step * dt;

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

end