seisrectime[1] 地震模型正演程序

地震模型正演程序，特别适合物探专业学生学习。

function varargout =seisrectime[1]（varargin）

%z denotes shot number
%x denotes time t  （x1的时间大，看作多次波）
z=-3:0.025:3;
deltx=0.1;
x=0:deltx:30;
v1=1;
v2=1.4;
v3=1.8;
m=length（x）;
n=length（z）;
%反射系数
%r1temp=（v2-v1）/（v2+v1）;
%r2temp=（（v3-v2）/（v3+v2））*（1-r1temp^2）;
r1temp=0.15;
r2temp=0.3;
%时距曲线方程（双曲线）
x1=abs（sqrt（28.74*z.^2+297.29））;
x2=abs（sqrt（33.13*z.^2+132.89））;
%把时间分配到网格上，与deltx有关
for i=1:n
    if （x1（i）/deltx-fix（x1（i）/deltx））>=0.5
        x1（i）=fix（x1（i）/deltx）*deltx+deltx;
    else x1（i）=fix（x1（i）/deltx）*deltx;
    end
end

for i=1:n
    if （x2（i）/deltx-fix（x2（i）/deltx））>=0.5
        x2（i）=fix（x2（i）/

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1357  2007-03-20 14:27  seisrectime[1].m

----------- ---------  ---------- -----  ----

                 1357                    1