NivaarExam PrepOfficial exam papers ↗

22-Elec-B2 Advanced Control Systems · May 2015

Question 4 of 6: Pole placement by state feedback with reference scaling

Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)

Notes on this paper

Paper format. National Exams, May 2015, 07-Elec-B2 Advanced Control Systems — three hours, closed book. Six questions are set; the rubric states that “any four questions constitute a complete paper” and that “all questions are of equal value”, so each carries 25 marks. Tables of inverse Laplace and inverse z-transforms are appended as pages 4 and 5, and only an approved Casio or Sharp calculator is permitted. All six questions are solved here, because this set is a study resource rather than a timed sitting. The paper deliberately mixes continuous-time frequency-response work (Questions 1 and 6), state-space methods (Questions 2 and 4) and sampled-data design (Questions 3 and 5).

Reference texts. G. F. Franklin, J. D. Powell and A. Emami-Naeini, Feedback Control of Dynamic Systems, 7th ed., Pearson (frequency response, stability margins, steady-state error, lead compensation); K. J. Åström and R. M. Murray, Feedback Systems: An Introduction for Scientists and Engineers, 2nd ed., Princeton (sensitivity functions and loops containing transport delay); K. Ogata, Modern Control Engineering, 5th ed., Pearson (state-space realisations, controllability, observability and pole placement); G. F. Franklin, J. D. Powell and M. L. Workman, Digital Control of Dynamic Systems, 3rd ed., Addison-Wesley (zero-order-hold equivalents, the Jury test, discrete steady-state error); L. Ljung, System Identification: Theory for the User, 2nd ed., Prentice Hall (least-squares estimation of difference-equation models). These are the works listed by Engineers Canada and EGBC for the Elec-B2 syllabus.

Two readings matter for the marks: in Question 1 the disturbance d enters the plant-input summing junction through a minus sign while u enters through a plus, so the plant sees $u - d$; and the Question 3 table is $y = \{0,\,10,\,4,\,3,\,1.6,\,0.4,\,0.3\}$ against $u = \{1,\,0,\,0,\,0,\,0,\,0,\,0\}$.

