Старый 22.11.2009, 15:28   #141
Edik
Собеседник
 
Регистрация: 19.12.2008
Сообщения: 7
Edik репутация выше +10
По умолчанию

Зраствуйте,

Прочитал тему, ну нереально интересная и очень полезная. Только я не нашёл как посчитать деклинации и широту. Думаю такая возможность в ШЕ ест, надо только декларироват в функции PLC ?
LordWilex, можете добавит деклинации ?

Заранее спасибо,

Едуард



If CType = "STrop" Then iflag = SEFLG_TRUEPOS + SEFLG_SWIEPH
If CType = "SSid" Then iflag = SEFLG_SIDEREAL + SEFLG_SWIEPH
If CType = "SHel" Then iflag = SEFLG_HELCTR + SEFLG_SWIEPH
If CType = "SXYZ" Then iflag = SEFLG_XYZ + SEFLG_SWIEPH
If CType = "SRad" Then iflag = SEFLG_RADIANS + SEFLG_SWIEPH
Edik вне форума   Ответить с цитированием
Старый 22.11.2009, 16:14   #142
LordWilex
В отпуске
 
Аватар для LordWilex
 
Регистрация: 01.06.2008
Адрес: Таганрог
Сообщения: 28,983
LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000
По умолчанию

Цитата:
Сообщение от Edik
Зраствуйте,

Прочитал тему, ну нереально интересная и очень полезная. Только я не нашёл как посчитать деклинации и широту. Думаю такая возможность в ШЕ ест, надо только декларироват в функции PLC ?
LordWilex, можете добавит деклинации ?

Заранее спасибо,

Едуард



If CType = "STrop" Then iflag = SEFLG_TRUEPOS + SEFLG_SWIEPH
If CType = "SSid" Then iflag = SEFLG_SIDEREAL + SEFLG_SWIEPH
If CType = "SHel" Then iflag = SEFLG_HELCTR + SEFLG_SWIEPH
If CType = "SXYZ" Then iflag = SEFLG_XYZ + SEFLG_SWIEPH
If CType = "SRad" Then iflag = SEFLG_RADIANS + SEFLG_SWIEPH
Конечно есть, все типы координат выводятся через функции PLC, и ни чего дополнительно декларировать не нужно!!!

Plc(Юлианская_дата, Номер_планеты, Тип_Зодиака, Координата);

Если в последний парраметр (в данном примере я его обозвал "координата"), поставить 0, то функция выдаст обычную долготу, если поставить 1 - то выдаст широту, если 2, - радиус-вектор (расстояние до планеты), если поставить 3,4,5 соответствующие скорости.

Ну это конечно, при условии, что вы выбирите в предпоследнем парраметре, STrop/MTrop/JTrop или SSid/MSid/JSid (эклиптические координаты, сидерический или тропический зодиак).
LordWilex вне форума   Ответить с цитированием
Старый 22.11.2009, 16:18   #143
Edik
Собеседник
 
Регистрация: 19.12.2008
Сообщения: 7
Edik репутация выше +10
По умолчанию

извиняюсь, поспешил ... всё получилось :) большое спасибо LordWilex за программу :)
Edik вне форума   Ответить с цитированием
Старый 23.11.2009, 22:21   #144
nsk-1
Собеседник
 
Регистрация: 23.11.2009
Сообщения: 1
nsk-1 отключил(а) отображение уровня репутации
По умолчанию

А как бы обратиться к этой библиотеке из PHP во FreeBSD?
nsk-1 вне форума   Ответить с цитированием
Старый 23.11.2009, 22:46   #145
LordWilex
В отпуске
 
Аватар для LordWilex
 
Регистрация: 01.06.2008
Адрес: Таганрог
Сообщения: 28,983
LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000
По умолчанию

Цитата:
Сообщение от nsk-1
А как бы обратиться к этой библиотеке из PHP во FreeBSD?
Нет, PHP пока не поддерживается, но есть 2 способа выкрутиться.
1. Если вам без разницы PHP или Perl.

