📚 الرئيسية

مسابقة دكتوراه 2023Université 8 Mai 1945 - Guelma — الموضوع 03

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

Concours d'accès à la formation doctorale de 3ème cycle (LMD) 2022/2023, Spécialité : Probabilités et équations différentielles stochastiques, Épreuve écrite 2 : Modélisation stochastique, Université 8 Mai 1945, Faculté de Mathématiques, de l'Informatique et des Sciences de la Matière, Département de Mathématiques, 21 janvier 2023, durée 2h00mn.

التمرين 1

Exercice 1 — Couple aléatoire de densité $f_{X,Y}=4y(x-y)e^{-(x+y)}$ : marginale, conditionnelle, espérance

#joint-density#conditional-density#conditional-expectation#marginal-density

Soit (X,Y)(X, Y) un couple de variables aléatoires dont la loi admet la densité conjointe suivante :

fX,Y(x,y)=4y(xy)exp{(x+y)}1D(x,y),f_{X,Y}(x, y) = 4y(x - y)\exp\{-(x + y)\}\mathbf{1}_D(x, y), D={(x,y):0yx}.D = \{(x, y) : 0 \leq y \leq x\}.

  1. Déterminer la densité marginale de YY, puis la densité de la loi conditionnelle de XX sachant YY.
  2. Calculer E[XY]E[X|Y].
  3. Calculer P(X<1Y=y)P(X \lt 1|Y = y), selon les cas : y1y \leq 1 et y>1y \gt 1.
الحل

1. Marginale de Y et conditionnelle de X|Y

fY(y)=y4y(xy)e(x+y)dxf_Y(y) = \int_y^\infty 4y(x-y)e^{-(x+y)}dx. Changement u=xyu=x-y : 04yue(u+2y)du=4ye2yΓ(2)=4ye2y\int_0^\infty 4y \cdot u \cdot e^{-(u+2y)}du = 4ye^{-2y}\Gamma(2) = 4ye^{-2y}. Donc YΓ(2,1/2)Y \sim \Gamma(2, 1/2).

Conditionnelle : fXY(xy)=4y(xy)e(x+y)4ye2y=(xy)e(xy)f_{X|Y}(x|y) = \frac{4y(x-y)e^{-(x+y)}}{4ye^{-2y}} = (x-y)e^{-(x-y)} pour xyx \geq y. Donc XY=yy+Γ(2,1)X|Y=y \sim y + \Gamma(2,1).

2. E[X|Y]

E[XY=y]=y+2E[X|Y=y] = y + 2 (espérance de Γ(2,1)\Gamma(2,1) vaut 2).

E[XY]=Y+2\boxed{E[X|Y] = Y + 2}

3. P(X < 1 | Y = y)

Si y>1y \gt 1 : P(X<1Y=y)=0P(X \lt 1|Y=y) = 0. Si y1y \leq 1 : P(X<1Y=y)=y1(xy)e(xy)dxP(X \lt 1|Y=y) = \int_y^1 (x-y)e^{-(x-y)}dx. Posons u=xyu=x-y : 01yueudu=1(1+(1y))e(1y)=1(2y)e(1y)\int_0^{1-y} ue^{-u}du = 1-(1+(1-y))e^{-(1-y)} = 1-(2-y)e^{-(1-y)}.

