CONSIDÉRATIONS SUR LES PERTURBATIONS PLANETAIRES
de ce Volume, que la fonction perturbatrice R peut en général être développée
on une série convercçente de la forme
> A rfi'-C l-«*»i
Lt^-u^Ai *'V « (L€ /'^u-*'-'"
R=:2]CcosD,
ou 1 on a
D —i{nt^t) -hi'{n't + £') + km + k'xs' +je +jB' \
{, i', k, k',j etj' sont des nombres entiers quelconques, positifs, nuls ou négatifs. Les coefficients G sont des fonctions de a, a', e, e' , ^ et o', qui diminuent en général assez rapidement quand les valeurs absolues des nombres entiers i, i' , k, k' , j, f augmentent.
Dans l'expression (4) devraient figurer aussi des termes analogues à ceux que nous avons mis en évidence, et dans lesquels n' , t' , ... seraient remplacés par n" , £", .... On voit bien ainsi de quelle manière entrent les divers éléments des planètes P, P', P", . . . dans le développement de R.
Les dérivées partielles de R par rapport à l'un quelconque des cinq éléments c, ^, 0, m, £ seront exprimées par des développements de même forme que (4). sauf que les cosinus pourront être remplacés par des sinus.
Il en va tout autrement de la sixième dérivée partielle -v-j elle se compose, en
da
effet, de deux parties : la première, que nous représenterons par (-p )' s'obtient
en faisant varier a seulement dans les coefficients G; la seconde provient de la variation de a dans n sous les signes cosinus. D'après la formule (5), les arguments D dépendent der-t, par suite de a, d'après la relation
On aura donc
ii'-a} :=!. f(i H- m).
ail /^\ ^R da dn^ \da J du da
ou bien, en remarquant que n n'entre dans R que par nt qui accompagne toujours £,
da
dK\ ÔK da
da J de da
On trouvera ainsi, en se reportant aux formules (4) et (5),
dR da
dC da
cos
da
D-^|^2^'Cs'»I>-
On voit que le temps t est sorti des signes cosinus; de là un grave inconvé-
de
nient que présenterait l'emploi de la valeur (2) de -j-^ la seule des six dérivées
192 CnAPITPxE XI.
des éléments qui contienne j-- Malgré la petitesse du facteur m' qui entre dans
le coefficient G, le terme G^sinD pourrait prendre des valeurs très grandes, et serait gênant de toutes façons. Voilà l'inconvénient dont on a parlé; on l'évite comme il suit :
Si l'on a égard à la valeur (6) de -r-> la dernière des formules (2) donne, en n écrivant pas, pour abréger, les termes en -p et -j-^
de 2 /(?R\ 2 dK dn
^ dt na\da ) na dt da
Or on a
dn dn da '>. (?R dn
dt da dt na de da
ce qui permet d'écrire comme il suit l'équation (7),
de dn 2 /dR
dt dt na \da
On est ainsi conduit à prendre, au lieu de £, un nouvel élément i^*\ tel que l'on ait
d£^^'> dt dn
(^^ ~dr~di^^Tt'
On trouvera immédiatement, en écrivant maintenant les termes en -v- et -c-j
ae ocf)
\\f ' rit ' n/^\ An I 9. /~~ I9 r)m '
dt na\da) nà^Ji — e'^ ^9 '^^'^ ^^
Or on tire de (8)
e(') — £ -{- I t -T- dt =^ £ -\- nt — I ndt,
(lo) nt -h £ =z I ndt -+- £''^K
On voit donc que le changement de variable sera bien facile à faire, puisqu'il se bornera à remplacer dans les expressions (i) de oo, y, z, nt -]- z = l par
I ndt -h £''^ La formule (lo) montre d'ailleurs que l'on aura
dR dR
PERTURBATIONS DES DIVERS ORDRES. IqS
Si l'on remplace dans la première, la quatrième et la cinquième des formules (2), -p- par ^-^j et si on les rapproche ensuite de (9), on voit que les nouvelles équations différentielles ne différeront des anciennes qu'en ce que £ et -r- auront été remplacés respectivement par £^'' et (-j-j- Il convient, pour ne pas multiplier les notations, de supprimer l'indice de £^'' et la parenthèse de (-p); cela permettra de conserver les équations (2) sous leur première forme. Seulement il sera bien entendu que la dérivée -p y devra être prise sans faire varier a sous les signes cosinus, et que, dans les formules (i), qui donnent X, y, z, nt -\- 1 devra être remplacé par / ndt -\- i.
Nous ferons, pour abréger.
ndt = p, d'où n=~^;
quand n sera connu, en effectuant la quadrature (ndt, nous n'ajouterons pas de constante d'intégration, parce qu'elle irait se fondre avec £ qui accompagne toujours In dt. Enfin nous ferons remarquer qu'ici, comme partout ailleurs, la
lettre /i n'a d'autre sens que celui qui est défini par la formule n = t/ /^^ ; de sorte que l'on a
Pour déterminer le mouvement de la planète P', il y a lieu de considérer des équations toutes pareilles à (2), qu'on obtiendra en accentuant les lettres, et mettant au lieu de R la fonction perturbatrice R'. On devra former la dérivée
-yr sans faire varier a' sous les signes cosinus; mais, dans les formules qui donnent ^', y, z' en fonction de t, a, 0', ..., il faudra remplacer n't-\-t' par in'dt-hi'; nous poserons aussi In'dt^p'. La considération des équations différentielles en -1-, -7-, •• • est indispensable, même pour déterminer le mouvement de la planète P, quand on va au delà de la première approximation.
71. Pour fixer les idées, ne considérons que deux planètes P et P'; nous aurons à intégrer par approximations successives le système des douze équations
différentielles simultanées suivantes :
l da 2 ^R
l dt na di
(ij) \^ I ^
da' 1 dV^
dt ~ n' a' de'
(12) \ dd' I dR'
dt
«'a'Vi— e"sin9' ^?'
OU 1 on a
(r^^GcosD, (i3) I D = ï(p + £) + «'(?' + £') + /*^^ + ^'^'+y^+./'^'>
1 p ::= I ndt, P'^^ f fi'dt;
R'=2]C'cosD';
C" et D' sont de même forme que C et D.
Nous avons déjà fait observer que les seconds membres des équations (ii) et (12) sont de petites quantités du premier ordre, à cause des facteurs m' et m qui entrent dans les coefficients C et C.
Nous allons chercher à développer les expressions des éléments variables sous la forme
9 = 9, + di9,-hà,9o+...,
«0, ôo, . • . , «„, 6„, . . . sont douze constantes arbitraires dont on trouvera les valeurs numériques en comparant la théorie à l'observation; les quantités représentées d'une manière générale par la caractéristique 0^ sont des fonctions inconnues du temps t, des constantes ci-dessus et des masses m et m'; relativement à ces masses, tous les ^i seront de l'ordre i; ceux des O/ qui se rapportent à la planète P devront s'annuler avec m', et contenir m' en facteur, tandis que pour la planète P' ils auront le facteur m.
PERTURBATIONS DES DIVERS ORDRES. IqS
On mettra ainsi en évidence les quantités §,«o» ^2^0» • • -, ou les perturbations du premier ordre, du second ordre, etc., de l'élément a, et de même pour les autres éléments.
Il s'agit de calculer ces perturbations des divers ordres. Nous poserons aussi
(15) /i = Wq-H ô, /*0+ «^î/to -f-. . .
et nous prendrons
En substituant les valeurs (i4) et (i5) de a et ^ dans la relation
il viendra
«o-H<5i«o+<52«o + - .. = nAi+ -!— 5 4- ^-^ + . .. )
d'où, en égalant de part et d'autre les quantités de même ordre.
^ 3 ôjao
Oi«o = '«0 i
On posera ensuite
(18) po=«gi, dipo = j di/iodt, àipQ= I difi^dt, ...,
et la formule p =/ /i </^ combinée avec la relation (i 5) donnera
(^9) P = Po+^iPo+ «^aPoH-... ;
S,po sera du premier ordre, B^po du second, etc. On aura des formules toutes pareilles pour la planète P'.
Il faut substituer dans les équations diiTérentielles (11) et (12) les expressions (i4), (i5) et (19).
72. Perturbations du premier ordre. — Pour commencer, nous allons faire la substitution indiquée, en ne considérant que les quantités du premier ordre, et négligeant celles du second. On pourra donc, dans les seconds membres des équations (11) et (12) qui sont déjà du premier ordre, remplacer a, 0, . .., a', 0', . . . par «0» ^0» • • •' ^0' ^0' •••» et aussi p et p' i^slt n^^t et n'^t. On trouvera ainsi
(20) Ro = 2] ^oCOSDo,
(21) J)^= i (n^t -\- £0) + i' {n'^t -^ e'J -^ kT^o+ k' ^'0 -^J Oo -hj%,
d^^a^) 2 ÔRq
dt «0^0 ^£0
dàiB^ I ^Ro
Les seconds membres de ces formules sont des fonctions connues de/; on aura donc
0, an := f -^i — at,
n^al \J i — el siiK^o
On est ainsi ramené à des quadratures; il est inutile d'ajouter des constantes
qui, dans les expressions (i4) de a, 6, . . . , iraient se fondre avec a^, Go, Au
point de vue analytique, toutes ces quadratures dépendent d'une seule,
I^Qdt; car on a, par exemple, On aura de même
«0 a^ si
I — e'n^ sin cp' t/ <^9o
Ro et R'o sont des fonctions très compliquées du temps t et des constantes «„» Go, ...; de telle sorte qu'il ne faut pas songer à effectuer rigoureusement les quadratures qui figurent dans les formules (a) et (a').
On pourrait bien avoir recours aux quadratures mécaniques; c'est ce qu'on
