%CONE dessine des cones de diffraction dont l'enveloppe est une reflexion
% voir : http://jmmeost.free.fr/Sismique/isis.html
% fig1 = http://jmmeost.free.fr/Sismique/cone.gif  
% fig3 = http://jmmeost.free.fr/Sismique/pente.gif
% fig4 = http://jmmeost.free.fr/Sismique/syncli.gif
% fig5 = http://jmmeost.free.fr/Sismique/copt.gif

clear all

V0 = 2000;

XM = 2000; xx = 0:50:XM; ZM = 2000; zz = 0:50:ZM; 
[X,Z] = meshgrid(xx,zz);

figure(1), clf, hold on

% fig 1 déport nul point diffractant (enlever % sur 2 lignes suivantes)
xd = XM/2; t0 = 0; TM = ZM/V0;
for zd = xd
plot3(xd,zd,0,'ko','MarkerFaceColor','k')
plot3([xd,xd],[0,zd],[0, 0],'k',[0,xd],[zd,zd],[0, 0],'k',[xd,xd],[zd,zd],[0, t0],'k',[xd,xd],[0,0],[0, zd/V0],'k',[0,xd],[0,0],zd/V0*[1,1],'k')
title(['CONE DE DIFFRACTION  XD = ' num2str(xd) ' M, ZD = ' num2str(zd) ' M, V = ' num2str(V0) ' M/S, DT/DX = SIN{\Theta}/V (couleur)']);

% fig 2 reflecteur horizontal (enlever % sur 2 lignes suivantes)
% xd = XM/2; zd = xd; t0 = 0; TM = ZM/V0;
% title(['POINTS DIFFRACTANTS HORIZONTAUX, ZD = ' num2str(zd) ' M, V = ' num2str(V0) ' M/S, DT/DX = SIN{\Theta}/V (couleur)']);
% for xd = xd*.5:50:xd*1.5 
% plot3(xd,zd,0,'ko','MarkerFaceColor','k')

% fig 3 reflecteur pente (enlever % sur 2 lignes suivantes)
% al = 30*pi/180; xd = XM/2; t0 = 0; TM = ZM/V0;
% title(['POINTS DIFFRACTANTS PENTES {\alpha} = ' num2str(al*180/pi) ' °, V = ' num2str(V0) ' M/S, DT/DX = SIN{\Theta}/V (couleur)']);
% for xd = xd*.5:50:xd*1.5   zd = xd*tan(al); 
% plot3(xd,zd,0,'ko','MarkerFaceColor','k')

% fig 4 synclinal (enlever % sur 2 lignes suivantes)
% t0 = 0; TM = ZM/V0;
% title(['POINTS DIFFRACTANTS SUR SYNCLINAL V = ' num2str(V0) ' M/S, DT/DX = SIN{\Theta}/V (couleur)']);
% for xd = 0:50:XM  zd=1000+500*sin(pi*xd/2000); 
% plot3(xd,zd,0,'ko','MarkerFaceColor','k')

% fig 5 point de tir  point diffractant (enlever % sur 2 lignes suivantes)
% xd = XM/2; 
% for zd = xd
% t0 = sqrt(xd^2+zd^2)/V0; TM = 2*ZM/V0;
% plot3(xd,zd,0,'ko','MarkerFaceColor','k')
% plot3([xd,xd],[0,zd],[0, 0],'k',[0,xd],[zd,zd],[0, 0],'k',[xd,xd],[zd,zd],[0, t0],'k',[xd,xd],[0,0],[0, t0+zd/V0],'k',[0,xd],[0,0],t0+zd/V0*[1,1],'k')
% plot3([0,xd],[0,zd],[0, t0],'r','Linewidth',2)
% title(['CONE DE DIFFRACTION XS = ZS = 0,  XD = ' num2str(xd) ' M, ZD = ' num2str(zd) ' M, V = ' num2str(V0) ' M/S, DT/DX = SIN{\Theta}/V (couleur)']);

% dessine les cones T(X,Z) en coloriant avec DTDX
R = sqrt((X-xd).^2 +(Z-zd).^2)+eps; 
T = R/V0+t0; DTDX = ((X-xd))./(R*V0); 
surf(X,Z,T,DTDX)

end
view(30,20)
axis([0 XM 0 ZM 0 TM]), colormap('hsv'), box on, grid on
xlabel('X(M)'); ylabel('Z(M)'); zlabel('T(S)'); 
set(gca, 'ZDir', 'reverse'), colorbar