one_layer.m

代码为本人为写毕业论文

clc;
clear;
%define variable
depth=300;
velocity_up=2800;
velocity_down=3500;
sample_inter=1;
trace_num=49;
trace_length=600;
density=1;
wave_t=160; %define the wavelet length 160ms
fp=45; %the wavelet main frequence is 45hz
%-----------------------
t=2*depth/velocity_up;%自激自收时间
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
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;
%disp(syn_l);
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);