В последней nix'ской версии есть расширение для Perl'a.
Просто компилите под FreeBSD тем компилятором, что у вас есть nix'скую версию ШЭ (на сколько я знаю, под FreeBSD "родной" компилятор тоже gcc ?), и устанавливайте Perl'овское расширение ШЭ себе в систему.
Только для этого нужны права администратора, чего вам вряд ли дадут на обычном хостинге, а я так понимаю вас оно интересует именно в качестве web-приложения?

2. Компилите под FreeBSD программу "swetest" (есть в комплекте никс'овскаой версии ШЭ, или если она вас не устраивает, - пишите свою), и вызываете её через PHP, кажется в PHP для этого есть функция "SysExec()" или что-то в этом роде, не помню точно....
Но для этого обязательно нужно, чтобы ваш хостинг поддерживал запуск скомпилированных исполняемых файлов.

***

Больше ни как
LordWilex вне форума   Ответить с цитированием
Старый 19.12.2009, 12:37   #146
Tom Sher
Собеседник
 
Регистрация: 19.12.2009
Сообщения: 2
Tom Sher отключил(а) отображение уровня репутации
По умолчанию

Здравствуйте!

Скажите пожалуйста, а есть ли возможность решения обратной задачи: вычислить точный момент времени по долготе той или иной планеты? Особенно интересует Солнце. Быть может, кто-то уже написал макрос для решения этой проблемы? Заранее благодарю за ответ...
Tom Sher вне форума   Ответить с цитированием
Старый 19.12.2009, 16:40   #147
LordWilex
В отпуске
 
Аватар для LordWilex
 
Регистрация: 01.06.2008
Адрес: Таганрог
Сообщения: 28,983
LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000
По умолчанию

Цитата:
Сообщение от Tom Sher
Здравствуйте!

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


*****
Кстати, анонс:

В ближайшее время (до нового года) выйдет обновление с исправлением некоторых ошибок, а так же, новым набором функций для работы с управителями/диспозиторами, а еще появятся функции для работы с ГША и жребиями.
LordWilex вне форума   Ответить с цитированием
Старый 19.12.2009, 16:45   #148
Алексей
Собеседник
 
Аватар для Алексей
 
Регистрация: 01.11.2007
Сообщения: 8,521
Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000
По умолчанию

Цитата:
Сообщение от LordWilex
Кстати, анонс:

В ближайшее время (до нового года) выйдет обновление с исправлением некоторых ошибок, а так же, новым набором функций для работы с управителями/диспозиторами, а еще появятся функции для работы с ГША и жребиями.
Вау!!!
Алексей вне форума   Ответить с цитированием
Старый 20.12.2009, 16:22   #149
Tom Sher
Собеседник
 
Регистрация: 19.12.2009
Сообщения: 2
Tom Sher отключил(а) отображение уровня репутации
По умолчанию

Цитата:
Сообщение от LordWilex
Пока нет, но в.. нет не в следующей версии... а через версию или две, конечно же планируется.


*****
Кстати, анонс:

В ближайшее время (до нового года) выйдет обновление с исправлением некоторых ошибок, а так же, новым набором функций для работы с управителями/диспозиторами, а еще появятся функции для работы с ГША и жребиями.

Спасибо за ответ. Пока приходится решать эту проблему через .... подручными средствами Excel. Будет хорошо, если эта опция появится. Тем более, что программка здесь не слишком сложная должна быть...
Tom Sher вне форума   Ответить с цитированием
Старый 20.12.2009, 19:21   #150
LordWilex
В отпуске
 
Аватар для LordWilex
 
Регистрация: 01.06.2008
Адрес: Таганрог
Сообщения: 28,983
LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000
По умолчанию

Цитата:
Сообщение от Tom Sher
Спасибо за ответ. Пока приходится решать эту проблему через .... подручными средствами Excel. Будет хорошо, если эта опция появится. Тем более, что программка здесь не слишком сложная должна быть...

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

В скором времени (через версию или две, как получиться) я планирую координальным образом изменить архитектуру проекта.
Все дополнительные выичсления будут проводиться не в VB, как сейчас, а в отдельном dll модуле, в эксель будут подоваться только окончательные результаты вычислений, а в среде VB будут находиться только декларации функций и ни чего больше.
Причем, часть кода будет реализованна на ассемблере.

Благодаря этим мерам, как я надеюсь, скорость появления результата в экселевских ячейках увеличиться раза в три (в смысле скорости, увы, это все, что можно теоретически выжать из среды эксель), а код внутри самой dll будет выполняться раз в 10 быстрей, а то и больше.
Это позволит не слишком беспокоиться о скорости, и реализовывать внутри этой dll довольно сложные алгоритмы с иттерациями, и вообще моделировать такие эксперименты, которые раньше было реализовать почти невозможно. Вот тогда и сделаем то, что вы хотите.

К тому же, такая архитектура позволит сделать аналогичные модули для mathcad и mathlab, с минимальными издержками, а там возможностей гораздо больше, чем в excel.
LordWilex вне форума   Ответить с цитированием
Старый 29.12.2009, 08:02   #151
LordWilex
В отпуске
 
Аватар для LordWilex
 
Регистрация: 01.06.2008
Адрес: Таганрог
Сообщения: 28,983
LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000
Arrow Вышла новая версия Швейцарских эфемерид по-русски. Ура, товарищи!!!

Швейцарские эфемериды по-русски.
Версия 1.2 (от 31 декабря 2009 года)

[img]http://s4.**********/f9d464fe51269a475f8c33303e9ccb64.gif[/img] [img]http://s4.**********/f9d464fe51269a475f8c33303e9ccb64.gif[/img] [img]http://s4.**********/f9d464fe51269a475f8c33303e9ccb64.gif[/img]

Тока шо вышла в свет новая версия Швейцарских эфемерид по-русски.

А шо новго?

[img]http://s17.**********/7c94993375e6ff8fd7b6d9ed5cb60188.gif[/img]
  • Исправлен ряд ошибок
  • Появились новые функции, которые обеспечивают полноценную работут с системами управителей/диспозиторов (класска и модерн), термов, триплицитетов, а так же функции для работы с планетными картинами ГША и жребиями:

RulerNum() - Возвращает номер планеты (в системе ШЭ), которая управляет, знаком Пример, если ввести “5” (лев), то функция возвратит “0” (Солнце)

RulerName() - Делает все тоже, что и RulerNum, только вместо номера, возвращает английское название планеты-управителя

На пример, если ввести 5 (Лев), то функция возвратит “Sun”

Dispositor() - Возвращает долготу диспозитора или долготу планеты, которая в этом знаке в изгнании.

На пример, если ввести 40 (10 градус Тельца), то функция вернет 3 (Венера), или 4 (Марс), если запросить планету в изгнании.

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

Пример. Допустим, что мы ввели 40 (10й градус Тельца), Тельцом управляет Венера, которая находится в 95 градусе, диспозитор которой Луна в 120 градусе. Если в параметрах задать управителя, то функция выдаст долготу Венеры (95), если задать диспозитора, то функция возвратит Луну (120).

RulerHouse() – Возвращает зодиакальную долготу управителя дома.

Пример. Нас интересует управитель 3 дома. 3й дом в Весах, управитель Весов Венера, функция возвратит долготу Венеры. Так же можно получить номер планеты, которая в этом знаке в изгнании.

RulerHouseInHouse() - Возвращает номер дома, в котором находится управитель дома

Пример. Нас интересует II дом. Второй дом в Весах. управитель II дома Венера в III доме, номер этого дома функция и вернет.

DayHight() – Возвращает 1, если карта дневная, и 0, если карта ночная.

jdayLT() – Делает тоже самое, что и функция jdayLT, т.е. выводит юлианскую дату, только вводить нужно местное время + поправка, а не время по гринвичу. Думаю так многим это облегчит работу.

Triplitsitet() – Возвращает номер управителя по триплицитету зодиакальной координаты Point,, в зависимости от того дневная карта или ночная (классика). Пример. Планета находиться в 10 градусе Стрельца. Карта дневная. Дневной управитель огненного триплицитета Солнца, следовательно его номер функция и выдаст.

TriplitsitetName() – Делает тоже самое, что и функция Triplitsitet (), только вместо номера возвращает латинское название управителя по триплицитету.

Пример. Планета находиться в 10 градусе Стрельца. Карта дневная. Дневной управитель огненного триплицитета Солнца, следовательно его имя (Sun) функция и выдаст.

Term() – Возвращает номер планету управителя термом.

Пример. Нас интересует 8 градус Овна. По Птолемею это терм Венеры. Следовательно функция выдаст 3 (номер Венеры в системе ШЭ).

TermName() – Делает то же самое, что и функция Term(), только вместо номера возвращает имя планеты.

Пример. Нас интересует 8 градус Овна. По Птолемею это терм Венеры. Следовательно функция выдаст “Venus”

Midp() – Возвращает среднюю точку между А и Б

MidpHar() - Возвращает среднюю точку между А и Б в заданной гармонике. В градации 0-360

MidpHarN() - Возвращает среднюю точку между A и B в заданной гармонике в градации 0-360, деленной на номер гармоники

Пример, для 4 гармоники (90 круг), максимально возможная координата = 90 градусов). Для того, чтобы проверить картину вида: A/B=C/D, просто вычислите: MidpHarN(A,B) и MidpHarN(C,D), если разница меньше 1 градуса, значит картина есть (орб. 1 градус).

