3.5.42 · D3 · Physics › Guidance, Navigation & Control (GNC) › Gain margin, phase margin — stability margins
Yeh page stability margins ke liye drill hall hai. Parent note ne do rulers banaye the — gain margin (GM) aur phase margin (PM) — critical point − 1 se. Yahan hum har tarah ke cases un rulers pe throw karte hain taaki koi bhi exam question, koi bhi sign, koi bhi degenerate loop tumhe surprise na kar sake.
Shuru karne se pehle, teen reminders plain words mein taaki koi bhi line one se lost na ho:
Definition Woh vocabulary jo hum use karte rahenge
s — Laplace variable. Isse "frequency behaviour" ka stand-in samjho. Ek transfer function L ( s ) ek formula hai s mein. Yeh poochna ke "loop kisi pure oscillation ka ω frequency par kaise respond karta hai?", hum set karte hain s = j ω , jahan j = − 1 algebra ko ek rotating wave mein badal deta hai. Toh L ( j ω ) bas L ( s ) hai jo oscillation axis par evaluate kiya gaya ho.
L ( j ω ) — open-loop response : ek complex number jisme ek size ∣ L ∣ hai (loop kitna amplify karta hai) aur ek angle ∠ L hai (wave ko kitna delay karta hai, degrees mein).
ω g c — gain crossover : woh ω jahan size 1 ke barabar ho.
ω p c — phase crossover : woh ω jahan angle − 18 0 ∘ ke barabar ho.
GM = 1/∣ L ( j ω p c ) ∣ aur PM = 18 0 ∘ + ∠ L ( j ω g c ) .
Is topic ka har question in cells mein se kisi ek mein aata hai. Neeche har worked example us cell ke saath tagged hai jo woh cover karta hai.
Cell
Kya cheez isse alag banati hai
Example
A. Clean GM
Loop ka ek real phase crossover hai; GM compute karo
Ex 1
B. Clean PM
Loop $
L
C. Dono margins saath
Usi loop par, GM aur PM dono nikalো
Ex 3
D. Infinite GM (degenerate)
Phase kabhi − 18 0 ∘ tak nahi pahunchti ⟹ GM = ∞
Ex 4
E. Negative margin (unstable)
Gain bahut zyada push ho gayi ⟹ GM < 0 dB, PM < 0
Ex 5
F. Time-delay twist
Delay phase margin kha jaati hai; τ ma x dhundho
Ex 6
G. Design / word problem
Target PM hit karne ke liye K choose karo (overshoot se link)
Ex 7
H. Exam trap
Conditionally stable / sign confusion
Ex 8
Worked example Ex 1 — Cell A: clean gain margin
Statement. L ( s ) = s ( s + 1 ) ( s + 2 ) 4 ke liye, gain margin dB mein nikalो.
Forecast: Andaaza lagao — kya GM yahan parent ke Worked Example 1 se bada hai ya chota (jahan K = 1 ne GM = 6 diya tha)? Aage padhne se pehle apna guess likhो.
ω p c dhundho (jahan ∠ L = − 18 0 ∘ ho).
∠ L = − 9 0 ∘ − arctan ω − arctan ( ω /2 ) . = − 18 0 ∘ set karne par arctan ω + arctan ( ω /2 ) = 9 0 ∘ milta hai.
Yeh step kyun? GM exactly wahan padhna padta hai jahan phase already danger angle par baith chuki ho — tab bas extra size chahiye − 1 tak pahunchne ke liye.
Jab do arctans 9 0 ∘ mein add hote hain tab sum ka tangent blow up ho jaata hai (infinity tak); tan ( a + b ) = 1 − t a n a t a n b t a n a + t a n b use karte hue, denominator 1 − ω ⋅ 2 ω = 0 ho jaata hai, toh ω p c = 2 .
Wahan size. ∣ L ∣ = ω ω 2 + 1 ω 2 + 4 4 . ω = 2 par: denominator = 2 ⋅ 3 ⋅ 6 = 6 , toh ∣ L ∣ = 6 4 = 3 2 .
Yeh step kyun? ω p c par ∣ L ∣ woh ek number hai jis par GM depend karta hai.
GM = ∣ L ∣ 1 = 2 3 = 1.5 , yaani 20 log 10 1.5 ≈ 3.52 dB.
Yeh step kyun? Hum ∣ L ∣ ko invert karte hain kyunki GM woh factor hai jisse tum gain multiply kar sakte ho before ∣ L ∣ danger phase par 1 tak pahunche: agar ∣ L ∣ = 3 2 hai, toh 2 3 se multiply karne par woh exactly 1 par land karta hai (yaani − 1 par). Bada ∣ L ∣ ⟹ disaster ke zyada paas ⟹ chota GM, aur yahi divide karna karta hai.
Verify: Parent ke K = 1 (GM 6 ) se gain 4 × badh gayi. K multiply karne se ∣ L ∣ scale up hota hai, toh GM ko 4 se divide hona chahiye: 6/4 = 1.5 . ✓ Match karta hai — aur woh chota hai, kyunki badi gain GM ko shrink karti hai.
Worked example Ex 2 — Cell B: clean phase margin
Statement. L ( s ) = s ( s + 1 ) 1 ke liye, phase margin nikalо (parent Ex 2 ka repeat hai lekin figure se verify karte hain).
Forecast: Parent ne PM ≈ 51. 8 ∘ nikala tha. Kya yeh "healthy" hai? Yes/no guess karo.
Gain crossover. ∣ L ∣ = ω ω 2 + 1 1 = 1 ⇒ ω 4 + ω 2 − 1 = 0 . u = ω 2 let karo: u = 2 − 1 + 5 ≈ 0.618 , toh ω g c = 0.618 ≈ 0.786 .
Yeh step kyun? PM wahan padhte hain jahan size exactly 1 ho; tab bas extra phase lag chahiye − 18 0 ∘ tak pahunchne ke liye.
Wahan phase. ∠ L = − 9 0 ∘ − arctan ( 0.786 ) = − 9 0 ∘ − 38. 2 ∘ = − 128. 2 ∘ .
PM = 18 0 ∘ + ( − 128. 2 ∘ ) = 51. 8 ∘ .
Figure dekho: loop curve unit circle (size = 1 ) ko amber dot par cross karti hai; us ray se − 18 0 ∘ axis tak woh amber wedge hi phase margin hai.
Verify: tan ( 38. 2 ∘ ) = 0.787 ≈ ω g c . ✓ Aur 3 0 ∘ –6 0 ∘ healthy hai, toh haan. ✓
Worked example Ex 3 — Cell C: ek hi loop par dono margins
Statement. L ( s ) = s ( s + 1 ) ( s + 2 ) 2 ke liye, dono GM (dB) aur PM nikalо.
Forecast: K = 2 ke saath (parent ke K = 1 ka double), guess karo kya dono margins abhi bhi positive hain.
GM. Wahi ω p c = 2 (phase K par depend nahi karti). ∣ L ∣ = 6 2 = 3 1 , toh GM = 3 = 20 log 10 3 ≈ 9.54 dB.
Yeh step kyun? K sirf size curve ko slide karta hai; danger frequency unchanged rehti hai.
PM — ω g c nikalо. Solve karo ∣ L ∣ = ω ω 2 + 1 ω 2 + 4 2 = 1 . Numerically ω g c ≈ 0.749 .
Yeh step kyun? Size= 1 wahi jagah hai jahan PM rehta hai.
Wahan phase. ∠ L = − 9 0 ∘ − arctan ( 0.749 ) − arctan ( 0.749/2 ) = − 9 0 ∘ − 36.8 3 ∘ − 20.5 3 ∘ = − 147.3 6 ∘ .
PM = 18 0 ∘ − 147.3 6 ∘ = 32. 6 ∘ .
Verify: Dono positive ⟹ stable, GM = 3 > 1 ke saath consistent. Step 2 ka sanity check: ω = 0.749 ko ∣ L ∣ mein plug karo ⟹ ≈ 1.00 . ✓ (Yeh dono numbers sidha do stacked graphs se padhne ke liye Bode plots dekho.)
Worked example Ex 4 — Cell D: infinite gain margin (degenerate)
Statement. First-order loop L ( s ) = s + 1 5 ke liye GM aur PM nikalо.
Forecast: Instability se pehle kitni extra gain? "Koi finite factor" vs "unlimited" guess karo.
Phase range. ∠ L = − arctan ω . Jaise ω : 0 → ∞ hota hai, ∠ L 0 ∘ se neeche − 9 0 ∘ ki taraf jaata hai lekin kabhi − 18 0 ∘ tak nahi pahunchta .
Yeh step kyun? GM − 18 0 ∘ par rehta hai — agar phase wahan kabhi nahi pahunchi, toh koi ω p c nahi hai.
GM = ∞ (ya "∞ dB"): koi bhi finite gain multiple ek aisi curve ko — jo − 9 0 ∘ par ruk jaati hai — negative-real axis − 1 par nahi pahuncha sakta.
PM. Size ∣ L ∣ = ω 2 + 1 5 = 1 ⇒ ω 2 + 1 = 25 ⇒ ω g c = 24 ≈ 4.899 . Phase = − arctan ( 4.899 ) = − 78.4 6 ∘ . PM = 18 0 ∘ − 78.4 6 ∘ = 101. 5 ∘ .
Verify: Ek single pole kabhi bhi ek simple loop ko gain akele unstable nahi kar sakta — classic result, GM = ∞ . ✓ PM bahut bada (> 9 0 ∘ ) ⟹ extremely damped, Root locus ke saath match karta hai jo sab K ke liye left half-plane mein rehta hai. ✓
Worked example Ex 5 — Cell E: negative margin (pehle se unstable)
Statement. Ex 1 ka plant lo lekin gain crank karo: L ( s ) = s ( s + 1 ) ( s + 2 ) 12 . GM aur PM nikalо; kya loop stable hai?
Forecast: Hum jaante hain K = 6 exact edge hai (parent Ex 1: K = 1 par GM = 6 ). Humne K = 12 choose kiya. GM ka sign guess karo.
GM. ω p c = 2 unchanged. ∣ L ∣ = 6 12 = 2 , toh GM = 2 1 = 0.5 = 20 log 10 0.5 ≈ − 6.02 dB.
Yeh step kyun? GM < 1 (negative dB) instability ka algebraic signature hai: curve ne − 1 ko pehle se swallow kar liya hai.
PM. ∣ L ∣ = 1 solve karo: ω ω 2 + 1 ω 2 + 4 12 = 1 ⇒ ω g c ≈ 1.96 . Phase = − 9 0 ∘ − arctan ( 1.96 ) − arctan ( 0.98 ) = − 9 0 ∘ − 62.9 7 ∘ − 44.4 2 ∘ = − 197. 4 ∘ . PM = 18 0 ∘ − 197. 4 ∘ ≈ − 17. 4 ∘ .
Yeh step kyun? Negative PM confirm karta hai ke phase gain ke 1 tak girne se pehle hi danger angle cross kar chuki hai.
Verdict. GM < 0 dB aur PM < 0 ⟹ unstable .
Verify: Critical gain K = 6 hai (GM= 1 ). Humne 12 = 2 × 6 use kiya, toh GM 1/2 hona chahiye. ✓ Dono margins negative ⟹ Nyquist stability criterion − 1 ko encircled dikhata. ✓
Worked example Ex 6 — Cell F: time-delay twist
Statement. Stable loop L ( s ) = s ( s + 1 ) 1 (Ex 2: PM = 51. 8 ∘ , ω g c = 0.786 ) ko ek transport delay τ milta hai (jaise ek communication lag). Woh maximum τ kya hai jo ise stable rakhega?
Forecast: Delay sirf phase lag add karta hai, size mein koi change nahi. Guess karo kya τ ma x second ka kuch fraction hai ya kaafi seconds.
Delay kya karta hai. Ek delay e − s τ ki size 1 hoti hai (magnitude curve ko kabhi nahi chhuta) lekin phase − ω τ radians add karta hai.
Yeh step kyun? Kyunki size untouched hai, ω g c 0.786 par rehta hai — extra lag bas PM ko wahan khaata hai.
Saari margin kha jao. Instability us moment aati hai jab added lag poori PM ke barabar ho jaata hai:
τ ma x = ω g c PM in radians .
PM = 51. 8 ∘ = 51.8 × 180 π = 0.9040 rad.
Compute karo. τ ma x = 0.786 0.9040 = 1.150 s.
Figure mein amber PM wedge dikhti hai jo delay ke growing lag se nibble hoti rehti hai jab tak woh khatam nahi ho jaati.
Verify: Units: rad/s rad = s . ✓ Numeric: 0.786 × 1.150 = 0.904 rad = 51. 8 ∘ extra lag = exactly PM. ✓ e − s τ ko rational function ke roop mein approximate karne ke liye Time delay and Padé approximation dekho.
Worked example Ex 7 — Cell G: design / word problem (gain choose karo)
Statement. Ek satellite reaction-wheel loop model hai L ( s ) = s ( s + 2 ) K . Mission spec: PM = 4 5 ∘ (ek well-damped, low-overshoot response ke liye). K nikalо.
Forecast: Ek handy rule hai (neeche justify ki gayi) ke PM degrees mein 100 se divide karne par damping ratio ζ ka estimate milta hai, toh PM = 4 5 ∘ targets ζ ≈ 0.45 . Guess karo kya required K 4 se upar hai ya neeche.
PM = 4 5 ∘ impose karo. PM = 18 0 ∘ + ∠ L ( j ω g c ) = 4 5 ∘ ⇒ ∠ L ( j ω g c ) = − 13 5 ∘ .
∠ L = − 9 0 ∘ − arctan ( ω /2 ) = − 13 5 ∘ ⇒ arctan ( ω /2 ) = 4 5 ∘ ⇒ ω g c = 2 .
Yeh step kyun? Phase K par depend nahi karta, toh target angle pehle ω g c fix karta hai.
Us ω g c par size = 1 impose karo. ∣ L ∣ = ω ω 2 + 4 K = 1 at ω = 2 : 2 ⋅ 8 K = 1 ⇒ K = 2 8 = 4 2 ≈ 5.657 .
Yeh step kyun? Ab K (K hi ek free knob hai) use karo size curve ko slide karne ke liye taaki unity gain exactly ω g c = 2 par land kare.
ζ ≈ PM ∘ /100 " kahan se aata hai?
Damping ratio ζ (dekho Damping ratio and overshoot ) measure karta hai ke second-order response kitna oscillatory hai: ζ = 0 hamesha ring karta hai, ζ = 1 just-non-oscillatory hai. Ek standard second-order loop ke liye PM aur ζ ko link karne wala ek exact formula hai; useful range PM ≲ 6 0 ∘ mein woh curve origin se nearly ek straight line hai jiska slope roughly 1/100 per degree hai, jo pocket rule ζ ≈ PM ∘ /100 deta hai. Yeh ek approximation hai, law nahi — gain choose karne ke liye kaafi hai, phir simulation se confirm karo.
Verify: K = 5.657 par PM check karo: ω g c = 2 , phase = − 9 0 ∘ − arctan ( 1 ) = − 13 5 ∘ , PM = 4 5 ∘ . ✓ ζ ≈ 0.45 se expected overshoot roughly 20% hai (dekho Damping ratio and overshoot ); overshoot trim karne ke liye bina speed khoye, PID controller tuning ya Loop shaping se reshape karo. K ≈ 5.66 > 4 , forecast direction se match karta hai.
Worked example Ex 8 — Cell H: exam trap (sign / "which frequency" confusion)
Statement. Ek student likhta hai: "L ( s ) = s ( s + 1 ) ( s + 2 ) 10 ke liye, ∠ L ( j ω g c ) = − 15 0 ∘ hai, isliye PM = − 15 0 ∘ hai aur system unstable hai." Sahi PM nikalо aur galti expose karo.
Forecast: Guess karo kya − 15 0 ∘ ki phase ek positive ya negative PM matlab hai.
Trap pakdo. PM phase khud nahi hai ; yeh − 18 0 ∘ tak ka gap hai: PM = 18 0 ∘ + ∠ L ( j ω g c ) .
Yeh step kyun? Sabse common mistake (parent ka Mistake #4) ∠ L ko PM samajhna hai.
Sahi value. PM = 18 0 ∘ + ( − 15 0 ∘ ) = + 3 0 ∘ . Positive ⟹ stable (barely; ringing likely).
Doosra trap — which frequency? Agar student ne PM padhne ke liye ω p c (jahan ∠ = − 18 0 ∘ ) use kiya hota, toh unhe galat 0 ∘ milta. Mnemonic "GaP" yaad karo: margin wahan measure hota hai jahan doosri quantity apni critical value hit kare.
Verify: 180 + ( − 150 ) = 30 > 0 ⟹ stable, student se contradict karta hai. ✓ − 15 0 ∘ ki phase − 18 0 ∘ danger point se kam lag hai, toh abhi bhi positive margin honi chahiye — geometry algebra confirm karta hai.
Recall Rapid self-test
Phase crossover par tum kaun sa margin padhte ho? ::: Gain margin (GM = 1/∣ L ( j ω p c ) ∣ ).
Ek first-order loop s + 1 K ka gain margin kya hota hai? ::: Infinite — phase kabhi − 18 0 ∘ tak nahi pahunchti.
Ek delay size add karta hai ya phase? ::: Sirf phase, − ω τ ; size 1 rehti hai.
Agar PM ≈ − 17. 4 ∘ aur GM = − 6 dB, kya loop stable hai? ::: Nahi — dono negative matlab pehle se unstable.
PM = 4 5 ∘ hit karne ke liye tum pehle kaun si quantity fix karte ho? ::: Gain-crossover frequency, target phase se (phase K se independent hai).
Student kehta hai PM = ∠ L ( j ω g c ) . Fix karo? ::: PM = 18 0 ∘ + ∠ L ( j ω g c ) — gap, angle nahi.
"GaP" — G M at P hase crossover, P M at G ain crossover. Angle mein 180 add karo PM paane ke liye. Delay phase khaata hai, size kabhi nahi.