P(X<1Y=y)={1(2y)e(1y)y10y>1\boxed{P(X \lt 1|Y=y) = \begin{cases}1-(2-y)e^{-(1-y)} & y \leq 1 \\\\ 0 & y \gt 1\end{cases}}

التمرين 2

Exercice 2 — EDS : processus $Y_t = tB_t$, formule d'Itô, espérance, covariance

#stochastic-differential-equations#ito-formula#brownian-motion#covariance

Soient (Ω,A,P,(Ft)t0)(\Omega, \mathcal{A}, \mathbb{P}, (\mathcal{F}_t)_{t\geq 0}) un espace probabilisé filtré, et (Bt)t0(B_t)_{t\geq 0} un (Ft)(\mathcal{F}_t)-mouvement brownien. On définit le processus Yt=tBtY_t = tB_t.

  1. Calculer dYtdY_t.
  2. Calculer l'espérance de YtY_t.
  3. Calculer E(YtYs)E(Y_t Y_s).
الحل

1.

Par la formule d'Itô appliquée à f(t,x)=txf(t,x) = tx : dYt=Btdt+tdBt\boxed{dY_t = B_t\,dt + t\,dB_t}

2.

E[Yt]=tE[Bt]=0E[Y_t] = t \cdot E[B_t] = 0.

3.

Pour sts \leq t : E[YtYs]=tsE[BtBs]=tss=ts2E[Y_tY_s] = ts\,E[B_tB_s] = ts \cdot s = ts^2.

E[YtYs]=tsmin(t,s)\boxed{E[Y_tY_s] = ts\min(t,s)}

التمرين 3

Exercice 3 — Coefficient de sécurité, sinistres géométriques et exponentiels, charge totale

#geometric-distribution#exponential-distribution#compound-sum#actuarial-science

Soit la formule du coefficient de sécurité : β=K+npμσn\beta = \frac{K + np\mu}{\sigma\sqrt{n}}.

On suppose qu'un risque peut être modélisé par un nombre de sinistres NN obéissant à la loi : P(N=k)=p(1p)kP(N=k) = p(1-p)^k, 0<p<10 \lt p \lt 1. Les montants de sinistres YY ont la densité : f(y)=λexp(λy)f(y) = \lambda\exp(-\lambda y), y>0y \gt 0.

  1. Quels sont l'espérance et la variance de NN ? Même question pour YY ?
  2. En déduire l'espérance et la variance de X=i=1NYiX = \sum_{i=1}^N Y_i.
  3. Déterminer les paramètres pp et λ\lambda pour que E(N)=0.1E(N) = 0.1 et E(Y)=9750E(Y) = 9750. Calculer E(X)E(X) et V(X)V(X).
  4. Quelles conditions sur nn pour que β4\beta \geq 4, avec K=500000K = 500000 DZD, σ=4468\sigma = 4468, ρμ=145\rho\mu = 145 ?
الحل

1.

NN suit une loi géométrique (à partir de 0) : E[N]=(1p)/pE[N] = (1-p)/p, Var(N)=(1p)/p2\text{Var}(N) = (1-p)/p^2. YE(λ)Y \sim \mathcal{E}(\lambda) : E[Y]=1/λE[Y] = 1/\lambda, Var(Y)=1/λ2\text{Var}(Y) = 1/\lambda^2.

2. Formule de Wald

E[X]=E[N]E[Y]=1ppλE[X] = E[N]\cdot E[Y] = \frac{1-p}{p\lambda}. Var(X)=E[N]Var(Y)+Var(N)(E[Y])2=(1p)(1+p)p2λ2\text{Var}(X) = E[N]\text{Var}(Y) + \text{Var}(N)(E[Y])^2 = \frac{(1-p)(1+p)}{p^2\lambda^2}.

3. Avec E(N)=0.1, E(Y)=9750

E[N]=(1p)/p=0.1p=10/11E[N] = (1-p)/p = 0.1 \Rightarrow p = 10/11. E[Y]=1/λ=9750λ=1/9750E[Y] = 1/\lambda = 9750 \Rightarrow \lambda = 1/9750. E[X]=0.1×9750=975E[X] = 0.1 \times 9750 = 975. Var(X)0.1×97502×(1+1/0.909)10,226,250\text{Var}(X) \approx 0.1 \times 9750^2 \times (1+1/0.909) \approx 10{,}226{,}250.

4. Condition sur n

β4\beta \geq 4 : K+nρμ4σnK + n\rho\mu \geq 4\sigma\sqrt{n}. 500000+145n17872n500000 + 145n \geq 17872\sqrt{n}. Posant u=nu=\sqrt{n} : 145u217872u+5000000145u^2 - 17872u + 500000 \geq 0. Δ=1787224×145×50000029,410,000\Delta = 17872^2 - 4\times145\times500000 \approx 29{,}410{,}000, Δ5423\sqrt{\Delta} \approx 5423. u42.9u \leq 42.9 ou u80.3u \geq 80.3, donc n1840n \leq 1840 ou n6450n \geq 6450.