CALCUL DE L ERREUR MOYENNE
On trouverait de même e' (valeur moyenne) pour (y — B)|x', [z — C){a" Il en résulte qu'on a :
eî ^ it* 1^- [A2] d'où ,, [A^]
formule qui s'accorde avec le § 113; car, dans les observations directes, on n'avait qu'une seule inconnue à déterminer, d'où i = \.
Dans l'exemple précédent, on a donc
I /0, 08041 d'où £ = ± 0,283.
§ 133. — Il nous reste à résoudre une question trèsimportante. L'utilité de la méthode des moindres carrés ne réside pas seulement dans la précision du résultat qu'elle fournit : un de ses avantages les plus précieux est de permettre d'apprécier numériquement le degré d'exactitude que ce résultat comporte. Nous allons donc chercher le poids de nos inconnues.
Dans ce but, rappelons-nous (§ 125) qu'après avoir tiré de
l'équation différentielle
du = V do' + l" do" + /'" do'" ...
nous en avons conclu
P étant le poids de la fonction u; p celui des quantités
observées, o'. o" , o'" .... Ici, comme toutes les observations sont supposées également exactes, nous prendrons l'une d'entre elles pour unité de poids, et nous
aurons ^
Voyons ce qui arrive dans le cas actuel. — Puisque nos i inconnues ont été déterminées en fonction des e quantités observées 0, 0', 0" .... (§ 129), le problème serait résolu, d'après ce qui précède, si nous possédions i équations de la forme
a; = /-"(0, 0', 0" ...)elc.;
car nous en déduirions immédiatement
dx=aidO + a, dO' + asdO"...
(q) dî/=pidO + P»dO' + ?3dO"...
d;5 = YidO+Y,dO' + Y.dO"...etc.;
d'Où i. = t=c,3;^ = [p*];i = OT.
Nous n'avons pas directement ces i équations (q); mais nous allons les déduire des e équations (p). En effet, puisque 0 + A = V, on a ^V = c^O + rf^; différentiant les équations (p) et substituant aux dN leurs
valeurs, on aura
dA = — dO +a da7 + fc dî/ -f c d;5 4- ••• (q') dA' = — dO' -]ra'dx-\-b'ày^c'dz-\-...
^l^n do" + a"dx + b"dy + c"dz + ... etc. ;
équations qui, comparées à celles du groupe (p'), pourraient s'en déduire en changeant dans celles-ci les A en rfA; les X, 2/, z .... en dx, dy, dz..... et les n en dO. Si donc on traite le groupe (q') par la méthode des
moindres carrés, on obtiendra i équations normales qui
ne différeront du groupe (p^") que par les changements
que nous venons d'indiquer. Elles seront donc
\aa] d£ + [a&] dy-\-[ac\dz...^ [a dO] = 0 (q") • • • W àx + [W] dy ^- [6c] dz... — \bd(y\ = 0
[oc] dx + \pc\ dy + [ce] dz ... — [cdO']=0 ... etc.
Théoriquement, le problème est résolu; car l'élimination entre ces i équations permet d'exprimer les dx, dy, dz en fonction des dO, et de trouver par conséquent les valeurs des e X i coefficients (a, j3, y...) des équations (q).
Mais les calculs peuvent être singulièrement simplifiés par des considérations particulières, résultant de ce que nous n'avons besoin, en définitive, que de connaître les [a^j\ [p*^, [y*]....
Adoptons, dans ce but, la méthode d'élimination par les coefficients indéterminés, et multiplions chacune de nos équations (q") par les coefficients q^, Çt, q^... Il vient ainsi
[aa] q\dx -f- [at] qxdy -j- [ex] qidz ... — [a dO] qi =0 (q'") ... [aft] q^dx + [bb] q.dy + [bc] q^dz ... — [6 dO] ^j =■ 0
[ac] qzdx -f- [bc] q^y -{- [ce] q^dz [c dO] Çj = 0 . . . etc .
Supposons que nous voulions trouver en premier lieu dœ : il faudra assigner aux coefficients d'élimination des valeurs telles, que la somme des coefficients de dx soit égale à l'unité, et les sommes analogues égales à zéro, pour les coefficients de dy, dz.... Cela nous donne les i équations de e termes chacune :
[aa] gi + [afc] gj 4- [ac] çs ... = < (q'v) ... [a*] g, + [frft] q, + [bc] q, ... = 0
lac] qi + \bc] q^ + [ce] qs ... = 0 ... etc.
Nous les nommons équations du poids pour la valeur de œ. Éliminant par une méthode quelconque entre ces équations, nous pourrons trouver les valeurs de Qx, q^, q^, qui, substituées dans (q'"), donneront par voie d'addition
dx = [a dO] qi + \b dO] 72 + [c dO] q%...,
c'est-à-dire une équation entre dœ, et les différents dO accompagnés de coefficients numériques. Développons-la de manière à séparer les c?0, nous aurons
dx= (0^1+ ^?s + cqz-[- ...)dO (qv) ... -f {a'qi + h'q^ -f c'q^ + . . .) dO'
+ (a"qi + V'q^ + c"qz + ...) dO" ... etc.
Comparant cette équation avec la première du groupe (q), nous en déduisons
ai = aqi -\- hqi + cqz + ...
(q^'j ... aj = a'qi + t'çj -4- c'qz + ...
as = a"gi + b"qi -\- c"qz + ...
Par suite, sans calculer les valeurs particulières des a, on peut, avec la plus grande facilité, former [a^J dont nous avons uniquement besoin. Il suffit, en effet, de multiplier la première des équations (q^') par aj, la seconde par a^, la troisième par oLi... et d'ajouter. On obtient ainsi
(qvii) [a2] = [aa] qi + [afe] qi + [ac] qz ...
Cette valeur de [aj] est bien plus simple encore qu elle ne le paraît ici ; car nous allons faire voir que le coefficient de qy est égal à l'unité, et que ceux de g^, ^3... sont égaux à zéro.
En effet, nous avons déjà fait remarquer que les e7'reu7^s A peuvent être considérées comme des diffé-
rentielles des quantités observées : A est donc du même ordre que dO, et la difiPérentielle o?A peut être négligée vis-à-vis àe dO, comme étant une quantité du second ordre. Cette considération réduit nos équations (q') à la
forme
dO = a dx -\- b dy -{- c dz -\- ... (Q')'.-- dO' =a' dx-\-h' dy -\-c' dz -\- ...
dO" = a"dx + V'dy + (f'dz -\- ... etc.
Multipliant la 1" par a^, la 2^ par ag, la 3® par ag.... et ajoutant les e équations résultantes, le premier membre ainsi obtenu sera identique avec le second membre de la l^des équations (q); par suite, le second membre de l'équation résultante devra être égal à dx; donc le facteur de dx, dans ce second membre, devra être égal à l'unité, et les facteurs de dy, dz... égaux à zéro.
Ainsi on a
[aa] = l; [aft] = 0 ; [ac] = 0.
Notre équation (q^'j devient donc enfin
Il est clair, d'ailleurs, qu'en agissant d'une manière analogue, on trouverait pour les poids de y, z....
Remarquons que la recherche des poids des inconnues, dont nous venons d'expliquer assez longuement la théorie, exige dans la pratique peu de calculs nou-
veaux ; car les coefficients des équations du poids sont les mêmes que ceux des équations normales ; et on n'a d'autre travail que celui de l'élimination, pour trouver Çi, ou g'2 ou ^'z,..
Tous les développements qui précèdent se résument donc dans cette règle pratique bien simple :
« Pour trouver le poids P„ d'une inconnue, œ, on « remplace, dans l'équation normale relative à a?, la « quantité toute connue par — 1 ; dans toutes les autres « équations normales, on remplace la quantité toute « connue par zéro. — En même temps, on met q^ à la «< place de a;, et on élimine toutes les inconnues res- « tantes. » On trouve ainsi
Connaissant le poids d'une inconnue, on trouvera son erreur moyenne par la relation
yp
Dans le cas de trois inconnues, les équations du poids présentent donc les trois systèmes suivants, dans lesquels Qi, Qj, Q3 sont les véritables inconnues à déterminer, et les q des coefficients arbitraires à éliminer.
[ 0 = — 1 + [aa] Qi + [aè] çj + \ac] qs
(1) ) 0 = 0 4- [a*] Qi + m Qi + Ibc} qz
0 = 0 + lac] Q, + [bc] qt + [ce] ?3.
CALCUL DU POfDà DES INCOiNNUES. 381
( 0 = 0 + [aa] g/ + [ah-] Q, + [oc] q^'
(II) 0 = ~ 1 + [«fr] î/ + \bU] Q. + [èc] q^'
( 0 = 0 + [ac] ?/ + \bc] Q, 4- M q^ •
10= 0 4- [««] 9i" + [a^] 9î" 4- lac-] Qs
(III) ) 0 = 0 + [afc] î/' + [W] g/' + [k] Qî ( 0 = - 1 + [ac] qr + [k] îi" + [c-c] Q,.
Il est facile de voir que ces équations sont les mêmes que celles que nous avons déjà établies, dans un cas de substitution analogue, au paragraphe précédent. Si, en effet, dans les équations (p^'') on remplace les m par leur valeur en q, on retrouve identiquement les équations du poids.
Remarque. — Les équations du poids ne différant des équations normales que par les termes tout connus que renferment ces dernières, il s'ensuit que si, dans chacun de ces deux systèmes, on éliminait successivement toutes les inconnues sauf une seule, œ, le coefficient de X serait identiquement le même dans les deux équations finales : or, d'après ce que viennent de nous montrer les équations du poids, ce coefficient est précisément le poids de œ.
Il suit de là qu'on peut se dispenser d'opérer sur les équations du poids, et éliminer simplement entre les équations normales (en ayant soin de n'introduire aucun facteur d'élimination). Le coefficient de l'inconnue qui reste la dernière est le poids de cette inconnue. Ainsi, le poids de v (§ 129) est [dd. 3] ; et pour obtenir ceux de z, y, X, il suffirait de former les coefficients [ce. 3], [hh. 3], [aa. 3].
Cette méthode, qui est celle de Gauss [Theor. motus corpor. cœL, § 182), fournit donc à la fois les inconnues et leurs poids. Cette particularité peut s'expliquer, en dehors de tout calcul, par la remarque que les coefficients de toutes les inconnues dans chacune des équations normales sont formés symétriquement au moyen de ceux des équations primitives. Il en résulte qu'après l'élimination, la dernière inconnue restante est déterminée par une équation M a; = N , ou M' 2/ = N', ou ... etc., dans laquelle M et N, M' et N', etc., sont aussi symétriques l'un par rapport à l'autre, et ne diffèrent que par la substitution des coefficients àa y k ceux de x ou inversement. Cette relation de symétrie est évidemment en harmonie avec ce que nous nommons le poids d'une valeur par rapport à une autre. ^.
§ 134. — Comme application, calculons les poids des inconnues dont nous avons déterminé (§ 131) les valeurs les plus probables.
Les équations du poids de œ sont ici 27Qi 4- 6îî — 1=0
qi 4- 54^3 ==0.
On en déduit par l'élimination ordinaire
19899 d'où Pi =24,596 ... poids de X.
On trouverait, par une marche analogue,
Ps = 13,649, ... poids de j/.
Ps = 53,927, ... — 2.
Cette dernière valeur n'est autre chose que celle du facteur [ce. 2].
Quant à l'erreur moyenne de chaque inconnue, elle se calculera en fonction de son poids, et de l'erreur moyenne (§ 132) de l'unité de poids.
On aura donc (k"")
e' = 0,0S7 erreur moyenne de x s" = 0,077 — y
§ 135. — Ajoutons quelques exemples pour familiariser le lecteur avec l'application de la méthode des moindres carrés.
1° On doit faire passer une droite par quatre points dont les coordonnées n'ont pu être mesurées qu'approximativement : quelle est l'équation la plus probable de la droite ?
Les coordonnées approximatives des quatre points
sont
Puisque la ligne est droite, elle doit satisfaire à l'équation
Y = aX + P ,
dans laquelle a et [3 sont les inconnues à déterminer. La mettant sous la forme
0 = aX + ? — Y ,
et la comparant à la forme générale • ^ = ax -\- hxj ■\- ... -\- n ,
384 THÉORIE DES PROBADILITÉS.
on a pour équations normales :
[XX1a + [X.i]P-[XY] = 0, [X]a-f [1] P- [Y] =0.
Substituant aux quatre groupes de coordonnées leurs
valeurs observées, on trouve
