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 />
Let us set x0 1 = z(t − d), <strong>and</strong> x02 = x(t − d) then Lemma 33 gives:<br />
�z(t − d) − x(t − d)� 2 ≤ 1<br />
λ0 � t<br />
�y (τ) − y (t − d, z (t − d) , τ)�<br />
d t−d<br />
2 dτ,<br />
or, equivalently,<br />
3.3 Innovation<br />
�z(t − d) − x(t − d)� 2 ≤ 1<br />
λ0 Id(t).<br />
d<br />
This is to say, that up to a multiplicative constant, innovation at time t upper bounds the<br />
estimation error at time t − d.<br />
Remark 34<br />
One could think that the adaptation scheme is likely to react with a delay time d when<br />
the estimation error is large. However, as is explained in Chapter 4, Remark 45, it may not<br />
always be the case in practice.<br />
Proof.<br />
Let x1 (t) =x x 0 1 ,u (t) <strong>and</strong> x2 (t) =x x 0 2 ,u (t) be the solutions of (3.2) with xi (0) = x 0 i ,<br />
i =1, 2. For any a ∈ [0, 1] ,<br />
b (ax2 + (1 − a) x1,u)<br />
� a ∂<br />
= b (x1,u)+<br />
0 ∂α b (α x2 + (1 − α) x1,u) dα<br />
� a ∂<br />
= b (x1,u)+<br />
∂x b (α x2 + (1 − α) x1,u) dα (x2 − x1) .<br />
Hence for a =1<br />
b (x2,u) − b (x1,u) =<br />
�� 1 ∂b<br />
0 ∂x (α x2<br />
=<br />
�<br />
+ (1 − α) x1,u) dα (x2 − x1)<br />
B (t)(x2 − x1) ,<br />
where B (t) = (bi,j) (i,j)∈{1,..,n} is a lower triangular matrix since<br />
Set ε = x1 − x2, <strong>and</strong> consider the system:<br />
⎧<br />
⎨<br />
⎩<br />
0<br />
b (x, u) =(b (x1,u) ,b(x1,x2,u) , . . . , b (x, u)) ′<br />
.<br />
˙ε = A (u) x1 + b (x1,u) − A (u) x2 − b (x2,u)<br />
= [A (u)+B (t)] ε<br />
yɛ = C (u) ε = a1 (u) ε1.<br />
It is uniformly observable5 as a result of the structure of B (t). Let us consider Ψ( t), the<br />
resolvent of the system, <strong>and</strong> the Gramm observability matrix Gd:<br />
� d<br />
Gd = Ψ (v) ′ C ′ CΨ (v) dv.<br />
5<br />
See [57] for instance, or compute the observability matrix<br />
�<br />
φO = C ′<br />
|(CA) ′<br />
| . . . |(CA n ) ′�′<br />
,<br />
<strong>and</strong> check the full rank condition for all inputs.<br />
0<br />
42