📚 الرئيسية

مسابقة دكتوراه 2017Université Mohamed Khider - Biskra — الموضوع 06

مسابقة تخصص · Probabilités & Statistiques

Concours d'accès à la formation 3ème Cycle LMD, Épreuve de Processus Stochastique, Université Mohamed Khider - Biskra, Faculté des Sciences Exactes et des Sciences de la Nature et de la Vie, Département de Mathématiques, 21 octobre 2017 (version B : brownien, Poisson, modèle de Vasicek, BM 2D).

التمرين 1

Exercice 1 — Mouvement brownien : dérivabilité, martingales et intégrales

#brownian-motion#martingale#stochastic-integral#ito-formula

Soit (Bt)t0(B_t)_{t\geq 0} un mouvement brownien standard sur l'espace de probabilité filtré (Ω,F,(Ft),P)(\Omega,\mathcal{F},(\mathcal{F}_t),\mathbb{P}).

  1. Étudier la dérivabilité de Yt=BtY_t=B_t pour tout tRt\in\mathbb{R}.
  2. Pour tout αR\alpha\in\mathbb{R}, montrer que le processus Kt=eα2t/2cosh(αBt)K_t=e^{-\alpha^2 t/2}\cosh(\alpha B_t) est une Ft\mathcal{F}_t-martingale. Indication : cosh(x)=ex+ex2\cosh(x)=\dfrac{e^x+e^{-x}}{2}.
  3. Calculer E ⁣[(0t1dBs) ⁣(0tBs2dBs)]\mathbb{E}\!\left[\left(\int_0^t 1\,dB_s\right)\!\left(\int_0^t B_s^2\,dB_s\right)\right]. Indication : xy=14[(x+y)2(xy)2]xy=\dfrac{1}{4}[(x+y)^2-(x-y)^2].
  4. Montrer que E(Bt4)=3t2\mathbb{E}(B_t^4)=3t^2.
  5. Par deux méthodes différentes, montrer que (Bt33tBt)(B_t^3-3tB_t) est une martingale.
الحل

1. Dérivabilité

Le BM n'est p.s. nulle part dérivable (théorème classique), car ses trajectoires sont continues mais de variation totale infinie et de variation quadratique tt.

2. K_t martingale

Kt=12(eαBtα2t/2+eαBtα2t/2)K_t = \frac{1}{2}(e^{\alpha B_t-\alpha^2 t/2}+e^{-\alpha B_t-\alpha^2 t/2}). Chaque terme est l'exponentielle stochastique de ±αBt\pm\alpha B_t, qui est une martingale bornée dans L2L^2. La somme est aussi une martingale.

3. Isométrie d'Itô

Par l'indication : E[(0tdBs)(0tBs2dBs)]=14{E[(Bt+0tBs2dBs)2]E[(Bt0tBs2dBs)2]}E[(\int_0^t dB_s)(\int_0^t B_s^2 dB_s)]=\frac{1}{4}\{E[(B_t+\int_0^t B_s^2 dB_s)^2]-E[(B_t-\int_0^t B_s^2 dB_s)^2]\}

Par l'isométrie : E[(0t1dBs)(0tBs2dBs)]=0tBs21dsE[(\int_0^t 1 dB_s)(\int_0^t B_s^2 dB_s)]=\int_0^t B_s^2\cdot 1\,ds... plus précisément =0tE[Bs2]ds=0tsds=t22=\int_0^t E[B_s^2]\,ds=\int_0^t s\,ds=\frac{t^2}{2}.

4. E(B_t^4) = 3t²

Par la formule d'Itô : d(Bt4)=4Bt3dBt+6Bt2dtd(B_t^4)=4B_t^3 dB_t+6B_t^2 dt. En prenant l'espérance : ddtE[Bt4]=6E[Bt2]=6t\frac{d}{dt}E[B_t^4]=6E[B_t^2]=6t, donc E[Bt4]=3t2E[B_t^4]=3t^2.

5. B_t³ - 3tB_t martingale

Méthode 1 (Itô) : d(Bt3)=3Bt2dBt+3Btdtd(B_t^3)=3B_t^2 dB_t+3B_t dt, d(3tBt)=3Btdt+3tdBtd(3tB_t)=3B_t dt+3t dB_t. Donc d(Bt33tBt)=(3Bt23t)dBtd(B_t^3-3tB_t)=(3B_t^2-3t)dB_t, intégrale stochastique de carré intégrable : martingale.

Méthode 2 (critère) : E[Bt33tBtFs]=(Bs+ξ)33t(Bs+ξ)E[B_t^3-3tB_t|\mathcal{F}_s]=(B_s+\xi)^3-3t(B_s+\xi) avec ξ=BtBs\xi=B_t-B_s indépendant de Fs\mathcal{F}_s. En développant et utilisant E[ξ]=0E[\xi]=0, E[ξ2]=tsE[\xi^2]=t-s, E[ξ3]=0E[\xi^3]=0 : on obtient Bs33sBsB_s^3-3sB_s.

التمرين 2

Exercice 2 — Processus de Poisson et martingale exponentielle

#poisson-process#martingale#exponential-martingale#stochastic-process

Soit NtN_t un processus de Poisson d'intensité λ\lambda (λ>0\lambda\gt 0) défini sur l'espace de probabilité filtré (Ω,F,(Ft),P)(\Omega,\mathcal{F},(\mathcal{F}_t),\mathbb{P}), et soit Mt=exp(aNt+bt)M_t=\exp(aN_t+bt)aa et bb sont des constantes.

