Statistics for Inverse Problems 101
June 5, 2011 Introduction
I Inverse problem: to recover an unknown object from indirect noisy observations Example: Deblurring of a signal f R 1 Forward problem: µ(x) = 0 A(x − t)f (t) dt Noisy data: yi = µ(xi ) + εi , i = 1,..., n Inverse problem: recover f (t) for t ∈ [0, 1] General Model
I A = (Ai ), Ai : H → IR, f ∈ H
I yi = Ai [f ] + εi , i = 1,..., n
I The recovery of f from A[f ] is ill-posed
I Errors εi modeled as random
I Estimate bf is an H-valued random variable
Basic Question (UQ) How good is the estimate bf ?
I What is the distribution of bf ? Summarize characteristics of the distribution of bf E.g., how ‘likely’ is it that bf will be ‘far’ from f ? Review: E, Var & Cov
(1) Expected value of random variable X with pdf f (x): Z E(X) = µx = x f (x) dx