Méthode des moindres carrés/Supplément

Carl Friedrich Gauss

Méthode des moindres carrés

Traduction par Joseph Bertrand.
Mallet-Bachelier, 1855 (p. 70-112).

SUPPLÉMENT
À LA
THÉORIE DE LA COMBINAISON DES OBSERVATIONS
QUI EXPOSE AUX MOINDRES ERREURS.

---

PRÉSENTÉ LE 16 SEPTEMBRE 1826 À LA SOCIÉTÉ ROYALE DE GOTTINGUE.

---

Dans le Mémoire précédent, on a supposé que les quantités à déterminer, à l’aide d’observations imparfaites, fussent dépendantes de certains éléments inconnus, en fonction desquels on sût les exprimer : le problème consistait alors à déduire des observations, le plus exactement possible, la valeur de ces éléments.

Dans la plupart des cas, c’est effectivement ainsi que la question se présente ; mais quelquefois il en est un peu autrement, et l’on pourrait même se demander, au premier abord, si le problème se ramène au précédent. Il n’est pas rare, en effet, que les quantités auxquelles se rapportent les observations ne soient pas explicitement exprimées en fonction de certains éléments, et qu’elles ne paraissent réductibles à une telle forme que par des opérations difficiles ou ambiguës. Il arrive souvent, d’autre part, que la nature du problème fournit certaines conditions auxquelles les valeurs observées doivent rigoureusement satisfaire.

Cependant, en y regardant de plus près, on aperçoit que ce cas ne diffère pas essentiellement du précédent, et qu’il peut s’y ramener. Si l’on nomme, en effet, $\varpi$ le nombre des grandeurs observées, et $\sigma$ celui des équations de condition, en choisissant, parmi les premières, $\varpi -\sigma$ quantités, rien ne nous empêche de les considérer comme nos seules inconnues, les autres étant regardées comme des fonctions de celles-là, que les équations de condition définissent. Par cet artifice, nous rentrerons dans le cas du Mémoire précédent.

Néanmoins, quoique cette marche conduise souvent au résultat d’une manière assez commode, ou ne peut nier qu’elle ne soit moins naturelle, et il est, par conséquent, désirable de traiter le problème sous l’autre forme, qui admet, du reste, une solution très-élégante. Il y a plus : cette solution nouvelle, conduisant à des calculs plus rapides que la précédente toutes les fois que $\sigma$ est moindre que ${\frac {1}{2}}\,\varpi$, ou, ce qui revient au même, toutes les fois que le nombre des éléments que nous avons désigné par $\rho$ dans le Mémoire précédent, est plus grand que ${\frac {\varpi }{2}}$, il faudra, dans ce cas, préférer la solution nouvelle, quand bien même il serait facile, par la nature du problème, de faire disparaître sans ambiguïté les équations de condition.

Désignons par $v$, $v'$, $v''$, etc., les quantités, en nombre $\varpi$, dont les valeurs sont fournies par l’observation. Supposons qu’une inconnue dépende de celles-là et soit exprimée par une fonction connue $u$, de $v$, $v'$, $v''$, etc. Soient $l$, $l'$, $l''$, etc., ce que deviennent les quotients différentiels

