一种导航卫星精密定轨的任务成功性确定方法与流程

本发明提供了一种导航卫星精密定轨的任务成功性确定方法，属于卫星导航系统可用性技术领域。

背景技术：

定轨的任务成功性是指导航卫星在任务开始时处于可用状态的情况下，在规定的任务剖面中的任一时刻，能够使用且能完成规定功能的能力。与定轨的任务成功性相关因素包括：运行时间、硬件短期故障、硬件长期故障、软件短期故障、软件长期故障、故障检测时间、软硬件耦合故障等因素。

技术实现要素：

本发明的技术解决问题是：克服现有技术的不足，提供一种导航卫星精密定轨的任务成功性确定方法。

本发明的技术解决方案是：一种导航卫星精密定轨的任务成功性确定方法，步骤如下：

(1)根据导航监测站中断发生的原因，确定导航监测站的各类相关状态；

(2)根据上述各类相关状态构建基于广义随机petri网的导航监测站可用性模型，确定模型的稳定状态及状态转换关系；

(3)根据广义随机petri网同构马尔科夫链，确定导航监测站微分方程，进而确定每个导航监测站可用性pupi；

(4)根据导航卫星精密定轨中的动力学方程和观测方程，确定导航卫星的平均用户距离精度ure计算公式；

(5)根据监测站分布及监测站中断情况，对在所有监测站均正常情况下、当任意一个监测站发生中断、任意两个监测站同时发生中断、……、任意n-4个监测站同时发生中断时各种情况进行仿真，根据仿真结果结合步骤(5)的计算公式确定上述各种中断组合情况下的实际平均用户距离精度；

(6)将步骤(5)中不同中断组合情况下的实际平均用户距离精度进行归一化处理；

(7)根据归一化结果，构建基于贝叶斯网络的导航卫星精密定轨的任务成功性模型：其中，pai是xi的父节点集，其输入值为可用性为xi表示第i个监测站，x＝{x1,x2,...,xn}；导航卫星精密定轨的任务成功性pd即为pd＝p(x)。

进一步的，当导航卫星精密定轨的任务成功性不满足要求时，识别薄弱导航监测站，对薄弱监测站进行优化后重新执行步骤(2)-(7)，直至任务成功性满足要求。

进一步的，步骤(2)中的导航监测站的各类相关状态包括可用状态pup，故障状态pdown；硬件或软件故障选择状态phos，硬件短期或长期故障选择状态psolh，硬件短期故障状态psh，硬件长期故障状态plh，备件选择状态psonh，备件等待状态pswh，硬件长期故障维修等待状态prwh；软件短期或长期故障选择状态psols，软件故障转移选择状态psths，软件长期故障状态pls，软件短期故障选择状态prois，软件重启故障状态prs，软件自修复状态pis；软件故障向硬件故障转移状态psth。

进一步的，步骤(2)中模型的稳定状态包括：可用状态pup、故障状态pdown，硬件短期故障状态psh，硬件长期故障状态plh，备件等待状态pswh，硬件长期故障维修等待状态prwh，软件长期故障状态pls，软件重启故障状态prs，软件自修复状态pis，软件故障向硬件故障转移状态psth。

进一步的，所述pup的值在0～1之间。

进一步的，同构马尔科夫链的原则如下：

将基于广义随机petri网的导航监测站可用性模型中的库所状态分为稳定状态和不稳定状态，稳定状态是指库所输出有向弧指向时间变迁，不稳定状态是指库所输出有向弧指向瞬时变迁；

将所有稳定状态库所设置为马尔科夫状态，各稳定状态库所之间的状态转换为马尔科夫转移，相关转移值称为转移概率；

转移概率的确定过程为：a)若稳定状态库所之间只存在时间变迁，则转移概率为时间变迁对应故障率，b)若稳定状态库所之间存在时间变迁和瞬时变迁，则转移概率为时间变迁对应故障率与瞬时变迁对应概率的乘积。

进一步的，步骤(4)中的导航监测站微分方程：

