现代化的天文定位——计算机代替天文定位中的查表的论文.docVIP

　　现代化的天文定位——计算机代替天文定位中的查表的论文 摘要：利用计算机完成天文定位中的太阳赤经、赤纬、太阳格林时角、春分点格林时角、测者地方时角的计算。 关键词：天文定位太阳赤经太阳赤纬太阳格林时角 　　近年来，电子导航系统不断发展，特别是gps导航仪的普及，使得船舶驾驶人员不重视天文定位。但是，天文导航有其独特的优越性:设备简单、可靠，观测的目标是自然天体而不受人为控制，不发射任何声、光和电波而具有隐蔽性等。因此，船舶驾驶员适任证书评估考试中仍然考核天文定位和天测罗经差，这些内容包含在天文导航中。但是，天文计算内容繁琐、劳动量大、容易出现错误。为了解决这些问题，用计算机实现天文定位就成为必然。以前的天文定位是用六分仪观测数据，然后查表、计算得出船位，这样费时费力。而现在我们可以通过射电六分仪准确快速的观测天体高度，然后通过计算机计算得出船位。整个过程在一分钟之内完成，可以说是实现了时时定位。 射电六分仪它工作在微波波段，自动跟踪太阳、月球或人造卫星，把天体（或人造卫星）地平坐标数据连同时间信息输入电子计算机，处理后给出船只或飞机所在地的地理位置。射电六分仪观测不受云雾甚至暴风雪的影响，因此可以全天候工作。假如在战争期间为了不发射信号，我们可以把射电六分仪换成普通六分仪，虽然操作需要人工，但也并不复杂。例如：应用太阳中天高度求纬度，在太阳中天时间，前后只需一分钟就可以测得太阳的高度，把数值输入计算机后马上得出测者的纬度，前后所需时间最多两分钟。. 传统天文定位费时费力的关键是查表和计算。计算可以轻松的通过计算机来实现，天文定位要实现计算机化，查表求天文数据就成为关键。下面主要介绍如何利用计算机求取天文数据来代替查表。 一、时间计算 1.求出儒略日jd 即jd=367*year-int（7*（year+int（（month+9）/12））/4）+int（275*month/9）+date+1721013.5+gmt/24 其中：year——年数（如2006）； month——月数（如10）； 　　date——日期（如16）； 　　gmt——世界时（如21）； 　　int——取整数。 2.从j2000.0起算的儒略世纪数t t=（jd-2451545.0）/36525 二、太阳平均轨道根数 1.太阳赤道偏心率 es=0.016708617-0t-0.0000001236*t2 2.太阳平均平近点角 ms0=357°31′39#733;.804+（99r+359°03′01#733;.224）*t-0#733;.577*t2-0#733;.012*t3 3.平黄赤交角ε0 ε0=23°26′21#733;.448-46#733;.815*t-0#733;.00059*t2+0#733;.001813*t3 4.白道升交点的平黄经ωn ωn=125°02′40#733;.28-(5r+134°08′10#733;.539)*t+7#733;.455*t2+0#733;.008*t3 5.对于当天平春分点的太阳几何平黄经ls0 ls0=280°27′57#733;.85+（100r+0°46′11#733;.27）*t+1#733;.089*t2 其中1r=360°=1296000#733; 三、其他所需的天体轨道根数 1.日月距角d d=277°51′01#733;.307+（1236r+307°06′41#733;.328）*t-6#733;.891*t2+0#733;.011*t3 2.金星的平近点角mv mv=50°24′57#733;.937+(162r+197°48′38#733;.915)*t+4#733;.98*t2+0#733;.019*t3 3.火星的平近点角mm mm=19°22′22#733;.947+(53r+59°51′19#733;.55)*t+0#733;.6331*t2+0#733;.0009*t3 四、太阳轨道根数的改正 1.太阳平近点角的摄动改正δms δms=6#733;.4*sin(251°.4+20°.2*t)+1#733;.9*sin(207°.5+150°.3*t) 2.黄赤交角章动改正δε δε=9#733;.2*cosωn 3.太阳黄经章动改正δφ δφ=-17#733;.2*sinωn 4.中心差v v=(2es-(es)3/4)sinms+(5(es)2/4-11(es)4/24)sin2ms+13(es)3sin(3ms)/12+103(es)4sin(4ms)/96 5.太阳平黄经的主要摄动项 （1）长周期项 δl