📚 الرئيسية

مسابقة دكتوراه 2023Université Ahmed Draia d'Adrar — الموضوع 02

مسابقة تخصص · Probabilités & Statistiques · المدة: 2سا

Concours d'accès à la formation 3e cycle, filière Mathématiques, spécialité Processus aléatoire et Statistique, année universitaire 2022-2023, Université Ahmed Draia d'Adrar, Faculté des Sciences et de Technologie, Département de Mathématiques et Informatique, Épreuve : Processus stochastique et équations différentielles, 07 février 2023, durée 02h00.

التمرين 1

Exercice 1 — Processus gaussien construit à partir d'un mouvement brownien et d'une variable gaussienne

#stochastic-processes#brownian-motion#gaussian-process#conditional-expectation#covariance

Soient B=(Bt)tB=(B_t)_t un mouvement brownien sur (Ω,F,(Ft)t0,P)(\Omega,\mathcal F,(\mathcal F_t)_{t\ge 0},\mathbb P) et η\eta une variable aléatoire distribuée suivant une loi gaussienne N(m,ρ2)\mathcal N(m,\rho^2) indépendante de (Bt)t(B_t)_t. Soit

Yt=ηt+τBt,Gt=σ(Ys,st).Y_t=\eta t+\tau B_t,\qquad \mathcal G_t=\sigma(Y_s,\, s\le t).

  1. (1,5 pt) Calculer Cov(η,Yt)\mathrm{Cov}(\eta,Y_t) et Cov(Ys,Yt)\mathrm{Cov}(Y_s,Y_t).
  2. (1,5 pt) Montrer que (Yt)t(Y_t)_t est un processus gaussien.
  3. (2 pts) Montrer que, pour tout t0t\ge 0, il existe λ\lambda (dépendant de tt) tel que η=λYt+Z\eta=\lambda Y_t+Z, avec ZZ indépendant de Gt\mathcal G_t.
  4. (1 pt) Calculer E[ηGt]\mathbb E[\eta\mid \mathcal G_t] et la variance de ηE[ηGt]\eta-\mathbb E[\eta\mid \mathcal G_t], puis montrer que limt+E[ηGt]=ηp.s.\lim_{t\to+\infty}\mathbb E[\eta\mid \mathcal G_t]=\eta\quad \text{p.s.}
الحل

1. Covariances

On a

Yt=ηt+τBt.Y_t=\eta t+\tau B_t.

Comme η\eta est indépendante de BtB_t et E(Bt)=0\mathbb E(B_t)=0 :

Cov(η,Yt)=Cov(η,ηt+τBt)=tVar(η)+τCov(η,Bt)=tρ2.\mathrm{Cov}(\eta,Y_t)=\mathrm{Cov}(\eta,\eta t+\tau B_t)=t\,\mathrm{Var}(\eta)+\tau\,\mathrm{Cov}(\eta,B_t)=t\rho^2.

Donc

Cov(η,Yt)=tρ2.\boxed{\mathrm{Cov}(\eta,Y_t)=t\rho^2.}

Pour s,t0s,t\ge 0,

Cov(Ys,Yt)=Cov(ηs+τBs,ηt+τBt).\mathrm{Cov}(Y_s,Y_t)=\mathrm{Cov}(\eta s+\tau B_s,\eta t+\tau B_t).

En développant et en utilisant l'indépendance,

Cov(Ys,Yt)=stρ2+τ2Cov(Bs,Bt).\mathrm{Cov}(Y_s,Y_t)=st\,\rho^2+\tau^2\mathrm{Cov}(B_s,B_t).

Or pour un mouvement brownien,

Cov(Bs,Bt)=min(s,t).\mathrm{Cov}(B_s,B_t)=\min(s,t).

Ainsi,

Cov(Ys,Yt)=stρ2+τ2min(s,t).\boxed{\mathrm{Cov}(Y_s,Y_t)=st\rho^2+\tau^2\min(s,t).}

2. (Yt)t(Y_t)_t est un processus gaussien

Considérons t1,,tn0t_1,\dots,t_n\ge 0 et a1,,anRa_1,\dots,a_n\in\mathbb R. Alors

i=1naiYti=ηi=1naiti+τi=1naiBti.\sum_{i=1}^n a_iY_{t_i}=\eta\sum_{i=1}^n a_it_i+\tau\sum_{i=1}^na_iB_{t_i}.

