|
Получив моменты астрономических весенних равноденствий, когда долгота солнца в системе координат ECLIPJ2000 равна 0, я предполагаю, что в системах координат FK4 и GALACTIC долгота солнца на эти моменты будет изменяться в положительную сторону из-за прецессии, которая равна 1-му обороту точки равноденствия в эклиптике за ~ 25 765 лет
, но этого не наблюдается (см. график)!
Получается, что системы координат FK4 и GALACTIC в программном инструменте SPICE NASA (http://bible-exodus.narod2.ru/articles/astro_ephemeris/jpl_ephemeris/jpl_ephemeris.html )
привязаны к прецессии ECLIPJ2000 и не соответствуют своему названию?
http://fastpic.ru/view/46/2013/0409/650cb9f1548181a634e3658d0c5a472d.png.html>![]() 
============================================
Прилагаю скрипт 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;
= cspice_spkpos ( 'SUN', st_win(ii), ' FK4 ', abcorr, observer); % позиция Солнца в момент весеннего равноденствия по FK4
= 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
|
.
