Skip to main content

State space model of circuit with current source

Short answer

Choose states: vC1, iL, vC2=vo. Then state equation: dot{x} = A x + B u, with A = [[-1,-4,0],[0.5,0,-0.5],[0,2,-2]], B = [[1,0],[0,0],[0,2]], C = [0,0,1], D=0. Inputs: Vi and Is.

Find state space representation of the system. [Figure: voltage source Vi in series with 4 Ω to node Vc1; 0.25 F capacitor from Vc1 to ground; 2 H inductor (current iL) from Vc1 to node Vc2; 0.5 F capacitor (voltage Vo) from Vc2 to ground; 1 Ω resistor from Vc2 to ground; current source is injecting into node Vc2 from ground]

For an electric circuit the natural state variables are the capacitor voltages and inductor currents, since they describe the stored energy. Here the states are vC1v_{C1}, iLi_L and vC2=vov_{C2} = v_o; the inputs are ViV_i and IsI_s.

Circuit

Vi --[4 ohm]--+--[2 H, iL ->]--+-------+
              |                |       |
           0.25 F            0.5 F   1 ohm   Is (up into
              |                |       |     node Vc2)
GND ----------+----------------+-------+

State variables and inputs

x1=vC1,x2=iL,x3=vC2=vo,u1=Vi,u2=Isx_1 = v_{C1},\quad x_2 = i_L,\quad x_3 = v_{C2} = v_o, \qquad u_1 = V_i,\quad u_2 = I_s

KCL at node VC1V_{C1}

Current in through 4 Ω = capacitor current + inductor current:

Vi−x14=0.25 x˙1+x2  ⇒  x˙1=−x1−4x2+Vi\frac{V_i - x_1}{4} = 0.25\,\dot{x}_1 + x_2 \;\Rightarrow\; \dot{x}_1 = -x_1 - 4x_2 + V_i

KVL around the inductor

2 x˙2=x1−x3  ⇒  x˙2=0.5x1−0.5x32\,\dot{x}_2 = x_1 - x_3 \;\Rightarrow\; \dot{x}_2 = 0.5x_1 - 0.5x_3

KCL at node VC2V_{C2}

Inductor current + source current = capacitor current + resistor current:

x2+Is=0.5 x˙3+x31  ⇒  x˙3=2x2−2x3+2Isx_2 + I_s = 0.5\,\dot{x}_3 + \frac{x_3}{1} \;\Rightarrow\; \dot{x}_3 = 2x_2 - 2x_3 + 2I_s

State equation

[v˙C1i˙Lv˙C2]=[−1−400.50−0.502−2][vC1iLvC2]+[100002][ViIs]\begin{bmatrix}\dot{v}_{C1}\\ \dot{i}_L\\ \dot{v}_{C2}\end{bmatrix} = \begin{bmatrix}-1 & -4 & 0\\ 0.5 & 0 & -0.5\\ 0 & 2 & -2\end{bmatrix} \begin{bmatrix}v_{C1}\\ i_L\\ v_{C2}\end{bmatrix} + \begin{bmatrix}1 & 0\\ 0 & 0\\ 0 & 2\end{bmatrix} \begin{bmatrix}V_i\\ I_s\end{bmatrix}

Output equation

vo=[001][vC1iLvC2]+[00][ViIs]v_o = \begin{bmatrix}0 & 0 & 1\end{bmatrix}\begin{bmatrix}v_{C1}\\ i_L\\ v_{C2}\end{bmatrix} + \begin{bmatrix}0 & 0\end{bmatrix}\begin{bmatrix}V_i\\ I_s\end{bmatrix}

Check

  • Characteristic equation: ∣sI−A∣=s3+3s2+5s+5|sI - A| = s^3 + 3s^2 + 5s + 5.
  • Vo/Vi=1s3+3s2+5s+5V_o/V_i = \dfrac{1}{s^3 + 3s^2 + 5s + 5}. At DC this is 1/51/5, which matches the resistive divider 1/(4+1)1/(4+1).
  • Vo/Is=2(s2+s+2)s3+3s2+5s+5V_o/I_s = \dfrac{2(s^2 + s + 2)}{s^3 + 3s^2 + 5s + 5}. At DC this is 0.8 Ω0.8\ \Omega, which matches 1∥4=0.8 Ω1 \parallel 4 = 0.8\ \Omega.

Answer: x˙=Ax+Bu\dot{x} = Ax + Bu, vo=Cxv_o = Cx with A=[−1−400.50−0.502−2]A = \begin{bmatrix}-1 & -4 & 0\\ 0.5 & 0 & -0.5\\ 0 & 2 & -2\end{bmatrix}, B=[100002]B = \begin{bmatrix}1 & 0\\ 0 & 0\\ 0 & 2\end{bmatrix}, C=[0    0    1]C = [0\;\;0\;\;1], D=0D = 0. (If the current source points the other way, the sign of the IsI_s column in BB changes.)

✓ Answer checked by Claude Opus 5.5 (high effort)

All Control Systems chapter 8 questions →