Где прецессия Земли в FK4 и GALACTIC?

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

Модераторы: Ulmo, Булдаков Сергей

Где прецессия Земли в FK4 и GALACTIC?

Непрочитанное сообщение vasnas » 09 апр 2013 19:35

FK4 vs ECLIPJ2000 где прецессия земной оси?
Получив моменты астрономических весенних равноденствий, когда долгота солнца в системе координат ECLIPJ2000 равна 0, я предполагаю, что в системах координат FK4 и GALACTIC долгота солнца на эти моменты будет изменяться в положительную сторону из-за прецессии, которая равна 1-му обороту точки равноденствия в эклиптике за ~ 25 765 лет

, но этого не наблюдается (см. график)!

Получается, что системы координат FK4 и GALACTIC в программном инструменте SPICE NASA (http://bible-exodus.narod2.ru/articles/ ... meris.html )

привязаны к прецессии ECLIPJ2000 и не соответствуют своему названию?

Изображение

============================================
Прилагаю скрипт MatLab, которым получал график за 6000 лет.

% долгота и широта Солнца в FK4, во время весенних равноденствий ECLIPJ2000
% LON and LAT of sun in FK4 at equinoxes
cspice_kclear; clear all
cspice_furnsh( 'd:\Documents and Settings\user\Мои документы\MATLAB\astro\mice\data\standard.tm' );
TFMT='DD.MM.YYYY HR:MN:SC ERA::RND::MCAL::UTC'; TFMTH='HR:MN:SC.## ::RND::UTC'; a_lon=zeros(0); a_lat=zeros(0);

target = 'Sun'; % объект наблюдения солнце
abcorr = 'LT+S'; % поправка на, видимое наблюдателем, положение солнца в небе

st1= cspice_str2et ( '01-01-2999 BC' ); % начало периода
st2 = cspice_str2et( '01-01-3000 AD' ); % конец периода
st_win = cspice_wninsd ( st1, st2 ); % период



% загрузка координат наблюдателя в систему координат
observer='Cairo' ; % место наблюдения Каир (Египет) CAIRO
crds = 'LATITUDINAL'; % crds = 'CYLINDRICAL'; %
crd = 'LONGITUDE' ; % ищется эклиптическая долгота

adj = 0.0; % допуск
MAXWIN = 40000; MAXIVL = MAXWIN / 2;
step = cspice_spd * 180;

% получаем моменты астрономических равноденствий - geting of moments of equinoxes for the given period
st_win = cspice_gfposc ( 'SUN', 'ECLIPJ2000', abcorr, observer, crds, crd,'=', 0 * cspice_rpd , adj, step, MAXIVL, st_win);
st_win=(reshape(st_win,2,[]))';
st_win = st_win(:,1);
st_win= st_win (:,1); % т.к. начала и конец периодов моментов равны, то отбросим пары - get off the ends or the starts, because they are =

%fprintf ( '\n\n\n The LON of the sun in equinoxes moments in FK4 - Долгота солнца в весенние равноденствия по FK4 : \n' );
n=numel(st_win);
for ii=1:n;
[rectan, lt] = cspice_spkpos ( 'SUN', st_win(ii), ' FK4 ', abcorr, observer); % позиция Солнца в момент весеннего равноденствия по FK4
[r, lat, lon] = cspice_reclat(rectan); % преобразовываем позицию Солнца из прямоугольной системы координат в прямое восхождение и склонение
lon =cspice_convrt(lon,'radians','degrees'); a_lon(ii)=lon;
lat=cspice_convrt(lat,'radians','degrees'); a_lat(ii)=lat;
%disp(strcat( sprintf( 'Долгота, ° : % 12.25f\n', lon ), ' (', angl2minsec(lon), ')'));
%disp(strcat( sprintf( 'Широт а, ° : % 12.25f\n', lat ), ' (', angl2minsec(lat), ')'));
end

p=figure('Name','Долгота и высота Солнца в FK4 в весенние равноденствия ','NumberTitle','off');
subplot(2,1,1); p=plot(a_lon);
set(p,'Color','red','LineWidth',1)
title('Lon ° of Sun in FK4 frame at ECLIPJ2000''s vernal equinoxes')
subplot(2,1,2); p=plot(a_lat);
title('Lat ° of Sun in FK4 frame at ECLIPJ2000''s vernal equinoxes')
set(p,'Color','green','LineWidth',1)

cspice_kclear
Последний раз редактировалось vasnas 11 апр 2013 10:44, всего редактировалось 3 раз(а).
vasnas
Статус: Новичок
Статус: Новичок
 
Сообщения: 35
Зарегистрирован: 10 апр 2012 13:55
Благодарил (а): 0 раз.
Поблагодарили: 0 раз.

Re: Где нутация Земли в FK4 и GALACTIC?

Непрочитанное сообщение vasnas » 11 апр 2013 10:27

никто не поможет?
vasnas
Статус: Новичок
Статус: Новичок
 
Сообщения: 35
Зарегистрирован: 10 апр 2012 13:55
Благодарил (а): 0 раз.
Поблагодарили: 0 раз.


Вернуться в Общение астрономов и любителей

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 7