ВНИМАНИЕ! На форуме началось голосование в конкурсе - астрофотография месяца АПРЕЛЬ!
0 Пользователей и 1 Гость просматривают эту тему.
$Year = 2000; $Mon = 1; $Day = 20; $hous_zona = 6; $min = 0; $sec = 0; $lon = 73.3682; $Var1 = 10000 * $Year + 100 * $Mon + $Day; if($Mon <= 2 ) { $Mon = $Mon + 12; $Year = $Year - 1; } if($Var1 <= 15821004){ $Var2 = -2 + floor(($Year + 4716) / 4) - 1179; }else{ $Var2 = floor($Year /400) - floor($Year / 100) + floor($Year / 4); } $Var3 = 365 * $Year - 679004; // MJD - Модифицированная Юлианская дата $MJD = $Var3 + $Var2 + floor(306001 * ($Mon + 1)/ 10000) + $Day; // Вычисление Гринвеческого звездного времени $T0 = ($MJD - 51544.5) / 36525; // мод.юл.дата на начало суток в юлианских столетиях echo $T0; /* */ $a1 = 24110.54841; $a2 = 8640184.812; $a3 = 0.093104; $a4 = 0.0000062; $S0 = $a1 + $a2 * $T0 + $a3 * $T0*$T0- $a4 * $T0 *$T0*$T0;// звездное время в Гринвиче на начало суток в секундах //UT - всемирное время в часах, момент расчета $UT = $hous_zona + $min/60 + $sec/3600; if($UT>24)$UT=$UT-24; if($UT<0)$UT=$UT+24; $Nsec = $UT * 3600; // количество секунд, прошедших от начала суток до момента наблюдения $NsecS = $Nsec * 366.2422 / 365.2422; //количество звездных секунд $GMT = ($S0 + $NsecS) /3600 * 15;//гринвическое среднее звездное время в градусах SG while($GMT>360)$GMT=$GMT-360; $GST = $GMT + $lon;// местное звездное время ST //Lon – долгота наблюдателя // Вычисление эклиптических координат Солнца $T0 = ($MJD - 51544.5) / 36525; // мод.юл.дата на начало суток в юлианских столетиях $M = 357.528 + 35999.05 * $T0 + 0.04107 * $UT;// средняя аномалия while($M>360)$M=$M-360; $L0 = 280.46 + 36000.772 * $T0 + 0.04107 * $UT; $L = $L0 + (1.915 - 0.0048 * $T0) * sin(rad2deg($M)) + 0.02 * sin(rad2deg(2 *$M));//долгота Солнца while($L>360)$L=$L-360; echo "\n"; echo $L;
А почему угол к примеру -30 градусов не должен считаться легитимным?Что такое периодическая функция помните?
Но вводя года менее 2000 года получаем отрицательный результат. А при дате 20 января должно получаться где-то 300 градусов. Почему?