Menu
  • branch
    • thinklab
    • insights
      • automating a lunar process plant
      • the future of process control
  • loops
    • Small Modular Reactor
    • Reduction Cell
    • Crude Oil Distillation Column
    • Adaptive Control of Ball Mill Grinding Circuits
    • Autoclave Control – a regulatory approach
    • Sliding Mode Control – a SAG Mill application
  • library
    • MODE
  • branch
    • thinklab
    • insights
      • automating a lunar process plant
      • the future of process control
  • loops
    • Small Modular Reactor
    • Reduction Cell
    • Crude Oil Distillation Column
    • Adaptive Control of Ball Mill Grinding Circuits
    • Autoclave Control – a regulatory approach
    • Sliding Mode Control – a SAG Mill application
  • library
    • MODE
home/Knowledge Base/State Space/SS.04 / Transformation from Classical Form to State-Space Representation +

SS.04 / Transformation from Classical Form to State-Space Representation +

71 views 0

The block diagram provides a convenient method for deriving a set of state equations for a system that is specified in terms of a single input/output differential equation. A set of n state variables can be identified as the outputs of integrators in the diagram, and state equations can be written from the conditions at the inputs to the integrator blocks (the derivative of the state variables). There are many methods for doing this; we present here one convenient state equation formulation that is widely used in control system theory.

Let the differential equation representing the system be of order n, and without loss of generality assume that the order of the polynomial operators on both sides is the same:

(1)   \begin{equation*} (a_{n}s^{n} + a_{n-1}s^{n-1} + \dots + a_{0}) Y(s) = (b_{n}s^{n} + b_{n-1}s^{n-1} + \dots + b_{0}) U(s). \end{equation*}

We may multiply both sides of the equation by s^{−n} to ensure that all differential operators have been eliminated:

(2)   \begin{equation*} (a_{n} + a_{n-1}s^{-1} + \dots + a_{1}s^{-(n-1)} + a_{0}s^{-n}) Y(s) = (b_{n} + b_{n-1}s^{-1} + \dots + b_{1}s^{-(n-1)} + b_{0}s^{-n}) U(s). \end{equation*}

from which the output may be specified in terms of a transfer function. If we define a dummy variable Z(s), and split Eq. (2) into two parts

(3)   \begin{equation*} Z(s) = \frac{1}{a_{n} + a_{n-1}s^{-1} + \dots + a_{1}s^{-(n-1)} + a_{0}s^{-n}} U(s) \end{equation*}

(4)   \begin{equation*} Y(s) = (b_{n} + b_{n-1}s^{-1} + \dots + b_{1}s^{-(n-1)} + b_{0}s^{-n}) Z(s), \end{equation*}

Eq. (3) may be solved for U(s),

(5)   \begin{equation*} U(s) = (a_{n} + a_{n-1}s^{-1} + \dots + a_{1}s^{-(n-1)} + a_{0}s^{-n}) X(s) \end{equation*}

and rearranged to generate a feedback structure that can be used as the basis for a block diagram:

(6)   \begin{equation*} Z(s) = \frac{1}{a_{n}} U(s) - (\frac{a_{n-1}}{a_{n}}\frac{1}{s} + \dots + \frac{a_{1}}{a_{n}}\frac{1}{s^{-n}}) Z(s). \end{equation*}

The dummy variable Z(s) is specified in terms of the system input u(t) and a weighted sum of successive integrations of itself. Figure 1.1 shows the overall structure of this direct-form block diagram. A string of n cascaded integrator (1/s) blocks, with Z(s) defined at the input to the first block, is used to generate the feedback terms, Z(s)/s^{−i}, i = 1, \dots n, in Eq. (6).

Eq. (4) serves to combine the outputs from the integrators into the output y(t).

Figure 1.1 [Block diagram of a system represented by a classical differential equation].

A set of state equations may be found from the block diagram by assigning the state variables x_{i}(t) to the outputs of the n integrators. Because of the direct cascade connection of the integrators, the state equations take a very simple form. By inspection:

(7)   \begin{equation*} \begin{split} & \dot{x}_{1} = x_{2} \\ & \dot{x}_{2} = x_{3} \\ & \vdots  \\ & \dot{x}_{n-1} = x_{n} \\ & \dot{x}_{n} = -\frac{a_{0}}{a_{n}}x_{1} - \frac{a_{1}}{a_{n}}x_{2} - \dots - \frac{a_{n-1}}{a_{n}} x_{n} + \frac{1}{a_{n}}u(t). \end{split} \end{equation*}

In matrix form these equations are

