GATE GUIDE

GATE EE Control Systems Formula Sheet: Key Formulas and Traps

By MD ANISH AHAMADUpdated 4 Oct 202610 min read
GATE EE Control Systems Formula Sheet: Key Formulas and Traps

The key GATE EE control systems formulas are the closed-loop gain G/(1+GH)G/(1+GH), Mason's formula, the second order specifications, the error constants, the Routh conditions, the root-locus rules, the margins, Z=N+PZ = N + P, the compensator results and C(sI−A)−1B+DC(sI-A)^{-1}B + D. Each comes with a condition that decides whether it applies. This sheet teaches them in learning order, with a worked example and a trap for each part.

In this guide
  1. Key takeaways
  2. The terms this sheet uses
  3. Transfer functions, feedback and Mason's gain formula
  4. First and second order time response
  5. Steady-state error and the error constants
  6. Routh-Hurwitz: counting unstable roots
  7. Root locus: where the closed-loop poles travel
  8. Bode plots, margins and frequency-domain specifications
  9. Nyquist: counting encirclements
  10. Compensators and P, PI and PID controllers
  11. State space and the state transition matrix
  12. Quick revision

Key takeaways

The terms this sheet uses

The open-loop transfer function G(s)H(s)G(s)H(s) is the gain around the broken loop, with GG forward and HH in the feedback path. The closed-loop transfer function T(s)T(s) is output over input with the loop closed. Unity feedback means H=1H = 1.

The type is the number of poles of G(s)H(s)G(s)H(s) at the origin. A second order system has a natural frequency ωn\omega_n and a damping ratio ζ\zeta; it is underdamped when 0<ζ<10 < \zeta < 1.

The gain crossover frequency ωgc\omega_{gc} is where the open-loop magnitude is 0 dB. The phase crossover frequency ωpc\omega_{pc} is where the open-loop phase is −180∘-180^\circ.

Transfer functions, feedback and Mason's gain formula

A transfer function is output over input in the Laplace domain with zero initial conditions. Negative feedback divides the forward gain by one plus the loop gain.

Formula Watch out for
T=G1+GH\displaystyle T = \dfrac{G}{1 + GH} (negative feedback); T=G1−GH\displaystyle T = \dfrac{G}{1 - GH} (positive) Follow the summing junction sign
SGT=11+GH\displaystyle S_G^T = \dfrac{1}{1 + GH}; SHT=−GH1+GH\displaystyle S_H^T = \dfrac{-GH}{1 + GH} High loop gain makes TT sensitive to HH
Cascade multiplies, parallel adds Moving a take-off point past a block changes that branch
T=∑kPkΔkΔ\displaystyle T = \dfrac{\sum_k P_k \Delta_k}{\Delta} Δk\Delta_k drops every loop touching path kk
Δ=1−∑Li+∑LiLj−∑LiLjLk+…\Delta = 1 - \sum L_i + \sum L_i L_j - \sum L_i L_j L_k + \dots Products of non-touching loops only

Worked example. With G=20/(s+4)G = 20/(s+4) and H=1H = 1, T=20/(s+24)T = 20/(s+24), the DC gain is 20/24=0.83320/24 = 0.833 and the DC sensitivity 1/(1+5)=0.1671/(1+5) = 0.167.

Worked example (Mason). One path of gain 6 touches two non-touching loops of gains −2-2 and −3-3: Δ=1−(−5)+6=12\Delta = 1 - (-5) + 6 = 12, Δ1=1\Delta_1 = 1, so T=6/12=0.5T = 6/12 = 0.5.

The trap is the sign: loop gains carry their own minus signs, so subtracting them adds them back.

First and second order time response

A first order system rises with one time constant τ\tau. A second order system's speed is set by ωn\omega_n and its overshoot by ζ\zeta.