${\begin{aligned}{\frac {\mathrm {d} u}{\mathrm {d} v}}&,&{\frac {\mathrm {d} u}{\mathrm {d} v'}}&,&{\frac {\mathrm {d} u}{\mathrm {d} v''}},\ldots ,\end{aligned}}$

lorsqu’on y substitue à $v$, $v'$, $v''$, etc., leurs valeurs véritables. Si l’on substituait à $v$, $v'$, $v''$, etc., dans la fonction $u$, leurs valeurs véritables, on obtiendrait aussi la véritable valeur de $u$ ; mais si les observations sont affectées d’erreurs $e$, $e'$, $e''$, etc., il en résultera pour $u$ une erreur totale représentée par

$le+l'e'+l''e''+\ldots ,$

pourvu que l’on puisse, comme nous le supposerons toujours (lorsque $u$ n’est pas linéaire), négliger les carrés et les produits des erreurs $e$, $e'$, $e''$, etc.

Quoique la grandeur des erreurs $e$, $e'$, $e''$, etc., soit incertaine, l’incertitude attachée à la valeur trouvée pour $u$ peut généralement se mesurer par l’erreur moyenne à craindre dans la détermination adoptée. D’après les principes développés dans le premier Mémoire, cette erreur moyenne est

${\sqrt {{l}^{2}{m}^{2}+{l'}^{2}{m'}^{2}+{l''}^{2}{m''}^{2}+\ldots }},$

$m$, $m'$, $m''$, etc., étant les erreurs moyennes des diverses observations. Si toutes les observations sont affectées du même degré d’incertitude, cette expression devient

$m\,{\sqrt {{l}^{2}+{l'}^{2}+{l''}^{2}+\ldots }}.$

Il est clair d’ailleurs qu’au degré d’approximation auquel nous nous arrêtons, on peut remplacer $l$, $l'$, $l''$, etc., par les valeurs que prennent les coefficients différentiels

${\begin{aligned}{\frac {\mathrm {d} u}{\mathrm {d} v}}&,&{\frac {\mathrm {d} u}{\mathrm {d} v'}}&,&{\frac {\mathrm {d} u}{\mathrm {d} v''}},\ldots ,\end{aligned}}$

lorsqu’on y remplace $v$, $v'$, $v''$, etc., par leurs valeurs observées.

3.

Lorsque les quantités $v$, $v'$, $v''$, etc., sont indépendantes, l’inconnue ne peut se déterminer que d’une seule manière, et l’incertitude attachée au résultat ne peut être ni évitée ni diminuée. Les observations fournissent une valeur de l’inconnue qui n’a, dans ce cas, rien d’arbitraire. Il en est tout autrement lorsque les quantités $v$, $v'$, $v''$, etc., sont liées par des relations nécessaires, que nous supposerons exprimées par $\sigma$ équations,

${\begin{aligned}\mathrm {X} &=0,&\mathrm {Y} &=0,&\mathrm {Z} &=0,\ldots ,\end{aligned}}$

$\mathrm {X}$, $\mathrm {Y}$, $\mathrm {Z}$, etc., désignant des fonctions données des indéterminées $v$, $v'$, $v''$, etc. ; car, à la fonction $u$, on peut, dans ce cas, substituer toute autre fonction $\mathrm {U}$ telle, que la différence $\mathrm {U} -u$ s’évanouisse identiquement en vertu des équations

${\begin{aligned}\mathrm {X} &=0,&\mathrm {Y} &=0,&\mathrm {Z} &=0,\ldots .\end{aligned}}$

Si les observations étaient rigoureusement exactes, cette substitution ne changerait en rien le résultat ; mais, en raison des erreurs inévitables, à chaque forme adoptée pour $u$ correspondra un résultat différent, et l’erreur commise, au lieu d’être

$le+l'e'+l''e''+\ldots ,$

deviendra

$\mathrm {L} e+\mathrm {L} 'e'+\mathrm {L} ''e''+\ldots ,$

en désignant par $\mathrm {L}$, $\mathrm {L} '$, $\mathrm {L} ''$, etc., les quotients différentiels

${\begin{aligned}{\frac {\mathrm {d} \mathrm {U} }{\mathrm {d} v}}&,&{\frac {\mathrm {d} \mathrm {U} }{\mathrm {d} v'}}&,&{\frac {\mathrm {d} \mathrm {U} }{\mathrm {d} v''}},\ldots .\end{aligned}}$

Quoiqu’il soit impossible d’assigner la valeur des diverses erreurs, nous pouvons cependant comparer les erreurs moyennes à craindre dans les diverses combinaisons. La combinaison la plus avantageuse sera celle qui donnera, à l’erreur moyenne, la valeur minimum. Cette erreur étant

${\sqrt {\mathrm {L} ^{2}{m}^{2}+\mathrm {L''} ^{2}{m'}^{2}+\mathrm {L''} ^{2}{m''}^{2}+\,\ldots }},$

nous devons chercher à rendre la somme

$\mathrm {L} ^{2}{m}^{2}+\mathrm {L''} ^{2}{m'}^{2}+\mathrm {L''} ^{2}{m''}^{2}+\ldots ,$

aussi petite que possible.

4.

Les fonctions $\mathrm {U}$, en nombre infini, par lesquelles on peut remplacer $u$, ne différeront les unes des autres, dans nos recherches, que par les valeurs qu’elles fourniront pour $\mathrm {L}$, $\mathrm {L} '$, $\mathrm {L} ''$, etc. : nous devons donc, avant tout, chercher les relations qui existent entre les systèmes de valeurs que peuvent prendre ces coefficients. Désignons par

${\begin{array}{c}{\begin{array}{lll}a,&a',&a'',\ldots ,\\b,&b',&b'',\ldots ,\\c,&c',&c'',\ldots ,\end{array}}\\\cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

les valeurs que prennent les coefficients

${\begin{array}{c}{\begin{array}{lll}{\dfrac {\mathrm {dX} }{\mathrm {d} v}},&{\dfrac {\mathrm {dX} }{\mathrm {d} v'}},&{\dfrac {\mathrm {dX} }{\mathrm {d} v''}},\ldots ,\\[0.75ex]{\dfrac {\mathrm {dY} }{\mathrm {d} v}},&{\dfrac {\mathrm {dY} }{\mathrm {d} v'}},&{\dfrac {\mathrm {dY} }{\mathrm {d} v''}},\ldots ,\\[0.75ex]{\dfrac {\mathrm {dZ} }{\mathrm {d} v}},&{\dfrac {\mathrm {dZ} }{\mathrm {d} v'}},&{\dfrac {\mathrm {dZ} }{\mathrm {d} v''}},\ldots ,\end{array}}\\\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

si l’on y substitue pour $v$, $v'$, $v''$, etc., leurs valeurs véritables. Il est clair que si l’on donne à $v$, $v'$, $v''$, etc., des accroissements $\mathrm {d} v$, $\mathrm {d} v'$, $\mathrm {d} v''$, etc., qui ne changent pas $\mathrm {X}$, $\mathrm {Y}$, $\mathrm {Z}$, etc., et leur laissent, par conséquent, la valeur zéro, ces accroissements, qui satisferont aux équations

${\begin{array}{l}{\begin{alignedat}{4}0{}={}&a\,&\mathrm {d} v&{}+{}a'&\mathrm {d} v'&{}+{}a''&\mathrm {d} v''&{}+{}\ldots ,\\0{}={}&b\,&\mathrm {d} v&{}+{}b'&\mathrm {d} v'&{}+{}b''&\mathrm {d} v''&{}+{}\ldots ,\\0{}={}&c\,&\mathrm {d} v&{}+{}c'&\mathrm {d} v'&{}+{}c''&\mathrm {d} v''&{}+{}\ldots ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

ne changeront rien à la valeur de $\mathrm {U} -u$, et l’on aura, par conséquent,

$0=(l-\mathrm {L} )\,\mathrm {d} v+(l'-\mathrm {L} ')\,\mathrm {d} v'+(l''-\mathrm {L} '')\,\mathrm {d} v''+\ldots .$

On en conclut facilement que $\mathrm {L}$, $\mathrm {L} '$, $\mathrm {L} ''$, etc., doivent avoir la forme

${\begin{array}{l}{\begin{alignedat}{6}&\mathrm {L} &{}={}&l&{}+{}&a&x{}+{}&b&y{}+{}&c&z+\ldots ,\\&\mathrm {L} '&{}={}&l'&{}+{}&a'&x{}+{}&b'&y{}+{}&c'&z+\ldots ,\\&\mathrm {L} ''&{}={}&l''&{}+{}&a''&x{}+{}&b''&y{}+{}&c''&z+\ldots ,\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}$

$x$, $y$, $z$, etc., désignant des multiplicateurs déterminés. Réciproquement, il est clair que, pour toutes les valeurs de $x$, $y$, $z$, etc., on pourra former une fonction $\mathrm {U}$, pour laquelle les valeurs $\mathrm {L}$, $\mathrm {L'}$, $\mathrm {L''}$, etc., seront précisément celles que fournissent ces équations, et cette fonction pourra, d’après ce qui précède, être substituée à $u$. La forme la plus simple qu’on puisse lui donner est

$\mathrm {U} =u+x\,\mathrm {X} +y\,\mathrm {Y} +z\,\mathrm {Z} +\ldots +u',$

$u'$ désignant une fonction de $v$, $v'$, $v''$, etc. qui s’annule identiquement lorsque $\mathrm {X}$, $\mathrm {Y}$, $\mathrm {Z}$, etc., sont nuls, et dont la valeur, dans le cas actuel, sera maximum ou minimum (car ses dérivées devront toutes s’annuler). Mais l’introduction de cette fonction n’apporte aucune différence dans les résultats.

5.

Il est maintenant facile d’assigner à $x$, $y$, $z$, etc., des valeurs telles, que la somme

$\mathrm {L} ^{2}m^{2}+\mathrm {L'} ^{2}{m'}^{2}+\mathrm {L''} ^{2}{m''}^{2}+\ldots$

soit un minimum.

Il est clair que pour atteindre ce but, la connaissance des erreurs moyennes absolues $m$, $m'$, $m''$, etc., n’est pas nécessaire ; il suffit de connaître leurs rapports. Introduisons, en effet, au lieu de ces quantités, les poids des observations, $p$, $p'$, $p''$, etc., c’est-à-dire des nombres réciproquement proportionnels à $m^{2}$, ${m'}^{2}$, ${m''}^{2}$, etc. Les quantités $x$, $y$, $z$, etc., devront être déterminées de telle sorte que le polynôme

${\frac {\left(ax+by+cz+\ldots +l\right)^{2}}{p}}+{\frac {\left(a'x+b'y+c'z+\ldots +l'\right)^{2}}{p'}}+\ldots$

acquière une valeur minimum. Supposons que $x^{0}$, $y^{0}$, $z^{0}$, etc., soient les valeurs déterminées de $x$, $y$, $z$, etc., auxquelles répond ce minimum, et adoptons les notations suivantes :

${\begin{array}{l}{\begin{alignedat}{4}&{\frac {a^{2}}{p}}&{}+{}&{\frac {{a'}^{2}}{p'}}&{}+{}&{\frac {{a''}^{2}}{p''}}&{}+{}&\ldots =(aa),\\[0.5ex]&{\frac {ab}{p}}&{}+{}&{\frac {a'b'}{p'}}&{}+{}&{\frac {a''b''}{p''}}&{}+{}&\ldots =(ab),\\[0.5ex]&{\frac {ac}{p}}&{}+{}&{\frac {a'c'}{p'}}&{}+{}&{\frac {a''c''}{p''}}&{}+{}&\ldots =(ac),\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \\{\begin{alignedat}{4}&{\frac {b^{2}}{p}}\,&{}+{}&\;{\frac {{b'}^{2}}{p'}}\,&{}+{}&\;{\frac {{b''}^{2}}{p''}}\,&{}+{}&\ldots =(bb),\\[0.5ex]&{\frac {c^{2}}{p}}&{}+{}&\;{\frac {{c'}^{2}}{p'}}\,&{}+{}&\;{\frac {{c''}^{2}}{p''}}\,&{}+{}&\ldots =(cc),\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \\{\begin{alignedat}{4}&{\frac {a\,l}{p}}\,&{}+{}&{\frac {a'l'}{p'}}\,&{}+{}&{\frac {a''l''}{p''}}&{}+{}&\ldots =(al),\\[0.5ex]&{\frac {b\,l}{p}}&{}+{}&{\frac {b'l'}{p'}}&{}+{}&{\frac {b''l''}{p''}}&{}+{}&\ldots =(bl),\\[0.5ex]&{\frac {c\,l}{p}}&{}+{}&{\frac {c'l'}{p'}}&{}+{}&{\frac {c''l''}{p''}}&{}+{}&\ldots =(cl),\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

La condition de minimum exige évidemment que l’on ait

(1)

$\left\{{\begin{array}{l}{\begin{alignedat}{4}0{}={}&(aa)\,&x^{0}{}+{}&(ab)\,&y^{0}{}+{}&(ac)\,&z^{0}+\ldots {}+{}&(al),\\0{}={}&(ba)\,&x^{0}{}+{}&(bb)\,&y^{0}{}+{}&(bc)\,&z^{0}+\ldots {}+{}&(bl),\\0{}={}&(ca)\,&x^{0}{}+{}&(cb)\,&y^{0}{}+{}&(cc)\,&z^{0}+\ldots {}+{}&(cl),\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right.$

Après que ces équations auront fourni les valeurs de $x^{0}$, $y^{0}$, $z^{0}$, etc., on posera

(2)

$\left\{{\begin{array}{l}{\begin{alignedat}{5}&a&x^{0}{}+{}&b&y^{0}{}+{}&c&z^{0}+\ldots {}+{}&l&{}={}&\mathrm {L} ,\\&a'&x^{0}{}+{}&b'&y^{0}{}+{}&c'&z^{0}+\ldots {}+{}&l'&{}={}&\mathrm {L} ',\\&a''&x^{0}{}+{}&b''&y^{0}{}+{}&c''&z^{0}+\ldots {}+{}&l''&{}={}&\mathrm {L} '',\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}\right.$

et la fonction la plus propre à déterminer notre inconnue, à laquelle correspond l’erreur moyenne minimum, sera celle dont les coefficients différentiels, pour les valeurs considérées des variables, seront égaux à $\mathrm {L}$, $\mathrm {L'}$, $\mathrm {L''}$, etc. Le poids de la détermination ainsi obtenue, sera

(3)

$\mathrm {P} ={\frac {1}{{\dfrac {\mathrm {L} ^{2}}{p}}+{\dfrac {\mathrm {L'} ^{2}}{p'}}+{\dfrac {\mathrm {L''} ^{2}}{p''}}+\ldots }}$

c’est-à-dire que ${\frac {1}{\mathrm {P} }}$ sera précisément la valeur que prend le polynôme considéré plus haut, pour les valeurs des variables $x$, $y$, $z$, etc., qui satisfont aux équations (1).

6.

Dans l’article précédent, nous avons montré à déterminer la fonction $\mathrm {U}$ qui fournit la détermination la plus convenable de l’inconnue $u$. Examinons actuellement quelle est la valeur qui en résulte Désignons-la par $\mathrm {K}$ : on l’obtiendra en substituant, dans $\mathrm {U}$, les valeurs observées des quantités $v$, $v'$, $v''$, etc. Soit $k$ la valeur que prend $u$ lorsqu’on y fait les mêmes substitutions, soit enfin $\varkappa$ la véritable valeur de cette inconnue, telle qu’on l’obtiendrait par la substitution des valeurs véritables de $v$, $v'$, $v''$, etc., soit dans $u$, soit dans $\mathrm {U}$. On aura

${\begin{alignedat}{4}k&{}={}\varkappa {}+{}&l\,&e{}+{}&l'&e'{}+{}&l''&e''\ldots ,\\\mathrm {K} &{}={}\varkappa {}+{}&\mathrm {L} \,&e{}+{}&\mathrm {L} '&e'{}+{}&\mathrm {L} ''&e''\ldots ,\end{alignedat}}$

et, par conséquent,

$\mathrm {K} =k+(\mathrm {L} -l)\,e+(\mathrm {L} '-l')\,e'+(\mathrm {L} ''-l'')\,e''+\ldots .$

Substituant, dans cette équation, à la place de $\mathrm {L} -l$, $\mathrm {L'} -l'$, $\mathrm {L''} -l''$, etc., leurs valeurs fournies par (2), et posant

(4)

$\left\{{\begin{array}{l}{\begin{alignedat}{4}&a\,&e{}+{}&a'&e'{}+{}&a''&e''+\ldots &{}={}{\mathcal {A}},\\&b\,&e{}+{}&b'&e'{}+{}&b''&e''+\ldots &{}={}{\mathcal {B}},\\&c\,&e{}+{}&c'&e'{}+{}&c''&e''+\ldots &{}={}{\mathcal {C}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}\right.$

nous aurons

(5)

$\mathrm {K} =k+{\mathcal {A}}\,x^{0}+{\mathcal {B}}\,y^{0}+{\mathcal {C}}\,z^{0}+\ldots .$

Il n’est pas possible de calculer ${\mathcal {A}}$, ${\mathcal {B}}$, ${\mathcal {C}}$, etc., par le moyen des formules (4), car les erreurs $e$, $e'$, $e''$, etc., qui y figurent, ont des valeurs inconnues, mais on voit facilement que ces quantités ${\mathcal {A}}$, ${\mathcal {B}}$, ${\mathcal {C}}$, etc., ne sont autre chose que les valeurs de $\mathrm {X}$, $\mathrm {Y}$, $\mathrm {Z}$, etc., qui correspondent aux valeurs observées de $v$, $v'$, $v''$, etc., et alors le système des équations (1), (3), (5), forme la solution complète de notre problème. Il est clair, en effet, que l’on peut appliquer au calcul de $a$, $a'$, $a''$, $\ldots$, $b$, $b'$, $b''$, etc., la remarque faite à la fin de l’article 2, à l’occasion des quantités $l$, $l'$, $l''$, etc., c’est-à-dire remplacer les valeurs véritables de $v$, $v'$, $v''$, etc., par les valeurs observées.

7.

On peut substituer à la formule (3), qui représente le poids de la détermination la plus probable, plusieurs expressions qu’il est utile d’indiquer ; remarquons, d’abord, qu’en ajoutant les équations (2) après les avoir multipliées par ${\frac {a}{p}}$, ${\frac {a'}{p'}}$, ${\frac {a''}{p''}}$, etc., on aura

$(aa)\,x^{0}+(ab)\,y^{0}+(ac)\,z^{0}+\ldots ={\frac {a\,\mathrm {L} }{p}}+{\frac {a'\mathrm {L} '}{p'}}+{\frac {a''\mathrm {L} ''}{p''}}+\ldots .$

Le premier membre est nul ; en désignant donc, d’après la notation adoptée, le second membre par $(a\,\mathrm {L} )$, on aura

$(a\,\mathrm {L} )=0,$

et de même

${\begin{aligned}(b\,\mathrm {L} )&=0,&(c\,\mathrm {L} )&=0,\ldots .\end{aligned}}$

Multiplions ensuite les équations (2) par ${\frac {\mathrm {L} }{p}}$, ${\frac {\mathrm {L} '}{p'}}$, ${\frac {\mathrm {L} ''}{p''}}$, etc. ; nous aurons, en les ajoutant,

${\frac {l\,\mathrm {L} }{p}}+{\frac {l'\mathrm {L} '}{p'}}+{\frac {l''\mathrm {L} ''}{p''}}+\ldots ={\frac {\mathrm {L} ^{2}}{p}}+{\frac {\mathrm {L'} ^{2}}{p'}}+{\frac {\mathrm {L''} ^{2}}{p''}}+\ldots ,$

et, par suite, nous avons cette seconde expression du poids,

$\mathrm {P} ={\frac {1}{{\dfrac {l\,\mathrm {L} }{p}}+{\dfrac {l'\mathrm {L} '}{p'}}+{\dfrac {l''\mathrm {L} ''}{p''}}+\ldots }}\cdot$

Si, enfin, nous ajoutons les mêmes équations (2) après les avoir multipliées par ${\frac {l}{p}}$, ${\frac {l'}{p'}}$, ${\frac {l''}{p''}}$, etc., nous obtiendrons une troisième expression du poids,

$\mathrm {P} ={\frac {1}{(al)\,x^{0}+(bl)\,y^{0}+(cl)\,z^{0}+\ldots +(ll)}},$

où l’on a posé, conformément à la notation précédente,

${\frac {{l}^{2}}{p}}+{\frac {{l'}^{2}}{p'}}+{\frac {{l''}^{2}}{p''}}+\ldots =(ll).$

On passera facilement de là à la quatrième expression du poids,

${\begin{aligned}{\frac {1}{\mathrm {P} }}&=(ll)-(aa)\,{x^{0}}^{2}-(bb)\,{y^{0}}^{2}-(cc)\,{z^{0}}^{2}-\ldots \\&{}-{}2\,(ab)\,x^{0}y^{0}-2\,(ac)\,x^{0}z^{0}-2\,(bc)\,y^{0}z^{0}-\ldots .\end{aligned}}$

8.

La solution générale que nous venons d’exposer s’applique principalement au cas où l’on n’a qu’une seule inconnue à déterminer. Lorsqu’on cherche, au contraire, les valeurs les plus plausibles de plusieurs inconnues, dépendant des mêmes observations, ou lorsqu’on ignore quelles sont les inconnues qu’il faut, de préférence, déduire des observations, il convient de procéder d’une manière différente, dont nous allons actuellement nous occuper.

Considérons $x$, $y$, $z$, etc., comme des indéterminées, et posons

(6)

$\left\{{\begin{array}{l}{\begin{alignedat}{4}&(aa)\,&x{}+{}&(ab)\,&y{}+{}&(ac)\,&z+\ldots &{}={}\xi ,\\&(ab)\,&x{}+{}&(bb)\,&y{}+{}&(bc)\,&z+\ldots &{}={}\eta ,\\&(ac)\,&x{}+{}&(bc)\,&y{}+{}&(cc)\,&z+\ldots &{}={}\zeta ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}\right.$

Supposons que l’on en déduise, par l’élimination,

(7)

$\left\{{\begin{array}{l}{\begin{alignedat}{4}&(\alpha \alpha )\,&\xi {}+{}&(\alpha \beta )\,&\eta {}+{}&(\alpha \gamma )\,&\zeta +\ldots &{}={}x,\\&(\beta \alpha )\,&\xi {}+{}&(\beta \beta )\,&\eta {}+{}&(\beta \gamma )\,&\zeta +\ldots &{}={}y,\\&(\gamma \alpha )\,&\xi {}+{}&(\gamma \beta )\,&\eta {}+{}&(\gamma \gamma )\,&\zeta +\ldots &{}={}z,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right.$

et remarquons, avant tout, que les coefficients placés symétriquement sont nécessairement égaux, c’est-à-dire que

${\begin{array}{l}{\begin{aligned}(\beta \alpha )&=(\alpha \beta ),\\(\gamma \alpha )&=(\alpha \gamma ),\\(\gamma \beta )&=(\beta \gamma ),\end{aligned}}\\\;\cdots \cdots \cdots \cdots \end{array}}$

comme cela résulte de la théorie de l’élimination et comme d’ailleurs nous le démontrons plus loin ; nous aurons

(8)

$\left\{{\begin{array}{l}{\begin{alignedat}{4}x^{0}{}={}&{}-{}(\alpha \alpha )\,&(al){}-{}&(\alpha \beta )\,&(bl){}-{}&(\alpha \gamma )\,&(cl)-\ldots ,\\y^{0}{}={}&{}-{}(\beta \alpha )\,&(al){}-{}&(\beta \beta )\,&(bl){}-{}&(\beta \gamma )\,&(cl)-\ldots ,\\z^{0}{}={}&{}-{}(\gamma \alpha )\,&(al){}-{}&(\gamma \beta )\,&(bl){}-{}&(\gamma \gamma )\,&(cl)-\ldots ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right.$

et en posant

(9)

$\left\{{\begin{array}{l}{\begin{alignedat}{4}&(\alpha \alpha )\,&{\mathcal {A}}{}+{}&(\alpha \beta )\,&{\mathcal {B}}{}+{}&(\alpha \gamma )\,&{\mathcal {C}}+\ldots &{}={}\mathrm {A} ,\\&(\beta \alpha )\,&{\mathcal {A}}{}+{}&(\beta \beta )\,&{\mathcal {B}}{}+{}&(\beta \gamma )\,&{\mathcal {C}}+\ldots &{}={}\mathrm {B} ,\\&(\gamma \alpha )\,&{\mathcal {A}}{}+{}&(\gamma \beta )\,&{\mathcal {B}}{}+{}&(\gamma \gamma )\,&{\mathcal {C}}+\ldots &{}={}\mathrm {C} ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}\right.$

nous obtiendrons

$\mathrm {K} =k-\mathrm {A} \,(al)-\mathrm {B} \,(bl)-\mathrm {C} \,(cl)-\ldots ;$

et si nous posons de plus,

(10)

$\left\{{\begin{array}{l}{\begin{alignedat}{5}&a&\mathrm {A} {}+{}&b&\mathrm {B} {}+{}&c&\mathrm {C} +\ldots &{}={}p&&\varepsilon ,\\&a'&\mathrm {A} {}+{}&b'&\mathrm {B} {}+{}&c'&\mathrm {C} +\ldots &{}={}p'&&\varepsilon ',\\&a''&\mathrm {A} {}+{}&b''&\mathrm {B} {}+{}&c''&\mathrm {C} +\ldots &{}={}p''&&\varepsilon '',\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right.$

il viendra

(11)

$\mathrm {K} =k-l\,\varepsilon -l'\varepsilon '-l''\varepsilon ''-\ldots .$

9.

La comparaison des équations (7) et (9) nous apprend que les quantités auxiliaires $\mathrm {A}$, $\mathrm {B}$, $\mathrm {C}$, etc., sont les valeurs que prennent les indéterminées $x$, $y$, $z$, etc., lorsque l’on suppose

${\begin{aligned}\xi &={\mathcal {A}},&\eta &={\mathcal {B}},&\zeta &={\mathcal {C}},\ldots ;\end{aligned}}$

on en conclut évidemment,

(12)

$\left\{{\begin{array}{l}{\begin{alignedat}{4}&(aa)\,&\mathrm {A} {}+{}&(ab)\,&\mathrm {B} {}+{}&(ac)\,&\mathrm {C} +\ldots &{}={}{\mathcal {A}},\\&(ba)\,&\mathrm {A} {}+{}&(bb)\,&\mathrm {B} {}+{}&(bc)\,&\mathrm {C} +\ldots &{}={}{\mathcal {B}},\\&(ca)\,&\mathrm {A} {}+{}&(cb)\,&\mathrm {B} {}+{}&(cc)\,&\mathrm {C} +\ldots &{}={}{\mathcal {C}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}\right.$

En ajoutant les équations (10) après les avoir multipliées par ${\frac {a}{p}}$, ${\frac {a'}{p'}}$, ${\frac {a''}{p''}}$, etc., on obtient

(13)

$\left\{{\begin{array}{l}{\begin{alignedat}{4}{\mathcal {A}}{}&={}a\,\varepsilon {}&+\;a'\varepsilon '{}&+{}&a''&\varepsilon ''+\ldots ,\\&\,\mathrm {et} \;\mathrm {de} &\mathrm {m{\hat {e}}me} &&&\\{\mathcal {B}}{}&={}b\,\varepsilon {}&+\;b'\varepsilon '{}&+{}&b''&\varepsilon ''+\ldots ,\\{\mathcal {C}}{}&={}c\,\varepsilon {}&+\;c'\varepsilon '{}&+{}&c''&\varepsilon ''+\ldots ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right.$

${\mathcal {A}}$ étant, comme on l’a dit, la valeur que prend $\mathrm {X}$ lorsqu’on substitue à $v$, $v'$, $v''$, etc., leurs valeurs observées ; on aperçoit facilement que si l’on applique à chacune de ces quantités les corrections $-\varepsilon$, $-\varepsilon '$, $-\varepsilon ''$, etc., la valeur de $\mathrm {X}$ devient égale à zéro, et que de même $\mathrm {Y}$, $\mathrm {Z}$, etc., s’évanouissent par cette hypothèse. L’équation (11) prouve aussi que $\mathrm {K}$ est la valeur que prend $u$ par suite des mêmes substitutions.

En nommant compensation des observations, l’application des corrections $-\varepsilon$, $-\varepsilon '$, $-\varepsilon ''$, etc., aux grandeurs directement observées, nous sommes conduit à une conséquence très-importante :

Les observations compensées comme nous l’avons indiqué, satisfont exactement à toutes les équations de condition, et font prendre à toute fonction des quantités observées la valeur qui résulterait de la combinaison la plus convenable des observations non modifiées ; et puisque les équations de condition sont trop peu nombreuses pour qu’on puisse en déduire les valeurs exactes des erreurs, nous aurons au moins trouvé, par ce qui précède, les erreurs les plus plausibles. C’est sous ce nom que les quantités $\varepsilon$, $\varepsilon '$, $\varepsilon ''$, etc., seront désormais désignées.

10.

Le nombre des observations étant plus grand que celui des équations de condition, outre le système des corrections les plus plausibles, on peut en trouver un nombre infini qui rendent les équations de condition exactement satisfaites.

Il importe d’examiner les relations qui lient entre eux ces divers systèmes. Soit $-\mathrm {E}$, $-\mathrm {E} '$, $-\mathrm {E} ''$, etc., un pareil système de corrections, autre que le système le plus plausible ; nous aurons

${\begin{array}{l}{\begin{alignedat}{4}&a\,&\mathrm {E} &{}+{}a'&\mathrm {E} '&{}+{}a''&\mathrm {E} ''+\ldots &{}={}{\mathcal {A}},\\&b\,&\mathrm {E} &{}+{}b'&\mathrm {E} '&{}+{}b''&\mathrm {E} ''+\ldots &{}={}{\mathcal {B}},\\&c\,&\mathrm {E} &{}+{}c'&\mathrm {E} '&{}+{}c''&\mathrm {E} ''+\ldots &{}={}{\mathcal {C}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

Multipliant ces équations par $\mathrm {A}$, $\mathrm {B}$, $\mathrm {C}$, etc., et ajoutant, il vient, en ayant égard aux équations (10),

$p\,\varepsilon \,\mathrm {E} +p'\varepsilon '\mathrm {E} '+p''\varepsilon ''\mathrm {E} ''+\ldots =\mathrm {A} \,{\mathcal {A}}+\mathrm {B} \,{\mathcal {B}}+\mathrm {C} \,{\mathcal {C}}+\ldots .$

Mais les équations (13) combinées de la même manière donnent

(14)

$p\,{\varepsilon }^{2}+p'{\varepsilon '}^{2}+p''{\varepsilon ''}^{2}+\ldots =\mathrm {A} \,{\mathcal {A}}+\mathrm {B} \,{\mathcal {B}}+\mathrm {C} \,{\mathcal {C}}+\ldots ;$

de la combinaison de ces résultats on déduit facilement

${\begin{aligned}p\,\mathrm {E} ^{2}&+p'\mathrm {E'} ^{2}+p''\mathrm {E''} ^{2}+\ldots =p\,{\varepsilon }^{2}+p'{\varepsilon '}^{2}+p''{\varepsilon ''}^{2}+\ldots \\&+p\,\left(\mathrm {E} -\varepsilon \right)^{2}+p'\left(\mathrm {E} '-\varepsilon '\right)^{2}+p''\left(\mathrm {E} ''-\varepsilon ''\right)^{2}+\ldots ,\end{aligned}}$

et, par suite, la somme

$p\,\mathrm {E} ^{2}+p'\mathrm {E'} ^{2}+p''\mathrm {E''} ^{2}+\ldots$

est, nécessairement, plus grande que

$p\,{\varepsilon }^{2}+p'{\varepsilon '}^{2}+p''{\varepsilon ''}^{2}+\ldots ,$

ce que l’on peut énoncer de la manière suivante :

Théorème. — Les carrés des corrections qui peuvent concilier les observations avec les équations de condition étant respectivement multipliés par les poids des observations correspondantes, donnent une somme minimum quand on adopte les corrections les plus plausibles.

On reconnaît précisément le principe des moindres carrés dont, au reste, les équations (12) et (10) peuvent facilement se déduire. La somme minimum, que nous désignerons désormais par $\mathrm {S}$, est égale, d’après l’équation (14), à

$\mathrm {A} \,{\mathcal {A}}+\mathrm {B} \,{\mathcal {B}}+\mathrm {C} \,{\mathcal {C}}+\ldots .$

11.

La détermination des erreurs les plus plausibles étant indépendante de $l$, $l'$, $l''$, etc., fournit la préparation la plus commode, quel que soit l’usage que l’on veuille faire ultérieurement des observations. En outre, on voit sans peine que, pour atteindre ce but, il n’est pas nécessaire d’effectuer l’élimination indéfinie, c’est-à-dire de calculer $(\alpha \alpha )$, $(\alpha \beta )$, etc., il suffit de déduire des équations (12), par une élimination définie, les quantités auxiliaires $\mathrm {A}$, $\mathrm {B}$, $\mathrm {C}$, etc., que nous nommerons, dans ce qui va suivre, les corrélatifs des équations de condition

${\begin{aligned}\mathrm {X} &=0,&\mathrm {Y} &=0,&\mathrm {Z} &=0,\ldots ;\end{aligned}}$

ces quantités seront ensuite substituées dans l’équation (10).

Cette méthode ne laisse rien à désirer lorsqu’on demande seulement les valeurs les plus plausibles des quantités fournies par l’observation. Mais il en est autrement lorsqu’on désire, en outre, le poids de chacune des valeurs trouvées. Quelle que soit alors celle des quatre formules précédentes que l’on veuille employer, il est indispensable de connaître $\mathrm {L}$, $\mathrm {L'}$, $\mathrm {L''}$, etc., ou, ce qui revient au même, $x^{0}$, $y^{0}$, $z^{0}$, etc. ; par cette raison, il sera utile d’étudier, de plus près, l’élimination qui fournit ces quantités et d’obtenir une méthode plus commode pour la détermination des poids.

12.

Les relations qui lient les quantités dont nous nous occupons sont notablement simplifiées par la considération de la fonction indéfinie du second degré

${\begin{aligned}(aa)\,x^{2}&+2\,(ab)\,xy+2\,(ac)\,xz+\ldots \\&+(bb)\,y^{2}+2\,(bc)\,yz+\ldots +(cc)\,z^{2}+\ldots ,\end{aligned}}$

que nous désignerons par $\mathrm {T}$.

Cette fonction est évidemment égale à

(15)

${\frac {\left(ax+by+cz+\ldots \right)^{2}}{p}}+{\frac {\left(a'x+b'y+c'z+\ldots \right)^{2}}{p'}}+\ldots$

De plus, on a évidemment

(16)

$\mathrm {T} =\xi \,x+\eta \,y+\zeta \,z+\ldots ;$

et si, enfin, ou exprime $x$, $y$, $z$, etc., au moyen des équations (7), en fonction de $\xi$, $\eta$, $\zeta$, etc., on aura

${\begin{aligned}\mathrm {T} =(\alpha \alpha )\,\xi ^{2}&+2\,(\alpha \beta )\,\xi \eta +2\,(\alpha \gamma )\,\xi \zeta +\ldots +(\beta \beta )\,\eta ^{2}\\&+2\,(\beta \gamma )\,\eta \zeta +\ldots +(\gamma \gamma )\,\zeta ^{2}+\ldots .\end{aligned}}$

La théorie développée plus haut fournit deux systèmes de valeurs déterminées pour les quantités $x$, $y$, $z$, etc., $\xi$, $\eta$, $\zeta$, etc. Le premier est

${\begin{aligned}x&=x^{0},&y&=y^{0},&z&=z^{0},\ldots ,\\\xi &=-(al),&\eta &=-(bl),&\zeta &=-(cl),\ldots ;\end{aligned}}$

à ce système correspond la valeur

$\mathrm {T} =(ll)-{\frac {1}{\mathrm {P} }},$

ainsi qu’on le voit en comparant à l’équation (16), la troisième forme du poids $\mathrm {P}$, ou par la considération directe de la forme (4).

Le second système de valeurs est

${\begin{aligned}x&=\mathrm {A} ,&y&=\mathrm {B} ,&z&=\mathrm {C} ,\ldots ,\\\xi &={\mathcal {A}},&\eta &={\mathcal {B}},&\zeta &={\mathcal {C}},\,\ldots ;\end{aligned}}$

la valeur correspondante de $\mathrm {T}$ est

$\mathrm {T} =\mathrm {S} ,$

comme cela est évident par les formules (10) et (15), et encore par les formules (14) et (16).

13.

Nous devons, avant tout, faire subir à la fonction $\mathrm {T}$ une transformation semblable à celle qui a été indiquée (Theoria Motus, art. 182), et, avec plus de développements, dans les Recherches sur Pallas.

Posons, à cet effet,

${\begin{array}{l}{\begin{alignedat}{2}(bb,1)&{}={}(bb)&{}-{}&{\frac {(ab)^{2}}{(aa)}},\\[0.75ex](bc,1)&{}={}(bc)&{}-{}&{\frac {(ab)\,(ac)}{(aa)}},\\[0.75ex](bd,1)&{}={}(bd)&{}-{}&{\frac {(ab)\,(ad)}{(aa)}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \\{\begin{alignedat}{2}(cc,2)&{}={}(cc)&{}-{}&{\frac {(ac)^{2}}{(aa)}}-{\frac {(bc,1)^{2}}{(bb,1)}},\\[0.75ex](cd,2)&{}={}(cd)&{}-{}&{\frac {(ac)\,(ad)}{(aa)}}-{\frac {(bc,1)\,(bd,1)}{(bb,1)}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \\{\begin{alignedat}{2}(dd,3)&{}={}(dd)&{}-{}&{\frac {(ad)^{2}}{(aa)}}-{\frac {(bd,1)^{2}}{(bb,1)}}-{\frac {(cd,2)^{2}}{(cc,2)}},\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

et, ensuite^[1]

${\begin{array}{l}{\begin{alignedat}{4}&(bb,1)\,&y{}+{}&(bc,1)\,&z{}+{}&(bd,1)\,&w+\ldots &{}={}\eta ',\\&&&(cc,2)\,&z{}+{}&(cd,2)\,&w+\ldots &{}={}\zeta '',\\&&&&&(dd,3)\,&w+\ldots &{}={}\varphi ''',\end{alignedat}}\\\qquad \qquad \qquad \qquad \quad \;\cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}$

on aura

$\mathrm {T} ={\frac {\xi ^{2}}{(aa)}}+{\frac {{\eta '}^{2}}{(bb,1)}}+{\frac {{\zeta ''}^{2}}{(cc,2)}}+{\frac {{\varphi '''}^{2}}{(dd,3)}}+\ldots ,$

et $\eta '$, $\zeta ''$, $\varphi '''$, etc., se déduiront du $\xi$, $\eta$, $\zeta$, $\varphi$, etc., par les équations suivantes :

${\begin{array}{l}{\begin{alignedat}{6}&\eta '&{}={}&\eta &{}-{}&{\frac {(ab)}{(aa)}}&\,\xi &,\\[0.75ex]&\zeta ''&{}={}&\zeta &{}-{}&{\frac {(ac)}{(aa)}}&\,\xi &{}-{}{\frac {(bc,1)}{(bb,1)}}&\,\eta '&,\\[0.75ex]&\varphi '''&{}={}&\varphi &{}-{}&{\frac {(ad)}{(aa)}}&\,\xi &{}-{}{\frac {(bd,1)}{(bb,1)}}&\,\eta '&{}-{}{\frac {(cd,2)}{(cc,2)}}&\,\zeta '',\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}$

et l’on en tirera facilement toutes les formules utiles à notre but. Ainsi, pour déterminer les corrélatifs $\mathrm {A}$, $\mathrm {B}$, $\mathrm {C}$, etc., nous poserons

(18)

$\left\{{\begin{array}{l}{\begin{alignedat}{6}&{\mathcal {B}}'&{}={}&{\mathcal {B}}&{}-{}&{\frac {(ab)}{(aa)}}&\,{\mathcal {A}}&,\\[0.75ex]&{\mathcal {C}}''&{}={}&{\mathcal {C}}&{}-{}&{\frac {(ac)}{(aa)}}&\,{\mathcal {A}}&{}-{}{\frac {(bc,1)}{(bb,1)}}&\,{\mathcal {B}}'&,\\[0.75ex]&{\mathcal {D}}'''&{}={}&{\mathcal {D}}&{}-{}&{\frac {(ad)}{(aa)}}&\,{\mathcal {A}}&{}-{}{\frac {(bd,1)}{(bb,1)}}&\,{\mathcal {B}}'&{}-{}{\frac {(cd,2)}{(cc,2)}}&\,{\mathcal {C}}'',\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}\right.$

et enfin $\mathrm {A}$, $\mathrm {B}$, $\mathrm {C}$, $\mathrm {D}$, etc., s’obtiendront par les formules suivantes, en commençant par la dernière :

(19)

$\left\{{\begin{array}{l}{\begin{aligned}(aa)\,\mathrm {A} +(ab)\,\mathrm {B} +(ac)\,\mathrm {C} +(ad)\,\mathrm {D} +\ldots &={\mathcal {A}},\\(bb,1)\,\mathrm {B} +(bc,1)\,\mathrm {C} +(bd,1)\,\mathrm {D} +\ldots &={\mathcal {B}}',\\(cc,2)\,\mathrm {C} +(cd,2)\,\mathrm {D} +\ldots &=\mathrm {C} '',\\(dd,3)\,\mathrm {D} +\ldots &={\mathcal {D}}''',\end{aligned}}\\\quad \qquad \qquad \qquad \qquad \qquad \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}\right.$

Pour exprimer $\mathrm {S}$, nous aurons la formule nouvelle

(20)

$\mathrm {S} ={\frac {{\mathcal {A}}^{2}}{(aa)}}+{\frac {{{\mathcal {B}}'}^{2}}{(bb,1)}}+{\frac {{{\mathcal {C}}''}^{2}}{(cc,2)}}+{\frac {{{\mathcal {D}}'''}^{2}}{(dd,3)}}+\ldots ;$

enfin, le poids $\mathrm {P}$, qu’il faut attribuer à la détermination la plus plausible de la quantité $u$, sera donné par la formule

(21)

${\frac {1}{\mathrm {P} }}=(ll)-{\frac {{(al)}^{2}}{(aa)}}-{\frac {{(bl,1)}^{2}}{(bb,1)}}-{\frac {{(cl,2)}^{2}}{(cc,2)}}-{\frac {{(dl,3)}^{2}}{(dd,3)}}-\ldots ;$

dans cette formule, on fait

(22)

$\left\{{\begin{alignedat}{5}&(bl,1)&{}={}&(bl)&{}-{}&{\frac {(ab)\,(al)}{(aa)}},\\&(cl,2)&{}={}&(cl)&{}-{}&{\frac {(ac)\,(al)}{(aa)}}&{}-{}&{\frac {(bc,1)\,(bl,1)}{(bb,1)}},\\&(dl,3)&{}={}&(dl)&{}-{}&{\frac {(ad)\,(al)}{(aa)}}&{}-{}&{\frac {(bd,1)\,(bl,1)}{(bb,1)}}&{}-{}&{\frac {(cd,2)\,(cl,2)}{(cc,2)}}\cdot \end{alignedat}}\right.$

Les formules (17), …, (21), dont la simplicité ne laisse rien à désirer, fournissent la solution complète de notre problème.

14.

Après avoir résolu le problème que nous avions en vue, nous allons aborder quelques questions secondaires qui éclaireront davantage l’ensemble de cette théorie.

Nous chercherons, en premier lieu, s’il peut arriver que l’élimination qui fournit $x$, $y$, $z$, etc., en fonction de $\xi$, $\eta$, $\zeta$, etc., devienne, dans certains cas, impossible. Cela aurait évidemment lieu si les fonctions $\xi$, $\eta$, $\zeta$, etc., n’étaient pas indépendantes les unes des autres. Supposons, pour un instant, qu’il en soit ainsi, et que l’une d’elles puisse s’exprimer en fonction des autres, de telle sorte que l’on ait la relation identique

$\alpha \xi +\beta \eta +\gamma \zeta +\ldots =0,$

$\alpha$, $\beta$, $\gamma$, etc., désignant des nombres déterminés.

On aura alors

${\begin{array}{l}{\begin{alignedat}{5}\alpha &\,(aa)&{}+{}\beta &\,(ab)&{}+{}\gamma &\,(ac)+\ldots &{}={}&0,\\\alpha &\,(ab)&{}+{}\beta &\,(bb)&{}+{}\gamma &\,(bc)+\ldots &{}={}&0,\\\alpha &\,(ac)&{}+{}\beta &\,(bc)&{}+{}\gamma &\,(cc)+\ldots &{}={}&0,\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}$

et si nous posons

${\begin{array}{l}{\begin{alignedat}{5}\alpha &\,a&{}+{}\beta &\,b&{}+{}\gamma &\,c&{}+{}\ldots {}={}&p\,\Theta ,\\\alpha &\,a'&{}+{}\beta &\,b'&{}+{}\gamma &\,c'&{}+{}\ldots {}={}&p'\Theta ',\\\alpha &\,a''&{}+{}\beta &\,b''&{}+{}\gamma &\,c''&{}+{}\ldots {}={}&p''\Theta '',\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}$

on en déduira

${\begin{array}{l}{\begin{alignedat}{5}a&\,\Theta {}+{}&a'&\Theta '{}+{}&a''&\Theta ''+\ldots &{}={}&0,\\b&\,\Theta {}+{}&b'&\Theta '{}+{}&b''&\Theta ''+\ldots &{}={}&0,\\c&\,\Theta {}+{}&c'&\Theta '{}+{}&c''&\Theta ''+\ldots &{}={}&0,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

et, par suite,

$p\,{\Theta }^{2}+p'{\Theta '}^{2}+p''{\Theta ''}^{2}+\ldots =0\;;$

$p$, $p'$, $p''$, etc. étant, par leur nature, tous positifs, cette équation exige

${\begin{aligned}\Theta &=0,&\Theta '&=0,&\Theta ''&=0,\ldots .\end{aligned}}$

Si nous considérons les différentielles complètes $\mathrm {dX}$, $\mathrm {dY}$, $\mathrm {dZ}$, etc., répondant aux valeurs de $v$, $v'$, $v''$, etc., immédiatement fournies par les observations, ces différentielles

${\begin{array}{l}{\begin{alignedat}{3}a&\,\mathrm {d} v{}+{}&a'&\mathrm {d} v'{}+{}&a''&\mathrm {d} v''+\ldots ,\\b&\,\mathrm {d} v{}+{}&b'&\mathrm {d} v'{}+{}&b''&\mathrm {d} v''+\ldots ,\\c&\,\mathrm {d} v{}+{}&c'&\mathrm {d} v'{}+{}&c''&\mathrm {d} v''+\ldots ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}$

d’après les résultats précédents, seront liées, les unes aux autres, de telle sorte qu’en les multipliant respectivement par $\alpha$, $\beta$, $\gamma$, etc., la somme des produits sera identiquement nulle, en sorte que, parmi les équations

${\begin{aligned}\mathrm {X} &=0,&\mathrm {Y} &=0,&\mathrm {Z} &=0,\ldots ,\end{aligned}}$

il en est une au moins que l’on peut regarder comme inutile, car elle sera satisfaite dès que les autres le seront.

En examinant la question de plus près, on voit que cette conclusion n’est applicable qu’à des valeurs des variables infiniment peu différentes de celles que fournit l’observation. Il y a, en effet, deux cas à distinguer : le premier est celui où l’une des équations

${\begin{aligned}\mathrm {X} &=0,&\mathrm {Y} &=0,&\mathrm {Z} &=0,\ldots ,\end{aligned}}$

est renfermée dans les autres d’une manière générale et absolue, et peut, par conséquent, être supprimée ; le second est celui où, pour les valeurs particulières de $v$, $v'$, $v''$, etc., auxquelles se rapportent les observations, l’une des fonctions $\mathrm {X}$, $\mathrm {Y}$, $\mathrm {Z}$, etc., $\mathrm {X}$ par exemple, acquiert une valeur maximum ou minimum, ou, plus généralement, une valeur dont la différentielle s’annule lorsque les autres équations restent satisfaites.

Mais comme nous ne considérons pour nos variables que des variations dont les carrés soient négligeables, ce second cas (qui dans les applications ne se présentera que bien rarement), pourra être assimilé au premier, et l’une des équations de condition pourra être supprimée comme surabondante.

Si les équations restantes sont indépendantes dans le sens que nous venons d’indiquer, on peut, d’après ce qui précède, être certain que l’élimination est possible. Nous nous réservons, du reste, de revenir sur cette matière qui mérite d’être examinée comme subtilité théorique plutôt que comme question d’une utilité pratique.

15.

Dans le premier Mémoire, art. 37 et suivants, nous avons montré le moyen de fixer, à posteriori, d’une manière très-approchée, le poids d’une détermination. Si les valeurs approchées de $\varpi$ quantités sont fournies par des observations également précises, et qu’on les compare avec les valeurs qui résultent pour elles des hypothèses les plus plausibles qu’on puisse faire sur les $\rho$ éléments dont elles dépendent, on a vu qu’il fallait ajouter les carrés des différences obtenues, diviser la somme par $\varpi -\rho$, et que le quotient pouvait être regardé comme une valeur approchée du carré de l’erreur moyenne inhérente à ce genre d’observations.

Si les observations sont inégalement précises, la seule modification que l’on doive apporter aux préceptes précédents, consiste en ce que l’on doit multiplier les carrés des différences par les poids respectifs des observations correspondantes, et l’erreur moyenne obtenue de cette manière se rapporte aux observations dont le poids est pris pour unité.

Dans le cas actuel, la somme des carrés dont nous parlons se confond évidemment avec la somme $\mathrm {S}$, et la différence $\varpi -\rho$ avec le nombre $\sigma$ des équations de condition. Par suite, pour l’erreur moyenne des observations dont le poids est 1, nous aurons l’expression ${\sqrt {\frac {\mathrm {S} }{\sigma }}}$, et la détermination sera d’autant plus digne de confiance que $\sigma$ sera plus considérable.

Mais il est bon d’établir ce résultat, indépendamment des raisonnements du premier Mémoire ; pour y parvenir, nous introduirons quelques notations nouvelles. Supposons qu’aux valeurs

${\begin{aligned}\xi &=a,&\eta &=b,&\zeta &=c,\ldots ,\end{aligned}}$

répondent

${\begin{aligned}x&=\alpha ,&y&=\beta ,&z&=\gamma ,\ldots ,\end{aligned}}$

de sorte que l’on ait

${\begin{alignedat}{5}&\alpha &{}={}a&\,(\alpha \alpha )&{}+{}b&\,(\alpha \beta )&{}+{}c&\,(\alpha \gamma )&{}+{}\ldots ,\\&\beta &{}={}a&\,(\alpha \beta )&{}+{}b&\,(\beta \beta )&{}+{}c&\,(\beta \gamma )&{}+{}\ldots ,\\&\gamma &{}={}a&\,(\alpha \gamma )&{}+{}b&\,(\beta \gamma )&{}+{}c&\,(\gamma \gamma )&{}+{}\ldots ;\end{alignedat}}$

et, en outre, qu’aux valeurs

${\begin{aligned}\xi &=a',&\eta &=b',&\zeta &=c',\ldots ,\end{aligned}}$

répondent

${\begin{aligned}x&=\alpha ',&y&=\beta ',&z&=\gamma ',\ldots ;\end{aligned}}$

enfin, qu’aux valeurs

${\begin{aligned}\xi &=a'',&\eta &=b'',&\zeta &=c'',\ldots ,\end{aligned}}$

répondent

${\begin{aligned}x&=\alpha '',&y&=\beta '',&z&=\gamma '',\ldots ,\end{aligned}}$

et ainsi de suite.

La combinaison des équations (4) et (9) fournit

${\begin{alignedat}{4}\mathrm {A} &{}={}\alpha &\,e{}+{}&\alpha '&e'{}+{}&\alpha ''&e''+\ldots ,\\\mathrm {B} &{}={}\beta &\,e{}+{}&\beta '&e'{}+{}&\beta ''&e''+\ldots ,\\\mathrm {C} &{}={}\gamma &\,e{}+{}&\gamma '&e'{}+{}&\gamma ''&e''+\ldots ,\end{alignedat}}$

et, comme on a

${\mathcal {S}}={\mathcal {A}}\,\mathrm {A} +{\mathcal {B}}\,\mathrm {B} +{\mathcal {C}}\,\mathrm {C} +\ldots ,$

on aura

${\begin{alignedat}{7}\mathrm {S} &{}={}(a&\,e{}+{}&a'&e'{}+{}&a''&e''+\ldots )\;&(\alpha &\,e{}+{}&\alpha '&e'{}+{}&\alpha ''&e''+\ldots )\\&{}+{}(b&\,e{}+{}&b'&e'{}+{}&b''&e''+\ldots )\;&(\beta &\,e{}+{}&\beta '&e'{}+{}&\beta ''&e''+\ldots )\\&{}+{}(c&\,e{}+{}&c'&e'{}+{}&c''&e''+\ldots )\;&(\gamma &\,e{}+{}&\gamma '&e'{}+{}&\gamma ''&e''+\ldots )&{}+{}\ldots .\end{alignedat}}$

16.

La série des observations qui fournissent les quantités $v$, $v'$, $v''$, etc., affectées des erreurs fortuites $e$, $e'$, $e''$, etc., peut être considérée comme une épreuve qui ne fait pas connaître, il est vrai, la grandeur de chaque erreur, mais qui, par le moyen des règles exposées plus haut, permet de déterminer la quantité $\mathrm {S}$, fonction connue de toutes les erreurs. Dans une telle épreuve, les erreurs peuvent être les unes plus grandes, les autres plus petites ; mais plus sera grand le nombre des erreurs employées, plus il y aura une grande probabilité que $\mathrm {S}$ diffère peu de sa valeur moyenne : la difficulté revient donc à trouver la moyenne de $\mathrm {S}$.

Par les principes exposés dans le premier Mémoire et qu’il est inutile de reproduire ici, on trouve pour cette valeur moyenne

$(a\,\alpha +b\,\beta +c\,\gamma +\ldots )\,{m}^{2}+(a'\alpha '+b'\beta '+c'\gamma '+\ldots )\,{m'}^{2}+\ldots .$

En nommant $\mu$ l’erreur moyenne qui correspond aux observations dont le poids est 1, de telle sorte que l’on ait

$\mu ^{2}=pm^{2}=p'{m'}^{2}=p''{m''}^{2}=\ldots ,$

l’expression précédente peut s’écrire comme il suit :

${\begin{aligned}&\left({\frac {a\,\alpha }{p}}\!+\!{\frac {a'\alpha '}{p'}}\!+\!{\frac {a''\alpha ''}{p''}}+\ldots \right)\mu ^{2}+\left({\frac {b\,\beta }{p}}+{\frac {b'\beta '}{p'}}+{\frac {b''\beta ''}{p''}}+\ldots \right)\mu ^{2}\\+{}&\left({\frac {c\,\gamma }{p}}+{\frac {c'\gamma '}{p'}}+{\frac {c''\gamma ''}{p''}}+\ldots \right)\mu ^{2}.\end{aligned}}$

Mais on a trouvé

${\frac {a\,\alpha }{p}}+{\frac {a'\alpha '}{p'}}+{\frac {a''\alpha ''}{p''}}+\ldots =(aa)(\alpha \alpha )+(ab)(\alpha \beta )+(ac)(\alpha \gamma )+\ldots ;$

or, le second membre est l’unité, comme on le reconnaît facilement par la comparaison des équations (6) et (7).

On trouvera de même

${\begin{alignedat}{4}&{\frac {b\,\beta }{p}}{}&&+{}{\frac {b'\beta '}{p'}}{}&&+{}{\frac {b''\beta ''}{p''}}{}&&+{}\ldots =1,\\&{\frac {c\,\gamma }{p}}{}&&+{}{\frac {c'\gamma '}{p'}}{}&&+{}{\frac {c''\gamma ''}{p''}}{}&&+{}\ldots =1,\end{alignedat}}$

et ainsi de suite.

D’après cela, la valeur moyenne de $\mathrm {S}$ devient $\sigma \mu ^{2}$, et si l’on juge permis de regarder la valeur fortuite de $\mathrm {S}$ comme égale à la valeur moyenne, on en conclut

$\mu ={\sqrt {\frac {\mathrm {S} }{\sigma }}}\cdot$

17.

On peut apprécier la confiance que mérite cette détermination en calculant l’erreur moyenne à craindre, soit pour sa valeur propre, soit pour celle de son carré. La seconde sera la racine carrée de la valeur moyenne de l’expression

$\left({\frac {\mathrm {S} }{\sigma }}-\mu ^{2}\right)^{2},$

dont le développement s’obtiendra par des raisonnements semblables à ceux qui ont été exposés dans le premier Mémoire (art. 39 et suivants). Nous les supprimons pour abréger, en nous contentant d’indiquer le résultat.

L’erreur moyenne à craindre dans la détermination du carré $\mu ^{2}$ s’exprime par

${\sqrt {{\frac {2\,\mu ^{4}}{\sigma }}+{\frac {\nu ^{4}-3\,\mu ^{4}}{\sigma ^{2}}}\,\mathrm {N} }},$

$\nu ^{4}$ étant la valeur moyenne des quatrièmes puissances des erreurs dont le poids est l’unité, et $\mathrm {N}$ la somme

${\begin{aligned}&(a\,\alpha +b\,\beta +c\,\gamma +\ldots )^{2}+(a'\alpha '+b'\beta '+c'\gamma '+\ldots )^{2}\\{}+{}&(a''\alpha ''+b''\beta ''+c''\gamma ''+\ldots )^{2}+\ldots .\end{aligned}}$

Cette somme ne peut pas, en général, se simplifier ; mais, par un procédé analogue à celui dont on a fait usage au paragraphe 40 du premier Mémoire, on peut montrer que sa valeur est comprise entre $\varpi$ et ${\frac {\sigma ^{2}}{\varpi }}$. Dans l’hypothèse sur laquelle nous avions primitivement établi la méthode des moindres carrés, le terme qui contient cette somme disparaît à cause de

$\nu ^{4}=3\,\mu ^{4},$

et la précision que l’on doit attribuer à la détermination

$\mu ={\sqrt {\frac {\mathrm {S} }{\sigma }}}$

est par conséquent la même que si l’on avait opéré sur $\sigma$ observations entachées d’erreurs exactement connues, conformément aux préceptes des art. 15 et 16 du premier Mémoire.

18.

Pour la compensation des observations, il y a, comme nous l’avons dit, deux opérations à exécuter : premièrement, il faut déterminer les corrélatifs des équations de condition, c’est-à-dire $\mathrm {A}$, $\mathrm {B}$, $\mathrm {C}$, etc., qui satisfont aux équations (12) ; secondement, substituer ces quantités dans l’équation (10). La compensation ainsi obtenue peut être dite parfaite et complète, par opposition à la compensation imparfaite ou incomplète. Nous désignerons de cette dernière manière celles qui résultent des mêmes équations (10), dans lesquelles on substituera des valeurs de $\mathrm {A}$, $\mathrm {B}$, $\mathrm {C}$ qui ne satisferont pas aux équations (12), c’est-à-dire qui satisferont à quelques-unes seulement, ou à aucune. Nous ne nous occuperons pas ici d’un tel système de corrections, et nous ne leur accorderons même pas le nom de compensation.

Lorsque les équations (10) sont satisfaites, les systèmes (12) et (13) deviennent équivalents, et la différence dont nous parlons peut alors s’énoncer comme il suit : Les observations complètement compensées satisfont aux équations de condition

${\begin{aligned}\mathrm {X} &=0,&\mathrm {Y} &=0,&\mathrm {Z} &=0,\ldots ;\end{aligned}}$

les observations incomplètement compensées ne satisfont qu’à une portion de ces équations, et peut-être à aucune ; la compensation à la suite de laquelle toutes les équations sont satisfaites, est nécessairement complète.

19.

Il résulte de la définition même des compensations, que la réunion de deux systèmes de compensations peut en fournir un troisième, et l’on voit qu’il importe peu que les règles données pour obtenir une compensation parfaite soient appliquées aux observations primitives ou aux observations déjà imparfaitement compensées.

Soient $-\Theta$, $-\Theta '$, $-\Theta ''$, etc., un système de compensations incomplètes, résultant des formules

(I)

$\left\{{\begin{array}{l}{\begin{alignedat}{5}&\Theta \,p&{}={}\mathrm {A} ^{0}&\,a&{}+{}\mathrm {B} ^{0}&\,b&{}+{}\mathrm {C} ^{0}&\,c&{}+{}\ldots ,\\&\Theta 'p'&{}={}\mathrm {A} ^{0}&\,a'&{}+{}\mathrm {B} ^{0}&\,b'&{}+{}\mathrm {C} ^{0}&\,c'&{}+{}\ldots ,\\&\Theta ''p''&{}={}\mathrm {A} ^{0}&\,a''&{}+{}\mathrm {B} ^{0}&\,b''&{}+{}\mathrm {C} ^{0}&\,c''&{}+{}\ldots ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right.$

Les observations ainsi changées ne satisfaisant pas à toutes les équations de condition, soient ${\mathcal {A}}^{\ast }$, ${\mathcal {B}}^{\ast }$, ${\mathcal {C}}^{\ast }$, etc., les valeurs que prennent $\mathrm {X}$, $\mathrm {Y}$, $\mathrm {Z}$, etc., quand on y substitue les valeurs ainsi obtenues pour $v$, $v'$, $v''$, etc. On devra chercher les valeurs $\mathrm {A} ^{\ast }$, $\mathrm {B} ^{\ast }$, $\mathrm {C} ^{\ast }$, etc., satisfaisant aux équations

(II)

$\left\{{\begin{array}{l}{\begin{alignedat}{5}{\mathcal {A}}^{\ast }&{}={}&\mathrm {A} ^{\ast }&\,(aa)&{}+{}\mathrm {B} ^{\ast }&\,(ab)&{}+{}\mathrm {C} ^{\ast }&\,(ac)&{}+{}\ldots ,\\{\mathcal {B}}^{\ast }&{}={}&\mathrm {A} ^{\ast }&\,(ab)&{}+{}\mathrm {B} ^{\ast }&\,(bb)&{}+{}\mathrm {C} ^{\ast }&\,(bc)&{}+{}\ldots ,\\{\mathcal {C}}^{\ast }&{}={}&\mathrm {A} ^{\ast }&\,(ac)&{}+{}\mathrm {B} ^{\ast }&\,(bc)&{}+{}\mathrm {C} ^{\ast }&\,(cc)&{}+{}\ldots ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right.$

et cela fait, la compensation complète des observations ainsi modifiées se fera par les nouveaux changements $-\varkappa$, $-\varkappa '$, $-\varkappa ''$, etc. ; $\varkappa$, $\varkappa '$, $\varkappa ''$, etc., se déduisant des formules

(III)

$\left\{{\begin{array}{l}{\begin{alignedat}{5}&\varkappa \,p&{}={}\mathrm {A} ^{\ast }&\,a&{}+{}\mathrm {B} ^{\ast }&\,b&{}+{}\mathrm {C} ^{\ast }&\,c&{}+{}\ldots ,\\&\varkappa 'p'&{}={}\mathrm {A} ^{\ast }&\,a'&{}+{}\mathrm {B} ^{\ast }&\,b'&{}+{}\mathrm {C} ^{\ast }&\,c'&{}+{}\ldots ,\\&\varkappa ''p''&{}={}\mathrm {A} ^{\ast }&\,a''&{}+{}\mathrm {B} ^{\ast }&\,b''&{}+{}\mathrm {C} ^{\ast }&\,c''&{}+{}\ldots ,\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{array}}\right.$

Cherchons comment ces corrections s’accordent avec la compensation complète des observations primitives. Il est clair d’abord que l’on a

${\begin{array}{l}{\begin{alignedat}{5}{\mathcal {A}}^{\ast }&{}={}{\mathcal {A}}&{}-{}&a&\,\Theta &{}-{}a'&\Theta '&{}-{}a''&\Theta ''{}-{}\ldots ,\\{\mathcal {B}}^{\ast }&{}={}{\mathcal {B}}&{}-{}&b&\,\Theta &{}-{}b'&\Theta '&{}-{}b''&\Theta ''{}-{}\ldots ,\\{\mathcal {C}}^{\ast }&{}={}{\mathcal {C}}&{}-{}&c&\,\Theta &{}-{}c'&\Theta '&{}-{}c''&\Theta ''{}-{}\ldots ,\end{alignedat}}\\\;\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \cdot \end{array}}$

En substituant dans ces équations, pour $\Theta$, $\Theta '$, $\Theta ''$, etc., leurs valeurs fournies par le système (I), pour ${\mathcal {A}}^{\ast }$, ${\mathcal {B}}^{\ast }$, ${\mathcal {C}}^{\ast }$, etc., celles que donne le système (II), il vient

${\begin{array}{l}{\begin{alignedat}{3}{\mathcal {A}}{}={}&(\mathrm {A} ^{0}+\mathrm {A} ^{\ast })\,(aa)&{}+{}&(\mathrm {B} ^{0}+\mathrm {B} ^{\ast })\,(ab)&{}+{}\ldots ,\\{\mathcal {B}}{}={}&(\mathrm {A} ^{0}+\mathrm {A} ^{\ast })\,(ab)&{}+{}&(\mathrm {B} ^{0}+\mathrm {B} ^{\ast })\,(bb)&{}+{}\ldots ,\\{\mathcal {C}}{}={}&(\mathrm {A} ^{0}+\mathrm {A} ^{\ast })\,(ac)&{}+{}&(\mathrm {B} ^{0}+\mathrm {B} ^{\ast })\,(bc)&{}+{}\ldots ,\end{alignedat}}\\\,\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdot \end{array}}$

et il suit de là que les corrélatifs des équations de condition (12) sont

${\begin{aligned}\mathrm {A} &=\mathrm {A} ^{0}+\mathrm {A} ^{\ast },&\mathrm {B} &=\mathrm {B} ^{0}+\mathrm {B} ^{\ast },&\mathrm {C} &=\mathrm {C} ^{0}+\mathrm {C} ^{\ast },\ldots ,\end{aligned}}$

et alors les équations (10), (I) et (III) montrent que l’on a

${\begin{aligned}\varepsilon &=\Theta +\varkappa ,&\varepsilon '&=\Theta '+\varkappa ',&\varepsilon ''&=\Theta ''+\varkappa '',\ldots ;\end{aligned}}$

et, par suite, la compensation parfaite a la même valeur pour chaque inconnue, soit qu’on la calcule directement, soit qu’on l’obtienne médiatement en partant d’une compensation incomplète.

20.

Lorsque les équations de condition sont très-nombreuses, la détermination des quantités corrélatives $\mathrm {A}$, $\mathrm {B}$, $\mathrm {C}$, etc., peut exiger des calculs tellement longs, que le calculateur en soit rebuté ; il pourra être avantageux d’obtenir, dans ce cas, une compensation complète, à l’aide d’une série d’approximations reposant sur le théorème de l’article précédent. On partagera, pour cela, les équations de condition en deux ou plusieurs groupes, et l’on cherchera d’abord une compensation qui rende satisfaites les équations du premier groupe. On traitera ensuite les valeurs modifiées par ce premier calcul, et on les corrigera de nouveau en ayant égard seulement aux équations du second groupe. Ce second calcul donnera des résultats qui, en général, ne satisferont plus aux équations du premier groupe, et il faudra, si l’on n’a formé que deux groupes, revenir alors au premier et y satisfaire à l’aide de nouvelles corrections. Les observations seront ensuite soumises à une quatrième compensation, dans laquelle on n’aura égard qu’aux conditions du second groupe ; et en opérant ainsi alternativement sur l’un et l’autre groupe d’équations, on formera des corrections qui seront nécessairement de plus en plus petites. Si le choix des groupes a été fait habilement, on arrivera bien vite à des valeurs que les corrections ultérieures ne changeront plus.

Quand on forme plus de deux groupes, on doit procéder de la même manière, les divers groupes étant employés successivement jusqu’au dernier, après quoi on revient au premier pour les reprendre dans le même ordre. Il nous suffit d’avoir indiqué ce procédé, dont la réussite dépendra beaucoup de l’habileté du calculateur.

21.

Il nous reste à donner la démonstration du lemme admis dans l’art. 8. Adoptons, pour plus de clarté, des notations plus propres à mettre la démonstration en lumière.

Soient $x^{0}$, $x'$, $x''$, etc., des indéterminées ; supposons que les équations