$\textstyle \parbox{10cm}{\scriptsize Universit\'e de Paris 5 -- Ren\'e Descart... ...matiques et Informatique}\ 45, rue des Saints-P\\lq eres 75270 Paris cedex 06 }$ \includegraphics[width=22mm]{/home/tcd/icones/p5_petit.ps}




Statistique des Processus - Maîtrise MASS


Examen du 10 juin 2003


9h - 12h


Les documents sont autorisés.


Le sujet comporte quatre pages.



Ex 1.Répondre aux questions suivantes par Vrai ou Faux, en justifiant votre choix :

  1. La somme de deux processus MA indépendants est un processus MA.
  2. La somme d'un processus AR et d'un processus MA indépendants est un processus ARMA.
  3. La somme de deux processus AR indépendants est un processus AR.
  4. Soit $ X$ un processus faiblement stationnaire. Si sa fonction d'autocorrélation s'annule à partir du rang $ 1$ inclus, alors sa fonction d'autocorrélation partielle est également nulle.
  5. Un processus ARMA caractérisé par un modèle non causal et non inversible ne possède pas de représentation de Wold.
  6. La densité spectrale d'un processus ARMA ne s'annule pas.
  7. Un modèle linéaire avec erreurs ARMA ne peut modéliser une série temporelle avec une tendance périodique.
  8. On peut valider un modèle linéaire avec erreurs corrélées en faisant un test de bruit blanc.
  9. Un processus ARIMA($ p,d,q$) est faiblement stationnaire si et seulement si $ d\leq \max(p,q)$.
  10. Le processus d'innovation d'un processus SARIMA est périodique.


Ex 2.Soit $ X$, $ Y$ et $ Z$ trois processus ARMA vérifiant pour tout $ t\in{\mathbb{Z}}$

$\displaystyle X_t$ $\displaystyle =$ $\displaystyle {\varepsilon}_t-1,6X_{t-1}-0,8X_{t-2}$  
$\displaystyle Y_t$ $\displaystyle =$ $\displaystyle {\varepsilon}_t+0,4{\varepsilon}_{t-1}+0,7Y_{t-1}$  
$\displaystyle Z_t$ $\displaystyle =$ $\displaystyle {\varepsilon}_t+0,4{\varepsilon}_{t-1}-0,5{\varepsilon}_{t-2}$  

avec $ {\varepsilon}$ bruit blanc de variance $ {\sigma}^2$.
  1. Préciser la nature de ces processus, déterminer leurs bruits blancs d'innovation, et calculer leurs fonctions d'autocovariance.
  2. On a tracé ci-dessous la densité spectrale de ces trois processus. En le justifiant, déterminer quel processus correspond à quelle densité.
    \includegraphics[scale=.31]{Dessins/SpeA.epsi-pg1} \includegraphics[scale=.31]{Dessins/SpeB.epsi-pg1} \includegraphics[scale=.31]{Dessins/SpeC.epsi-pg1}
  3. À partir d'une simulation du bruit blanc $ {\varepsilon}$, on a tracé les trois trajectoires suivantes, correspondant à une réalisation de chacun des processus $ X$, $ Y$ et $ Z$ entre les instants $ 1$ à $ 100$. En le justifiant, déterminer quel processus correspond à quelle trajectoire.
    \includegraphics[scale=.31]{Dessins/SimA.epsi-pg1} \includegraphics[scale=.31]{Dessins/SimB.epsi-pg1} \includegraphics[scale=.31]{Dessins/SimC.epsi-pg1}
  4. Pour chacune des réalisations précédentes, on a tracé l'autocorrélogramme correspondant (autocorrélation empirique et autocorrélation partielle empirique). En le justifiant, déterminer quel processus correspond à quel autocorrélogramme.
    \includegraphics[scale=.31]{Dessins/AcfA.epsi-pg1} \includegraphics[scale=.31]{Dessins/PacfA.epsi-pg1}
    \includegraphics[scale=.31]{Dessins/AcfB.epsi-pg1} \includegraphics[scale=.31]{Dessins/PacfB.epsi-pg1}
    \includegraphics[scale=.31]{Dessins/AcfC.epsi-pg1} \includegraphics[scale=.31]{Dessins/PacfC.epsi-pg1}


Ex 3.Soit $ {\varepsilon}=\left({\varepsilon}_t,t\in{\mathbb{Z}}\right)$ un bruit blanc de variance $ {\sigma}^2$. Soit $ X=(X_t,t\in{\mathbb{Z}})$ le processus AR($ 1$) tel que

$\displaystyle \forall t\in{\mathbb{Z}},\ \ X_t={\varepsilon}_t+{\varphi}X_{t-1}. $
avec $ \vert {\varphi}\vert<1$.
  1. On souhaite estimer $ {\varphi}$ et $ {\sigma}^2$ en fonction de $ (X_1,\ldots,X_N)$.
    1. Déterminer l'estimateur préliminaire de $ {\varphi}$.
    2. On pose pour tout $ t\geq 1$
      $\displaystyle \tilde {\varepsilon}_t=X_t-P^\bot_{V^2(1,X_1,\ldots,X_{t-1})}\left(X_t\right) $
      Calculer explicitement $ \tilde{\varepsilon}_t$ en fonction de $ (X_1,\ldots, X_t)$ pour tout $ t\geq 1$.
    3. Déterminer l'estimateur $ \hat{\varphi}_{mco}$ des moindres carrés ordinaires de $ {\varphi}$ et l'estimateur correspondant de $ {\sigma}^2$.
    4. Déterminer l'estimateur $ \hat{\varphi}_{mcc}$ des moindres carrés conditionnels de $ {\varphi}$. Comparer $ \hat{\varphi}_{mco}$ et $ \hat{\varphi}_{mcc}$.
    5. On suppose ici le bruit blanc $ {\varepsilon}$ gaussien. Proposer une équation caractérisant alors l'estimateur $ \hat{\varphi}_{mv}$ du maximum de vraisemblance de $ {\varphi}$. Montrer que
      $\displaystyle \left\vert\hat{\varphi}_{mco}-\hat{\varphi}_{mv}\right\vert\leq\frac{1}{\sqrt{\sum_{t=2}^{N-1}X_t^2}}. $
      Que peut-on en conclure ?
  2. En supposant $ {\varphi}$ et $ {\sigma}$ connus, déterminer pour tout $ h\geq 1$ le meilleur prédicteur linéaire de $ X_{N+h}$, sachant $ (X_1,\ldots,X_N)$, ainsi que la variance de l'erreur de prédiction.
  3. On pose pour tout $ t\in{\mathbb{Z}}$ $ Y_t=X_{-t}$.
    1. Montrer que $ Y=(Y_t,t\in{\mathbb{Z}})$ est un processus faiblement stationnaire, de même fonction d'autocovariance que $ X$.
    2. Préciser la nature du processus $ Y$, et calculer son bruit blanc d'innovation en fonction de $ {\varepsilon}$.
    3. On suppose dans toute la suite que
      $\displaystyle \forall t\in{\mathbb{Z}},\ \ \P\left({\varepsilon}_t=1\right)=\P\left({\varepsilon}_{t}=-1\right)=\frac{1}{2} $
      Déterminer la variance de $ {\varepsilon}_t$.
    4. On suppose encore que $ {\varphi}=1/3$, et que sont connus $ (Y_1,\ldots,Y_N)$. Sans calcul ou presque, déterminer la variance de l'erreur de prédiction commise au pas $ h$ (i.e. entre $ Y_{N+h}$ et le meilleur prédicteur linéaire de $ Y_{N+h}$), avec $ h\geq 1$.
    5. Montrer que pour tout $ h\geq 1$, $ Y_{N+h}$ est en fait une fonction de $ Y_N$, et qu'on peut donc le prédire sans erreur. Indication : représenter $ X_t$ sous forme MA($ \infty$).
    6. Expliquer l'apparente contradiction entre les deux questions précédentes.



Thierry Cabanal-Duvillard