Définition du gain

Un ansatz pour le vecteur d'analyse ${\mathbf x}^a$ est de le décomposer de la manière suivante :

$${\mathbf x}^a = \mathbf{L}{\mathbf x}^b + {\mathbf K}{\mathbf y}^o,$$ (4.1)

où $\mathbf{L}$ est une matrice de dimension $n \times n$ et ${\mathbf K}$ une matrice de dimension $n \times p$. L'état analysé est donc une combinaison linéaire des différentes informations disponibles.

Son erreur associée peut être obtenu facilement en soustrayant ${\mathbf x}^t$ à l'Eq. 4.1.

$${\mathbf x}^a - {\mathbf x}^t = \mathbf{L}({\mathbf x}^b-{\mathbf x}^t+{\mathbf x}^t) + {\mathbf K}({\mathbf H}{\mathbf x}^t+\boldsymbol {\epsilon}^o) - {\mathbf x}^t$$

$$\boldsymbol {\epsilon}^a = \mathbf{L}\boldsymbol {\epsilon}^b + {\mathbf K}\boldsymbol {\epsilon}^o + (\mathbf{L}+{\mathbf K}{\mathbf H}-{\mathbf I}){\mathbf x}^t$$ (4.2)

En recherchant une erreur d'analyse sans biais et comme les erreurs d'ébauche et d'observation sont, par hypothèse, non biaisées, l'Eq. 4.2 donne :

$$E(\boldsymbol {\epsilon}^a) = \mathbf{L}E(\boldsymbol {\epsilon}^b) + {\mathbf K}E(\boldsymbol {\epsilon}^o) + (\mathbf{L}+{\mathbf K}{\mathbf H}-{\mathbf I})E({\mathbf x}^t),$$

$$E(\boldsymbol {\epsilon}^a) = (\mathbf{L}+{\mathbf K}{\mathbf H}-{\mathbf I})E({\mathbf x}^t).$$ (4.3)

Il faut donc que $\mathbf{L}+{\mathbf K}{\mathbf H}-{\mathbf I}=0$. L'ansatz prend donc la forme simplifiée suivantes :

$${\mathbf x}^a = ({\mathbf I}- {\mathbf K}{\mathbf H}){\mathbf x}^b + {\mathbf K}{\mathbf y}^o$$

$${\mathbf x}^a = {\mathbf x}^b + {\mathbf K}({\mathbf y}^o-{\mathbf H}{\mathbf x}^b)$$ (4.4)

On retrouve donc l'innovation ${\mathbf y}^o-{\mathbf H}{\mathbf x}^b$ définit précédemment et un opérateur linéaire ${\mathbf K}$ allant de l'espace des observations à l'espace du modèle généralement appelé gain.

Comme ${\mathbf K}$ est linéaire, l'analyse est donc une régression linéaire. Néanmoins, on parle, en général, d'interpolation linéaire pour des raisons historiques. En effet, les premières méthodes d'analyse, comme l'analyse de Cressman présentée dans la section 3.1.3, étaient réellement des interpolations linéaires au sens mathématique du terme. Pour obtenir une bonne analyse, il suffit donc de déterminer un gain satisfaisant.