Formula Watch out for
1τs+1\displaystyle \dfrac{1}{\tau s + 1}: step response 1−e−t/τ1 - e^{-t/\tau} 63.2 per cent at t=τt = \tau
tr=2.2τt_r = 2.2\tau (10 to 90 per cent); ts=4τt_s = 4\tau (2 per cent) trt_r is τln⁡9\tau \ln 9, not τ\tau
ωn2s2+2ζωns+ωn2\displaystyle \dfrac{\omega_n^2}{s^2 + 2\zeta\omega_n s + \omega_n^2}; poles −ζωn±jωn1−ζ2-\zeta\omega_n \pm j\omega_n\sqrt{1-\zeta^2} Read ωn2\omega_n^2 from the constant term
ωd=ωn1−ζ2\omega_d = \omega_n\sqrt{1-\zeta^2}; tp=π/ωdt_p = \pi/\omega_d Uses ωd\omega_d, not ωn\omega_n
Mp=100 e−πζ/1−ζ2M_p = 100\,e^{-\pi\zeta/\sqrt{1-\zeta^2}} per cent Depends on ζ\zeta alone
ts=4/(ζωn)t_s = 4/(\zeta\omega_n) (2 per cent); 3/(ζωn)3/(\zeta\omega_n) (5 per cent) Check the stated band
tr=(π−cos⁡−1ζ)/ωdt_r = (\pi - \cos^{-1}\zeta)/\omega_d (0 to 100 per cent) A different rise-time definition

Worked example. T=25/(s2+6s+25)T = 25/(s^2 + 6s + 25) gives ωn=5\omega_n = 5, ζ=6/10=0.6\zeta = 6/10 = 0.6, ωd=5×0.8=4\omega_d = 5 \times 0.8 = 4 rad/s, tp=π/4=0.785t_p = \pi/4 = 0.785 s, Mp=100 e−0.75π=9.48M_p = 100\,e^{-0.75\pi} = 9.48 per cent and ts=4/3=1.33t_s = 4/3 = 1.33 s.

Remember: settling time is fixed by how far left the poles sit, and overshoot by the angle they make with the negative real axis. Equal real parts mean equal settling times.

The trap is a reduced model. You may drop a pole far to the left of the dominant pair, but keep the DC gain unchanged.

Steady-state error and the error constants

Steady-state error is what remains once the transient dies out. For unity feedback it follows from the type and three constants.

Formula Watch out for
Kp=lim⁡s→0G(s)K_p = \lim_{s \to 0} G(s); step error 11+Kp\displaystyle \dfrac{1}{1 + K_p} Finite only for type 0
Kv=lim⁡s→0sG(s)K_v = \lim_{s \to 0} sG(s); ramp error 1Kv\displaystyle \dfrac{1}{K_v} Type 1: zero step error, finite ramp error
Ka=lim⁡s→0s2G(s)K_a = \lim_{s \to 0} s^2G(s); parabola error 1Ka\displaystyle \dfrac{1}{K_a} Type 2: zero step and ramp errors
Input of size AA: multiply the unit error by AA A ramp 5t5t gives 5/Kv5/K_v
Sinusoid: output ∣T(jω)∣sin⁡(ωt+∠T(jω))\lvert T(j\omega) \rvert \sin(\omega t + \angle T(j\omega)) Use the closed-loop TT, not GG

Worked example. For G=40(s+2)/(s(s+4)(s+5))G = 40(s+2)/\big(s(s+4)(s+5)\big) the closed loop s3+9s2+60s+80s^3 + 9s^2 + 60s + 80 is stable since 9×60>809 \times 60 > 80, so Kv=80/20=4K_v = 80/20 = 4 and the unit ramp error is 0.250.25.

The trap is non-unity feedback. Reduce it to a unity feedback form before taking the limits.

Routh-Hurwitz: counting unstable roots

The Routh array counts right-half-plane roots without solving the polynomial. Each first-column sign change marks one.

