📚 الرئيسية

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

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

Concours d'accès à la formation de troisième cycle (Doctorat LMD) — Épreuve : Processus Stochastiques, Option Probabilités et EDS, Université Mohamed Khider - Biskra, Faculté des Sciences Exactes et des Sciences de la Nature et de la Vie, Département de Mathématiques — 21/10/2018.

التمرين 1

Exercice 1 — Intégrale stochastique d'Itô, loi de X_t et intégrale de Brownien

#stochastic-processes#ito-integral#brownian-motion#variance

Soient (Ω,F,(Ft)t0,P)(\Omega, \mathcal{F}, (\mathcal{F}_t)_{t \geq 0}, P) un espace de probabilité filtré et (Bt)t0(B_t)_{t \geq 0} un mouvement Brownien standard. Soit Xt=0tsdBsX_t = \displaystyle\int_0^t s\,dB_s.

  1. Calculer E(Xt)E(X_t) et Var(Xt)\mathrm{Var}(X_t). Quelle est la loi de XtX_t ?
  2. Calculer d(tBt)d(tB_t) à l'aide de la formule d'Itô.
  3. En déduire une relation entre XtX_t et Yt=0tBsdsY_t = \displaystyle\int_0^t B_s\,ds.
  4. Calculer la variance de YtY_t, a. directement à partir de sa définition, b. en calculant d'abord la covariance de BtB_t et XtX_t, à l'aide d'une partition de [0,t][0,t]. En déduire la loi de YtY_t.
الحل

1.

Par propriété de l'intégrale de Itô : E(Xt)=0E(X_t) = 0. L'isométrie d'Itô donne :

Var(Xt)=E(Xt2)=0ts2ds=t33.\mathrm{Var}(X_t) = E(X_t^2) = \int_0^t s^2\,ds = \frac{t^3}{3}.

Puisque XtX_t est une intégrale stochastique gaussienne (intégrande déterministe) :

XtN ⁣(0,t33).\boxed{X_t \sim \mathcal{N}\!\left(0,\,\frac{t^3}{3}\right).}

2.

On applique la formule d'Itô à f(t,Bt)=tBtf(t,B_t) = tB_t :

d(tBt)=Btdt+tdBt+0d(tBt)=Btdt+tdBt.d(tB_t) = B_t\,dt + t\,dB_t + 0 \quad \Longrightarrow \quad \boxed{d(tB_t) = B_t\,dt + t\,dB_t.}

3.

En intégrant la formule d'Itô de 00 à tt :

tBt=0tBsds+0tsdBs=Yt+Xt.tB_t = \int_0^t B_s\,ds + \int_0^t s\,dB_s = Y_t + X_t.

Donc Xt=tBtYt.\boxed{X_t = tB_t - Y_t.}

4a.

On calcule directement :

Var(Yt)=Cov ⁣(0tBsds,0tBudu)=0t0tCov(Bs,Bu)dsdu=0t0tmin(s,u)dsdu.\mathrm{Var}(Y_t) = \mathrm{Cov}\!\left(\int_0^t B_s\,ds, \int_0^t B_u\,du\right) = \int_0^t\int_0^t \mathrm{Cov}(B_s,B_u)\,ds\,du = \int_0^t\int_0^t \min(s,u)\,ds\,du.

Par symétrie :

=20t0usdsdu=20tu22du=0tu2du=t33.= 2\int_0^t\int_0^u s\,ds\,du = 2\int_0^t \frac{u^2}{2}\,du = \int_0^t u^2\,du = \frac{t^3}{3}.

4b.

De la relation Yt=tBtXtY_t = tB_t - X_t, on obtient Var(Yt)=t2Var(Bt)+Var(Xt)2tCov(Bt,Xt)\mathrm{Var}(Y_t) = t^2\mathrm{Var}(B_t) + \mathrm{Var}(X_t) - 2t\,\mathrm{Cov}(B_t,X_t). En calculant Cov(Bt,Xt)=0tsds=t2/2\mathrm{Cov}(B_t,X_t) = \int_0^t s\,ds = t^2/2 :

Var(Yt)=t2t+t332tt22=t3+t33t3=t33.\mathrm{Var}(Y_t) = t^2\cdot t + \frac{t^3}{3} - 2t\cdot\frac{t^2}{2} = t^3 + \frac{t^3}{3} - t^3 = \frac{t^3}{3}.

