source1.m

来自「时间域有限差分法（TD-FDM）源代码」· M 代码 · 共 11 行

%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