02 State Space Realizations
02 State Space Realizations
Introduction
This lecture will begin with a look at the concept of observability and controllability. A system is observable if the states can
be observed at the output, likewise a system is controllable if the states can be controlled from the input. The ideas of
observability and controllability are crucial to most state space control design methods.
In the previous topic, it was shown how the TF can be determined from the state-space matrices 𝐴̅, 𝐵̅, 𝐶̅ and 𝐷 ̅ . Realization is
the reverse process, i.e. determining 𝐴̅, 𝐵̅, 𝐶̅ and 𝐷
̅ from the TF. While the former process is a one-to-one mapping, realization
is a one-to-many mapping. This topic will also look at a number of different approaches to realization.
u 𝑥ሶ 1 𝑥1 +
y
+ ∫
- +
1
I
𝑥ሶ 2 𝑥2
∫
-
2
II
𝑥ሶ 3 𝑥3
+
∫
-
3
III
𝑥ሶ 4 𝑥4
∫
-
4
IV
Figure 2.1
(As an exercise confirm that this diagram does represent the above state space matrices).
States x1 and x2 are said to be observable in that they contribute to the output y. Likewise, states x1 and x3 are said to be
controllable in that they are influenced by the input u. Any state space system can be partitioned into four subspaces.
1. observable and controllable, e.g. space I in earlier diagram
2. observable and uncontrollable, e.g. II
3. unobservable and controllable, e.g. III
9
4. unobservable and uncontrollable, e.g. IV
Any system which has an uncontrollable subspace (II and/or IV) is said to be uncontrollable. Likewise, any system with an
unobservable subspace (III and/or IV) is said to be unobservable.
𝑠+1 0 0 0
̅ ̅
[𝑠𝐼 − 𝐴] = [ 0 𝑠 + 2 0 0 ]
0 0 𝑠+3 0
0 0 0 𝑠+4
(𝑠 + 2)(𝑠 + 3)(𝑠 + 4) 0 0 0
0 (𝑠 + 1)(𝑠 + 3)(𝑠 + 4) 0 0
𝑎𝑑𝑗[𝑠𝐼 ̅ − 𝐴̅] = [ 0 ]
0 0 (𝑠 + 1)(𝑠 + 2)(𝑠 + 4)
0 0 0 (𝑠 + 1)(𝑠 + 2)(𝑠 + 3)
1
0 0 0
𝑠+1
1 0 0
[𝑠𝐼 ̅ − 𝐴̅]−1 = 0 𝑠+2 1 0
0 0 𝑠+3 1
[0 0 0 𝑠+4]
𝑌̅ (𝑠)
𝐺 (𝑠 ) = ̅ (𝑠)
= 𝐶̅ (𝑠𝐼 ̅ − 𝐴̅)−1 𝐵̅ + 𝐷
̅
𝑈
1
𝑠+1
0 0 0 1
1 0 0 0
= [1 1 0 0] 0 𝑠+2 1 0 [1] + [0]
0 0 𝑠+3 1
[0 0 𝑠+4] 0
0
1
= 𝑠+1
This is a non-minimal system as it has 1 pole and 4 eigenvalues. The state space of a minimal system consists entirely of an
observable and controllable subspace (I).
Definitions
Controllability
A system is said to be controllable if and only if it is possible, by means of the input, to transfer the system form any initial
state 𝑥̅ (𝑡) = 𝑥̅ 𝑖 to any other state 𝑥̅ (𝑇) = 𝑥̅ 𝑇 in a finite time 𝑇 − 𝑡 ≥ 0.
Observability
An unforced system is said to be observable if and only if it is possible to determine any (arbitrarily initial) state 𝑥̅ (𝑡) = 𝑥̅ 𝑡 by
using only a finite record, 𝑦(𝜏) for 𝑡 ≤ 𝜏 ≤ −𝑡.
Tests
The controllability and observability of a system can be determined by the following two tests:
Controllability
The time-invariant system 𝑥ሶ ̅ = 𝐴̅𝑥̅ + 𝐵̅ 𝑢̅ is controllable if and only if the rank 𝑟(𝑄̅) of the controllability test matrix 𝑄̅ is equal
to n the order of the system.
If the test matrix is not of full rank, then the system is only partially state controllable. The rank deficiency of the test matrix is
equal to the number of uncontrollable modes in the system.
10
Observability
The (unforced) time invariant system 𝑥ሶ ̅ = 𝐴̅𝑥̅ with observable vector 𝑦̅(𝑡) = 𝐶̅ 𝑥̅ is observable if and only if the rank 𝑟(𝑁
̅ ) of
̅
the observability matrix 𝑁 is equal to n, the order of the system.
𝐶̅
𝐶̅ 𝐴̅
𝑄̅ = 𝐶̅ 𝐴2̅
⋯
[𝐶̅ 𝐴𝑛̅ −1 ]
If the test matrix is not of full rank, then the system is only partially state observable. The rank deficiency of the test matrix is
equal to the number of unobservable modes in the system.
Example 2.1
Consider the state space system
1.5 1] 1
𝐴̅ = [ 𝐵̅ = [ ], 𝐶̅ = [ 2 1 ] , ̅ = [ 0]
𝐷
1 0 0
1 1.5]
Controllability matrix: 𝑄̅ = [𝐵̅ 𝐴̅𝐵̅] = [
0 1
𝑑𝑒𝑡𝑄̅ = 1, and hence 𝑟(𝑄̅) = 2
Example2.2
Consider the state space system
−3 −4 ̅ 4
𝐴̅ = [ ] 𝐵 = [ ], 𝐶̅ = [−1 ̅ = [ 0]
−1], 𝐷
−1 0 1
The system is not observable. The rank deficiency is one; therefore, there is one unobservable mode.
The reason for this (possibly pole-zero cancellation) may be checked by deriving the LTF from the state space model.
𝑌̅ (𝑠)
𝐺 (𝑠 ) = ̅ (𝑠)
= 𝐶̅ (𝑠𝐼 ̅ − 𝐴̅)−1 𝐵̅ + 𝐷
̅
𝑈
The reason that the system is not completely observable is that if some disturbance signal excites the mode et arising from the
denominator term (s − 1), the dynamic effects will be cancelled out by the simultaneous response of the numerator term
(s − 1), thus leaving no visible effect at the output.
Checking the controllability of this system shows that it is not completely controllable. That is, one of the modes cannot be
excited form the input.
Stabilizability
If a system is completely controllable then any unstable eigenvalues may be moved to stable locations by state feedback. If the
system fails the controllability test, what then?
Tests (using the sensitivity matrix) may be carried out to find those modes that are uncontrollable. If the unstable pole that
needs to be moved is controllable, then the system is said to be stabilizable, and the design can proceed to move that pole.
The sensitivity matrix is formed corresponding to each eigenvalue of the system in turn (𝑖 = 1 𝑡𝑜 𝑛 for an nth-order system):
𝑆(𝜆𝑖 ) = [𝐴̅ − 𝜆𝑖 𝐼 ̅ 𝐵̅]
11
If this test is carried out on the eigenvalues of example 2.2, it is found that the unstable mode corresponding to 𝜆 = +1 is also
uncontrollable. Hence the system cannot be stabilized. The state space model must then be converted to some other form as
will discussed later.
Realization
We will start by looking at two simple realization methods namely the controllable canonical form and the observable canonical
form.
The 𝐴̅ matrix is in a canonical form, sometimes known as the companion form or the controllable canonical form as it satisfies
the controllability test matrix; that is, all the states can be modified by the control input.
y
+ + +
+ + +
bn bn b1 b0
u + 𝑥𝑛−1 𝑥2 𝑥1
𝑥𝑛 ∫ ∫
∫ ∫
+
an an a1 a0
+ + + +
+ +
Figure 2.2
Note that there is no unique state-space model of any given system. The canonical form given above has some mathematical
advantages. However, there is actually an infinite number of other possible state-space models of the same system. These are
obtained by choosing different sets of states, of which some are even more mathematically convenient (for example the 𝐴̅
matrix can be made diagonal by suitable choice of the states, thus making certain types of analysis relatively easy).
The main advantage of the controllable canonical form over other forms is that the differential equation describing the system
can be converted into this canonical form by simple substitution of the differential equation coefficients into the matrices 𝐴̅
and 𝐶̅ . The conversion can be performed by inspection, with none of the intervening mathematics being necessary.
Note
Some texts reverse the state numbering to give the following equivalent state matrices:
12
−𝑎𝑛−1 … −𝑎2 −𝑎1 −𝑎0 1
1 … 0 0 0 0
𝐴̅ = ⋮ ⋱ ⋮ ⋮ ⋮ , 𝐵̅ = ⋮ , 𝐶̅ = [𝑏𝑛−1 … 𝑏2 𝑏1 𝑏0 ]
0 … 1 0 0 0
[ 0 … 0 1 0 ] [ 0]
b0 b1 bn
+ +
+ 𝑥1 𝑥2 𝑥𝑛−1 𝑥𝑛 y
∫ ∫ ∫ ∫
+ + +
- - -
an an-1 a1
Figure 2.3
Note
Some texts reverse the state numbering to give the following equivalent state matrices:
−𝑎𝑛−1 1 … 0 0 𝑏𝑛−1
⋮ ⋮ ⋱ ⋮ ⋮ ⋮
𝐴̅ = −𝑎2 0 … 1 0, 𝐵̅ = 𝑏2 , 𝐶̅ = [ 1 … 0 0 0]
−𝑎1 0 … 0 1 𝑏1
[ −𝑎0 0 … 0 0] [ 𝑏0 ]
Notes
1. The controllable and observable canonical forms are easy to determine but are rather sensitive to parameter variations and
are thus of limited direct practical use.
2. The approach adopted also determines the choice of state variables.
Diagonal forms
This is when the state space matrices are transformed to give a diagonal or near diagonal 𝐴̅. Diagonal and near diagonal methods
have a number of practical advantages over the controllable and observable canonical forms already considered, and
consequently the extra effort in determining them is usually justified. Such diagonal matrices play an important role in the
13
analysis of MIMO systems, where decoupling of the states becomes very crucial. Three distinct types of 𝐴̅ matrix will be
considered:
1. All eigenvalues are real and distinct.
2. Some eigenvalues are complex but all are distinct.
3. Some eigenvalues are repeated.
Example 2.3
Consider the following transfer function:
𝑌(𝑠) 3𝑠 2 +7𝑠+15
𝐺 (𝑠 ) = =
𝑈(𝑠) 𝑠 3 +7𝑠 2 +14𝑠+14
Defining:
1
𝑋1 (𝑠) = 𝑠+1 𝑈(𝑠)
1
𝑋2 (𝑠) = 𝑠+2 𝑈(𝑠)
1
𝑋3 (𝑠) = 𝑠+4 𝑈(𝑠)
The output is then given by: 𝑌(𝑠) = 11⁄3 𝑋1 (𝑠) − 13⁄2 𝑋2 (𝑠) + 35⁄6 𝑋2 (𝑠). Figure 2.4 gives the Laplace block diagram.
1 𝑋1 (𝑠)
11⁄
𝑠+1 3
+
𝑈(𝑠) 1 𝑋2 (𝑠) -
𝑌(𝑠)
13⁄
𝑠+2 2
+
1 𝑋3 (𝑠) 35⁄
𝑠+4 6
Figure 2.4
1
Considering 𝑋2 (𝑠) = 𝑠+2 𝑈(𝑠). The expression can be written in the form: (𝑠 + 2)𝑋2 (𝑠) = 𝑈(𝑠)
𝑠𝑋2 (𝑠) + 2𝑋2 (𝑠) = 𝑈(𝑠)
Example 2.4
1 1
Given the TF 𝐺 (𝑠) = (𝑠+1)(𝑠+2) = 𝑠 2+3𝑠+2
2 1 0 1
𝐵̅𝐷 = [ ][ ] = [ ]
1 1 1 1
1 −1
𝐶𝐷̅ = [1 0] [ ] = [1 −1]
−1 2
Example 2.5
6𝑠 2 +26𝑠+8 −2 4+𝑗 4−𝑗
Given the TF 𝐺 (𝑠) = (𝑠+2)(𝑠 2+2𝑠+10) = 𝑠+2 + 𝑠+1+𝑗3 + 𝑠+1−𝑗3
Defining:
1
𝑋1 (𝑠) = 𝑈(𝑠)
𝑠+2
1
𝑋2 (𝑠) = 𝑠+1+𝑗3 𝑈(𝑠)
1
𝑋3 (𝑠) = 𝑠+1−𝑗3 𝑈(𝑠)
15
Then 𝑌(𝑠) = −2𝑋1 (𝑠) − (4 + 𝑗)𝑋2 (𝑠) + (4 − 𝑗)𝑋2 (𝑠)
In matrix form:
−2 0 0 1
𝐴̅ = [ 0 −1 − 𝑗3 0 ], 𝐵̅ = [1], 𝐶̅ = [−2 (4 + 𝑗) (4 − 𝑗)]
0 0 −1 + 𝑗3 1
Example 2.6
(4𝑠+2)(𝑠 2 +4𝑠+8) −5 5 4𝑠+12
Given the TF 𝐺 (𝑠) = (𝑠+3)(𝑠 2+2𝑠+2)(𝑠+1) = 𝑠+1 + 𝑠+3 + 𝑠 2+2𝑠+2
Defining:
1
𝑋1 (𝑠) = 𝑠+1 𝑈(𝑠)
1
𝑋2 (𝑠) = 𝑠+3 𝑈(𝑠)
Repeated eigenvalues
Again, a near diagonal; form is used called a Jordan canonical form. The basic idea is shown by the following example.
Example 2.7
3𝑠 2 +17𝑠+8 −3 6 −7
Given the TF 𝐺 (𝑠) = (𝑠+4)(𝑠+2)2 = + + (𝑠+2)2
𝑠+4 𝑠+2
Defining:
1
𝑋1 (𝑠) = 𝑈(𝑠)
𝑠+4
1
𝑋3 (𝑠) = 𝑠+2 𝑈(𝑠)
1 1 1 1
𝑋2 (𝑠) = (𝑠+2)2 𝑈 (𝑠) = 𝑠+2 × 𝑠+2 𝑈(𝑠) = 𝑠+2 𝑋3 (𝑠)
−4 0 0 1
𝐴̅ = [ 0 −2 1 ], 𝐵̅ = [0], 𝐶̅ = [−3 −7 6]
0 0 −2 1
Jordan blocks
16
1 𝑋1 (𝑠)
3
𝑠+4
-
𝑈(𝑠) 𝑋3 (𝑠) + 𝑌(𝑠)
6
-
1 1 𝑋2 (𝑠)
7
𝑠+2 𝑠+2
Figure 2.5
Basically, the Jordan canonical form (JCF) is a block diagonal matrix in which each n x n diagonal block matrix corresponds
to n-multiplicity eigenvalues. The main diagonal elements of the JCF are the eigenvalues (with their multiplicities) and for the
repeated eigenvalues the entry above the main diagonal is a 1 and the rest are 0s.
The idea can be extended to three or more repeated poles and multiple Jordan blocks can be included in the same system as
shown in the following example.
Example 2.8
−5 4 6 9 8 7
𝐺 (𝑠) = 𝑠+2 + (𝑠+2)2 + 𝑠+3 + 𝑠−4 − (𝑠−4)2 + (𝑠−2)3
This gives the following Jordan canonical form of the state matrices:
−2 0 0 0 0 0 0
0 −2 0 0 0 0 1
0 0 −3 0 0 0
𝐴̅ = 0 0 0 4 1 0, 𝐵̅ = 1 , 𝐶̅ = [4 −5 6 7 −8 9]
0
0 0 0 0 4 1 0
[0 0 0 0 0 4] [ 1]
Minimal realizations
If a transfer function has no common roots between the numerator and the denominator then it is a minimal system and its state
space realization is a minimal realization. If, however, there are common roots then the realization is a non-minimal realization
as shown in the following example.
Example 2.9
Consider the following transfer function:
𝑠 3 +4𝑠 2 +7𝑠+6
𝐺 (𝑠) = 𝑠 4+5𝑠 3+10𝑠 2+11𝑠+3
This is a non-minimal realization as 𝐺(𝑠) has the common factor 𝑠 2 + 2𝑠 + 3 in the numerator and denominator.
(𝑠+2)(𝑠 2 +2𝑠+3)
𝐺 (𝑠) = (𝑠 2+3𝑠+1)(𝑠 2+2𝑠+3)
17
Notes
𝑁(𝑠)
1. For any minimal system 𝐺(𝑠) = 𝐷(𝑠), an infinite number of non-minimal systems can be obtained simply by multiplying
numerator and denominator by 𝑃(𝑠) thus:
𝑁(𝑠)𝑃(𝑠)
𝐺(𝑠) = 𝐷(𝑠)𝑃(𝑠)
2. The use of minimal systems rather than non-minimal systems is generally preferred.
Similarity transforms
The state space representation is not unique and in general infinitely many representations are possible. It can be shown that
each state representation can be related to another by a simple linear transformation thus:
𝑧̅ = 𝑃̅𝑥̅ or 𝑥̅ = 𝑃̅ −1 𝑧̅
𝑦̅ = 𝐶̅ 𝑥̅ + 𝐷̅ 𝑢̅
𝑦̅ = 𝐶̅ 𝑃̅−1 𝑧̅ + 𝐷̅ 𝑢̅
Thus:
𝐴̅𝑧 = 𝑃̅𝐴̅𝑃̅−1 , 𝐵̅𝑧 = 𝑃̅𝐵̅, 𝐶𝑧̅ = 𝐶̅ 𝑃̅ −1 , ̅𝑧 = 𝐷
𝐷 ̅
Notes
1. Some texts use the convention: 𝑥̅ = 𝑃̅𝑧̅. Thus 𝑃̅ and 𝑃̅−1 are reversed in the above equations.
2. An important characteristic of similarity transformations is that they do not alter the eigenvalues of the underlying system,
i.e. the eigenvalues of 𝐴̅ and 𝑃̅𝐴̅𝑃̅−1 are the same. The various state descriptions of the system that preserve the eigenvalues
are known as canonical forms.
A simple example of the use of the similarity transformation is given in the following example:
Example 2.10
1 1 0
Given the system 𝑥ሶ ̅ = [ ] 𝑥̅ + [ ] 𝑢 , 𝑦 = [1 0]𝑥̅ :
1 1 1
0.5 −0.5]
Consider the transformation matrix 𝑃̅ = [ :
0.5 0.5
1 1
The inverse is 𝑃̅−1 = [ ]
−1 1
0.5 −0.5] [1 1 1 1 0 0
𝐴̅𝑧 = 𝑃̅𝐴̅𝑃̅−1 = [ ][ ]=[ ]
0.5 0.5 1 1 −1 1 0 2
18
Chapter 2 Tutorial Questions
1. For each of the following systems determine the TF, draw the state block diagram and check for observability and
controllability. In all cases D 0 .
2
a) A 1 B 1 C 2 ; C/O
s 1
0 0 1 2
b) A B C 2 2 ; C/NO
1 1 0 s 1
3 1 2 2
c) A B C 1 0 ; NC/O
2 0 4 s 1
3 1 0 2
2
d) A 2 0 0 B 4 C 1 0 0 ; NC/NO
s 1
1 0 0 1
4. Find the controllable canonical form realizations of the following, and draw the state block diagrams:
2 s 10 3s 2 s 2 2
(a) G ( s ) (b) G(s) (c) G ( s )
s2 2 s 3 2s 2 1 s s 1
3
5. Find the observable canonical form realizations of the following, and draw the state block diagrams:
5 3 4s 3 2s 1
(a) G ( s ) (b) G ( s ) (c) G( s)
2s 4s 3
2
2s 1
5
2s 3 3s 2 2
6. Using partial fractions determine the diagonal form of state matrices for the following:
7 s 4 2 s 2 3s 7 10
(a) G ( s) (b) G ( s ) (c) G ( s )
s 2 8s 12 ( s 2)(s 8)(s 5) s 3 8s 2 15s
7. For each of the following determine both diagonal form (complex) realizations and near diagonal form (real)
realizations:
4s 4
(a) G ( s) (b) G ( s)
s 2 2s 7 ( s 2)(s 2 2s 17)
8. For the following transfer functions determine the Jordan canonical form state matrices:
3s 1 s 3 4s 2 s 2
(a) G ( s) (b) G ( s )
s 4s 4
2
( s 2)(s 3) 3
19
9. Transform the controllable canonical form state matrices from Q4 using the similarity transformation matrix:
0 0 0 1
0 0 1 0
P
0 1 0 0
1 0 0 0
10. Find the controllable canonical minimal realizations of the following transfer function:
s 2 1
G( s)
s3 s2 s 3
20