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

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

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

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

Непрочитанное сообщение vasnas » 07 июн 2012 15:43

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

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

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

Ссылка на документацию SPICE NASA в т.ч. на русском http://www.bible-exodus.narod2.ru/artic ... meris.html
Вложения
}precession_earth_eclipj2000.png
Последний раз редактировалось vasnas 07 июн 2012 15:46, всего редактировалось 1 раз.
vasnas
Статус: Новичок
Статус: Новичок
 
Сообщения: 35
Зарегистрирован: 10 апр 2012 13:55
Благодарил (а): 0 раз.
Поблагодарили: 0 раз.

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

Непрочитанное сообщение vasnas » 07 июн 2012 15:44

Мой скрипт:

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

% 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
vasnas
Статус: Новичок
Статус: Новичок
 
Сообщения: 35
Зарегистрирован: 10 апр 2012 13:55
Благодарил (а): 0 раз.
Поблагодарили: 0 раз.

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

Непрочитанное сообщение sergey2011 » 09 июн 2012 12:15

Так много программ! Не скрывайте от нас конец графика! Это математическое доказательство, что мы все умрем???
sergey2011
Статус: Новичок
Статус: Новичок
 
Сообщения: 26
Зарегистрирован: 02 сен 2011 19:54
Благодарил (а): 0 раз.
Поблагодарили: 0 раз.


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

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

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

cron