留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

路由算法在LEO通信星座构型优化中的应用

孙远辉 韩潮

王尔申, 杨福霞, 庞涛, 等 . BDS/GPS组合导航接收机自主完好性监测算法[J]. 北京航空航天大学学报, 2018, 44(4): 684-690. doi: 10.13700/j.bh.1001-5965.2017.0277
引用本文: 孙远辉, 韩潮. 路由算法在LEO通信星座构型优化中的应用[J]. 北京航空航天大学学报, 2012, (11): 1435-1439.
WANG Ershen, YANG Fuxia, PANG Tao, et al. BDS/GPS combined navigation receiver autonomous integrity monitoring algorithm[J]. Journal of Beijing University of Aeronautics and Astronautics, 2018, 44(4): 684-690. doi: 10.13700/j.bh.1001-5965.2017.0277(in Chinese)
Citation: Sun Yuanhui, Han Chao. Application of routing algorithm in configuration optimization of LEO communication satellite constellation[J]. Journal of Beijing University of Aeronautics and Astronautics, 2012, (11): 1435-1439. (in Chinese)

路由算法在LEO通信星座构型优化中的应用

详细信息
  • 中图分类号: V412.41

Application of routing algorithm in configuration optimization of LEO communication satellite constellation

  • 摘要: 低轨(LEO,Low Earth Orbit)通信星座的特殊功能决定了其构型优化需要综合考虑成本、覆盖和路由性能等因素.以往的研究往往忽略了路由性能,导致优化得到的构型不能很好的满足要求.首先给出了两个意义明确的星座路由性能的评价指标,然后通过将它们作为目标函数进行仿真,验证了其优化必要性和评价有效性.最后利用多目标遗传算法实现了星座的构型优化.仿真结果显示,将路由算法加入到星座的构型优化中,在仅对星座构型参数进行微调的前提下,不仅满足了成本、覆盖等基本要求,而且星座的路由性能大大提高.

     

  • 接收机自主完好性监测(RAIM)算法源于民航用户对卫星导航系统可靠性的苛刻要求[1]。目前,卫星导航完好性监测技术主要有3种:卫星自主完好性监测(SAIM)、利用接收机内部冗余信息的RAIM和外部辅助增强系统[2]。相对于其他2种完好性监测技术,RAIM是利用接收机的超定解对卫星测量值进行一致性校验的技术,其更能保证用户对告警时间的要求。早期的RAIM算法大多基于GPS单星座,可以满足民航航路到非精密进近阶段的完好性要求。目前,全球卫星导航系统(GNSS)发展为多星座多频点,参与定位的卫星成倍增加,接收机具有更多的冗余观测信息,可有效地减少RAIM空洞,也使RAIM技术可应用于一些对完好性性能要求更高的应用,如民航飞行中垂直引导进近(APV)Ⅰ、Ⅱ等阶段。

    多星座RAIM算法是指使用2个或2个以上的卫星导航系统的冗余观测信息进行一致性校验,以检测和识别故障。单星座系统中至少观测到5颗卫星才能进行故障检测,考虑到不同星座的时间系统之间的偏差,对于多星座需要更多冗余观测量才能进行故障检测。文献[3]对RAIM的可用性预测进行了分析。文献[4-6]提出了针对多系统多故障情况下的完好性监测算法,其中文献[4]对GPS/Galileo/GLONASS组合导航在APV-Ⅰ、APV-Ⅱ阶段利用RAIM算法进行了分析。Ene等[7-8]分析了GPS/Galileo组合导航在APV阶段的RAIM可用性。文献[9]对GPS/Galileo组合导航RAIM算法中的截止高度角、用户等效距离误差等问题进行了分析。徐君毅等[10]对APV阶段的垂直保护水平在全球的分布情况进行了研究。Lee等[11]提出了最优加权平均解(OWAS)算法,国内的研究人员也对此算法进行了扩展,郭靖等[12]采用新的分组方法,满足双星座双故障的应用场景,王式太[13]设计了基于M估计的RAIM粗差探测和剔除算法。其中,OWAS算法能够方便地兼容不同的星座,避免了不同系统间时间系统、坐标系统转换带来的误差影响,具有处理速度快、可用性高的特点,因此,本文采用OWAS算法对BDS/GPS组合导航应用于APV阶段进行了研究。

    OWAS算法是将卫星按星座分为2组,监测同时发生在同一星座中的多星故障,但此算法不能有效地对双星座同时发生故障进行检测。本文采用OWAS算法,通过优化加权因子,提高算法的可用性,并对OWAS算法进行了改进。通过实测数据验证了OWAS算法对APV阶段的可用性和双星座双故障的检测。

    OWAS RAIM算法原理如图 1所示。

    图  1  OWAS RAIM算法原理
    Figure  1.  OWAS RAIM algorithm principle

    加权平均后的定位结果为

    (1)

    式中:加权因子0≤r≤1;xn(n=1, 2)分别为BDS和GPS的定位结果。

    图 1中:dn(n=1, 2)分别为xnxA的间隔距离。对于加权平均解,要求以0.999的概率满足:|xAx|≤|xAxn|+|xnx|≤T+an。其中:T(n=1, 2)为检验门限,an为无故障假设条件下定位误差上限。所以,无故障卫星时,每个导航系统内的保护水平为:PLn=T+an,PL=max(PL1, PL2)。OWAS算法在只考虑垂直方向上的定位精度与完好性要求时,可以归纳为如下问题:

    (2)

    式中:VPL为垂直保护级别。

    同时,满足以下条件时,RAIM算法可用:

    (3)

    (4)

    (5)

    (6)

    式中:σV为垂直方向上的定位误差标准差;σV-max为垂直方向上的定位误差标准差的最大值;dmajor为最终定位解水平方向上的定位误差标准差;dmajor-max为飞行阶段规定的95%水平方向上的定位误差标准差的最大值;VAL为垂直保护门限;HPL为水平保护级别;HAL为水平保护门限[10]

    σVdmajor的标准由以下2个操作需求推导:

    1) 95%的导航系统误差(Navigation System Errors,NSE)边界:1.96σV=4m,1.96dmajor=16m(1.96对应于高斯分布假设的95%分布)。

    2) NSE边界满足无故障条件完好性可用性时:KVσV≤VAL(垂直方向),KHdmajor≤HAL(水平方向),当Pr{HMI}=10-8/approach(Pr为概率, HMI为危险误导信息, approach为进近)时,KV=5.73,KH=6.47[10]

    σV-maxdmajor-max的值分别为

    OWAS算法为

    (7)

    式中:w1为GPS定位解在最终解中所占权重;w2为BDS定位解在最终解中所占权重;,且0≤r≤1。

    OWAS算法的位置误差与加权因子的关系如下:

    (8)

    式中:Δx1为GPS定位解与真实位置之间的误差;Δx2为BDS定位解与真实位置之间的误差;ΔxA为最终加权平均定位解与真实位置之间的误差。

    加权最小二乘定位解的位置误差[14]表达式如下:

    (9)

    (10)

    式中:G为伪距误差与接收机位置及钟差的关联矩阵;W为加权矩阵;ΔR为伪距误差矢量;Pn为定位误差的协方差矩阵(n=1, 2分别表示GPS和BDS)。

    最优解和各导航系统定位解误差互相关表达式如下:

    (11)

    (12)

    假设GPS与BDS伪距误差不相关,则ER2·ΔR1T)=0,ER1ΔR2T)=0。由此,得

    (13)

    (14)

    因此,ExAΔx1T)=rP1ExAΔx2T)=(1-r)P2,即

    (15)

    (16)

    同理,加权平均解xA的定位误差协方差矩阵PA如下:

    (17)

    协方差矩阵Jn

    (18)

    将式(10)、式(11)、式(15)、式(17)代入式(18)中,可得

    (19)

    式中:PS=P1+P2

    同理

    (20)
    1.1.1   加权因子确定保护级别

    (21)

    (22)

    式中:Q-1为标准单边累加分布函数的逆函数;Pfd为虚警率(参考国际民航不同进近阶段导航性能要求,虚警率与连续性保持一致。虚警率除以4是因为总的允许虚警率在垂直和水平方向上平均分配[10]);Pmd为最大允许漏警率(完好性风险等于漏警率与卫星故障率的乘积)。

    (23)

    (24)

    (25)

    由式(19)和式(20)可以得出

    (26)

    (27)

    VPL1和VPL2分别为

    (28)

    (29)

    图 2可以看出,VPL1是加权因子的一元线性递减函数;VPL2是加权因子的一元线性递增函数。在OWAS算法中,VPL=max(VPL1, VPL2),当VPL1=VPL2时,得到最小化的VPL,此时可得

    (30)

    (31)
    图  2  VPL1、VPL2与加权因子的函数
    Figure  2.  VPL1 and VPL2 as a function of weight factor
    1.1.2   加权因子确定误差精度σVdmajor

    OWAS位置误差的协方差矩阵PA=ExAΔxAT)=r2P1+(1-r)2P2

    σVdmajor

    (32)

    (33)

    由式(3)、式(4)可知,RAIM算法可用要满足:

    1.1.3   确定最优加权因子

    OWAS最重要的是在降低保护水平的条件下满足精度的要求。

    图 3所示,rmin-VPL表示垂直保护水平最小时的加权因子,但此时不满足σV的标准,而rmin-σV满足最小化σV,却不满足VPL≤VAL,rbest-int表示在保证精度要求下修正加权因子rmin-VPL来最小化VPL,rbest-acc表示通过修正最小化精度加权因子σV保证可用性的要求。

    图  3  加权因子的计算
    Figure  3.  Determination of weight factor

    APV飞行阶段时需更大可能的保证算法的可用性,故采用rbest-int作为最终加权因子,则

    (34)

    式(23)和式(24)代入到式(34)中,可得

    (35)

    由式(35)计算可得2个解,如下:

    (36)

    (37)

    由于r1rr2,且0≤r≤1,所以

    (38)

    式中:racc, min=max(r1, 0);racc, max=min(r2, 1)。

    racc, minrmin-VPLracc, max时,加权因子r如下:

    (39)

    rmin-VPLracc, min时,加权因子r如下:

    (40)

    racc, maxrmin-VPL时,加权因子r如下:

    (41)

    OWAS RAIM算法是以GPS、BDS到xA的距离作为故障检测的检验统计量[15-16],按照国际民航组织(ICAO)要求的性能标准得到检验门限。当双星座同时出现故障时,2个检验统计量均会在故障产生的时刻超过检验门限。当单星座产生故障时,故障星座的定位解使得OWAS与无故障星座的距离增加,容易产生误警,此时采用加权最小二乘残差的故障检测对OWAS定位域的故障检测进行辅助,完成单星座故障的检测。

    1.2.1   故障检测

    GNSS伪距观测线性化方程为

    (42)

    式中:ε为观测噪声矢量。

    对式(42)运用加权最小二乘法求解,得到的解为

    (43)

    式中:B=(GTWG)-1GTW

    伪距残差矢量v

    (44)

    矩阵C的协因数矩阵为

    伪距残差平方和为

    (45)

    在卫星导航系统正常工作的情况下,各卫星伪距残差都较小,伪距残差平方和也较小;当某个测量伪距存在较大偏差时,伪距残差平方和会变大。假设卫星导航系统正常工作时,观测噪声ε中的各个分量是均值为0、方差为σ2的正态分布随机误差,且它们相互独立。依据统计分布理论,若ε的均值为0,SSEWLS服从自由度为(n-4)的χ2分布;若ε的均值不为0,则SSEWLS服从自由度为(n-4)的非中心化χ2分布。对SSEWLS作二元假设:

    1) 无故障假设H0:E(ε)=0,则SSEWLS~χ2(n-4)。

    2) 有故障假设H1E(ε)≠0,则SSEWLS~χ2(λ, n-4)。

    (46)

    由式(46)可以解得tWLS值。

    1.2.2   故障识别

    根据残差和观测误差的关系式,定义故障识别的检验统计量为

    (47)

    式中:i=1, 2, …, n,并且n为可见卫星数目;vi为利用加权最小二乘法求得的伪距残差;Qii为矩阵Q的第i行第i列;di~N(0, 1)。

    对统计量di做二元假设:

    1) H0(无故障):E(εi)=0,则di~N(0, 1)。

    2) H1(有故障):E(εi)≠0,则di~N(δi, 1)。

    其中,δi为统计量偏移参数。给定总体的误警率PFA,可计算出对应的检验门限T1WLS。计算如下:

    (48)

    由式(48)可解得门限值T1WLS

    为检验统计量,对应的检验门限值为T1WLS。比较检验统计量di与检验门限T1WLS进行卫星故障识别。

    采用IGS(International GPS Service)2016年7月31日0时的BDS、GPS系统RINEX数据对本文算法进行仿真验证。

    在GPS可见星中的第16颗卫星、BDS可见星中的第7颗卫星中,分别在500~900s时刻人为加入30、50、70m伪距偏差。

    图 4(a)(b)分别为未加伪距偏差时GPS和BDS的垂向定位域检验统计量。图 5(a)(b)分别为在500 s开始加入30、50、70 m伪距偏差时GPS、BDS的检验统计量。可知,当加入伪距偏差时,检验统计量值会发生跳变; 未加入伪距偏差时,检验统计量值波动较小。

    图  4  BDS/GPS未加入伪距偏差的检验统计量
    Figure  4.  BDS/GPS test statistics with no pseudo-range bias
    图  5  BDS/GPS加入不同伪距偏差的检验统计量
    Figure  5.  BDS/GPS test statistics with different pseudo-range bias

    双故障卫星检测结果如图 6所示。

    图  6  BDS/GPS故障卫星检测
    Figure  6.  BDS fault satellite detection

    图 6可知,在GPS和BDS系统的第16颗卫星和第7颗卫星中加入50m伪距偏差,检验统计量大于检验门限,可以对故障进行检测。

    双星座双故障算法识别率结果如图 7所示。可知,当伪距偏差值较小时,不能识别出故障卫星;当故障卫星伪距偏差在30m以上时,BDS/GPS能识别全部故障。

    图  7  故障识别率
    Figure  7.  Fault identification rate

    本文研究了一种BDS/GPS组合导航RAIM算法,对算法进行了详细的理论分析,并利用IGS网站数据对算法进行验证。

    1) 本文算法对BDS/GPS双卫星故障检测与识别是可用的。

    2) 当伪距偏差在30m以上时可有效地检测与识别。

    研究结果对多星座组合导航接收机自主完好性监测算法研究有一定的参考价值。

  • [1] 郦苏丹,朱江,李广侠.采用遗传算法的低轨区域通信星座优化设计[J].通信学报,2005,26(8):122-128 Li Sudan,Zhu Jiang,Li Guangxia.Optimization of LEO r-egional communication satellite constellation with GA algor-ithm [J].Journal of Communications,2005,26(8):122-128(in Chinese)
    [2] 刘文,张育林,刘昆.基于多目标进化算法的卫星通信星座优化设计[J].宇航学报,2008,29(1):95-99 Liu Wen,Zhang Yulin,Liu Kun.Optimization of multi obj-ective evolutionary algorithm-based communication satellite constellation [J].Journal of Astronautics,2008,29(1):95-99(in Chinese)
    [3] 郦苏丹,朱江,李广侠.基于多目标进化算法的低轨区域通信星座优化设计[J].解放军理工大学学报(自然科学版),2005,6(1):1-6 Li Sudan,Zhu Jiang,Li Guangxia.Multi-object optimization of LEO regional communication satellite constellation [J].Journal of PLA University of Science and Technology (Natural Science Edition),2005,6(1):1-6(in Chinese)
    [4] 郦苏丹,朱江,李广侠.采用遗传算法的中轨区域通信星座优化设计[J].系统仿真学报,2005,17(6):1366-1369 Li Sudan,Zhu Jiang,Li Guangxia.Optimization of MEO regional communication satellite constellation with GA algo-rithm [J].Journal of System Simulation,2005,17(6):1366-1369(in Chinese)
    [5] 王桂华.军用通信卫星星座设计与仿真研究.长沙:国防科技大学电子科学与工程学院,2006 Wang Guihua.The design and simulation research in milit-ary communication satellite constellation.Changsha:School of Electronic Science and Engineering,National University of Defense Technology,2006(in Chinese)
    [6] 贺泉,韩潮.基于蚁群算法的多任务导航星座载荷配置[J].北京航空航天大学学报,2007,33(10):1154-1157 He Quan,Han Chao.Payload configuration for multi-mission navigation constellation based on ant colony algorithm [J].Journal of Beijing University of Aeronautics and Astronautics,2007,33(10):1154-1157(in Chinese)
    [7] Kelly C W.Minimizing the cost of availability of services from a constellation of satellites.Los Angeles:School of Engineering,University of Southern California,2003
    [8] 蒙波,伊成俊,韩潮.基于多目标粒子群算法的导航星座优化设计[J].航空学报,2009,30(7):1284-1291 Meng Bo,Yi Chengjun,Han Chao.Optimization of navigation satellite constellation by multi-objective particle swarm algorithm [J].Acta Aeronautics et Astronautics Sinica,2009,30(7):1284-1291(in Chinese)
    [9] 孙远辉,韩潮.星间链路及区域不均匀对路由性能的影响[J].系统仿真学报,2011,23(S1):219-223 Sun Yuanhui,Han Chao.Effects of inter-satellite link and regional heterogeneity on routing performance [J].Journal of System Simulation,2011,23(S1):219-233(in Chinese)
  • 加载中
计量
  • 文章访问数:  1394
  • HTML全文浏览量:  26
  • PDF下载量:  644
  • 被引次数: 0
出版历程
  • 收稿日期:  2011-06-16
  • 网络出版日期:  2012-11-30

目录

/

返回文章
返回
常见问答