28.01.2013 Views

Adaptative high-gain extended Kalman filter and applications

Adaptative high-gain extended Kalman filter and applications

Adaptative high-gain extended Kalman filter and applications

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!