%共炮點道集
clc;
clear;
%define?variable
depth=100;
velocity_up=2000;
velocity_down=3000;
offset_min=10;
sample_inter=1;
trace_num=50;
trace_length=600;
density=1;
wave_t=160;?%define?the?wavelet?length?160ms
fp=35;?%the?wavelet?main?frequence?is?45hz
%------------------------------------------------
ref=（density*velocity_down-density*velocity_up）/（density*velocity_down+density*velocity_up）;
sample_num=trace_length/sample_inter;
for?i=1:trace_num
for?j=1:sample_num
R（ji）=0;
t（i）=sqrt（（2*depth）^2+（offset_min*i）^2）/velocity_up;
k（i）=int32（t*1000/sample_inter）