图像矩阵matlab代码贝叶斯
MATLAB代码可运行维稳健的MCMC来进行分层贝叶斯反演,如Dunlop,Iglesias和Stuart撰写的《分层贝叶斯能级集反演》所述。
提供了三个示例正向模型:直接点观测,地下水流模型和电阻抗层析成像模型。
提供了以下文件:
run_mcmc.m
,在文件开头定义的参数(网格分辨率,样本数量,先前的平滑度等)下执行MCMC。
可选地,输出可以显示为数字。
长度比例参数和一些傅立叶模式的痕迹保存在内存中。
gaussrnd.m根据给定的平滑度参数alpha,反比例尺参数tau和网格大小N,从本文概述的高斯先验生成2D样本。该样本在傅立叶空间中提供,并重整为N
^
2
*
1向量。
make_lvl.m取一个代表连续函数的方阵,以及两个级别的阈值,返回代表分段常数函数的矩阵。
阈值函数采用的值在此文件中定义。
ell.m选择适当的正向模型,将阈值函数映射到输出测量值。
这三个模型包含在model_id.m,model_gwf.m和model_eit.m中。
model_id.m对分段常数字段进行直接观察。
观察是在J点的正方形网格上进行的。
J在文件
2022-11-10 15:10:08
19KB
系统开源
1