月球位置计算Word格式.docx

1、 近日点平均离角 M 利用级数展开的克卜勒方程式（Keplers Equation）求得真实的近日点离角 M360 E sin（M）/900 E2 sin（2M）/4180E3sin（M）/ 4 太阳的黄道经度 月球的轨道资料（以黄道面来定义） 月球的平均经度 L218.31664613.1* JD10.0014664 T2 月球的平均近日点离角 MmL（83.3532430.1* JD10.0103217 T2） 升交点引数 N125.0445550.1* JD10.0020756 T2 月面轨道倾角 I5.15668983 Ev1.2739 sin（ 2L-2Mm ） Ae0.1858 s

2、in（ M ） A30.37 sin（ M ） 修正月球的平均近日点离角 MmMmEvAeA3 Ec6.2886 sin（ Mm ） A40.214 sin（ 2Mm 修正月球的平均经度 LLEvEcAeA4 V0.6583 sin（ 2L2） 再次修正月球的平均经度 LLV 修正升交点引数 NN0.16 sin（ M ） 计算月球的黄道经度与纬度 sin（）sin（ LN ） cos（ I ） cos（）cos（ L 判断的象限后，求得。 道经度 mN 月球的黄道纬度 msin-1 sin（ L ） sin（ I ） 若需要月球的赤道经度与纬度，可参考座标系统的变换部分。将黄道座标转成赤道座

3、标。或者更进一步转成地平方为角。 2000年9月1日 AM10:30 月球的黄道经纬度 、赤道经纬度 与地平方为角 查得当年力学时与平均太阳时的差值 T64秒 ET=UT+64秒 计算ET儒略日得到 JD2451788.60492 JD1243.60492 儒略世纪 T0.00667 520.492123 282.948812 E0.01671 太阳黄道位置计算结果 M237.543311 158.8947 0 月球的平均经度 L188.3759761 月球的平均近日点离角 Mm77.88239044 升交点引数 N112.14385644 Ev0.41717013 Ae0.15695074

4、A30.3125499277.93472097 Ec6.14968536 A40.08748689194.35292896 V0.62138281194.97431177112.27901316 sin（）0.98786949 cos（）0.127146 判断的象限后，求得82.66593614 月球的黄道经度 m194.9449493 月球的黄道纬度 m5.11472622 经过黄道座标与赤道座标系统后 13h 2m56s 110 7 经过地平座标转换后 白天也可以看到月球！ 若是经过多项的摄动计算后，月球的位置会更精确，它的赤道座标位置： 13h 1m24s 2 31 计算月球位置的c语言

5、程序 / / 名称： / 作者：胡铂（http:/hubble.lamost.org） / 日期：20210929 / 说明：根据北京天文同好会提供的Astronomy Algrithms翻译版实现，并得到了 / zjuglr的帮助。 / #include math.hstdio.h #define DE 3.141592654/180 /计算儒略日历书时/ float jde（int Y,int M,int D,int hour,int min,int sec） int f,g; double mid1,mid2,J,JDE,A; if（M=3） if（M=1|M=2） g=M+12; mi

6、d1=floor（365.25*f）; mid2=floor（30.6001*（g+1）; A=2-floor（f/100）+floor（f/400）; J=mid1+mid2+D+A+1720994.5; JDE=J+hour/24+min/1440+sec/86400; return JDE; / void main（void） /变量定义/ int i,year,month,day,hour,min,sec; double JDE,T,L1,D,M,M1,F,A1,A2,A3,E,SUML,lamda,SUMB,beta,SUMR,SIN1,SIN2,COS1,Dist; / /数据/

7、static double La60=0,2,2,0,0,0,2,2,2,2,0,1,0,2,0,0,4,0,4,2,2,1,1,2,2,4,2,0,2,2,1,2,0,0,2,2,2,4,0,3,2,4,0,2,2,2,4,0,4,1,2,0,1,3,4,2,0,1,2,2; static double Lb60=0,0,0,0,1,0,0,-1,0,-1,1,0,1,0,0,0,0,0,0,1,1,0,1,-1,0,0,0,1,0,-1,0,-2,1,2,-2,0,0,-1,0,0,1,-1,2,2,1,-1,0,0,-1,0,1,0,1,0,0,-1,2,1,0,0; static do

8、uble Lc60=1,-1,0,2,0,0,-2,-1,1,0,-1,0,1,0,1,1,-1,3,-2,-1,0,-1,0,1,2,0,-3,-2,-1,-2,1,0,2,0,-1,1,0,-1,2,-1,1,-2,-1,-1,-2,0,1,4,0,-2,0,2,1,-2,-3,2,1,-1,3,-1; static double Ld60=0,0,0,0,0,2,0,0,0,0,0,0,0,-2,2,-2,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,0,0,0,-2,2,0,2,0,0,0,0,0,0,-2,0,0,0,0,-2,-2,0,0,0,0,0,0,0,-2

9、; static double Sl60=6288774,1274027,658314,213618,-185116,-114332,58793,57066,53322,45758,-40923,-34720,-30383,15327,-12528,10980,10675,10034,8548,-7888,-6766,-5163,4987,4036,3994,3861,3665,-2689,-2602,2390,-2348,2236,-2120,-2069,2048,-1773,-1595,1215,-1110,-892,-810,759,-713,-700,691,596,549,537,5

