NivaarExam PrepOfficial exam papers ↗

20-Bio-A2 Process Dynamics and Control · May 2017

Question 2 of 8: State-Space Model — Transfer Function and Step Response

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

Notes on this paper

National Exams / EGBC — May 2017 — 04-Bio-A2 Process Dynamics and Control. Three-hour open-book examination; any non-communicating calculator is permitted. The cover page states that only the first five questions in the answer book are marked, yet all eight problems are printed on the paper — all eight are solved below for completeness, since the full set is a study resource. Content spans sampled-data (digital) control, state-space-to-transfer-function modelling, Internal Model Control (IMC) of a non-minimum-phase dead-time process, Nyquist and Routh stability, frequency response (Bode/gain-margin) design and nonlinear-reactor linearisation.

Reference texts: D. E. Seborg, T. F. Edgar, D. A. Mellichamp & F. J. Doyle III, Process Dynamics and Control (4th ed., Wiley) — Laplace/z-domain modelling, transfer functions, Routh and Jury stability, Nyquist/Bode frequency response and IMC design; G. Stephanopoulos, Chemical Process Control: An Introduction to Theory and Practice (Prentice Hall) — Nyquist criterion, dead-time systems and sampled-data control. Standard control conventions (deviation variables; unity sensor/valve gains unless stated) are used throughout.

Problem 2: State-Space Model — Transfer Function and Step Response (20%)

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. Linear state-space model $\dot x_1=-2.4048x_1+7u$, $\dot x_2=0.8333x_1-2.2381x_2-1.117u$, $y=x_2$, zero initial conditions (deviation variables).

Find. (a) $Y(s)/U(s)$; (b) $y(t)$ for a unit step in $u$.

Approach. Laplace-transform each state equation from rest, eliminate $X_1(s)$ to get $Y/U$, then invert $Y(s)=[Y(s)/U(s)]\cdot(1/s)$ by partial fractions.

  1. (a) Eliminate $X_1$. $sX_1=-2.4048X_1+7U\Rightarrow X_1=\dfrac{7U}{s+2.4048}$. Substituting into $sX_2=0.8333X_1-2.2381X_2-1.117U$ and solving for $X_2=Y$: $$\boxed{\frac{Y(s)}{U(s)}=\frac{5.8331/(s+2.4048)-1.117}{s+2.2381}=\frac{3.1469-1.117\,s}{(s+2.2381)(s+2.4048)}.}$$ The numerator has a zero at $s=+2.818$ — a right-half-plane (RHP) zero, since $3.1469/1.117=2.818>0$.
  2. Recognise the inverse-response signature. An RHP zero makes $y(t)$ initially move opposite to its eventual direction before turning around — a classic "wrong-way" response seen in boiler drum level and other biomedical/process loops with competing fast and slow paths.
  3. (b) Step response by partial fractions. With $U(s)=1/s$, $Y(s)=\dfrac{3.1469-1.117s}{s(s+2.2381)(s+2.4048)}=\dfrac{A}{s}+\dfrac{B}{s+2.2381}+\dfrac{C}{s+2.4048}$. Evaluating each residue at its own pole gives $A=0.5847$, $B=-15.135$, $C=14.551$, so $$\boxed{y(t)=0.5847-15.135\,e^{-2.2381t}+14.551\,e^{-2.4048t}.}$$
  4. Check limits. $y(0)=0.5847-15.135+14.551\approx0$ (correct start), and $y(\infty)=A=0.5847$ (matches the final-value $Y(s)/U(s)$ evaluated at $s=0$). At $t=0.1\ \text{s}$, $y\approx-0.075$ — the response briefly dips negative before rising to $0.5847$, confirming the inverse response predicted by the RHP zero.
ResultValue
Transfer function$Y/U=(3.1469-1.117s)/[(s+2.2381)(s+2.4048)]$
RHP zero$s=+2.818$ (inverse response)
Step response$y(t)=0.5847-15.135e^{-2.2381t}+14.551e^{-2.4048t}$
Steady-state gain$y(\infty)=0.5847$