AandBsubC() - Вычисляет точку вида A+(B-C) в круге 360. Для перевода результата, на пример, в 90 круг, воспользуйтесь функцией Harmonica, а для того, чтобы перевести то, что получилось в градацию 0-90 и т.д., разделите результат на номер гармоники

AandBdivC() - Вычисляет точку вида A+B/C в круге 360.

AandB() - Вычисляет точку вида A+B.

Grebiy() - Вычисляет жребий (формула вида A=B-C для дневной карты и A=C-B для ночной карты). Для вычисления жребия без учтета дневной/ночной карты просто воспользуйтесь функцией AandBsubC.


Более подробное описание новых функций см. в документации, если так можно выразица


*****
Люди, документация версталась в спешке, очень хотелось успеть закончить это дело до нового года, поэтому за верстку не судите строго )) Обещаю на праздниках подготовить нормальную документацию, и структурировать её в более удобочитаемом виде. :) Очень большая просьба в случае глюков, багов, косяков и пр. и пр. отписаться на форуме, либо lordwilex@lordwilex.ru , а я побежал готовиться к НГ :)
******

Скачать последнюю версию Швейцарских эфемерид по-русски можно тут: http://sweru.lordwilex.ru/sweru.zip

Документация тут: http://sweru.lordwilex.ru/

