100030, RE: не вполне адекватный метод Послано guest, 21-04-2013 14:12
>С учётом разброса протяжённости созвездий избранный метод >является ошибочным. Он бездумно перенесён из анализа >знакового гороскопа (где все знаки имеют одинаковую >протяжённость 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
|