The discrepancy principle and parameter choice
Choose a regularization level using an explicit noise bound and distinguish a residual target from an object-error guarantee.
On this page
The question for today
What you’ll understand
- State sufficient a priori parameter limits for deterministic convergence.
- Solve a scalar Morozov discrepancy equation with a safety factor greater than one.
- Explain the consequences of an inaccurate noise bound or unattainable target.
A regularized estimator is a family, so selecting its parameter is part of the reconstruction problem. Too much regularization hides signal; too little lets weak-direction noise dominate. We use the zero-centered identity Tikhonov family for a compact linear Hilbert-space operator and data with ||yδ−Ax†||≤δ for δ>0. The discrepancy principle compares a residual with that known norm bound. It does not replace the need to state what counts as noise or what object is being recovered.
A priori rules balance two limits
Without a source-rate assumption, a sufficient deterministic convergence rule is α(δ)→0 together with δ/√α(δ)→0 as δ→0. The first limit removes bias for each fixed minimum-norm exact solution; the second controls worst-case noise amplification. Neither limit alone is sufficient to establish the usual convergence conclusion. Keeping α fixed preserves stability but may retain bias, whereas decreasing it too quickly can amplify a shrinking data error into a nonshrinking object error.
For example, α=δ satisfies both conditions, while α=δ² makes the displayed noise bound remain of constant size. That failure of a sufficient condition does not assert divergence for every special datum, but it prevents this general guarantee. An optimal rate-based choice also needs a source exponent and source bound, often unknown in practice. This motivates a posteriori rules that use the observed residual, though their guarantees still depend on accurate noise information and the properties of the chosen regularization family.
The residual is a monotone spectral function
For Tikhonov, each observed singular component of the residual is multiplied by α/(σ²+α). Its magnitude increases with α. Components orthogonal to the range closure are unchanged because no object can fit them. Consequently the residual norm is continuous and nondecreasing, with lower limit the distance from yδ to the range closure and upper limit ||yδ|| as α grows. If A*yδ is nonzero, the nontrivial component makes the growth strict. A finite target must lie between those limits to admit a positive solution.
Under a valid noise bound and exact data in the range, the irreducible residual is at most δ, so a target τδ with τ>1 lies above it. If ||yδ|| is already at most τδ, the zero object already meets the stopping tolerance and one should state that convention rather than force a positive-parameter equality with no solution. Otherwise an interior target determines a positive parameter. The residual equation can be solved by monotone search, but that numerical root finding is separate from the inverse problem's stability analysis.
A residual tolerance is conditional information
A discrepancy rule can deliver convergent regularization under the appropriate linear assumptions, but its object-error rate is not universal. Rate claims require source conditions and attention to the method's qualification; an a priori rate formula cannot simply be copied onto every discrepancy rule. A residual of size δ does not imply object error of size δ when the inverse is ill conditioned. It means the fit is commensurate with the stated data uncertainty, which is a different and useful conclusion.
Underestimating δ can force the algorithm to reproduce noise or model discrepancy; overestimating it can stop reconstruction before useful information has been extracted. Correlated statistical errors call for a justified weighted norm or a probabilistic tolerance, not substitution of an unexplained standard deviation into a deterministic theorem. L-curves, generalized cross-validation, and other heuristics can be useful choices, but they rely on different assumptions and do not acquire the discrepancy principle's guarantees merely by producing a visually appealing parameter.
Work it through
Solve a discrepancy target exactly
For scalar A=1 and observed yδ=2, use δ=0.1 and τ=2. Find the identity-Tikhonov parameter satisfying the discrepancy equation and the resulting estimate.
The estimator is xαδ=2/(1+α). Its residual magnitude is 2−2/(1+α)=2α/(1+α), increasing from zero toward two.
The target is τδ=0.2, inside that residual range. Solve 2α/(1+α)=0.2, giving 1.8α=0.2 and α=1/9.
The estimate is 2/(10/9)=1.8, and direct subtraction gives residual 0.2. This matches the chosen safety tolerance, not an assertion that the unknown truth is exactly 1.8.
The parameter follows from the estimator's residual curve and the declared uncertainty. The target does not itself identify the true object.
The living notebook
Explore the model
Interactive lab · root-finding
Inspect a certified scalar search
Reset, set target c=2 and initial upper bound U=3. Predict a bracketed root √2; one Step performs one bisection and gives bracket [0,1.5] with midpoint 0.75 and bound 0.75. Compare a fresh Reset, c=4 and U=2, where the root is the endpoint and the bracket collapses immediately. This instrument bisects x²−c, not the discrepancy residual. It illustrates how a scalar parameter search needs a valid bracket and an explicit error bound.
- Midpoint after 0 refinements
- 1.5
- Bracket error bound
- 1.5
- Computed residual · x² − c
- 0.25
Bisection of x² − c on [0,U], capped at 50 iterations for floating-point resolution. The error bound is the maximum distance from the displayed midpoint to either bracket endpoint (half the width in exact arithmetic). Compensated square evaluation preserves the bracket when ordinary multiplication rounds a residual to zero. A truly exact endpoint or midpoint collapses the bracket. This is one monotone quadratic, not a general root solver. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forSolving the parameter equation accurately is distinct from validating the noise model that defines it.
Pause & reason
Check your understanding
Your turn
Try it yourself
Keep the scalar datum 2 and τ=2, but incorrectly declare δ=0.001 when the actual error may be 0.1. Compute the resulting parameter and explain the limitation.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
The new target is 0.002; solve 2α/(1+α)=0.002.
Compare with the solution
The equation gives α=0.002/1.998=1/999 and estimate x=1.998. The much smaller parameter fits the observation much more closely. Because the assumed bound 0.001 does not cover an actual error as large as 0.1, a theorem using that bound no longer certifies the reconstruction. This example does not prove that every underestimated bound produces a bad answer; it shows precisely which premise has failed and why the reduced residual alone does not repair it.
Keep the idea
Take it with you
- Convergence needs a parameter that removes bias without uncontrolled noise gain.
- Morozov targets must be attainable and tied to a valid uncertainty bound.
- Residual agreement is not a universal object-error certificate.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.