Mohamad Magdy 0 Report post Posted November 30, 2005 · العمليات الحسابية Algorithms · قيم طول القمر Periodic Terms of longitude · قيم مسافة القمر Periodic Terms ( radius vector) · قيم عرض القمر Periodic Terms ( Latitude ) Beta · المسافة بين القمر والأرض بالكيلومتر Distance To Earth in km Δ · المسافة بين القمر والأرض بالوحده الفلكية R · المطلع المستقيم للقمر Alpha · الميل الزاوى للقمر Declination · إختلاف منظر القمر Parralax · قدر القمر Magnitude · قطر القمر الزاوى Diameter · زاوية الإستطالة Elongation · زاوية الطور Phase Angle · إختلاف منظر القمر الإستوائى Equatorial Horizontial Parallax · إختلاف منظر القمر الإستوائى المركزى EGeocentric Equatorial Horizontial Parallax · إرتفاع القمر بحسب موقع الراصد على سطح الأرض Topocentric ALtitude · الزاوية الساعية للقمر HA · إرتفاع القمر Altitude · إرتفاع القمر + عنصر الإنكسار Altitude with Refraction · زاوية سمت القمر مع مراعاة ان نقطةالجنوب = 0° Azimuth (S = 0°) · زاوية سمت القمر مع مراعاة ان نقطة الشمال = 0° Azimuth (S = 180°) · الزاوية الساعة للقمر Hour Angle [h:m:s] · إتجاه القمر فى السماء Compass Direction · إحداثيات القمر بحسب موقع الراصد على سطح الأرض Topocentr.Coordinates · المطلع المستقيم بحسب موقع الراصد على سطح الأرض Topo RA · الميل الزاوى للقمر بحسب موقع الراصد على سطح الأرض Topo Declination · إختلاف منظر القمر الحادث فى الإرتفاع Parallaxe in Altitude° · إرتفاع القمر بحسب موقع الراصد على سطح الأرض Topocentric Altitude · إرتفاع القمر بحسب موقع الراصد على سطح الأرض Topocentric Altitude + أثر الإنكسار · الزاوية الساعية بحسب موقع الراصد على سطح الأرض HA · طول وعرض القمر الجغرافى Geographic Longitude and latitude of The Moon · نسبة الضوء الحادث فى القمر Phase · العناصر المدارية للقمر Lunar Orbital Elements العناصر المطلوبة للحساب : الزمن الفلكى مصحح بقيمة Delta T JDE ميل دائرة البروج الحقيقى True Obliquity Of Ecliptic ورمزها e الترنح فى الطول Nutation in Long ورمزه ΔΨ الترنح فى ميل دائرة البروج Nutation in Obliquity ورمزهΔε الزمن النجمى العالمى الحقيقى بالكسور العشرية True GST decimal الزمن النجمى المحلى الحقيقى بالكسور العشرية True LST decimal الزمن النجمى الحقيقى العالمى للساعة 00:00 True GST UT العمليات الحسابية Algorithms TE = (JDE-2451545)/36525 L' = 218.3164591 + (481267.88134236*TE) - (0.0013268*TE^2)+ (TE^3/538841) - (TE^4/65194000) D = 297.8502042 + (445267.1115168 *TE) - (0.00163 *TE^2)+ (TE^3/545868) - (TE^4/113065) M = 357.5291092 + (35999.0502909 *TE) - (0.0001536*TE^2)+ (TE^3/24490000) M' = 134.9634114 + (477198.8676313 *TE) + (0.008997 *TE^2)+ (TE^3/69699) - (TE^4/14712000) F = 93.2720993 + (483202.0175273 *TE) - (0.0034029*TE^2)- (TE^3/3526000)+ (TE^4/863310000) A1 = 119.75 + 131.849 *TE A2 = 53.09 + 479264.29 *TE A3 = 313.45 + 481266.484 *TE E = 1-(0.002516*TE)-(0.0000074*TE^2) قيم طول القمر Periodic Terms of longitude Lambda = 6288774 * sin(0 * D + 0 * M + 1 * M' + 0 * F) + 1274027 * sin(2 * D + 0 * M - 1 * M' + 0 * F) + 658314 * sin(2 * D + 0 * M + 0 * M' + 0 * F) + 213618 * sin(0 * D + 0 * M + 2 * M' + 0 * F) + -185116 * E * sin(0 * D + 1 * M + 0 * M' + 0 * F) + -114332 * sin(0 * D + 0 * M + 0 * M' + 2 * F) + 58793 * sin(2 * D + 0 * M + -2 * M' + 0 * F) + 57066 * E * sin(2 * D + -1 * M + -1 * M' + 0 * F) + 53322 * sin(2 * D + 0 * M + 1 * M' + 0 * F) + 45758 * E * sin(2 * D + -1 * M + 0 * M' + 0 * F) + -40923 * E * sin(0 * D + 1 * M + -1 * M' + 0 * F) + -34720 * sin(1 * D + 0 * M + 0 * M' + 0 * F) + -30383 * E * sin(0 * D + 1 * M + 1 * M' + 0 * F) + 15327 * sin(2 * D + 0 * M + 0 * M' + -2 * F) + -12528 * sin(0 * D + 0 * M + 1 * M' + 2 * F) + 10980 * sin(0 * D + 0 * M + 1 * M' + -2 * F) + 10675 * sin(4 * D + 0 * M + -1 * M' + 0 * F) + 10034 * sin(0 * D + 0 * M + 3 * M' + 0 * F) + 8548 * sin(4 * D + 0 * M + -2 * M' + 0 * F) + -7888 * E * sin(2 * D + 1 * M + -1 * M' + 0 * F) + -6766 * E * sin(2 * D + 1 * M + 0 * M' + 0 * F) + -5163 * sin(1 * D + 0 * M + -1 * M' + 0 * F) + 4987 * E * sin(1 * D + 1 * M + 0 * M' + 0 * F) + 4036 * E * sin(2 * D + -1 * M + 1 * M' + 0 * F) + 3994 * sin(2 * D + 0 * M + 2 * M' + 0 * F) + 3861 * sin(4 * D + 0 * M + 0 * M' + 0 * F) + 3665 * sin(2 * D + 0 * M + -3 * M' + 0 * F) + -2689 * E * sin(0 * D + 1 * M + -2 * M' + 0 * F) + -2602 * sin(2 * D + 0 * M + -1 * M' + 2 * F) + 2390 * E * sin(2 * D + -1 * M + -2 * M' + 0 * F) + -2348 * sin(1 * D + 0 * M + 1 * M' + 0 * F) + 2236 * E ^ 2 * sin(2 * D + -2 * M + 0 * M' + 0 * F) + -2120 * E * sin(0 * D + 1 * M + 2 * M' + 0 * F) + -2069 * E ^ 2 * sin(0 * D + 2 * M + 0 * M' + 0 * F) + 2048 * E ^ 2 * sin(2 * D + -2 * M + -1 * M' + 0 * F) + -1773 * sin(2 * D + 0 * M + 1 * M' + -2 * F) + -1595 * sin(2 * D + 0 * M + 0 * M' + 2 * F) + 1215 * E * sin(4 * D + -1 * M + -1 * M' + 0 * F) + -1110 * sin(0 * D + 0 * M + 2 * M' + 2 * F) + -892 * sin(3 * D + 0 * M + -1 * M' + 0 * F) + -810 * E * sin(2 * D + 1 * M + 1 * M' + 0 * F) + 759 * E * sin(4 * D + -1 * M + -2 * M' + 0 * F) + -713 * E ^ 2 * sin(0 * D + 2 * M + -1 * M' + 0 * F) + -700 * E ^ 2 * sin(2 * D + 2 * M + -1 * M' + 0 * F) + 691 * E * sin(2 * D + 1 * M + -2 * M' + 0 * F) + 596 * E * sin(2 * D + -1 * M + 0 * M' + -2 * F) + 549 * sin(4 * D + 0 * M + 1 * M' + 0 * F) + 537 * sin(0 * D + 0 * M + 4 * M' + 0 * F) + 520 * E * sin(4 * D + -1 * M + 0 * M' + 0 * F) + -487 * sin(1 * D + 0 * M + -2 * M' + 0 * F) + -399 * E * sin(2 * D + 1 * M + 0 * M' + -2 * F) + -381 * sin(0 * D + 0 * M + 2 * M' + -2 * F) + 351 * E * sin(1 * D + 1 * M + 1 * M' + 0 * F) + -340 * sin(3 * D + 0 * M + -2 * M' + 0 * F) + 330 * sin(4 * D + 0 * M + -3 * M' + 0 * F) + 327 * E * sin(2 * D + -1 * M + 2 * M' + 0 * F) + -323 * E ^ 2 * sin(0 * D + 2 * M + 1 * M' + 0 * F) + 299 * E * sin(1 * D + 1 * M + -1 * M' + 0 * F) + 294 * sin(2 * D + 0 * M + 3 * M' + 0 * F) + 0 * sin(2 * D + 0 * M + -1 * M' + -2 * F) قيم مسافة القمر Periodic Terms ( radius vector) R = -20905355 * cos(0 * D + 0 * M + 1 * M' + 0 * F) + -3699111 * cos(2 * D + 0 * M + -1 * M' + 0 * F) + -2955968 * cos(2 * D + 0 * M + 0 * M' + 0 * F) + -569925 * cos(0 * D + 0 * M + 2 * M' + 0 * F) + 48888 * E * cos(0 * D + 1 * M + 0 * M' + 0 * F) + -3149 * cos(0 * D + 0 * M + 0 * M' + 2 * F) + 246158 * cos(2 * D + 0 * M + -2 * M' + 0 * F) + -152138 * E * cos(2 * D + -1 * M + -1 * M') + -170733 * cos(2 * D + 0 * M + 1 * M' + 0 * F) + -204586 * E * cos(2 * D + -1 * M + 0 * M' + 0 * F) + -129620 * E * cos(0 * D + 1 * M + -1 * M' + 0 * F) + 108743 * cos(1 * D + 0 * M + 0 * M' + 0 * F) + 104755 * E * cos(0 * D + 1 * M + 1 * M' + 0 * F) + 10321 * cos(2 * D + 0 * M + 0 * M' + -2 * F) + 0 * cos(0 * D + 0 * M + 1 * M' + 2 * F) + 79661 * cos(0 * D + 0 * M + 1 * M' + -2 * F) + -34782 * cos(4 * D + 0 * M + -1 * M' + 0 * F) + -23210 * cos(0 * D + 0 * M + 3 * M' + 0 * F) + -21636 * cos(4 * D + 0 * M + -2 * M' + 0 * F) + 24208 * E * cos(2 * D + 1 * M + -1 * M' + 0 * F) + 30824 * E * cos(2 * D + 1 * M + 0 * M' + 0 * F) + -8379 * cos(1 * D + 0 * M + -1 * M' + 0 * F) + -16675 * E * cos(1 * D + 1 * M + 0 * M' + 0 * F) + -12831 * E * cos(2 * D + -1 * M + 1 * M' + 0 * F) + -10445 * cos(2 * D + 0 * M + 2 * M' + 0 * F) + -11650 * cos(4 * D + 0 * M + 0 * M' + 0 * F) + 14403 * cos(2 * D + 0 * M + -3 * M' + 0 * F) + -7003 * E * cos(0 * D + 1 * M + -2 * M' + 0 * F) + 0 * cos(2 * D + 0 * M + -1 * M' + 2 * F) + 10056 * E * cos(2 * D + -1 * M + -2 * M' + 0 * F) + 6322 * cos(1 * D + 0 * M + 1 * M' + 0 * F) + -9884 * E ^ 2 * cos(2 * D + -2 * M + 0 * M' + 0 * F) + 5751 * E * cos(0 * D + 1 * M + 2 * M' + 0 * F) + 0 * E ^ 2 * cos(0 * D + 2 * M + 0 * M' + 0 * F) + -4950 * E ^ 2 * cos(2 * D + -2 * M + -1 * M' + 0 * F) + 4130 * cos(2 * D + 0 * M + 1 * M' + -2 * F) + 0 * cos(2 * D + 0 * M + 0 * M' + 2 * F) + -3958 * E * cos(4 * D + -1 * M + -1 * M' + 0 * F) + 0 * cos(0 * D + 0 * M + 2 * M' + 2 * F) + 3258 * cos(3 * D + 0 * M + -1 * M' + 0 * F) + 2616 * E * cos(2 * D + 1 * M + 1 * M' + 0 * F) + -1897 * E * cos(4 * D + -1 * M + -2 * M' + 0 * F) + -2117 * E ^ 2 * cos(0 * D + 2 * M + -1 * M' + 0 * F) + 2354 * E ^ 2 * cos(2 * D + 2 * M + -1 * M' + 0 * F) + 0 * E * cos(2 * D + 1 * M + -2 * M' + 0 * F) + 0 * E * cos(2 * D + -1 * M + 0 * M' + -2 * F) + -1423 * cos(4 * D + 0 * M + 1 * M' + 0 * F) + -1117 * cos(0 * D + 0 * M + 4 * M' + 0 * F) + -1571 * E * cos(4 * D + -1 * M + 0 * M' + 0 * F) + -1739 * cos(1 * D + 0 * M + -2 * M' + 0 * F) + 0 * E * cos(2 * D + 1 * M + 0 * M' + -2 * F) + -4421 * cos(0 * D + 0 * M + 2 * M' + -2 * F) + 0 * E * cos(1 * D + 1 * M + 1 * M' + 0 * F) + 0 * cos(3 * D + 0 * M + -2 * M' + 0 * F) + 0 * cos(4 * D + 0 * M + -3 * M' + 0 * F) + 0 * E * cos(2 * D + -1 * M + 2 * M' + 0 * F) + 1165 * E ^ 2 * cos(0 * D + 2 * M + 1 * M' + 0 * F) + 0 * E * cos(1 * D + 1 * M + -1 * M' + 0 * F) + 0 * cos(2 * D + 0 * M + 3 * M' + 0 * F) + 8752 * cos(2 * D + 0 * M + -1 * M' + -2 * F) قيم عرض القمر Periodic Terms ( Latitude ) Beta Beta = 5128122 * Sin(0 * D + 0 * M + 0 * M' + 1 * F) + 280602 * Sin(0 * D + 0 * M + 1 * M' + 1 * F) + 277693 * Sin(0 * D + 0 * M + 1 * M' + -1 * F) + 173237 * Sin(2 * D + 0 * M + 0 * M' + -1 * F) + 55413 * Sin(2 * D + 0 * M + -1 * M' + 1 * F) + 46271 * Sin(2 * D + 0 * M + -1 * M' + -1 * F) + 32573 * Sin(2 * D + 0 * M + 0 * M' + 1 * F) + 17198 * Sin(0 * D + 0 * M + 2 * M' + 1 * F) + 9266 * Sin(2 * D + 0 * M + 1 * M' + -1 * F) + 8822 * Sin(0 * D + 0 * M + 2 * M' + -1 * F) + 8216 * E * Sin(2 * D + -1 * M + 0 * M' + -1 * F) + 4324 * Sin(2 * D + 0 * M + -2 * M' + -1 * F) + 4200 * Sin(2 * D + 0 * M + 1 * M' + 1 * F) + -3359 * E * Sin(2 * D + 1 * M + 0 * M' + -1 * F) + 2463 * E * Sin(2 * D + -1 * M + -1 * M' + 1 * F) + 2211 * E * Sin(2 * D + -1 * M + 0 * M' + 1 * F) + 2065 * E * Sin(2 * D + -1 * M + -1 * M' + -1 * F) + -1870 * E * Sin(0 * D + 1 * M + -1 * M' + -1 * F) + 1828 * Sin(4 * D + 0 * M + -1 * M' + -1 * F) + -1794 * E * Sin(0 * D + 1 * M + 0 * M' + 1 * F) + -1749 * Sin(0 * D + 0 * M + 0 * M' + 3 * F) + -1565 * E * Sin(0 * D + 1 * M + -1 * M' + 1 * F) + -1491 * Sin(1 * D + 0 * M + 0 * M' + 1 * F) + -1475 * E * Sin(0 * D + 1 * M + 1 * M' + 1 * F) + -1410 * E * Sin(0 * D + 1 * M + 1 * M' + -1 * F) + -1344 * E * Sin(0 * D + 1 * M + 0 * M' + -1 * F) + -1335 * Sin(1 * D + 0 * M + 0 * M' + -1 * F) + 1107 * Sin(0 * D + 0 * M + 3 * M' + 1 * F) + 1021 * Sin(4 * D + 0 * M + 0 * M' + -1 * F) + 833 * Sin(4 * D + 0 * M + -1 * M' + 1 * F) + 777 * Sin(0 * D + 0 * M + 1 * M' + -3 * F) + 671 * Sin(4 * D + 0 * M + -2 * M' + 1 * F) + 607 * Sin(2 * D + 0 * M + 0 * M' + -3 * F) + 596 * Sin(2 * D + 0 * M + 2 * M' + -1 * F) + 491 * E * Sin(2 * D + -1 * M + 1 * M' + -1 * F) + -451 * Sin(2 * D + 0 * M + -2 * M' + 1 * F) + 439 * Sin(0 * D + 0 * M + 3 * M' + -1 * F) + 422 * Sin(2 * D + 0 * M + 2 * M' + 1 * F) + 421 * Sin(2 * D + 0 * M + -3 * M' + -1 * F) + -366 * E * Sin(2 * D + 1 * M + -1 * M' + 1 * F) + -351 * E * Sin(2 * D + 1 * M + 0 * M' + 1 * F) + 331 * Sin(4 * D + 0 * M + 0 * M' + 1 * F) + 315 * E * Sin(2 * D + -1 * M + 1 * M' + 1 * F) + 302 * E ^ 2 * Sin(2 * D + -2 * M + 0 * M' + -1 * F) + -283 * Sin(0 * D + 0 * M + 1 * M' + 3 * F) + -229 * E * Sin(2 * D + 1 * M + 1 * M' + -1 * F) + 223 * E * Sin(1 * D + 1 * M + 0 * M' + -1 * F) + 223 * E * Sin(1 * D + 1 * M + 0 * M' + 1 * F) + -220 * E * Sin(0 * D + 1 * M + -2 * M' + -1 * F) + -220 * E * Sin(2 * D + 1 * M + -1 * M' + -1 * F) + -185 * Sin(1 * D + 0 * M + 1 * M' + 1 * F) + 181 * E * Sin(2 * D + -1 * M + -2 * M' + -1 * F) + -177 * E * Sin(0 * D + 1 * M + 2 * M' + 1 * F) + 176 * Sin(4 * D + 0 * M + -2 * M' + -1 * F) + 166 * E * Sin(4 * D + -1 * M + -1 * M' + -1 * F) + -164 * Sin(1 * D + 0 * M + 1 * M' + -1 * F) + 132 * Sin(4 * D + 0 * M + 1 * M' + -1 * F) + -119 * Sin(1 * D + 0 * M + -1 * M' + -1 * F) + 115 * E * Sin(4 * D + -1 * M + 0 * M' + -1 * F) + 107 * E ^ 2 * Sin(2 * D + -2 * M + 0 * M' + 1 * F) AD Sl = 3958 * SIN(A1)+1962 * SIN(L'-F)+ 318 * SIN(A2) SL = Sum lambda terms + AD SL Lambda = ( L'+ SL /1000000 ) Lambda = lambda + Nutation in longitude AD SB = - SIN (L') *2235 + SIN (A3) *382 + SIN (A1-F) *175 + SIN (F+A1) *175 + SIN (L'-M')*127 + - SIN (L'+M')*115 SB = Sum Beta terms + AD SB Beta = SB /1000000 True Beta = Beta + Nutation in obliquity المسافة بين القمر والأرض بالكيلومتر Distance To Earth in km Δ SR = Sum R terms Distance To Earth in KM Delta = 385000.56 +( SR /1000) المسافة بين القمر والأرض بالوحده الفلكية R R = SQRT( r * r + (Delta*1/149597870.706008)*(Delta*1/149597870.706008) - COS (Elongation)*(Delta*1/149597870.706008)* r *2) المطلع المستقيم للقمر Alpha A = sin (lamda)* cos(ob) - tan(beta) * sin(ob) B = cos (lamda) Alpha = ATAN2( B ; A ) Alpha hrs = Alpha / 15 الميل الزاوى للقمر Declination Dec = sin (beta)* cos(ob)+ cos(beta)* sin(ob) * sin(lambda ) Dec = Asin (Dec) إختلاف منظر القمر Parralax P = 0.9983271 + 0.001674 * COS(2*lat)- 0.0000035 * COS(4*lat) قدر القمر Magnitude Mag = (-0.23 + 5 * LOG (R * Delta *((6.684587122 *10^(-12))*1000)) + 0.026 * Phase angle + 4 * 10^-9 * Phase angle ^4) + 0.44 قطر القمر الزاوى Diameter Diameter = 2*((1+(6378.14/Delta)*SIN(Alt))* 358473400 /Delta ) زاوية الإستطالة Elongation Elongation = ABS(ACOS(COS(lamda Moon - lambda Sun)*COS(Beta))) زاوية الطور Phase Angle Phase Angle = ACOS (-COS(Elongation)) إختلاف منظر القمر الإستوائى Equatorial Horizontial Parallax Equ.-Horizontalparallax p° = ASIN(SIN(8.794/3600)/(Delta * 0.000000006684587122)) إختلاف منظر القمر الإستوائى المركزى Geocentric Equatorial Horizontial Parallax Geocentric Equatorial-Horizontal-Parallax: = ASIN((SIN(8.794/3600))/(Delta/149597892)) إرتفاع القمر بحسب موقع الراصد على سطح الأرض Topocentric ALtitude Alt = ASIN(SIN(Lat) * SIN(Dec)+ COS(Lat)* COS(Dec) * COS(H)) P = 0.9983271 + 0.001674 * COS(2*Lat)- 0.0000035 * COS(4*Lat) Top Alt = Alt-p * COS(alt) * Equatorial-Horizontal-Parallaxe الزاوية الساعية للقمر HA Co Lat = 360-lat IF(((((AGST-Colat-alpha°)/360)-ROUNDDOWN((AGST-Colat-alpha°)/360;0))*360)<0 Then H = 360+((((AGST-Colat-alpha°)/360)-ROUNDDOWN((AGST-Colat-alpha°)/360;0))*360) Else H = ((((AGST-Colat-alpha°)/360)-ROUNDDOWN((AGST-Colat-alpha°)/360;0))*360)) إرتفاع القمر Altitude Alt =SIN(Lat)*SIN(dec)+COS(lat)*COS(dec)*COS(H) Alt = Asin ( Alt ) إرتفاع القمر + عنصر الإنكسار Altitude with Refraction =IF(ALT<=-1;ALT;(((((1.02/(TAN((ALT+(10.3/(ALT+5.11)))*PI()/180)))+0.0019279)*(1013/1010)*(283/(273+10)))/60)+ALT)) زاوية سمت القمر مع مراعاة ان نقطةالجنوب = 0° Azimuth (S = 0°) =ATAN2((COS(H)*SIN(lat)-TAN(dec)*COS(lat));(SIN(H))) زاوية سمت القمر مع مراعاة ان نقطة الشمال = 0° Azimuth (S = 180°) =IF(Azimuth>=180;Azimuth-180;180+Azimuth) الزاوية الساعة للقمر Hour Angle [h:m:s] HA = LST + (24 - RA) إتجاه القمر فى السماء Moon Compass Direction in The Sky Az = Azimuth ( N = 0° ) IF Az >=180 AND <=191.25 Then Moon Compass Direction is "South" IF Az >191.25 AND <=213.75 Then Moon Compass Direction is "South South West" IF Az >213.75 AND <=236.25 Then Moon Compass Direction is "South West " IF Az >236.25 AND <=258.75 Then Moon Compass Direction is "West South West " IF Az >258.75 AND <=281.25 Then Moon Compass Direction is "West " IF Az >281.25 AND <=303.75 Then Moon Compass Direction is "West North West " IF Az >303.75 AND <=326.25 Then Moon Compass Direction is "North West " IF Az >326.25 AND <=348.75 Then Moon Compass Direction is "North North West " IF Az >348.75 AND <=360 Then Moon Compass Direction is "North " IF Az >=0 AND <=11.25 Then Moon Compass Direction is "North " IF Az >11.25 AND <=33.75 Then Moon Compass Direction is "North North East" IF Az>33.75 AND <=56.25 Then Moon Compass Direction is "North East " IF Az>56.25 AND <=78.75 Then Moon Compass Direction is "East North East " IF Az>78.75 AND <=101.25 Then Moon Compass Direction is "East " IF Az>101.25 AND <=123.75 Then Moon Compass Direction is "East South East " IF Az>123.75 AND <=146.25 Then Moon Compass Direction is "South East " IF Az>146.25 AND <=168.75 Then Moon Compass Direction is "South South East " IF Az>168.75 AND <=180 Then Moon Compass Direction is "South " إحداثيات القمر بحسب موقع الراصد على سطح الأرض Topocentr.Coordinates Pcos H = m Lat = ° a = 6378.14 f =1/298.257 b = a *(1-f) b/a = b / a U =(b / a ) * TAN( lat ) U = ATAN(U) psin =(b/a) * SIN( U ) pcos =COS ( U ) Sin Pi = sin ( pi ) Pi = Parallax Parallax = 0.00244277777777778 / R X = 0.000156785520543656/1000 X = hight * x P = (0.9983271 + 0.0016764 * COS(2*lat)- 0.0000035 * COS(4*lat)) + X المطلع المستقيم بحسب موقع الراصد على سطح الأرض Topo RA Delta RA = (-pcos * sin pi *SIN(HA))/(COS(dec)-pcos*sin pi *COS(HA)) Delta RA = ATAN(Delta RA) Topo RA = Delta RA + RA الميل الزاوى للقمر بحسب موقع الراصد على سطح الأرض Topo Declination Delta dec = ((SIN(dec) – psin * sin(Pi)) * COS(delta RA ))/(COS(dec)- pcos* sin(Pi) * COS( HA )) Topo Dec = ATAN( Delta dec ) إختلاف منظر القمر الحادث فى الإرتفاع Parallaxe in Altitude° =(0.9983271+0.001674*COS(2*lat)-0.0000035*COS(4*Lat)) * Equ.Horizontalparallaxep * COS (geocentric altitude) إرتفاع القمر بحسب موقع الراصد على سطح الأرض Topocentric Altitude = -(Parallax in Altitude)+ geocentric Altitude إرتفاع القمر بحسب موقع الراصد على سطح الأرض Topocentric Altitude + أثر الإنكسار =IF(Alt<=-1;Alt;(((((1.02/(TAN(Alt+(10.3/(Alt+5.11)))))+0.0019279)*(1013/1010)*(283/(273+10)))/60)+Alt)) الزاوية الساعية بحسب موقع الراصد على سطح الأرض HA = LST + topoRA طول وعرض القمر الجغرافى Geographic Longitude and latitude of The Moon Geocentric Lat = dec Geo Lat = dec + 0.1924 * SIN(2* dec) for geoc & geo lat : Lat = ROUND(lat;3)&"°" direction = IF(lat>=0;"N";"S") Lat = LAt& direction Geo long = IF(( RA-GST)<-180;360+(RA-GST);IF((RA-GST)>180;-360+(RA-GST);(RA-GST))) Direction= IF (Geo long<0;"West";"East") Geo long = Geo long & direction Longitude Location = IF(Long>180;-360+Long;Long) Latitude Location = latitude نسبة الضوء الحادث فى القمر Phase RA Moon Dec Moon Dis Earth-Moon km RA Sun Dec Sun Distance Earth - Sun AU R = Dis Earth - Sun AU *(149971500/1.00249769) F = ACOS(SIN(Sun Dec)*SIN(Moon Dec) + COS(Sun Dec)*COS(Moon Dec)* COS(Sun RA-Moon RA)) I = ( R * SIN( F )) / ((Earth Moon Dis - R * COS (F))) I = Atan (I) K = (1 + COS(I)) /2 K = (1 - K ) * 100 A = 100 * k B = (1 – K ) * 100 Phase = IF ( F > 90 ; A ; B ) العناصر المدارية للقمر Lunar Orbital Elements T = JDE - 2451543.5 M = e = 0.0459 i = 5.1454 N = ascending Moon Node w = 318.0634 + 0.1643573223* t M = 115.3654 + 13.0649929509 *t a = 0.389 * 1000000 * 0.000000006684587122 pi = L = M + W + N دَعْوَاهُمْ فِيهَا سُبْحَانَكَ اللَّهُمَّ وَتَحِيَّتُهُمْ فِيهَا سَلامٌ وَآخِرُ دَعْوَاهُمْ أَنِ الْحَمْدُ لِلَّهِ رَبِّ الْعَالَمِينَ وصلىاللهم وسلم وبارك علىسيدنا محمد وعلى آله وصحبه وسلم Share this post Link to post Share on other sites
Mohamad Magdy 0 Report post Posted December 1, 2005 مرفق ملف txt به المعادلات مرتبه للفائدة وشكرا محمد مجدى Moon_Position.txt Share this post Link to post Share on other sites
العبادي 0 Report post Posted December 1, 2005 اخي الكريم محمد مجدي اننا نشكرك على هذه الدروس القيمة الشاملة ونأمل منكم ايضاح معنى الإختصارات الفلكية الموجودة في الدرس الأول الذي يدور عن حساب مواقع الكزاكب حول الشمس وفي البروج مثل Jd w a e v r وبقية الإختصارات الأخرى ولكم جزيل الشكر وجزاكم الله خيرا كثيرا Share this post Link to post Share on other sites
Mohamad Magdy 0 Report post Posted December 1, 2005 الأخ العزيز العبادى المصطلحات التى كتبتها هى العناصر المدارية للقمر على العموم اى مصطلح ارجوا ان تكتبه لأن المصطلحات كثيره جدا ولا ادرى ايهما تود إيضاحه وساقوم بالرد بعون الله وشكراً Share this post Link to post Share on other sites
elgavry osama 0 Report post Posted December 7, 2005 بسم الله الرحمن الرحيم استاذى العزيز: جزاكم الله خيرا عن مجهودكم الكبير الذي تبذله فى نشر العلم النافع وارجو منكم ان تلحق حسابات الكسوف والخسوف مع حسابات القمر اسامة فتحي Share this post Link to post Share on other sites
العبادي 0 Report post Posted December 8, 2005 اخي الفاضل/ محمد مجدي ..........................حفظه الله ورعاه لقد ذكرتم الآتي فما معنى تلك القيم مع توضيح ذلك بمثال من واقع اى تاريخ وليكن تاريخ هذا الشهر مثلا حتى تتضح الصورة Lambda = 6288774 * sin(0 * D + 0 * M + 1 * M' + 0 * F) + 1274027 * sin(2 * D + 0 * M - 1 * M' + 0 * F) + 658314 * sin(2 * D + 0 * M + 0 * M' + 0 * F) + 213618 * sin(0 * D + 0 * M + 2 * M' + 0 * F) + -185116 * E * sin(0 * D + 1 * M + 0 * M' + 0 * F) + -114332 * sin(0 * D + 0 * M + 0 * M' + 2 * F) + 58793 * sin(2 * D + 0 * M + -2 * M' + 0 * F) + 57066 * E * sin(2 * D + -1 * M + -1 * M' + 0 * F) + 53322 * sin(2 * D + 0 * M + 1 * M' + 0 * F) + 45758 * E * sin(2 * D + -1 * M + 0 * M' + 0 * F) + -40923 * E * sin(0 * D + 1 * M + -1 * M' + 0 * F) + -34720 * sin(1 * D + 0 * M + 0 * M' + 0 * F) + -30383 * E * sin(0 * D + 1 * M + 1 * M' + 0 * F) + 15327 * sin(2 * D + 0 * M + 0 * M' + -2 * F) + -12528 * sin(0 * D + 0 * M + 1 * M' + 2 * F) + 10980 * sin(0 * D + 0 * M + 1 * M' + -2 * F) + 10675 * sin(4 * D + 0 * M + -1 * M' + 0 * F) + 10034 * sin(0 * D + 0 * M + 3 * M' + 0 * F) + 8548 * sin(4 * D + 0 * M + -2 * M' + 0 * F) + -7888 * E * sin(2 * D + 1 * M + -1 * M' + 0 * F) + -6766 * E * sin(2 * D + 1 * M + 0 * M' + 0 * F) + -5163 * sin(1 * D + 0 * M + -1 * M' + 0 * F) + 4987 * E * sin(1 * D + 1 * M + 0 * M' + 0 * F) + 4036 * E * sin(2 * D + -1 * M + 1 * M' + 0 * F) + 3994 * sin(2 * D + 0 * M + 2 * M' + 0 * F) + 3861 * sin(4 * D + 0 * M + 0 * M' + 0 * F) + 3665 * sin(2 * D + 0 * M + -3 * M' + 0 * F) + -2689 * E * sin(0 * D + 1 * M + -2 * M' + 0 * F) + -2602 * sin(2 * D + 0 * M + -1 * M' + 2 * F) + 2390 * E * sin(2 * D + -1 * M + -2 * M' + 0 * F) + -2348 * sin(1 * D + 0 * M + 1 * M' + 0 * F) + 2236 * E ^ 2 * sin(2 * D + -2 * M + 0 * M' + 0 * F) + -2120 * E * sin(0 * D + 1 * M + 2 * M' + 0 * F) + -2069 * E ^ 2 * sin(0 * D + 2 * M + 0 * M' + 0 * F) + 2048 * E ^ 2 * sin(2 * D + -2 * M + -1 * M' + 0 * F) + -1773 * sin(2 * D + 0 * M + 1 * M' + -2 * F) + -1595 * sin(2 * D + 0 * M + 0 * M' + 2 * F) + 1215 * E * sin(4 * D + -1 * M + -1 * M' + 0 * F) + -1110 * sin(0 * D + 0 * M + 2 * M' + 2 * F) + -892 * sin(3 * D + 0 * M + -1 * M' + 0 * F) + -810 * E * sin(2 * D + 1 * M + 1 * M' + 0 * F) + 759 * E * sin(4 * D + -1 * M + -2 * M' + 0 * F) + -713 * E ^ 2 * sin(0 * D + 2 * M + -1 * M' + 0 * F) + -700 * E ^ 2 * sin(2 * D + 2 * M + -1 * M' + 0 * F) + 691 * E * sin(2 * D + 1 * M + -2 * M' + 0 * F) + 596 * E * sin(2 * D + -1 * M + 0 * M' + -2 * F) + 549 * sin(4 * D + 0 * M + 1 * M' + 0 * F) + 537 * sin(0 * D + 0 * M + 4 * M' + 0 * F) + 520 * E * sin(4 * D + -1 * M + 0 * M' + 0 * F) + -487 * sin(1 * D + 0 * M + -2 * M' + 0 * F) + -399 * E * sin(2 * D + 1 * M + 0 * M' + -2 * F) + -381 * sin(0 * D + 0 * M + 2 * M' + -2 * F) + 351 * E * sin(1 * D + 1 * M + 1 * M' + 0 * F) + -340 * sin(3 * D + 0 * M + -2 * M' + 0 * F) + 330 * sin(4 * D + 0 * M + -3 * M' + 0 * F) + 327 * E * sin(2 * D + -1 * M + 2 * M' + 0 * F) + -323 * E ^ 2 * sin(0 * D + 2 * M + 1 * M' + 0 * F) + 299 * E * sin(1 * D + 1 * M + -1 * M' + 0 * F) + 294 * sin(2 * D + 0 * M + 3 * M' + 0 * F) + 0 * sin(2 * D + 0 * M + -1 * M' + -2 * F) Share this post Link to post Share on other sites
سالم الجعيدي 0 Report post Posted May 10, 2006 يا جماعة ، هل تريدون أن أشرح لكم القانون الذي شرحه الأستاذ محمد مجدي في حساب القمر أنا مستعد هل لديكم همة ؟؟؟ أم أنكم تعتمدون على الأشياء الجاهزة للأستهلاك سنرى Share this post Link to post Share on other sites
حيدر الموسوي 0 Report post Posted September 25, 2006 حقيقة اخي العزيز انت ما قصرت في هذا العرض الشيق وانا ان شاء الله من المستفيدين وسوف اقدم لك اي سؤال في حالة حصول حاجة لذلك . تحياتي . Share this post Link to post Share on other sites
أحمد الأنصاري 0 Report post Posted October 6, 2006 والله يا سالم الجعيدي احنا ما نقول لأ Share this post Link to post Share on other sites
فخري حسن 0 Report post Posted November 22, 2006 بصراحه شكلك تعبت كثير و الله يعطيك العافية بس شكل الموضوع أعلى من المستوى و سؤالي هو : كيف تعلل رؤيتنا لوجه واحد من القمر و شكرا لكم Share this post Link to post Share on other sites
