|
>С учётом разброса протяжённости созвездий избранный метод
>является ошибочным. Он бездумно перенесён из анализа
>знакового гороскопа (где все знаки имеют одинаковую
>протяжённость 30%). То есть, у вас ошибки на разных
>созвездиях имеют разный вес, и, по сути, вы повторяете
>нелепицу Лжеархимеда, соврершенно незнакомого с математикой
>и статистикой.
>
>Таким образом, вам следует применить свой инструмент почище.
>Интересно - что он даст в этом случае?
всё это учтено!
, просто с переездом сайтов на UCOZ я не могу опубликовать полностью решение...
пока вот здесь без картинок http://bible-exodus.narod.ru/articles/astro_ephemeris/jpl_ephemeris/jpl_ephemeris.html#zodiac_sibli_ae
и вот программа (проверяйте):
% вычисление даты зодиака Христа из книги Эбенезера Сибли (зодиак "AE")
% calculating of date of ...
% разбивка, принятая ФиН, с 17.11.2007 г. см. http://www.chronologia.org/egruit_zodiaki/pril.html
% Fomenko and Nosovsky's separation since 11/17/2007
% РАЗБИВКА ЗОДИАКА% In ecliptics % In constellations
% СОЗВЕЗДИЕ/Const ОТ ЭКЛИПТИКИ В СОЗВЕЗДИЯХ
% 1 Овен 0-30 31-56
% 2 Телец 30-60 56-92
% 3 Близнецы 60-90 92-118
% 4 Рак 90-120 118-137
% 5 Лев 120-150 137-172
% 6 Дева 150-180 172-215
% 7 Весы 180-210 215-236
% 8 Скорпион 210-240 236-266
% 9 Стрелец 240-270 266-296
% 10 Козерог 270-300 296-326
% 11 Водолей 300-330 326-349
% 12 Рыбы 330-360 349-31
% во всех условиях интервал расширяется в обеи стороны на 5%
% in all conditions the interval wided onto 5 %
clear all;
cspice_furnsh( 'd:\Documents and Settings\user\Мои документы\MATLAB\astro\mice\data\standard.tm' );
TIMFMT = 'DD.MM.YYYY ERA HR:MN ::RND::JCAL::UTC'; % JCAL - Юлианский календарь, Рождество было до Григореанской реформы
TIMLEN = 35; % выводимый формат времени
d=0; % поправка по эклиптике для FK4
res_date = zeros(1);
% начальный интервал времени
et1=cspice_str2et ( '12-21-1000 BC' ); et2=cspice_str2et ( '01-01-1700 AD' ); res = cspice_wninsd (et1 ,et2 );
obsrvr='CAIRO' ; % место наблюдения Каир (Египет)
<~, lon, ~, ~, ~> = get_observer (obsrvr); % узнать долготу места для местного солн.врем. % Lon of pl.
crds = 'CYLINDRICAL' ; crd = 'LONGITUDE' ; % ищется эклиптическая долгота
adj = 0.0; MAXWIN = 20000; MAXIVL = MAXWIN / 2;
abcr = 'LT+S'; % учитывать оберацию и скорость света, хотя это излишне
Tol=5; % +- допуск к границам созвездий
% для ускорения расчетов ограничиваем период окнами времени, начиная с медленных планет
% если диапазон долготы планеты пересекает ноль то вводим две записи,
% с параметрами до нуля и после
% for increasing of speed of calulations arrange "slowest" planets to the top of data array
% if interval of position cross 0 then put 2 records (duplet) from begining till 360 degree and from 0 to the end
% % strong conditions +- 5°
% params = {
% 'JUPITER BARYCENTER', (215-Tol), (239+Tol), 30; % Юпитер в Весах |J. in Libra
% 'SATURN BARYCENTER', ( 92-Tol), (118+Tol), 30; % Сатурн в Близнецах | S. in Gemini
% 'VENUS', (296-Tol), (326+Tol), 5; % Венера в Козероге | Venus in Capricorn
%
% 'MARS', (349-Tol), (359.99999999), 2; % Марс в Козероге | Mars in Capricorn
% 'MARS', (0), (31+Tol), 2; % Марс в Козероге | Mars in Capricorn
%
% 'MERCURY', (266-Tol), (296+Tol), 1; % Меркурий в Стрельце | Mercury in Sagittarius
% 'SUN', (266-Tol), (296+Tol), 1/3; % Солнце в Стрельце | Sun in Sagittarius
%
% 'MOON', ( 56-Tol), ( 92+Tol), 1/12; % Луна в Тельце | Moon in Taurus
% };
% % soft conditions +-5° and +- 3/4 ecliptics longitude of constellations
% params = {
% 'JUPITER BARYCENTER', (182.75-Tol), (259.25+Tol), 30; % Юпитер в Весах |J. in Libra
% 'SATURN BARYCENTER', ( 65-Tol), (132.25+Tol), 30; % Сатурн в Близнецах | S. in Gemini
% 'VENUS', (273.5-Tol), (343.25+Tol), 5; % Венера в Козероге | Venus in Capricorn
%
% 'MARS', (331.75-Tol), (359.99999999), 2; % Марс в Козероге | Mars in Capricorn
% 'MARS', (0), ( 49.75+Tol), 2; % Марс в Козероге | Mars in Capricorn
%
% 'MERCURY', (245.75-Tol), (318.5+Tol), 1; % Меркурий в Стрельце | Mercury in Sagittarius
% 'SUN', (245.75-Tol), (318.5+Tol), 1/3; % Солнце в Стрельце | Sun in Sagittarius
%
% 'MOON', ( 37.25-Tol), (111.5+Tol), 1/12; % Луна в Тельце | Moon in Taurus
% };
%
% моя расстановка по созвездиям
params = {
'JUPITER BARYCENTER', (215-Tol), (236+Tol), 30; % Юпитер в Весах |J. in Libra
'SATURN BARYCENTER', ( 92-Tol), (118+Tol), 30; % Сатурн в Близнецах | S. in Gemini
'VENUS', (326-Tol), (349+Tol), 5; % Венера в Водолее | Venus in
'MARS', (31), ( 56+Tol), 2; % Марс в Овне | Mars in Ari
'MERCURY', (296-Tol), (326+Tol), 1; % Меркурий в Стрельце | Mercury in Capricorn
'SUN', (296-Tol), (326+Tol), 1/3; % Солнце в Стрельце | Sun in Capricorn
'MOON', ( 56-Tol), ( 92+Tol), 1/12; % Луна в Тельце | Moon in Taurus
};
fl=0; % akfu
n=numel(params(:,1)); % количество записей
for i=1:n
target=params(i,1); degree_l=cell2mat(params(i,2))+d; degree_h=cell2mat(params(i,3))+d; step = cell2mat(params(i,4))* cspice_spd;
if (i < n) && (isequal(params(i,1),(params(i+1,1)))) % если встретился дуплет, две записи одной планеты, что проходит через 0
fl=1; % флаг, чтобы обойти следующую запись, как пару дуплета
degree_l_2=cell2mat(params(i+1,2))+d; degree_h_2=cell2mat(params(i+1,3))+d;
res_1 = cspice_gfposc ( target, 'FK4', abcr, obsrvr, crds, crd,'>', degree_l * cspice_rpd , adj, step, MAXIVL, res);
res_1 = cspice_gfposc ( target, 'FK4', abcr, obsrvr, crds, crd,'<', degree_h * cspice_rpd , adj, step, MAXIVL, res_1);
res_2 = cspice_gfposc ( target, 'FK4', abcr, obsrvr, crds, crd,'>', degree_l_2 * cspice_rpd , adj, step, MAXIVL, res);
res_2 = cspice_gfposc ( target, 'FK4', abcr, obsrvr, crds, crd,'<', degree_h_2 * cspice_rpd , adj, step, MAXIVL, res_2);
res = cspice_wnunid( res_1, res_2 );
else
if fl==0
res = cspice_gfposc ( target, 'FK4', abcr, obsrvr, crds, crd,'>', degree_l * cspice_rpd , adj, step, MAXIVL, res);
res = cspice_gfposc ( target, 'FK4', abcr, obsrvr, crds, crd,'<', degree_h * cspice_rpd , adj, step, MAXIVL, res);
else fl=0;
end
end
end
res=cspice_wnfild(cspice_spd/(12/6), res); % склеивание 5 минутных интервалов из-за перехода нуля
fprintf ( '\n\nРезультат:\n' );
for i = 1 : cspice_wncard( res )
< start, finish > = cspice_wnfetd ( res, i );
<~, ~, ~, time_st, ~> = cspice_et2lst( start, 399, lon * cspice_rpd, 'PLANETOGRAPHIC'); % Local solar time for S.
<~, ~, ~, time_end, ~> = cspice_et2lst( finish, 399, lon * cspice_rpd, 'PLANETOGRAPHIC'); % Local solar time for End
output = cspice_timout( , TIMFMT );
fprintf ( 'с %s UTC (%s LST) по %s UTC (%s LST)\n', output(1  , time_st, output(2 , time_end );
end
cspice_kclear
|
