Барицентрическое движение Солнца

[граф мошкин]

Сейчас читаю про барицентрическое движение Солнца (движение относительно центра масс солнечной системы).
Интересные выводы в этой статье. Вычислены циклы, время которых можно сопоставить с фиктивными и др. точками.
Хотелось бы проверить кое-что, но для этого нужно сопоставить координаты ЦМСС ,
которые есть в разделе "разное" , со множеством карт.
Есть ли возможность вывести их на круг карты?
Задаю этот вопрос здесь, т.к. не нашёл подходящей темы .
Кому интересно, вот ещё на эту тему

[Юрий_1971]

Насколько я знаю, Zet не рассчитывает барицентрические позиции небесных объектов. Но сами швейцарские эфемериды это позволяют, с некоторыми оговорками.
Я изменил один из своих Python скриптов для этой цели:

Код: Выделить всё

#!/usr/bin/env python3
# -*- coding: utf-8 -*-


import swisseph
import math


def get_sky_object(jd, sky_object):
    """ Функция получает небесный объект и возвращает его эклиптикальную
    долготу. """
    result_calc = swisseph.calc_ut(
        jd,
        sky_object,
        swisseph.FLG_BARYCTR
    )
    return result_calc[0][0]


def convert_result(value):
    """ Функция преобразует переданное значение градусов с плавающей точкой в
    градусы, минуты, секунды, доли секунд. И возвращает результат как
    строку. """
    degrees = math.floor(value)
    remainder = value - degrees
    minutes = remainder * 3600 / 60
    remainder = minutes - math.floor(minutes)
    seconds = remainder * 60
    remainder = seconds - math.floor(seconds)
    part_seconds = remainder * 100
    return "{0:>3}\xB0{1:>2}'{2:>2}\".{3:>2}".format(
        math.floor(degrees),
        math.floor(minutes),
        math.floor(seconds),
        math.floor(part_seconds)
        )


if __name__ == "__main__":
    # Дата события.
    date_event = "2.3.2021"
    # Время события по Гринвичу.
    time_event = "12:00:00"
    # Путь к швейцарским эфемеридам.
    swisseph.set_ephe_path("/home/yuri/Exclude/swisseph/")
    date_event = date_event.split('.')
    time_event = time_event.split(':')
    flag = swisseph.GREG_CAL
    et_ut = swisseph.utc_to_jd(
        int(date_event[2]),
        int(date_event[1]),
        int(date_event[0]),
        int(time_event[0]),
        int(time_event[1]),
        int(time_event[2]),
        flag
    )
    jd = et_ut[1]
    sun = get_sky_object(jd, swisseph.SUN)
    earth = get_sky_object(jd, swisseph.EARTH)
    moon = get_sky_object(jd, swisseph.MOON)
    merkury = get_sky_object(jd, swisseph.MERCURY)
    venus = get_sky_object(jd, swisseph.VENUS)
    mars = get_sky_object(jd, swisseph.MARS)
    jupiter = get_sky_object(jd, swisseph.JUPITER)
    saturn = get_sky_object(jd, swisseph.SATURN)
    uranus = get_sky_object(jd, swisseph.URANUS)
    neptune = get_sky_object(jd, swisseph.NEPTUNE)
    pluto = get_sky_object(jd, swisseph.PLUTO)
    # Вывод результатов.
    string_date = '.'.join(date_event) + ' ' + ':'.join(time_event) + '\n'
    print(string_date)
    print("Солнце  :", convert_result(sun))
    print("Меркурий:", convert_result(merkury))
    print("Венера  :", convert_result(venus))
    print("Земля   :", convert_result(earth))
    print("Луна    :", convert_result(moon))
    print("Марс    :", convert_result(mars))
    print("Юпитер  :", convert_result(jupiter))
    print("Сатурн  :", convert_result(saturn))
    print("Уран    :", convert_result(uranus))
    print("Нептун  :", convert_result(neptune))
    print("Плутон  :", convert_result(pluto))

Результат на сегодня, на 12 часов по Гринвичу:

Код: Выделить всё

2.3.2021 12:00:00

Солнце  : 141°39'10".29
Меркурий: 223°53'31".86
Венера  : 328°10'22". 5
Земля   : 161°56'30".60
Луна    : 162° 2'11".25
Марс    :  96°49'18".16
Юпитер  : 312°10'30".61
Сатурн  : 305°28'28".14
Уран    :  40° 3'16".34
Нептун  : 350°35'41".42
Плутон  : 294°52' 8".46

Можете по аналогии воспользоваться Экселем, или чем там вам удобно.

[граф мошкин]

Насколько я знаю, Zet не рассчитывает барицентрические позиции небесных объектов.
""" Функция получает небесный объект и возвращает его эклиптикальную
долготу. """

""" Функция преобразует переданное значение градусов с плавающей точкой в
градусы, минуты, секунды, доли секунд. И возвращает результат как
строку. """


Можете по аналогии воспользоваться Экселем, или чем там вам удобно.

Для меня это тёмный лес, к сожалению..
Я думал, что если уже существует в ZET диаграмма, показывающая координаты барицентра на эклиптике (ЦМСС) ,
то можно их вывести в виде значка на круге карты.
Правильно ли я понял, что в результате Ваших преобразований - искомая точка - это барицентрическая позиция Солнца ?

[Юрий_1971]

Правильно ли я понял, что в результате Ваших преобразований - искомая точка - это барицентрическая позиция Солнца ?

Это проекции небесных объектов на эклиптику относительно барицентра солнечной системы. Вы же на карте хотели их увидеть, которая обычно и есть "круг" эклиптики.

P.S. Или вы хотели точку барицентра увидеть на эклиптике? Но барицентр находится в самом Солнце, хоть и не в самом его центре, это Zet показывает.

[граф мошкин]

Правильно ли я понял, что в результате Ваших преобразований - искомая точка - это барицентрическая позиция Солнца ?

Это проекции небесных объектов на эклиптику относительно барицентра солнечной системы. Вы же на карте хотели их увидеть, которая обычно и есть "круг" эклиптики.

P.S. Или вы хотели точку барицентра увидеть на эклиптике? Но барицентр находится в самом Солнце, хоть и не в центре самого Солнца, это Zet показывает.

Нет , не всегда в самом Солнце. Иногда и на расстоянии его диаметра.
В общем, это не имеет значения, потому что в любом случае Солнце должно как-то реагировать на движение барицентра.
Мне важна линия, соединяющая центр массы Солнца с барицентром, и её проекция на эклиптику.
Дело в том, что влияние планет земной группы на Солнце с периодом ~11лет - уже хорошо известно - это вспышечеая активность.
Но есть ещё влияние планет-гигантов - оно создаёт другой цикл , выражающийся в движении солнца вокруг барицентра.
Мне интересно сопоставить это движение с событиями и свойствами рождённых в разные фазы этого движения.