jOas wEb Sci

Aller au contenu | Aller au menu | Aller à la recherche

Assimilation

Fil des billets

mardi, novembre 13 2007

Interpolation statistique - Le problème de l'estimation

Le système étudié est décrit par \[ \mathbf{x}^t \]. La première estimation faite est \[ \mathbf{x}^b \] qui peut, par exemple, provenir d'une analyse antérieure. C'est la meilleure estimation du système en l'absence d'autres informations. Des observations \[ \mathbf{y}^o \] permettent d'obtenir des renseignements partiels au travers de l'opérateur d'observations non-linéaire \[ H \]. Dans la suite de ce chapitre, l'opérateur d'observations sera considéré linéaire et noté \[ \mathbf{H} \]. De plus, les erreurs d'ébauche \[ \epsilon^b \] et d'observation \[ \epsilon^o \] sont non-biaisées (ou débiaisées) et leurs statistiques sont connues.

L'objectif est alors, à l'aide des observations, d'améliorer l'estimation du système \[ \mathbf{x}^a \] vis-à-vis de la première estimation \[ \mathbf{x}^b \]. De plus, l'erreur sur l'état du système analysé \[ \epsilon^a \] est aussi recherchée.

Plusieurs possibilités existent pour obtenir l'analyse. Cependant, l'objectif est de réaliser une analyse aussi bonne qui possible (voire la meilleure). Il faut donc minimiser l'erreur commise a posteriori \[ \epsilon^a \], soit par exemple en minimisant \[ \mathrm{Tr}(\mathbf{A})\].

Interpolation statistique - Notations et hypothèses

Soit l'espace du modèle de dimensions n et l'espace des observations de dimensions p. En reprenant les définitions et notations précédentes :

Lire la suite...

Présentation du problème - Modélisation des erreurs

Pour prendre en compte les incertitudes dans l'ébauche, les observations et l'analyse, il faut faire des hypothèses sur la modélisation des erreurs entre ces vecteurs et leurs équivalents "vrais". L'utilisation des fonctions de densité de probabilité, nommées pdf, est une approche adaptée pour construire des modèles d'erreur. Les fonctions de densité de probabilité sont largement et rigoureusement décrites à travers des théories mathématiques. Une description simplifiée est donnée ci-dessous.

Lire la suite...

Présentation du problème - Description des vecteurs et espaces

Vecteur d'état

Pour formaliser mathématiquement le problème d'analyse, il faut définir un espace de travail. L'état du modèle est défini par une série de nombres ordonnés en une matrice colonne appelé vecteur d'état. Ce vecteur d'état \[ \mathbf{x} \] est le même que celui utilisé dans le modèle de prévision. Le choix de la discrétisation détermine comment les composantes vectorielles sont liées à l'état vrai du système. La discrétisation est ainsi équivalant au choix d'une base en mathématique.

Plusieurs vecteurs d'état peuvent être définis. Il n'est pas possible de décrire la réalité parfaitement car elle ne peut pas être représentée dans un vecteur d'état. Néanmoins, un vecteur \[ \mathbf{x}^t_c \] pourra permettre de décrire cet état vrai en continu. Plus classiquement, \[ \mathbf{x}^t" \] représentera la meilleure représentation possible de la réalité \[ \mathbf{x}^t_c" \] dans un vecteur d'état nommé état vrai. La première estimation de l'état analysé est appelé l'ébauche \[ \mathbf{x}^b" \], tandis que l'état analysé lui-même est noté \[ \mathbf{x}^a" \].

Lire la suite...

Présentation du problème - Concepts de base

L'analyse

L'analyse est une description fiable de l'état vrai du système à un instant donné. Elle est déjà utile par elle même en tant que représentation globale et consistante du système étudié. Elle peut aussi servir comme état initial pour une prévision du système à l'aide du modèle ou comme pseudo-observation. Elle peut aussi servir de référence afin de vérifier la qualité des observations.

Pour obtenir l'état analysé, les seules informations objectives sont les mesures des observations effectuées sur l'état vrai. Le système peut parfois être surdéterminé. Dans ce cas, l'analyse se résume à un problème d'interpolation. Il est, en général, sous-déterminé car les observations sont clairsemées et pas toujours liées directement aux variables du modèle. Ce qui n'empêche pas d'avoir des régions où les observations sont très denses et où le système est ainsi sur-déterminé. Afin de bien poser le problème, il est nécessaire de disposer d'une ébauche de l'état du modèle (c'est-à-dire une estimation a priori de l'état du modèle). Des contraintes physiques peuvent aussi permettre de mieux déterminer le système. Cette ébauche peut aussi bien être une climatologie, un état quelconque ou état obtenu à partir de précédentes analyses.  Dans ce cas, si le système efficace, l'information est sensée s'accumuler dans l'état du système et se propager entre les variables du modèle.

Lire la suite...

- page 6 de 7 -