WOW !! MUCH LOVE ! SO WORLD PEACE !
Fond bitcoin pour l'amélioration du site: 1memzGeKS7CB3ECNkzSn2qHwxU6NZoJ8o
  Dogecoin (tips/pourboires): DCLoo9Dd4qECqpMLurdgGnaoqbftj16Nvp


Home | Publier un mémoire | Une page au hasard

 > 

Conception d'un pro logiciel interactif sous r pour la simulation de processus de diffusion

( Télécharger le fichier original )
par Arsalane Chouaib GUIDOUM
Université des sciences et de technologie de Houari Boumedienne - Magister en mathématiques 2012
  

précédent sommaire suivant

Bitcoin is a swarm of cyber hornets serving the goddess of wisdom, feeding on the fire of truth, exponentially growing ever smarter, faster, and stronger behind a wall of encrypted energy

2.9 Pont brownien

Un pont brownien est un objet mathématique de la théorie des probabilités, également appelé mouvement brownien attaché ("tied down brownian motion" en anglais), mouvement brownien attaché en a ("brownian motion tied down at a" en anglais) ou mouvement brownien épinglé ("pinned brownian motion" en anglais).

Considérons un mouvement brownien (Wt)0=t=1 sur l'intervalle de temps [0,1]. On veut calculer la loi conditionnelle de Wt par rapport à la tribu engendrée par la variable aléatoire W1. Notons £(a,.) la loi conditionnelle de Wt sachant W1 = a. Si l'on modifie la loi du processus Wt, celui-ci changer de nom et de propriétés.

Définition 2.5 Le processus Wt sous la loi £(a,.), est appelé le pont brownien, et pont brownien standard si a = 0.

Définition 2.6 Un pont brownien Xt est un processus stochastique à temps continu dont la loi celle d'un processus de Wiener sachant l'événement W0 = W1 = a. Il s'agit d'un processus aléatoire gaussien, c'est à dire que la loi de probabilité de tout vecteur (Wt1,...,Wtn) conditionnellement à W1 = a est gaussienne. Il est alors caractérisé par sa moyenne et sa fonction de covariance, qui sont données par:

E(Wt|W1 = a) = at ?0 = t = 1

cov(Ws,Wt|W1 = a) = s(1-t) ?0 = s < t = 1

Observons que la moyenne dépend de a, mais pas la covariance. Et si a = 0 on a un pont brownien standard dont

E(Wt|W1 = 0) = 0 et cov(Ws,Wt|W1 = 0) = s(1-t) ?0 = s < t = 1

Son semi-groupe de transition est donné par la formule suivante, pour 0 = s < t < 1

( ~y - 1

1-s(x(1 -t) + a(t - s)))2 )

P(a) 1

s,t (x,dy) = q2ð(t-s)(1-t) exp dy

- 2(t-s)(1-t)

1-s 1-s

Le pont brownien peut être vu comme la définition de la loi d'un point Xè à la date è ? [s,t], intermédiaire sur la trajectoire d'un mouvement brownien dont Xs et Xt sont connus, un tel point suit une loi normale

( ~

Xs + è - s

t - s (Xt - Xs), (t - è)(è - s)

N t - s

Proposition 2.9 Si (Wt)0=t=1 est un mouvement brownien standard, alors on a

(1). Xt = Wt -tW1 est un pont brownien.

(2). Xt = (1 -t)Wt/1-t est un pont brownien.

Preuve Il suffit de vérifier que le processus est gaussien et de covariance égale à s(1-t) ?0 = s < t = 1. Vérifions (1), on a

E(XsXt) = E((Ws -sW1)(Wt -tW1))

= E(WsWt) -tE(WsW1) -sE(WtW1)+tsE(W21 )

= min(s,t) -ts = s(1 -t)

De même pour (2), on a

E(XsXt) = E((1 -s)(1 -t)Ws/1-sWt/1-t)

( s )

= (1 - s)(1 -t)min 1 - s, t

1 -t

= min(s(1 -t),t(1 -s)) = s(1 -t)

2.9.1 Construction par processus contraint

Soit (Wt)t=0 un mouvement brownien standard, a et b des réels quelconques. On définit le processus (Xè)s=è=t, comme la déformation de Wt, à partir de l'instant è = s forcée de passer par a à la date è = t. A tout instant, utilisant la nullité de l'espérance de dWt pour tout t, on écrit donc le drift du processus Xè comme la différence entre a et Xè, sur le temps restant avant t. On obtient donc une description du processus Xè comme :

a - Xè

dXè = t - è dè + dWè, avec Xs = b et Xt = a. (2.16)

L'équation différentielle stochastique (2.16) (voir le chapitre correspondant) admet pour solution:

è - s

Xt,a

s,b(è) = b +Wè-s - (Wt-s -a+b) (2.17)

t -s

Nous pouvons simuler une trajectoire d'un pont brownien directement à partir de l'équation (2.17). On considère la subdivision de l'intervalle de temps [0,T] suivante 0 = t1 <
·
·
· < tN < tN+1 = T, avec ti+1 --ti = Ät, pour i = 1 et W(0) = W(t1) = 0, on a l'algorithme suivant :

1. Générer une nouvelle variable aléatoire Z de distribution gaussienne N(0,1).

2. i = i+1.

V'

3. W(ti) = W(ti--1) + Z Ät.

4.X(ti)=b+Wti -- ti T (WT--a+b)

5. Si i N + 1, réitérez a l'étape 1.

Remarque 2.2 Si a = b = 0 on a un pont brownien standard.

La fonction BB (Brownian Bridge) permet de simuler un pont brownien sur l'intervalle de temps [t0,T] avec un pas Ät = (T --t0)/N, et la fonction BBF permet de simuler un flux de pont brownien.

R> b <- - 2; a <- 1;

R> BB(N = 1000, t0 = 0 , T = 1 , x0 = b, y = a)

R> BBF(N = 1000, M = 100, t0 = 0, T = 1, x0 = 0, y = 0)

FIGURE 2.15 - Trajectoire d'un pont brownien à partir de Xt0 = --2 et XT = 1.

FIGURE 2.16 - Flux de 100 trajectoires d'un pont brownien standard Xt0 = XT = 0.

2.9.2 Construction par le développement de Karhunen-Loève (D.K.L)

Pour un pont brownien standard {Xt,0 =t = 1} est un processus gaussien centré, à trajectoires continues (p.s), tel que X(0) = X(1) = 0. Le développement de Karhunen-Loève [15, 16] s'écrit pour une suite de variables aléatoires Zn de loi normale centrée réduite telles que E|Z2 n| = 1,

v

Xt = 2

8

?

n=1

sin(ðnt)

Zn

ðn

?t ? [0,1] (2.18)

Utilisant le code 5 pour simulée un pont brownien standard a partir de D.K.L. La figure 2.17 donne une représentation graphique d'une approximation d'un pont brownien standard par le D.K.L pour n = 10,n = 100 et n = 1000.

FIGURE 2.17 - Approximation d'un pont brownien standard par le D.K.L.

précédent sommaire suivant






Bitcoin is a swarm of cyber hornets serving the goddess of wisdom, feeding on the fire of truth, exponentially growing ever smarter, faster, and stronger behind a wall of encrypted energy








"Il ne faut pas de tout pour faire un monde. Il faut du bonheur et rien d'autre"   Paul Eluard