Formula Watch out for
Sign changes in column one == right-half-plane roots Positive coefficients are necessary, not sufficient
s3+a2s2+a1s+a0s^3 + a_2 s^2 + a_1 s + a_0 stable iff all positive and a2a1>a0a_2 a_1 > a_0 Build the full array for higher orders
Zero in column one: replace it by ε\varepsilon Let ε→0+\varepsilon \to 0^+ before counting
Row of zeros: auxiliary polynomial from the row above Its roots give the oscillation frequency

Worked example. For G=K/((s+1)(s+2)(s+3))G = K/\big((s+1)(s+2)(s+3)\big) the characteristic equation is s3+6s2+11s+6+K=0s^3 + 6s^2 + 11s + 6 + K = 0, so stability needs 66>6+K66 > 6 + K and 6+K>06 + K > 0, giving −6<K<60-6 < K < 60; at K=60K = 60, 6s2+66=06s^2 + 66 = 0 gives oscillation at 11=3.32\sqrt{11} = 3.32 rad/s.

Trap: a right-half-plane pole or zero in the plant can make closed-loop coefficients change sign with the gain. Always build the closed-loop polynomial, then the array.

Root locus: where the closed-loop poles travel

The root locus traces the closed-loop poles as KK rises from zero. Here nn counts open-loop poles and mm zeros.

Formula Watch out for
Starts at the poles (K=0K = 0), ends at the zeros or at infinity n−mn - m branches go to infinity
Real axis: left of an odd number of real poles and zeros Count only those to the right of the point
Asymptote angles (2q+1)180∘n−m\displaystyle \dfrac{(2q+1)180^\circ}{n - m}; centroid ∑pi−∑zjn−m\displaystyle \dfrac{\sum p_i - \sum z_j}{n - m} Use signed pole values
Breakaway where dKds=0\displaystyle \dfrac{dK}{ds} = 0 Keep only roots on the locus
Angle condition ∠G(s)H(s)=±180∘(2q+1)\angle G(s)H(s) = \pm 180^\circ(2q+1); gain K=∏∣s−pi∣∏∣s−zj∣\displaystyle K = \dfrac{\prod \lvert s - p_i \rvert}{\prod \lvert s - z_j \rvert} Factors must have unit leading coefficients
Departure from a complex pole =180∘−∑(other pole angles)+∑(zero angles)= 180^\circ - \sum \text{(other pole angles)} + \sum \text{(zero angles)} Arrival at a zero swaps the roles

Worked example. For K/(s(s+1)(s+5))K/\big(s(s+1)(s+5)\big) the centroid is (0−1−5)/3=−2(0 - 1 - 5)/3 = -2, the asymptotes lie at 60∘60^\circ, 180∘180^\circ and 300∘300^\circ, 3s2+12s+5=03s^2 + 12s + 5 = 0 gives breakaway at s=−0.473s = -0.473, and Routh gives the imaginary-axis crossing at K=30K = 30, ω=5=2.24\omega = \sqrt{5} = 2.24 rad/s.

The trap is the other breakaway root, s=−3.53s = -3.53. It lies between −1-1 and −5-5, where there is no locus, so reject it.

Bode plots, margins and frequency-domain specifications

A Bode plot draws magnitude in dB and phase against log frequency. Each factor adds a straight-line asymptote.