Puisque YtY_t est une intégrale d'une fonction déterministe d'un Brownien : YtN ⁣(0,t33).\boxed{Y_t \sim \mathcal{N}\!\left(0,\,\frac{t^3}{3}\right).}

التمرين 2

Exercice 2 — EDS linéaire : existence, unicité et moments

#stochastic-differential-equations#linear-sde#moments#ito-formula

Soient (Ω,F,(Ft)t0,P)(\Omega, \mathcal{F}, (\mathcal{F}_t)_{t \geq 0}, P) un espace de probabilité filtré et (Bt)t0(B_t)_{t \geq 0} un mouvement Brownien standard. On considère l'équation différentielle stochastique :

\begin{cases} dX_t = (a + \gamma X_t)\,dt + (b + \beta X_t)\,dB_t, \\ X_0 = x, \end{cases} \tag{1}

avec a,γ,b,βa, \gamma, b, \beta des constantes réelles et xRx \in \mathbb{R}.

  1. Montrer que l'équation (1) admet une solution unique.
  2. On note m(t)=E(Xt)m(t) = E(X_t) et M(t)=E(Xt2)M(t) = E(X_t^2). i. Montrer que m(t)m(t) est l'unique solution de l'équation différentielle ordinaire : {yγy=a,y(0)=x.\begin{cases} y' - \gamma y = a, \\ y(0) = x. \end{cases} (2) ii. Ecrire Xt2X_t^2XtX_t est la solution de l'équation (1). iii. En déduire que M(t)M(t) est l'unique solution de l'équation différentielle ordinaire : \begin{cases} y' - (2\gamma+\beta^2)y = 2(a+b\beta)m + b^2, \\ y(0) = x^2, \end{cases} \tag{3}mm est la solution de (2). iv. Calculer E(Xt)E(X_t) et E(Xt2)E(X_t^2).
الحل

1.

Les coefficients de dérive f(t,x)=a+γxf(t,x) = a+\gamma x et de diffusion σ(t,x)=b+βx\sigma(t,x) = b+\beta x sont lipschitziens en xx (uniformément en tt) et à croissance linéaire. Par le théorème d'existence et d'unicité des EDS lipschitziennes, l'équation (1) admet une unique solution forte.

2i.

En prenant l'espérance de (1) (l'intégrale stochastique est de moyenne nulle) :

m(t)=ddtE(Xt)=E(a+γXt)=a+γm(t),m(0)=x.m'(t) = \frac{d}{dt}E(X_t) = E(a + \gamma X_t) = a + \gamma m(t), \quad m(0) = x.

C'est exactement l'équation (2). Par unicité des EDO : m(t)=xeγt+aγ(eγt1)\boxed{m(t) = xe^{\gamma t} + \frac{a}{\gamma}(e^{\gamma t}-1)} (si γ0\gamma\neq 0).

2ii.

Par la formule d'Itô appliquée à f(Xt)=Xt2f(X_t) = X_t^2 :

dXt2=2XtdXt+(dXt)2=2Xt[(a+γXt)dt+(b+βXt)dBt]+(b+βXt)2dt.dX_t^2 = 2X_t\,dX_t + (dX_t)^2 = 2X_t[(a+\gamma X_t)dt + (b+\beta X_t)dB_t] + (b+\beta X_t)^2 dt.

dXt2=[2aXt+(2γ+β2)Xt2+2bβXt+b2]dt+2Xt(b+βXt)dBt.dX_t^2 = [2aX_t + (2\gamma+\beta^2)X_t^2 + 2b\beta X_t + b^2]\,dt + 2X_t(b+\beta X_t)\,dB_t.

2iii.

En prenant l'espérance (l'intégrale stochastique est une martingale, donc de moyenne nulle) :

M(t)=2aE(Xt)+(2γ+β2)M(t)+2bβE(Xt)+b2=(2γ+β2)M(t)+2(a+bβ)m(t)+b2.M'(t) = 2aE(X_t) + (2\gamma+\beta^2)M(t) + 2b\beta E(X_t) + b^2 = (2\gamma+\beta^2)M(t) + 2(a+b\beta)m(t) + b^2.

Cela donne l'équation (3). Avec M(0)=E(X02)=x2M(0) = E(X_0^2) = x^2, c'est une EDO du premier ordre à coefficients constants non homogène, dont la solution est unique.

2iv.

De (2) : E(Xt)=m(t)=xeγt+aγ(eγt1).\boxed{E(X_t) = m(t) = xe^{\gamma t} + \frac{a}{\gamma}(e^{\gamma t}-1).} (pour γ0\gamma \neq 0)

