七政四余三王的算法

《七政四余三王的算法》由会员分享，可在线阅读，更多相关《七政四余三王的算法（9页珍藏版）》请在装配图网上搜索。

1、七政四余三王的算法目录1,儒略数d的计算2，行星轨道要素。3, 根据儒略历计算太阳的平黄经。4, 根据儒略历年计算月球的平黄经。5, 天球坐标转换。6, 恒星时和地方时间7 一个精确计算的例子正文丄儒略数d的计算YYYY年MM月DD日TT时d = 367XY - (7X(Y + (M+9)/12)/4 + (275XM)/9 + D+TT/24 - 7305302,行星轨道要素：21)参数介绍a (半长轴)e (偏心率)I （轨道面和黄道交角），01180度，若I大于90度，则是逆轨天体, 如哈雷彗星。N （升交点黄经）（大写Q）从春分点到升交点的角度。W （升交点到近日点的角度一近日点黄经）

2、3）升交点到近日点角度。M （平近点角）L （平均黄经度）二M +w+ NE （偏近点角） M 二 E - e X sin（E）V（真近点角）r （向径）22）太阳系主要天体的轨道要素太阳轨道要素：w = 282.9404_deg + 4.70935E-5_deg * da =1.000000e = 0.016709 - 1.151E-9 * dM = 356.0470_deg + 0.9856002585_deg * d太阴轨道要素：N = 125.1228_deg - 0.0529538083_deg * di= 5.1454_degw = 318.0634_deg + 0.1643573

3、223_deg * d (a= 60.2666e = 0.054900水星轨道要素：N = 48.3313_deg + 3.24587E-5_deg * di = 7.0047_deg + 5.00E-8_deg * d w = 29.1241_deg + 1.01444E-5_deg * d a =0.387098e = 0.205635 + 5.59E-10M = 168.6562_deg + 4.0923344368_deg * d金星轨道要素：N = 76.6799_deg + 2.46590E-5_deg * di = 3.3946_deg + 2.75E-8_deg * dw =

4、54.8910_deg + 1.38374E-5_deg * da = 0.723330e = 0.006773 - 1.302E-9 * dM = 48.0052_deg + 1.6021302244_deg * d火星轨道要素：49.5574_deg + 2.11081E-5_deg * d i = 1.8497_deg - 1.78E-8_deg * dw = 286.5016_deg + 2.92961E-5_deg * da = 1.523688e = 0.093405 + 2.516E-9 * dM = 18.6021_deg + 0.5240207766_deg * d木星轨道要

5、素：N = 100.4542_deg + 2.76854E-5_deg * di = 1.3030_deg - 1.557E-7_deg * dw = 273.8777_deg + 1.64505E-5_deg * da = 5.20256e = 0.048498 + 4.469E-9 * dM = 19.8950_deg + 0.0830853001_deg * d土星轨道要素：N = 113.6634_deg + 2.38980E-5_deg * di = 2.4886_deg - 1.081E-7_deg * dw = 339.3939_deg + 2.97661E-5_deg * da

6、 = 9.55475e = 0.055546 - 9.499E-9 * dM = 316.9670_deg + 0.0334442282_deg * d3，根据公元纪年计算太阳的平黄经及黄赤交角即黄经的变化量与儒略数d的关系PI = 3.14159265358979323846lw = 282.9404_deg（度）+ 4.70935E-5_deg X da =1.000000e = 0.016709 - 1.151E-9X dlM = 356.0470_deg + 0.9856002585_deg X dl黄赤交角 oblecl = 23.4393_deg - 3.563E-7_deg X

7、dlL = M +w+ NN （升交点黄经）（大写Q）从春分点到升交点的角度因为对太阳来说，升 交点就是春分点，故N=0L=M+w2005年6月17日23： 00UT (UT：格林威治时间)d=1995.35, 取整，d=1995w = 282.9404 + 4.70935E-5 x d = 283.03437M = 356.0470 + 0.9856002585 x d =2322.6686即：M = 162.6686太阳的平黄经L=M+ w=445.7029度L=85.7029岁差：3.82394E-5 X ( 365.2422 X (公元纪年-2000.0 ) - d )=0.2 L=8

8、5.7029+0.2=85.9黄赤交角 oblecl = 23.4393 - 3.563E-7 X d=23.4393 度4,根据公元纪年计算月球的平黄经及轨道角度。月球的计算较复杂，因其轨道与黄道有一个交角（5 . 1 4 1 5度）,精确计算还要考虑太阳木星 金星的摄动，其黄经黄纬的校正项多达100多项.限于篇幅此处就不介绍了，只讲求近似值。轨道要素：N = 125.1228 - 0.0529538083 X di =5.1454度w = 318.0634 + 0.1643573223 Xa= 60.2666(Mean distance)（升交点经度）（轨道倾角）d （Arg. of pe

9、rigee（近地点角）e =0.054900 (Eccentricity)M = 115.3654 + 13.0649929509 X d (平近点角)l试计算1990年4月19日0:00的月球黄经4-1）,求得儒略数为：d=-3543,代入轨道要素得:312.7381i =5.1454w = -264.2546 修正为95.7454a = 60.2666e = 0.054900M = -46173.9046,用 360X 129 加上修正为 266.095442),根据开普勒第二(面积)定理推得E=M+exsi n(E)用迭代法求解此超越方程先求 E0, PI=3.1415926E0-M +

