3.5.37Guidance, Navigation & Control (GNC)

H∞ control — robust to uncertainty (intro)

2,043 words9 min readdifficulty · medium

WHY — the problem we are solving

WHAT goes wrong with classical/LQR design? LQR minimises a quadratic cost for the known model. It has great margins for that model, but says nothing about a plant that differs. H∞ turns "how different can the plant be, and still be OK?" into a solvable optimisation.

HOW we measure "worst case": with a norm on signals and systems.


Building block 1 — signal size (the 2-norm)


Building block 2 — the H∞ norm (system "gain")

Derivation: why the peak of the frequency response = worst-case energy gain

Start from what we actually want to bound — the energy gain: γworst=supw0z2w2=supw0Gw2w2.\gamma_{\text{worst}} = \sup_{w\neq 0}\frac{\|z\|_2}{\|w\|_2}=\sup_{w\neq 0}\frac{\|Gw\|_2}{\|w\|_2}.

Why this step? This is literally "worst output energy per unit input energy" — the physical meaning of robustness margin.

By Parseval's theorem, energy is preserved between time and frequency: w22=12πW(jω)2dω.\|w\|_2^2 = \frac{1}{2\pi}\int_{-\infty}^{\infty}|W(j\omega)|^2\,d\omega.

Why this step? It lets us work frequency-by-frequency, where GG acts as simple multiplication Z(jω)=G(jω)W(jω)Z(j\omega)=G(j\omega)W(j\omega).

Then (scalar case for clarity): z22=12πG(jω)2W(jω)2dω(supωG(jω)2)12πW2dω.\|z\|_2^2=\frac{1}{2\pi}\int|G(j\omega)|^2|W(j\omega)|^2 d\omega \le \Big(\sup_\omega |G(j\omega)|^2\Big)\frac{1}{2\pi}\int|W|^2 d\omega.

Why this step? Pull the biggest gain out of the integral — it upper-bounds every frequency's contribution.

So z22w22supωG(jω)2\dfrac{\|z\|_2^2}{\|w\|_2^2}\le \sup_\omega|G(j\omega)|^2, and by concentrating ww's energy near the peak frequency the bound is achieved. Hence  supwz2w2=G \boxed{\ \sup_w \frac{\|z\|_2}{\|w\|_2}=\|G\|_\infty\ }

That is why the H∞ norm is the worst-case energy amplification. Bounding it bounds the worst case.


The design goal

The Small-Gain Theorem — why this gives robustness

So designing KK to push M\|M\|_\infty below 1 (and, with weights, below γ\gamma) guarantees stability for every uncertainty in the ball. That is the whole point.


Worked example 1 — reading an H∞ norm from a Bode plot

Take G(s)=10s2+0.4s+4G(s)=\dfrac{10}{s^2+0.4s+4} (a lightly damped mode, ωn=2\omega_n=2, ζ=0.1\zeta=0.1).

  • Step: Find peak of G(jω)|G(j\omega)|. Why? G\|G\|_\infty is that peak.
  • Resonant peak magnitude 10/42ζ1ζ2=2.52(0.1)12.6\approx \dfrac{10/4}{2\zeta\sqrt{1-\zeta^2}}=\dfrac{2.5}{2(0.1)}\approx 12.6.
  • So G12.6\|G\|_\infty\approx 12.6: a resonant gust at ω2\omega\approx2 rad/s is amplified ~12.6×.
  • Interpretation: If this is a bending mode of a rocket, H∞ synthesis will insert a weight there to force the loop to attenuate it.

Worked example 2 — small-gain sizing

Suppose the true actuator gain is 1±30%1\pm30\%: Δ\Delta with Δ0.3\|\Delta\|_\infty\le 0.3.

  • Step: Normalise: write uncertainty as 0.3Δ~0.3\,\tilde\Delta, Δ~1\|\tilde\Delta\|_\infty\le 1. Why? Small-gain wants a unit ball.
  • Step: Robust stability requires 0.3M<1M<3.33\|0.3\,M\|_\infty<1\Rightarrow \|M\|_\infty<3.33.
  • So as long as the nominal loop's peak gain at the uncertain channel stays below 3.333.33, any actuator error in ±30% is tolerated.
  • Why this is powerful: one norm inequality certifies an infinite family of plants.

Worked example 3 — weight shaping intuition

We want good tracking (low error at low frequency) and gentle control (roll off high frequency).

  • Choose WS(s)=s/M+ωbs+ωbAW_S(s)=\frac{s/M+\omega_b}{s+\omega_b A} on the sensitivity SS (error/reference). Why? Large WSW_S at low ω\omega forces S|S| small there ⇒ good tracking.
  • Require WSS<1\|W_S S\|_\infty<1. Why? Since WSS<1S(jω)<1/WS(jω)\|W_S S\|_\infty<1\Rightarrow |S(j\omega)|<1/|W_S(j\omega)| at every frequency — you literally draw the ceiling 1/WS1/|W_S| that S|S| must stay under.


Recall Explain it to a 12-year-old (Feynman)

Imagine you build a robot to balance a broom. You practised with one broom, but tomorrow the broom might be heavier, or a fan might blow on it. A normal robot only knows its practice broom. The H∞ robot asks: "What is the meanest push the wind could give, and can I still stay standing?" It designs itself so that even the worst push only tips the broom a little. The "H∞ number" is just the biggest wobble the worst wind can cause — and we make that number small.


