Produit de la densité et fonction de repartition d'une loi normale


Colle de mathématiques

Sujet de colle de maths:

Énoncé du sujet

  1. Calculer $\dsp\int_{-\infty}^{+\infty}x^2e^{-x^2}\,dx$ sachant que $\dsp\int_{-\infty}^{+\infty}e^{-x^2}\,dx=\sqrt\pi$
  2. Soit $X$ une variable aléatoire qui suit la loi normale centrée réduite.
    Soit $\varphi$ la densité et $\Phi$ la fonction de répartition de $X$. Pour tout réel $x$, on pose $g(x)=c\,\varphi(x)\Phi(x)$.
    1. Quelle relation a-t-on entre $\Phi(x)$ et $\Phi(-x)$ ?
    2. Quelle est la valeur de $c$ pour que $g$ soit la densité d'une variable aléaltoire $Y$ ?
    3. Calculer l'espérance et la variance de $Y$.



Correction

Correction

  1. En intégrant par parties,
    \[\begin{array}{ll}\dsp\int_{-\infty}^{+\infty}x^2e^{-x^2}\,dx &\dsp=\int_{-\infty}^{+\infty}x \left( xe^{-x^2}\right) dx\\[1.4em] &\dsp=\left[ x\dfrac{e^{-x^2}}{-2}\rb_{-\infty}^{+\infty} -\int_{-\infty}^{+\infty}\dfrac{e^{-x^2}}{-2}\,dx\\[1.4em] &=\dfrac12\sqrt\pi\end{array} \]

    car $\dsp\lim_{x\to-\infty}xe^{-x^2}=\lim_{x\to+\infty}xe^{-x^2}=0$.
  2. Soit $X$ une variable aléatoire qui suit la loi normale centrée réduite.
    Soit $\varphi$ la densité et $\Phi$ la fonction de répartition de $X$. Pour tout réel $x$, on pose $g(x)=c\,\varphi(x)\Phi(x)$.
    1. On a $\Phi(x)=\dsp\int_{-\infty}^x\varphi(t)dt$ et alors, en posant $u=-t$ et donc $du=-dt$,
      \[\Phi(-x)=P\left( X<-x\rp=\int_{-\infty}^{-x}\varphi(t)dt =\int_{+\infty}^x\varphi(-u)\lp-du\rp\]

      puis, comme $\varphi$ est paire, $\varphi(-u)=\varphi(u)$,
      \[\begin{array}{ll}\Phi(-x)&\dsp=\int_x^{+\infty}\varphi(u)du\\[.8em] &\dsp=\int_{-\infty}^{+\infty}\varphi(u)du-\int_{-\infty}^x\varphi(u)du\\[1em] &\dsp=1-\Phi(x) \enar\]

    2. On a clairement $g(x)>0$ pour tout réel $x$, et on souhaite donc aussi que
      \[\int_{-\infty}^{+\infty}g(t)dt =c\int_{-\infty}^{+\infty}\varphi(t)\Phi(t)dt =1\]

      On pose $I=\dsp\int_{-\infty}^{+\infty}\varphi(t)\Phi(t)dt$.
      1ère méthode: avec un changement de variable.
      Comme précédemment, avec le changement de variable $u=-t$, on a

      \[I=\int_{+\infty}^{-\infty}\varphi(-u)\Phi(-u)\lp-du\rp\]

      et donc, en utilisant la parité de $\varphi$ et la relation de la question précédente
      \[\begin{array}{ll}I&=\dsp\int_{-\infty}^{+\infty}\varphi(u)\lp1-\Phi(u)\right) du\\[1.3em] &\dsp=\int_{-\infty}^{+\infty}\varphi(u)du -\int_{-\infty}^{+\infty}\varphi(u)\Phi(u) du\\[1.3em] &=1-I\enar\]

      et on trouve donc que   $2I=1\iff I=1/2$.
      On doit donc avoir   $cI=1\iff c=2$.

      2ème méthode: avec une intégration par parties.
      On a $\Phi'=\varphi$ et donc
      \[\lp\Phi^2\rp'=2\Phi'\Phi=2\varphi\Phi\]

      et alors
      \[\begin{array}{ll} I&\dsp=\int_{-\infty}^{+\infty}\varphi(t)\Phi(t)dt\\[1em] &=\Bigl[\,\dfrac12\Phi^2(t)\,\Bigr]_{-\infty}^{+\infty} =\dfrac12 \enar\]

      car $\dsp\lim_{t\to+\infty}\Phi(t)=1$ et $\dsp\lim_{t\to-\infty}\Phi(t)=0$.
      On trouve donc que nécessirement $c=2$.

    3. L'espérance de $Y$ est
      \[E(Y)=\int_{-\infty}^{+\infty}yg(y)dx =2\int_{-\infty}^{+\infty}y\varphi(y)\Phi(y)dy\]

      On se rappelle que $\varphi(y)=\dfrac1{\sqrt{2\pi}}e^{-y^2/2}$ et donc $\varphi'(y)=-y\varphi(y)$, et que les fonctions de densité et de répartition sont reliées par $\Phi'(y)=\varphi(y)$.
      On peut donc penser à intégrer par parties avec $u(y)=\Phi(y)$ et $v'(y)=y\varphi(y)$, donc $u'(y)=\varphi(y)$ et $v(y)=-\varphi(y)$, et on a donc
      \[E(Y)=2\Bigl[-\Phi(y)\varphi(y)\Bigr]_{-\infty}^{+\infty} +2\int_{-\infty}^{+\infty}\varphi^2(y)dy\]

      Comme $\dsp\lim_{y\to\pm\infty}\varphi(y)=0$ et $\dsp\lim_{y\to-\infty}\Phi(y)=0$ et $\dsp\lim_{y\to+\infty}\Phi(y)=1$, on a
      \[\Bigl[-\Phi(y)\varphi(y)\Bigr]_{-\infty}^{+\infty}=0\]

      et donc
      \[\begin{array}{ll}E(Y) &\dsp=2\int_{-\infty}^{+\infty}\lp\dfrac1{\sqrt{2\pi}}e^{-y^2/2}\rp^2dy\\[1.5em] &\dsp=\dfrac1\pi\int_{-\infty}^{+\infty}e^{-y^2}dy\\[1.4em] &=\dfrac1\pi\tm\sqrt{\pi}=\dfrac1{\sqrt\pi}\enar\]



      La variance de $Y$ est
      \[V(Y)=E\left( Y^2\rp-\left( E(Y)\rp^2\]

      avec
      \[\begin{array}{ll}E\left( Y^2\rp&\dsp=\int_{-\infty}^{+\infty}x^2g(x)dx\\[1.3em] &\dsp=2\int_{-\infty}^{+\infty}x^2\varphi(x)\Phi(x)dx\\[1.3em] &\dsp=2\int_{+\infty}^{-\infty}(-u)^2\varphi(-u)\Phi(-u)\lp-du\rp\\[1.3em] &\dsp=2\int_{-\infty}^{+\infty}u^2\varphi(u)\lp1-\Phi(u)\right) du\\[1.3em] &=2E\left( X^2\rp-E\left( Y^2\rp\enar\]

      or
      \[\begin{array}{ll}V\left( X^2\rp=1&=E\left( X^2\rp-\left( E(X)\rp^2=E\left( X^2\rp-0\\[1em] &\iff E\left( X^2\rp=1\enar\]

      et alors
      \[E\left( Y^2\rp=2-E\left( Y^2\rp\iff E\left( Y^2\rp=1\]

      et donc, finalement,
      \[V(Y)=1-\dfrac1\pi\]



Tag:Variables aléatoires continues

Autres sujets au hasard: Lancer de dés

Voir aussi:
LongPage: h2: 3 - h3: 0