Lecture 13
Lecture 13
Lecture 13
Introduction
The white box view of dynamical systems
K. J. strm
1. Introduction
2. The Observer Problem
3. Combining with State Feedback
4. Integral action
6. Summary
dy
CA
dt
x
=
W
x
=
.
o
n1
n1
d y
CA
dtn1
Wo observability matrix
y
dy
dt
1
x = Wo
..
dn1 y
dtn1
An Other Attempt
The problem can be solved if the observability matrix has full
rank. Can we avoid differentiations?
dx
= Ax + Bu
dt
y = Cx
Introduce error x = x x
d x
= A( x x ) = A x
dt
An Alternative Solution
dx
= Ax + Bu
dt
y = Cx
dt
This method has the drawback that it only uses the output and
that it is differentiated many times which is highly sensitive to
noise.
By integrating the input
d x
= A x + Bu
dt
d x
= A x + Bu + K ( y C x )
dt
Introduce error x = x x
d x
= A( x x ) K ( y C x ) = ( A K C) x
dt
This method has the drawback that it only uses the input and
that it only works if the matrix A is stable.
Determine matrix K so that the matrix A K C has all its eigenvalues in the left half plane, then x will go to zero!
0
x+
u
1
A KC =
0)x
Observability matrix
Wo =
The matrix
k1
k2
k2
Choosing
k1 = 2 o
k2 = 2o
Input-output relations
X 1 ( s) =
X 2 ( s) =
U ( s) +
2
s2 + 2 o s + o
s2
2 o s + 2o
Y ( s)
s2 + 2 o s + 2o
s + 2 o
2o s
U ( s) + 2
Y ( s)
2
+ 2 o s + o
s + 2 o s + 2o
dy
dt
a1
1
a2
dz
..
.
dt
an1 0
an
0
y = 1 0 0...
...
0
1
0
0
0 z
1
0
0 ... 0
a1
1
0 ... 0
a
1
.
.
.
0
a
2
1
1
,
Wo =
.
.
...
0
b1
0
b2
..
z
+
0
bn
1
W
o
a1
a2
=
..
an
a1
...
...
...
an1
...
0
1
a1
1
a2
dz
..
.
dt
an1 0
an
0
y = 1 0 0...
...
0
0
0 z
0
b1
0
b2
..
z
+
0
bn
a
b
k
1 0 ... 0
1
1
1
a
0
1
0
b
k
2
2
2
d z
.
.
.
( y z1 )
.
.
.
=
+
u
+
z
.
.
.
dt
1
bn1
kn1
an1 0 0
an
bn
kn
Summary
Output Feedback
If the system
dx
= Ax
dt
y = Cx
dx
= Ax + Bu
dt
y = Cx
u = L x
d x
= A x + Bu + K ( y C x )
dt
dx
= Ax + Bu = Ax B L x = Ax B L( x x ) = ( A B L) x + B L x
dt
d x
= ( A K C) x
dt
K
u
Characteristic equation
det (sI A + B L) det (sI A + K C) = 0
c K. J. strm August, 2001
&
Process
dx
= Ax + Bu,
dt
y = Cx
Integral action
How to introduce reference signals
Modeling Disturbances
Classical disturbance models:
dv
= 0, v(t) = a = constant
Step
dt
d2 v
Ramp 2 = 0, v(t) = a + bt
dt
d2 v
Sinusoid 2 + 2 v = 0, v(t) = a sin( t + b)
dt
Disturbances!
How to describe disturbances
dx
= Ax + B (u + v)
dt
dv
=0
dt
0
dx
= 0 0 1 x + 1u
dt
0 0 0
0
y = Cx = ( 1 0 0 ) x
Observer
d
dt
x
v
=
x
0
= ( y C x )
+
u+
K
Kv
k1
k2
k3
A KC =
det (sI A + K C) =
1
s
k3
k3
To introduce the reference values we use the standard configuration with two degrees of freedom.
ur
0 !
k2
Model and
Feedforward
xm
Generator
u ff
u fb
Process
= s3 + k1 s2 + k2 s + k3
Observer
A nice structure
0
dx
= 0
dt
0
0 1 x + 1u
1
0
y = Cx = ( 1
0 1 0
0
k1
d x
= 0 0 1 x + 1 u + k2 ( y x 1 )
dt
k3
0 0 0
0
0
0
0)x
2m
R,
s2 + 2 m s + 2m
X m2 = s X m1
Hence
2m s2
R
+ 2 m s + 2m
2 m s + 2m
2
1 2
R = 2m R
sYm Ym
s + 2 m s + 2m
m
U m = P 1 M y R =
= 2m
s2
A Simulation
Summary
2 y och ym
0
0
10
15
0
0
1
state feedback
observer
reference trajectory generator
10
15
10
15
v and ve
0.5
0
0