Le vecteur (Bt1,,Btn)(B_{t_1},\dots,B_{t_n}) est gaussien, donc aiBti\sum a_iB_{t_i} est gaussienne. La variable η\eta est gaussienne et indépendante de cette combinaison linéaire. Une somme de variables gaussiennes indépendantes est gaussienne.

Donc toute combinaison linéaire finie de (Yti)(Y_{t_i}) est gaussienne, d'où

(Yt)t est un processus gaussien.\boxed{(Y_t)_t\text{ est un processus gaussien}.}

3. Décomposition orthogonale de η\eta

Comme (η,Yt)(\eta,Y_t) est un vecteur gaussien, on peut projeter η\eta sur YtY_t :

λ=Cov(η,Yt)Var(Yt).\lambda=\frac{\mathrm{Cov}(\eta,Y_t)}{\mathrm{Var}(Y_t)}.

Or

Var(Yt)=t2ρ2+τ2t.\mathrm{Var}(Y_t)=t^2\rho^2+\tau^2 t.

Donc

λ=tρ2t2ρ2+τ2t=ρ2tρ2+τ2.\lambda=\frac{t\rho^2}{t^2\rho^2+\tau^2 t}=\frac{\rho^2}{t\rho^2+\tau^2}.

Posons

Z=ηλYt.Z=\eta-\lambda Y_t.

Alors

Cov(Z,Yt)=Cov(η,Yt)λVar(Yt)=0.\mathrm{Cov}(Z,Y_t)=\mathrm{Cov}(\eta,Y_t)-\lambda\mathrm{Var}(Y_t)=0.

Comme le couple est gaussien, l'absence de covariance implique l'indépendance entre ZZ et YtY_t. Plus généralement, dans ce modèle gaussien linéaire, ZZ est indépendante de Gt\mathcal G_t.

Ainsi,

η=ρ2tρ2+τ2Yt+Z,Z ⁣ ⁣ ⁣Gt.\boxed{\eta=\frac{\rho^2}{t\rho^2+\tau^2}Y_t+Z,\qquad Z\perp\!\!\!\perp \mathcal G_t.}

4. Espérance conditionnelle et variance résiduelle

Par la propriété précédente,

E[ηGt]=λYt+E[ZGt]=λYt+E[Z].\mathbb E[\eta\mid\mathcal G_t]=\lambda Y_t+\mathbb E[Z\mid\mathcal G_t]=\lambda Y_t+\mathbb E[Z].

Comme E(η)=m\mathbb E(\eta)=m et E(Yt)=mt\mathbb E(Y_t)=mt, la formule de régression gaussienne donne plus directement

E[ηGt]=m+Cov(η,Yt)Var(Yt)(Ytmt).\mathbb E[\eta\mid\mathcal G_t]=m+\frac{\mathrm{Cov}(\eta,Y_t)}{\mathrm{Var}(Y_t)}(Y_t-mt).

Donc

E[ηGt]=m+ρ2tρ2+τ2(Ytmt).\boxed{\mathbb E[\eta\mid\mathcal G_t]=m+\frac{\rho^2}{t\rho^2+\tau^2}(Y_t-mt).}

La variance conditionnelle résiduelle vaut

Var(ηE[ηGt])=ρ2Cov(η,Yt)2Var(Yt).\mathrm{Var}(\eta-\mathbb E[\eta\mid\mathcal G_t])=\rho^2-\frac{\mathrm{Cov}(\eta,Y_t)^2}{\mathrm{Var}(Y_t)}.

Ainsi,

Var(ηE[ηGt])=ρ2t2ρ4t2ρ2+τ2t=ρ2τ2tρ2+τ2.\mathrm{Var}(\eta-\mathbb E[\eta\mid\mathcal G_t])=\rho^2-\frac{t^2\rho^4}{t^2\rho^2+\tau^2 t}=\rho^2\frac{\tau^2}{t\rho^2+\tau^2}.

Donc

Var(ηE[ηGt])=ρ2τ2tρ2+τ2.\boxed{\mathrm{Var}(\eta-\mathbb E[\eta\mid\mathcal G_t])=\frac{\rho^2\tau^2}{t\rho^2+\tau^2}.}

