Analysis of regularized inversion of data corrupted by white Gaussian noise

Abstract

Tikhonov regularization is studied in the case of linear pseudodifferential operator as the forward map and additive white Gaussian noise as the measurement error. The measurement model for an unknown function u(x) ism(x)=Au(x)+δϵ(x),?>where δ > 0 is the noise magnitude. If ε was an L2-function, Tikhonov regularization gives an estimateTα(m)= arg min u∈Hr∥Au−m∥L22+α∥u∥Hr2?>for u where α = α(δ) is the regularization parameter. Here penalization of the Sobolev norm covers the cases of standard Tikhonov regularization (r = 0) and first derivative penalty (r = 1). Realizations of white Gaussian noise are almost never in L2, but do belong to Hs with probability one if s < 0 is small enough. A modification of Tikhonov regularization theory is presented, covering the case of white Gaussian measurement noise. Furthermore, the convergence of regularized reconstructions to the correct solution as δ → 0 is proven in appropriate function spaces using microlocal analysis. The convergence of the related finite-dimensional problems to the infinite-dimensional problem is also analysed.