图像矩阵matlab代码-bayes-hier:用于贝叶斯分层反演的MATLABMCMC

图像矩阵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在文件