matlab瑞利波频散曲线代码-Rayleigh-Wave-Monte-Carlo-Inversion:一种联合反演R波频散曲线的代码

时间:2024-06-21 15:29:17
【文件属性】:

文件名称:matlab瑞利波频散曲线代码-Rayleigh-Wave-Monte-Carlo-Inversion:一种联合反演R波频散曲线的代码

文件大小:126KB

文件格式:ZIP

更新时间:2024-06-21 15:29:17

系统开源

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)

网友评论