GLONASS广播星历轨道计算方法的研究.pdfVIP

GLONASS广播星历轨道计算方法研究 张永军徐绍铨张小红李英冰 (武汉测绘科技大学地测学院，武汉市珞瑜路129号，430079) 【捅要1GLONASS系统建成后，它与GPS系统的观测数据联合处理所提供的48早鹰组 合定位系统，已经引起了曲’多现彳{GPS川户的广泛关注。本文庄简单介}“ 种 GLONASS系统的基础上，阐述了常川的J。善罱甩历轨道积分模州，”提出了 新的积分方法。该方法贝订编科简单，运算速度快等优点。 一、 引言 前苏联的全球导航p星系统(Global Satellite Navigation 1996年初建成，f=投入运行。该系统由24倾1J吊细成，均匀分布在夹角为1200的二个轨道 面I～．其轨道倾角为64．80。GLONASS系统建成并投入整体运行以来．该系统尤其是它与 GPS系统的组合定位系统，受剑了国内外彳『天专家和学者lfJftqJ“泛关注。 GLONASS系统同GPS系统扫：设计原理上很相似，二者都I竹在L波段上，分别发射 C／A码平¨P码。原始的码观测值年¨载波相位观测值可转换为卫星剑接收机的距离，『J星的 位置则通过广j罱星历获取，然屙通过距离交会法确定接收机的位置。但由于GLONASS系 统所采_HJ的时间基准利空间基准与GPS系统不同，其广播星历的数据格式也不同丁．GPS 系统，而且GLONASS系统采川频分多址技术，即每颗■吊都使川不同的i‘作频率发射信 号． 这就使得现有GPS数据处理软什无法对GLONASS系统的观测数据进行处理。解决 这一问题的常川方法是利川GLONASS系统的广播星历进行轨道积分，计算小GLONASS 中，这样GLONASS码观测数据就可以直接与GPS码观测数据进行联合处理。本文在阐述 积分时间较短时的轨道积分模耻的基础上．提山采川编程简单，运算速度快的定步长四阶 龙格一库塔方法进行轨道积分．并分别对积分区间雨l积分步K的不同情况编群进行了r计算 平¨It较，撮后讨论了适合丁实时定位和后处理定位的GLONASSp星坐标的计筇方法。 二、 GLONASS广播星历 GLONASS广播星历每隔、r小时更新‘。‘次，屉历中包含着与GPSJ“插屉历相似的信息． 但其数锯格式及星历内容并不相同。GLONASS甲尿的位置需要依据其广橘星历通过轨道 积分进行计算，另外，与GPSp星星历每小时更新一次并：可向左右外推两小时所不同的是． GLONASS TJ．昂星历只有、15,J,Bt的有效性。GLONASSJ“}播屉历中与轨道确定有芙的参数 士要有参考时刻t、、相对论效戍政止、n，a钟筹、p艰忙置、■星速度、太日1和月亮摄动 加速度之币I等，其中t。为J“播晕历的参考时间基准，其余所有相关信息均以tb为参考值。 GLONASS系统的稳定性还有待进一步提高，其广播星历的可靠性会有所变化，偶尔 可能会有某颗卫是的导航信息部分或全部置零，致使接收机虽然观测剑了该星的码平¨载波 相位数据，却冈得不到甲星星历而无法参加计算。一般而青，GLONASSJ“播星历中甲昂 的位置年Il速度的误筹在沿轨道方向分别为20m和0．05em／s，与轨道Ⅱ三交方向为10m和 0．1cm／s，而径向则为5m和O．3em／s左右。 三、 轨道模型的数值积分 153 p星在空问的延行轨道主要是由作_L}j在其上的巫力所决定的，璁星的重力位可表示为 r．_坐+u (I) 其中GM为常量，U为重力场的离心力位部分，r为从地球儿何¨0到1]肇质心的距离。 将U埘球谐函数展开井忽略地球重力场中非带谐部分和太Rf辐射成’譬次要冈素的影响，可 以得剑如r的二阶{。数表达式： 孚一警x+五!丁；Ma z,z[t一爿- 幢， 华=一等，+五等y(1_等H基 ’