Adaptative high-gain extended Kalman filter and applications
Adaptative high-gain extended Kalman filter and applications
Adaptative high-gain extended Kalman filter and applications
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
tel-00559107, version 1 - 24 Jan 2011<br />
5.1 Multiple Inputs, Multiple Outputs Case<br />
The matrices A(u) <strong>and</strong> C(u) are given by:<br />
⎛<br />
⎞<br />
A1(u) 0 . . . 0<br />
⎜ 0 A2(u) . . . 0 ⎟<br />
A = ⎜<br />
⎝<br />
.<br />
. .. . ..<br />
⎟<br />
. ⎠<br />
0 . . . 0 Any(u)<br />
, Ai(u)<br />
⎛<br />
0 α<br />
⎜<br />
= ⎜<br />
⎝<br />
2 i . . . 0<br />
0 . ⎞<br />
. .<br />
. ..<br />
⎟<br />
0 ⎟<br />
.<br />
. .. . ..<br />
⎟ ni α ⎠<br />
i<br />
0 . . . 0 0<br />
<strong>and</strong> the matrix C(u) is the generalization:<br />
⎛<br />
C1(u)<br />
⎜ 0<br />
C = ⎜<br />
⎝ 0<br />
0<br />
C2(u)<br />
...<br />
...<br />
...<br />
. ..<br />
0<br />
0<br />
0<br />
⎞<br />
⎟<br />
⎠<br />
0 0 ... Cny(u)<br />
, Ci = � α1 i (u) 0 ... 0 � .<br />
Finally the vector field b(x, u) is defined as:<br />
⎛ ⎞<br />
b1(x, u)<br />
⎜<br />
b(x, u) = ⎜ b2(x, u) ⎟<br />
⎝ ... ⎠<br />
bny(x, u)<br />
, bi(x,<br />
⎛<br />
b<br />
⎜<br />
u) = ⎜<br />
⎝<br />
1 i (x1i ,u)<br />
b2 i (x1i ,x2i ,u)<br />
...<br />
b ni<br />
i (x,u)<br />
⎞<br />
⎟<br />
⎠ .<br />
Remark 47<br />
The very last component of each element bi(., .) of the vector field b is allowed to depend<br />
on the full state. As one can see, the linear part is a Brunovsky form of observability: this<br />
is clearly a generalization of system (3.2). Nevertheless not every observable system can be<br />
transformed into this form.<br />
5.1.2 Definition of the Observer<br />
In order to preserve the convergence result, apply a few modifications to the observer. There<br />
are mainly three points to consider<br />
− the definition of innovation is adapted, rather trivially, to a multidimensional output<br />
space;<br />
− the matrix∆ , generated with the parameter θ is not the generalization one would<br />
imagine initially;<br />
− <strong>and</strong> the definition of the matrix Rθ must remain compatible with the matrix∆.<br />
Let us first recall the equations of the observer:<br />
⎧<br />
⎨<br />
⎩<br />
dz<br />
dt = A(u)z + b(z, u) − S−1C ′<br />
dS<br />
dt<br />
(Cz − y(t))<br />
= −(A(u)+b∗ (z, u)) ′<br />
S − S(A (u)+b∗ (z, u)) + C ′<br />
R −1<br />
dθ<br />
dt<br />
θ C − SQθS<br />
= F(θ, Id (t))<br />
where Q <strong>and</strong> R are symmetric positive definite matrices of dimension (n × n) <strong>and</strong> (ny × ny)<br />
respectively. The innovation at time t is:<br />
� t<br />
Id(t) = �y(s) − y (t − d, z(t − d),s) � 2 ny R ds<br />
where:<br />
t−d<br />
R −1<br />
θ<br />
94