(8)   \begin{equation*} \begin{bmatrix} \dot{x}_{1} \\ \dot{x}_{2} \\ \vdots \\ \dot{x}_{n-2} \\ \dot{x}_{n-1} \\ \dot{x}_{n} \end{bmatrix} = \begin{bmatrix} 0 & 1 & \dots & 0 & 0 \\ 0 & 0 & \dots & 0 & 0 \\ \vdots & \vdots & \ddots & \vdots & \vdots \\ 0 & 0 & \dots & 1 & 0 \\ 0 & 0 & \dots & 0 & 1 \\ -a_{0}/a_{n} & -a_{1}/a_{n} & \dots & -a_{n-2}/a_{n} & -a_{n-1}/a_{n} \\\end{bmatrix} \begin{bmatrix} x_{1} \\ x_{2} \\ \vdots \\ x_{n-2} \\ x_{n-1} \\ x_{n}\end{bmatrix} + \begin{bmatrix} 0 \\ 0 \\ \vdots \\ 0 \\ 0 \\ 1/a_{n}\end{bmatrix} u(t). \end{equation*}

The \mathrm{A} matrix has a very distinctive form. Each row, except the bottom one, is filled of zeroes except for a one in the position just above the leading diagonal. Eq. (8) is a common form of the state equations, used in control system theory and known as the phase variable or companion form. This form leads to a set of state variables which may not correspond to any physical variables within the system.

The corresponding output relationship is specified by Eq. (4) by noting that X_{i}(s) = Z(s)/s^{(n+1−i)}.

(9)   \begin{equation*} y(t) = b_{0}x_{1} + b_{1}x_{2} + b_{2}x_{3} + \dots + b_{n-1}x_{n} + b_{n}z_{t}. \end{equation*}

But z(t) = dx_{n}/dt, which is found from the nth state equation in Eq. (7). When substituted into Eq. (9) the output equation is:

(10)   \begin{equation*} Y(s) = \begin{bmatrix} (b_{0} - \frac{b_{n}a_{0}}{a_{n}}) & (b_{1} - \frac{b_{n}a_{1}}{a_{n}}) & \dots & (b_{n-1} - \frac{b_{n}a_{n-1}}{a_{n}}) \end{bmatrix} \begin{bmatrix} x_{1} \\ x_{2} \\ \vdots \\ x_{n} \end{bmatrix} + \frac{b_{n}}{a_{n}}u(t). \end{equation*}

\rule{20cm}{0.02cm}

Example

Draw a direct form realization of a block diagram, and write the state equations in phase variable form, for a system with the differential equation.

(11)   \begin{equation*} \frac{d^{3}y}{dt^{3}} + 7 \frac{d^{2}y}{dt^{2}} + 19 \frac{dy}{dt} + 13y = 13 \frac{du}{dt} + 26u \end{equation*}

Solution: The system order is 3, and using the structure shown in Fig. 1.1 the block diagram is as shown in Fig. 1.2.

The state and output equations are found directly from Eqs. (8) and (10):

(12)   \begin{equation*} \begin{bmatrix} \dot{x}_{1} \\ \dot{x}_{2} \\ \dot{x}_{3} \end{bmatrix} = \begin{bmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -13 & -19 & -7 \end{bmatrix} \begin{bmatrix} x_{1} \\ x_{2} \\ x_{3} \end{bmatrix} + \begin{bmatrix} 0 \\ 0 \\1 \end{bmatrix} u(t), \end{equation*}

(13)   \begin{equation*} y(t) = \begin{bmatrix} 26 & 13 & 0\end{bmatrix} \begin{bmatrix} x_{1} \\ x_{2} \\ x_{3} \end{bmatrix} + \begin{bmatrix} 0 \end{bmatrix} u(t). \end{equation*}

Figure 1.2 [Block diagram of the transfer operator of a third-order system found by a direct realization].

\rule{20cm}{0.02cm}

State Space
  • SS.04 / Transformation from Classical Form to State-Space Representation +
  • SS.03 / The Matrix Transfer Function +
  • SS.05 / Transformation From State-Space Equations to Classical Form +
View All 5  
KB Categories
  • PID
  • MPC
  • State Space
  • Optimisation
  • Feedback Control
  • Feedforward Control
  • Statistics
  • Control Theory
  • Genetic Algorithms
  • Control Performance
  • Models
State Space
  • SS.04 / Transformation from Classical Form to State-Space Representation +
  • SS.03 / The Matrix Transfer Function +
  • SS.05 / Transformation From State-Space Equations to Classical Form +
  • SS.01 / Linear Systems Described by State Equations +
  • SS.02 / State-Space Representation of Linear Time-Invariant Systems +