П. В. Артемий
|
|
Помогите или объясните.
|
26.02.2008 16:49 |
|
|
Пишу диплом по мореходной астрономии, дошёл до части где необходимо написать самому программу для
вычисления положения планет (Венера, Марс, Юпитер, ,Сатурн). Но возникла проблема.
В книге Миуса "Астрономические алгоритмы", есть вычисления гелиоцентрических широт, долгот и радиус-
векторов планет.
Там имеется такая формула:
L = (L0 + L1*t +L2*t*t +L3*t*t*t +L4*t*t*t*t +L5*t*t*t*t*t)/100 000 000 (рад)
В конце имеются приложения где есть все эти коэффициенты (L0,L1,L2,L3,L4,L5), но я не могу понять по
какому принципу их выбирать и для чего дается формула
L0=A*Cos(B+C*t)
L0: A B C
1 1.0564 3.6933 89624555.3
2 1803.0 4.1033 5661.3320
3 . . .
3 . . .
4 . . .
5 . . .
и т.д
Так понимаю все эти коэф-ты подставляются в формулу L0=A*Cos(B+C*t) т.е.:
L0=1.0564 *cos(3.6933+89624555.3*t)
+1803.0*cos(4.1033+5661.3320*t)
+........ и т.д
Это так?
|
|
В. В. Чазов
|
|
Re: Помогите или объясните.
|
27.02.2008 9:10 |
|
Уважаемый коллега. Вы всё правильно написали. В моделях движения планет присутствуют
многочлены по времени (коэффициенты L0,L1,...), но есть также и периодические изменения
долготы (нулевого коэффициента L0), широты и расстояния от Солнца. Периодические изменения
надо учитывать по приводимой Вами формуле.
Особенного физического смысла в этих числовых коэффициентах нет, особенно в коэффициентах
многочленов по степеням времени. Так принято для удобства вычислений: формулы простые,
короткие и обеспечивают приличную точность на интервалах времени 100-200 лет.
С поклоном, Вадим.
|
|
