马尔可夫链蒙特卡洛算法中提议密度的归一化
数值分析
2015-03-31 v2
摘要
我们在根据边界点温度 重构二维热方程中电导率项 的背景下,探讨了在马尔可夫链蒙特卡洛 (MCMC) 算法中对提议密度进行归一化的影响。我们通过实现 Metropolis-Hastings 马尔可夫链蒙特卡洛算法来解决这一非线性反问题。马尔可夫链生成了以观测数据为条件的可能解的概率分布。我们生成候选解 并求解正问题,得到 。在第 步,以概率 ,我们设定 。我们指出了该构造中存在的一些问题,这些问题源于数据项的大幅值和波动性。利用此框架,我们开发了归一化项 和参数 ,以保留我们所能利用的固有稀疏信息。我们考察了该 MCMC 算法变体在重构多个二维电导率函数时的结果。
引用
@article{arxiv.1410.2972,
title = {Normalizations of the Proposal Density in Markov Chain Monte Carlo Algorithms},
author = {Antoine E. Zambelli},
journal= {arXiv preprint arXiv:1410.2972},
year = {2015}
}
备注
6 pages, 6 figures