clear?all;
clc;
L=71;
fm=40;%采樣頻率
dt=0.001;%時間間隔
R=zeros（1L）;%雷克子波長度
%雷克子波
B=（L-1）/2;
for?m=1:1:L????
??
????t=dt*（m-B）;
????R（m）=（1-2*（t*pi*fm）^2）*exp（-（t*pi*fm）^2）;
end


subplot（241）;
?plot（R）;hold?on;

%?%卷積
%?RE=zeros（1128）;
%?RE（164）=0.5;%反射系數
%?Result=zeros（128+L-1）;%卷積長度
%?Result=conv（RRE）;
%??subplot（242）;
%?plot（Result）;hold?on;
%?title（‘一維卷積‘）

%時深轉換
c=256;
M=load（‘Model_V.txt‘）;%輸入數據
V=（reshape（Mcc））‘;%轉換成256*256的矩陣
dh=10;%深度間隔

t=zeros（2c）;??%時間域分界面矩陣??????
T=zeros（cc）;%時間矩陣

???
????????
for?j=1:1:c??%按列循環
????T（1j）=0;??%第一行時間為零
????d=1;????%深度分界面矩陣行數
????for?i=2:1:c
????????T（ij）=2*dh/V（ij）+T（i-1j）;%時間矩陣
?????????if?V（i-1j）??????????h（dj）=（i-1）*dh;%深度分界面矩陣
??????????d=d+1;
?????????end
?????????t（dj）=ceil（T（ij）/0.001）;
???????????
?????????
????end
????
????????
end
%時間域賦速度值
for?j=1:1:c
????
????n=0;
????for?i=2:1:c
????????for?k=n+1:1:ceil（T（ij）/0.001）
????????????
????????????Vt（kj）=V（i-1j）;
????????????n=ceil（T（ij）/0.001）;
????????
???????????????
????????end
?????end
????
end
I=size（Vt1）;%取行數
for?j=1:1:c
????for?k=1:1:I
?????????if?k>t（1j）?+t（2j）
?????????????Vt（kj）=V（cc）;
?????????end???
????end
end
%反射系數

????????????
for?j=1:1:c??%按列循環
???
????for?i=2:1:I
????????
?????????Re（ij）=（Vt（ij）-Vt（