NivaarExam PrepOfficial exam papers ↗

22-Elec-B2 Advanced Control Systems · December 2017

Question 2 of 5: Root locus from a pole-zero map

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

Notes on this paper

Paper format. 16-Elec-B2 Advanced Control Systems, December 2017 — a three-hour open-book examination; any non-communicating calculator is permitted. The cover page states “Any four questions constitute a complete paper. Only the first four questions as they appear in your answer book will be marked” and “All questions are of equal value (25%)”. The sitting prints five questions, so each carries 25 marks and a candidate answers four. Question 1 is a seventeen-item multiple-choice block whose per-item weights are printed in the margin as [1], [2] or [3] and total exactly 25. A table of inverse Laplace transforms and a table of Laplace/z transforms are appended. All five questions are worked below, because this set is a study resource rather than a timed sitting.

Reference texts. N. S. Nise, Control Systems Engineering, 8th ed. (Ch. 4 time response, percent overshoot, peak time and settling time; Ch. 6 Routh–Hurwitz stability; Ch. 7 steady-state error and system type; Ch. 8 root-locus sketching rules including break-away/break-in points and angles of departure; Ch. 10 frequency response and stability margins); R. C. Dorf and R. H. Bishop, Modern Control Systems, 13th ed. (Ch. 5, 6, 7, 9); K. Ogata, Modern Control Engineering, 5th ed. (Ch. 5 transient response, Ch. 6 root-locus analysis, Ch. 7 frequency-response analysis); G. F. Franklin, J. D. Powell and A. Emami-Naeini, Feedback Control of Dynamic Systems, 8th ed. (Ch. 3, 5, 6). These are the references listed for exam code 16-Elec-B2 in the Engineers Canada syllabus, and this paper’s vocabulary (“break-in point”, “intercept of the asymptotes”, “angle of departure”) follows Nise closely.

Check — three source readings checked against the printed paper.

(1) Question 1, item 2. The printed sentence reads “where \(\theta(s)\), \(T(s)\), \(J\) and \(K\) represent torque, angular displacement, moment of inertia and spring constant, respectively”. The paper really does print the first two descriptions in the wrong order. The drawing beside it settles the physics — \(T(t)\) is the applied torque acting on the inertia \(J\), and \(\theta\) is the resulting angular displacement, which is also the only reading under which \(\theta(s)/T(s)\) has the units of a compliance. The answer does not depend on the naming at all, because only the denominator is used.

(2) Question 1, item 16. The margins are read off a printed Bode pair, so they are inherently approximate. The printed Bode pair gives a gain crossover near \(\omega \approx 3\) rad/s where the phase curve sits near \(-116^\circ\), and a phase crossover near \(\omega \approx 10\) rad/s where the magnitude curve is near \(-21\) dB. Both readings point at the same option. The graph is coarse enough that the honest statement is “approximately”, exactly as the question words it.

(3) Question 2. The pole-zero map prints its two crosses to the right of the imaginary axis, at \(s = 1 \pm j1\), and its two circles on the negative real axis at \(-2\) and \(-3\). This is the whole point of the question, because an open-loop plant with two right-half-plane poles is unstable until enough gain is applied, which is the opposite of the usual root-locus habit.