******

А я всем желаю веселого новогодопрепровожждения, да прибудет с нами сила!!!

[img]http://s16.**********/f59f71b9e322cb609b9823bf6a81f4a4.gif[/img][img]http://s18.**********/291af55c62cdb3261bb82c7ed2c19432.gif[/img] [img]http://s16.**********/f59f71b9e322cb609b9823bf6a81f4a4.gif[/img]

LordWilex вне форума   Ответить с цитированием
Старый 29.12.2009, 08:22   #152
Алексей
Собеседник
 
Аватар для Алексей
 
Регистрация: 01.11.2007
Сообщения: 8,521
Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000Алексей репутация выше +2000
По умолчанию

Молодчина, Виталий!!!

Алексей вне форума   Ответить с цитированием
Старый 29.12.2009, 08:28   #153
LordWilex
В отпуске
 
Аватар для LordWilex
 
Регистрация: 01.06.2008
Адрес: Таганрог
Сообщения: 28,983
LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000
По умолчанию

Цитата:
Сообщение от Алексей
Молодчина, Виталий!!!

Ну там еще наверняка полно ошибок, после НГ придеться исправлять, как обычно все делалось в попыхах
LordWilex вне форума   Ответить с цитированием
Старый 30.12.2009, 09:05   #154
vitek_tlt
Собеседник
 
