Questions tagged [inverse-problems]

Inverse problems involve for example reconstruction of an object based on physical measurements and finding a best model/parameters out of a family given observed data. Typically the corresponding "forward" problems are well-posed and can be solved straightforwardly, while the inverse problems are often ill-posed. Not to be confused with the (inverse) tag.

Typically a forward problem is a family of well-posed problems that are parametrized by some set $\mathcal{P}$. We can often write the forward problems as a nonlinear transformation $G_p$ indexed by $p\in \mathcal{P}$ such that given initial data $d$ the transformation returns some final data (or solution) $s_p$. Or, in notations: $$s_p = G_p[d]$$

The forward problem then is the problem of finding $s_p$ for a given parameter $p$ and a given initial data $d$. Some examples are:

  • Given a family of wave equations parametrised by the potential function $V$ $$ -\partial_t^2 u(x,t) + \partial_x^2 u(x,t) = V(x,t)u(x,t) $$ one possible forward problem would be solving for the function $u$ and its velocity $\partial_t u$ at time $t = T$ given the initial data $(u,\partial_t u)_{t = 0} = (f,g)$.

  • We let the parameter space $\mathcal{P}$ be the space of compactly supported smooth functions on Euclidean space $\mathbb{R}^n$, and for any initial data $d$ we let $G_p[d] = s_p$ be the X ray transform of $p$.

  • We let the parameter space $\mathcal{P}$ be the density distribution in the bedrock of a region of land. We let that transformation $G_p$ be the mapping that sends the initial data of the strength of controlled explosion to the final data which is the observed surface seismic wave. The operator $G_p$ can be (in theory) computed from $p$ based on accepted model of the elastodynamics of the interior of the earth.

In principle the forward problems can be solved, at least numerically, by straightforward methods.

The corresponding inverse problems are the problems of finding the best (or some) parameter $p$ such that given some initial data elicits some observed solution. Quite often the inverse problems are ill-posed: there may not be admissible parameters at all that reproduce the transformation from initial to final data (in this case the observed data are believed to contain errors, or that our a priori assumptions on what the operators are can be fallacious); or there could be multiple admissible parameters that produce the same observed results.

The inverse problem corresponding to the examples above are

  • Solving for the potential $V(x,t)$ between time $t = 0$ and $t = T$ given the solution $u$ and its time derivative at those two boundary times.
  • Finding the compactly support function $p$ given its X ray transform
  • Solving for the density of the rock by seismic sounding.

Some often-studied inverse problems include the inverse scattering problem in partial differential equations and the inverse Sturm-Liouville problem in ordinary differential equations. Typical applications include medical imaging (X rays, CAT scans), seismic sounding, various tomography methods, machine learning, statistical analysis and more.

295 questions
5
votes
2 answers

What is an Inverse problem in Mathematics?

I have come accross a lot of articles that talk about inverse problems. However, I dont really appreciate the uses due to my poor understanding of the notion. From the mathematics point of view, when does a problem qualify to be called an inverse…
OKPALA MMADUABUCHI
  • 1,063
  • 6
  • 16
4
votes
2 answers

Given a vector field $f:\mathbb R^3\to\mathbb R^3$, is there a mass distribution that generates $f$ as its Newtonian gravity?

Let $f:\mathbb R^3\to\mathbb R^3$ be a smooth bounded vector field. I want to produce a density $\rho:\mathbb R^3\to\mathbb R$ such that the Newtonian acceleration experienced by a particle in the space $\mathbb R^3$ is given by $f$. Actually it's…
Derivative
  • 1,566
4
votes
2 answers

What is the meaning of variational regularization methods for inverse problems?

I know if a inverse problem is ill posed, then we can partially obtain the information about the solution by applying regularization techniques. The most commonly used regularization method is Tikhonov regularization. We can also use iterative…
kapil
  • 709
2
votes
1 answer

Confusion related to inverse problems in statistics

I am getting started with inverse problems in statistics. However, I didn't something related to it. I was reading this paper http://math.uni-heidelberg.de/studinfo/reiss/CavalierInvProb.pdf. It says The classical problem is the following : let A…
user34790
  • 4,192
2
votes
0 answers

Inverse conductivity problem - how to obtain the reconstruction numerically?

I'm working on an Electrical Impedance Tomography (or EIT) problem i.e. I wish to reconstruct an image based on information obtained at the boundary. The problem consists of the generalised LaPlace equation with Neumann boundaries…
ANYN11
  • 365
2
votes
1 answer

Regularization and Optimization

What is the difference between regularization and optimization? I keep reading these terms in various papers on solutions of inverse problems but none of them describe what these terms physically mean.
1
vote
0 answers

What is the main difference between interior and exterior inverse problem?

I wonder whether there is any difference between interior inverse problem and exterior inverse problem with respect to stability, numerical solution or anything.. Is it harder to solve interior inverse problem than exterior inverse problem ? What…
izaag
  • 908
1
vote
0 answers

Why do we choose $\lambda,\mu \in C^1(\partial D)$ and $\partial D $ in the class of $C^4$ for this problem?

I am reading article with title " Integral equation methods in Inverse Obstacle scattering with a Generalized Impedance Boundary Condition" written by Rainer Kress. The problem is fomulated in the follwing. Let $D$ be simply connected domain in…
izaag
  • 908
1
vote
1 answer

Is there a reason for why some problems reduce to Caldéron's inverse conductivity equation?

Is there a reason for why some problems reduce to Caldéron's inverse conductivity equation? Or is it "coincindence"? Such problems include, optical tomography, inverse scattering. E.g. in https://ims.nus.edu.sg/events/2018/theo/files/tutn3.pdf
mavavilj
  • 7,270
1
vote
1 answer

Well-posed Inverse Problems Definition

$X,Y$ Banach spaces, $F: X \to Y$. Then the operator-equation "$Fx = y$" is well-posed, if for all $y \in Y$ (1) there exists $x \in X: Fx = y$ (2) the solution is unique (3) the solution continuously depends on y Q: What exactly is meant by (3)?…
Pazu
  • 1,077
1
vote
0 answers

Well posed or ill posed inverse problem?

Consider the problem of finding a control parameter $p(t)$ in the following \begin{equation} \frac{\partial u}{\partial t}= \frac{\partial^2 u}{\partial x^2}+p(t)u(x,t)+f(x,t), \ ~~\ 0\leq x \leq1 ,\ ~~\ 0< t \leq1 \end{equation} with the initial…
0
votes
2 answers

Roughning matrix format for the first-order Tikhonov Regularization (inverse problem)

I have been trying to solve the regularized least square problem of min||Gm-d||^2 + a ||Lm||^2 using first order Tikhonov regularization method. the general form of L for calculating the first derivative of m…
K-1
  • 101
  • 2
0
votes
1 answer

how can we solve this equation without mellin transform?

given the equation (functional equation) $$ f(x)+f(2x)+f(3x)+.... =g(x) $$ we can use the Mobius tranform to obtain $$ f(x)=\sum_{n=1}^{\infty}g(nx)\mu(n) $$ however, what can we do with the equation $$ f(x)-f(2x)+f(3x)-f(4x)+... =g(x) $$ ? how can…
Jose Garcia
  • 8,506