Formula Watch out for
Each pole −20-20 dB/decade and −90∘-90^\circ; each zero +20+20 dB/decade and +90∘+90^\circ −45∘-45^\circ and a 3 dB error at a first-order corner
Quadratic pole pair −40-40 dB/decade, resonant peak if ζ<0.707\zeta < 0.707 The asymptote hides the peak
Type 1: initial slope meets 0 dB at ω=Kv\omega = K_v; type 0: level 20log⁡10Kp20\log_{10} K_p Extend the initial line past any corner
e−sTe^{-sT}: magnitude 1, phase −ωT-\omega T rad Convert to degrees before adding
GM=1∣GH(jωpc)∣\displaystyle \text{GM} = \dfrac{1}{\lvert GH(j\omega_{pc}) \rvert}; in dB −20log⁡10∣GH(jωpc)∣-20\log_{10}\lvert GH(j\omega_{pc}) \rvert Read at the phase crossover
PM=180∘+∠GH(jωgc)\text{PM} = 180^\circ + \angle GH(j\omega_{gc}) Read at the gain crossover
Mr=12ζ1−ζ2\displaystyle M_r = \dfrac{1}{2\zeta\sqrt{1-\zeta^2}} at ωr=ωn1−2ζ2\omega_r = \omega_n\sqrt{1 - 2\zeta^2} Valid only for ζ<0.707\zeta < 0.707
ωb=ωn1−2ζ2+4ζ4−4ζ2+2\omega_b = \omega_n\sqrt{1 - 2\zeta^2 + \sqrt{4\zeta^4 - 4\zeta^2 + 2}} Closed-loop bandwidth

Worked example (gain margin). For 10/(s(s+1)(s+5))10/\big(s(s+1)(s+5)\big) the phase is −180∘-180^\circ at ωpc=5\omega_{pc} = \sqrt{5}, where ∣G∣=10/30\lvert G \rvert = 10/30, so GM=3\text{GM} = 3, or 9.549.54 dB, matching the Routh limit K=30K = 30.

Worked example (phase margin). For 2/(s(s+1))2/\big(s(s+1)\big), ω4+ω2−4=0\omega^4 + \omega^2 - 4 = 0 gives ωgc=1.25\omega_{gc} = 1.25 rad/s, so PM=180∘−90∘−tan⁡−11.25=38.7∘\text{PM} = 180^\circ - 90^\circ - \tan^{-1}1.25 = 38.7^\circ.

Worked example (specifications). With ζ=0.5\zeta = 0.5 and ωn=10\omega_n = 10 rad/s, Mr=1/(2×0.5×0.866)=1.155M_r = 1/(2 \times 0.5 \times 0.866) = 1.155 at ωr=7.07\omega_r = 7.07 rad/s, and ωb=101.618=12.72\omega_b = 10\sqrt{1.618} = 12.72 rad/s.

These answers are decimals, so read the virtual calculator guide before the exam.

Nyquist: counting encirclements

The Nyquist criterion counts closed-loop right-half-plane poles from the open-loop frequency response, even when the open loop is unstable.

Formula Watch out for
Z=N+PZ = N + P NN counts clockwise encirclements of −1+j0-1 + j0
Stable iff Z=0Z = 0, so anticlockwise encirclements must equal PP Anticlockwise counts as negative NN
Strictly proper GHGH: the infinite arc maps to the origin A proper GHGH maps it to lim⁡s→∞GH\lim_{s \to \infty} GH
Pole of order qq at the origin: an arc of q×180∘q \times 180^\circ at infinite radius, clockwise Depends on the side of the indentation

Worked example. For G=K/(s−1)G = K/(s-1), P=1P = 1, and for K>1K > 1 the plot circles −1-1 once anticlockwise, so N=−1N = -1 and Z=−1+1=0Z = -1 + 1 = 0: stable, matching the closed loop s−1+Ks - 1 + K.

In one line: Nyquist is bookkeeping, so write down PP first, count NN with its sign, and only then decide stability.

Compensators and P, PI and PID controllers

A lead compensator adds phase near crossover. A lag compensator raises low-frequency gain relative to crossover.

