Расчет движения планет.

[NiOl]

Хочется сделать свою примитивную программку (скрин-савер), где бы наглядно крутились основные "камушки" Солнечной системы.

Достаточно-ли для расчетов использовать закон Всемирного тяготения? (без учета торможения вращения вокруг своей оси спутниками и т.д.).  Как я понимаю, можно воспользоваться обычной формулой "F=G*(m1*m2)/r^2" для расчета ускорения, которым модифицировать значения dX,dY,dZ, а далее - X,Y,Z. Только нужно знать некие начальные значения dX,dY,dZ,X,Y,Z, а также массу каждого участника "движения".

Где можно взять подобную информацию о положении и характере движения планет?

Если у кого есть e-теория и нежалко поделиться - кинте, плиз, на мыло: ll8vl@fromru.com

[kis]

Тут проблема не в формулах а в том что если на экран засунуть дальние планеты Уран Нептун Плутон то внутренние планеты будут буквально в одной точке. Ну и собственно говоря видимый на экране размер планеты (так что бы это было симпатично ) съест любые погрешности расчетов. Можно даже не мудрствуя просто посчитать по периодам обращения смещения за нужный период времени и все.

[Анатолий Волчков]

NiOl: «Достаточно-ли для расчетов использовать закон Всемирного тяготения»?

Все зависит от необходимой точности. Если ограничиться точностью в секунду дуги – эффекты Теории относительности можно не учитывать.

Фактически Вы говорите о так называемом численном интегрировании орбит планет. Суть такого метода в том, что если в начальный момент Т известны положения и скорости планет, можно рассчитать грпвитационное влияние всех планет на каждую и вычислить ускорения всех планет. Считая, что на протяжении некоторого интервала времени dT ускорение не меняется, можно найти положения планет в момент Т+dT. Далее процесс повторяется.
В зависимости от требуемой точности координат выбирается интервал dT.
Размеры планет тут роли не играют, так как приливные процессы существенны лишь для Луны.

По такому примерно принципу построена теория движения планет DE200/LE200. Теория эта содержит полиномы Чебышева, позволяющие на некотором интервале с точностью 0.001 сек дуги найти положения и скорости планет, Солнца и Луны. Все это вычисляется относительно центра масс солнечной системы (барицентр).

Если желаете, можете взять файл на сайте SIMFOV.RU в разделе SOFT. Там этот файл в двоичной системе, охватывает интервал 1975-2025 годы (в распакованном виде это 5.4 Мб). Форматы, если будете использовать файл, опишу.

[panov]

1) Для начала пусть NiOl напишет сверх-примитивный "скрин-савер", чтобы, к примеру, вертелся хотя бы 1 pixel. Нечто вроде:
a=0: while true: x=r*cos(a): y=r*sin(a): a=a+0.01: if a>=6.28 then a=0: wend
Затем уже можно будет этот "савер" наворачивать по Всемирному тяготению.

           "The orbit of  any  one  planet  depends  on  the
           combined  motion  of  all  the  planets,  not  to
           mention the action of all of these on each other.
           But  to  consider simultaneously all these causes
           of motion and to define these  motions  by  exact
           laws   allowing   of   conventional   calculation
           exceeds, unless I am mistaken, the force  of  the
           entire human intellect."
                      -- Sir Isaac Newton, Principia

Но можно, сделав энное количество упрощений, и одному программисту справиться с задачей, если почитать теорию, изложенную в:

A.K.   Dewdney,  "The  Armchair Universe", W.  H.  Freeman: New York, 1988.
        Chapter  "A  Cosmic  Ballet",  pp.   229-238.

2) Кто там оформляет голосования? Нельзя ли поаккуратней с языком!