3.5.27Guidance, Navigation & Control (GNC)

Transfer function — Laplace domain, poles and zeros

1,809 words8 min readdifficulty · medium1 backlinks

WHY do we even need this?

A GNC system (a rocket's attitude loop, an autopilot) is described by linear differential equations. Solving them by hand for every new command is painful. We want a shortcut that answers: "if I command X, what does the vehicle do?"

The Laplace transform is that shortcut. It converts ddt\frac{d}{dt} into multiplication by ss, so a differential equation becomes a polynomial equation. The ratio of output-polynomial to input-polynomial is the transfer function.


WHAT is the Laplace transform (first principles)

WHY the kernel este^{-st}? Because este^{-st} is the eigenfunction of differentiation: ddtest=sest\frac{d}{dt}e^{-st} = -s\,e^{-st}. Multiplying by such an exponential and integrating "sorts" the signal by how fast it grows/decays (σ\sigma) and how fast it oscillates (ω\omega).

The property that does all the work — HOW ss replaces ddt\frac{d}{dt}

Derive the derivative rule from the definition using integration by parts: L{f(t)}=0f(t)estdt\mathcal{L}\{f'(t)\} = \int_0^\infty f'(t)e^{-st}\,dt Let u=estu=e^{-st}, dv=f(t)dtdu=sestdtdv=f'(t)dt \Rightarrow du=-se^{-st}dt, v=f(t)v=f(t): =[f(t)est]0+s0f(t)estdt=f(0)+sF(s)= \Big[f(t)e^{-st}\Big]_0^\infty + s\int_0^\infty f(t)e^{-st}\,dt = -f(0) + sF(s)

Why zero initial conditions? A transfer function describes the system itself, not one particular start state. We deliberately zero the ICs so the answer depends only on the machine.


Building the transfer function

HOW to get it from an ODE — worked example

System: y¨+3y˙+2y=u(t)\ddot{y} + 3\dot{y} + 2y = u(t) (a damped mass on a spring; uu = force).

Step 1 — Laplace each term, zero ICs. Why? Derivatives → powers of ss. s2Y+3sY+2Y=Us^2 Y + 3sY + 2Y = U

Step 2 — Factor out YY. Why? We want the ratio Y/UY/U. (s2+3s+2)Y=U(s^2+3s+2)Y = U

Step 3 — Form the ratio. G(s)=YU=1s2+3s+2=1(s+1)(s+2)G(s) = \frac{Y}{U} = \frac{1}{s^2+3s+2} = \frac{1}{(s+1)(s+2)}


Poles and Zeros — the roots that decide behaviour

WHY does the pole become epte^{pt}? Because L1{1sp}=ept\mathcal{L}^{-1}\{\frac{1}{s-p}\}=e^{pt}. Partial fractions split GG into a sum of Akspk\frac{A_k}{s-p_k}, and each piece inverts to AkepktA_k e^{p_k t}.

Figure — Transfer function — Laplace domain, poles and zeros

Example: read the poles of our system

G(s)=1(s+1)(s+2)G(s)=\frac{1}{(s+1)(s+2)} → poles at s=1s=-1 and s=2s=-2, no finite zeros. Both poles are in the left half-plane → response c1et+c2e2t\sim c_1 e^{-t}+c_2 e^{-2t}stable, no oscillation. The slow pole s=1s=-1 (closer to imaginary axis) dominates (time constant τ=1\tau=1 s vs 0.50.5 s).


Steel-man the common mistakes


Forecast-then-Verify checkpoint


Flashcards

What operation does multiplying by ss represent in the Laplace domain?
Differentiation w.r.t. time, L{f}=sF(s)f(0)\mathcal{L}\{f'\}=sF(s)-f(0).
Define a transfer function.
Ratio Y(s)/U(s)Y(s)/U(s) of output to input Laplace transforms, with zero initial conditions.
What are poles?
Roots of the denominator; frequencies where G(s)G(s)\to\infty; they set the system's natural modes epte^{pt}.
What are zeros?
Roots of the numerator; where G(s)=0G(s)=0; they shape the transient but not stability.
Stability condition in terms of poles?
All poles in the left half-plane (negative real part).
Physical meaning of a pole's real part?
Growth/decay rate (σ\sigma); negative = decaying/stable.
Physical meaning of a pole's imaginary part?
Oscillation frequency (ω\omega) of that mode.
Poles of 1(s+1)(s+2)\frac{1}{(s+1)(s+2)} and behaviour?
s=1,2s=-1,-2; stable, non-oscillatory, response c1et+c2e2tc_1e^{-t}+c_2e^{-2t}.
Why este^{-st} in the Laplace integral?
It's the eigenfunction of d/dtd/dt, turning derivatives into multiplication by ss.
Poles of s2+2ζωns+ωn2s^2+2\zeta\omega_n s+\omega_n^2?
s=ζωn±jωn1ζ2s=-\zeta\omega_n\pm j\omega_n\sqrt{1-\zeta^2}.
Recall Feynman: explain to a 12-year-old

Imagine a swing. If you push it once, it wiggles a bit and stops. The way it wiggles — how fast it swings and how quickly it calms down — is decided by the swing itself, not by how you pushed it. The "transfer function" is a card that tells you the swing's personality. The special numbers on that card are the poles: they say "this swing calms down in 2 seconds and swings twice a second." If a pole ever says "this thing gets bigger forever," the swing is broken (unstable). Zeros are like a mute button for certain pushes — they change the shape of the wiggle but can't break the swing.

Connections

  • Laplace transform — the tool that maps time → frequency domain
  • Block diagrams and feedback loops — transfer functions multiply along signal paths
  • Stability - Routh-Hurwitz criterion — checks pole locations without factoring
  • Root locus — how poles move as gain changes
  • Second-order systems - damping ratio and natural frequency
  • Bode plot — evaluate G(jω)G(j\omega) along the imaginary axis
  • State-space representation — equivalent time-domain description; eigenvalues = poles

Concept Map

hard to solve repeatedly

use

kernel e^-st sorts signal

converts d/dt

zero initial conditions

ratio output over input

equals

roots of numerator

roots of denominator

decide

shape

Linear differential equation

Need a shortcut

Laplace transform

Complex frequency s

Derivative rule s F s minus f0

ODE becomes polynomial equation

Transfer function G s

N s over D s

Zeros

Poles

System behaviour and stability

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, transfer function ka basic idea simple hai. Kisi bhi system (jaise rocket ka attitude control ya autopilot) ko hum differential equations se describe karte hain, lekin har baar equation solve karna bahut painful hai. Isliye hum Laplace transform ka jugaad use karte hain — yeh time domain se nikaal ke ek "frequency domain" mein le jaata hai jahan ddt\frac{d}{dt} (differentiation) sirf ss se multiply karne jaisa ho jaata hai. Matlab calculus ban gaya simple algebra! Phir output ka Laplace divide by input ka Laplace = transfer function G(s)G(s), jo poore system ki fingerprint hai.

Ab G(s)G(s) hamesha do polynomials ka ratio hota hai: N(s)/D(s)N(s)/D(s). Zeros matlab numerator zero, aur poles matlab denominator zero (jahan GG infinity ho jaata hai). Yaad rakho: poles hi asli boss hain. Har pole pp time response mein ek term epte^{pt} deta hai. Agar pole ka real part negative hai (left half-plane), toh ete^{-t} jaisa decay hoga — system stable. Agar real part positive (right half-plane), toh e+te^{+t} badhta jaayega — system unstable, phat jaayega. Imaginary part ho toh oscillation (jhoolna) aata hai.

Ek common galti: log sochte hain zeros bhi stability decide karte hain. Nahi bhai! Zeros sirf response ka shape badalte hain, stability sirf poles se aati hai. Isliye mantra yaad rakho: "Left is alright" — jab tak saare poles left half-plane mein hain, system theek hai. GNC mein yeh sab isliye important hai kyunki hum poles ko dekh ke turant bata sakte hain vehicle command ke baad settle karega ya oscillate karega ya diverge — bina puri equation solve kiye. Yeh 80/20 shortcut hai: pole locations dekho, kahani samajh jaao.

Go deeper — visual, from zero

Test yourself — Guidance, Navigation & Control (GNC)

Connections