Расчет видимых координат Луны и Солнца.

[Yaropolkow]

Здравствуйте, дорогие небожильцы. Подскажите где взять методику расчетов положения Луны и Солнца в нужный момент времени. Только никаких программ калькуляторов ит.д
Порекомендуйте литературу, учебник, статью, ссылку. Буду благодарен. Очень. Желательно, что б самому эфемериды-считать.
Астрономия с какулятором не-подходит, нет времени еще и вникать-в програмирования калькуляторов.

Нужны просто формулы.И последовательность расчетов.
Если нужны стронние данные например эферы, то подскмжите где их брать на просторах интернета.-а то эферы я находил только до 10 года, а другие предоставляют инф. Только с сайта, что тоже не подходит.

[Pluto]

http://stjarnhimlen.se/comp/ppcomp.html

[Astro_Coder]

В этой книге есть - Монтенбрук, Пфлегер "Астрономия на персональном компьютере".
+ исходники для ПК на С++ и Паскале.

[Toth]

Если нужны стронние данные например эферы, то подскмжите где их брать на просторах интернета

Для проверки можно вот тут например - http://vo.imcce.fr/webservices/miriade/?forms
PS Если нужны видимые, то надо будет учитывать рефракцию в атмосфере, которая зависит от погоды. В IMCCE есть наблюдаемые ( Apparent ).

[Yaropolkow]

Всем привет!
И так начал изучать методику данную Pluto тут (  http://stjarnhimlen.se/comp/ppcomp.html)


И так вырезка - перевод от туда. И прошу помочь разобратся.

"Этот «номер дня» d вычисляется следующим образом (y = год, m = месяц,  D = дата, UT = UT в часах + десятичные знаки):
d= 367 * y - 7 * (y + (m + 9) / 12) / 4 + 275 * m / 9 + D - 730530"

Коротко - суть проблемы в том, как правильно считать, что на что делить и когда.
Как у Задорнова:
2+2*2
Можно получить 6  = 2+(2*2)
А можно восемь  8 = (2+2)*2
(да такой я неуч)



Кому не лень и кто правильно понимает эту формулу перепроерте меня пожалуйста на дату

29.09.2018 у меня вышло число - 6875,875

Меня смущает несколько моентов.
1. Что как и где делить например
 /4 + 275 * m/   - делить на 4 или на 4+275 и т.д.

2. Его пояснения

"Обратите внимание, что здесь ВСЕ разделы должны быть подразделениями INTEGER. В Паскале используйте «div» вместо «/», в MS-Basic используйте «\» вместо «/». В Fortran C и C ++ «/» могут использоваться, если оба y и m являются целыми числами. "


Его формулу я просто скопипастил в свою программу.

d= 367 * y - 7 * (y + (m + 9) / 12) / 4 + 275 * m / 9 + D - 730530

Она высчитала на дату 29.09.2018 число  6875,875

Я правильно все подсчитал ?

На всякий случай приложу полный перевод статьи.


Просто Pluto подскажет может или кто занимался этим.

Скажу тоже - я пишу свою программу. Но она имеет исключительно свои цели и задачи, поэтому множество программ готовых мне не подходят.
Приходится свою мудрить.
Это не для астрологии и я не студент, а исключительно в научных целях.




Ну или просто если не правильно рассчитал поставте где надо скобки пожалуйста. С меня математик никакой, сори.

[Toth]

Как я понял, d- число дней от 2000 Jan 0.0 UT.
Значит, у вас должно получиться целое число. Там же написано В Паскале используйте «div» вместо «/». Div - это деление без учета остатка. Например, 7 div 4 получится 1.
И вообще, на чём собираетесь считать ?
В Excel можно так - см. рис.
В Паскаль (VCL Delphi, LCL Lazarus ) есть тип TDateTime и подпрограммы DecodeDateTime, EncodeDateTime. С ними просто.

[Yaropolkow]

Спасибо за ответ.
Я уже пишу на C#
Но я и запнулся на этом вопросе. Вот у вас уже другое число вышло.
Я использовал тип Double то, есть с запятыми.
Дальше и не могу продолжать, потому что накопится потом ряд ошибок и я так и не смогу рассчитать координаты на нужную дату.
С уважением Ярослав.


P/S

А понял   "d- число дней от 2000 Jan 0.0 UT."

[xd]

Целочисленное деление должно быть.

[Toth]

Я уже пишу на C#

Вот, посмотрите https://metanit.com/sharp/tutorial/19.1.php

[Yaropolkow]

Спасибо за ответ. Этим я и занимался как раз. С этим разобрался проблем нет.

Проблема то просто в этом, как правильно работать с этой формулой. Что сперва умножать или делить.
Может Pluto ответит раз ссылку дал. как считать. 

видимо придется дальше считать сравнивая (если ни кто не подскажет) с уже проверенными программами, Орбитрон на пример. И СтарКальк и т.д.

[Yaropolkow]

Целочисленное деление должно быть.

Хороше. Int -ами буду считать (целыми).
Спасибо большее всем. Завтра перепишу.
У меня там будет еще ряд вопросов.
Взамен выложу программу которая будет считать гео-координаты где Солнце и Луна на какой либо период времени находятся в Зените.
Кому это понадобится.... уже другой вопрос )))
(я не плоско земельщик, а то вдруг подумаете ))))
С уважением Ярослав.

