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 />

N(0,1)<br />

Normally distributed<br />

r<strong>and</strong>om generator<br />

zn Xn<br />

Xn = µXn−1 + σ � 1 − µ 2 zn<br />

C.5 Ornstein-Ulhenbeck Process<br />

Colored<br />

noise<br />

Figure C.1: Simulation of colored noise via block wise programing.<br />

The law of the r<strong>and</strong>om variable at time 0 is supposed the invariant probability associated<br />

to (C.1) such that Xt is a stationary process.<br />

Equation (C.1) can be rewritten:<br />

<strong>and</strong> then<br />

d � e ρt � ρt<br />

Xt = αe dBt,<br />

Xt = e −ρt<br />

� � t<br />

X0 + αe ρs �<br />

dBs .<br />

The stochastic process Xt is such that [28]:<br />

1. If X0 is a gaussian variable with zero mean <strong>and</strong> variance equals to α2<br />

2ρ , then X(t) is<br />

gaussian <strong>and</strong> stationary of covariance:<br />

2. X(t) is a markovian process.<br />

0<br />

E[X(t)X(s + t)] = α2<br />

2ρ e−rρ|s| .<br />

3. When X(0) = c, the stochastic law of X(t) is a normal law with mean e−ρtc, <strong>and</strong><br />

variance α2<br />

�<br />

2ρ 1 − e2ρt � .<br />

As explained in [59], Part 3, equation (3.15), for a sample time∆ t, the exact updating<br />

formulas of the Ornstein-Uhlenbeck process is given by:<br />

where:<br />

− µ = e −ρ∆t ,<br />

− γ 2 = � 1 − e −2ρ∆t� � α 2<br />

2µ<br />

�<br />

,<br />

X(t + ∆t) =X(t)µ + γzn<br />

− zn is the realization of a st<strong>and</strong>ard gaussian distribution (i.e. N(0, 1)).<br />

163<br />

(C.2)

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

Saved successfully!

Ooh no, something went wrong!