Skip to main content

Transfer function from third-order state model

Short answer

The transfer function is Y(s)/U(s) = 10(s^2+3s+2)/(s^3+3s^2+2s+1), found from C adj(sI-A)B / det(sI-A).

For given state equation and output equation, find transfer function Y(s)/U(s). Ẋ = [0 1 0; 0 0 1; −1 −2 −3]X + [10; 0; 0]u and y = [1 0 0]X

For x˙=Ax+Bu\dot{x} = Ax + Bu, y=Cx+Duy = Cx + Du:

Y(s)U(s)=C(sI−A)−1B+D=C adj⁡(sI−A) B∣sI−A∣+D\frac{Y(s)}{U(s)} = C(sI - A)^{-1}B + D = \frac{C\,\operatorname{adj}(sI-A)\,B}{|sI-A|} + D A=[010001−1−2−3],B=[1000],C=[100],D=0A = \begin{bmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -1 & -2 & -3 \end{bmatrix},\quad B = \begin{bmatrix} 10 \\ 0 \\ 0 \end{bmatrix},\quad C = \begin{bmatrix} 1 & 0 & 0 \end{bmatrix},\quad D = 0

Step 1: sI−AsI - A and its determinant

sI−A=[s−100s−112s+3]sI - A = \begin{bmatrix} s & -1 & 0 \\ 0 & s & -1 \\ 1 & 2 & s+3 \end{bmatrix} ∣sI−A∣=s[s(s+3)+2]−(−1)[0⋅(s+3)−(−1)(1)]+0=s3+3s2+2s+1\begin{aligned} |sI - A| &= s\left[s(s+3) + 2\right] - (-1)\left[0\cdot(s+3) - (-1)(1)\right] + 0 \\ &= s^3 + 3s^2 + 2s + 1 \end{aligned}

Step 2: Required part of the adjoint

Because BB has only its first element non-zero and CC picks only the first row, only the (1,1) element of the adjoint is needed:

adj⁡(sI−A)11=∣s−12s+3∣=s2+3s+2\operatorname{adj}(sI-A)_{11} = \begin{vmatrix} s & -1 \\ 2 & s+3 \end{vmatrix} = s^2 + 3s + 2

Full first column of the adjoint: [s2+3s+2−1−s]\begin{bmatrix} s^2+3s+2 \\ -1 \\ -s \end{bmatrix}, so

adj⁡(sI−A)B=10[s2+3s+2−1−s]\operatorname{adj}(sI-A)B = 10\begin{bmatrix} s^2+3s+2 \\ -1 \\ -s \end{bmatrix}

Step 3: Transfer function

C adj⁡(sI−A)B=10(s2+3s+2)C\,\operatorname{adj}(sI-A)B = 10(s^2 + 3s + 2) Y(s)U(s)=10(s2+3s+2)s3+3s2+2s+1=10(s+1)(s+2)s3+3s2+2s+1\frac{Y(s)}{U(s)} = \frac{10(s^2+3s+2)}{s^3+3s^2+2s+1} = \frac{10(s+1)(s+2)}{s^3+3s^2+2s+1}

Stability check

Routh array for s3+3s2+2s+1s^3 + 3s^2 + 2s + 1:

Row
s3s^312
s2s^231
s1s^13×2−13=1.67\frac{3\times2 - 1}{3} = 1.670
s0s^01

No sign change, so the system is stable. The poles are −2.325-2.325 and −0.338±j0.562-0.338 \pm j0.562; the zeros are −1-1 and −2-2.

Answer: Y(s)U(s)=10(s2+3s+2)s3+3s2+2s+1\dfrac{Y(s)}{U(s)} = \dfrac{10(s^2+3s+2)}{s^3+3s^2+2s+1}.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Control Systems chapter 8 questions →