低轨导航增强卫星广播星历拟合模型研究.pdf

V o l N o 航天器工程S P A C E C R A F TE N G I N E E R I N G第 卷第期 年月低轨导航增强卫星广播星历拟合模型研究杨开伟刘涵宇孙腾达孙秀聪徐明(中国电子科技集团公司第五十四研究所,石家庄 )(北京航空航天大学 宇航学院,北京 )收稿日期:;修回日期:基金项目:中央高校基本科研业务费专项资金资助(YWF L )作者简介:杨开伟,男,硕士,高级工程师,从事卫星导航研究工作.通讯作者:孙秀聪,男,博士,副研究员,从事航天器导航与控制研究工作.摘要传统的广播星历模型(如G P S的 参数模型和 参数模型)是针对中高轨卫星设计的,而中高轨卫星与低轨卫星所受摄动力有所区别,如果直接将传统模型应用到低轨卫星上,拟合精度会降低,因此需要针对低轨卫星所受摄动的特点设计更加合理的广播星历模型.考虑到大气阻力主要改变平面内的轨道要素,文章在 参数模型的基础上增加偏心率变化率这一参数,设计基于轨道根数的 参数模型,并分别对不同轨道高度和倾角的卫星进行位置拟合试验.试验结果表明:在相同仿真场景下,参数模型的位置拟合精度高于 参数模型;同时,对于大倾角卫星的长弧段位置,其拟合精度与 参数模型近似.关键词低轨卫星;广播星历;导航增强;参数模型中图分类号:V 文献标志码:AD O I:/j i s s n R e s e a r c ho nB r o a d c a s tE p h e m e r i sF i t t i n gM o d e l o fL o wO r b i tN a v i g a t i o nA u g m e n t a t i o nS a t e l l i t eYANGK a i w e iL I U H a n y uS UNT e n g d aS UNX i u c o n gXU M i n g(T h e t hR e s e a r c hI n s t i t u t eo fC h i n aE l e c t r o n i c sT e c h n o l o g yG r o u pC o r p o r a t i o n,S h i j i a z h u a n g ,C h i n a)(S c h o o l o fA s t r o n a u t i c s,B e i h a n gU n i v e r s i t y,B e i j i n g ,C h i n a)A b s t r a c t:T r a d i t i o n a lb r o a d c a s t e p h e m e r i sp a r a m e t e rm o d e l s,s u c ha s t h eG P S p a r a m e t e ra n d p a r a m e t e rm o d e l s,a r ed e s i g n e d f o rm e d i u mt oh i g ho r b i t s a t e l l i t e s H o w e v e r,t h e r ea r ed i f f e r e n c e s i nt h ep e r t u r b a t i o nf o r c e sb e t w e e nm e d i u mt oh i g ho r b i t s a t e l l i t e sa n d l o wo r b i t s a t e l l i t e s I f t h e t r a d i t i o n a lm o d e l i sd i r e c t l ya p p l i e dt ol o wo r b i ts a t e l l i t e s,t h ef i t t i n ga c c u r a c yw i l lb er e d u c e d T h e r e f o r e,i t i sn e c e s s a r y t od e s i g nam o r e r e a s o n a b l eb r o a d c a s t e p h e m e r i sm o d e l b a s e do nt h ec h a r a c t e r i s t i c so ft h ep e r t u r b a t i o no f l o wo r b i ts a t e l l i t e s C o n s i d e r i n gt h a ta t m o s p h e r i cd r a gm a i n l yc h a n g e s t h eo r b i t a l e l e m e n t s i nt h eo r b i t a lp l a n e,t h ec h a n g er a t eo fe c c e n t r i c i t yi sa d d e dt ob r o a d c a s t e p h e m e r i sp a r a m e t e r sb a s e do nt h e p a r a m e t e rm o d e l,l e a d i n gt oa p a r a m e t e rm o d e lb a s e do no r b i t a l e l e m e n t s T h e n,p o s i t i o nf i t t i n gt e s t sa r ec o n d u c t e do ns a t e l l i t e sw i t hd i f f e r e n to r b i t a l t i t u d e s a n d i n c l i n a t i o na n g l e s T h e t e s t r e s u l t s s h o wt h a tu n d e r t h e s a m e s i m u l a t i o ns c e n a r i o,t h e f i t t i n ga c c u r a c yo f t h e p a r a m e t e rm o d e l i sh i g h e rt h a nt h a to f t h e p a r a m e t e rm o d e l,a n dt h ep o s i t i o nf i t t i n ga c c u r a c yo f t h e p a r a m e t e rm o d e l f o r t h e l o n ga r co f l a r g e i n c l i n a t i o na n g l es a t e l l i t e s i ss i m i l a r t ot h a to f t h e p a r a m e t e rm o d e l K e yw o r d s:l o wo r b i t s a t e l l i t e;b r o a d c a s t e p h e m e r i s;n a v i g a t i o na u g m e n t a t i o n;p a r a m e t e rm o d e l近年来,全球卫星导航系统(GN S S)迅速发展,在人类社会生活的方方面面发挥着不可或缺的作用,涉及到测绘、交通、救援等诸多方面.随着自动驾驶、物联网、G等新兴产业的发展,社会生产生活对卫星导航服务提出了更高的要求,对时空信息的需求发展为精准、实时、动态和全球 .GN S S的不足之处也逐渐显露出来.GN S S在常规情况下只能提供 m左右的定位精度,而无法满足用户高精度的需求;此外,导航卫星信号到达地面时极其微弱,不但导致其极易受到干扰和欺骗,而且导致其难以穿透物理遮蔽,无法在室内、水下等复杂地形中提供高精度的定位服务,.导航增强技术是解决上述问题、提升导航服务性能的一种可行方案,能利用低轨导航增强卫星产生测距信号并与中高轨GN S S系统协同提供导航服务.传统的导航增强系统采用的是地球静止轨道(G E O)卫星,但G E O卫星轨道资源有限,而且用户接收到的信号弱,容易受到干扰,无法满足用户日益增长的高精度定位需求.低轨卫星则没有上述问题,其轨道高度低,信号的自由空间损耗低,落地信号更强,抗干扰性能更好.此外,由于低轨导航增强卫星的飞行速度较快,相等时间的位置变化较大,从而使得历元间量测信息的相关性下降,能提高系统的可观测度,加快载波相位模糊度的收敛速度.例如,美国的铱下一代(I r i d i u m N E X T)星座提供的卫星授时与定位服务,已经能够脱离G P S单独提供定位服务,因而低轨导航增强卫星可以在G P S受到干扰的情况下提供备份的导航服务.我国的珞珈一号卫星搭载有导航增强载荷,并开展了一系列低轨卫星增强相关试验.高精度确定用户的位置需要导航卫星高精度的位置信息,导航卫星的位置信息由广播星历计算,因此广播星历的高精度计算方法得到了广泛的研究.目前,应用广泛的广播星历参数模型有种,主要是针对中高轨卫星,一种是以G P S为代表的基于轨道根数的广播星历参数模型;一种是以格洛纳斯(G L ONA S S)为代表的基于轨道状态的广播星历参数模 型.G P S广 播 星 历 参 数 模 型 的 精 度 和G L ONA S S相当,但前者的预报能力更强,用户算法更简单;后者的优势是其参数更少,占用导航信号的资源少.卫星设备、通信链路和用户性能等多方面的限制,要求广播星历模型具有参数少、外推能力强、用户算法简单等特点,因此G P S广播星历参数模型应用更加广泛.传统的G P S系统采用的是包括个开普勒轨道参数、个轨道参数变化项、个调和项系数和个参考历元的 参数模型.在原来 参数的基础上,G P S系统后续又发展了 参数模型.参数模型中考虑卫星所受摄动力对卫星轨道半长轴、平均角速度及升交点赤经的影响,相较于 参数模型更加注重卫星轨道的瞬时变化,拟合精度更高.参数模型和 参数模型都是针对中高轨卫星设计的,而低轨卫星与中高轨卫星所受摄动力有所区别,如果将这种参数模型直接应用在低轨卫星上,将无法满足高精度的定位需求,因此需要根据低轨卫星的特点重新设计广播星历参数模型.文献 中利用轨道高度约 k m的低轨卫星在 参数模型的基础上设计了 参数的广播星历拟合算法,增加了轨道偏心率的变化率、升交角线性变化率和个调和改正项振幅共 个参数,其h轨道弧段的均方差小于 m,局部弧段拟合误差小于 m,但参数过多,因此占用的导航信号资源多.文献 中考虑了低轨卫星摄动力的短期变化快,在G L ONA S S广播星历参数模型的基础上设计了基于轨道状态的 参数模型.其研究显示:当低轨卫星轨道高度大于 k m时,用户距离误差(F UR E)精度在 c m以内;当轨道高度降低至 k m时,F UR E精度优于 m.虽然相对于 参数模型来说 参数模型参数少,但是用于确定卫星位置的用户算法外推效率较低,因此限制了其被广泛应用.针对 参数和 参数占用导航信号资源多、用户算法复杂的问题,本文提出了一种兼具拟合精度、计算效率和占用导航信号资源少的模型.考虑到大气阻力主要改变平面内的轨道要素,在 参数模型的基础上仅增加偏心率变化率这一参数,设计了 参数模型,并针对 k m高度的低轨卫星进行了 参数模型最小二乘拟合,检验了这种星历模型的有效性,验证了其拟合精度要优于 参数模型.本文提出的 参数模型对大倾角(不小于 )卫星的长时弧段(h)的位置拟合精度与 参数模型近似,但所需参数量大大减少,可为低轨导航增强卫星的广播星历提供一种兼顾精度与计算资源的模型.第期杨开伟 等:低轨导航增强卫星广播星历拟合模型研究 参数模型 广播星历参数 参数模型采用轨道根数加摄动参数的形式表示,包括个开普勒轨道参数、个轨道根数变化项、个调和项系数和个参考历元.参数模型在 参数模型的基础上添加偏心率变化率这一参数.表给出了参考历元to e时刻 参数广播星历的符号说明.表 参数广播星历参数定义T a b l e p a r a m e t e rb r o a d c a s t e p h e m e r i sp a r a m e t e rd e f i n i t i o n参数符号参数意义t星历参考历元A参考历元半长轴A相对于半长轴参考值Ar e f的偏差A半长轴的变化率e参考历元t的偏心率i参考时刻的轨道倾角参考时刻的升交点赤经近地点幅角M参考时刻的平近点角n相对于参考时刻卫星平均角速度n的偏差n卫星平均运动角速度的变化率i轨道倾角i的变化率相对于参考时刻升交点赤经变化率r e f的偏差cu s纬度幅角的正弦调和改正数的振幅cu c纬度幅角的余弦调和改正数的振幅cr s轨道向径的正弦调和改正数的振幅cr c轨道向径的余弦调和改正数的振幅ci s轨道倾角的正弦调和改正数的振幅ci c轨道倾角的余弦调和改正数的振幅e偏心率e的变化率 参数广播星历用户算法用户算法是指用户使用广播星历计算卫星位置的公式.在短时间内,使用广播星历可以较精确地描述卫星的位置.在之后的一段时间里,广播星历仍能描述卫星的位置,只是误差会随着外推时间变长而增大.对于一般轨道,参数广播星历用户算法如下.)计算观测历元t到参考历元t的时间差tktktt()式中:k为第k时刻.)计算观测历元半长轴AAAr e fAAAA tk()计算卫星平均角速度n/Annntknnn()式中:为地心引力常数;n为观测历元相比参考历元的卫星平均角速度变化量;n为观测历元卫星平均角速度.)计算观测瞬间的平近点角MkMkMn tk()计算偏心率eeee tk()迭代计算偏近点角EkEkMkes i nEk()计算真近点角VkVka r c t a nes i nEkc o sEke()计算摄动改正前的纬度辐角kkVk()计算摄动改正后的纬度辐角uk、卫星矢量半径rk和轨道倾角ik的摄动改正项 uk,rk,ik ukcu ss i n(k)cu cc o s(k)rkcr ss i n(k)cr cc o s(k)ikci ss i n(k)ci cc o s(k)()计算经过摄动改正的uk,rk,ikukk ukrkA(ec o sEk)rkikii tk ik()计算卫星在轨道平面上的位置(xk,yk)xkrkc o sukykrks i nuk()计算观测时刻的升交点经度(格林威治子午圈到卫星轨道升交点)r e fket(e)tk()式中:为升交点赤经变化率;k为第k时刻的升交点赤经;e为地球自转角速率.)计算卫星在地心固连坐标系中的位置(X,Y,Z)Xxkc o skykc o siks i nkYxks i nkykc o sikc o skZyks i nik()航天器工程 卷当直接利用常规的轨道根数方法对轨道倾角近似为 的卫星进行广播星历参数拟合时,会导致法方程系数矩阵成为病态矩阵,使方程组的解误差变大,进而造成拟合精度降低.因此,为了解决小轨道倾角造成的卫星广播星历拟合精度降低甚至拟合失败的问题,本文采用轨道旋转法,先将地心惯性坐标系绕X轴(由地心指向春分点)旋转一个角度,在新坐标系计算广播星历,在用户预测卫星位置时将惯性坐标系旋转相反的角度,得到位置矢量在原惯性坐标系下的表达,进一步得到地心固连坐标系下的位置矢量.对于轨道倾角近似为 的情况,坐标转换的过程见图.图轨道倾角为 时坐标转换过程F i g C o o r d i n a t ec o n v e r s i o np r o c e s sw h e ni n c l i n a t i o na n g l e i s 参数广播星历拟合算法 参数广播星历拟合算法流程如图所示.首先,输入被拟合卫星一段时间的精密位置和初始时刻的精密速度,设置拟合时长,并对小倾角情况进行坐标系旋转.然后,设计 参数的迭代初值,除轨道六要素对应的参数外其他参数置为零.最后,利用最小二乘拟合算法迭代计算 参数,直至迭代改正量足够小时停止.其中,当倾角不大于时,认为属于小倾角情况,需要进行坐标系旋转,本文取为.通过判断改正量的模是否足够小(可以取阈值为 )或已经达到最大迭代次数(可以取 次)作为迭代终止条件.下面对拟合算法进行详细介绍.图广播星历拟合算法流程F i g B r o a d c a s t e p h e m e r i s f i t t i n ga l g o r i t h mf l o w 最小二乘拟合算法广播星历参数的设计是实现高精度拟合的前提条件,而拟合算法的设计则是拟合精度和稳定性的重要保证.拟合算法的核心内容是参数估计,经典的参数估计法有最小二乘估计、极大似然估计、极大验后估计、贝叶斯估计等.其中,经典的最小二乘估计法是比较简单且容易理解的估计方法,在参数模型合理、观测值误差较小的情况下,采用最小二乘拟合算法就可以保证较高的拟合精度.批处理最小二乘拟合算法的基本原理为是结合动力学模型与观测模型对状态量进行迭代平差,使得给定时间段内测量估计值与实际测量值之差的平第期杨开伟 等:低轨导航增强卫星广播星历拟合模型研究方和最小.对于一个非线性系统,有zh(x)()式中:z为低轨卫星的精确位置信息;x为 参数广播星历;h(x)为通过用户算法计算出的低轨卫星的位置信息;为精确位置和计算位置的误差.最小二乘拟合算法是使残差的平方和最小,即让评价函数J最小.J(zh(x)T(zh(x)()由于最小二乘拟合算法只能用于线性系统,因此需要对系统线性化,在x处线性化后的残差可以表示为zh(x)zh(x)hx(xx)zHx()式中:测量模型在状态量初始点处的雅各比矩阵Hhxxx.可以使用线性最小二乘拟合算法求解式(),得到x(HTH)(HTz)()进而,原问题的最小二乘解为xxx.在线性化过程中舍弃了系统二阶及以上的信息,所以不可避免地引入误差,造成精度损失,即所得到的估计值并不是原问题的最优估计值.因此,还需要以所得估计值作为新的初值,重新迭代计算,直到算法收敛,即可得到状态量的最优估计值.文献 中 定量研究了 调和项 对A,e,M,n,n这个参数的偏导数精确值的影响的数量级为 ,而其他参数不受调和项的影响,因此可以忽略调和项对偏导数的影响,从而简化偏导数的计算.在文献 中还详细给出了各时刻位置对参考时刻 参数广播星历各参数的偏导数.本文在此基础上进一步给出各时刻位置对偏心率变化率e的偏导数.rkerkeeerketk()式中:rk为第k时刻的卫星位置;rk/e由文献 给出.数值稳定修正数值稳定修正是为了确保拟合算法能得到正确的结果,下面对修正方法进行具体介绍.()偏心率.在批处理最小二乘拟合算法迭代过程中,可能会出现偏心率大于的情况,导致后续运算出错(如e),因此,当某次迭代改正后的偏心率大于时,应将其置为,可以取 .()升交点赤经、近地点幅角、平近点角.为保证数值稳定性,需要在每次迭代施加改正量之后对升交点赤经、近地点幅角、平近点角进行归一化,也就是让其落在,之后再继续进行下一次迭代.例如,当Mk过大时,可能会因为数值问题无法通过迭代求出Ek.()解线性方程组.最小二乘拟合算法问题需要求解线性方程组(HTH)xHTz,但HTH条件数很小,甚至接近病态,直接求逆或者直接使用高斯消元法求解误差较大.因此,本文使用下三角上三角(L U)分解法求解线性方程组.L U分解法可以将个矩阵分解为个下三角矩阵和个上三角矩阵的乘积.如果使用杜利特尔(D o o l i t t l e)分解,得到的下三角矩阵为单位下三角矩阵;如果使用克劳特(C r o u t)分解,得到的上三角矩阵为单位上三角矩阵.首先,对(HTH)xHTz中的HTH进行L U分解,即L UHTH,则(HTH)xHTz则可以写为L UxHTz,令yUx,可以先解线性方程组L yHTz,得到y,再解线性方程组yUx,得到x,从而相比于直接解线性方程组能够降低系数矩阵的病态程度,误差更小.仿真结果及分析本文使用高精度轨道预报器预报的卫星位置作为精密星历,高精度轨道预报器中非球形引力摄动考虑 阶 次 全 球 超 高 阶 地 球 重 力 场 模 型(E GM ),大气阻力模型使用质谱仪非相关散射模型(N R LM S I S E ),同时考虑日月三体引力;参考时刻为 T :(UT C),参考时刻的偏心率取,升交点赤经、近地点幅角和平近点角均取.下面对不同轨道高度(k m,k m,k m,k m)及不同轨道倾角(,)的轨道进行 参数广播星历最小二乘拟合.设置拟合初值如下.参考时刻半长轴取为参考时刻卫星的半长轴,参考时刻升交点赤经变化率取为,参数中的参考时刻、偏心率、升交点赤经、近地点幅角和平近点角按照卫星参考时刻的实际数据设置,其余参数设为.卫星位置采样时间间隔为m i n,即每分钟获取组位置,拟合时长分别考虑 m i n,m i n,航天器工程 卷 m i n,h,h,h的情况.绘制上述仿真场景下的位置均方根误差(R M S E),计算公式为RMS E(xyz)/()式中:x,y,z为地心固连坐标系下的位置误差.表中给出了 参数和 参数模型拟合的RM S E曲线对各个时刻的均方根,图中给出了同一轨道倾角、同一拟合时长时不同轨道高度对拟合效果的影响.从表和图可以看出:对于任何拟合时长,拟合RM S E随轨道高度的升高而减小,原因是轨道高度越高,同等时间下轨道的非线性越弱,最小二乘法更能发挥优势;同时,随着轨道高度的升高及非球形引力摄动和大气阻力摄动的下降,轨道更接近二体轨道,因此拟合误差会有所下降.参数模型拟合时长为 m i n时精度达到厘米级,m i n时精度达到分米级,m i n时精度达到米级.随着拟合时间的增长,RM S E有所增大,这是因为轨道根数预报是基于二体模型的,即使广播星历增加了摄动修正项,但其只有个系数,远远不足以描述全部的摄动力,因此也只能解决短时间内的拟合问题,随着时间增长,误差仍然会越来越大.总体上来说,参数模型的拟合精度要高于 参数模型.从表和图可以看出:轨道倾角为 时同样可以得到与轨道倾角为 时相同的结论,即对于任何拟合时长,拟合RM S E随轨道高度升高而减小,整体上 参数模型的精度要高于 参数模型.对比轨道倾角为 和 的种场景,可以发现:当轨道高度和拟合时长相等时,轨道的拟合误差略大于 轨道的拟合误差.从表和图可以看出:轨道倾角为 时同样可以得到与轨道倾角为 和 时相同的结论.轨道倾角 时的拟合误差大于轨道倾角 时的拟合误差,与 倾角时近似.从表和图可以看出:轨道倾角为 时同样可以得到与轨道倾角为、和 时相同的结论.总的来看,参数模型对不同轨道高度(k m)、不同轨道倾角的卫星在不同拟合弧段时长下的拟合均具有很好的适用性;同时,参数模型的精度要高于 参数模型.表轨道倾角为 时不同轨道高度拟合结果的RM S ET a b l eRM S E so f f i t t i n gr e s u l t s f o rd i f f e r e n t o r b i t a l t i t u d e sa t a n i n c l i n a t i o na n g l eo f m轨道高度/k m m i n m i n m i n m i n m i n m i n 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 表轨道倾角为 时不同轨道高度拟合结果的RM S ET a b l eRM S E so f f i t t i n gr e s u l t s f o rd i f f e r e n t o r b i t a l t i t u d e sa t a n i n c l i n a t i o na n g l eo f m轨道高度/k m m i n m i n m i n m i n m i n m i n 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 表轨道倾角为 时不同轨道高度拟合结果的RM S ET a b l eRM S E so f f i t t i n gr e s u l t s f o rd i f f e r e n t o r b i t a l t i t u d e sa t a n i n c l i n a t i o na n g l eo f m轨道高度/k m m i n m i n m i n m i n m i n m i n 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 第期杨开伟 等:低轨导航增强卫星广播星历拟合模型研究表轨道倾角为 时不同轨道高度拟合结果RM S ET a b l eRM S E so f f i t t i n gr e s u l t s f o rd i f f e r e n t o r b i t a l t i t u d e sa t a n i n c l i n a t i o na n g l eo f m轨道高度/k m m i n m i n m i n m i n m i n m i n 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 参数 图轨道倾角、拟合时长 m i n时不同轨道高度对 参数和 参数模型的拟合效果F i g F i t t i n gr e s u l t so fd i f f e r e n to r b i t a l t i t u d e so n p a r a m e t e ra n d p a r a m e t e rm o d e l sw h e ni n c l i n a t i o na n g l e i s a n df i t t i n gt i m e i s m i n u t e s图轨道倾角 、拟合时长 m i n时不同轨道高度对 参数和 参数模型的拟合效果F i g F i t t i n gr e s u l t so fd i f f e r e n to r b i t a l t i t u d e so n p a r a m e t e ra n d p a r a m e t e rm o d e l sw h e ni n c l i n a t i o na n g l e i s a n df i t t i n gt i m e i s m i n u t e s图轨道倾角 、拟合时长 m i n时不同轨道高度对 参数和 参数模型的拟合效果F i g F i t t i n gr e s u l t so fd i f f e r e n to r b i t a l t i t u d e so n p a r a m e t e ra n d p a r a m e t e rm o d e l sw h e ni n c l i n a t i o na n g l e i s a n df i t t i n gt i m e i s m i n u t e s图轨道倾角 、拟合时长 m i n时不同轨道高度对 参数和 参数模型的拟合效果F i g F i t t i n gr e s u l t so fd i f f e r e n to r b i t a l t i t u d eo n p a r a m e t e ra n d p a r a m e t e rm o d e l sw h e ni n c l i n a t i o na n g l e i s a n df i t t i n gt i m e i s m i n u t e s航天器工程 卷进一步对比 参数模型和 参数模型的长弧段(h)位置拟合精度,结果如表所示.从表可以看出:当轨道倾角不小于 时,在 种不同轨道高度和倾角的情况中,有种情况是 参数模型精度更高,另外种情况是 参数模型精度更高.这证明了增加轨道倾角变化率这一参数对长弧段拟合精度的提升是有显著帮助的,使得 参数模型对于大倾角卫星的长弧段拟合达到了与 参数模型近似的精度.同时,参数模型在有些情况下拟合精度低于 参数模型,这是因为 参数模型是由 参数模型发展而来的,其没有考虑 参数模型相较 参数模型增加的轨道高度变化率和角速度变化率.表不同轨道高度、不同倾角的 h拟合结果RM S ET a b l eRM S E so f f i t t i n gr e s u l t s f o rd i f f e r e n t o r b i t a l t i t u d e sa n d i n c l i n a t i o na n g l e sw h e nf i t t i n g t i m e i sh o u r sm倾角/()高度/k m 参数 参数 参数 本文在G P S的 参数模型基础上增加了偏心率变化率这一参数,仿真结果证明这一参数对于提高低轨导航增强卫星的位置拟合精度具有重要作用.如何在本文提出的 参数模型基础上,在不增加过多参数前提下进一步提高位置拟合精度,值得进一步研究.表为此提供了重要方向,即:对于小倾角()卫星的长弧段拟合,参数模型的位置拟合精度远高于 参数模型;而对于大倾角(不小于 )卫星,二者精度则近似.因此,导致 参数模型对小倾角卫星的位置拟合精度较高的主要参数十分值得研究,进而可以将这些起主要作用的参数与 参数模型结合,以提高后者对于小倾角卫星的位置拟合精度.结论本文针对轨道高度为 k m的低轨卫星,在G P S的 参数模型基础上仅增加偏心率变化率这一参数,设计了 参数模型,并与 参数模型的精度进行了比较,试验分析结论如下.()对于轨道高度为 k m的低轨卫星,采用 参数模型拟合精度高.()随着低轨卫星轨道高度的降低,参数模型的拟合精度降低,这是因为轨道高度越低,大气阻力的影响越显著,导致拟合精度降低.()总体来说,在相同轨道高度、轨道倾角和拟合时长的条件下,参数模型的拟合精度要高于G P S的 参数模型.此外,本文通过仿真验证了 参数模型对于大轨道倾角卫星的长弧段位置拟合精度与 参数模型近似,但所需参数量大大减少,因此可以提高计算效率和减少存储资源的需求.因此,本文的研究为低轨导航增强卫星的广播星历提供了一种兼顾精度高与资源需求少的模型.参考文献(R e f e r e n c e s)张小红,马福建低轨导航增强G N S S发展综述J测绘学报,():Z HANGX i a o h o n g,MAF u j i a n O v e r v i e wo fGN S Sd e v e l o p m e n t e n h a n c e db yl o we a r t ho r b i tn a v i g a t i o nJA c t aG e o d a e t i c ae tC a r t o g r a p h i c aS i n i c a,():(i nC h i n e s e)徐伟证基于低轨导航增强的全球实时精准定位研究D南京:南京航空航天大学,X U W e i z h e n g R e s e a r c ho ng l o b a l r e a l t i m e a c c u r a t ep o s i t i o n i n gb a s e do ne n h a n c e dl o we a r t ho r b i tn a v i g a t i o nD N a n j i n g:N a n j i n g U n i v e r s i t yo fA e r o n a u t i c sa n dA s t r o n a u t i c s,(i nC h i n e s e)梅林,鲍尧 G N S S干扰分析及抗干扰技术综述J中国无线电,():ME IL i n,B AO Y a o GN S Si n t e r f e r e n c ea n a l y s i sa n da n t i i n t e r f e r e n c et e c h n o l o g yr e v i e wJ C h i n aR a d i o,():(i nC h i n e s e)田润,崔志颖,张爽娜,等基于低轨通信星座的导航增强技术发 展 概 述 J导 航 定 位 与 授 时,():T I AN R u n,C U IZ h i y i n g,Z HAN G S h u a n g n a,e ta l O v e r v i e wo fn a v i g a t i o ne n h a n c e m e n tt e c h n o l o g yd e v e l o p m e n tb a s e do nl o we a r t ho r b i tc o mm u n i c a t i o nc o n 第期杨开伟 等:低轨导航增强卫星广播星历拟合模型研究s t e l l a t i o n sJN a v i g a t i o n P o s i t i o n i n g a n d T i m i n g,():(i nC h i n e s e)王磊,李德仁,陈锐志,等低轨卫星导航增强技术 机遇与挑战J中国工程科学,():WAN GL e i,L ID er e n,CHE NR u i z h i,e t a l E n h a n c e dl o we a r t ho r b i ts a t e l l i t en a v i g a t i o nt e c h n o l o g yo p p o r t u n i t i e sa n dc h a l l e n g e sJ S t r a t e g i cS t u d yo fC A E,():(i nC h i n e s e)方善传,杜兰,周佩元,等低轨导航增强卫星的轨道状态型星历参数设计J测绘学报,():F AN GS h a n c h u a n,DUL a n,Z HOUP e i y u a n,e t a l D e s i g no fo r b i ts t a t ee p h e m e r i sp a r a m e t e r sf o rl o wo r b i tn a v i g a t i o ne n h a n c e m e n ts a t e l l i t e sJ A c t aG e o d a e t i c ae tC a r t o g r a p h i c aS i n i c a,():(i nC h i n e s e)许钡榛基于北斗系统的珞珈一号低轨卫星近实时精密定轨方法研究D武汉:武汉大学,XUB a z h e n R e s e a r c ho nn e a rr e a lt i m ep r e c i s i o no r b i td e t e r m i n a t i o nm e t h o df o rL u o j i al o we a r t ho r b i ts a t e l l i t eb a s e do nt h eB e i d o us y s t e mD Wu h a n:W u h a nU n i v e r s i t y,(i nC h i n e s e)方善传L E O导航增强卫星广播星历模型设计D郑州:解放军信息工程大学,F AN GS h a n c h u a n D e s i g no fL E On a v i g a t i o ne n h a n c e ds a