Вход

Просмотр полной версии : FK4 vs ECLIPJ2000 где прецессия земной оси?


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

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

Получается, что системы координат FK4 и GALACTIC в программном инструменте SPICE привязаны к прецесии ECLIPJ2000 и не соответствуют своему названию?

http://i46.fastpic.ru/thumb/2013/0409/2d/650cb9f1548181a634e3658d0c5a472d.jpeg (http://fastpic.ru/view/46/2013/0409/650cb9f1548181a634e3658d0c5a472d.png.html)

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

% долгота и широта Солнца в 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

астровадян
09.04.2013, 22:07
Ещё один такой бессмысленный пост - отпралю на месяц отдохнуть от форума. Уже не в первой делаю замечания о недопустимости бессмысленный выкриков! :mad:
Грин.

vasnas
11.04.2013, 11:23
Я задал теоретический вопрос, т..к. системы координат FK4 или ECLIPJ2000 международные и относятся к любой программе ...

Объясните Ваши слова, чтобя я мог соответсвующе "себя вести".

Нельзя задавать астрономические вопросы по теории... Или что?

Грин
11.04.2013, 15:47
Моё замечание было не Вам, а пользователю астровадян.

астровадян
11.04.2013, 22:14
Ещё один такой бессмысленный пост - отпралю на месяц отдохнуть от форума. Уже не в первой делаю замечания о недопустимости бессмысленный выкриков! :mad:
Грин.ладно согласен люблю честно говоря пофлудить, раскаиваюсь:cry:
если по теме лично для меня это слишком нагруженная не понятным теоретическим материалом статья с момента про скрипты к мат лаб, это я собственно и хотел сказать.

vasnas
11.04.2013, 22:29
можно отбросить загруженность и ограничиться

В моменты весенних равноденствий в FK4 и GALACTIC долгота солнца не изменяется, что говорит о том, что эти системы координат не соответствуют своим номиналам и привязаны к ECLIPJ2000 системе координат в плоскости эклиптики.

FK4- фундаментальный каталог (см. Интернет)