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.

306 questions
35
votes
3 answers

Why are additional constraint and penalty term equivalent in ridge regression?

Tikhonov regularization (or ridge regression) adds a constraint that $\|\beta\|^2$, the $L^2$-norm of the parameter vector, is not greater than a given value (say $c$). Equivalently, it may solve an unconstrained minimization of the least-squares…
25
votes
3 answers

Can we recover a space from its continuous functions?

Let $X$ be a topological space and let $\mathcal{F}(X,\Bbb R)$ be the set of continuous function from $X$ to $\Bbb R$. Can we recover the topology of $X$ by only the knowledge of $\mathcal{F}(X,\Bbb R)$? That is, can we determinate whether or not…
18
votes
1 answer

Can you hear the pins fall from bowling game scores?

Let $\mathbb T=\{1,\dotsc,10\}$ represent the ten pins in a standard game of bowling. Given two sets of pins $T\subseteq S\subseteq \mathbb T$, let's write $p_{S\to T}$ to represent the conditional probability that given the current pins up are $S$,…
15
votes
0 answers

Recovering a binary function on a lattice by studying its sum along closed paths

I have a binary function $f:\mathbb N^2\rightarrow\{0,1\}$. While I do not known $f$ explicitly, I have a "device" located at the origin $(1,1)$ which can do the following: Given an even number $m$, the device runs over all closed walks of length…
14
votes
2 answers

Is there a Lagrangian that produces these equations? How can I find one if it exists?

Consider the two differential equations \begin{align*} \ddot{x}_{A} - \gamma(x_{A} + x_{B}) &= 0, \\ \ddot{x}_{B} + \gamma(x_{A} + x_{B}) &= 0 \end{align*} where $\gamma$ is a constant. I am looking either for a Lagrangian $L = L(x, \dot{x}, t)$…
12
votes
7 answers

A good book on inverse problems for engineers

I'm looking for a book on inverse problems which is suitable for engineers; both introduction and practical applications are required. Currently I'm looking to Parameter Estimation and Inverse Problems by Richard C. Aster, Brian Borchers and…
10
votes
0 answers

A question about the article 'You can't hear the shape of a drum'

I have read the article http://www.math.upenn.edu/~kazdan/425S11/Drum-Gordon-Webb.pdf, , where Gordon and Webb describe in a simple a way the contruction of a pair of isospectral but non isometric plane domains. At pag $52$ they mention the group…
10
votes
0 answers

Inverse Galois problem

The inverse Galois problem conjectures that every finite group is (isomorphic to) the Galois group of some Galois extension of $\mathbb Q$, however it is not known. My question is: what is the smallest finite group such that it is not known…
6
votes
1 answer

When is $R \, A^{-1} \, R^t$ invertible?

In the context of a Gaussian model, I came across a matrix product $R \, A^{-1} \, R^t$ where $R$ is a $m \times n$ rectangular matrix and as implied $A$ is $n \times n$ and invertible. On which properties of $R$ does the existence of $(R \, A^{-1}…
user69958
6
votes
0 answers

Inverse vs. adjoint operators

I hope this is enough of a question (and less of a start of a discussion) to be allowed here. I am trying to get my head around the ideas of inverse and adjoint operators. To keep it simple, let's stay with real, linear ones. Let's say we have $$ y…
bers
  • 562
6
votes
0 answers

How to infer properties of a function from an integral transform?

I have a function $$ f(x) = \int_{0}^{\infty} \exp\left\{-\iint_{R(x,y)} g(x',y') \, dx' \, dy'\right\} \, dy, $$ where $R(x,y) \subset \mathbb{R}^2$ is a given family of closed regions in the $(x', y')$-plane, parametrized by $(x, y)$. Suppose $f$…
6
votes
1 answer

Inverse propagation of information from the PDF of $Y=f(X)$ to the PDF of $X$

Assume a non-linear relation between the random variables $\mathbf{Y} = f(\mathbf{X})$, where $\mathbf{Y}\sim p_Y$ takes values $\mathbf{y} \in \mathbb{R}^M$ and $\mathbf{X}\sim p_X$ takes values $\mathbf{x} \in \mathbb{R}^N$, with $M\leq N$. My…
6
votes
0 answers

How is the Jacobian matrix computed in finite difference problems?

I have come across many papers which reference the Jacobian when solving certain finite difference inverse problems. And I have seen many articles and textbooks which discuss the mathematical properties of the Jacobian in an abstract sense. I have…
6
votes
0 answers

Inverse spectral problem: How to recover the function $ q(x) $?

The forward problem is a second order Sturm-Liouville operator $$ - \frac{d^{2}}{dx^{2}}y(x)+q(x)y(x)=zy(x) $$ with the boundary conditions $ y(0)=0=y(\infty) $. If I know the spectral measure function $ \sigma (x) =\sum_{\lambda_{n} \le x}1 $, then…
6
votes
1 answer

Please, help to identify this numerical constant

I'm trying to find an answer to this question. Let $K(k)$ be the elliptic integral of the first kind and $K'=K(\sqrt{1-k^2})$. According to Abel's theorem (see this link) we know that if $\frac{K'}{K}=\frac{a+b\sqrt{n}}{c+d\sqrt{n}}$ where…
1
2 3
20 21