%計算并畫出單發雙收聲息在砂巖層（厚度h=?5m，Vp?=?3850m/s）
%以及上下圍巖為泥巖（厚度為10m，Vp?=?1800m/s）的
%時差曲線示意圖
%?函數主體?%
function?Time_curve（）
h?=?0.5:0.01:25;%深度
L?=?length（h）;
delta_t?=?zeros（L1）;%初始化?聲波時差
for?i?=?1:1:L
????if?h（i）?<=?10
????????delta_t（i）?=?1000000/1800;
????end
????if?h（i）?>?10?&&?h（i）?????????delta_t（i）?=?2*（10^6）*（?（10.5-h（i））?/?1800?+?（h（i）-10）?/?3850?）;?
????end
????if?