Choosing α: the L-Curve
You have a powerful cure in hand —
Tikhonov regularization —
but it comes with a dial, the parameter \alpha, and the whole thing only
works if you turn that dial to the right spot. Turn it too little and noise floods
your reconstruction, drowning the answer in wobbles. Turn it too much and you
smear away the real signal, leaving a smooth, featureless blur. There is a Goldilocks amount — and
finding it, without ever knowing the true answer, is the central practical decision in
regularization.
The good news: you are not just guessing. Several clever, semi-automatic methods exist to point you
at a good \alpha. Three are classics:
- The L-curve. Plot the model size \|m_\alpha\| against the data misfit \|Gm_\alpha - d\| on log–log axes, as \alpha sweeps. The graph is L-shaped; the best \alpha sits at the corner.
- The discrepancy principle. If you know the noise level, pick the \alpha that makes the misfit equal to that noise — fit the data exactly as well as the noise allows, no better.
- Generalized cross-validation (GCV). Choose the \alpha that best predicts left-out data, needing no noise estimate.
The trade-off underneath it all
Every one of these methods is really navigating the same tension you met with Tikhonov: the
bias–variance trade-off that runs through all of statistics and machine learning.
- Small \alpha: small bias (you barely bend the answer)
but large variance — the reconstruction jitters wildly as the noise wiggles.
- Large \alpha: small variance (a calm, stable answer)
but large bias — you have smoothed real features out of existence.
The total error is bias-squared plus variance, and it has a minimum somewhere in the
middle. Every selection method below is, at heart, a different way of hunting for that
minimum — the same balancing act as choosing the cutoff in
truncated SVD.
Why the corner is the sweet spot
The L-curve has two arms. The steep vertical arm (small
\alpha) is the under-regularized regime: the misfit barely improves
while the model norm rockets up — you are fitting noise. The flat horizontal arm
(large \alpha) is over-regularized: the model is tiny but the misfit
grows as real signal is thrown away. The corner is the balance point — the
largest drop in misfit for the smallest growth in model size.
Read it like any effort-versus-reward curve. On the vertical arm you are paying a fortune
(a huge, noisy model) for almost no reward (a barely-smaller misfit). On the horizontal arm you are
paying almost nothing but getting nothing back. The corner is the point of diminishing
returns — exactly where a little more fit stops being worth the price in solution size.
Worked example — the discrepancy principle by hand
Suppose you actually know your noise level: each of m = 100 data points
carries independent noise of standard deviation \eta = 0.02. Then a fit
that is "as good as the noise allows" should leave a residual of about
\|Gm_\alpha - d\| \approx \eta\sqrt{m} = 0.02\sqrt{100} = 0.2.
So the recipe is: start with a large \alpha (big misfit) and
decrease it, watching the residual fall. The moment the residual drops to the
noise floor 0.2, stop. Pushing
\alpha any lower would drive the residual below
0.2 — but that means you are now fitting the noise itself, chasing
wobbles that carry no real information. Fit the data exactly as well as the noise allows, and not
one bit better.
Worked example — reading the corner off an L-curve
No noise estimate? Use the L-curve. Compute the Tikhonov solution for a whole ladder of
\alpha values — say \alpha = 10^{-3}, 10^{-2.5}, \dots, 1 —
and for each one record the pair
\big(\log\|Gm_\alpha - d\|,\; \log\|m_\alpha\|\big). Plot those pairs.
Trace the curve from one end: coming down the near-vertical arm, the model norm plunges while the
misfit hardly moves — good progress, keep going. Then the curve bends and flattens
into the horizontal arm, where now the misfit starts climbing fast for little shrinkage in the
model. The elbow between the two — the point of tightest curvature — is the
L-curve's pick. In the interactive figure below, slide \alpha and watch
the marker walk along the curve; park it right at the bend.
Find the corner
The curve is the L-curve for a noisy test problem (axes are log model-norm vs log misfit). Slide
\alpha and the marker travels along it: down the vertical arm for small
\alpha (noise amplified), along the horizontal arm for large
\alpha (over-smoothed). Park it at the corner — that is the
\alpha the L-curve criterion picks.
- L-curve: the corner of log \|m_\alpha\| vs log \|Gm_\alpha - d\|.
- Discrepancy principle: make the misfit match the known noise level.
- GCV: minimise predicted error on left-out data — no noise estimate needed.
It is tempting to treat these methods as black boxes that spit out the answer. They don't.
The truly correct \alpha depends on the noise level — which you often
don't know — and every selection method can fail on some problems:
- The discrepancy principle needs a good noise estimate. Feed it a
wrong noise level and it will confidently over- or under-smooth.
- The L-curve's corner can be ambiguous, rounded, or on some problems
absent altogether — no sharp bend to point at.
- GCV is elegant but occasionally picks a wildly-too-small
\alpha, leaving noise in.
So choosing \alpha stays part science, part art. The single best habit:
look at the reconstruction. Does it show plausible features, or noise speckle, or
an over-smoothed blur? Your eye, informed by what the answer should look like, is a
sanity check no automatic criterion can replace.
Plot solution size against residual size on log–log axes and, for a well-behaved problem, the trace
comes out looking like a capital letter L (or a backwards-L, depending which way
\alpha runs): a steep near-vertical stroke, a sharp turn, then a long
flat foot. The name is just that shape.
But the corner of that L is a genuinely beautiful idea. It is the point of diminishing
returns — the exact place where fitting the data any better only buys you a wildly larger, noisier
solution. That is the same sweet spot you look for on any effort-versus-reward curve:
training time versus accuracy, engineering tolerance versus cost, study hours versus exam marks.
Find the corner, and you have found where extra effort stops paying for itself.