NivaarExam PrepOfficial exam papers ↗

22-Elec-B2 Advanced Control Systems · December 2013

Question 4 of 6: Deadbeat design of a digital PI controller for an integrating plant

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 4: Deadbeat design of a digital PI controller for an integrating plant (25 marks)

Question text not reproduced: the examination questions are © Engineers and Geoscientists BC. Open the official past paper (linked at the top of this page) to read the question, then follow the worked solution below.

Given. A digital PI controller driving a pure integrator through a zero-order hold, with everything sampled at the same period.

Given data — Question 4
QuantityValue
Controller$C(z) = K_{1} + \dfrac{K_{2}}{1-z^{-1}}$ (proportional plus discrete integral)
Continuous plant$P(s) = 1/s$
Holdzero-order, output constant over each period
Sample period$h$ (symbolic; $h=0.5$ s is used for the plotted response)
Requirementall closed-loop poles at $z=0$ (deadbeat)

Find. (a) $K_{1}$ and $K_{2}$ as functions of $h$; (b) the closed-loop transfer function $T(z)=Y(z)/R(z)$; (c) the unit-step response of $y(t)$, including what happens between the sampling instants.

rh+−C(z)uhZOHP(s)ysampled feedback path
Figure Q4.1 — the sampled-data loop. Everything inside the digital part advances once per period h; the ZOH holds u constant between samples, so the continuous integrator sees a staircase input.

Approach. Replace the hold-plus-plant combination by its exact discrete equivalent, form the closed-loop characteristic polynomial in $z$, force it to equal $z^{2}$, and then reconstruct the continuous output by integrating the resulting staircase of control values.

  1. Discretise the hold and plant together. The zero-order-hold equivalent is $$G(z) = \left(1-z^{-1}\right)\mathcal{Z}\!\left\{\frac{P(s)}{s}\right\} = \left(1-z^{-1}\right)\mathcal{Z}\!\left\{\frac{1}{s^{2}}\right\} = \frac{z-1}{z}\cdot\frac{hz}{(z-1)^{2}} = \frac{h}{z-1},$$ using the $t \leftrightarrow hz/(z-1)^{2}$ entry of the supplied transform table. Physically this says the obvious thing: an integrator fed a constant $u$ for one period gains exactly $hu$.
  2. Put the controller over a common denominator. $$C(z) = K_{1} + \frac{K_{2}}{1-z^{-1}} = K_{1} + \frac{K_{2}z}{z-1} = \frac{(K_{1}+K_{2})z - K_{1}}{z-1}.$$ The controller pole at $z=1$ is the discrete integrator, which is what will guarantee zero steady-state step error.
  3. Form the closed-loop characteristic polynomial. The open loop is $L(z)=C(z)G(z)$, so $1+L(z)=0$ gives $$(z-1)^{2} + h\big[(K_{1}+K_{2})z - K_{1}\big] = 0,$$ which expands to $$z^{2} + \big[h(K_{1}+K_{2}) - 2\big] z + \big(1 - hK_{1}\big) = 0 .$$
  4. Impose the deadbeat condition. Both poles at the origin means the characteristic polynomial must be exactly $z^{2}$, so both remaining coefficients vanish: $$1 - hK_{1} = 0 \;\Rightarrow\; K_{1} = \frac{1}{h}; \qquad h(K_{1}+K_{2}) - 2 = 0 \;\Rightarrow\; K_{1}+K_{2} = \frac{2}{h} \;\Rightarrow\; K_{2} = \frac{1}{h}.$$ $$\boxed{K_{1}=K_{2}=\frac{1}{h}, \qquad C(z) = \frac{1}{h}\left[1 + \frac{z}{z-1}\right] = \frac{1}{h}\cdot\frac{2z-1}{z-1}}$$ At $h = 0.5$ s this is $K_{1}=K_{2}=2$. The gains scale as $1/h$: a faster sampler needs proportionally more gain to move the integrator the same distance in one period.
  5. Assemble the closed-loop transfer function. With the denominator now equal to $z^{2}$, $$T(z) = \frac{C(z)G(z)}{1+C(z)G(z)} = \frac{h\big[(K_{1}+K_{2})z-K_{1}\big]}{z^{2}} = \frac{h\left[\tfrac{2}{h}z - \tfrac{1}{h}\right]}{z^{2}} \quad\Longrightarrow\quad \boxed{T(z) = \frac{2z-1}{z^{2}} = 2z^{-1}-z^{-2}}$$ Note that $T(z)$ is independent of $h$: the sample period is absorbed entirely into the controller gains. Checking the DC gain, $T(1) = (2-1)/1 = 1$, so the steady-state step error is zero as the integrator promised.
  6. Evaluate the sampled step response. Because $T(z)=2z^{-1}-z^{-2}$ is a two-term finite impulse response, the unit-step response is just the running sum of its coefficients $\{0,\,2,\,-1\}$: $$y(0)=0,\quad y(h)=2,\quad y(2h)=2-1=1,\quad y(kh)=1 \ \ \text{for all } k\ge 2 .$$ The output reaches its final value in two steps and stays there exactly — the defining property of a deadbeat design, and the reason the poles were placed at the origin.
  7. Recover the inter-sample behaviour. The control sequence follows from $U(z)/R(z) = C(z)\big[1-T(z)\big]$ with $R(z)=z/(z-1)$, which reduces to $U(z)=\dfrac{1}{h}\left(2 - z^{-1}\right)$, that is $$u(0)=\frac{2}{h},\qquad u(h)=-\frac{1}{h},\qquad u(kh)=0 \ \ (k\ge2).$$ Since the hold keeps each value constant for a full period and the plant integrates, $y(t)$ is piecewise linear, not a staircase: it ramps up at slope $2/h$ from $0$ to $2$ over $0\le t\lt h$, ramps down at slope $-1/h$ from $2$ to $1$ over $h\le t\lt 2h$, and is then exactly constant at 1. At $h = 0.5$ s the two control values are $+4$ and $-2$. $$\boxed{y(t)\ \text{rises linearly to } 2 \text{ at } t=h,\ \text{falls linearly to } 1 \text{ at } t=2h,\ \text{and stays at } 1}$$