Cette quantité tend vers 00 quand t+t\to +\infty, donc

ηE[ηGt]0dans L2.\eta-\mathbb E[\eta\mid\mathcal G_t]\to 0\quad \text{dans }L^2.

Comme Gt\mathcal G_t est croissante et que E[ηGt]\mathbb E[\eta\mid\mathcal G_t] est une martingale bornée dans L2L^2, elle converge presque sûrement et dans L2L^2 vers E[ηG]\mathbb E[\eta\mid\mathcal G_\infty]. L'unicité de la limite L2L^2 donne

limt+E[ηGt]=ηp.s.\boxed{\lim_{t\to+\infty}\mathbb E[\eta\mid\mathcal G_t]=\eta\quad \text{p.s.}}

التمرين 2

Exercice 2 — Construction hilbertienne d'un mouvement brownien

#brownian-motion#hilbert-space#orthonormal-basis#gaussian-series#holder-regularity

Soient I=[0,T]I=[0,T] et (φn)n(\varphi_n)_n une base hilbertienne orthonormale de l'espace de Hilbert L2(I)L^2(I) des fonctions φ:IR\varphi:I\to\mathbb R de carré intégrable pour la mesure de Lebesgue sur II. On note

φ,ψ=0Tφ(t)ψ(t)dt\langle \varphi,\psi\rangle=\int_0^T \varphi(t)\psi(t)\,dt

le produit scalaire de L2(I)L^2(I). Soit (Xn)nN(X_n)_{n\in\mathbb N} une suite de variables aléatoires i.i.d. de loi N(0,1)\mathcal N(0,1) définie sur un espace de probabilités (Ω,F,P)(\Omega,\mathcal F,\mathbb P). Pour tout tIt\in I, on pose

Bt=nN1[0,t],φnXn.B_t=\sum_{n\in\mathbb N}\langle \mathbf 1_{[0,t]},\varphi_n\rangle X_n.

  1. (1,5 pt) Montrer que la série nN1[0,t],φn2\sum_{n\in\mathbb N}|\langle \mathbf 1_{[0,t]},\varphi_n\rangle|^2 est convergente.
  2. (1,5 pt) Montrer que le processus B=(Bt)tIB=(B_t)_{t\in I} est bien défini.
  3. (2 pts) Montrer que B=(Bt)tIB=(B_t)_{t\in I} est un mouvement brownien naturel sur (Ω,F,P)(\Omega,\mathcal F,\mathbb P).
  4. (2 pts) Montrer que B=(Bt)tIB=(B_t)_{t\in I} admet une modification dont les trajectoires sont localement α\alpha-höldériennes pour tout 0<α<120\lt \alpha\lt \frac12.
الحل

1. Convergence de la série des coefficients

Comme (φn)n(\varphi_n)_n est une base hilbertienne orthonormale de L2(I)L^2(I), l'identité de Parseval donne, pour tout tIt\in I,

nN1[0,t],φn2=1[0,t]L2(I)2.\sum_{n\in\mathbb N}|\langle \mathbf 1_{[0,t]},\varphi_n\rangle|^2=\|\mathbf 1_{[0,t]}\|_{L^2(I)}^2.

Or

1[0,t]L2(I)2=0T1[0,t](s)ds=t.\|\mathbf 1_{[0,t]}\|_{L^2(I)}^2=\int_0^T \mathbf 1_{[0,t]}(s)\,ds=t.

Donc

nN1[0,t],φn2=t<+.\boxed{\sum_{n\in\mathbb N}|\langle \mathbf 1_{[0,t]},\varphi_n\rangle|^2=t\lt +\infty.}

2. Bonne définition de BtB_t

Posons an(t)=1[0,t],φna_n(t)=\langle \mathbf 1_{[0,t]},\varphi_n\rangle. La série définissant BtB_t est

nNan(t)Xn.\sum_{n\in\mathbb N} a_n(t)X_n.

Comme les XnX_n sont centrées, indépendantes, de variance 11, on a

E[(n=1Nan(t)Xn)2]=n=1Nan(t)2.\mathbb E\left[\left(\sum_{n=1}^N a_n(t)X_n\right)^2\right]=\sum_{n=1}^N a_n(t)^2.

