发明名称 一种卫星导航系统的完好性监测可用性确定方法
摘要 本发明提供了一种卫星导航系统的完好性监测可用性确定方法,综合利用马尔科夫链(MC)、贝叶斯网络(BN)、中断分析、应力——强度分析,通过分析单机设备可用性、监测站可用性以及空间信号监测精度(SISMA)仿真,确定系统级完好性监测的可用性,为系统级完好性监测设计优化、以及单机与监测站可用性的提高提供依据。该方法可进一步推广于卫星轨道精密测定与长期预报、时间同步与钟差预报、导航信息注入、星座运行管理等任务的可用性确定,为卫星导航系统的研制建设提供重要专业技术支持。
申请公布号 CN104215974B 申请公布日期 2017.01.11
申请号 CN201410469448.8 申请日期 2014.09.15
申请人 中国航天标准化研究所 发明人 杨卓鹏;郑恒;顾长鸿;王晋婧;薛恩;龚佩佩;李海生;周波
分类号 G01S19/20(2010.01)I 主分类号 G01S19/20(2010.01)I
代理机构 中国航天科技专利中心 11009 代理人 安丽
主权项 一种卫星导航系统的完好性监测可用性确定方法,其特征在于步骤如下:(1)根据卫星导航系统的完好性监测过程中单机设备的运行与故障情况,构建基于马尔科夫链的单机设备可用性模型,确定单机设备可用性;(2)确定监测站的可用性,具体为:(2.1)根据所有监测站内的单机设备的可用性以及单机设备在监测站中的逻辑关系,构建基于贝叶斯网络的监测站可用性模型,所述逻辑关系包括串联、并联、k/n、桥联和网络;基于贝叶斯网络的监测站可用性模型为:<img file="FDA0001084760800000011.GIF" wi="422" he="116" />其中,pa<sub>i</sub>是x<sub>i</sub>的父节点集,x<sub>i</sub>表示第i个单机设备,其可用性为A<sub>ei</sub>(t);p(X)为单机设备的联合概率分布,X={x<sub>1</sub>,x<sub>2</sub>,...,x<sub>n</sub>};(2.2)根据步骤(2.1)中确定的基于贝叶斯网络的监测站可用性模型,确定监测站的可用性A<sub>m</sub>(t),该可用性A<sub>m</sub>(t)=p(X);(3)确定卫星导航系统的完好性监测可用性,具体为:(3.1)根据监测站分布及监测站中断情况,仿真确定不同中断情况下的空间信号监测精度;所述中断是指监测站整体停止工作;(3.1.1)利用监测站的观测数据计算卫星轨道误差和卫星钟差误差组成整体的误差矢量x,历元t时刻卫星j可视域内有n个监测站接收卫星信号,利用星历和监测站坐标计算站星几何距离,确定伪距观测方程并将其线性化得:<img file="FDA0001084760800000012.GIF" wi="1077" he="286" />式中,<img file="FDA0001084760800000013.GIF" wi="547" he="63" />为监测站的测量残差向量;<img file="FDA0001084760800000014.GIF" wi="741" he="83" />为卫星j到监测站i的单位矢量;x=[Δx<sub>s</sub>,Δy<sub>s</sub>,Δz<sub>s</sub>,Δt<sub>s</sub>]<sup>T</sup>为卫星在地心坐标下的轨道误差Δx<sub>s</sub>、Δy<sub>s</sub>、Δz<sub>s</sub>和钟差Δt<sub>s</sub>,<img file="FDA0001084760800000015.GIF" wi="420" he="70" />为n个监测站在伪距域上产生的误差,服从正态分布;误差矢量x及其协方差误差阵<img file="FDA0001084760800000021.GIF" wi="54" he="63" />为:<img file="FDA0001084760800000022.GIF" wi="557" he="70" /><img file="FDA0001084760800000023.GIF" wi="387" he="62" />式中,D<sub>b</sub>为卫星导航系统正常状态下伪距域噪声的标准差,为先验信息;将<img file="FDA0001084760800000024.GIF" wi="58" he="63" />投影到用户位置得:<img file="FDA0001084760800000025.GIF" wi="443" he="71" />式中,<img file="FDA0001084760800000026.GIF" wi="128" he="70" />表示用户估计误差,a<sub>uS</sub>=(e<sub>uX</sub>,e<sub>uY</sub>,e<sub>uZ</sub>,1)<sup>T</sup>表示用户到卫星方向的单位向量,e<sub>uX</sub>,e<sub>uY</sub>,e<sub>uZ</sub>表示用户与卫星的向量信息,取最差用户位置处投影值为:<img file="FDA0001084760800000027.GIF" wi="422" he="63" />(3.1.2)依照步骤(3.1.1)的计算最差用户位置处投影值SISMA的方式,在所有监测站均正常、任意一个监测站发生监测中断、任意两个监测站同时发生监测中断、……、任意n‑4个监测站同时发生中断的情况下,分别进行中断仿真,确定不同中断组合情况下的空间信号监测精度,即SISMA,n为监测站数量;(3.2)通过公式<img file="FDA0001084760800000028.GIF" wi="1261" he="87" />将步骤(3.1)中确定的不同中断组合情况下的空间信号监测精度进行归一化处理;其中,SISMA<sub>BN</sub>为各监测站之间的逻辑关系,<img file="FDA0001084760800000029.GIF" wi="379" he="87" />为应力函数,σ为应力,g(δ)为卫星导航系统要求的空间信号监测精度值,δ为强度;(3.3)构建基于贝叶斯网络的卫星导航系统的完好性监测可用性模型:<img file="FDA00010847608000000210.GIF" wi="430" he="115" />其中,ma<sub>i</sub>是y<sub>i</sub>的父节点集,且ma<sub>i</sub>=A<sub>mi</sub>(t),即ma<sub>i</sub>为第i个监测站的可用性,y<sub>i</sub>表示第i个监测站;m(Y)为监测站的联合概率分布,Y={y<sub>1</sub>,y<sub>2</sub>,...,y<sub>n</sub>};(3.4)根据步骤(3.3)中确定的基于贝叶斯网络的完好性监测可用性模型,确定完好性监测的可用性A<sub>integrity</sub>(t),该可用性即为A<sub>integrity</sub>(t)=m(Y)。
地址 100071 北京市丰台区小屯路89号