Lecture-06 Design of Control Systems in State Space Nust Masters 2022
Lecture-06 Design of Control Systems in State Space Nust Masters 2022
Lecture-6
Design of Control Systems in Sate Space
Dr. Imtiaz Hussain
Assistant Professor (Control Systems),
Department of Electronic and Power Engineering
PNEC-NUST, Karachi, Pakistan
email: [email protected]
Fall 2022
1
Lecture Outline
• Introduction
𝒙=
˙ 𝑨𝒙 +𝑩𝑢
𝑦 =𝑪𝒙
Topology of Pole Placement
• In a typical feedback control system, the output, , is fed
back to the summing junction.
𝒙=
˙ 𝑨𝒙 + 𝑩(𝑟 − 𝑲𝒙 )
𝒙=
˙ 𝑨𝒙 +𝑩𝑟 − 𝑩𝑲𝒙 𝑦 =𝑪𝒙
𝒙=(
˙ 𝑨− 𝑩𝑲 ) 𝒙 + 𝑩𝑟
Topology of Pole Placement
• For example consider a plant signal-flow graph in phase-
variable form
Topology of Pole Placement
• Each state variable is then fed back to the plant’s input, u,
through a gain, ki, as shown in Figure
Pole Placement
• We will limit our discussions to single-input, single-output
systems (i.e. we will assume that the control signal and
output signal to be scalars).
• We will also assume that the reference input is zero.
𝒙=(
˙ 𝑨 − 𝑩𝑲 )𝒙 + 𝑩𝑟
𝒙=(
˙ 𝑨− 𝑩𝑲 )𝒙 𝑢=− 𝑲𝒙
Pole Placement
𝒙=(
˙ 𝑨− 𝑩𝑲 )𝒙
• The stability and transient response characteristics are
determined by the eigenvalues of matrix .
• Ackermann’s formula
𝐶𝑀 = [ 𝐵 𝐵]
2 𝑛− 1
𝐴𝐵 𝐴 𝐵⋯ 𝐴
Pole Placement (Using Transformation Matrix P)
• Following are the steps to be followed in this particular method.
s+
Pole Placement (Using Transformation Matrix P)
• Following are the steps to be followed in this particular
method.
𝑲 =[ 𝛼 𝑛 − 𝑎 𝑛 𝛼 𝑛 −1 − 𝑎 𝑛 −1 ⋯ 𝛼2 − 𝑎 2 𝛼1 − 𝑎 1 ]
Pole Placement (Using Transformation Matrix P)
• Example-1: Consider the regulator system shown in following figure.
The plant is given by
[ ][ ][ ] [ ]
𝑥1 0 1 0 𝑥1 0
𝑥2 = 0 0 1 𝑥 2 + 0 𝑢 (𝑡 )
𝑥3 −1 −5 −6 𝑥3 1
• The system uses the state feedback control u=-Kx. The desired
eigenvalues are , ,. Determine the state feedback gain matrix K.
Pole Placement (Using Transformation Matrix P)
• Example-1: Step-1
[ ][ ][ ] [ ]
𝑥1 0 1 0 𝑥1 0
𝑥2 = 0 0 1 𝑥 2 + 0 𝑢 (𝑡 )
𝑥3 −1 −5 −6 𝑥3 1
[ ]
0 0 1
𝐶𝑀 = [ 𝐵 𝐴 𝐵 ]= 0
2
𝐴𝐵 1 −6
1 −6 31
[ ][ ][ ] [ ]
𝑥1 0 1 0 𝑥1 0
𝑥2 = 0 0 1 𝑥 2 + 0 𝑢 (𝑡 )
𝑥3 −1 −5 −6 𝑥3 1
[ ][ ][ ] [ ]
𝑥1 0 1 0 𝑥1 0
𝑥2 = 0 0 1 𝑥 2 + 0 𝑢 (𝑡 )
𝑥3 −1 −5 −6 𝑥3 1
• Hence
𝑎 1 =6 , 𝑎 2= 5 , 𝑎 3=1
Pole Placement (Using Transformation Matrix P)
• Example-1: Step-4
• Hence
𝛼1 =14 , 𝛼 2=60 , 𝛼 3=200
Pole Placement (Using Transformation Matrix P)
• Example-1: Step-4
𝑲 =[ 𝛼 3 − 𝑎 3 𝛼 2 − 𝑎2 𝛼1 − 𝑎 1 ]
𝑲 =[ 199 5 5 8]
Pole Placement (Using Transformation Matrix P)
[ ][ ][ ] [ ]
𝑥1 0 1 0 𝑥1 0
𝑥2 = 0 0 1 𝑥 2 + 0 𝑢 (𝑡 )
𝑥3 −1 −5 −6 𝑥3 1
+¿ 𝑥˙ 3 𝑥3 𝑥˙ 2 𝑥2 𝑥˙ 1 𝑥1
𝑢(𝑡)
+¿
∫ ∫ ∫
+¿ -6
+¿
+¿
-5
+¿
+¿ -1
[]
𝑥1
𝑢=− 𝑲 𝑥 𝑲 =[ 199 5 5 8] 𝑢= − [ 199 55 8 ] 𝑥2
𝑥3
199
+¿
+¿
55
+¿
+¿
8
+¿
𝑢(𝑡) +¿ 𝑥˙ 3 𝑥3 𝑥˙ 2 𝑥2 𝑥˙ 1 𝑥1
-1 ∫ ∫ ∫
+¿
+¿ -6
+¿
+¿
-5
+¿
+¿ -5
Pole Placement (Direct Substitution Method)
• Following are the steps to be followed in this particular
method.
𝐶𝑀 = [ 𝐵 𝐵]
2 𝑛− 1
𝐴𝐵 𝐴 𝐵⋯ 𝐴
Pole Placement (Direct Substitution Method)
• Following are the steps to be followed in this particular
method.
[ ][ ][ ] [ ]
𝑥1 0 1 0 𝑥1 0
𝑥2 = 0 0 1 𝑥 2 + 0 𝑢 (𝑡 )
𝑥3 −1 −5 −6 𝑥3 1
• The system uses the state feedback control u=-Kx. The desired
eigenvalues are , ,. Determine the state feedback gain matrix K.
Pole Placement (Using Transformation Matrix P)
• Example-1: Step-1
[ ][ ][ ] [ ]
𝑥1 0 1 0 𝑥1 0
𝑥2 = 0 0 1 𝑥 2 + 0 𝑢 (𝑡 )
𝑥3 −1 −5 −6 𝑥3 1
[ ]
0 0 1
𝐶𝑀 = [ 𝐵 𝐴 𝐵 ]= 0
2
𝐴𝐵 1 −6
1 −6 31
|[ ][ ][] |
𝑠 0 0 0 1 0 0
|𝑠𝐼 − 𝐴 + 𝐵𝐾 |= 0 𝑠 0 − 0 0 1 + 0 [ 𝑘1 𝑘2 𝑘3 ]
0 0 𝑠 −1 −5 −6 1
𝐶𝑀 = [ 𝐵 𝐵]
2 𝑛− 1
𝐴𝐵 𝐴 𝐵⋯ 𝐴
Pole Placement (Ackermann’s Formula)
• Following are the steps to be followed in this particular
method.
∅ ( 𝐴 ) = 𝐴 𝑛 +𝛼1 𝐴 𝑛 − 1 +⋯ + 𝛼𝑛 − 1 𝐴+ 𝛼𝑛 𝐼
Pole Placement (Ackermann’s Formula)
• Example-1: Consider the regulator system shown in following figure.
The plant is given by
[ ][ ][ ] [ ]
𝑥1 0 1 0 𝑥1 0
𝑥2 = 0 0 1 𝑥 2 + 0 𝑢 (𝑡 )
𝑥3 −1 −5 −6 𝑥3 1
• The system uses the state feedback control u=-Kx. The desired
eigenvalues are , ,. Determine the state feedback gain matrix K.
Pole Placement (Using Transformation Matrix P)
• Example-1: Step-1
[ ][ ][ ] [ ]
𝑥1 0 1 0 𝑥1 0
𝑥2 = 0 0 1 𝑥 2 + 0 𝑢 (𝑡 )
𝑥3 −1 −5 −6 𝑥3 1
[ ]
0 0 1
𝐶𝑀 = [ 𝐵 𝐴 𝐵 ]= 0
2
𝐴𝐵 1 −6
1 −6 31
∅ ( 𝐴 ) = 𝐴 3 + 𝛼 1 𝐴 2 +𝛼 2 𝐴 + 𝛼 3 𝐼
[ ][ ][ ] [ ]
𝑥1 0 1 0 𝑥1 0
𝑥2 = 0 0 1 𝑥 2 + 0 𝑢 (𝑡 )
𝑥3 −1 −5 −6 𝑥3 1
∅ ( 𝐴 ) = 𝐴 3 +14 𝐴 2 +6 0 𝐴 +200 𝐼
[ ] [ ] [ ] [ ]
3 2
0 1 0 0 1 0 0 1 0 1 0 0
∅ ( 𝐴 )= 0 0 1 +14 0 0 1 +6 0 0 0 1 +2 00 0 1 0
−1 −5 −6 −1 −5 −6 −1 −5 −6 0 0 1
[ ]
1 99 55 8
∅ ( 𝐴 )= −8 1 59 7
−7 − 34 1 17
Pole Placement (Ackermann’s Formula)
[ ] [ ]
0 0 1 1 99 55 8
[𝐵 𝐴𝐵 𝐴 𝐵 ]= 0
2
1 −6 ∅ ( 𝐴 )= −8 1 59 7
1 −6 31 −7 −34 1 17
2 −1
𝐾 =[ 0 0 1] [ 𝐵 𝐴𝐵 𝐴 ] ∅( 𝐴 )
[ ][ ]
−1
0 0 1 199 55 8
𝐾 =[ 0 0 1] 0 1 −6 −8 159 7
1 −6 31 −7 − 34 117
𝐾 = [ 199 55 8 ]
Pole Placement
• Example-2: Consider the regulator system shown in following figure. The
plant is given by
[ ][ ][ ] [ ]
𝑥1 1 2 1 𝑥1 1
𝑥2 = 0 1 3 𝑥2 + 0 𝑢(𝑡 )
𝑥3 1 1 1 𝑥3 1
• Determine the state feedback gain for each state variable to place the poles
at -1+j, -1-j,-3. (Apply all methods)
Observer Based Control
Introduction
• In the pole-placement approach to the design of control systems,
we assumed that all state variables are available for feedback.
𝒙 ˇ˙ = 𝑨𝒙 + 𝑩 𝑢 − 𝑨 𝒙ˇ − 𝑩 𝑢− 𝑲 𝑒 ( 𝑪𝒙 − 𝑪 𝒙
˙ −𝒙 ˇ)
Full Order State Observer
𝒙 ˇ˙ = 𝑨𝒙 + 𝑩 𝑢 − 𝑨 𝒙ˇ − 𝑩 𝑢− 𝑲 𝑒 ( 𝑪𝒙 − 𝑪 𝒙
˙ −𝒙 ˇ)
• Simplifications in above equation yields
𝒙 ˇ˙ = 𝑨( 𝒙 − 𝒙
˙ −𝒙 ˇ )− 𝑲 𝑒 𝑪 ( 𝒙 − 𝒙
ˇ) (3)
𝒆=𝒙 − 𝒙ˇ
• Equation (3) can now be written as
𝒆=
˙ 𝑨 𝒆 − 𝑲 𝑒 𝑪𝒆
˙
𝒆=( 𝑨 − 𝑲 𝑒 𝑪) 𝒆
Full Order State Observer
˙
𝒆=( 𝑨 − 𝑲 𝒆 𝑪 )𝒆
• From above equation we see that the dynamic behavior of the
error vector is determined by the eigenvalues of matrix A-KeC.
𝑣=− 𝑲𝒛
Duality Property
• If the dual system is completely state controllable, then the
state feedback gain matrix K can be determined such that
matrix A*-C*K will yield a set of the desired eigenvalues.
𝑲 ∗= 𝑲 𝑒
• Since 𝑲 𝑒=𝑲 ∗
[ ]
𝛼𝑛 − 𝑎𝑛
𝛼𝑛 − 1 − 𝑎𝑛 −1
∗
𝑲 𝑒= 𝑲 = ⋮
𝛼 2 − 𝑎2
𝛼 1 − 𝑎1
Observer Gain Matrix
• Direct Substitution Method
[ ]
𝑘1
𝑘2
𝑲 𝒆 = ⋮
𝑘3
𝑘𝑛
Observer Gain Matrix
• Ackermann’s Formula
−1
𝐾 =[ 0 1][ 𝐵 𝐵]
2 𝑛− 1
0 ⋯0 𝐴𝐵 𝐴 𝐵⋯ 𝐴 ∅( 𝐴 )
• For the dual system
𝒛˙ = 𝑨∗ 𝒛 +𝑪 ∗ 𝑣
𝑛= 𝑩 ∗ 𝒛
∗ −1
𝑲 =[ 0 1 ] [𝐶 ]
∗ ∗ ∗ ∗ 2 ∗ ∗ 𝑛 −1 ∗
0 ⋯0 𝐴 𝐶 (𝐴 ) 𝐶 ⋯ (𝐴 ) 𝐶 ∅(𝐴 )
• Since
𝑲 𝑒= 𝑲 = [ 0
∗
{ 0 ⋯0 1] [ 𝐶
∗ ∗
𝐴 𝐶
∗ ∗ 2
(𝐴 ) 𝐶 ⋯
∗ ∗ 𝑛 −1
(𝐴 ) 𝐶
∗ −1
] ∗
∅(𝐴 ) ∗}
Observer Gain Matrix
{
𝑲 𝑒= 𝑲 ∗ = [ 0 0 ⋯0 1] [ 𝐶∗ 𝐴∗ 𝐶 ∗ ( 𝐴∗ )2 𝐶 ∗ ⋯ ( 𝐴 ∗ )𝑛 −1 𝐶 ∗ ]
−1
}
∅ ( 𝐴∗ ) ∗
• Simplifying it further
{[ 𝐶 }
∗
∗ −1
∗ ∗
𝐾 𝑒= ∅ ( 𝐴 )
∗ ∗
𝐴 𝐶
∗ ∗
(𝐴 ) 𝐶
2 ∗
⋯ (𝐴 )
∗ 𝑛− 1
𝐶 ] [0 0 ⋯0 1]
∗
[ ][]
−1
𝐶 0
𝐶𝐴 0
𝐾 𝑒 = ∅( 𝐴 ) ⋮ ⋮
𝑛 −2
𝐶𝐴 0
𝐶𝐴
𝑛 −1
1
Observer Gain Matrix
• The observer poles must be two to five times faster than the
controller poles to make sure the observation error (estimation
error) converges to zero quickly.
[ ][
𝑥
˙1
𝑥
˙2
=
0
1
20.6
0 ][ ] [ ]
𝑥1
𝑥2
+
0
1
𝑢(𝑡 )
𝑦 =[ 0 1 ]
[ ]
𝑥1
𝑥2
[ ][
𝑥
˙1
𝑥
˙2
=
0
1
20.6
0 ][ ] [ ]
𝑥1
𝑥2
+
0
1
𝑢(𝑡 )
𝑦 =[ 0 1 ]
[ ]
𝑥1
𝑥2
[ ][
𝑥
˙1
𝑥
˙2
=
0
1
20.6
0 ][ ] [ ]
𝑥1
𝑥2
+
0
1
𝑢(𝑡 )
𝑦 =[ 0 1 ]
[ ]
𝑥1
𝑥2
[ ][
𝑥
˙1
𝑥
˙2
=
0
1
20.6
0 ][ ] [ ]
𝑥1
𝑥2
+
0
1
𝑢(𝑡 )
𝑦 =[ 0 1 ]
[ ]
𝑥1
𝑥2
[ ][
𝑥
˙1
𝑥
˙2
=
0
1
20.6
0 ][ ] [ ]
𝑥1
𝑥2
+
0
1
𝑢(𝑡 )
𝑦 =[ 0 1 ]
[ ]
𝑥1
𝑥2
• We have
𝛼1 =20 , 𝛼2 =100
Example-1 (Method-1)
• Observer gain matrix Ke can be calculated using following formula
𝑲 𝑒=
𝛼 2 − 𝑎2
𝛼 1 − 𝑎1[ ]
• Where
𝑲 𝑒=
[ 1 00 − (− 20.6)
20 − 0 ]
𝑲 𝑒= [ 1 20.6
20 ]
Example-1 (Method-2)
[ ][
𝑥
˙1
𝑥
˙2
=
0
1
20.6
0 ][ ] [ ]
𝑥1
𝑥2
+
0
1
𝑢(𝑡 )
𝑦 =[ 0 1 ]
[ ]
𝑥1
𝑥2
|𝑠 𝑰 − 𝑨+ 𝑲 𝑒 𝑪|= |[ 𝑠
0 𝑠] [
0 − 0
1 0 ][ ]
20.6 + 𝑘 𝑒 1 [ 0
𝑘𝑒 2 |
1]
¿ 𝑠 2 +𝑘𝑒 2 𝑠 −2 0 .6+ 𝑘𝑒 1
Example-1 (Method-2)
• The desired characteristic polynomial is
𝑲 𝑒=
[ 1 20.6
20 ]
Example-1 (Method-3)
[ ][
𝑥
˙1
𝑥
˙2
=
0
1
20.6
0 ][ ] [ ]
𝑥1
𝑥2
+
0
1
𝑢(𝑡 )
𝑦 =[ 0 1 ]
[ ]
𝑥1
𝑥2
• Where ∅ ( 𝑨 ) = 𝐴 2+ 𝛼1 𝐴 +𝛼 2 𝐼
𝛼1 =20 , 𝛼2 =100
∅ ( 𝑨 ) = 𝐴 2+ 20 𝐴+100 𝐼
Example-1 (Method-3)
∅ ( 𝑨 ) = 𝑨 2+20 𝑨 +100 𝑰
[ ] [ ] [ ]
2
0 20.6 0 20.6 1 0
∅ ( 𝑨) = + 20 +100
1 0 1 0 0 1
∅ ( 𝑨) =
[ 1 2 0 .6
20
412
12 0 .6 ]
Example-1 (Method-3)
• Using Ackermann’s formula
[ ] [ ]
−1
𝑪 0
𝑲 𝑒= ∅ ( 𝑨 )
𝑪𝑨 1
[ ][ ] []
−1
1 20.6 412 0 1 0
𝑲 𝑒=
20 120.6 1 0 1
𝑲 𝑒 =
[ 1 20.6
20 ]
Example-1
• We get the same Ke regardless of the method employed.
[ ]
ˇ˙ 𝟏
𝒙
ˇ˙ 𝟐
𝒙
= [
0
1
20.6
0 ][ ] [ ]
ˇ1
𝒙
ˇ2
𝒙
+
0
1
𝑢 (𝑡 )+ [
1 20.6
20 ]
(𝑦 −[ 0 1]
[ ]
ˇ1
𝒙
ˇ2
𝒙
)
[ ][
ˇ˙ 𝟏
𝒙
ˇ˙ 𝟐
𝒙
= 0
1
20.6
0 ][ ] [ ]
ˇ1
𝒙
ˇ2
𝒙 [ ] [
+ 0 𝑢 (𝑡 )+ 1 20.6 𝑦 − 1 20.6 [ 0
1 20 20 ] 1]
[ ]
ˇ1
𝒙
ˇ2
𝒙
[ ] ][ ] [ ]
˙ 𝟏
[ [ ]𝑦
ˇ
𝒙 0 − 100 ˇ 1
𝒙 0 1 20.6
= + 𝑢 (𝑡 )+
˙
ˇ
𝒙 1 − 20 ˇ 2
𝒙 1 20
𝟐
Example-2
• Design a regulator system for the following plant:
[ ][
𝑥
˙1
𝑥
˙2
=
0
1
20.6
0 ][ ] [ ]
𝑥1
𝑥2
+
0
1
𝑢(𝑡 )
𝑦 =[ 1 0 ]
[ ]
𝑥1
𝑥2
• The desired closed-loop poles for this system are at , . Compute the
state feedback gain matrix K to place the poles of the system at
desired location.
• Obtain the observer gain matrix Ke and draw a block diagram for the
observed-state feedback control system.
Introduction
• The state observers discussed in previous lecture was
designed to estimate all the state variables.
• Or
• The output equations for the full order and minimum order
observers are
Minimum Order State Observer
• List of Necessary Substitutions for Writing the Observer Equation for
the Minimum-Order State Observer
Table–1
Minimum Order State Observer
• The observer equation for the full-order observer is given by :
END OF LECTURE-06