22-Elec-B2 Advanced Control Systems · December 2013
Question 6 of 6: Proportional design to a specified gain margin for a plant with transport delay
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
Paper format. National Examinations, 07-Elec-B2 Advanced Control
Systems, December 2013 — 3 hours, closed book. Six questions; the rubric states that
“any four questions constitute a complete paper” and that “all questions are of
equal value”, so each is worth 25 marks. Tables of Laplace and z-transforms are appended to
the paper. Only a Casio FX-991 or Sharp EL-540 calculator is permitted.
Reference texts. G. F. Franklin, J. D. Powell and A. Emami-Naeini,
Feedback Control of Dynamic Systems, 7th ed., Pearson (Routh–Hurwitz, stability
margins, steady-state error, frequency response); K. J. Åström and R. M. Murray,
Feedback Systems: An Introduction for Scientists and Engineers, 2nd ed., Princeton
(sensitivity function, loops with transport delay, non-minimum-phase limitations); K. Ogata,
Modern Control Engineering, 5th ed., Pearson (state-space realisations, controllability
and observability, pole placement); G. F. Franklin, J. D. Powell and M. L. Workman,
Digital Control of Dynamic Systems, 3rd ed., Ellis-Kagle (zero-order-hold equivalents,
deadbeat design, discrete identification). These are the works listed by Engineers Canada / EGBC
for the Elec-B2 syllabus.
Scope of this solution. Although the paper is marked on four questions, all
six are worked here so that the set serves as a complete study resource. Every boxed number is
recomputed from the question’s own data for this paper.
Check: two printing defects in the source paper. In Question 3 the leading
entry of the second row of the system matrix is missing from the print — the row reads
“0 −1” with the first column blank — and the
equation is typeset as “$\dots x + B = [0\ 1\ 1]^{\mathsf T}u$” rather than
$\dot x = Ax + Bu$. Reading the column alignment against rows 1 and 3, the missing entry is
$a_{21}=0$; that reading is adopted throughout and makes the system controllable, so the design
asked for is well posed. In Question 4(b) the paper asks for “$T(z)$ that relates $X(z)$ to
$R(z)$”, but no signal $X$ is defined anywhere in the question or its figure; $X(z)$ is read
as $Y(z)$.
Question 6: Proportional design to a specified gain margin for a plant with transport delay (25 marks)
Given. A type-1 plant with one lag and a pure transport delay, to be closed
with a proportional gain against a gain-margin specification.
Given data — Question 6
Quantity
Value
Plant
$P(s) = \dfrac{e^{-0.2s}}{s(2s+1)}$
Transport delay
$\theta = 0.2$ s
Lag time constant
$T = 2$ s (pole at $s=-0.5$)
Controller
$C(s) = K$, proportional only
Specification
gain margin $= 8$ dB
Feedback
unity negative
Find. (a) the gain $K$ meeting the 8 dB gain margin; (b) the phase margin that
results; (c) the steady-state output to a unit step set-point.
Approach. Locate the phase crossover from the phase equation alone — the
delay contributes phase but no attenuation, so the crossover frequency does not depend on $K$
— then scale the magnitude there to the required margin, and finally find the new gain
crossover to read the phase margin.
Write the frequency response. A pure delay has unit magnitude at every
frequency and a phase that grows linearly with frequency, so
$$|K P(j\omega)| = \frac{K}{\omega\sqrt{1+4\omega^{2}}},\qquad
\angle KP(j\omega) = -\,\theta\omega - 90^\circ - \arctan(2\omega),$$
with $\theta\omega$ expressed in radians before conversion. Because the delay is transcendental,
the closed-loop characteristic equation $s(2s+1)+Ke^{-0.2s}=0$ is not a polynomial and the
Routh test is unavailable — the frequency-domain route is the only closed-form
one.
Locate the phase crossover. Setting the phase to $-180^\circ$ and working in
radians,
$$0.2\,\omega + \arctan(2\omega) = \frac{\pi}{2}.$$
The left side rises monotonically from 0, so the root is unique; solving numerically,
$$\omega_{c\varphi} = 1.5553\ \text{rad/s}.$$
Crucially this frequency is set by the phase alone and is therefore independent of
$K$ — a proportional gain slides the magnitude curve up and down but never moves the
$-180^\circ$ crossing.
Evaluate the plant magnitude there.
$$|P(j\omega_{c\varphi})| = \frac{1}{1.5553\sqrt{1+4(1.5553)^{2}}}
= \frac{1}{1.5553 \times 3.2674} = 0.19679,$$
so the uncompensated loop ($K=1$) already has a gain margin of
$-20\log_{10}(0.19679) = 14.12$ dB. Meeting an 8 dB specification therefore means
increasing the gain.
Solve for the required gain. A gain margin of 8 dB means the loop magnitude at
the phase crossover must equal $10^{-8/20}$:
$$K\,|P(j\omega_{c\varphi})| = 10^{-8/20} = 0.39811
\quad\Longrightarrow\quad K = \frac{0.39811}{0.19679}$$
$$\boxed{K = 2.0230}$$
Checking, $20\log_{10}(2.0230 \times 0.19679) = -8.00$ dB exactly.
Find the new gain crossover. Setting $|KP(j\omega)|=1$,
$$\omega\sqrt{1+4\omega^{2}} = K = 2.0230 \quad\Longrightarrow\quad
\omega_{cg} = 0.9456\ \text{rad/s}.$$
Read the phase margin. Evaluating the phase at that frequency, in degrees,
$$\angle KP(j\omega_{cg}) = -\frac{180}{\pi}(0.2)(0.9456) - 90^\circ - \arctan\!\big(2(0.9456)\big)
= -10.84^\circ - 90^\circ - 62.13^\circ = -162.97^\circ,$$
$$\boxed{\text{PM} = 180^\circ - 162.97^\circ = 17.0^\circ \ \text{at } \omega_{cg}=0.946\
\text{rad/s}}$$
The delay costs almost 11° of that phase on its own. A 17° margin is lightly damped
— roughly $\zeta \approx 0.17$ by the usual $\zeta \approx \text{PM}/100$ rule of thumb
— which is the price of insisting on an 8 dB rather than a larger gain margin on a
delayed plant.
Determine the steady-state output. The plant contains a free integrator, so
the loop is type 1 and the position error constant is infinite:
$$K_{pos} = \lim_{s\to0} KP(s) = \lim_{s\to0}\frac{K e^{-0.2s}}{s(2s+1)} = \infty
\quad\Longrightarrow\quad e_{ss} = \frac{1}{1+K_{pos}} = 0 .$$
Therefore the output tracks the unit step exactly in steady state:
$$\boxed{y(\infty) = 1, \qquad e_{ss} = 0}$$
This holds for any stable $K$, so it is a property of the plant structure rather than of the
design just completed.
Figure Q6.1 — Bode plot of the compensated loop. The delay contributes phase only, so the −180° crossing sits at ω = 1.555 rad/s regardless of K; choosing K = 2.0230 pulls the magnitude there down to −8 dB, which is the required gain margin.
Check: the delay is handled exactly, not by approximation. All numbers above
use $e^{-j0.2\omega}$ directly. A first-order Padé substitution
$e^{-0.2s}\approx(1-0.1s)/(1+0.1s)$ would make the characteristic equation a polynomial and allow
a Routh solution, but it shifts the phase crossover by a few per cent and with it the design gain;
the exact frequency-domain calculation is preferable and is what the mark scheme expects here. The
$17^\circ$ phase margin also means the closed-loop step response will be markedly oscillatory: if
the specification really is on transient quality rather than on gain margin, this plant wants a
lead or a PI-lead compensator rather than a bare proportional gain.