Регистрация: 29.12.2009
Сообщения: 10
vitek_tlt отключил(а) отображение уровня репутации
По умолчанию

Цитата:
Сообщение от LordWilex
Ну там еще наверняка полно ошибок, после НГ придеться исправлять, как обычно все делалось в попыхах


Спасибо!!!Титанический труд.....
vitek_tlt вне форума   Ответить с цитированием
Старый 30.12.2009, 09:25   #155
LordWilex
В отпуске
 
Аватар для LordWilex
 
Регистрация: 01.06.2008
Адрес: Таганрог
Сообщения: 28,983
LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000
По умолчанию

Цитата:
Сообщение от vitek_tlt
Спасибо!!!Титанический труд.....
Та ни, там все легко, на несколько ночей работы, если бы не лень, давно бы сделал из этого проекта астрологическогоо монстра
LordWilex вне форума   Ответить с цитированием
Старый 05.01.2010, 18:27   #156
LordWilex
В отпуске
 
Аватар для LordWilex
 
Регистрация: 01.06.2008
Адрес: Таганрог
Сообщения: 28,983
LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000
По умолчанию

Небольшой ап.


****
По просбам трудящихся добавленна функция для работы с дирекциями.

Вы можите скопировать от сюда:

Цитата:
Public Function DirectArc(ByVal BirthJD As Double, ByVal SobJD As Double, ByVal PCoeff As Double) As Double
'
' Выводит юлианскую дату соответствующую дате события, соответствующему положению дуги прогрессивного Солнца
' Где
' BirthJD - Дата рождения
' SobJD - Дата события
' PCoeff - Прогрессивный коэффициент (продолжительность года 365.25, 365.2425 или 265.2422 суток)
'
'
'

Arc = (SobJD - BirthJD) / PCoeff
days = 360 / PCoeff
DirectArc = BirthJD + (Arc * days)

End Function


Или скачать модуль целиком со всеми причиндалами, от предыдущей версии отличие только в наличии данной функции.


Как ей пользоваться?

Допустим, вам нужно расчитать дирекцию Солнечной дуги для человека, который родился 10-08-1980 17:00:35 по Московскому времени на дату 31:12:2009 12:09:07 По московскому времени.

1. Вычисляем юлианскую дату рождения человека, при помощи функции jdayLT():
10-08-1980 17:00:35 = 2444462,084

2. Вычисляем юлианскую дату события, при помощи функции jdayLT():
31:12:2009 12:09:07 = 2455196,881

3. В качестве дирекционного коэффицента берем 365.2422 (величину солнечного года в днях, хотя это завист от типа дирекции, можно, на пример, взять и 365.25).

4. Вводим данные в функцию
DirectArc(2444462.084, 2455196.881 , 365.2422)
= 2444491.051


В результате получаем дату: 2444491.051 , которую подставляем в функцию PLC() для расчета координат планет.

5. Чтобы перевести эту дату в обычную воспользуемся функцией JDToDay(), в результате получим: 8-9-1980 13:13:45 .


****
Описание новой функции в документацию пока не включено, так что считайте это неофициаальным релизом, документацией займусь по позже
LordWilex вне форума   Ответить с цитированием
Старый 14.01.2010, 17:42   #157
vitek_tlt
Собеседник
 
