Lecture 12
Lecture 12
Introduction
K. J. strm
1. Introduction
2. An Example
4. Integral Action
5. Summary
State Feedback
dx
= Ax + Bu
dt
y = Cx
Matlab
dx
= Ax + Bu = Ax + B ( Lx + lr r) = ( A B L) x + Blr r
dt
Ae = e
( A I)e = 0
det A I = 0
det(sI A + B L) = 0
Simple Examples
a22
b21
b22
c21
c22
c23
dx
=
dt
Choosing
y = Cx = ( 1
x+
0
1
0)x
u = l1 x1 l2 x2 + lr r
l1 = 2 0, l2 = 02
Control law
a21
Process model
Characteristic equation
s 2 + l1 s + l2 = 0
lr
R( s )
s2 + 2 0 s + 02
Y ( s) =
02
R( s )
s2 + 2 0 s + 02
s a22
Example
Controllability
dx
= Ax + Bu
dt
Find control signal that moves the system from x(0) = a to
x(t) = b. Hence
Z t
At
b = x(t) = e a +
eA(ts) Bu(s)ds
y = Cx = ( 1
x+
1
0
0) x
This implies
Characteristic equation
s ( s + l1 ) = s 2 + l1 s = 0
e As Bu(s)ds = a + be At
Examples
Controllability ...
Z
As
dx
=
dt
At
Bu(s)ds = a + be
n1
B . Equation
Wc = ( B , AB ) =
dx
=
dt
we have
x+
Wc = ( B , AB ) =
Wc = ( B AB . . . An1 B )
we have
k=1
, Controllable!
x+
1
0
Not Controllable!
A Special Case
The system
a1 a2 . . . an1 an
1
1
0
0
0
0
dz
0
1
0
0
0
=
z
+
dt
..
..
.
.
is controllable
1 a1 a21 a2 . . .
0
1
a1
...
0
0
1
...
Wc =
.
0
...
The system
a1 a2 . . . an1 an
1
1
0
0
0
0
dz
0
1
0
0
0
=
z
+
dt
..
..
.
.
A Special Case
a1 a2 . . . an1 an
1
1
0
0
0
0
dz
0
1
0
0
0
=
z
+
u
dt
.
.
.
.
.
.
a1 a2 . . . an1 an
1
1
0
0
0
0
dz
0
1
0
0
0
=
z
+
u
dt
.
.
.
.
.
.
a1 l1 a2 l2 . . . an ln
lr
1
0
0
0
dz
0
1
0
0
=
z
+
r
dt
.
.
.
.
.
0
u = l1 z1 l2 z2 . . . ln z n + lr r
l1 = p1 a1 ,
l2 = p2 a2 , . . . , ln = pn an
A system
dx
= Ax + Bu
dt
which is controllable can be transformed to
a1
dz
0
=
dt
The feedback
n1
B)
Change coordinates z = T x
dz
dx
+ Bu
=T
= T Ax + T Bu = T AT 1 z + T Bu = Az
dt
dt
a2
. . . an
...
...
...
z
+
u
.
0
p2 a2
...
pn an )
c = ( B A B . . . A n1 B ) = T ( B AB . . . An1 B ) = T Wc
W
c K. J. strm August, 2001
&
An Algorithm
dx
= Ax + Bu to controllable canonical
dt
dz
+ Bu
using the transformation z = T x = W
c W 1 x .
form
= Az
c
dt
Feedback for transformed system u = z where
= ( p1 a1 p2 a2 . . . pn an )
L
We have u = Lx = LT
L = (0
...
1)( B
AB
...
An1 B )1 P( A)
... 0 1)
Wc = ( B AB . . . An1 B )
en = ( 0
, hence
z = Lz
L = en Wc1 P( A)
=L
W
c W 1
L = LT
c
dx
0 1
= Ax + Bu =
x+
1 0
dt
The system is controllable. The control law
Using Matlab
0
1
u = l1 x1 l2 x2
l1 1
1 l1
l2
1
s + l2
x
= s2 + l2 s + l1 1
Process model
dx
= Ax + Bu
dt
y = Cx
x
I
=
+
u+
0
1
u = ( L lI )
x
I
+ lr r
CB
...
CA
Summary
A nice way to use process dynamics to predict
A nice alternative to pole placement
Actually works for systems with many inputs and outputs
Nice computational methods matrix calculations
Works for high order systems
The idea of augmentation
Integral action still ad hoc will be dealt with later
Reference values will be done later
The control problem is thus easy if all state variables are
measured
What to do if all states are not measured? (Stay tuned for
next lecture!)