10、 (180/pi) x e x sin(M) x (1 + e x cos(M)度E1 = E0 - (E0 -(180/pi) X e X sin(E0) - M) / (1 -e X cos(E0)En+1=En-(En-(180/pi) x e x sin( En)-M)/(1- e x cos (En)直到En+1和En的差非常小为止代入数值，得到:E0u262.9689度迭代后解得E=262.973543)，求轨道面上的直角坐标x = r X cos(v) = a X (cos(E) - e)y = r X sin(v) = a X sqrt(1 - eXe) X sin(E) 得到

11、 x = -10.68095y = -59.723774-4)，求真近点角v及向径rr = sqrt( xXx + yXy ) = 60.67134 地球半径V= 180/3.1415926*atan(-59.72377/-10.68095)-180=-100.1395 加上360度的倍数,使之大于0,得到V=259.8605_deg4-5) ,化为黄道直角坐标:公式: xeclip = r * ( cos(N) * cos(v+w) - sin(N) * sin(v+w) * cos(i) ) yeclip = r * ( sin(N) * cos(v+w) + cos(N) * sin(v

12、+w) * cos(i) ) zeclip = r * sin(v+w) * sin(i) 代入数值,注意度要化弧度,N= 312.7381r=60.67134 ,v=259.8605*PI/180 w=95.7454*PI/180I=5.1454*PI/180得到日心三维直角坐标：xeclip = +37.65311 yeclip = -47.57180 zeclip = -0.4168746),化为黄道坐标系6.1公式:向径:r =sqr t(x *x+y *y+z *z)黄经 L: RA=atan(y/x)黄纬 Lat : Decl=atan(z/sqrt(x*x+y*y)代入数值,注意

13、求反正切时要用第4中判断角的象限:方法得到:ra = 308.3616_deglat = -0.3937_degr = 60.67137，恒星时与时角、地方时SIDTIME = GMST0 + UT + LON/15 地方恒星时=格林尼治00：00恒星时+世界时(格林尼治时间)+经度(东经为正 西经为负)/15GMSTO (格林尼治 0 点恒星时)二(L + 180_deg ) / 15 = L/15 + 12hL：太阳平黄经，L=26.8388 :经度为东经15度的地方恒星时为：GMST0 = 26.8388_deg/15 + 12h = 13.78925 hoursSIDTIME = GM

14、ST0 + UT + LON/15 = 13.78925h + 0 + 15_deg/15 = 14.78925hoursSIDTIME = 14h 47m 21.3s计算地平坐标需要知道时角HA:HA = SIDTIME - RAHA = 14.78925h - 1.77720h = 13.01205h = 195.1808_deg坐标转换，求出七政黄道后可以做出七政静盘，弧角星盘的算法，要用到坐标 变换，公式如下，有心人可以自己研究推算，另有黄道十二宫的详细算法，因 为笔者在其他网站已经发表过，此处就不写了地平坐标系：设天体P点的高度为a,是从天文地平沿通过P点和天文天顶的大圆向天文天 顶

15、（或天底）计量的角距离，向天顶为正（向天底为负），P点和天文天顶（天底）重合，则a =十90。（一 90。）.高度的余数称为天项距，北天极的高度是观测者的天文纬度P点的方位角A,是从南点向西沿着天文地平到通过P点和z点的大圆与天文地平的交点所量度的角距离.&赤纬，h时角，设a为赤经，h=T-a, t恒星时，天文纬度，A是地平方位角B黄纬，入黄经，e黄赤交角。a是地平高度7.1，从地平坐标到赤道坐标cos&*cosh=cos *sin a+s 吨* cosa *cosAcos&*sinh= cosa *sinAsin 6=s ingsin a-coscosa *cosA7. 2从时角和赤纬坐标系

16、到地平坐标系cosa *cosA= sing cos6 *cos h-coss in&cosa *sinA= cos6 *sinhsin a=s ingsin 6+coscos6 *cosh7.3 从黄道坐标系到赤道坐标系cos6cosa=cosgcosAcos&s ina= cosecosgs in入-sinzsinB sin 6=coses in g+s in ecosgs in入7. 4 从赤道坐标系到黄道坐标系cosgcosA=cos6cosacosgs inA=sinesin 6+cosecos6s inasin g=coses in 6-s in ecos&s ina举一个精确计算

17、的例子：昆明市七政四余三王真躔经度102.7纬度25年月日时分北京时间200742812 0200742810 50.8到方 高 度黄经赤 经 赤 纬日 124.96672.23571 37.83792.36503714.12268044月 261.4755 -54.7374241.045115.85572 -25.0010048金 76.4509537.2137878.751235.171086 25.1310498木259.7284 -35.9382259.090517.21336 -22.2646662水140.488872.5296531.81744 2.0008711.1401967

18、4火227.201346.76688 346.78 23.2265 -6.53865875土 59.57946 -20.1115138.26629.40809216.59498244计都 224.7706 53.2332353.275923.58842 -2.66942139 罗睺 44.77057 -53.2332173.275911.588422.669421395 紫炁 218.231559.581560.5750210.0351720.228718286 月孛 38.23147 -59.5816 180.57512.03517 -0.22871829天王 227.674247.6974

19、7347.489223.25243 -5.63336956海王 239.512723.85027321.8535 21.6204 -14.4575283冥王 256.2059 -29.5288266.116217.71795 -23.34818042007年4 月28日 12：00，昆明的十二宫黄道如下宫数 黄道1/7 宫 117.4619172/8 宫 142.7954333/9 宫 169.5201364/10 宫200.9969015/11 宫236.4052266/12 宫269.5243837/1 宫 297.4619178/2 宫 322.7954339/3 宫 349.52013610/4 宫 20.99690111/5 宫 56.40522612/6 宫 89.524383