Регистрация: 29.12.2009
Сообщения: 10
vitek_tlt отключил(а) отображение уровня репутации
По умолчанию

Я посмотрел в архиве есть под Делфи файлы и они вполне работоспособны. Зачем заморачиваться с офисом, не понял теперь. Все прекрасно работает. Визуально.


Последний раз редактировалось vitek_tlt, 14.01.2010 в 17:47.
vitek_tlt вне форума   Ответить с цитированием
Старый 14.01.2010, 18:44   #158
LordWilex
В отпуске
 
Аватар для LordWilex
 
Регистрация: 01.06.2008
Адрес: Таганрог
Сообщения: 28,983
LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000
По умолчанию

Цитата:
Сообщение от vitek_tlt
Я посмотрел в архиве есть под Делфи файлы и они вполне работоспособны. Зачем заморачиваться с офисом, не понял теперь. Все прекрасно работает. Визуально.


Вообще-то этот пакет предназначен в первую очередь для целей статистической обработки больших масивов информации простых генераторов эфемерид в текстовом виде в интернете и так полно, почти любая астропрограмма это умеет делать, тока как вы при помощи этой проги собираетесь работать со статистикой?
Если это вас устраивает, то пользуйтесь на здоровье
LordWilex вне форума   Ответить с цитированием
Старый 17.01.2010, 15:44   #159
Edik
Собеседник
 
Регистрация: 19.12.2008
Сообщения: 7
Edik репутация выше +10
По умолчанию

Здраствуйте,

Тут такой вопросик выскочил: как можно посчитать координаты Прозерпины и Селены (Белой Луны) - есть такие планеты по Глобе :)

Эдуард
Edik вне форума   Ответить с цитированием
Старый 17.01.2010, 19:43   #160
LordWilex
В отпуске
 
Аватар для LordWilex
 
Регистрация: 01.06.2008
Адрес: Таганрог
Сообщения: 28,983
LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000LordWilex репутация выше +1000
По умолчанию

Цитата:
Сообщение от Edik
Здраствуйте,

Тут такой вопросик выскочил: как можно посчитать координаты Прозерпины и Селены (Белой Луны) - есть такие планеты по Глобе :)

Эдуард

Это есть в документации и в этой теме.


Вот процедура подключения подобного рода объектов на примере Прозерпины:

Цитата:
Работа с нестандартными фиктивными планетами в ШЭ?



Что делать, если вашей любимой фиктивной планеты нет в Швейцарских Эфемеридах?



Две Прозерпины….

Вчера один из посетителей форума обратился ко мне с вопросом, «а куда вы, собственно говоря, дели мою любимую Прозерпину?». :)

Я с сей экзотичной планиДой ни когда не работал, но будучи уверен, что в ШЭ все есть, сначала посоветовал почитать внимательней форум, но когда внимательней поискал, все оказалось не так просто, как казалось, я обнаружил, что Прозерпины там действительно нет…

Верней я её конечно нашел, но эта Прозерпина оказалась какой-то не такой, во первых, её элементы оказались геоцентрическими, в то время, как мой знакомый уверял меня, что они должны быть геолиоцентрическими, во вторых координаты той Прозерпины, которая была в ШЭ (Ско 016° 22' 02,6953''), градусов на 5 отличались от того, что выдавала программа моего знакомого (Ско 011° 31' 50,3279'').



В общем-то если у вас имеются орбитальные элементы нужной вам фиктивной планеты, - это не проблема!

Достаточно создать файл seorbel.txt, добавить в него орбитальные элементы нужной «планеты» и поместить этот файл в папку со швейцарскими эфемеридами (или дописать элементы нужной фиктивной планеты в файл, если такой уже существует) и.. и все…



На пример, для этой самой Прозерпины, строка с орбитальными элементами будет выглядеть следующим образом (орбитальные элементы Прозерпины позаимствованы мной из программы ZET):



