Équivalence avec le BLUE

En reprenant exactement les mêmes hypothèses que pour le BLUE, il est possible de résoudre le problème par une approche variationnelle. Pour cela, il faut définir une fonctionnelle :
$\displaystyle J({\mathbf x})$ $\textstyle =$ $\displaystyle \frac{1}{2}({\mathbf x}-{\mathbf x}^b){\mathbf B}^{-1}({\mathbf x...
...\mathbf H}{\mathbf x})^T{\mathbf R}^{-1}({\mathbf y}^o-{\mathbf H}{\mathbf x}),$ (4.18)

appelée fonction coût et qui a pour caractéristique d'être quadratique en ${\mathbf x}$. Comme, de plus, les matrices ${\mathbf B}$ et ${\mathbf R}$ sont définies positives, alors cette fonction coût est convexe et possède un seul minimum qui peut être estimé par son gradient :
$\displaystyle \nabla J({\mathbf x})$ $\textstyle =$ $\displaystyle {\mathbf B}^{-1}({\mathbf x}-{\mathbf x}^b)-{\mathbf H}^T{\mathbf R}^{-1}({\mathbf y}^o-{\mathbf H}{\mathbf x}).$ (4.19)

Le point selle ${\mathbf x}^*$ est donc logiquement obtenu pour $\nabla J({\mathbf x})=0$ :
$\displaystyle {\mathbf B}^{-1}({\mathbf x}^*-{\mathbf x}^b)-{\mathbf H}^T{\mathbf R}^{-1}({\mathbf y}^o-{\mathbf H}{\mathbf x}^*)$ $\textstyle =$ $\displaystyle 0,$  
$\displaystyle {\mathbf B}^{-1}({\mathbf x}^*-{\mathbf x}^b)+{\mathbf H}^T{\mathbf R}^{-1}{\mathbf H}{\mathbf x}^* - {\mathbf H}^T{\mathbf R}^{-1}{\mathbf y}^o$ $\textstyle =$ $\displaystyle 0,$  
$\displaystyle {\mathbf B}^{-1}({\mathbf x}^*-{\mathbf x}^b)+{\mathbf H}^T{\mathbf R}^{-1}{\mathbf H}({\mathbf x}^*-{\mathbf x}^b)$ $\textstyle =$ $\displaystyle {\mathbf H}^T{\mathbf R}^{-1}{\mathbf y}^o - {\mathbf H}^T{\mathbf R}^{-1}{\mathbf H}{\mathbf x}^b,$  
$\displaystyle ({\mathbf B}^{-1}+{\mathbf H}^T{\mathbf R}^{-1}{\mathbf H})({\mathbf x}^*-{\mathbf x}^b)$ $\textstyle =$ $\displaystyle {\mathbf H}^T{\mathbf R}^{-1}({\mathbf y}^o-{\mathbf H}{\mathbf x}^b),$  

et la forme obtenue (Eq. 4.20) en écrivant cette égalité par rapport à $x^*$ est identique à l'Eq. 4.11 :
$\displaystyle {\mathbf x}^*$ $\textstyle =$ $\displaystyle {\mathbf x}^b + ({\mathbf B}^{-1}+{\mathbf H}^T{\mathbf R}^{-1}{\...
... H})^{-1}{\mathbf H}^T{\mathbf R}^{-1}({\mathbf y}^o-{\mathbf H}{\mathbf x}^b).$ (4.20)

La formulation variationnelle est donc strictement identique au BLUE à l'optimalité.

Nicolas Daget 2007-11-16