Linear methods for non-linear inverse problems
Abstract
We consider the recovery of an unknown function f from a noisy observation of the solution uf to a partial differential equation that can be written in the form L uf=c(f,uf), for a differential operator L that is rich enough to recover f from L uf. Examples include the time-independent Schr\"odinger equation uf = 2uff, the heat equation with absorption term (∂t -x/2) uf=fuf, and the Darcy problem ∇· (f ∇ uf) = h. We transform this problem into the linear inverse problem of recovering L uf under the Dirichlet boundary condition, and show that Bayesian methods with priors placed either on uf or L uf for this problem yield optimal recovery rates not only for uf, but also for f. We also derive frequentist coverage guarantees for the corresponding Bayesian credible sets. Adaptive priors are shown to yield adaptive contraction rates for f, thus eliminating the need to know the smoothness of this function. The results are illustrated by numerical experiments on synthetic data sets.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.