2445375.5, 2445375.5, 187.502320, 77.926077235, 0.08550236, 319.538004, 63.696644, 3.000799, Прозерпина # 19 (PLC = 59)



Где каждый элемент отделяется друг от друга запятой:



1. Эпоха относительно которой вычисляем координаты (Юлианская дата)

2. Эквоникс (Юлианская дата, либо один из этих: "J1900" или "B1950" или "J2000")

3. Средняя аномалия

4. Большая полуось орбиты

5. Экцентриситет

6. Аргумент перегелия

7. Долгота восходящего узла

8. Наклон орбиты к эклиптики

9. Название планеты

10. Если элементы геоцентрические, после названия, через запятую пишем еще и geo



Знак - '#' обозначает комментарий, после этого знака можите писать что угодно, программа это все равно не увидет :)



Для вычисления координат при помощи функции swe_calc() (в эксели функция PLC() ), используйте номер планеты

Который вычисляется по формуле:

ipl = SE_FICT_OFFSET_1 + номер элементов в этом файле. На пример, номер Крона ipl = 39 + 4 = 43



А вот и пример файла с орбитальными элементами (перевод с английского мой, прошу больно не пинать, если что не так перевелось :) ):




Цитата:
Сообщение от Орбитальные элементы фиктивных планет
# Орбитальные элементы фиктивных планет
# 27 Jan. 2000

#

# Этот файл является частью Швейцарских Эфемерид, версии 1.52.

# Пользователь может добавить свои собственные фиктивные планеты (максимум 960 планет).

#

# Порядок элементов в строке следующий:

# 1. Эпоха относительно которой вычисляем координаты (Юлианская дата)

# 2. Эквоникс (Юлианская дата, либо один из этих: "J1900" или "B1950" или "J2000")

# 3. Средняя аномалия

# 4. Большая полуось орбиты

# 5. Экцентриситет

# 6. Аргумент перегелия

# 7. Долгота восходящего узла

# 8. Наклон орбиты к эклиптики

# 9. Название планеты

# 10. Если элементы геоцентрические, после названия, через запятую пишем еще и geo

#

# Знак - '#' обозначает комментарий, после этого знака можите писать что угодно, программа это все равно не увидет :)

# Для вычисления координат при помощи функции swe_calc() (в эксели функция PLC() ), используйте номер планеты

# Который вычисляется по формуле:

# ipl = SE_FICT_OFFSET_1 + номер элементов в этом файле. На пример, номер Крона ipl = 39 + 4 = 43

#

# Транснептуны ГША, по данным Джейма Нили (James Neely)

J1900, J1900, 163.7409, 40.99837, 0.00460, 171.4333, 129.8325, 1.0833, Купидон # 1 (PLC = 40)

J1900, J1900, 27.6496, 50.66744, 0.00245, 148.1796, 161.3339, 1.0500, Гадес # 2 (PLC = 41)

J1900, J1900, 165.1232, 59.21436, 0.00120, 299.0440, 0.0000, 0.0000, Зевс # 3 (PLC = 42)

J1900, J1900, 169.0193, 64.81960, 0.00305, 208.8801, 0.0000, 0.0000, Крон # 4 (PLC = 43)

J1900, J1900, 138.0533, 70.29949, 0.00000, 0.0000, 0.0000, 0.0000, Апплалон # 5 (PLC = 44)

J1900, J1900, 351.3350, 73.62765, 0.00000, 0.0000, 0.0000, 0.0000, Адмет # 6 (PLC = 45)

J1900, J1900, 55.8983, 77.25568, 0.00000, 0.0000, 0.0000, 0.0000, Вулкан # 7 (PLC = 46)

J1900, J1900, 165.5163, 83.66907, 0.00000, 0.0000, 0.0000, 0.0000, Посейдон # 8 (PLC = 47)

#

# Изида трансплутоновая "Die Sterne" 3/1952, p. 70ff.