Formula Watch out for
Lead 1+aTs1+Ts\displaystyle \dfrac{1 + aTs}{1 + Ts}, a>1a > 1: sin⁡ϕm=a−1a+1\displaystyle \sin\phi_m = \dfrac{a-1}{a+1} at ωm=1Ta\displaystyle \omega_m = \dfrac{1}{T\sqrt{a}} ωm\omega_m is the geometric mean of the corners
Lead gain at ωm\omega_m: 10log⁡10a10\log_{10} a dB Half the high-frequency gain 20log⁡10a20\log_{10} a
Lag 1+Ts1+βTs\displaystyle \dfrac{1 + Ts}{1 + \beta Ts}, β>1\beta > 1: attenuation 20log⁡10β20\log_{10}\beta dB Cuts error, slows the response
PI Kp+Ki/sK_p + K_i/s: raises the type by one Removes step error, adds phase lag
PD Kp+KdsK_p + K_d s: phase lead tan⁡−1(Kdω/Kp)\tan^{-1}(K_d\omega/K_p) Amplifies noise
PID Kp+Ki/s+KdsK_p + K_i/s + K_d s: two zeros and a pole at the origin Lead and lag together form lead-lag

Worked example. With a=4a = 4 and T=0.05T = 0.05 s, sin⁡ϕm=3/5\sin\phi_m = 3/5 gives ϕm=36.9∘\phi_m = 36.9^\circ at ωm=1/(0.05×2)=10\omega_m = 1/(0.05 \times 2) = 10 rad/s, where the gain is 10log⁡104=6.0210\log_{10}4 = 6.02 dB.

The trap is the P controller alone. A higher KpK_p shrinks a type 0 step error but never removes it.

State space and the state transition matrix

The poles and the transfer function of a state-space model both come from sI−AsI - A.

Formula Watch out for
x˙=Ax+Bu\dot{x} = Ax + Bu, y=Cx+Duy = Cx + Du; G(s)=C(sI−A)−1B+DG(s) = C(sI - A)^{-1}B + D Keep the DD term
Poles == roots of det⁡(sI−A)=0\det(sI - A) = 0 Eigenvalues of AA
Φ(t)=eAt=L−1[(sI−A)−1]\Phi(t) = e^{At} = \mathcal{L}^{-1}\big[(sI - A)^{-1}\big] Check Φ(0)=I\Phi(0) = I
Φ(t)−1=Φ(−t)\Phi(t)^{-1} = \Phi(-t); Φ(t1+t2)=Φ(t1)Φ(t2)\Phi(t_1 + t_2) = \Phi(t_1)\Phi(t_2) A product, not a sum
x(t)=Φ(t)x(0)+∫0tΦ(t−τ)Bu(τ) dτx(t) = \Phi(t)x(0) + \int_0^t \Phi(t - \tau)Bu(\tau)\,d\tau Zero-input plus zero-state
z=Txz = Tx: Az=TAT−1A_z = TAT^{-1}, Bz=TBB_z = TB, Cz=CT−1C_z = CT^{-1} Eigenvalues and G(s)G(s) unchanged
Controllable iff rank [B   AB   ⋯   An−1B]=n\text{rank}\,[B \;\, AB \;\, \cdots \;\, A^{n-1}B] = n; observable iff rank [C; CA; ⋯ ; CAn−1]=n\text{rank}\,[C;\, CA;\, \cdots;\, CA^{n-1}] = n Full rank is the test

Worked example. Below, det⁡(sI−A)=s2+5s+6\det(sI - A) = s^2 + 5s + 6, so the poles are −2-2 and −3-3 and G(s)=1/(s2+5s+6)G(s) = 1/(s^2 + 5s + 6); partial fractions of (sI−A)−1(sI-A)^{-1} give Φ(t)\Phi(t).

A=[01−6−5],  B=[01],  C=[10],Φ(t)=[3e−2t−2e−3te−2t−e−3t−6e−2t+6e−3t−2e−2t+3e−3t]A = \begin{bmatrix} 0 & 1 \\ -6 & -5 \end{bmatrix},\; B = \begin{bmatrix} 0 \\ 1 \end{bmatrix},\; C = \begin{bmatrix} 1 & 0 \end{bmatrix},\qquad \Phi(t) = \begin{bmatrix} 3e^{-2t} - 2e^{-3t} & e^{-2t} - e^{-3t} \\ -6e^{-2t} + 6e^{-3t} & -2e^{-2t} + 3e^{-3t} \end{bmatrix}