pup+pdown+psh+plh+pswh+prwn+psth+pls+prs+pis＝1

式中，为导航监测站平均故障率；

λfd为导航监测站故障检测率；

为导航监测站平均管理延误故障率；

为导航监测站平均备件供应反应率；

λtsth为导航监测站软件故障传播至硬件故障的故障率；

为导航监测站硬件短期平均修复率；

为导航监测站平均重启修复率；

为导航监测站平均自修复率；

为导航监测站硬件长期平均修复率；

为导航监测站软件长期平均修复率；

fh为导航监测站硬件故障概率；

fs为导航监测站软件故障概率；

fsh为导航监测站硬件短期故障概率；

flh为导航监测站硬件长期故障概率；

fsph为导航监测站备件保障率；

fnsph为导航监测站备件非保障率；

fls为导航监测站软件长期故障概率；

fss为导航监测站软件短期故障概率；

fsths为导航监测站发生软硬件耦合故障概率；

fnsths为导航监测站不发生软硬件耦合故障概率；

frs为导航监测站软件重启类概率；

fis为导航监测站软件自修复类概率。

进一步的，fh+fs＝1，fsh+flh＝1，fsph+fnsph＝1，fls+fss＝1，fsths+fnsths＝1，frs+fis＝1。

进一步的，中地球轨道卫星、对地静止轨道卫星、倾斜地球同步轨道卫星的平均用户距离精度uremeo、uregeo、ureigso计算公式如下：

其中，卫星轨道与标准轨道比对时，轨道在径向r误差为dr、切向t误差为dt和法向n误差为dn；

步骤(5)中的导航卫星的实际平均用户距离精度ure计算选择上述任意一种计算公式进行计算。

进一步的，步骤(6)中采用如下的公式进行归一化处理：

其中cure为平均用户距离精度转换值，表示各监测站之间的逻辑关系，f(σst)为步骤(5)中计算的实际平均用户距离精度，σst为应力，g(δin)为导航卫星要求的平均用户距离精度值，δin为强度。

本发明与现有技术相比的优点在于：

(1)本发明综合综合利用广义随机petri网和贝叶斯网络开展精密定轨的任务成功性建模与分析。考虑监测站硬件故障、软件故障、软硬件耦合故障等多种因素，利用广义随机petri网构建监测站可用性模型，有效反映监测站各状态之间的转换关系；结合精密定轨算法、中断分析、以及应力——强度算法，将定轨性能数据归一化处理，构建基于贝叶斯网的精密定轨的任务成功性模型；识别模型薄弱监测站。

(2)本发明综合考虑了导航监测站运行时间，中断状态关注硬件短期故障、硬件长期故障、软件短期故障、软件长期故障、故障检测时间、软硬件耦合故障等因素，并建立了各因素与导航卫星精密定轨的任务成功性pd的映射关系。

(3)利用本方法和模型，不仅可以确定导航卫星精密定轨的任务成功性，还可进一步推广于导航信息上行注入、完好性监测等关键任务的任务成功性确定。

附图说明

图1是基于广义随机petri网的导航监测站可用性模型示意图；

图2是导航监测站各类有效状态转换示意图；

图3是导航监测站可用性仿真结果示意图；

图4是构建基于贝叶斯网络的导航卫星精密定轨任务成功性模型示意图；

图5是导航卫星精密定轨任务成功性可用性结果示意图；

图6是薄弱导航监测站识别结果示意图；

图7是改进后的导航卫星精密定轨任务成功性可用性结果示意图；

图8是本发明方法实施过程示意图。

具体实施方式

本发明提供了一种导航卫星精密定轨的任务成功性确定方法；所述定轨的任务成功性是指导航卫星在任务开始时处于可用状态的情况下，在规定的任务剖面中的任一时刻，能够使用且能完成规定功能的能力；本发明方法如图8所示，步骤如下：

(1)确定导航监测站可用性