Question 4: Pole placement by state feedback with reference scaling (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 third-order single-input single-output state model and three required closed-loop pole locations.

QuantityValue
State matrix$A = \begin{pmatrix}-1 & -2 & -2\\ 0 & -1 & 1\\ 1 & 0 & -1\end{pmatrix}$
Input matrix$B = \begin{pmatrix}2 & 0 & 1\end{pmatrix}^{T}$
Output matrix$C = \begin{pmatrix}1 & 0 & 1\end{pmatrix}$, $D = 0$
Control law$u = L\,r - K x$, with $K = \begin{pmatrix}k_1 & k_2 & k_3\end{pmatrix}$
Required closed-loop poles$s = -10$, $s = -3 \pm j4$
Steady-state requirement$e = r - y \to 0$ for a step reference

Find. The feedback row vector $K$ that places the three poles, and the scalar reference gain $L$ that removes the steady-state tracking error.

-10-8-6-4-20-4-2024Real axis (1/s)Imaginary axis (1/s)grey x = open loop red x = placed
Question 4: the three open-loop poles (grey) and the three placed closed-loop poles (red) at s = -10 and s = -3 +/- j4.

Approach. Confirm the pair $(A,B)$ is controllable, expand the closed-loop characteristic polynomial symbolically in the three unknown gains, match it term by term against the desired polynomial, and finally choose $L$ so that the closed-loop DC gain from $r$ to $y$ is exactly one.

  1. Verify that arbitrary pole placement is possible. Successive multiplications give $$AB = \begin{pmatrix}-4\\ 1\\ 1\end{pmatrix}, \qquad A^2B = \begin{pmatrix}0\\ 0\\ -5\end{pmatrix}, \qquad \mathcal{C} = \begin{pmatrix}2 & -4 & 0\\ 0 & 1 & 0\\ 1 & 1 & -5\end{pmatrix},$$ and expanding along the second row, $\det\mathcal{C} = -10 \neq 0$. The pair is controllable, so any self-conjugate set of three poles is reachable by state feedback and the design is guaranteed to have a solution.
  2. Write the desired characteristic polynomial. Multiplying out the required factors, $$\Delta_{d}(s) = (s+10)\bigl[(s+3)^2 + 16\bigr] = (s+10)(s^2 + 6s + 25),$$ $$\Delta_{d}(s) = s^3 + 16s^2 + 85s + 250 .$$
  3. Expand the closed-loop polynomial in the unknown gains. With $u = -Kx$ the state matrix becomes $A - BK$, and expanding $\det\bigl(sI - A + BK\bigr)$ symbolically gives $$\Delta_{cl}(s) = s^3 + \bigl(2k_1 + k_3 + 3\bigr)s^2 + \bigl(2k_1 + k_2 + 4k_3 + 5\bigr)s + \bigl(-2k_1 + 3k_2 + 3k_3 + 5\bigr).$$ Note that the open-loop polynomial $s^3 + 3s^2 + 5s + 5$ is recovered when $k_1 = k_2 = k_3 = 0$, which is a useful check on the expansion before going further.
  4. Match the coefficients. Equating like powers of $s$ produces three linear equations in three unknowns: $$\begin{aligned} s^2: &\quad 2k_1 + k_3 + 3 = 16\\ s^1: &\quad 2k_1 + k_2 + 4k_3 + 5 = 85\\ s^0: &\quad -2k_1 + 3k_2 + 3k_3 + 5 = 250 \end{aligned}$$
  5. Solve the linear system. From the first equation $k_3 = 13 - 2k_1$; substituting into the second gives $k_2 = 80 - 2k_1 - 4k_3 = 28 + 6k_1$, and the third then reduces to $-2k_1 + 3(28 + 6k_1) + 3(13 - 2k_1) = 245$, i.e. $10k_1 = 122$. Hence $$\boxed{\;K = \begin{pmatrix}k_1 & k_2 & k_3\end{pmatrix} = \begin{pmatrix}12.2 & 101.2 & -11.4\end{pmatrix} = \tfrac{1}{5}\begin{pmatrix}61 & 506 & -57\end{pmatrix}\;}$$ The third gain is negative, which is perfectly admissible — state feedback gains carry no sign restriction, and here the third state must be fed back positively to pull the lightly damped open-loop pair leftwards.
  6. Confirm the placement. Substituting $K$ back and factoring, $$\det\bigl(sI - A + BK\bigr) = s^3 + 16s^2 + 85s + 250 = (s+10)(s^2 + 6s + 25),$$ whose roots are $s = -10$ and $s = -3 \pm j4$ as required. The dominant pair has $\zeta = 3/5 = 0.6$ and $\omega_n = 5\ \text{rad/s}$, with the real pole placed twice as far out so that the pair genuinely dominates.
  7. Set up the reference gain. The feedback gain $K$ fixes the poles but also changes the DC gain, so a step reference would otherwise be tracked with an offset. The closed-loop system is $$\dot x = (A - BK)x + BL\,r, \qquad y = Cx,$$ and at steady state $\dot x = 0$ gives $x_{ss} = -(A-BK)^{-1}BL\,r$, so $$y_{ss} = -C\,(A-BK)^{-1}B\,L\,r \stackrel{!}{=} r \;\Longrightarrow\; L = \frac{1}{-C\,(A-BK)^{-1}B}.$$
  8. Evaluate the reference gain. Inverting the closed-loop matrix and forming the scalar product gives $-C(A-BK)^{-1}B = 1/250$, so $$\boxed{\;L = 250\;}$$ and the complete control law is $u(t) = 250\,r(t) - \begin{pmatrix}12.2 & 101.2 & -11.4\end{pmatrix}x(t)$.
  9. Check the finished design. The closed-loop transfer function is $$\frac{Y(s)}{R(s)} = C\bigl(sI - A + BK\bigr)^{-1}BL = \frac{250\,(3s^2 + 6s + 1)}{s^3 + 16s^2 + 85s + 250},$$ whose value at $s = 0$ is $250/250 = 1$. The step response therefore settles on the reference exactly, and $e = r - y \to 0$ as required.

Check: the zero steady-state error here is a calibration, not an integrator. The gain $L = 250$ cancels the DC error of the nominal model only. Any drift in $A$, $B$ or $C$ reintroduces an offset, because there is no integral action to absorb it, and a load disturbance entering at the plant input is not rejected at all. If robust zero error were required the design would have to be extended with an integrator state and the poles re-placed for the augmented fourth-order system. The question as set asks only for the nominal condition, which is what is delivered above.

Question 4 — results
QuantitySymbolValue
Controllability determinant$\det\mathcal{C}$$-10$ (controllable)
Open-loop characteristic polynomial$\Delta_{ol}$$s^3 + 3s^2 + 5s + 5$
Desired characteristic polynomial$\Delta_{d}$$s^3 + 16s^2 + 85s + 250$
Feedback gains$K$$\begin{pmatrix}12.2 & 101.2 & -11.4\end{pmatrix}$
Reference gain$L$250
Dominant pair damping and frequency$\zeta,\ \omega_n$0.6, 5 rad/s
Closed-loop DC gain$Y(0)/R(0)$1 (zero step error)