2368547.66, 2431456.5, 0.0, 77.775, 0.3, 0.7, 0, 0, Isis-Transpluto # 9 (PLC = 48)

# Нибиру элементы по Криятьяну Уолтегу (Christian Woeltge)

1856113.380954, 1856113.380954, 0.0, 234.8921, 0.981092, 103.966, -44.567, 158.708, Нибиру # 10 (PLC = 49)

# 10-я планета по Харингтону (Harrington), по данным Astronomical Journal 96(4), Oct. 1988

2374696.5, J2000, 0.0, 101.2, 0.411, 208.5, 275.4, 32.4, Harrington # 11 (PLC = 50)

# Другие попытки найти 10-ю "планету - X" по данным W.G. Hoyt, "Planets X and Pluto", Tucson 1980, p. 63

2395662.5, 2395662.5, 34.05, 36.15, 0.10761, 284.75, 0, 0, Leverrier (Neptune) # 12 (PLC = 51)

2395662.5, 2395662.5, 24.28, 37.25, 0.12062, 299.11, 0, 0, Adams (Neptune) # 13 (PLC = 52)

2425977.5, 2425977.5, 281, 43.0, 0.202, 204.9, 0, 0, Lowell (Pluto) # 14 (PLC = 53)

2425977.5, 2425977.5, 48.95, 55.1, 0.31, 280.1, 100, 15, Pickering (Pluto) # 15 (PLC = 54)

# Вулкан (интромеркуриарный :) ) по данным Л. Уестона (L.H. Weston)

J1900,JDATE, 252.8987988 + 707550.7341 * T, 0.13744, 0.019, 322.212069+1670.056*T, 47.787931-1670.056*T, 7.5, Vulcan # 16 (PLC = 55)

# Селена / Белая луна

J2000,JDATE, 242.2205555 + 5143.5418158 * T, 0.05280098949, 0.0, 0.0, 0.0, 0.0, Selena/White Moon, geo # 17 (PLC = 56)

# Прозерпина, по данным http://www.geocities.com/Hollywood/A...roserpina.html

# J1900, 170.73 + 51.05 * T

J1900,JDATE, 170.73, 79.225630, 0, 0, 0, 0, Proserpina (гео) #18 (PLC = 57)

# Вторая Луна по Уальдемату (Waldemath's Second Earth Moon)

2414290.95827875,2414290.95827875, 70.3407215 + 109023.2634989 * T, 0.0068400705250028, 0.1587, 8.14049594 + 2393.47417444 * T, 136.24878256 - 1131.71719709 * T, 2.5, Waldemath, geo # 19 (PLC = 58)

######################################## ######

### Дополнительные элементы

# Прозерпина по элементам, позаимствованным из программы ZET

2445375.5, 2445375.5, 187.502320, 77.926077235, 0.08550236, 319.538004, 63.696644, 3.000799, Прозерпина # 19 (PLC = 59)

# Дхарма по элементам, позаимствованным из программы ZET

2448830.500, 2448830.500, 050.7800, 57.839761679, 0.11800000, 288.017908, 047.587053, 2.500595, Дхарма # 20 (PLC = 60)


Что касается Селены, то процедура та же, её орбитальные элементы есть, на пример, в ЗЕТ.

Последний раз редактировалось LordWilex, 17.01.2010 в 19:50.
LordWilex вне форума   Ответить с цитированием
Ответ


Опции темы
Опции просмотра

Ваши права в разделе
Вы не можете создавать темы
Вы не можете отвечать на сообщения
Вы не можете прикреплять файлы
Вы не можете редактировать сообщения

BB коды Вкл.
Смайлы Вкл.
[IMG] код Вкл.
HTML код Выкл.
Быстрый переход


Часовой пояс GMT +1, время: 21:09.


Powered by vBulletin Version 3.5.4
Copyright ©2000 - 2024, Jelsoft Enterprises Ltd.
© 1995-2024, ARGO