将导航监测站mi(i＝1,2,...n，n为监测站数量)分为运行和中断两种状态，其中运行状态关注运行时间因素，中断状态关注硬件短期故障、硬件长期故障、软件短期故障、软件长期故障、故障检测时间、软硬件耦合故障等因素。

(1.1)确定导航监测站各类相关状态和状态参数。

(1.1.1)导航监测站各类相关状态包括：可用状态pup，故障状态pdown；硬件或软件故障选择状态phos，硬件短期或长期故障选择状态psolh，硬件短期故障状态psh，硬件长期故障状态plh，备件选择状态psonh，备件等待状态pswh，硬件长期故障维修等待状态prwh；软件短期或长期故障选择状态psols，软件故障转移选择状态psths，软件长期故障状态pls，软件短期故障选择状态prois，软件重启故障状态prs，软件自修复状态pis；软件故障向硬件故障转移状态psth。

(1.1.2)导航监测站各类状态参数包括，运行时间对应参数：平均故障间隔时间mtbfm(为监测站平均故障率)；故障检测对应参数：故障检测时间tfd(tfd＝1/λfd，λfd为故障检测率)，硬件故障概率fh、软件故障概率fs；硬件短期故障对应参数：硬件短期故障概率fsh，硬件短期平均修复时间mttrsh(为硬件短期平均修复率)；硬件长期故障对应参数：硬件长期故障概率flh，平均管理延误时间madtlh(为平均管理延误故障率)，备件保障率fsph(fnsph＝1-fsph)，平均备件供应反应时间msrtlh(为平均备件供应反应率)，硬件长期平均修复时间mttrlh(为硬件长期平均修复率)；软件短期故障对应参数：软件短期故障概率fss，软件重启类概率frs，软件自修复类概率fis，平均重启修复时间mttrrs(为平均重启修复率)，平均自修复时间mttris(为平均自修复率)；软件长期故障对应参数：软件长期故障概率fls，软件长期平均修复时间mttrls(为软件长期平均修复率)；软硬件耦合故障对应参数：软硬件耦合故障发生概率fsths(fnsths＝1-fsths，fnsths为不发生软硬件耦合故障概率)，软件故障传播至硬件故障时间tsth(tsth＝1/λtsth，λtsth为软件故障传播至硬件故障的故障率)。

(1.2)根据上述各类相关状态构建基于广义随机petri网的导航监测站可用性模型，确定模型的稳定状态及状态转换关系；

模型的稳定状态包括：可用状态pup，故障状态pdown，硬件短期故障状态psh，硬件长期故障状态plh，备件等待状态pswh，硬件长期故障维修等待状态prwh，软件长期故障状态pls，软件重启故障状态prs，软件自修复状态pis，软件故障向硬件故障转移状态psth。短期故障是指导航监测站可自行修复的故障，所述长期故障是指导航监测站需人工修复的故障。

上述状态转换关系为petri网相关公知常识，此处不过多进行赘述。

(1.3)根据广义随机petri网同构马尔科夫链，确定导航监测站微分方程，进而确定导航监测站可用性pup及模型各稳定状态值。

同构马尔科夫链的原则如下：

(a)将基于广义随机petri网的导航监测站可用性模型中的库所状态分为稳定状态和不稳定状态，稳定状态是指库所输出有向弧指向时间变迁，不稳定状态是指库所输出有向弧指向瞬时变迁；

(b)将所有稳定状态库所设置为马尔科夫状态，各稳定状态库所之间的状态转换为马尔科夫转移，相关转移值称为转移概率；

转移概率的确定过程为：a)若稳定状态库所之间只存在时间变迁，则转移概率为时间变迁对应故障率，b)若稳定状态库所之间存在时间变迁和瞬时变迁，则转移概率为时间变迁对应故障率与瞬时变迁对应概率的乘积。

建立的导航监测站微分方程为：

pup+pdown+psh+plh+pswh+prwn+psth+pls+prs+pis＝1

(2)确定导航卫星精密定轨的任务成功性，具体为：