La série des carrés converge d'après la question précédente, donc les sommes partielles forment une suite de Cauchy dans L2(Ω)L^2(\Omega). Il existe donc une variable aléatoire carrée intégrable BtB_t telle que

n=1Nan(t)XnNL2Bt.\sum_{n=1}^N a_n(t)X_n\xrightarrow[N\to\infty]{L^2} B_t.

Ainsi BtB_t est bien définie pour tout tIt\in I.

Bt est bien deˊfinie pour tout tI.\boxed{B_t\text{ est bien définie pour tout }t\in I.}

3. BB est un mouvement brownien

3.1. Gaussianité

Chaque BtB_t est limite en L2L^2 de combinaisons linéaires de variables gaussiennes indépendantes, donc est gaussienne. De même, tout vecteur fini (Bt1,,Btk)(B_{t_1},\dots,B_{t_k}) est gaussien. Le processus est donc gaussien centré.

3.2. Covariance

Pour s,tIs,t\in I,

E[BsBt]=nN1[0,s],φn1[0,t],φn,\mathbb E[B_sB_t]=\sum_{n\in\mathbb N}\langle \mathbf 1_{[0,s]},\varphi_n\rangle\langle \mathbf 1_{[0,t]},\varphi_n\rangle,

car les XnX_n sont orthogonales dans L2(Ω)L^2(\Omega). Par Parseval,

E[BsBt]=1[0,s],1[0,t]=0T1[0,s](u)1[0,t](u)du=min(s,t).\mathbb E[B_sB_t]=\left\langle \mathbf 1_{[0,s]},\mathbf 1_{[0,t]}\right\rangle=\int_0^T \mathbf 1_{[0,s]}(u)\mathbf 1_{[0,t]}(u)\,du=\min(s,t).

Donc

E[BsBt]=min(s,t).\boxed{\mathbb E[B_sB_t]=\min(s,t).}

3.3. Accroissements indépendants

Pour 0t0<t1<<tkT0\le t_0\lt t_1\lt \cdots \lt t_k\le T,

BtiBti1=n1(ti1,ti],φnXn.B_{t_i}-B_{t_{i-1}}=\sum_n \langle \mathbf 1_{(t_{i-1},t_i]},\varphi_n\rangle X_n.

Deux accroissements disjoints correspondent à des indicatrices orthogonales dans L2(I)L^2(I), donc leurs covariances sont nulles. Comme le processus est gaussien, ils sont indépendants.

3.4. Loi des accroissements

Var(BtBs)=E[(BtBs)2]=ts,\mathrm{Var}(B_t-B_s)=\mathbb E[(B_t-B_s)^2]=t-s,

pour 0s<tT0\le s\lt t\le T. Donc

BtBsN(0,ts).B_t-B_s\sim \mathcal N(0,t-s).

Enfin B0=0B_0=0 car 1[0,0]=0\mathbf 1_{[0,0]}=0.

Tous les axiomes du mouvement brownien sont satisfaits :

B=(Bt)tI est un mouvement brownien sur [0,T].\boxed{B=(B_t)_{t\in I}\text{ est un mouvement brownien sur }[0,T].}

4. Modification höldérienne

Pour s,tIs,t\in I,

BtBsN(0,ts).B_t-B_s\sim \mathcal N(0,|t-s|).

Pour tout entier p1p\ge 1, il existe une constante Cp>0C_p\gt 0 telle que

E(BtBs2p)=Cptsp.\mathbb E\big(|B_t-B_s|^{2p}\big)=C_p |t-s|^p.

Choisissons p2p\ge 2. Alors

E(BtBs2p)Cpts1+(p1).\mathbb E\big(|B_t-B_s|^{2p}\big)\le C_p |t-s|^{1+(p-1)}.

On applique le critère de Kolmogorov : il existe une modification de BB dont les trajectoires sont localement α\alpha-höldériennes pour tout

0<α<p12p.0\lt \alpha\lt \frac{p-1}{2p}.

En faisant tendre p+p\to +\infty, on obtient toute valeur

0<α<12.0\lt \alpha\lt \frac12.

Par conséquent,

B admet une modification aˋ trajectoires localement α-ho¨ldeˊriennes pour tout 0<α<12.\boxed{B\text{ admet une modification à trajectoires localement }\alpha\text{-höldériennes pour tout }0\lt \alpha\lt \frac12.}