With some basic indexing, rescaling of the gaussian length and line equation maths, we can do this by:
L=5;
x=1:10;
t=(1:10)/10;
[X,Y]=meshgrid(linspace(0,x(end)+L,500),linspace(0,t(end),500));
%slope of both lines
m=(t(2)-t(1))/(x(2)-x(1));
for ii=1:size(X,1)
z(ii,:)=exp(-((X(ii,:)-Y(ii,1)/m)/(L/2)).^2);
end
% basic line maths an dindexing
z(Y./X>m)=NaN;
z((Y+m*L)./X<m)=NaN;
surf(X,Y,z,'linestyle','none')
view(2)
与恶龙缠斗过久,自身亦成为恶龙;凝视深渊过久,深渊将回以凝视…