10、20,-487,-399,-381,351,-340,330,327,-323,299,294,0; static double Sr60=-20905355,-3699111,-2955968,-569925,48888,-3149,246158,-152138,-170733,-204586,-129620,108743,104755,10321,0,79661,-34782,-23210,-21636,24208,30824,-8379,-16675,-12831,-10445,-11650,14403,-7003,0,10056,6322,-9884,5751,0,-4950,4130

11、,0,-3958,0,3258,2616,0,-2117,2354,0,0,0,0,0,0,0,-4421,0,0,0,0,1165,0,0,8752; static double Sb60=5128122,280602,277693,173237,55413,46271,32573,17198,9266,8822,8216,4324,4200,-3359,2463,2211,2065,-1870,1828,-1794,-1749,-1565,-1491,-1475,-1410,-1344,-1335,1107,1021,833,777,671,607,596,491,-451,439,422

12、,421,-366,-351,331,315,302,-283,-229,223,223,-220,-220,-185,181,-177,176,166,-164,132,-119,115,107; static double Ba60=0,0,0,2,2,2,2,0,2,0,2,2,2,2,2,2,2,0,4,0,0,0,1,0,0,0,1,0,4,4,0,4,2,2,2,2,0,2,2,2,2,4,2,2,0,2,1,1,0,2,1,2,0,4,4,1,4,1,4,2; static double Bb60=0,0,0,0,0,0,0,0,0,0,-1,0,0,1,-1,-1,-1,1,0

13、,1,0,1,0,1,1,1,0,0,0,0,0,0,0,0,-1,0,0,0,0,1,1,0,-1,-2,0,1,1,1,1,1,0,-1,1,0,-1,0,0,0,-1,-2; static double Bc60=0,1,1,0,-1,-1,0,2,1,2,0,-2,1,0,-1,0,-1,-1,-1,0,0,-1,0,1,1,0,0,3,0,-1,1,-2,0,2,1,-2,3,2,-3,-1,0,0,1,0,1,1,0,0,-2,-1,1,-2,2,-2,-1,1,1,-1,0,0; static double Bd60=1,1,-1,-1,1,-1,1,1,-1,-1,-1,-1,

14、1,-1,1,1,-1,-1,-1,1,3,1,1,1,-1,-1,-1,1,-1,1,-3,1,-3,-1,-1,1,-1,1,-1,1,1,1,1,-1,3,-1,-1,1,-1,-1,1,-1,1,-1,-1,-1,-1,-1,-1,1; / /计算日期和时间/ / year=2021; month=9; hour=0; / /计算时间日期的儒略日历书时/ JDE=jde（year,month,day,hour,min,sec）; / /计算自J2000.0开始的儒略世纪数/ T=（JDE-2451545）/36525; / / L1=218.3164477+481267.8812342

15、1*T-0.0015786*T*T+T*T*T/538841-T*T*T*T/65194000;/月亮的平黄经 D=297.8501921+445267.1114034*T-0.0018819*T*T+T*T*T/545868-T*T*T*T/113065000;/月亮的平均太阳距角 M=357.5291092+35999.0502909*T-0.0001536*T*T+T*T*T/24490000;/太阳的平近点角 M1=134.9633964+477198.8675055*T+0.0087414*T*T+T*T*T/69699-T*T*T*T/14712000;/月亮的平近点角 F=93.

16、2720950+483202.0175233*T-0.0036539*T*T-T*T*T/3526000+T*T*T*T/863310000;/月亮的黄纬参量（由升交点起算的月球平均距离） A1=119.75+131.849*T;/金星的摄动 A2=53.09+479264.290*T;/木星的摄动 A3=313.45+481266.484*T; E=1-0.002516*T-0.0000074*T*T;/计算反映地球轨道偏心率变化 / /计算月球地心黄经周期项；/ SUML=0; for（i=0;i SIN1=Lai*D+Lbi*M+Lci*M1+Ldi*F; SUML=SUML+Sli*0

17、.000001*sin（SIN1*DE）*pow（E,fabs（Lbi）; /计算月球地心黄经/ lamda=L1+SUML+（3958*sin（A1*DE）+1962*sin（L1-F）*DE）+318*sin（A2*DE）/1000000; lamda=fmod（lamda,360）; /计算月球地心黄纬周期项/ SUMB=0; SIN2=Bai*D+Bbi*M+Bci*M1+Bdi*F; SUMB=SUMB+Sbi*0.000001*sin（SIN2*DE）*pow（E,fabs（Lbi）; /计算月球地心黄纬/ beta=SUMB+（-2235*sin（L1*DE） /（modifie

18、d） +382*sin（A3*DE）+175*sin（A1-F）*DE） +175*sin（A1+F）*DE）+127*sin（L1-M1）*DE） -115*sin（L1+M1）*DE）/1000000; /计算月球地心距离周期项/ SUMR=0; COS1=Lai*D+Lbi*M+Lci*M1+Ldi*F; SUMR=SUMR+Sri*0.001*cos（COS1*DE）*pow（E,fabs（Lbi）; /计算月球地心距离/ /Dist=385000.56+SUMR/1000; Dist=385000.56+SUMR; printf（n%d-%d-%d %d:%d:%dn,year,month,day,hour,min,sec）;JDE=%fn,JDE）;T=%fn,T）;L1=%fn,L1）;D=%fn,D）;M=%fn,M）;M1=%fn,M1）;F=%fn,F）;A1=%fn,A1）;A2=%fn,A2）;A3=%fn,A3）;E=%fn,E）;SUML=%fn,SUML）; prin