Math 2.0

Numerical analysis Definition

Regularization of ill-posed problems

Builds on Hilbert spaceMoore–Penrose inverseWell-posed problemCompact operatorTikhonov regularization

Definition

Let X,YX, Y be Hilbert spaces, let T ⁣:X→YT \colon X \to Y be a bounded linear operator, and let T†T^\dagger be its Moore–Penrose generalized inverse, defined on 𝒟(T†)=ℛ(T)+ℛ(T)⊥\mathcal{D}(T^\dagger) = \mathcal{R}(T) + \mathcal{R}(T)^\perp. A regularization of T†T^\dagger is a family (Rα)α>0(R_\alpha)_{\alpha > 0} of maps Y→XY \to X together with a parameter choice rule α ⁣:(0,∞)×Y→(0,∞)\alpha \colon (0, \infty) \times Y \to (0, \infty) such that:

  1. each RαR_\alpha is continuous;
  2. for every y∈𝒟(T†)y \in \mathcal{D}(T^\dagger), sup⁡{α(δ,yδ):yδ∈Y, ∥yδ−y∥≤δ}→0\sup \set{ \alpha(\delta, y^\delta) : y^\delta \in Y, \ \norm{y^\delta - y} \le \delta } \to 0 as δ→0\delta \to 0;
  3. for every y∈𝒟(T†)y \in \mathcal{D}(T^\dagger), sup⁡{∥Rα(δ,yδ)yδ−T†y∥:yδ∈Y, ∥yδ−y∥≤δ}→0\sup \set{ \norm{R_{\alpha(\delta, y^\delta)} y^\delta - T^\dagger y} : y^\delta \in Y, \ \norm{y^\delta - y} \le \delta } \to 0 as δ→0\delta \to 0.

The pair is then called a convergent regularization method.

Remarks

Example: Tikhonov regularization Rα=(T∗T+αI)−1T∗R_\alpha = (T^* T + \alpha I)^{-1} T^* with an a priori rule α=α(δ)\alpha = \alpha(\delta) such that α(δ)→0\alpha(\delta) \to 0 and δ2/α(δ)→0\delta^2 / \alpha(\delta) \to 0. Each RαR_\alpha is bounded with ∥Rα∥≤1/(2α)\norm{R_\alpha} \le 1/(2\sqrt{\alpha}), which gives clause 1. The error splits as ∥Rαyδ−T†y∥≤∥Rαy−T†y∥+δ/(2α)\norm{R_\alpha y^\delta - T^\dagger y} \le \norm{R_\alpha y - T^\dagger y} + \delta/(2\sqrt{\alpha}). The first term tends to 00 as α→0\alpha \to 0 by spectral calculus, and the second by the choice of α\alpha, which gives clause 3. Non-example: Rα=T†R_\alpha = T^\dagger for all α\alpha fails clause 1 whenever ℛ(T)\mathcal{R}(T) is not closed, because T†T^\dagger is then unbounded. This happens for every injective compact operator with infinite-dimensional range, the typical model of an ill-posed problem. Other standard methods are truncated singular value decomposition, Landweber iteration with the iteration index as parameter, and the discrepancy principle as a parameter choice rule.

Reading path

To understand this page from the ground up, read these 6 pages in order. The foundations come first.

Show the path as a graph
Each page uses the pages above it. The colour of the bar on a box shows how strong the evidence is, and the box also gives it as text. The list has the same pages.
  • no evidence level
  1. Compact operator concept
  2. Hilbert space concept
  3. Moore–Penrose inverse concept
  4. Tikhonov regularization concept
  5. Well-posed problem concept
  6. Regularization of ill-posed problems definition

Further reading

  1. Web https://en.wikipedia.org/wiki/Regularization_(mathematics)
About this page
TagX2AMH38NCR
URIm2://X2AMH38NCR
Kinddefinition
WritingW0 stub
MSC65J20, 65F22, 47A52
DefinesRegularization of ill-posed problems
UsesHilbert space, Moore–Penrose inverse, Well-posed problem, Compact operator, Tikhonov regularization

Markdown · JSON · Edit on GitHub

What this page needs for W1

  • an inventory row: the row of the inventory that the node accounts for (0 of 1)
  • a source with a locator (0 of 1)

The writing levels say how complete the text of a page is. They do not measure its truth.

What links here

No page links to this page yet.