HJ-sampler: A Bayesian sampler for inverse problems of a stochastic process by leveraging Hamilton-Jacobi PDEs and score-based generative models

Tingwei Meng, Zongren Zou, Jérôme Darbon, George Em Karniadakis
{"title":"HJ-sampler: A Bayesian sampler for inverse problems of a stochastic process by leveraging Hamilton-Jacobi PDEs and score-based generative models","authors":"Tingwei Meng, Zongren Zou, Jérôme Darbon, George Em Karniadakis","doi":"arxiv-2409.09614","DOIUrl":null,"url":null,"abstract":"The interplay between stochastic processes and optimal control has been\nextensively explored in the literature. With the recent surge in the use of\ndiffusion models, stochastic processes have increasingly been applied to sample\ngeneration. This paper builds on the log transform, known as the Cole-Hopf\ntransform in Brownian motion contexts, and extends it within a more abstract\nframework that includes a linear operator. Within this framework, we found that\nthe well-known relationship between the Cole-Hopf transform and optimal\ntransport is a particular instance where the linear operator acts as the\ninfinitesimal generator of a stochastic process. We also introduce a novel\nscenario where the linear operator is the adjoint of the generator, linking to\nBayesian inference under specific initial and terminal conditions. Leveraging\nthis theoretical foundation, we develop a new algorithm, named the HJ-sampler,\nfor Bayesian inference for the inverse problem of a stochastic differential\nequation with given terminal observations. The HJ-sampler involves two stages:\n(1) solving the viscous Hamilton-Jacobi partial differential equations, and (2)\nsampling from the associated stochastic optimal control problem. Our proposed\nalgorithm naturally allows for flexibility in selecting the numerical solver\nfor viscous HJ PDEs. We introduce two variants of the solver: the\nRiccati-HJ-sampler, based on the Riccati method, and the SGM-HJ-sampler, which\nutilizes diffusion models. We demonstrate the effectiveness and flexibility of\nthe proposed methods by applying them to solve Bayesian inverse problems\ninvolving various stochastic processes and prior distributions, including\napplications that address model misspecifications and quantifying model\nuncertainty.","PeriodicalId":501215,"journal":{"name":"arXiv - STAT - Computation","volume":null,"pages":null},"PeriodicalIF":0.0000,"publicationDate":"2024-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - STAT - Computation","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2409.09614","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0

Abstract

The interplay between stochastic processes and optimal control has been extensively explored in the literature. With the recent surge in the use of diffusion models, stochastic processes have increasingly been applied to sample generation. This paper builds on the log transform, known as the Cole-Hopf transform in Brownian motion contexts, and extends it within a more abstract framework that includes a linear operator. Within this framework, we found that the well-known relationship between the Cole-Hopf transform and optimal transport is a particular instance where the linear operator acts as the infinitesimal generator of a stochastic process. We also introduce a novel scenario where the linear operator is the adjoint of the generator, linking to Bayesian inference under specific initial and terminal conditions. Leveraging this theoretical foundation, we develop a new algorithm, named the HJ-sampler, for Bayesian inference for the inverse problem of a stochastic differential equation with given terminal observations. The HJ-sampler involves two stages: (1) solving the viscous Hamilton-Jacobi partial differential equations, and (2) sampling from the associated stochastic optimal control problem. Our proposed algorithm naturally allows for flexibility in selecting the numerical solver for viscous HJ PDEs. We introduce two variants of the solver: the Riccati-HJ-sampler, based on the Riccati method, and the SGM-HJ-sampler, which utilizes diffusion models. We demonstrate the effectiveness and flexibility of the proposed methods by applying them to solve Bayesian inverse problems involving various stochastic processes and prior distributions, including applications that address model misspecifications and quantifying model uncertainty.
查看原文
分享 分享
微信好友 朋友圈 QQ好友 复制链接
本刊更多论文
HJ-采样器:利用汉密尔顿-雅可比 PDE 和基于分数的生成模型解决随机过程逆问题的贝叶斯采样器
文献中对随机过程与最优控制之间的相互作用进行了广泛的探讨。随着最近扩散模型使用的激增,随机过程越来越多地被应用于样本生成。本文以对数变换(在布朗运动中称为 Cole-Hopft 变换)为基础,在包含线性算子的更抽象框架内对其进行了扩展。在这一框架内,我们发现科尔-霍普夫变换与最优传输之间的著名关系是线性算子充当随机过程无限小生成器的一个特殊实例。我们还介绍了一种新颖的情况,即线性算子是生成器的邻接,在特定的初始和终端条件下与贝叶斯推理相联系。利用这一理论基础,我们开发了一种新算法,命名为 HJ-取样器,用于给定终端观测值的随机微分方程逆问题的贝叶斯推断。HJ 采样器包括两个阶段:(1)求解粘性汉密尔顿-雅可比偏微分方程;(2)从相关的随机最优控制问题中采样。我们提出的算法自然允许灵活选择粘性 HJ 偏微分方程的数值求解器。我们介绍了求解器的两种变体:基于里卡提方法的里卡提-HJ 采样器和利用扩散模型的 SGM-HJ 采样器。我们将所提出的方法应用于解决涉及各种随机过程和先验分布的贝叶斯逆问题,包括解决模型错误规范和量化模型不确定性的应用,从而证明了这些方法的有效性和灵活性。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 去求助
来源期刊
自引率
0.00%
发文量
0
期刊最新文献
Model-Embedded Gaussian Process Regression for Parameter Estimation in Dynamical System Effects of the entropy source on Monte Carlo simulations A Robust Approach to Gaussian Processes Implementation HJ-sampler: A Bayesian sampler for inverse problems of a stochastic process by leveraging Hamilton-Jacobi PDEs and score-based generative models Reducing Shape-Graph Complexity with Application to Classification of Retinal Blood Vessels and Neurons
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
已复制链接
已复制链接
快去分享给好友吧!
我知道了
×
扫码分享
扫码分享
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1