At t=0t = 0 this is the identity. The trap is the direction of the transformation: x=Tzx = Tz gives T−1ATT^{-1}AT, while z=Txz = Tx gives TAT−1TAT^{-1}.

The last-minute sheet in the GATE EE 2027 book covers every syllabus topic this way, each result with its conditions. The book has 643 pages, 914 questions with worked solutions and 10 full mock tests, built on all 5 GATE EE papers from 2022 to 2026, counted question by question. For MCQ, MSQ and NAT marking, see the exam pattern and marking scheme, common to every paper, and the timetable in GATE 2027 exam dates.

More formula sheets: all of GATE EE · Electrical Machines · Power Electronics · Power Systems

Quick revision

  1. Negative feedback gives G/(1+GH)G/(1+GH); Mason subtracts the loop sum and adds non-touching products.
  2. Second order: tp=π/ωdt_p = \pi/\omega_d, overshoot from ζ\zeta alone, ts=4/(ζωn)t_s = 4/(\zeta\omega_n) at 2 per cent.
  3. Check stability, then use 1/(1+Kp)1/(1+K_p), 1/Kv1/K_v and 1/Ka1/K_a.
  4. A cubic is stable when all coefficients are positive and a2a1>a0a_2 a_1 > a_0.
  5. Root locus: centroid (∑p−∑z)/(n−m)(\sum p - \sum z)/(n-m), breakaway from dK/ds=0dK/ds = 0.
  6. Gain margin at the phase crossover, phase margin at the gain crossover.
  7. Nyquist: Z=N+PZ = N + P, clockwise NN positive, stable when Z=0Z = 0.
  8. Lead: sin⁡ϕm=(a−1)/(a+1)\sin\phi_m = (a-1)/(a+1); state-space poles from det⁡(sI−A)\det(sI - A).

Frequently asked questions

What is the formula for peak overshoot in a second order system?

For the standard underdamped second order system the peak overshoot is 100 e−πζ/1−ζ2100\,e^{-\pi\zeta/\sqrt{1-\zeta^2}} per cent. It depends on the damping ratio alone, not on the natural frequency. A damping ratio of 0.6 gives about 9.5 per cent, and 0.5 gives about 16.3 per cent. Use it only for a pure second order system with no zero, or when one pair of poles clearly dominates.

How do you find the steady-state error for a ramp input?

Compute the velocity error constant KvK_v as the limit of sG(s)sG(s) as ss tends to zero, for a unity feedback loop. The error to a unit ramp is 1/Kv1/K_v. A type 0 loop has infinite ramp error, type 1 a finite error and type 2 zero error. Check first that the closed loop is stable, because the error constants mean nothing for an unstable loop.

What is the Nyquist stability criterion formula?

The criterion is Z=N+PZ = N + P. Here PP is the number of open-loop poles in the right half plane, NN is the number of clockwise encirclements of the point minus one by the plot of G(jω)H(jω)G(j\omega)H(j\omega), and ZZ is the number of closed-loop right half plane poles. The closed loop is stable only when ZZ is zero. Anticlockwise encirclements count as negative.

What is the difference between gain margin and phase margin?

Gain margin is read at the phase crossover frequency, where the phase is minus 180 degrees, and equals one over the open-loop magnitude there. Phase margin is read at the gain crossover frequency, where the magnitude is one, and equals 180 degrees plus the phase there. For a minimum phase system the closed loop is stable when both margins are positive.

How do I revise control systems formulas before GATE EE?

Read a compact sheet daily in the last two weeks and attach one trap to each formula. Work one small example for each sub-topic by hand, such as a Routh array or a margin calculation, so the steps stay quick. Keep full precision on the virtual calculator and round only the final numerical answer.

Sources

Dates, fees and the syllabus are set by the GATE 2027 organising institute and can change. Always confirm at gate2027.iitm.ac.in.

Keep reading

GATE EE 2027 book643 pages · ₹250 ₹300
Buy now — ₹250