如何用两行在 plot 上应用功能 heatmap 或 imagesc？

Question

我需要像这样2画我的 plot 1 ，但不知道我该怎么做。 当曲线接触红线时，曲线的一部分变为红色[包络区域的总振幅变大] 3 。 Colors（和彩条）应该显示幅度的变化。 据我了解，颜色正在发生变化，因为抛物线中有冲动的积累（例如，Ricker 的冲动）。

clear all
phi=pi/6;
b=4000;
 if (phi>pi/4)|(phi<-pi/4)
    warndlg('Invalid value');
else
a=tan(phi);
x=-4000:100:4000;
y=zeros(1,length(x));
z=a*x+b;
plot(x,z,'r','LineWidth',3);xlim([-5000 9000]);ylim([0 7000]);
set(gca,'YDir','reverse');
xlabel('X,m');
zoom on
title({'Решение задачи сейсмического сноса' , 'для плоской наклонной границы методом огибающей'});
daspect([1 1 1])
zeta=z./cos(phi);
hold on[enter image description here][1]
dz=z*tan(phi)+x;
pause
plot(dz,zeta,'b','LineWidth',3);
pause
for i=1:4:length(x)
xk=x(i)-4000:12.5:x(i)+4000;    
d(i)=z(i)*tan(phi)+x(i);
hyp=sqrt((xk-x(i)).^2+z(i).^2);
p=plot(dz(i),zeta(i),'.k','MarkerSize',25);
plot(x(i),z(i),'.r','MarkerSize',25);
lz=plot([x(i) d(i)],[z(i)  0],'Color',[0.5 0.5 0.5]);
lzeta=plot([d(i) d(i)],[zeta(i) 0],'Color',[0.5 0.5 0.5]);
h=plot(xk,hyp,'k');
pause
delete(p)
delete(lz)
delete(lzeta)
plot(dz(i),zeta(i),'.c','MarkerSize',25);
delete(h)
plot(xk,hyp,'c');
end
end
clear all
close all