ВНИМАНИЕ! На форуме началось голосование в конкурсе - астрофотография месяца СЕНТЯБРЬ!
0 Пользователей и 1 Гость просматривают эту тему.
гравитационную постоянную на массу не умножать.
GM Солнца + GM Марса = 1.3271244001799E+20 + 4.2828314258067E+13
Через 59353325.08 секунд Марс должен вернуться на прежнее место.
GM Солнца + GM Марса = 1.3271244001799E+20 + 4.2828314258067E+13Нач. вектор положения Марса =(206648658093.026;0;0), скорость =(0;26498.9066619695;0) - в перигелии. У Солнца - все 0.Через 59353325.08 секунд Марс должен вернуться на прежнее место. В гелиоцентрических координатах ( то есть радиус-вектор Марса минус радиус-вектор Солнца. Хотя Марс легкий, но все же .. )
Так вы что, собираетесь разоблачать гнусных буржуазных лжеученых.Нет, здесь я вам не помощник, с моими программными поделками я не в силах ни опровергнуть, ни подтвердить. А если хотите просто проверить качество моих расчетов, то пожалуйста - задайте исходные данные , просчитаю на какой-то период, потом назад, сравним с исходными и т.д.Исходные данные - в виде векторов скорости и положения и масс ( или произведений GM ). Т е. (X,Y,Z,vX,vY,vZ)1, (X,Y,Z,vX,vY,vZ)2, ..., GM1,GM2.
Для этого воспользуемся моей программой Solsys7m5 (скачать вместе с исходниками и описанием можно на моем сайте - смотрите адрес в подписи)
но и то на каком языке написана ваша программа.
Вам надо будет скачать библиотеку MSVBVM60.DLL, а затем поместить ее в папку C: \ WINDOWS \ SYSTEM32.
А на сайте -Вам надо будет скачать библиотеку MSVBVM60.DLL, а затем поместить ее в папку C: \ WINDOWS \ SYSTEM32. Потом можете скачивать только заархивированные исполняемые файлы нужных Вам программ, например, Solsys4exe.zip или Solsys5exe.rar..Мне как-то не хочется вставлять какую-то .dll в папку C: \ WINDOWS \ SYSTEM32. А потом чего-то запускать. Зачем такие сложности.
Сделайте это сами. И выдайте здесь - всего 10 строк - массы ( или G*M ) Солнца, планет и Луны в каких-либо единицах . И 10 строк vPos,vVel ( по 6 чисел каждая строка ).
Можно не 10, можно какую-то упрощенную систему.
будем ближе к реальности и использовать все 10-ть масс.
Вот - в пред. посте результат . 10 масс. Какие исходные данные брать - не принципиально. Думаю ваши не сильно отличаются от DE405.
а здесь эфемериды DE405 не просто дают какую-то погрешность, а вообще не лезут ни в какие ворота.
Начнем с того, что у вас немного другая исходная система, т.к. вы вместо 10-го тела (Плутон) используете Луну,
Например, по методу Эрмита у вас погрешность расчета для Земли получилась больше, чем для Вернеры, чего не должно быть, а для Марса у вас получилась минимальная погрешность при шаге 0,1, а при меньшем и большем шаге погрешность растет, хотя для остальных внутренних планет она постоянно увеличивается вместе с увеличением шага решения, т.е. объяснить это накоплением "мусора" при мелких шагах решения нельзя.
И вообще, у вас какие-то подозрительные погрешности расчета.
Поэтому я бы для начала попросил вас выложить где ни будь (или прислать мне на емайл) исходный код вашей программы (или хотя бы только код методов численного интегрирования дифференциальных уравнений и самих диффуров)
Вы только JPL "опровергли" или и других тоже - вот например, наш ИПА РАН - http://iaaras.ru/dept/ephemeris/epm/
Плутон совсем легкий. Ну, можно было вместо Земли и Луны взять одно тело - Земля+Луна, и считать движение его барицентра.
Во первых - этот метод " вперед-назад" определяет только точность возвращения в исходные координаты. Он не определяет, насколько точные были пути туда и сюда. Главное, чтобы они совпали.
Потом - погрешности с течением времени растут не монотонно, а с колебаниями, что естественно для замкнутых орбит.И при малых шагах , я думаю начинает сказываться уже т.н. машинный эпсилон.
Для иллюстрации этих парадоксов - см. рис.Первый - СС методом "вперед-назад" на 2 года. По оси X - шаг, в днях, забыл подписать.Второй - система 4 тел Солнце, Земля, Луна и недавний близкий астероид 2021 RS2 , на 5 дней, но другим методом - сравнение с неким шагом и с половинным шагом. Красным пунктиром - момент сближения. Довольно близко, 20 тыс км от геоцентра, поэтому шаги такие мелкие брал. Там по оси X - время от начала.
Посмотрите в этой теме Ответ #8 - я все делал, как в том источнике.
тем более, что планеты движутся не по замкнутым орбитам, т.к. движется вся Солнечная система
а на несколько тысяч лет получается полный бардак. Вот поэтому то и моя теория Ser0 и теория Птолемея ... дают приличные результаты для очень старых наблюдений,
Хотя орбиты строго говоря не замкнутые, но не потому что система движется, они незамкнутые ( не совсем эллипсы) и в гелио- и барицентрических координатах - изза взаимных возмущений.
Вы как это проверяли? По каким-то древнегреческим наблюдениям ?