0.00.51.01.52.02.50.000.501.001.502.00output y(t)y(t) ramps in a straight line between samples; dots mark y(kh)0.00.51.01.52.02.5-2.500.002.505.00time t (s)control u(t)u(t) is held constant by the ZOH over each period
Figure Q4.2 — unit-step response with h = 0.5 s. Between samples y(t) is a straight ramp because P(s) = 1/s integrates the held control: it climbs to 2 at t = h, falls back to 1 at t = 2h, then stays flat. The sampled sequence is deadbeat (final value reached at k = 2) but the continuous output overshoots to 2.0 in between — that overshoot is the price of placing both poles at z = 0.

Check: the sampled response hides a 100 % overshoot. The sampled sequence $\{0,2,1,1,\dots\}$ settles in two steps, but the continuous output passes through 2 on the way — twice the set-point. Any specification written on $y(t)$ rather than on $y(kh)$ would fail this design. This is the standard objection to deadbeat control: it also demands large control amplitudes ($2/h$, unbounded as the sampler is speeded up) and it is exquisitely sensitive to plant gain error, because a pole at $z=0$ has no stability margin to lose. A more conservative pole placement inside the unit circle trades finite settling for a gentler response. The question also asks for $T(z)$ relating $X(z)$ to $R(z)$; $X$ is not defined anywhere in the paper and is read here as the output $Y$.

Final results — Question 4
QuantityResult
ZOH equivalent of the plant$G(z) = h/(z-1)$
Characteristic polynomial$z^{2}+[h(K_{1}+K_{2})-2]z+(1-hK_{1})$
Deadbeat gains$K_{1}=K_{2}=1/h$  (both $=2$ at $h=0.5$ s)
Controller$C(z) = \dfrac{2z-1}{h(z-1)}$
Closed-loop transfer function$T(z) = \dfrac{2z-1}{z^{2}}$, $T(1)=1$
Sampled step response$y(kh) = 0,\ 2,\ 1,\ 1,\ 1,\dots$
Control sequence$u(kh) = 2/h,\ -1/h,\ 0,\ 0,\dots$
Continuous outputpiecewise linear; peak 2 at $t=h$, constant 1 for $t\ge 2h$