Question 2: Root locus from a pole-zero map (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 unity-feedback loop whose open-loop pole-zero map prints two crosses in the right half-plane and two circles on the negative real axis, together with the break-in coordinate for part (d).

Given data read from the printed s-plane map
QuantitySymbolValue
Open-loop poles (crosses)\(p_{1,2}\)\(1 + j1\) and \(1 - j1\)
Open-loop zeros (circles)\(z_{1,2}\)\(-2\) and \(-3\)
Feedback\(H(s)\)unity, negative
Break-in coordinate (given in part d)\(\sigma_b\)\(-2.43\)
-3-2-112-11σjω× open-loop poles   ○ open-loop zeross-plane
The printed open-loop pole-zero map: two right-half-plane poles at \(1 \pm j1\) and two real zeros at \(-2\) and \(-3\).

Find. (a) the open-loop transfer function, (b) the angle of departure from the complex poles, (c) the closed-loop transfer function and, from a Routh table, the imaginary-axis crossing and the gain that produces it, and (d) the completed root-locus sketch.

Approach. Write \(G(s)\) directly from the map, apply the angle condition at a point infinitesimally close to a pole to get the departure angle, form the characteristic polynomial and use the Routh array to locate the gain at which its \(s^{1}\) row vanishes, then assemble the sketch from the real-axis rule, the departure angles, the crossing and the given break-in.

  1. Part (a) — read the singularities off the map and assemble \(G(s)\). A cross marks a pole and a circle marks a zero, so the open-loop function has poles at \(s = 1 \pm j1\) and zeros at \(s = -2\) and \(s = -3\). Building the factored form and expanding, $$G(s) = \frac{K(s+2)(s+3)}{(s-1-j1)(s-1+j1)} = \frac{K(s+2)(s+3)}{(s-1)^{2}+1} = \boxed{\frac{K\left(s^{2}+5s+6\right)}{s^{2}-2s+2}}$$ Both open-loop poles lie in the right half-plane, so the plant itself is unstable; the loop will have to supply enough gain to pull them across.
  2. Part (b) — angles from each singularity to the departing pole. The angle condition states that every point on the locus satisfies \(\sum \angle(\text{zeros}) - \sum \angle(\text{poles}) = (2k+1)180^\circ\). Applied at a test point a hair away from \(p_1 = 1+j1\), every vector is measured from the singularity to \(p_1\): $$\angle(p_1+2) = \arctan\frac{1}{3} = 18.435^\circ, \qquad \angle(p_1+3) = \arctan\frac{1}{4} = 14.036^\circ$$ $$\angle(p_1-p_2) = \angle(j2) = 90^\circ$$
  3. Solve the angle condition for the unknown departure angle. Writing the condition with \(\theta_d\) as the contribution of the departing pole itself, $$\theta_d = 180^\circ + \sum \angle(\text{zeros}) - \sum \angle(\text{other poles}) = 180^\circ + 18.435^\circ + 14.036^\circ - 90^\circ$$ $$\boxed{\theta_{d} = +122.47^\circ \ \text{from } 1+j1, \qquad \theta_{d} = -122.47^\circ \ \text{from } 1-j1}$$ The conjugate pole departs at the mirror-image angle, as it must, because the locus is symmetric about the real axis. Taking \(k = \pm 1, \pm 2, \dots\) merely adds multiples of \(360^\circ\) and returns the same direction, which is what the note in the question is pointing at.
  4. Part (c) — form the closed-loop transfer function. For unity negative feedback, $$T(s) = \frac{G(s)}{1+G(s)} = \frac{K\left(s^{2}+5s+6\right)}{\left(s^{2}-2s+2\right)+K\left(s^{2}+5s+6\right)}$$ Collecting powers of \(s\) in the denominator gives the characteristic polynomial $$\boxed{(1+K)s^{2} + (5K-2)s + (2+6K) = 0}$$
  5. Build the Routh array and find where the \(s^{1}\) row vanishes. For a quadratic the array is short, and the first column is simply the three coefficients: $$\begin{array}{c|cc} s^{2} & 1+K & 2+6K \\ s^{1} & 5K-2 & 0 \\ s^{0} & 2+6K & \end{array}$$ An entire row of zeros — the signature of a pair of roots on the imaginary axis — occurs when \(5K-2 = 0\), i.e. $$\boxed{K = 0.4}$$ Positive \(K\) below 0.4 makes the \(s^{1}\) entry negative, so the loop is unstable there: it is the lower bound on gain, not an upper one, because the open-loop poles start in the right half-plane.
  6. Extract the crossing frequency from the auxiliary polynomial. The row above the vanishing row forms the auxiliary polynomial, whose roots are the imaginary-axis pair. With \(K = 0.4\), $$(1+K)s^{2} + (2+6K) = 1.4\,s^{2} + 4.4 = 0 \;\Longrightarrow\; s^{2} = -\frac{4.4}{1.4} = -3.1429$$ $$\boxed{s = \pm j1.7728 \ \text{at } K = 0.4}$$ Substituting \(K = 0.4\) back into the quadratic and solving numerically returns \(s = \pm j1.7728\) exactly, with a real part of zero to machine precision — the cheapest possible check on the Routh work.
  7. Part (d) — assemble the sketch. Four construction facts fix the whole picture. There are \(n = m = 2\), so no branch runs to infinity and there are no asymptotes; every branch starts on a pole and ends on a zero. On the real axis the locus exists only where an odd number of real singularities lie to the right, which is the segment \(-3 \le \sigma \le -2\) between the two zeros. The branches leave \(1 \pm j1\) at \(\pm 122.47^\circ\), sweep leftward across the imaginary axis at \(\pm j1.7728\) when \(K = 0.4\), and meet the real axis at the given break-in point \(\sigma_b = -2.43\) — recomputing it from \(dK/ds = 0\) with \(K = -(s^{2}-2s+2)/[(s+2)(s+3)]\) gives \(7s^{2}+8s-22 = 0\), whose admissible root is \(-2.4341\), confirming the printed value. The gain there is \(K = 52.08\); beyond it the two roots separate along the real axis and run into the zeros at \(-2\) and \(-3\).
  8. State the stability conclusion the sketch implies. Because both branches start in the right half-plane and cross into the left at \(K = 0.4\), and no branch ever returns, $$\boxed{\text{the closed loop is stable for } K \gt 0.4}$$ This is the mark-earning observation: for this plant more gain is stabilising, which reverses the habit built on minimum-phase plants where increasing gain eventually destabilises the loop.
-3-2-112-2-112Re(s)Im(s)jω crossing ±j1.773 (K = 0.4)break-in −2.434 (K = 52.08)122.5°departure ±122.5°; stable for K > 0.4
Question 2(d): the completed root locus. Branches depart the right-half-plane poles at \(\pm 122.5^\circ\), cross the imaginary axis at \(\pm j1.773\) when \(K = 0.4\), break in on the real axis at \(-2.434\) when \(K = 52.08\), and terminate on the zeros at \(-2\) and \(-3\).
Question 2 — results
PartQuantityResult
(a)Open-loop poles\(s = 1 \pm j1\) (both unstable)
(a)Open-loop zeros\(s = -2\), \(s = -3\)
(a)Transfer function\(G(s) = K(s^{2}+5s+6)/(s^{2}-2s+2)\)
(b)Angle of departure\(\pm 122.47^\circ\)
(c)Characteristic polynomial\((1+K)s^{2}+(5K-2)s+(2+6K)\)
(c)\(j\omega\)-axis crossing\(s = \pm j1.7728\)
(c)Gain at the crossing\(K = 0.4\)
(d)Break-in point (and its gain)\(\sigma = -2.434\), \(K = 52.08\)
(d)Stability range\(K \gt 0.4\)