📚 الرئيسية

مسابقة دكتوراه 2023Université Hassiba Benbouali - Chlef — الموضوع 03

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

Concours d'accès à la formation doctorale 3ème cycle LMD, Mathématiques Appliquées, Université Hassiba Benbouali - Chlef, Faculté des Sciences Exactes et Informatique, Département de Mathématiques, Épreuve de Spécialité (Troisième Variante), 21 janvier 2023 (Coefficient 03, Durée 02 heures).

التمرين 1

Exercice 1 — Statistique inférentielle : estimation et information de Fisher

#inferential-statistics#method-of-moments#maximum-likelihood#fisher-information

On considère XX une variable aléatoire continue, de densité :

fX(x)=12(1+αx)1]1,1[(x)f_X(x) = \frac{1}{2}(1 + \alpha x)\mathbb{1}_{]-1,1[}(x)

α]1,1[\alpha \in ]-1, 1[ est un paramètre inconnu. Soit X1,,XnX_1, \ldots, X_n un échantillon pour cette variable aléatoire.

  1. Trouver un estimateur α^n\hat{\alpha}_n de α\alpha par la méthode des moments. Déterminer la loi asymptotique de α^n\hat{\alpha}_n.
  2. Peut-on donner une expression explicite de l'estimateur du maximum de vraisemblance de α\alpha.
  3. Calculer l'information de Fisher. Quelle est la loi limite de l'EMV α^n\hat{\alpha}_n.
الحل

1.

On calcule E(X)=11x12(1+αx)dx=α3E(X) = \int_{-1}^{1} x \cdot \frac{1}{2}(1+\alpha x)dx = \frac{\alpha}{3}. L'estimateur des moments est α^n=3Xˉn\hat{\alpha}_n = 3\bar{X}_n. Par le TCL, n(α^nα)LN(0,9Var(X))\sqrt{n}(\hat{\alpha}_n - \alpha) \xrightarrow{\mathcal{L}} \mathcal{N}(0, 9\text{Var}(X)).

2.

La log-vraisemblance est (α)=nln2+ln(1+αXi)\ell(\alpha) = -n\ln 2 + \sum \ln(1+\alpha X_i). L'équation Xi1+αXi=0\sum \frac{X_i}{1+\alpha X_i} = 0 n'a pas de solution explicite en general.

3.

I(α)=E[(X1+αX)2]I(\alpha) = E\left[\left(\frac{X}{1+\alpha X}\right)^2\right]. L'EMV suit asymptotiquement N(α,1nI(α))\mathcal{N}(\alpha, \frac{1}{nI(\alpha)}).

التمرين 2

Exercice 2 — Programmation linéaire : simplexe et analyse de sensibilité

#linear-programming#simplex-method#sensitivity-analysis

Considérons le problème linéaire suivant :

maxz=x1\max z = x_1

s.c.{x1+x22x1+x28x1+x24x10,x20\text{s.c.} \begin{cases} -x_1 + x_2 \leq 2 \\\\ x_1 + x_2 \leq 8 \\\\ -x_1 + x_2 \geq -4 \\\\ x_1 \geq 0, x_2 \geq 0 \end{cases}

  1. Mettre le programme linéaire sous forme standard.
  2. Résoudre ce problème linéaire à l'aide de la méthode du simplexe.
  3. Résoudre ce problème géométriquement et tracez également les étapes de la procédure simplexe graphiquement.
  4. Supposons que la fonction objectif soit changée en z=x1+αx2z = x_1 + \alpha x_2, αR\alpha \in \mathbb{R}. a. Déterminer graphiquement les valeurs de α\alpha pour lesquelles la solution trouvée dans les questions (2.) et (3.) reste optimale. b. Déterminer graphiquement la valeur de α\alpha pour que ce problème ait une infinité de solutions optimales. Quelle est la valeur de zz^* dans ce cas ?
الحل

1.

On introduit les variables d'écart s1,s2,s30s_1, s_2, s_3 \geq 0: x1+x2+s1=2-x_1+x_2+s_1=2, x1+x2+s2=8x_1+x_2+s_2=8, x1x2+s3=4x_1-x_2+s_3=4.

2.

Le simplexe donne la solution optimale x1=6,x2=2x_1^*=6, x_2^*=2 avec z=6z^*=6.

3.

Graphiquement, le domaine réalisable est un polygone. Le maximum de z=x1z=x_1 est atteint au sommet le plus à droite.

4.

a.

La solution reste optimale pour 1α1-1 \leq \alpha \leq 1.

b.

Pour α=1\alpha=1, la fonction objectif est parallèle à la contrainte x1+x2=8x_1+x_2=8, donnant une infinité de solutions.

التمرين 3

Exercice 3 — Processus stochastiques : martingales et mouvement brownien

#stochastic-processes#brownian-motion#martingales#stopping-times

On travaille sur un espace de probabilité filtré (Ω,F,(Ft)t0,P)(\Omega, \mathcal{F}, (\mathcal{F}_t)_{t\geq 0}, \mathbb{P})(Ft)t0(\mathcal{F}_t)_{t\geq 0} est une filtration P\mathbb{P}-complète et continue à droite. Soit (Bt)t0(B_t)_{t\geq 0} un (Ft)t0(\mathcal{F}_t)_{t\geq 0}-mouvement brownien standard.

  1. Montrer que le processus (Mt)t0=(Bt2t)t0(M_t)_{t\geq 0} = (B_t^2 - t)_{t\geq 0} est une (Ft)t0(\mathcal{F}_t)_{t\geq 0}-martingale.
  2. Soit τa,b=inf{t>0:Bt]a,b[}\tau_{a,b} = \inf\{t \gt 0 : B_t \notin ]a, b[\} avec a<0<ba \lt 0 \lt b. a. Montrer que τa,b\tau_{a,b} est un (Ft)t0(\mathcal{F}_t)_{t\geq 0}-temps d'arrêt fini p.s. b. Montrer que E(B1τa,b)=0\mathbb{E}(B_{1\wedge\tau_{a,b}}) = 0 puis en déduire que E(Bτa,b)=0\mathbb{E}(B_{\tau_{a,b}}) = 0 et P(Bτa,b=a)=bba\mathbb{P}(B_{\tau_{a,b}} = a) = \frac{b}{b-a}. c. Montrer que E(M1τa,b)=0\mathbb{E}(M_{1\wedge\tau_{a,b}}) = 0 et en déduire que E(τa,b)=ab\mathbb{E}(\tau_{a,b}) = -ab. d. Soit Xt=12Ba2tX_t = \frac{1}{2}B_{a^2 t}, t0t \geq 0; et σa,b=inf{t>0:Xt]a,b[}\sigma_{a,b} = \inf\{t \gt 0 : X_t \notin ]a, b[\}, où a<0<ba \lt 0 \lt b. Comparer entre E(τa,b)\mathbb{E}(\tau_{a,b}) et E(σa,b)\mathbb{E}(\sigma_{a,b}). Justifier !
  3. Supposons maintenant que (Xt)t0(X_t)_{t\geq 0} est un processus càdlàg adapté à la filtration (Ft)t0(\mathcal{F}_t)_{t\geq 0}. a. Montrer que si pour tout (Ft)t0(\mathcal{F}_t)_{t\geq 0}-temps d'arrêt borné τ\tau; XτL1X_\tau \in L^1 et E(Xτ)=E(X0)\mathbb{E}(X_\tau) = \mathbb{E}(X_0) alors (Xt)t0(X_t)_{t\geq 0} est une (Ft)t0(\mathcal{F}_t)_{t\geq 0}-martingale. b. Pourquoi la filtration (Ft)t0(\mathcal{F}_t)_{t\geq 0} est supposée complète et continue à droite ainsi que le processus (Xt)t0(X_t)_{t\geq 0} est supposé càdlàg ?
  4. Soit (Mt)t0(M_t)_{t\geq 0} une (Ft)t0(\mathcal{F}_t)_{t\geq 0}-martingale et τ\tau un (Ft)t0(\mathcal{F}_t)_{t\geq 0}-temps d'arrêt. Soit (Mtτ)t0(M_t^\tau)_{t\geq 0} le processus arrêté défini par :

Mtτ:=Mtτ;t0( deˊsigne le minimum).M_t^\tau := M_{t\wedge\tau}; t \geq 0 \quad (\wedge \text{ désigne le minimum}).

a. Pourquoi on peut supposer que le processus (Mtτ)t0(M_t^\tau)_{t\geq 0} est càdlàg ? b. En déduire de ce qui précède que (Mtτ)t0(M_t^\tau)_{t\geq 0} est une (Ft)t0(\mathcal{F}_t)_{t\geq 0}-martingale.

الحل

1.

Par Itô: d(Bt2)=2BtdBt+dtd(B_t^2) = 2B_t dB_t + dt, donc Mt=Bt2t=20tBsdBsM_t = B_t^2 - t = 2\int_0^t B_s dB_s qui est une martingale.

2a.

{τa,bt}={supstBsb}{infstBsa}Ft\{\tau_{a,b} \leq t\} = \{\sup_{s\leq t} B_s \geq b\} \cup \{\inf_{s\leq t} B_s \leq a\} \in \mathcal{F}_t.

2b.

Par théorème d'arrêt: E(Bτa,b)=0\mathbb{E}(B_{\tau_{a,b}}) = 0. Comme Bτa,b{a,b}B_{\tau_{a,b}} \in \{a,b\}: P(Bτa,b=a)=bba\mathbb{P}(B_{\tau_{a,b}}=a) = \frac{b}{b-a}.

2c.

E(Bτa,b2τa,b)=0\mathbb{E}(B_{\tau_{a,b}}^2 - \tau_{a,b}) = 0 donne E(τa,b)=a2bba+b2aba=ab\mathbb{E}(\tau_{a,b}) = a^2\frac{b}{b-a} + b^2\frac{-a}{b-a} = -ab.

2d-4.

Les propriétés classiques des martingales arrêtées s'appliquent directement.