Active recall

What does the H∞ norm G\|G\|_\infty represent physically?
The worst-case energy gain supwz2/w2\sup_w \|z\|_2/\|w\|_2, equal to the peak over frequency of the largest singular value supωσˉ(G(jω))\sup_\omega \bar\sigma(G(j\omega)).
Which theorem links H∞ synthesis to robustness?
The Small-Gain Theorem: MMΔ\Delta loop is stable for all Δ1\|\Delta\|_\infty\le1 iff M<1\|M\|_\infty<1.
Why does bounding Tzw<γ\|T_{zw}\|_\infty<\gamma help?
It guarantees no exogenous disturbance is amplified into the performance outputs by more than γ\gamma, for the whole uncertainty set.
What mathematical tool converts time-domain energy to frequency domain in the derivation?
Parseval's theorem, w22=12πW(jω)2dω\|w\|_2^2=\frac{1}{2\pi}\int|W(j\omega)|^2 d\omega.
What do weighting functions WiW_i encode?
The designer's frequency-dependent specs (tracking, control effort, noise rejection); e.g. WSS<1\|W_S S\|_\infty<1 forces S<1/WS|S|<1/|W_S|.
Why can't sensitivity S|S| be small at all frequencies?
Bode/waterbed integral lnSdω=0\int\ln|S|\,d\omega=0: reducing it in one band raises it in another.
For a plant with resonance ζ=0.1\zeta=0.1, roughly how big is the resonant peak factor?
About 1/(2ζ)=51/(2\zeta)=5 times the low-frequency gain, dominating G\|G\|_\infty.
How is a ±30% gain uncertainty turned into a small-gain requirement?
Write it as 0.3Δ~0.3\tilde\Delta, Δ~1\|\tilde\Delta\|_\infty\le1; robust stability needs 0.3M<1\|0.3 M\|_\infty<1, i.e. M<3.33\|M\|_\infty<3.33.

Connections

  • LQR and LQG control — optimal for the known model; H∞ adds worst-case robustness.
  • Sensitivity and Complementary Sensitivity (S+T=1) — the objects H∞ weights shape.
  • Small-Gain Theorem · μ-synthesis (structured uncertainty)
  • Singular Value Decomposition — where σˉ\bar\sigma comes from (MIMO gain).
  • Parseval's Theorem · Bode Sensitivity Integral (waterbed)
  • Model Uncertainty in GNC — gusts, flex modes, mass variation for launch vehicles.

Concept Map

causes

motivates

guarantees over

handles only

weakness

implemented by

minimises

measured with

equals

defined as

linked via

derives

Model uncertainty

Nominal controller may go unstable

Robust control

Plant set P

LQR design

Nominal plant P0

H-infinity control

Worst-case energy gain

Signal 2-norm energy

H-infinity norm

Peak singular value over frequency

Parseval theorem

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, real rocket ya drone kabhi bhi apne math model ke exactly barabar nahi hota. Wind gust aayega, mass change hoga, kuch bending modes model me the hi nahi. Agar tumne controller sirf "nominal" (perfect) model pe tune kiya, toh real hardware pe woh unstable ho sakta hai. H∞ control ka idea simple hai: sirf average case mat sochoin, balki worst-case — sabse kameena disturbance — pe bhi system safe aur accha rahe. Yehi robustness hai.

Core quantity hai H∞ norm, G\|G\|_\infty. Yeh koi average gain nahi hai — yeh saare frequencies me se sabse bada amplification hai (peak of the Bode magnitude, MIMO me largest singular value). Kyun? Kyunki worst disturbance apni saari energy us frequency pe daal degi jahan system ka gain sabse zyada hai. Parseval theorem se hum time-domain energy ko frequency me le jaate hain, aur prove karte hain ki worst-case energy gain exactly G\|G\|_\infty ke barabar hai.

Robustness ka guarantee Small-Gain Theorem deta hai: agar uncertainty ko ek bounded block Δ\Delta (Δ1\|\Delta\|_\infty\le1) maano jo loop me lagi hai, toh loop tabhi tak stable rahega jab tak M<1\|M\|_\infty<1. Loop ke around signal har baar MΔM\Delta se multiply hota hai; agar yeh product 1 se chhota hai toh geometric series converge kar jaati hai, matlab bounded, matlab stable. Isliye H∞ design ka goal: Tzw\|T_{zw}\|_\infty ko γ\gamma ke neeche push karo — tab koi bhi disturbance γ\gamma se zyada amplify nahi hogi, aur infinite plants ka poora family ek hi inequality se certify ho jaata hai.

Practical me hum weighting functions lagate hain — yeh batate hain ki kaunsi frequency pe kya important hai (tracking low frequency pe, control effort high pe). WSS<1\|W_S S\|_\infty<1 ka matlab hai S<1/WS|S| < 1/|W_S| — matlab tum ek ceiling draw karte ho jiske neeche error ko rehna padega. Bas yaad rakho: gain infinite nahi bana sakte, kyunki waterbed effect (Bode integral) hai — ek jagah dabaoge toh doosri jagah uthega. H∞ isi trade-off ka best solution nikaalta hai.

Go deeper — visual, from zero

Test yourself — Guidance, Navigation & Control (GNC)

Connections