Trouver une relation entre aa et bb qui implique que MtM_t est une Ft\mathcal{F}_t-martingale.

الحل

Condition martingale

Pour que MtM_t soit une martingale, on a besoin de E[MtFs]=MsE[M_t|\mathcal{F}_s]=M_s pour sts\leq t.

E[MtFs]=eaNs+bsE[ea(NtNs)+b(ts)]E[M_t|\mathcal{F}_s] = e^{aN_s+bs}\cdot E[e^{a(N_t-N_s)+b(t-s)}]

Comme NtNsPoisson(λ(ts))N_t-N_s\sim\text{Poisson}(\lambda(t-s)) indépendant de Fs\mathcal{F}_s :

E[ea(NtNs)]=eλ(ts)(ea1)E[e^{a(N_t-N_s)}] = e^{\lambda(t-s)(e^a-1)}

Donc E[MtFs]=Mse(ts)[b+λ(ea1)]E[M_t|\mathcal{F}_s]=M_s\cdot e^{(t-s)[b+\lambda(e^a-1)]}.

Pour que cela soit =Ms=M_s, il faut :

b=λ(ea1)\boxed{b = -\lambda(e^a-1)}

التمرين 4

Exercice 4 — Mouvement brownien 2D : semimartingale et probabilité d'atteinte

#brownian-motion-2d#semimartingale#doob-decomposition#harmonic-function#hitting-probability

Soit (Bt)t0(B_t)_{t\geq 0} un (Ft)(\mathfrak{F}_t)-mouvement brownien en dimension deux, issu d'un point z=(x,y)z=(x,y) distinct de 0. On écrit Bt=(Xt,Yt)B_t=(X_t,Y_t)(Xt)(X_t) et (Yt)(Y_t) sont deux mouvements browniens de dimension 1 indépendants.

  1. On pose Rt=Bt2=Xt2+Yt2R_t=|B_t|^2=X_t^2+Y_t^2. Vérifier que (Rt)t0(R_t)_{t\geq 0} est une semimartingale et écrire sa décomposition. Vérifier que (R;R)t=40tRsds(R;R)_t=4\int_0^t R_s\,ds.
  2. On fixe ε]0,z[\varepsilon\in]0,|z|[ et on pose Tε=inf{t0:Btε}\mathcal{T}_\varepsilon=\inf\{t\geq 0:|B_t|\leq\varepsilon\} puis Ut=RtTεU_t=R_{t\wedge\mathcal{T}_\varepsilon}. Montrer que log(Ut)\log(U_t) est une martingale locale.
  3. Soit A]z,+[A\in]|z|,+\infty[ et τA=inf{t0:BtA}\tau_A=\inf\{t\geq 0:|B_t|\geq A\}. Montrer que P(Tε<τA)=log(A)log(z)log(A)log(ε)P(\mathcal{T}_\varepsilon\lt\tau_A)=\frac{\log(A)-\log(|z|)}{\log(A)-\log(\varepsilon)}
الحل

1. Semimartingale

Par la formule d'Itô : dRt=2XtdXt+2YtdYt+2dtdR_t=2X_t\,dX_t+2Y_t\,dY_t+2\,dt. Donc :

Rt=z2+20tXsdXs+20tYsdYs+2tR_t = |z|^2 + 2\int_0^t X_s\,dX_s+2\int_0^t Y_s\,dY_s+2t

Partie martingale : Mt=20tXsdXs+20tYsdYsM_t=2\int_0^t X_s dX_s+2\int_0^t Y_s dY_s. Partie VF : At=2tA_t=2t.

(R;R)t=40t(Xs2+Ys2)ds=40tRsds(R;R)_t=4\int_0^t(X_s^2+Y_s^2)\,ds=4\int_0^t R_s\,ds. \blacksquare

2. log(U_t) martingale locale

Sur [0,Tε)[0,\mathcal{T}_\varepsilon), Rtε2>0R_t\geq\varepsilon^2\gt 0. Par Itô : d(logRt)=1RtdRt12Rt2d(R;R)t=dMtRt+2dtRt2RtdtRt2=dMtRtd(\log R_t)=\frac{1}{R_t}dR_t-\frac{1}{2R_t^2}d(R;R)_t=\frac{dM_t}{R_t}+\frac{2dt}{R_t}-\frac{2R_t dt}{R_t^2}=\frac{dM_t}{R_t}

Donc logRt\log R_t est une intégrale stochastique, donc martingale locale. En arrêtant en Tε\mathcal{T}_\varepsilon, logUt\log U_t est aussi une martingale locale.

3. Probabilité d'atteinte

Par le théorème d'arrêt optionnel appliqué à la martingale locale logUt\log U_t aux temps d'arrêt TετA\mathcal{T}_\varepsilon\wedge\tau_A :

log(z2)=E[logUTετA]\log(|z|^2) = E[\log U_{\mathcal{T}_\varepsilon\wedge\tau_A}] =2log(ε)P(Tε<τA)+2log(A)P(τA<Tε)= 2\log(\varepsilon)\cdot P(\mathcal{T}_\varepsilon\lt\tau_A) + 2\log(A)\cdot P(\tau_A\lt\mathcal{T}_\varepsilon)

Avec P(Tε<τA)+P(τA<Tε)=1P(\mathcal{T}_\varepsilon\lt\tau_A)+P(\tau_A\lt\mathcal{T}_\varepsilon)=1, on résout :

P(Tε<τA)=logAlogzlogAlogε\boxed{P(\mathcal{T}_\varepsilon\lt\tau_A)=\frac{\log A-\log|z|}{\log A-\log\varepsilon}}