Math 2.0

Integral equations Inverse problems for integral equations Theorem

Picard's criterion for first-kind equations

Builds on Hilbert spaceCompact operatorSingular values of a compact operatorHilbert space adjoint

Statement

Let XX and YY be Hilbert spaces and let K ⁣:X→YK \colon X \to Y be a compact linear operator with adjoint K∗K^*. Let (σn,vn,un)n∈N(\sigma_n, v_n, u_n)_{n \in N} be a singular system of KK, with NN finite or N={1,2,… }N = \set{1, 2, \dots}. This means that the σn>0\sigma_n > 0 are the nonzero singular values of KK, and Kvn=σnunK v_n = \sigma_n u_n and K∗un=σnvnK^* u_n = \sigma_n v_n. It also means that (vn)(v_n) is an orthonormal basis of ker⁡(K)⊥\ker(K)^\perp and (un)(u_n) is an orthonormal basis of the closure of K(X)K(X). Let y∈Yy \in Y. Then the equation Kx=yK x = y has a solution x∈Xx \in X if and only if

  1. y∈ker⁡(K∗)⊥y \in \ker(K^*)^\perp, and
  2. ∑n∈N∣⟨y,un⟩∣2σn2<∞\displaystyle \sum_{n \in N} \frac{\abs{\ip{y}{u_n}}^2}{\sigma_n^2} < \infty.

In this case x†=∑n∈Nσn−1⟨y,un⟩vnx^\dagger = \sum_{n \in N} \sigma_n^{-1} \ip{y}{u_n} v_n is a solution, and it is the unique solution of smallest norm.

Remarks

Proof idea: if Kx=yK x = y, then y⊥ker⁡(K∗)y \perp \ker(K^*) and ⟨y,un⟩=⟨x,K∗un⟩=σn⟨x,vn⟩\ip{y}{u_n} = \ip{x}{K^* u_n} = \sigma_n \ip{x}{v_n}, so the sum in clause 2 equals ∑n∣⟨x,vn⟩∣2≤∥x∥2\sum_n \abs{\ip{x}{v_n}}^2 \le \norm{x}^2 by Bessel's inequality. Conversely, clause 2 makes the series for x†x^\dagger converge. Then Kx†=∑n⟨y,un⟩unK x^\dagger = \sum_n \ip{y}{u_n} u_n is the orthogonal projection of yy onto the closure of K(X)K(X), which is ker⁡(K∗)⊥\ker(K^*)^\perp. Clause 1 makes it equal to yy. The solution x†x^\dagger lies in ker⁡(K)⊥\ker(K)^\perp, so it has the smallest norm among all solutions.

Example: the Volterra operator (Vx)(t)=∫0tx(s) ds(V x)(t) = \int_0^t x(s) \, ds on L2(0,1)L^2(0, 1) is injective with dense range. Its singular system is σn=2(2n−1)π\sigma_n = \frac{2}{(2n - 1)\pi}, vn(s)=2cos⁡((n−12)πs)v_n(s) = \sqrt{2} \cos\big((n - \tfrac{1}{2}) \pi s\big) and un(t)=2sin⁡((n−12)πt)u_n(t) = \sqrt{2} \sin\big((n - \tfrac{1}{2}) \pi t\big). So Vx=yV x = y is solvable if and only if ∑n(2n−1)2∣⟨y,un⟩∣2<∞\sum_n (2n - 1)^2 \abs{\ip{y}{u_n}}^2 < \infty, a condition on the decay of the coefficients of yy that expresses the differentiability of yy.

Going deeper

For a student

The criterion is the infinite-dimensional form of solving a linear system by the singular value decomposition, see singular value decomposition. It shows the ill-posedness of first-kind equations in a quantitative way: an error of size δ\delta in the coefficient ⟨y,un⟩\ip{y}{u_n} becomes an error δ/σn\delta / \sigma_n in the solution. The faster the singular values decay, the smoother the data must be, and the more ill posed the problem. Regularization methods such as truncated singular value expansion and Tikhonov regularization replace the factors 1/σn1/\sigma_n by bounded filter factors. The map y↦x†y \mapsto x^\dagger is the Moore–Penrose inverse of KK on its natural domain.

Reading path

To understand this page from the ground up, read these 5 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. Hilbert space adjoint concept
  4. Singular values of a compact operator concept
  5. Picard's criterion for first-kind equations statement

Further reading

  1. Web https://en.wikipedia.org/wiki/Singular_value_decomposition
About this page
TagMPB3QF09E5
URIm2://MPB3QF09E5
Kindstatement, theorem
WritingW0 stub
MSC45Q05, 47A52, 45B05
StatusThis copy of the site has no build data. Run m2 status MPB3QF09E5.
AboutSingular values of a compact operator
UsesHilbert space, Compact operator, Singular values of a compact operator, Hilbert space adjoint

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.