[xd]

Умножение и деление - левоассоциативные операции равного приоритета. Где надо поднять приоритет, есть скобки.

[Toth]

гео-координаты где Солнце и Луна на какой либо период времени находятся в Зените

Это называется географическое место светила. Солнца или Луны например.

[Yaropolkow]

Все сходится теперь. 6848 на сегодня или 6847 на вчера. Если считать целыми. Если это правильно (а я вам доверяю больше чем себе в этом вопросе) , то все должно работать в будущем.
Проект не коммерческий.
Спасибо, продолжаю....

[Yaropolkow]

Ну вот и приплыли
Считаем наклон эклиптики:

Другая необходимая нам величина - это ecl, наклон эклиптики, т. Е. «Наклон» оси вращения Земли (в настоящее время 23,4 градуса и медленно уменьшающийся). Сначала вычислите «d» ( 6848 на сегодня) момента интереса (раздел 3). Затем вычислим наклон эклиптики:
    ecl = 23,4393 - 3,563E-7 * d

Проблема в том, что я получаю одно и тоже число 23.436....
На любую дату что в Январе, что в Июне, что в марте - а ведь так не должно быть как я понимаю

вот как я работал с его формулой

Сперва я преобразовал число  3,563E-7 в нормально назовем его g у меня получилось

   3.563 >> g1 = 0.0000003563                      --- Это правильно?? 

Потом
g1 * d
 0.0000003563 * 6848
ну а потом
23,4393 - (0.0000003563 * 6848) и поллучаем практически вечное 23.4369

Что не так? Подскажите плииз

[Папа и Сын]

Ну вот и приплыли
Считаем наклон эклиптики:

Другая необходимая нам величина - это ecl, наклон эклиптики, т. Е. «Наклон» оси вращения Земли (в настоящее время 23,4 градуса и медленно уменьшающийся). Сначала вычислите «d» ( 6848 на сегодня) момента интереса (раздел 3). Затем вычислим наклон эклиптики:
    ecl = 23,4393 - 3,563E-7 * d

Проблема в том, что я получаю одно и тоже число 23.436....
На любую дату что в Январе, что в Июне, что в марте - а ведь так не должно быть как я понимаю

вот как я работал с его формулой

Сперва я преобразовал число  3,563E-7 в нормально назовем его g у меня получилось

   3.563 >> g1 = 0.0000003563                      --- Это правильно?? 

Потом
g1 * d
 0.0000003563 * 6848
ну а потом
23,4393 - (0.0000003563 * 6848) и поллучаем практически вечное 23.4369

Что не так? Подскажите плииз

Все правильно, скорость изменения угла наклона эклиптики порядка нескольких десятков тысяч лет (можно нагуглить картинку, как перемещается ось вращения - направление на север - по небу, а тут есть желание увидеть это изменение на нескольких десятках дней  . Если бы небо скакало с такой скоростью, то, боюсь, древние ни календарь, ни астрономию бы не придумали

[Yaropolkow]

Спасибо ВсЕМ за ответы.

Но вот идя по этому пособию я встал в тупик.

И так там есть некая
v = истинная аномалия (угол между положением и перигелием)

Далее дается формула для расчета Солнца

"Затем вычислите расстояние Солнца r и его истинную аномалию v:
    xv = r * cos (v) = cos (E) - e
    yv = r * sin (v) = sqrt (1.0 - e * e) * sin (E)
"


Обратите внимание что уже в этой формуле V используется.

Но чуть ниже написанно как ее вычислить.

    v = atan2 (yv, xv)


Вот и получается что б вычилить расстояние до Солнца нужно V

А что б вычислить V нужно вычислить YV и XV

Круг замкнулся.....

Что делать?


Прикрепляю ниже перевод этой стать. Пункт расчет Солнца

[Yaropolkow]

Вот

[AN14440]

16 лет назад на этом Форуме уважаемые Karen (до сих пор действующий модератор) и Анатолий Волчков (почетный член форума, которого с нами уже нет) довольно терпеливо разъяснили мне, дилетанту, аналогичную(ые) задачу(и):
https://astronomy.ru/forum/index.php/topic,1128.0.html
Там от Анатолия вроде были готовые обкатанные процедуры на Паскале, очень сильно помогающие в написании программ.
Там же рядом аналогичная тема про Луну: https://astronomy.ru/forum/index.php?board=9;action=display;threadid=1176
---------
Ужас! Как давно это было...

[xd]

Рекомендую что ли почитать что-нибудь по теме... Чтобы кроме математики было ещё и представление о физике явления.
Хорошее чтиво - Астрономический календарь, постоянная часть. Там первый раздел начинается с отличного введения в сферическую астрономию и небесную механику.