STATE SPACE ANALYSIS

Introduction
• State space analysis is an powerful and modern approach for the design
and analysis of control systems.
• The conventional or old methods for the design and analysis of control
systems is based on transfer function method.
• The transfer function method for design and analysis had many drawbacks
such as..
– Transfer function is defined under zero initial conditions
– Applicable to LTI systems
– SISO systems
– Does not provide the information regarding internal state of
the system
• Initial conditions can be incorporated in the system design
• State equations are highly compatible for simulation on analog or digital
computers
 Advantages of state variable analysis.
This can be applicable to
• Linear systems
• Non-linear system
• Time variant systems
• Time invariant systems
• Multiple input multiple output systems
• This gives idea about the internal state of the
system
 Concept
Input Output
SYSTEM
r(t) y(t)

The output not only depends on the input applied to the

system for t > t0, but also on the initial conditions at time t = t0
y(t) = y(t)│ t=t0 + y(t)│ t≥t0
t t
` = ‫׬‬−∞
0
y t + ‫׬‬t y(t)
0
t
= y(t0) + ‫׬‬0 y(t)
The term y(t0) is called the state of the system.
The variable that represents this state of the system is called
the state variable
 Definitions
State: The state of a dynamic system is the smallest set of
variables called state variables such that the knowledge of
these variables at time t = to (Initial condition), together
with the knowledge of input for t ≥ 𝑡0 , completely
determines the behavior of the system for any time 𝑡 ≥ 𝑡0 .
State Variables: A set of variables which describes the
system at any time instant are called state variables
State vector: A vector whose elements are the state
variables
State space: The n-dimensional space whose co-ordinate
axes consists of the x1 axis, x2 axis,.... xn axis, (where x1 , x2
,..... xn are state variables:) is called a state space.
 Illustration
Consider the circuit shown in figure
dvc(t)
i(t) = C
dt
1 t
Vc(t) = ‫׬‬−∞ i t dt
𝐶
1 t0 1 t
= ‫׬‬−∞ i t dt + ‫׬‬ i t dt
𝐶 C t0
= Vc(t0) + Vc(t)
Vc(t0) = initial voltage across capacitor
The voltage across capacitor can be taken as a
state variable
Consider the circuit shown in figure
i(t)
di(t)
VL(t) = L
dt
1 t L VL(t)
iL(t) = ‫׬‬−∞ VL(t)dt
𝐿
1 t0 1 t
= ‫׬‬−∞ VL(t) dt + ‫׬‬t
VL(t) dt
𝐿 L 0
= VL(t0) + VL(t)
VL(t0) = initial voltage across Inductor
The current through inductor can be taken as a
state variable
 State Model
consider a multi-input & multi-output system is having
m inputs 𝑢1(𝑡) , 𝑢2(𝑡) ,…….𝑢m(𝑡)
p outputs 𝑦1(𝑡) , 𝑦2(𝑡) ,…….𝑦p(𝑡)
n state variables 𝑥1(𝑡) , 𝑥2(𝑡),…….𝑥𝑛(𝑡)

The different variables may be represented by the vectors as shown below

u1(t) y1(t)
u (t) y (t)
Input vector U(t) = 2 ; Output vector Y(t) = 2
.. ..
um(t) yp(t)

x1(t)
x (t)
State variable vector X(t) = 2
..
xn(t)
 State Equations
The state variable representation can be arranged
in the form of n number of first order differential
equations as shown below:
dx1
= x1ሶ = f1(x1,x2,….xn ; u1,u2…..um)
dt
dx2
= x2ሶ = f1(x1,x2,….xn ; u1,u2…..um)
dt
……
……
dxn
= xnሶ = f1(x1,x2,….xn ; u1,u2…..um)
dt
ሶ f(X(t),U(t))
In vector notation, X(t)=
Similarly the output vector Y(t) = f(X(t),U(t))
 State Model of Linear System
The state model of a system consist of state equation and output
equation.
The state equation of a system is a function of state variables
and inputs.
For LTI systems, the first derivatives of state variables can be
expressed as a linear combination of sate variables and inputs
x1ሶ = a11x1+a12x2+----+a1nxn + b11u1+b12u2+---+b1mum
x2ሶ = a21x1+a22x2+----+a2nxn + b21u1+b22u2+---+b2mum
---
---
xnሶ = an1x1+an2x2+----+annxn + bn1u1+bn2u2+---+bnmum
where the coefficients aij and bij are constants
In the matrix form,
x1ሶ a11 a1n x1 b11 b1m u1
x2ሶ a21 ⋯ a2n x2 b21 ⋯ b2m u2
. = . . .. + . . ..
. ⋮ ⋱ ⋮ .. ⋮ ⋱ ⋮ ..
xnሶ an1 ⋯ ann xn bn1 ⋯ bnm um
Ẋ(t) = A X(t) + B U(t)……… state equation
where, A is state matrix of size (n×n)
B is the input matrix of size (n×m)
X(t) is the state vector of size (n×1)
U(t) is the input vector of size (m×1)
 Output equation
The output at any time are functions of state
variables and inputs.
output vector, Y(t) = f(x(t), U(t))
Hence the output variables can be expressed as a
linear combination of state variables and inputs.
y1= c11x1+c12x2+----+c1nxn + d11u1+d12u2+---+d1mum
y2 = c21x1+c22x2+----+c2nxn + d21u1+d22u2+---+d2mum
---
---
yp = cp1x1+cp2x2+----+cpnxn + dp1u1+dp2u2+---+dpmum
where the coefficients cij and dij are constants
In the matrix form,
𝑦1 c11 c1n x1 d11 d1m u1
𝑦2 c21 ⋯ c2n x2 d21 ⋯ d2m u2
.. = . . .. + . . ..
.. ⋮ ⋱ ⋮ .. ⋮ ⋱ ⋮ ..
𝑦𝑝 cp1 ⋯ cpn xn dp1 ⋯ dpm um

Y(t) = C X(t) + D U(t)………….. output equation

where, C is the output matrix of size (p×n)
D is the transmission matrix of size (p×m)
X(t) is the state vector of size (n×1)
Y(t) is the output vector of size (p×1)
U(t) is the input vector of size (m×1)
 State Model
Ẋ(t) = A X(t) + B U(t) state equation
Y(t) = C X(t) + D U(t) output equation
 Selection of state variables
• The state variables of a system are not unique.
• There are many choices for a given system
Guide lines:
1. For a physical systems, the number of state variables
needed to represent the system must be equal to the
number of energy storing elements present in the system
2. If a system is represented by a linear constant coefficient
differential equation, then the number of state variables
needed to represent the system must be equal to the
order of the differential equation
3. If a system is represented by a transfer function, then the
number of sate variables needed to represent the system
must be equal to the highest power of s in the
denominator of the transfer function.
State space Representation using Physical variables
• In state-space modeling of the systems, the
choice of sate variables is arbitrary.
• One of the possible choice is physical
variables.
• The state equations are obtained from the
differential equations governing the system
 State Space Model
Consider the following series of the RLC circuit.
It is having an input voltage vi(t) and the current
flowing through the circuit is i(t).
• There are two storage elements (inductor and
capacitor) in this circuit. So, the number of the
state variables is equal to two.
• These state variables are the current flowing
through the inductor, i(t) and the voltage
across capacitor, vc(t).
• From the circuit, the output voltage, v0(t) is
equal to the voltage across capacitor, vc(t).
Y(t) = v0(t) = vc(t)
Apply KVL around the loop,
di(t)
Vi(t) = R i(t) + L + vc(t)
dt
di(t) R 1 1
𝑖 𝑡ሶ = = - i(t) - vc(t) + Vi(t)
dt L L L
The voltage across the capacitor is
1
vc(t) = ‫ ׬‬i t dt
c
Differentiate the equation with respect to time,
dv (t) i(t)
𝑣 𝑡ሶ = c =
𝑐 dt c

di(t)
i(t) dt
State vector, X = ; Differential state vector, Xሶ =
vc(t) dvc(t)
dt
Arrange the differential equations and output equation into standard
form of state space model as,
di(t) −R −1
1
dt L L i(t)
Xሶ = = 1 + L vi(t)
dvc(t)
0 vc(t) 0
dt C

i(t)
Y= 0 1
vc(t)

Ẋ(t) = A X(t) + B U(t)

Y(t) = C X(t) + D U(t)
−R −1
1
L L
Here A = 1 ; B= L ; C= 0 1 ; D = [0]
0 0
C
 Problem
Represent the electrical circuit shown by a state
model
 Solution
Since there are three energy storing elements, choose
three state variables to represent the systems
The current through the inductors i1,i2 and voltage
across the capacitor vc are taken as state variables
Let the three sate variables be x1, x2 and x3 be related
to physical quantities as shown
Let, i1 = x1,
i2 = x2,
vc = x3
Applying KVL to loop 1,
di1(t)
L1 + R1i1(t) + u(t) – vc(t) = 0
dt
dx1(t)
⟹ L1 + R1x1(t) + u(t) – x3(t) = 0
dt
⟹ L1x1ሶ (t) + R1x1(t) + u(t) – x3(t) = 0
R1 1 1
⟹ x1ሶ (t) = - x1(t) + x3(t) - u(t)-------(1)
L1 L1 L1
Applying KVL to loop 2,
di2(t)
L2 + R2i2(t) – vc(t) = 0
dt
dx2(t)
⟹ L2 + R2x2(t) – x3(t) = 0
dt
⟹ L2 x2ሶ (t) + R2x2(t) – x3(t) = 0
R2 1
⟹ x2ሶ (t) = - x2(t) + x3(t) ------------------(2)
L2 L2
Applying KCL at node vc(t),
dvc(t)
i1(t) + i2(t) + C =0
dt
dx3(t)
⟹ x1(t) + x2(t) + C =0
dt
1 1
⟹ x3ሶ (t) = - x1(t) - x2(t) --------------(3)
C C
Putting 1, 2and 3 in matrix form,
R1 1
− 0 1
x1ሶ (t) L 1 L 1 x 1(t) −
R2 1 x (t) + L1
x2ሶ (t) = 0 − 2 0 u(t)
x3ሶ (t) L2 L2
1 1 x3(t) 0
− − 0
C C

This is State Equation

y(t) = R2i2(t) = R2x2(t)
x1(t)
y(t) = 0 R2 0 x2(t)
x3(t)
This is output equation
 Problem
Obtain the state model for a system represented
by an electrical system as shown in figure
 Solution
Since there are two energy storage elements
present in the system, assume two state
variables to describe the system behavior.
Let the two state variables be x1 and x2 be
related to physical quantities as shown

Let v1(t) = x1(t)

v2(t) = x2(t)
Applying KCL at node v1(t),
v1 t −u(t) dv1(t) v1 t −v2(t)
+C + =0
R dt R
x1 t −u(t) dx1(t) x1 t −x2(t)
⟹ +C + =0
R dt R
2x1 t u(t) x2(t)
⟹ - + C x1ሶ (t) - =0
R R R
2x1 t x2(t) u(t)
⟹ x1ሶ (t) = - + + ------------(1)
RC RC RC
Applying KCL at node v2(t),
dv2(t) v2 t −v1(t)
C + =0
dt R
dx2(t) x2 t −x1(t)
⟹ C + =0
dt R
x1(t) x2 t
⟹ C x2ሶ (t) - + =0
R R
x1 t x2(t)
⟹ xሶ 2(t)= - ----------(2)
RC R𝐶
putting 1 and 2 in matrix form,
−2 1
1
x1ሶ (t) RC RC x1(t)
= 1 −1 + RC u(t)
x2ሶ (t) x2(t) 0
RC RC
This is the state equation
The output of the circuit is given by
y(t) = v2(t)
= x2(t)
x1(t)
= 0 1
x2(t)
This is the output equation
 Problem
Represent the electrical network by a state
equation
State representation using Phase variables
• The phase variables are defined as those particular
state variables which are obtained from one of the
system variables and its derivatives.
• Usually the variables used is the system output and
the remaining state variables are then derivatives of
the output.
• The state model using phase variables can be easily
determined if the system model is already known in
the differential equation or transfer function form.
Consider the following nth order linear differential equation relating the
output y(t) to the input u(t) of a system.
dn y d𝑛−1𝑦 d𝑛−2𝑦 dy
+ 𝑎 1 dt𝑛−1 + 𝑎2 dt𝑛−2 +-----+ 𝑎𝑛−1 dt + 𝑎𝑛 𝑦 = u
dtn
Let us define the state variables as
𝑥1 = 𝑦
dy 𝑑𝑥1
𝑥2 = =
dt 𝑑𝑡

d2y dẏ 𝑑𝑥2
𝑥3 = = =
dt2 dt 𝑑𝑡
⋮⋮⋮
d𝑛−1𝑦 𝑑𝑥𝑛−1
𝑥𝑛 = 𝑛−1 =
dt 𝑑𝑡
From the above equations we can write
x1ሶ = 𝑥2
xሶ 2 = 𝑥3
⋮⋮
xሶ 𝑛−1 = 𝑥𝑛
xሶ 𝑛 + 𝑎1𝑥𝑛 + …..+ 𝑎 𝑛−1𝑥2 +𝑎𝑛𝑥1 = 𝑢

xሶ 𝑛 = − 𝑎𝑛𝑥1 − 𝑎 𝑛−1𝑥2……..−𝑎1𝑥𝑛 + 𝑢
writing the above state equation in vector matrix form,
ሶ = AX(t) + Bu(t)
X(t)

x1ሶ 0 0 x1 0
x2ሶ ⋯ x2 0
= 0 0 + [u]
⋮ ⋱ ⋮ .
.
.
x𝑛ሶ −𝑎𝑛 ⋯ −𝑎1 x𝑛 1

Output equation can be written as

x1
x2
Y(t) = C X(t)= 1 0 0 . .
.
x𝑛
 Problem
Construct a state model for a system characterized
by the differential equation,
d3y d2y dy
+6 + 11 + 6y + u = 0
dt3 dt2 dt
Also give the block diagram representation of the
state model
Solution: Let us choose y and their derivatives as
state variables. The system is governed by third
order differential equation, so the number of state
variables required are three.
Let the state variables x1, x2 and x3 are related to
phase variables as follows.
x1 = y
dy
x2 = = x1ሶ
dt
d2y dx2
x3 = 2 = = x2ሶ
dt dt
dy d2y d3y
Put y = x1, = x2, 2 = x3 and = xሶ 3 in the
dt dt dt3
given equation
∴ xሶ 3+ 6x3 + 11x2 + 6x1 + u = 0
⟹ xሶ 3 = - 6x1- 11x2 - 6x3 - u
The state equations are
xሶ 1 = x2
xሶ 2 = x3
xሶ 3 = - 6x1- 11x2 - 6x3 - u
Arranging the state equations in the matrix form,
x1ሶ 0 1 0 x1 0
x2ሶ = 0 0 1 x2 + 0 [u]
x3ሶ −6 −11 −6 x3 −1
Here y = output
But y = x1
x1
∴ The output equation is , y = [1 0 0] x2
x3
The block diagram for the state model is
 Problem
Represent the differential equation given below
in a state model
d3y d2y dy
+ +6 +7y = 2u(t)
dt3 dt2 dt
Solution:
Since, the given equation is a third-order
differential equation, choose three state
variables to represent the system
 Let y(t) = x1(t)
y(t)
ሶ = x2(t)
y(t)
ሷ = x3(t)
where x1(t), x2(t), x3(t) are the state variables of
the system.
y(t) = x1(t)
y(t)ሶ = x2(t) = x1ሶ (t) -----------(1)
y(t) ሷ = x3(t) = x2ሶ (t) -----------(2)
From the given diff. equation,
d3y d2y dy
+ +6 +7y = 2u(t)
dt3 dt2 dt
𝑦ഺ (t) + 𝑦(t)
ሷ + 6𝑦(t) ሶ +7y = 2u(t)
x3ሶ (t) + x3(t) + 6 x2(t) + 7 x1(t) = 2u(t)
x3ሶ (t) = -7 x1(t) - 6 x2(t) - x3(t) + 2u(t) ------(3)
x2ሶ (t) = x3(t) ------(2)
x1ሶ (t) = x2(t) ------(1)
Putting the above equations in matrix form,
x1ሶ 0 1 0 x1 0
x2ሶ = 0 0 1 x2 + 0 u(t)
x3ሶ −7 −6 −1 x3 2
This is the state equation
The output expression is y(t) = x1(t)
x1
= 1 0 0 x2
x3
This is the output equation
 Problem
Obtain the state model for the signal flow graph
given below:
 Problem
Obtain the state model of the system whose transfer
Y(s) 24
function is given by = 3 2
U(s) 𝑠 +9𝑠 +26𝑠+24
Solution:
Y(s) 24
=
U(s) 𝑠3+9𝑠2+26𝑠+24
Cross-multiplying yields
[s3 + 9s2 + 26s + 24] Y(s) = 24 U(s)
s3 Y(s) + 9s2 Y(s) + 26sY(s) + 24 Y(s) = 24U(s)
Taking inverse Laplace transforms,
d3y(t) d2y(t) dy(t)
+9 + 26 + 24 y(t) = 24u(t)
dt3 dt2 dt
ഺ
y(t) + 9y(t)
ሷ + 26y(t)
ሶ + 24 y(t) = 24u(t)
Choosing the state variables as successive derivatives
x1(t) = y(t)
x2(t) = y(t)
ሶ
x3(t) = y(t)
ሷ

x1(t) = y(t)
x1ሶ (t) = x2(t) = y(t)
ሶ -----(1)
x2ሶ (t) = x3(t) = 𝑦(t)
ሷ -----(2)
x3ሶ (t) = 𝑦ഺ(t)
ഺ
y(t) + 9y(t)ሷ + 26y(t) ሶ + 24 y(t) = 24u(t)
x3ሶ (t) + 9x3(t) + 26x2(t) + 24x1(t) = 24u(t)
x3ሶ (t) = - 24x1(t) - 26x2(t) - 9x3(t) + 24u(t)------(3)
Putting equations 1, 2 and 3 in matrix form,
x1ሶ 0 1 0 x1 0
x2ሶ = 0 0 1 x2 + 0 [u]
x3ሶ −24 −26 −9 x3 24
The output expression is y(t) = x1(t)
x1
= 1 0 0 x2
x3
 Matlab
Obtain the state model of the system whose
Y(s) 24
transfer function is given by =
U(s) 𝑠3+9𝑠2+26𝑠+24
using Matlab
 Problem
Obtain the state model of the system whose transfer function is
Y(s) 1
given by =
U(s) s2+s+1
Y(s) 1
Solution: =
U(s) s2+s+1

⟹ (s2 + s + 1)Y(s) = U(s)

s2 Y(s) + s Y(s) + Y(s) = U(s)
Taking inverse Laplace transform on both sides,
d2y(t) dy(t)
+ + y(t) = u(t)
dt2 dt

Let y(t) = x1
𝑑𝑦(𝑡)
= x2 = x1ሶ and u(t) = u
𝑑𝑡
Then the state equation is, x2ሶ = -x1-x2+u
The output equation is, y(t) = y = x1
The state space model is
x1ሶ 0 1 x1 0
xሶ = == + [u]
x2ሶ −1 − 1 x2 1

x1
Y= 1 0
x2
 Problem
Obtain the state model of the system whose
s2+7s+2
transfer function is given by 3 2
s +9s +26s+24
Solution:
Y(s) s2+7s+2
=
U(s) s3+9s2+26s+24
Y(s) Y(s) C(s)
= x
U(s) C(s) U(s)
Y(s)
= s2 + 7s + 2 --------(1)
C(s)
C(s) 1
= -------(2)
U(s) s3+9s2+26s+24
 C(s) 1
Consider equation (2), =
U(s) s3+9s2+26s+24

Cross-multiplying on both sides,

[s3 + 9s2 + 26s + 24] C(s) = U(s)
s3 C(s) + 9s2 C(s) + 26s C(s) + 24 C(s) = U(s)
Taking inverse Laplace transform,
d3c(t) d2c(t) dc(t)
+ 9 + 26 + 24 c(t) = u(t)
dt3 dt2 dt
ഺc (t) + 9 c(t)
ሷ + 26 c(t)
ሶ +24 c(t) = u(t)
x1(t) = c(t)
x1ሶ (t) = x2(t) = c(t)
ሶ -----(3)
x2ሶ (t) = x3(t) = c(t)
ሷ -----(4)
x3ሶ (t) = ഺc(t)
ഺc (t) + 9c(t)
ሷ + 26c(t)
ሶ + 24c(t) = u(t)
x3ሶ (t) + 9 x3(t) + 26x2(t) + 24x1(t) = u(t)
x3ሶ (t) = - 24x1(t) - 26x2(t) - 9 x3(t) + u(t) -----(5)
Putting equations 3, 4 and 5 in matrix form,
x1ሶ 0 1 0 x1 0
x2ሶ = 0 0 1 x2 + 0 [u]
x3ሶ −24 −26 −9 x3 1
Consider equation (1),
Y(s)
= s2 + 7s + 2
C(s)

Y(s) = [s2 + 7s + 2]C(s)

Y(s) = s2C(s) + 7sC(s) + 2C(s)
Taking inverse Laplace transform,
y(t) = c(t)
ሷ + 7c(t)
ሶ +2c(t)
y(t) = 2x1(t) + 7x2(t) + x3(t)
x1
y(t) = 2 7 1 x2
x3
 Matlab
Obtain the state model of the system whose
s2+7s+2
transfer function is given by 3 2 using
s +9s +26s+24
Matlab
 Problem
A feedback system has a closed-loop transfer
Y(s) 2(s+5)
function =
U(s) (s+2)(s+3)(s+4)
Solution:
Y(s) 2(s+5)
=
U(s) (s+2)(s+3)(s+4)
By partial fraction expansion,
Y(s) 2(s+5) A B C
= = + +
U(s) (s+2)(s+3)(s+4) (s+2) (s+3) (s+4)
Solving for A, B and C
A = 3; B = - 4; C=1
Y(s) 2(s+5) 3 4 1
= = - + -----(1)
U(s) (s+2)(s+3)(s+4) (s+2) (s+3) (s+4)
3 4 1
= - +
s(1+2/s) s(1+3/s) s(1+4/s)
1 1 1
s s s
= 1 x3 - 1 x4+ 1
(1+ ∗2) (1+ ∗3) (1+ ∗4)
s s s
1 1 1
s s s
∴Y s =[ 1 x3 - 1 x4+ 1 ]u(s)
(1+ ∗2) (1+ ∗3) (1+ ∗4)
s s s
1 1 1
=[ s
1
(1+s ∗2)
x 3]U(s) - [ s
1
(1+s ∗3)
x 4] U(s) + [ s
1
(1+s ∗4)
]u(s)
The equation can be represented by the block
diagram as shown
Assign state variables at the output of the
integrators as shown.
At the input of the integrators , first derivative of
the state variables are present.
The state equations are formed by adding all the
incoming signals to the integrator and equating
to the corresponding first derivative of state
variables.
The state equations are
xሶ 1 = -2x1 + u
xሶ 2 = -3x2 + u
xሶ 3 = -4x3 + u
The output equation is
y = 3x1- 4x2 + x3
The State model is given by,
xሶ 1 −2 0 0 x1 1
xሶ 2 = 0 −3 0 x2 + 1 u
xሶ 3 0 0 −4 x3 1
x1
Y = 3 − 4 1 x2
x3
 Matlab
Find the state model for the transfer function
2(s+5)
using Matlab
(s+2)(s+3)(s+4)
 Solution of State Equation
S-Domain:
The State equation of nth order system is given by,
x(t)
ሶ = A X(t) + BU(t); X(0) = X0 = initial condition vector
Taking Laplace transforms on both sides,
S X(s)-X(0) = A X(s) + B U(s)
X(s)[sI-A] = X(0) + B U(s) where I = unit matrix
X(s) = [sI-A]-1 X(0) + [SI-A]-1 B U(s) --------------(1)
Taking inverse Laplace transforms on both sides,
X(t) = L-1 [sI-A]-1 X(0) + L-1 [SI-A]-1 B U(s)
where L-1 [sI-A]-1 = ∅(t) = state transition matrix
[sI-A]-1 = ∅(s) = Resolvent matrix
The solution of state equation is,
X(t) = ∅(t) X(0) + L-1 [∅(s). B U(s)]
The output equation is,
y(t) = C x(t) + D u(t)
Taking Laplace transforms on both sides,
Y(s) = C X(s) + D U(s)
From equation (1), X(s) = [sI-A]-1 X(0) + [SI-A]-1 B U(s)
= C{[sI-A]-1 X(0) + [SI-A]-1 B U(s)} + D U(s)
= C[sI-A]-1 X(0) + C [SI-A]-1 B U(s) + D U(s)
For zero initial conditions, X(0) = 0
∴ Y(s) = C [SI-A]-1 B U(s) +D U(s)
= {C [SI-A]-1 B +D} U(s)
Y(s)
The transfer function = = C [SI-A]-1 B + D
U(s)
 Problem
A state variable description of a system is given
by the matrix equation,
ሶX = −1 0 X + 1 u
1 −2 0
Y = [ 1 1] X
Find (i) The Transfer function
(ii) The State transition matrix
(iii) State diagram
 Solution
The state model is given by
Xሶ = A X + B U
Y=CX+DU
From the given problem,
−1 0 1
A= B= C = [ 1 1]
1 −2 0
Y(s)
(i) The transfer function = = C [SI-A]-1 B + D
U(s)
Here D =0
Y(s)
∴ = C [SI-A]-1 B
U(s)
 1 0 −1 0
[sI-A] = s -
0 1 1 −2
s 0 −1 0
= -
0 s 1 −2
s+1 0
=
−1 s + 2
-1 Adj A s+2 0 1
[SI-A] = =
Det A 1 s+1 (s+1)(s+2)
Y(s)
= C [SI-A]-1 B
U(s)
s+2 0 1 1
= [1 1] { }
1 s+1 (s+1)(s+2) 0
1 s+2 0 1
={ } [1 1]
(s+1)(s+2) 1 s+1 0
 Y(s) 1 s+2
= { } [1 1]
U(s) (s+1)(s+2) 1
1
= { } [s+3]
(s+1)(s+2)
s+3
=
(s+1)(s+2)

(ii) State transition matrix = ∅(t) = L-1 [sI-A]-1

-1 s+2 0 1
[SI-A] =
1 s + 1 (s+1)(s+2)
s+2 0
(s+1)(s+2) (s+1)(s+2)
= 1 s+1
(s+1)(s+2) (s+1)(s+2)
 1
0
(s+1)
[SI-A]-1 = 1 1
(s+1)(s+2) (s+2)

∅(t) = L-1 [sI-A]-1

1
0
-1 (s+1)
= L 1 1
(s+1)(s+2) (s+2)
1
0
(s+1)
= L-1 1 1 1
−
(s+1) (s+2) (s+2)

= e−t 0
e−t − e−2t e−2t
(iii) State Diagram:
The state equation is,
xሶ 1 = - x1 + u(t)
xሶ 2 = x1 – 2x2
The output equation is y = x1 + x2
 Matlab
Find the transfer function using Matlab
 Problem
The state equation of a LTI system is given as
0 5 1
xሶ = X+ u and y = [1 1] X
−1 −2 1
Determine (i) State transition matrix
(ii) The transfer function
(iii) State diagram
 Solution
From the given system,
0 5 1
A= B= C= [1 1]
−1 −2 1
(i) The State transition matrix ,
∅(t) = L-1 [sI-A]-1
1 0 0 5 s −5
[sI – A] = s - =
0 1 −1 −2 1 s+2
s+2 5
s+2 5 1 s s+2 +5 s s+2 +5
[sI-A]-1 = = −1 s
−1 s s s+2 +5
s s+2 +5 s s+2 +5
 s+1+1 5
s+1 2+22 s+1 2+22
= −1 s
s+1 2+22 s+1 2+22

∅(t) = L−1 [sI−A]−1

s+1+1 5
s+1 2+22 s+1 2+22
= L-1 −1 s+1−1
s+1 2+22 s+1 2+22

1 5
e−t cos 2t + e−t sin 2t e−t sin 2t
2 2
= 1 1
− e−t sin 2t e−t cos 2t − e−t sin 2t
2 2
(ii) The transfer function
Y(s)
= C [SI-A]-1 B
U(s)
s+2 5 1 1
= [1 1]
−1 s s s+2 +5 1
1 S+7
= [1 1]
s s+2 +5 S−1
2S+6
=
s s+2 +5
 Matlab
Find the transfer function using Matlab
 Problem
Find the transfer function of a state model of a
system given by ,
0 1 0 0 0 u
xሶ = 0 0 1 X + 1 0 u1
2
−1 −2 −3 0 1

y1 1 0 0
and y = X
2 0 0 1
 Solution
From the given system,
0 1 0 0 0
1 0 0
A= 0 0 1 ; B= 1 0 ;C =
0 0 1
−1 −2 −3 0 1
1 0 0 0 1 0
[sI-A] = s 0 1 0 - 0 0 1
0 0 1 −1 −2 −3
s −1 0
= 0 s −1
1 2 s+3
 (s + 2)(s + 1) s+3 1
1
∴ [sI-A]-1 = −1 s(s + 3) s
s3+3s2+2s+1
−s −(2s + 1) s2
The transfer function of the system is given by
Y(s)
= C [SI-A]-1 B
U(s)

(s + 2)(s + 1) s+3 1 0 0
1 0 0 −1 s(s + 3) s
1
= 1 0
0 0 1 s3+3s2+2s+1
−s −(2s + 1) s2 0 1
s+3 1
1 0 0 s(s + 3) s 1
=
0 0 1 s3+3s2+2s+1
− 2s + 1 s2
s+3 1 1
=
− 2s + 1 s2 s3+3s2+2s+1
 Matlab
Find the transfer function using Matlab
 Problem
A linear time-invariant system is characterized
xሶ 1 1 0 x1
by state equation = . Compute
xሶ 2 1 1 x2
the solution of the state equation , assuming the
1
initial vector X0 =
0
 Solution
1 0
From the given system, A =
1 1
The solution of state equation is,
X(t) = L-1 [sI-A]-1 X(0) + L-1 [SI-A]-1 B U(s)
Here U= 0
∴ X(t) = L-1 [sI-A]-1 X(0)
1 0 1 0 s−1 0
[sI – A] = s - =
0 1 1 1 −1 s−1

s−1 0 1
[sI-A]-1 = 2
1 s−1 s−1
 s−1 0 1
[sI-A]-1 = 2
1 s−1 s−1
1
0
s−1
= 1 1
s−1 2 s−1

X(t) = L-1 [sI-A]-1 X(0)

1
0 1
s−1
= L-1 1 1 0
s−1 2 s−1
t 1 t
e
= t 0 e
= t
te et 0 te
Solution of state equation (Time Domain)
x(t)
ሶ = A x(t) + B u(t) X(0) = x0
x(t)
ሶ - A x(t) = B u(t)
pre-multiplying both sides by e−At
e−At [x(t)
ሶ - A x(t)] = e−At B u(t) ----------------(1)
𝑑
Consider, {e−At x(t)} = e−At x(t)
ሶ - A e−At x(t)]
𝑑𝑡
= e−At [x(t)
ሶ - A x(t)]
𝑑
∴ equation (1) = {e−At x(t)} = e−At B u(t)
𝑑𝑡
Integrating with respect to t
t d t −Aτ
‫׬‬0 dt {e−At
x(t)} dt = ‫׬‬0 [e B u(τ)]dτ
−At 𝑡 −Aτ
e x(t) – x(0) = ‫׬‬0 [e B u(τ)]d𝜏
Pre-multiplying both sides by eAt ,
At −At At 𝑡 −Aτ
e [e x(t) – x(0)] = e [‫׬‬0 e B u(τ)d𝜏]
At t −Aτ
x(t) = e [x(0) + ‫׬‬0 e B u(τ)dτ]
At t A(t−τ )
= e x(0) + ‫׬‬0 [e B u(τ)]dτ
t
x(t) = ∅ t x(0) + ‫׬‬0 ∅(t−τ) B u(τ) dτ
if the initial time is t = t0, the solution of state equation
becomes,
t
x(t) = ∅ t − t0 x(t0) + ‫׬‬0 ∅(t−τ) B u(τ)dτ
 Properties of state transition matrix
∅(t) = eAt = L-1 [sI-A]-1

1. ∅(0) = I
2. ∅-1(t) = ∅(-t)
3. ∅(t2-t1) ∅(t1-t0) = ∅(t2-t0) for any t2, t1, t0
4. [∅(t)]k = ∅(kt)
5. ∅(t1+t2) = ∅(t1) ∅(t2) = ∅(t2) ∅(t1)
 Problem
Compute the State transition matrix by infinite
0 1
series method A =
−1 −2
Solution: For the given system matrix A, the state
transition matrix is,
2 3
At At
∅ (t) = eAt = I + At + + +------
2! 3!
0 1
A=
−1 −2
2 0 1 0 1 −1 −2
A = A. A = . =
−1 −2 −1 −2 2 3
3 2 −1 −2 0 1 2 3
A = A .A = . =
2 3 −1 −2 −3 −4
 2 3
At At
∅(t) = I + At + + +---
2! 3!
1 0 0 1 −1 −2 t2 2 3 t3
= + t+ + +- -
0 1 −1 −2 2 3 2! −3 −4 3!

t2 t3 t3
1 − + +⋯ t − t2 + + ⋯
2 3 2
= t 3
3t2
−t + t2 − + ⋯ 1 − 2t + + ⋯
2 2

−t −t
= e + te te−t
−te−t e−t − te−t
 Problem
Find the state transition matrix by infinite series
1 1
method for the system matrix A =
0 1
Solution: For the given system matrix A, the
state transition matrix is,
2 3
At At
∅ (t) = eAt = I + At + + +------
2! 3!
1 1
A=
0 1
1 1 1 1 1 2
A2 = A.A = . =
0 1 0 1 0 1
 1 2 1 1 1 3
A3 = A2. A= . =
0 1 0 1 0 1
2 3
At At
∅(t) = I + At + + + ---
2! 3!
1 0 1 1 1 2 t2 1 3 t3
= + t+ + + ---
0 1 0 1 0 1 2! 0 1 3!

t2 t3 t3
1+t+ + + ⋯ t + t2 + +⋯
2! 3! 2
= t2 t3
0 1+t+ + +. .
2! 3!
t
= e tet
0 et
 References
1. Control Engineering by Nagrath & Gopal, New
Age International Publishers
2. Engineering control systems - Norman S. Nise,
John WILEY & sons , fifth Edition
3. Modern control Engineering-Ogata, Prentice
Hall
4. Automatic Control Systems- B.C Kuo, John
Wiley and Sons