Et E(Xt2)=M(t)E(X_t^2) = M(t) est la solution de (3), que l'on obtient par variation de la constante.

التمرين 3

Exercice 3 — Martingales, surmartingales, sous-martingales et inégalité de Doob

#martingales#submartingale#doob-inequality#iid-variables
  1. Soient Y1,Y2,Y_1, Y_2,\ldots des variables aléatoires indépendantes identiquement distribuées définies sur un espace probabilisé filtré (Ω,F,{Fn}n,P)(\Omega, F, \{\mathcal{F}_n\}_n, P){Fn}n\{\mathcal{F}_n\}_n est la filtration canonique, et soit Xn=i=1nYiX_n = \prod_{i=1}^n Y_i. Sous quelle condition la suite XnX_n est-elle une surmartingale ? Une sous-martingale ? Une martingale ?

  2. On suppose que Y1,Y2,Y_1, Y_2,\ldots sont d'espérance nulle et variance finie σ2\sigma^2. Montrer que Zn=i=1nYiZ_n = \sum_{i=1}^n Y_i est une martingale. En utilisant l'inégalité de Doob, montrer que :

λ>0,P ⁣(max1inZiλ)1λ2Var(Zn)=nσ2λ2.\forall \lambda \gt 0, \quad P\!\left(\max_{1 \leq i \leq n} |Z_i| \geq \lambda\right) \leq \frac{1}{\lambda^2}\,\mathrm{Var}(Z_n) = \frac{n\sigma^2}{\lambda^2}.

  1. Améliorer la borne précédente en appliquant l'inégalité de Doob à (Zn+c)2(Z_n+c)^2 et en optimisant sur cc.
الحل

1.

Si E(Y1)=μE(Y_1) = \mu, alors E(XnFn1)=Xn1E(Yn)=μXn1E(X_n|\mathcal{F}_{n-1}) = X_{n-1}\cdot E(Y_n) = \mu X_{n-1}.

  • Martingale     μ=1\iff \mu = 1, c'est-à-dire E(Y1)=1E(Y_1) = 1.
  • Surmartingale     μ1\iff \mu \leq 1, c'est-à-dire E(Y1)1E(Y_1) \leq 1.
  • Sous-martingale     μ1\iff \mu \geq 1, c'est-à-dire E(Y1)1E(Y_1) \geq 1.

2.

On a Zn=Zn1+YnZ_n = Z_{n-1} + Y_n et YnY_n est indépendante de Fn1\mathcal{F}_{n-1}, donc E(ZnFn1)=Zn1+E(Yn)=Zn1E(Z_n|\mathcal{F}_{n-1}) = Z_{n-1} + E(Y_n) = Z_{n-1}. Donc ZnZ_n est une martingale.

Puisque ZnZ_n est une martingale, Zn2Z_n^2 est une sous-martingale. L'inégalité de Doob pour les sous-martingales donne :

P ⁣(max1inZiλ)=P ⁣(max1inZi2λ2)E(Zn2)λ2=nσ2λ2.P\!\left(\max_{1\leq i\leq n} |Z_i| \geq \lambda\right) = P\!\left(\max_{1\leq i\leq n} Z_i^2 \geq \lambda^2\right) \leq \frac{E(Z_n^2)}{\lambda^2} = \frac{n\sigma^2}{\lambda^2}.

3.

On applique l'inégalité de Doob à (Zn+c)2(Z_n+c)^2 (sous-martingale pour tout cc) :

P ⁣(max1inZi+cλ)E[(Zn+c)2]λ2=nσ2+c2λ2.P\!\left(\max_{1\leq i\leq n}|Z_i+c| \geq \lambda\right) \leq \frac{E[(Z_n+c)^2]}{\lambda^2} = \frac{n\sigma^2 + c^2}{\lambda^2}.

Nota : P(maxZiλ)P(maxZi+cλc)P(\max|Z_i|\geq \lambda) \leq P(\max|Z_i+c|\geq\lambda-|c|)... En optimisant directement sur cc dans la borne E[(Zn+c)2]λ2\frac{E[(Z_n+c)^2]}{\lambda^2}, le minimum sur cc est atteint en c=E(Zn)=0c = -E(Z_n) = 0. On obtient la même borne. Pour une amélioration effective, on peut optimiser P(maxi(Zi+c)2(λ+c)2)P(\max_{i}(Z_i+c)^2 \geq (\lambda+|c|)^2), ce qui donne des raffinements dépendant de la situation.