Re: Вопросы Markab-у


Автор сообщения: zakharov
Дата и время сообщения: 28 June 2006 at 20:23:06:

В ответ на сообщение: Re: Вопросы Markab-у

Уважаемый Markab!

Пересчет координат по разным формулам на весьма удаленные эпохи может давть отличия до десятков секунд. Это связано:

1. С учетом собственных движений.
а) Звезды движутся в пространстве, поэтому надо учитывать не только тангециальные движения, но идвижение по лучу зрения (т.е. учитывать расстояние - параллакс и лучевые скорости). За счет этого, кстати, у близких звезд быстро приближающихся или удаляющихся может немного меняться блеск.
б) Звезды могут двигаться относительно Солнца не по прямой. Например они могут составлять пары с периодами тысячи лет; и, если звезды достаточно близки к нам, это может искажать собственные движения. Поэтому чем больше расстояние между эпохами при вычислении собственного движения, тем точнее оно определяется.
г) Только для некоторых (около двух десятков) звезд даются орбиты и собственные движения центра тяжести (например, для alp Cen). Для остальных в собственное движение входят и орбиты.

2. С учетом прецессии.
а) Современные формулы для учета прецессии содержат высокие гармоники, поэтому они ограничено годны во времени (например, формулы прецессии, использованные в большинстве программ, рекомендованы МАС (IAU) к использованию начиная с 1980 года, а с 2003 года вводятся другие).
б) Поэтому на 1-2 тысячи лет лучше использовать формулы без учета высших членов (до кубического члена вполне достаточно). А на десятки тысяч лет и этого много.

3. С учетом нутации (т.е. с приведением на видимое место). Нутацию можно не учитывать
а) Потому, что размах изменений за ее счет не превышает 18.5 угл. секунд.
б) Потому, что мы не знаем точной эпохи каталога Птолемея, а период нутации равен 18.6 лет

4. С учетом рефракции.
а) Я бы ее не учитывал, потому, что для большинства звезд мы не знаем высоты наблюдения над горизонтом.
б) Возможно ее надо учитывать для южных звезд, хотя от 5 градусов над горизонтом рефракция не превышает 10 угловых минут. Для звезд на истинном зенитном расстоянии 90 градусов рефракция примерно 1/2 градуса. При этом ошибка долготе может составить более градуса.

Программа пересчета координат, которой пользовался я (http://www.newchrono.net/almagest/zakharov.pdf) производит учет пpeцeccии в cooтвeтcтвии c пocтoяннoй пpeцeccии MAC 1976 гoдa. Она точна лишь при использовании координат стандартной эпохи и равноденствия J2000,0.

Если Вы ею воспользуетесь, то хочу предупредить об очепятке. Строку
epsilon = (84381.47 - 46.8156*dt - 0.00041*dt2 + 0.00183*dt2)/sectorad;
следует читать
epsilon = (84381.47 - 46.8156*dt - 0.00041*dt2 + 0.00183*dt3)/sectorad;
Я не понимаю откуда она взялась. В моих версиях все верно.

Если Вы захотите использовать Нькомбовские постоянные, то исправте вычисление углов прецессии


HY=36524.2198781d0; /* продолжительность бесселева столетия точно*/
T0=2415020.31352d0 /* эпоха B1900.0 */

T = (JDE1 - T0)/HY /* JDE1 - эпоха тех координат, которые Вы имеете,
например, если вы взяли координаты эпохи 1950.0 года,
то: JDE1=(1950-1900)*365.242198781+2415020.31352.
Аналогично JDE2 - эпоха, на которую рассчитываются
координаты */
dT = (JDE2 - JDE1)/HY;
T2 = T*T;
dT2 = dT*dT;
dT3 = dT2*dT;
/* углы пpецессии в секундах дуги */
ss2 = (2304.253 + 1.3973*T + 0.00006*T2)*dT;
ss1 = ss2 + (0.3023 - 0.0027*T)*dT2 + 0.01800*dT3;
ss3 = ss2 + (1.0950 + 0.0039*T)*dT2 + 0.01832*dT3;
ss2 = (2004.685 - 0.8533*T - 0.00037*T2)*dT -(0.4267 + 0.00037*T)*dT2 - 0.0418*dT3;


С уважением
Захаров Андрей



2340. Новохронологический анонс :)))) - А.Калашников 18:04 13.06.06 (245)
К списку тем на странице