Почему дрейфует широта Солнца в эклиптике в момент В.Р.?

[vasnas]

Почему дрейфует широта Солнце в эклиптике в момент В.Р.?

В момент весеннего равноденствия "я из центра Земли" измеряю долготу и широту Солнца в эклиптике.

Пользуюсь программным инструментарием SPICE NASA в Matlab.

Долгота, как и следовало ожидать ~0, из определения системы координат ECLIPJ2000, что она всегда привязана к В.Р.
Но почему дрейфует широта? Ведь широта солнца должна быть в плоскости эклиптики, т.е. ~0

Ссылка на документацию SPICE NASA в т.ч. на русском http://www.bible-exodus.narod2.ru/artic ... meris.html


Мой скрипт:

============

% Precession of Earth
clear all;
cspice_furnsh( 'd:\Documents and Settings\user\Мои документы\MATLAB\astro\mice\data\standard.tm' );
TIMFMT = 'DD.MM.YYYY HR:MN ::RND::MCAL::UTC';

et = cspice_str2et ( '1966-03-15' ); % Start at ...
years = 68; data= zeros( years, 4); % through years



for ii=1:years;
% approximately the date of the vernal equinox
st_win = cspice_wninsd( et - cspice_spd * 10, et + cspice_spd * 10);

% exactly an ET of the vernal equinox
abcorr = 'LT'; crds = 'LATITUDINAL' ; crd = 'LONGITUDE'; adj = 0.0; MAXWIN = 20000; MAXIVL=MAXWIN/2; step= cspice_spd * 9;
res = cspice_gfposc ( 'Sun', 'ECLIPJ2000', abcorr, '399', crds, crd,'=', 0 , adj, step, MAXIVL, st_win);
et_equinox= res(1);

% position of Sun concerning stars at the vernal equinox
[pos, ~] = cspice_spkpos ( 'Sun', et_equinox , 'ECLIPJ2000', abcorr, '399');
[~, lon, lat] = cspice_reclat(pos); lon = lon * cspice_dpr; lat = lat * cspice_dpr;
data(ii,1) = str2double( cspice_timout( et, 'YYYY')); data(ii,2) = et;
data(ii,3) = lon; data(ii,4) = lat;

fprintf('% s\n', cspice_timout( et, TIMFMT));
et = et + cspice_tyear;
end


figure1 = figure('Name','Precession of Earth');
title({'Precession of Earth',''}); xlabel('Year');
[AX,H1,H2] = plotyy(data(:,1), data(:,3), data(:,1), data(:,4),'plot');
set(get(AX(1),'Ylabel'),'String','lon')
set(get(AX(2),'Ylabel'),'String','lat')


cspice_kclear