22-Mec-A3 System Analysis and Control · December 2013
Question 6 of 6: Asymptotic Bode magnitude plot, phase margin and gain margin
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
Paper format. National Exams, December 2013 — 07-Mec-A3
System Analysis and Control. Three hours, closed book; a Casio or
Sharp approved calculator and semi-logarithmic graph paper are the only permitted aids.
Six questions are printed; the rubric states that any four questions constitute a
complete paper and that all questions are of equal value. A two-page table of Laplace
transform pairs is appended to the examination. Because this set is a study resource,
all six questions are solved in full below.
Reference texts.
K. Ogata, Modern Control Engineering, 5th ed. — Ch. 5 (transient response),
Ch. 5.8 (steady-state errors), Ch. 6 (root-locus analysis and design), Ch. 7 (frequency-response
methods, Bode diagrams, stability margins).
N. S. Nise, Control Systems Engineering, 8th ed. — Ch. 4 (time response),
Ch. 6 (Routh–Hurwitz stability), Ch. 7 (steady-state errors, system type),
Ch. 8 (root locus), Ch. 10 (gain and phase margins).
R. C. Dorf and R. H. Bishop, Modern Control Systems, 13th ed. — Ch. 5, 7, 9.
G. F. Franklin, J. D. Powell and A. Emami-Naeini, Feedback Control of Dynamic
Systems, 8th ed. — Ch. 3 (Laplace methods), Ch. 5 (root locus), Ch. 6 (frequency
response).
Question 6: Asymptotic Bode magnitude plot, phase margin and gain margin
Given. Loop gain function
$G(s)=20/[(s+5)(0.1s+1)(0.025s+1)]$, assumed to be closed with unity negative
feedback. Note that the first factor is not in time-constant form.
Find. (a) the asymptotic (straight-line) Bode magnitude plot; (b) the phase
margin and gain margin of the loop.
Figure 6.1 — Bode magnitude of
$G(s)$. The dashed red trace is the straight-line asymptotic construction
(12.04 dB flat, then $-20$, $-40$ and $-60$ dB/decade
after the corners at 5, 10 and 40 rad/s); the solid blue trace is the exact magnitude. Gain
crossover occurs at 11.5 rad/s.
Approach. Put every factor into time-constant form so the d.c. gain and
the corner frequencies can be read directly, build the asymptotic magnitude by accumulating
$-20$ dB/decade at each corner, then find the gain-crossover and
phase-crossover frequencies and evaluate the two margins there.
Convert to time-constant (Bode) form. The factor
$(s+5)$ must be written $5(0.2s+1)$:
$$G(s)=\frac{20}{5(0.2s+1)(0.1s+1)(0.025s+1)}
=\frac{4}{(0.2s+1)(0.1s+1)(0.025s+1)}$$
The d.c. gain is therefore $K=4$, i.e.
$20\log_{10}4 = 12.04\ \text{dB}$, not 26 dB as the un-normalised numerator
might suggest.
List the corner frequencies. Each factor
$(\tau s+1)$ corners at $\omega=1/\tau$:
$$\omega_1=\frac{1}{0.2}=5,\qquad
\omega_2=\frac{1}{0.1}=10,\qquad
\omega_3=\frac{1}{0.025}=40\ \ \text{rad/s}$$
All three are simple poles, so each adds $-20$ dB/decade.
Build the asymptotic magnitude. Starting flat at 12.04 dB (there is no
integrator, so the low-frequency asymptote is horizontal), the slope steepens at each corner:
Frequency band (rad/s)
Slope (dB/decade)
Asymptote value at the upper corner
$\omega \lt 5$
0
12.04 dB (flat)
$5 \lt \omega \lt 10$
$-20$
$12.04-20\log_{10}2=6.02$ dB at 10
$10 \lt \omega \lt 40$
$-40$
$6.02-40\log_{10}4=-18.06$ dB at 40
$\omega \gt 40$
$-60$
$-60$ dB at 200
The straight-line plot crosses 0 dB on the $-40$ dB/decade segment: solving
$6.02-40\log_{10}(\omega/10)=0$ gives $\omega\approx 14.1\ \text{rad/s}$.
The exact curve lies below the asymptotes near the corners (each simple pole costs 3 dB at
its own corner), so the true crossover is somewhat lower, as the next step shows.
Find the exact gain-crossover frequency. Setting
$$|G(j\omega)|=\frac{4}{\sqrt{1+(0.2\omega)^2}\sqrt{1+(0.1\omega)^2}\sqrt{1+(0.025\omega)^2}}=1$$
and solving numerically gives
$$\omega_{gc}=11.54\ \text{rad/s}$$
which is about 18 % below the straight-line estimate of 14.1 rad/s — the expected
penalty for reading a crossover off asymptotes alone when two corners are close together.
Evaluate the phase margin. The phase of three simple poles is
$$\angle G(j\omega)=-\left[\tan^{-1}(0.2\omega)+\tan^{-1}(0.1\omega)+\tan^{-1}(0.025\omega)\right]$$
At $\omega_{gc}=11.54\ \text{rad/s}$ the three terms are
$66.6^\circ$, $49.1^\circ$ and $16.1^\circ$, totalling
$131.8^\circ$, so
$$\text{PM}=180^\circ+\angle G(j\omega_{gc})=180^\circ-131.8^\circ
\quad\Longrightarrow\quad
\boxed{\;\text{PM}=48.2^\circ\;}$$
Find the phase-crossover frequency. Solving
$\angle G(j\omega)=-180^\circ$, i.e.
$\tan^{-1}(0.2\omega)+\tan^{-1}(0.1\omega)+\tan^{-1}(0.025\omega)=180^\circ$, gives
$$\omega_{pc}=25.50\ \text{rad/s}$$
which lies on the $-40$ dB/decade segment, well above gain crossover.
Evaluate the gain margin. The magnitude there is
$$|G(j\omega_{pc})|=0.2370
\quad\Longrightarrow\quad
\text{GM}=\frac{1}{0.2370}=4.22
\quad\Longrightarrow\quad
\boxed{\;\text{GM}=12.5\ \text{dB}\;}$$
Both margins are positive and the open loop has no right-half-plane poles, so the closed
loop is stable. As a check, multiplying the loop gain by exactly 4.22 (raising the numerator
from 20 to 84.4) places the closed-loop poles precisely on the imaginary axis, confirming the
gain margin independently.
Assess the design. A phase margin near
$48^\circ$ and a gain margin above 12 dB is a comfortable, well-damped design by
the usual industrial guideline (PM 45–60°, GM 6–12 dB). The closed-loop
response would show a modest overshoot — the rule of thumb
$\zeta\approx\text{PM}/100$ gives $\zeta\approx 0.48$, close to the
$\zeta=0.5$ targeted in Question 4 — and a bandwidth of roughly
$\omega_{gc}$, about 12 rad/s.