By formulating the inverse problem of partial differential equations (PDEs) as a statistical inference problem, the Bayesian approach provides a general framework for quantifying uncertainties. In the inverse problem of PDEs, parameters are defined on an infinite-dimensional function space, and the PDEs induce a computationally intensive likelihood function. Additionally, sparse data tends to lead to a multi-modal posterior. These features make it difficult to apply existing sequential Monte Carlo (SMC) algorithms. To overcome these difficulties, we propose new conditions for the likelihood functions, construct a Gaussian mixture based preconditioned Crank-Nicolson transition kernel, and demonstrate the universal approximation property of the infinite-dimensional Gaussian mixture probability measure. By combining these three novel tools, we propose a new SMC algorithm, named SMC-GM. For this new algorithm, we obtain a convergence theorem that allows Gaussian priors, illustrating that the sequential particle filter actually reproduces the true posterior distribution. Furthermore, the proposed new algorithm is rigorously defined on the infinite-dimensional function space, naturally exhibiting the discretization-invariant property. Numerical experiments demonstrate that the new approach has a strong ability to probe the multi-modality of the posterior, significantly reduces the computational burden, and numerically exhibits the discretization-invariant property (important for large-scale problems).
翻译:通过将偏微分方程反问题表述为统计推断问题,贝叶斯方法为量化不确定性提供了通用框架。在偏微分方程反问题中,参数定义在无限维函数空间上,且偏微分方程导出了计算密集的似然函数。此外,稀疏数据往往导致多峰后验分布。这些特征使得现有序贯蒙特卡洛算法难以直接应用。为克服这些困难,我们提出了似然函数的新条件,构建了基于高斯混合分布预处理Crank-Nicolson转移核,并证明了无限维高斯混合概率测度的通用逼近性质。通过整合这三种新工具,我们提出了一种名为SMC-GM的新序贯蒙特卡洛算法。针对该新算法,我们获得了适用于高斯先验的收敛定理,证明序贯粒子滤波器能够准确再现真实后验分布。此外,所提算法在无限维函数空间上严格定义,天然具备离散化不变性。数值实验表明:新方法具备强大的后验多峰探测能力,显著降低了计算负担,并在数值上展现出离散化不变性(这对大规模问题至关重要)。