程序说明

程序说明

一、直射线追踪及LSQR算法

SL文件夹下：

1、可执行文件

（1）shijianshuju.m 合成时间数据T

可修改参数M (横向网格数），N (纵向网格数)，P（横向网格间距），

Q（纵向网格间距），S0（给定模型）

（2）fanyan.m 运用LSQR进行反演

可修改参数M (横向网格数），N (纵向网格数)，P（横向网格间距），

Q（纵向网格间距）

注：当用文件（1）合成的数据时，M，N，P，Q需与文件1中的保持一致。（3）chengtu.m 用得到的数据画出图像

可修改参数M (横向网格数），N (纵向网格数)，P（横向网格间距），

Q（纵向网格间距），range1，range2（要显示的数据范围）注：这里的M，N，P，Q需与文件（2）中的保持一致。

2、函数

（1）HCSJ.m 合成时间数据T

（2）ZuoBiao.m射线与网格垂直线的交点坐标及射线的与水平方向夹角的正弦值

（3）JieDian.m 射线与网格交点上方的节点

（4）XSJZ 求系数矩阵

（5）convert2file.m 保存数据（肖锋老师编译的函数）

二、线性插值射线追踪及SIRT算法

LTI-SIRT文件夹下

1、可执行文件

（1）WSXZZ.m 线性插值射线追踪

可修改参数M（横向网格数），N（纵向网格数），P（横向网格单位），

Q（纵向网格单位），S（慢度模型）

（2）shuju.m 合成时间数据T

可修改参数M（横向网格数），N（纵向网格数），P（横向网格单位），

Q（纵向网格单位），S（慢度模型）

（3）diedai.m SIRT 迭代反演

可修改参数M（横向网格数），N（纵向网格数），P（横向网格单位），

Q（纵向网格单位），iter（迭代次数）

注：当用文件（2）合成的数据时，M，N，P，Q需与文件（2）中的保持一致。

2、函数

（1）Zuobiao.m 各计算点的坐标

（2）XQCL 线性插值射线追踪法中的向前处理过程

（3）LTI.m 线性插值

（4）HCSJ.m 利用向前处理的得到的各计算点合成时间数据T

（5）SXJZ.m 求取系数矩阵

（6）XHCL.m 线性插值射线追踪方法的向后处理过程

（7）ZXZS.m 选出接收点所在单元内个计算点到接收点的最小旅行时

（8）XSJZ.m 系数矩阵中的一行

（9）SIRT.m 联合迭代重建算法

三、最短路径射线追踪法

ZDLJ文件夹下

1、可执行文件

（1）zuiduanlujing.m 节点设置在单元角点上是的最短路径法

可修改参数Length（剖面横向长度），Width（剖面垂直长度），p（横向网格单位），q（纵向网格单位），VDOTMN（给定速度模型）

（2）GZDLJ1.m 在单元边界上设置一个节点时的最短路径法

可修改参数Length（剖面横向长度），Depth（剖面垂直长度），jsd（边界上设置的计算点个数），p（横向网格单位），q（纵向网格单位），SDOTMN

（给定慢度模型）

（3）GZDLJ.m 在单元边界上设置三个节点时的最短路径法

可修改参数Length（剖面横向长度），Depth（剖面垂直长度），jsd（边界上设置的计算点个数），p（横向网格单位），q（纵向网格单位），SDOTMN

（给定慢度模型）

2、函数

（1）gVnew_jisuan.m 节点设置在单元角点上时，与任意计算点相连的计算点的集合

（2）Vnew_jishuan.m 节点设置在单元边界上时（一个节点），与任意计算点相连的计算点的集合

（3）velocity_jisuan_1.m 节点设置在单元边界上时（一个节点）,任意两点间的慢度值

（4）Vnew_jishuan_3.m 节点设置在单元边界上时（三个节点），与任意计算点相连的计算点的集合

（5）velocity_jisuan.m 节点设置在单元边界上时（三个节点）,任意两点间的慢度值

四、论文中图的参数

运用程序：shijianshuju.m，fanyan.m，chengtu.m

1、图3.2

M=8 % (横向网格数）

N=15 % (纵向网格数)

P=2; %（横向网格间距）

Q=2; %（纵向网格间距）

S0=ones(N,M)*(1/0.13); %给定的模型

图3.7-图3.9 运用程序WSXZZ.m

2、图3.7

M=8; %横向网格数

N=15; %纵向网格数

P=1; %横向网格单位

Q=1; %纵向网格单位

S=ones(N,M)*(1/0.13); %慢度模型

3、图3.8

M=8; %横向网格数

N=15; %纵向网格数

P=1; %横向网格单位

Q=1; %纵向网格单位

S=ones(N,M)*(1/0.13); %慢度模型

S(7:8,4:5)=1/0.3;

4、图3.9

M=8; %横向网格数

N=15; %纵向网格数

P=1; %横向网格单位

Q=1; %纵向网格单位

S=ones(N,M)*(1/0.13); %慢度模型

S(7:8,4:5)=1/0.033;

图3.15-3.17 运用程序zuiduanlujing.m

5、图3.15

Length=8;% 剖面横向长度

Width=15; %剖面垂直长度

p=1; %横向网格单位

q=1;%纵向网格单位

VDOTMN=ones(m,n)*0.13; % 给定速度模型6、图3.16

Length=8;% 剖面横向长度

Width=15; %剖面垂直长度

p=1; %横向网格单位

q=1;%纵向网格单位

VDOTMN=ones(m,n)*0.13; % 给定速度模型VDOTMN([7:9],[4:6])=0.3;

7、图3.17

Length=8;% 剖面横向长度

Width=15; %剖面垂直长度

p=1; %横向网格单位

q=1;%纵向网格单位

VDOTMN=ones(m,n)*0.13; % 给定速度模型VDOTMN([7:9],[4:6])=0.033;