|
Получив моменты астрономических весенних равноденствий, когда долгота солнца в системе координат 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
|