-
近年来,全球卫星导航系统日趋完善,应用不断拓展,在各领域为用户提供高质量的定位、导航和授时服务[1]。另一方面,现有导航卫星大多部署在中高轨道,给导航服务质量带来不利因素。如卫星运行在较高轨道时,空间几何构型短时间不会显著变化,导致精密单点定位收敛时间较长[2-3];此外,高运行轨道必然引起传输链路的长时延和高损耗,导致接收端信号微弱且容易受到干扰和欺骗,影响服务性能等。
针对上述问题,研究人员提出了部署低轨卫星的解决方案并受到广泛关注[4]。低轨卫星星座的空间几何构型变化快,能够辅助提升定位收敛时间及定位精度[5-6];由于传输距离更短,低轨卫星播发的导航信号在接收端质量更好,能够大幅提升导航信号的精度和覆盖能力。目前,低轨卫星是对现有卫星导航系统的有效补充,是实现导航增强的重要手段。
在低轨卫星的部署过程中,采用何种星座构型是影响其覆盖能力和导航增强效能的重要因素。针对全球覆盖的要求,文献[7-9]归纳分析了低轨卫星轨道设计的主要准则和权衡策略。文献[7]给出了一种几何精度因子(geometric dilution of precision, GDOP)和全球可见星数量均匀分布的星座设计方案;文献[8]将常见低轨卫星星座方案与北斗3号系统(BeiDou-3 system, BDS-3)结合,针对位置精度因子和全球可见星数量进行了分析对比,论证了低轨卫星对BDS-3系统的增强效果;文献[10]分析了“鸿雁”低轨星座对现有导航系统在GDOP和精密单点定位首次收敛时间等方面的增强效果。
上述关于低轨卫星星座设计的研究都是针对全球导航增强的需求,主要采用了基于近圆形轨道的星座构型。另一方面,出于部署难度、成本以及应用场景的考虑,很多时候需要研究如何部署低轨卫星对局部重点区域进行导航增强。针对该问题,文献[11]设计了一种由256颗低轨近圆形轨道卫星组成的Flower星座构型方案。该方案仍然使用近圆形轨道,卫星对轨道下的所有区域有着几乎相同的覆盖和导航增强效能,当对覆盖和导航增强的需求有区域限制时,这种无区别覆盖的方案难以达到最优效率。与近圆形轨道不同,椭圆轨道更适合于实现局部区域的覆盖。文献[12]引入了回归系数,可以定量地计算椭圆卫星轨道。文献[13]给出了一种限制椭圆轨道远地点的星下点与目标点重叠的星座设计方法,验证了椭圆轨道星座对特定区域的覆盖性能优于近圆轨道星座。但是,该方案只保证了最初始设计的那颗卫星的远地点限制条件,不能保证同一轨道面的所有卫星在远地点时都能经过目标点上空,仍然会造成观测时间和空间上的浪费。
在上述研究基础上,本文提出了一种新型的面向区域导航增强的低轨卫星星座方案,和现有方案相比,在低轨卫星数量相同时,能够有效提升覆盖重数、过顶时间、GDOP等影响导航性能的关键指标,具有更好的导航增强效果。
为了更清楚地说明本文讨论的低轨区域导航增强星座的作用与必要性,将其与北斗2号(BDS-2)区域导航星座[14]进行比较。BDS-2采用地球静止轨道(geostationary earth orbit, GEO)卫星、倾斜地球同步轨道(inclined geosynchronous orbit, IGSO)卫星、与中圆地球轨道(medium earth orbit, MEO)卫星实现对亚太地区的覆盖。由于采用中高轨卫星,存在卫星运动动态小,精密单点定位过程中载波整周解模糊过程收敛时间长,卫星轨道高,需要大功率发射机才能保证信号落地电平等不足。本文的低轨局部星座与BDS-2星座主要异同如下:
1) 二者都是利用卫星广播导航信号,在用户容量上不受限。
2) BDS-2采用近圆形轨道,本文采用椭圆形轨道。椭圆形轨道能够显著增大设定区域的过顶时间;另一方面,卫星在椭圆轨道中运行受地球摄动影响较大,轨道预报难度提高,轨道预报精度相比圆轨道略差,同时,在椭圆轨道运行时,卫星离地表距离变化较大,地面接收信号动态大,需要接收机具有较大的信号接收动态范围。
3) BDS-2是针对亚太地区设计的独立导航系统,其位于中高轨,单星一般在1 000 kg量级,且星座和卫星均独立设计、独立运行。
4) 本文利用低轨星座配合现有全球导航系统,实现局部导航增强,理论上可以针对任意区域进行导航增强,且通过较少卫星数量和单星就可达到目的,单星一般在100 kg量级,同时具有易于和通信卫星一体化设计等优点。
-
本文以BDS-3全球导航系统为对象研究面向区域导航增强的低轨卫星星座设计。如表1所示,BDS-3系统标称空间星座由3颗GEO卫星、3颗IGSO卫星、24颗MEO卫星组成[15]。针对BDS-3系统的导航增强问题,文献[16]给出了一种GDOP极小化的Walker组合构型计算方法,即:
$$ \sum\limits_{j = 1}^K {{N_j}{{\cos }^2}{i_j}} = \frac{1}{3}\sum\limits_{j = 1}^K {{N_j}} $$ (1) 式中,Nj为第j个单Walker构型的卫星数;K为构型数量;ij为第j个单Walker构型的轨道倾角。常见的设置是K=2,即两种倾角的Walker构型组合。如表1所示,一个为倾角45°、48颗星的Walker星座;另一个为倾角为75°、30颗星的Walker星座。
表 1 GDOP极小化的导航增强低轨星座
轨道类型 LEO LEO 卫星数量/颗 30 48 星座构型 Walker(30/6/2) Walker(48/6/1) 轨道倾角/(o) 75 45 轨道高度/km 1 058.7 1 058.7 表1中的低轨导航增强星座采用近圆形轨道,且Walker构型星座在空间上均匀分布,此类星座对南北半球有对称的观测效果,当要求覆盖区域范围较小时会造成观测空间和时间的浪费。基于以上考虑,本文引入椭圆轨道,利用其远地点特性,使其覆盖重心聚焦在需要导航增强的局部区域。
-
参考文献[12-13]中的椭圆轨道设计流程,结合区域导航增强需求,对星座设计流程进行改进,使其聚焦在需要增强的区域,在部署相同数量的低轨卫星时,达到更好的导航增强效果。
现有椭圆轨道设计流程主要满足轨道远地点与增强区域相对应这一基本要求,并未从细节上考虑如何尽可能地优化导航增强性能。如现有设计流程往往只保证最初始设计卫星的远地点限制条件,由于地球自转、轨道回归特性等因素的影响,难以保证其余卫星的星下点轨迹也能经过目标区域中心点,这给区域覆盖和导航增强性能带来损失。
基于上述分析,考虑Flower星座构型[17],它是由一组共地面轨迹[18-19]的卫星构成的星座。利用该特性进行低轨卫星椭圆轨道设计,获得的星座也具有共地面轨迹特性,能够提升导航增强性能。
-
卫星椭圆轨道运动可用开普勒参数来描述,包括:轨道升交点赤经Ω(或升交点经度λ0)、轨道倾角i、近地点幅角ω、长半径a、偏心率e和卫星的真近点角v。围绕这6个参数,面向区域导航增强的椭圆轨道低轨卫星星座设计流程如下。
1) 确定目标点经纬度(λ,
$\varphi $ )和轨道倾角i。一般选取导航增强区域的中心点为目标点。轨道倾角主要影响覆盖的纬度范围[8],根据区域的纬度跨度确定,也可根据特定的需求确定,如临界轨道等。2) 确定近地点幅角ω和真近点角v。当目标点纬度
$\varphi $ 、轨道倾角i和近地点幅角ω满足:$$ \omega = {180^ \circ }{\text{ + }}\arcsin \left( {\frac{{\sin \varphi }}{{\sin i}}} \right) $$ (2) 可以保证卫星过远地点时刻位于目标点上空。真近点角v为:
$$ v = {360^ \circ } - \omega $$ (3) 3) 确定回归系数Q。回归系数是描述轨道回归特性的参数,回归轨道是指星下点轨迹经过一段时间后重复出现的轨道。根据轨道高度所处的范围,结合式(3),中长半轴与偏心率的关系为:
$$ a = \frac{{{h_{\rm{a}}} + {h_{\rm{p}}}}}{2} + {R_{\rm{e}}},\quad e = \frac{{{h_{\rm{a}}} - {h_{\rm{p}}}}}{{2a}} $$ (4) 式中,ha和hp分别为椭圆轨道的远地点和近地点高度;Re为地球平均半径。将相关参数范围带入式(5),可计算Q的范围:
$$ \begin{split} & \qquad\qquad\qquad\quad Q = \frac{R}{X} = I + \frac{C}{X} = \\ & \qquad \frac{{n \left[ {1 - \dfrac{3}{4}{J_2}{{\left( {\dfrac{{{R_{\rm{e}}}}}{p}} \right)}^2}\left[ {( {1 - 5{c^2}} ) + \eta ( {1 - 3{c^2}} )} \right]} \right]}}{{\left[ {{\omega _{\rm{e}}} + \dfrac{3}{2}{J_2}{{\left( {\dfrac{{{R_{\rm{e}}}}}{p}} \right)}^2}nc} \right]}} \\ & n = \sqrt {\frac{\mu }{{{a^3}}}} ,\;\;\eta = \sqrt {1 - {e^2}} ,\;\;p = a (1 - {e^2}),\;\;c = \cos i \end{split} $$ (5) 式中,μ为地球引力常数;J2=0.001 082 63为地球引力带二阶谐项系数;ωe为地球自转角速度;回归系数Q由整数部分I和真分数部分C/X组成;I、C、X均为正整数,C与X互质。式(5)中的回归系数Q表示卫星在X天运行R圈后完成一次回归。
在轨道设计过程中,远地点高度ha主要限制卫星在远地点时的覆盖范围。回归系数Q可参考卫星的重访特性选择,根据文献[20]的结论,一般选取回归天数较小、回归圈数较大所对应的Q值。
4) 确定长半轴a和偏心率e。将流程3)选定的回归系数Q和近地点高度hp重新代入式(4)和式(5),可确定远地点高度ha,并得到轨道长半轴a和偏心率e。
5) 确定升交点经度λ0。根据目标点经度λ,轨道倾角i和真近点角v,利用式(6)可以计算升交点经度λ0和偏近点角E:
$$ \left\{ \begin{gathered} {\lambda _0} = \lambda - \arcsin \left( {\frac{{\tan \varphi }}{{\tan i}}} \right) + {\omega _{\rm{e}}} \frac{T}{2} \left( {1 - \frac{{\left( {E - e\sin E} \right)}}{{\text π} }} \right) \hfill \\ \tan \frac{v}{2} = \sqrt {\frac{{1 + e}}{{1 - e}}} \tan \frac{E}{2} \hfill \\ \end{gathered} \right. $$ (6) 6) 构建星座。首先,利用式(7)计算初始卫星的平近点角M0:
$$ {M_0} = E - e\sin E $$ (7) 然后,设定Flower星座的每个轨道面卫星数等于回归天数X。基于初始卫星的升交点经度λ0和平近点角M0,均匀分布的Flower星座相邻轨道面的升交点赤经差Δλ和每颗卫星的平近点角M(m, n)计算如下:
$$ \left\{ \begin{gathered} \Delta \lambda = \frac{{2{\text π} }}{P} \hfill \\ M(m,n) = {M_0} - m \Delta \lambda Q + n \frac{{2{\text π} }}{X} \hfill \\ \end{gathered} \right. $$ (8) 式中,P为轨道面数;X为回归天数;
$m=0,1, \cdots, P- 1 $ ,为轨道面编号;$n=0,1, \cdots, X-1 $ ,为轨道面内卫星编号。 -
设计一个低轨卫星星座对我国中西部地区进行导航增强,区域中心点位于(30.666 7°N, 104.067°E),区域经纬度跨度均为20°。
-
基于临界轨道方案确定轨道倾角i,以避免轨道平面拱线漂移。考虑地球引力位函数对卫星轨道的摄动时,近地点幅角、偏心率的变化率分别为[21]:
$$ \begin{split} & \;\;\;\;\;\;\;\;\;\;\;\; \dot \omega = - \frac{{3n{J_3}R_{\rm{e}}^3}}{{2{a^2}{{\left( {1 - {e^2}} \right)}^2}}}( {\frac{5}{2}{{\sin }^2}i - 2} ) \times \hfill \\ & \left[ {1 + \frac{{{J_3}{R_{\rm{e}}}}}{{2{J_2}\left( {1 - {e^2}} \right)}}( {\frac{{{{\sin }^2}i - e{{\cos }^2}i}}{{\sin i}}} )\frac{{\sin \omega }}{e}} \right] \hfill \end{split} $$ (9) $$ \dot e = \frac{{3n{J_3}R_{\rm{e}}^3\sin i}}{{4{a^3}{{\left( {1 - {e^2}} \right)}^2}}}\left( {\frac{5}{2}{{\sin }^2}i - 2} \right)\cos \omega $$ (10) 令各变化率为零,求得临界轨道倾角为i=63.43°。
-
由于范·艾伦带中的高能粒子穿透力强,对星载电子设备的损害极大,在设计卫星轨道高度时应避开范·艾伦带的内带(1 500~8 000 km)和外带(15 000~25 000 km)。选取椭圆轨道近地点高度hp=1 000 km,位于内带以下;远地点高度位于内外带之间,即8 000~15 000 km范围内。
-
1) 选取临界轨道倾角i=63.43°;
2) 根据导航增强区域中心点的经纬度坐标,即(30.666 7°N, 104.067°E),由式(2)计算出近地点幅角ω=214.769°,由式(3)计算出真近点角v=145.231°;
3) 远地点高度范围为8 000~15 000 km,由式(4)和式(5)计算回归系数Q的范围[5.015 2, 7.612 0],由于Q取7时远地点高度降低过多,影响覆盖性能,选取Q=6,即卫星运行6圈后回归,回归周期为1天。
4) 根据回归系数Q和近地点高度ha的值,利用式(5)确定远地点高度hp=11 754 km,再利用式(4)可得长半轴a=12 755 km,偏心率e=0.421 5。
5) 直接将相关参数带入式(6),计算出升交点经度为λ0=98.83°和偏近点角E=127.72°。
6) 计算初始卫星的平近点角M0=108.61°,在Flower星座框架下,选择P=9个轨道面,每个轨道面的卫星数为X=1颗。利用式(8)可得相邻轨道面升交点经度差为40°,初始卫星平近点角差
$ - {240^ \circ } $ 。由此完成了面向区域导航增强的低轨卫星星座设计,其主要参数如表2所示。
表 2 面向区域导航增强的低轨卫星星座
卫星轨道设计参数 星座设计参数 长半轴
a/km偏心率
e轨道倾角
i/(°)升交点
经度λ0/(°)近地点
幅角ω/(°)初始平近
点角M0/(°)轨道平面
数量/个轨道面内卫
星数量/个升交点
经度差/(°)平近点
角差/(°)回归
系数回归周期/
d12 755 0.4215 63.43 98.83 214.77 108.61 9 1 40 −240 6 1 关于星座设计的结果有两点需要说明。
1) 椭圆轨道高度的选择是多个因素综合考虑折中的结果。一般而言,轨道高度越低,在载波快速收敛、地面信号强度方面的增强效果越好;但同时单星覆盖范围也越小,覆盖相同区域需要的卫星数也越多。本文在选择轨道高度时,综合考虑了导航增强效果和单星覆盖区域,同时尽量避开了对卫星寿命损害极大的范·艾伦带。基于上述考虑,本文设计在增强效果、卫星数量和卫星寿命方面进行了折中,椭圆轨道远地点高度为11 754 km。根据电波在自由空间传输的损耗模型[22]:
$$ \begin{split} &\;\;\;\;\; {\text{Lfs}} = - 20\lg a(r) = \\ & 92.5 + 20{\rm{lg}}F + 20{\rm{lg}}r \end{split} $$ (11) 式中,F为载频,单位为GHz;r表示传播距离,单位为km。在传输距离损耗方面,椭圆轨道远地点信号传输损耗为81.4 dB,与之相对,BDS-3的GEO/IGSO卫星轨道高度为35 786 km,对应的传输损耗为91.1 dB;MEO卫星轨道高度为21 528 km,传输损耗为86.7 dB。由此可见,如果星上发射功率相同,在地面用户端低轨卫星的信号强度最高有10 dB提升。
2) 本文设计方案具有通用性,当重点区域的位置发生变化后,利用本文的星座设计方法可以快速计算出新的低轨卫星星座参数。至于实际星座调整过程中涉及到的指令发送、姿态控制、变速变轨等操作的原理和具体实施过程,则不在本文讨论的范围内。
-
利用覆盖重数(可见星个数)、过顶时间、GDOP等指标验证本文设计的低轨卫星星座在区域导航增强上的性能优势。覆盖重数和GDOP的仿真软件为MATLAB,过顶时间的仿真软件为STK。
-
通过覆盖重数和过顶时间验证星座设计中的共地面轨迹条件带来的覆盖优势。为便于对比,在低轨卫星数量不变的条件下,设计了基于经典Walker构型的椭圆轨道星座,它具有3个均匀分布的轨道面,每个轨道面有3颗卫星,真近点角分别为145.231°、202.556°和329.725°,其他椭圆轨道参数不变。在给定区域上对两种星座的覆盖重数和过顶时间进行仿真。仿真开始时刻为世界协调时2020年5月15日4时,最低仰角设为5°。
覆盖重数仿真时长48 h,步进为300 s。图1为两种星座的全球平均覆盖重数分布情况,图中虚线框表示以目标点为中心、经纬度跨度为20°构成的区域。图1a和图1b分别为经典Walker星座和本文设计星座的全球覆盖重数,可以发现,Walker星座的覆盖特性在经度上比较均匀,而本文设计星座则呈现出明显的区域覆盖特性。图1c给出了在设定区域上,本文设计星座相对于Walker星座在覆盖重数上的提升,说明在给定区域上,本文设计星座的覆盖重数都高于经典Walker星座。
表3为两种星座对区域中心点的过顶时间对比。本文的星座在平均过顶时间和过顶时间占比方面都优于经典Walker星座,且过顶时间更为稳定和集中。
表 3 过顶时间对比
参数 经典Walker星座 本文设计星座 平均过顶时间/s 6 118.592 6 467.094 最小过顶时间/s 1 374.066 3 025.387 最大过顶时间/s 9 026.357 8 935.949 过顶时间占比/% 28.40 30.02 -
为了验证椭圆轨道方案的优势,将其与近圆形轨道方案对比。在面向全球导航增强的经典近圆形轨道Walker星座基础上,进一步增加本文设计的椭圆轨道星座进行辅助。具体而言,选取两种Walker星座,分别为表1中由78颗低轨卫星构成的星座,即Walker 30/6/2(75°)+Walker 48/6/1(45°)和由116颗低轨卫星构成的Walker星座,即Walker 36/6/2(75°)+Walker 80/8/1(45°)。
图2为在78星Walker星座上叠加本文设计的星座后,在给定区域上相对于78星Walker星座和116星Walker星座的覆盖重数提升情况,结果验证了椭圆轨道低轨星座在区域导航增强方面的高效性:在78星近圆形轨道Walker星座上叠加本文设计的9星椭圆轨道星座,就能达到与116星近圆形轨道Walker星座等同的区域导航增强效果。
最后,结合BDS-3系统,对比本文椭圆轨道低轨星座和其他相关星座在给定区域上的GDOP。同样部署相同数量的低轨卫星,即在表1中78星近圆形轨道Walker星座上分别叠加4.1节中两种9星椭圆轨道星座。图3a在全球范围内的经、纬两个维度上对比了上述两种星座的GDOP。本文星座设计方案联合考虑了椭圆轨道和共地面轨迹两个条件,在关注的重点区域,无论是经度维度还是纬度维度,本文方法都具有更低的GDOP,表明在给定区域内星座空间构型更好,定位导航精度更高,达到了区域导航增强的目标。关注区域之外的GDOP对比情况,不是本文关注的重点(本文研究的应用场景为区域导航增强、而非全球导航增强)。如果关注的区域发生了变化,则需要更新区域信息,重新计算新的低轨卫星星座。另一方面,从图中可以看出,相对于经典Walker星座,本文低轨星座在经纬度上的GDOP增强效果并不均匀。经度方面,重点区域中心增强效果最好,两侧的增强效果逐渐减弱;纬度方面,低纬度区域的增强效果优于高纬度区域。图3b从二维平面上给出重点区域内的GDOP增强情况,结果与上面的分析一致。在重点区域的下边界中心,增强效果最明显,GDOP值相对降低达2.5%;重点区域右上角的增强效果最弱,GDOP值降低0.5%。这种增强效果不均匀的现象是由椭圆轨道形状、近地点远地点高度限制、重点区域的选取、地表曲率等多种因素造成的,同时也与星座设计方案需要兼顾覆盖重数、过顶时间、GDOP等多方面的导航增强需求有关。尽管存在增强效果不均匀的情况,但在整个区域内GDOP值都有所减小,不同程度地提高了导航定位精度。
LEO Satellite Constellation Design for Regional Navigation Augmentation
-
摘要: 部署低轨导航卫星是对现有卫星导航系统的有效补充,是实现导航增强的重要途径。研究了低轨卫星的星座设计问题,面向区域导航增强的应用需求,采用基于椭圆轨道的星座设计方案,利用卫星在椭圆轨道远地点运动速度慢、覆盖区域大、相对地面滞留时间长等特点,提升其对局部区域的覆盖能力。此外,在Flower星座框架下优化每颗卫星的轨道参数,确保其星下点轨迹重复经过重点区域的中心,延长卫星相对于给定区域的过顶时间。有效改善了过顶时间、覆盖重数、空间几何分布等影响导航性能的关键指标。在以北斗3号系统为导航增强对象的仿真对比中发现,在低轨卫星数量相同的条件下,该文的卫星星座设计方案具有更好的导航增强性能。Abstract: The deployment of low earth orbit (LEO) navigation satellites is an important supplement to the existing satellite navigation systems and also an effective approach to navigation augmentation. The design of LEO satellite constellation is studied in this paper, in order to optimize its performance in regional navigation augmentation. Different from the near-circular-orbit-based scheme adopted by most LEO constellations for global navigation augmentation, an elliptical-orbit-based constellation scheme is employed in our design. Owing to the slow-motion speed, the large coverage area, and the long residence time of satellites in the apogee of elliptical orbit, the elliptical-orbit-based scheme is specifically suitable for regional navigation augmentation. In addition, under the framework of Flower constellation, the orbit parameter of each satellite is optimized to ensure that the trajectory of its sub-satellite point repeatedly passes through the center of the region, thus extending the zenith crossing time of the satellite. By the aforementioned constellation design, the key navigation performance metrics, such as zenith crossing time, coverage multiplicity, and geometric dilution of precision, are effectively improved. In the simulations with BeiDou-3 system (BDS-3) being the navigation system to be augmented, it is validated that with the same number of LEO satellites, the proposed design achieves better navigation augmentation performance.
-
表 1 GDOP极小化的导航增强低轨星座
轨道类型 LEO LEO 卫星数量/颗 30 48 星座构型 Walker(30/6/2) Walker(48/6/1) 轨道倾角/(o) 75 45 轨道高度/km 1 058.7 1 058.7 表 2 面向区域导航增强的低轨卫星星座
卫星轨道设计参数 星座设计参数 长半轴
a/km偏心率
e轨道倾角
i/(°)升交点
经度λ0/(°)近地点
幅角ω/(°)初始平近
点角M0/(°)轨道平面
数量/个轨道面内卫
星数量/个升交点
经度差/(°)平近点
角差/(°)回归
系数回归周期/
d12 755 0.4215 63.43 98.83 214.77 108.61 9 1 40 −240 6 1 表 3 过顶时间对比
参数 经典Walker星座 本文设计星座 平均过顶时间/s 6 118.592 6 467.094 最小过顶时间/s 1 374.066 3 025.387 最大过顶时间/s 9 026.357 8 935.949 过顶时间占比/% 28.40 30.02 -
[1] 杨元喜. 综合PNT体系及其关键技术[J]. 测绘学报, 2016, 45(5): 505-510. doi: 10.11947/j.AGCS.2016.20160127 YANG Y X. Concepts of comprehensive PNT and related key technologies[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(5): 505-510. doi: 10.11947/j.AGCS.2016.20160127 [2] PAZIEWSKI I, WIELGOSZ P. Investigation of some selected strategies for multi-GNSS instantaneous RTK positioning[J]. Advances in Space Research, 2017, 59(1): 12-23. doi: 10.1016/j.asr.2016.08.034 [3] 王磊, 李德仁, 陈锐志, 等. 低轨卫星导航增强技术——机遇与挑战[J]. 中国工程科学, 2020, 22(2): 144-152. WANG L, LI D R, CHEN R Z, et al. Low earth orbiter (LEO) navigation augmentation ―opportunities and challenges[J]. Strategic Study of CAE, 2020, 22(2): 144-152. [4] 李德仁, 沈欣, 龚健雅, 等. 论我国空间信息网络的构建[J]. 武汉大学学报, 2015, 40(6): 711-715. LI D R, SHEN X, GONG J Y, et al. On construction of China ’s space information network[J]. Geomatics and Information Science of Wuhan University, 2015, 40(6): 711-715. [5] LI X, MA F, LV H, et al. LEO constellation-augmented multi-GNSS for rapid PPP convergence[J]. Journal of Geodesy, 2019, 93(5): 749-764. doi: 10.1007/s00190-018-1195-2 [6] 马福建. 低轨星座增强GNSS精密定位关键技术研究[D]. 武汉: 武汉大学, 2018. MA F J. Research on the key technologies for GNSS precise positioning augmented with LEO constellation[D]. Wuhan: Wuhan University, 2018. [7] 田野, 张立新, 边朗. 低轨导航增强卫星星座设计[J]. 中国空间科学技术, 2019, 39(6): 55-61. TIAN Y, ZHANG L X, BIAN L. Design of LEO satellites augmented-constellation for navigation[J]. Chinese Space Science and Technology, 2019, 39(6): 55-61. [8] 王慧珍, 王广兴, 牛飞, 等. 低轨卫星增强北斗三号导航星座性能分析[J]. 测绘科学, 2021, 46(1): 76-83. WANG H Z, WANG G X, NIU F, et al. Performance analysis of LEO satellites augmented BDS-3 navigation constellation[J]. Science of Surveying and Mapping, 2021, 46(1): 76-83. [9] TAINI G, PIETROPAOLO A, NOTARANTONIO A. Criteria and trade-offs for LEO orbit design[C]//Aerospace Conference. MT, USA: IEEE, 2008: 1-11. [10] 沈大海, 蒙艳松, 边朗, 等. 基于低轨通信星座的全球导航增强系统[J]. 太赫兹科学与电子信息学报, 2019, 17(2): 209-215. doi: 10.11805/TKYDA201902.0209 SHEN D H, MENG Y S, BIAN L, et al. A global navigation argumentation system based on LEO communication constellation[J]. Journal of Terahertz Science and Electronic Information Technology, 2019, 17(2): 209-215. doi: 10.11805/TKYDA201902.0209 [11] 徐哲宇, 杜兰, 刘泽军, 等. 基于Flower星座的区域导航增强低轨卫星星座设计与优化[C]//第十一届中国卫星导航年会论文集-S07卫星导航增强技术. 成都: [s.n.], 2020: 116-120. XU Z Y, DU L, LIU Z J, et al. Design and optimization of LEO regional navigation augmentation satellite constellation based on flower constellation[C]//Proceedings of the 11st China Satellite Navigation Conference. Chengdu: [s.n.], 2020: 116-120. [12] 沈欣, 姚璜, 张靖. 一种区域覆盖遥感卫星椭圆轨道设计方法[J]. 导航定位学报, 2015, 3(3): 100-104. SHEN X, YAO H, ZHANG J. An elliptical orbit design method of remote sensing satellite for regional coverage task[J]. Journal of Navigation and Positioning, 2015, 3(3): 100-104. [13] 张雅声, 张育林. 一种有效的中椭圆轨道卫星星座设计与分析[J]. 装备指挥技术学院学报, 2006, 17(4): 42-47. ZHANG Y S ZHANG Y L. Design and analysis of an effective medium ellipse orbit satellites constellation[J]. Journal of the Academy of Equipment Command & Technology, 2006, 17(4): 42-47. [14] 中国卫星导航系统管理办公室. 北斗卫星导航系统空间信号接口控制文件: GBS-SIS-ICD-2.1[S]. 北京: 中国卫星导航系统管理办公室, 2016.11. China Satellite Navigation System Management Office. BeiDou navigation satellite system signal in space interface control document (ver. 2.1): GBS-SIS-ICD-2.1[S]. Beijing: China Satellite Navigation System Management Office, 2016.11. [15] 中国卫星导航系统管理办公室. 北斗卫星导航系统公开服务性能规范: GB/T 39473-2020[S]. 北京: 中国卫星导航系统管理办公室, 2021: 5 China Satellite Navigation System Management Office. BeiDou navigation satellite system open service performance standard (ver. 3.0): GB/T 39473-2020[S]. Beijing: China Satellite Navigation System Management Office, 2021: 5. [16] 薛树强, 杨元喜. 最小GDOP组合Walker星座构型[J]. 武汉大学学报(信息科学版), 2016, 41(3): 380-387. XUE S Q, YANG Y X. Combined Walker configuration with minimal GDOP[J]. Geomatics and Information Science of Wuhan University, 2016, 41(3): 380-387. [17] WILKINS M P. The flower constellations: Theory, design process, and applications[D]. Texas: Texas A&M University, 2004. [18] JAMALIPOUR A, TUNG T. The role of satellites in global IT: Trends and implications[J]. IEEE Personal Communications, 2001, 8(3): 5-11. doi: 10.1109/98.930091 [19] 吴廷勇, 朱立东, 吴诗其. 区域覆盖卫星系统的共地面轨迹星座实现[J]. 通信学报, 2006(8): 25-30. doi: 10.3321/j.issn:1000-436X.2006.08.005 WU T Y, ZHU L D, WU S Q. Design of regional coverage satellite system with common-track constellation[J]. Journal on Communications, 2006(8): 25-30. doi: 10.3321/j.issn:1000-436X.2006.08.005 [20] 贾向华, 徐明. 太阳同步回归轨道卫星的重访特性研究[J]. 兵器装备工程学报, 2017, 38(5): 159-163. doi: 10.11809/scbgxb2017.05.037 JIA X H, XU M. Research on revisit properties for sun-synchronous RGT orbits’ satellite[J]. Journal of Ordnance Equipment Engineering, 2017, 38(5): 159-163. doi: 10.11809/scbgxb2017.05.037 [21] 杨维廉. 资源一号卫星轨道: 理论与实践[J]. 航天器工程, 2001, 10(1): 30-43. YANG W L. Orbit of CBERS-1: Theory and practice[J]. Spacecraft Engineering. 2001, 10(1): 30-43. [22] GOLDSMITH A. Wireless communications[M]. Cambridge: Cambridge University Press, 2005.