common_sourse.m

代码为本人为写毕业论文

%共炮点道集
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(j,i)=0;
end
t=sqrt((2*depth)^2+(offset_min*i)^2)/velocity_up;
k=int32(t*1000/sample_inter);
R(k,i)=ref;
end
%------------------------------------------------------------------------------------
figure(1);
wigb(R);
w_t=-wave_t/(2*1000):0.001:wave_t/(2*1000);
rick=(1-2*(pi*fp*w_t).^2).*exp(-(pi*fp*w_t).^2);
figure(2);
plot(w_t,rick);
for i=1:trace_num
syn(:,i)=conv(R(:,i),rick);
end
wave_l=wave_t/sample_inter;
syn_l=wave_l+sample_num-1;
for i=1:trace_num
for j=(wave_l/2):1:(syn_l-wave_l/2)
syn_n((j+1-wave_l/2),i)=syn(j,i);%整体提升了wave_l/2,使得振幅最大值对应了反射系数值
end
end
figure(3);
wigb(syn_n,1,[10:10:500],[1:600]);
title('共炮点道集','fontsize',20);
xlabel('offset(m)');
ylabel('time(ms)');