diffraction_wave.zip_wave_地距_地震勘探_时距曲线_绕射波

clc,clear h = 1000; d = 200; v = 1500; x = -1000:2000; t = abs(x)/v;%直达波时距曲线 td = (sqrt(h^2 + d^2) + sqrt(h^2 + (x-d).^2))/v;%绕射波时距曲线 plot(x,t,'k-.'); hold on plot(x,td,'r') set(gca,'ydir','reverse') xlabel('Short-receiver distance(m)') ylabel('Two-way traveltime(s)') legend('直达波','绕射波')