(2.1)构建精密定轨算法的动力学方程和观测方程，分别为：动力学方程为式中：x为状态矢量，包括卫星在直角坐标系的三维位置分量、三维速度分量、卫星动力学参数以及不在卫星运动方程中出现的待估几何参数。卫星状态向量对时间的偏导数是t时刻卫星状态x的函数，称为状态函数，记作f(x,t)。卫星在t0时刻的状态量x(t0)记为x0；观测方程为yi＝g(xi,ti)+εi，式中：yi为第i个观测量，g(xi,ti)为ti时刻状态量的函数，即观测量理论值计算公式，εi为观测噪声。对动力学方程和观测方程进行线性化，确定卫星轨道，

与标准轨道比对，分别得出轨道在径向r误差dr、切向t误差dt和法向n误差dn。根据选择模型的参数，确定中地球轨道卫星meo、静止轨道卫星geo、倾斜地球同步轨道卫星igso的实际平均用户距离精度ure为

式中uremeo、uregeo和ureigso分别表示中地球轨道卫星、对地静止轨道卫星和倾斜地球同步轨道卫星的平均用户距离精度。实际应用过程中择一计算实际平均用户距离精度ure即可。

(2.2)根据监测站分布及监测站中断情况，仿真确定在所有监测站均正常、任意一个监测站发生中断、任意两个监测站同时发生中断、……、任意n-4个监测站同时发生中断的情况下的平均用户距离精度，n为监测站数量。

(2.3)通过公式(cure的值在0～1之间)将步骤(2.1)中确定的不同中断组合情况下的平均用户距离精度进行归一化处理；其中cure为平均用户距离精度转换值，表示各监测站之间的逻辑关系，f(σst)为导航卫星实际平均用户距离精度(即2.1中计算)，σst为应力，g(δin)为导航卫星要求的平均用户距离精度值，δin为强度。

(2.4)构建基于贝叶斯网络的导航卫星精密定轨的任务成功性模型：其中，pai是xi的父节点集，即可用性xi表示第i个监测站，x＝{x1,x2,...,xn}。

(2.5)根据步骤(2.4)确定的基于贝叶斯网络的导航卫星精密定轨的任务成功性模型，确定导航卫星精密定轨的任务成功性pd，该任务成功性即为pd＝p(x)。

(3)构建导航卫星精密定轨的任务成功性故障诊断模型：p(mi＝0|pd＝0)＝p(mi＝0)p(pd＝0|mi＝0)/p(pd＝0)，其中，“0”表示故障状态，mi表示第i个监测站，p(mi＝0|pd＝0)表示监测站mi的故障概率。通过各监测站的故障概率，识别薄弱监测站，当导航卫星精密定轨的任务成功性不满足要求时，对薄弱监测站进行优化，针对故障概率高的监测站通过延长运行时间、减少故障时间等措施进行设计优化。重新执行步骤(1.2)～(2.5)直至任务成功性满足要求。

实施示例：仿真构建卫星导航系统模型，模型中共7个监测站，选择meo卫星进行仿真。

(a)确定导航监测站可用性；根据步骤(1)构建基于广义随机petri网的导航监测站可用性，如图1所示。确定导航监测站各有效状态转换关系，如图2所示。表1所示为导航监测站m1的各参数及指标，由此确定导航监测站可用性仿真结果如图3所示。同理可确定导航监测站m2～m7的可用性。

表1监测站参数及指标

(b)确定导航卫星精密定轨的任务成功性；根据步骤(2)构建基于贝叶斯网络的导航卫星精密定轨的任务成功性模型，如图4所示；确定定轨的任务成功性pd，仿真结果如图5所示。开展导航卫星精密定轨的任务成功性故障诊断，识别薄弱监测站，如图6所示。针对薄弱监测站m1和m7进行设计优化，例如提升平均故障间隔时间、降低故障检测时间等，重新计算导航卫星精密定轨的任务成功性，如图7所示。由图7可知，相比于原方案，优化后的导航卫星精密定轨的任务成功性存在较为显著的提升。

本发明未详细说明部分属于本领域技术人员公知常识。