合成地震记录的matlab程序

博客石油www. blogoil .com%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)/deltx)*deltx+deltx; else x2(i)=fix(x2(i)/deltx)*deltx; end

end

%初始化、转换反射系数矩阵 r1=zeros(m,n);

for i=1:n

r1(fix(x1(i)*10),i)=r1temp; end

r2=zeros(m,n);

for i=1:n

r2(fix(x2(i)*10),i)=r2temp; end

%打开子波文件，去掉文件头