【文件属性】:
文件名称:matlab瑞利波频散曲线代码-Rayleigh-Wave-Monte-Carlo-Inversion:一种联合反演R波频散曲线的代码
文件大小:126KB
文件格式:ZIP
更新时间:2021-06-01 21:42:37
系统开源
matlab瑞利波频散曲线代码瑞利波蒙特卡罗反演
一种联合反演R波频散曲线的代码
MC_R_CyU_VariableThickness.m
此代码执行瑞利波频散曲线(相位和群联合)的蒙特卡罗反演(马尔可夫链
MC)。
它搜索适合观察数据的一维高概率模型。
该代码使用
Herrmann
地震学代码中的
surf96
()。
这组代码应该在尝试运行反演之前安装。
与任何
MC
代码一样,运行它需要时间!
速度取决于您的计算机和
Matlab
版本。
根据迭代次数,反演可能需要长达一天的时间,但经验表明
100000
次迭代就可以解决问题。
要有耐心!
请阅读:Bosch
(1999;
2001
and
2005)
和
Mosegaard
&
Sambridge
(2002)
了解方法和背景。
代码从一个模型开始(要修改的层和要考虑的层数应按照说明进行设置),并且很难找到解决方案。
这是通过选择一层或一组层(随机)并改变层的速度或厚度来完成的。
图层不允许出现或消失。
如果“扰动”减少了错误,则它被接受,否则它有被接受或拒绝的概率。
如果被拒绝,则重复先前的模型。
参数:
IT:迭代次数
M
【文件预览】:
Rayleigh-Wave-Monte-Carlo-Inversion-master
----R_U.dat(485B)
----MC_R_CyU_VariableThickness.m(11KB)
----R_C.dat(485B)
----Example2.png(31KB)
----README.md(3KB)
----Example1.png(78KB)
----bin()
--------make_par_surf96.m(2KB)
--------tempe.m(795B)
--------writedisp_surf96.m(1KB)
--------dispR_surf96.m(847B)
--------disp_obs.dsp(9KB)
--------readdisp_surf96.m(612B)
--------writemod_surf96.m(1KB)
--------viridis.m(12KB)
--------sobs.d(75B)
--------readmod_surf96.m(715B)
----model0.5(949B)
----model0.Crust2(928B)