一种基于北斗全链条故障激励的完好性测试方法与流程

1.本发明涉及完好性风险测试技术领域，尤其涉及一种基于北斗全链条故障激励的完好性测试方法。


背景技术：

2.北斗导航卫星系统(global navigation satellite system，gnss)涵盖我国海、陆、空等多个军民融合典型应用场景，例如海洋生态应用、智慧农机协同作业、舰载机批次全自动引导等，为保障北斗导航服务可信性，星基增强系统(satellite based augmentation system，sbas)、地基增强系统(ground based augmentation system，gbas)等兼顾精度与完好性的gnss增强系统，被广泛应用于民航等与生命安全相关的导航应用领域。
3.目前我国的北斗星基增强系统(bdsbas)正处于关键建设时期，亟需可靠测试方案，以开展系统完好性性能评估，为系统服务性能优化提供参考。北斗导航定位涉及从卫星到信号到应用平台的整条服务链路，鉴于北斗导航“卫星-信号-应用平台”全链条服务体系潜在面临多类型故障场景，同时为满足北斗导航定位服务高可信性需求。


技术实现要素：

4.为解决上述技术问题，本发明提出了一种基于北斗全链条故障激励的完好性测试方法，所述方法针对北斗全链条开展压力测试，遍历北斗定位服务潜在故障场景，实现高可靠完好性能测试。
5.本发明采用以下技术方案：
6.一种基于北斗全链条故障激励的完好性测试方法，所述方法包括以下步骤：
7.步骤1，依据应用平台所规定的安全性能要求，设定应用平台的待测试完好性性能指标，将待测试故障场景初始化；选定待测试的性能指标和操作指标，设置应用平台的分布参数，确定故障场景需求；向模拟器输入北斗卫星历书文件；
8.步骤2，开展卫星-信号传输-应用平台中全链路的应用平台测试区域范围划定，所述测试区域包括多个测试地点，每个所述测试区域确定待测试地点的经度和纬度坐标；
9.步骤3，选定故障北斗卫星；根据所述待测地点的经度和纬度以及所述卫星历书文件，逐个计算待测试地点可视北斗卫星的位置，并从待测试地点所有可视北斗卫星中，通过所述模拟器选定当前历元几何构型下发生完好性事件概率最大的北斗卫星，该卫星即被选为故障北斗卫星；所述可视指可接收到北斗卫星信号；
10.步骤4，模拟北斗卫星异常信号；所述模拟器根据所述故障北斗卫星的信号，进行模拟，依据所述故障北斗卫星的测距，添加一定误差值，得到模拟的异常信号；
11.步骤5，评估应用平台的完好性性能；模拟器根据步骤1输入的北斗卫星历书，模拟北斗卫星导航系统运行过程，按照步骤3确定故障北斗卫星，按照步骤4模拟北斗卫星异常；在北斗信号异常条件下，依次记录应用平台在步骤2所划分测试地点的完好性输出，计算完
好性风险；最后与步骤1所设定的待测试完好性性能指标比较，当计算完好性风险小于待测试完好性性能指标时，证明应用平台的完好性性能合格。
12.通过设置经度和纬度的区间，确定所述待测区域；对所述待测区域进行网格点划分，确定待测试地点的经度和纬度的坐标。
13.所述待测试区域的范围划定包括以下步骤：
14.进行测试区域大小划定：设置所述待测试区域经纬度范围lon
°
e-lon
°
w、lat
°
s-lat
°
n；
15.进行测试区域内网格点划分：将所述待测试区域划分为k
°×k°
的网格，选所述网格的网格线交点为所述待测地点的经、纬度参数；
16.计算所述待测地点的总数：station＝[(lon
°
e－lon
°
w)+1]
×
[(lat
°
s－lat
°
n)+1]/(k
°×k°
)；
[0017]
其中，lat为经度、lon为纬度、k为网格点间隔。
[0018]
所述模拟器通过数学模型选定所述故障北斗卫星，所述数学模型为：
[0019][0020]
其中a表示量测域投影至定位域转换矩阵，s表示量测域投影至检测统计域转换矩阵，q表示定位误差方向，j表示卫星星号；q在a中表示矩阵的第q行，j表示a的第j列；j在s中表示矩阵的j行j列；设置当前几何构型下slopej最大的北斗卫星为故障卫星，所述故障卫星的星号为k，在所述故障卫星k上模拟北斗卫星异常信号。
[0021]
所述步骤4还包括：向所述故障卫星k，添加一定误差值，所述误差为测距误差，通过概率模型搜索获取，所述搜索的目标函数为：
[0022]
max p
nd
·
p
pf
·
p
hi
\*mergeformat(2)；
[0023]
式中，p
nd
为检测失败概率、p
pf
为定位超限概率、p
hi
为先验故障概率，所述p
hi由
初始化参数给定，p
nd
和p
pf
通过计算得到：
[0024][0025]
其中，η表示检测阈值；n表示可视卫星个数；λ表示非中心化参数；b表示测距误差；σ表示定位误差标准差；q表示标准正态分布尾部描述函数；p
ncx
为非中心卡方分布概率分布函数；q表示定位误差方向，q在a中表示矩阵的第q行。
[0026]
所述待测试性能指标包括：完好性风险指标ir
req
，告警极限l，告警时间t；所述待测试性能指标用于完好性性能评估。
[0027]
所述操作指标包括：测试时长t，测试采样间隔ts；所述性能操作指标用于规定完好性测试时长与间隔。
[0028]
所述应用平台的分布参数包括：经、纬度lat、lon和分布密度k；
[0029]
所述故障场景需求包括：待测试故障类型hi，i表示故障类型编号，故障类型包括：电离层梯度异常hi、对流层梯度异常h
t
、卫星轨道机动ho、卫星钟跳变hs及周跳hc；待测试故
障对应先验故障概率p
hi
；
[0030]
所述故障场景包括重故障并发情形。
[0031]
所述步骤5包括：
[0032]
步骤51、定义测试样本即历元数目：根据待测试完好性性能指标ir
req
，确定待测试样本量n＝p
hi
/ir
req
，单个测试地点测试样本数ns＝n/station；
[0033]
步骤52、根据测试地点位置、故障卫星号k以及测距误差取值b，组成故障场景；
[0034]
步骤53、在故障场景下，激励完好性监测算法，并统计每个测在t时间内输出的完好性参数，完好性输出参数包括：保护水平p，告警间隔s，待测定位误差e；
[0035]
步骤54、定义完好性丢失：当定位误差e大于告警限制l，但保护水平p未发出告警或告警超时，即{|e|＞l＞p∩s＞t}，则判定完好性丢失，统计完好性丢失总次数n
ir
；
[0036]
步骤55、计算完好性风险：n
ir
/n。
[0037]
本发明实施例提供的一种基于北斗全链条故障激励的完好性测试方法，该方法包括：应用平台测试地点规划操作，支撑不同地点下的应用平台完好性性能测试；故障卫星选取操作，支撑不同故障卫星下的应用平台完好性性能测试；故障信号异常模拟操作，支撑不同故障源下的应用平台完好性性能测试。通过应用平台测试地点规划、故障卫星选定以及信号异常模拟操作，严格测试北斗定位服务体系所涉及的卫星-信号传输-应用平台三类关键环节，实现可靠的完好性性能测试。
附图说明
[0038]
图1为本发明实施例提供的一种基于北斗全链条故障激励的完好性测试方法流程图。
具体实施方式
[0039]
以下结合附图对本发明的具体实施方式作出详细说明。
[0040]
参见图1，本发明实施例提供了一种基于北斗全链条故障激励的完好性测试方法，该方法可应用于在建或现有支持定位服务应用平台的完好性性能测试。具体步骤包括：
[0041]
步骤1，待测试完好性性能指标以及待测试故障场景初始化。
[0042]
步骤2，开展全链条“卫星-信号-应用平台”中，应用平台测试区域范围划定；
[0043]
步骤3，开展全链条“卫星-信号-应用平台”中，故障北斗卫星选定；
[0044]
步骤4，开展全链条“卫星-信号-应用平台”中，北斗卫星信号异常模拟；
[0045]
步骤5，评估应用平台的完好性性能。
[0046]
本发明实施例提供的针对北斗系统“卫星-信号-应用平台”全链条测试方法，该方法包括：应用平台测试地点规划操作，支撑不同地点下的应用平台完好性性能测试；故障卫星选取操作，支撑不同故障卫星下的应用平台完好性性能测试；故障信号异常模拟操作，支撑不同故障源下的应用平台完好性性能测试。通过应用平台测试地点规划、故障卫星选定以及信号异常模拟操作，严格测试北斗定位服务体系所涉及的卫星-信号传输-应用平台三类关键环节，实现可靠的完好性性能测试。
[0047]
本发明的基于北斗全链条故障激励的完好性测试方法，该方法通过模拟北斗定位服务“卫星-信号传输-应用平台”全链路所面临的潜在故障场景，激励并记录完好性监测算
法输出响应值，统计完好性性能。
[0048]
其中，涵盖“卫星-信号传输-应用平台”三个方面的全链条故障场景是由应用平台测试地点选定、北斗故障卫星选定以及北斗信号异常模拟环节构成；完好性性能评估是通过记录完好性监测算法在故障场景下的输出响应值，计算完好性丧失的次数在总测试历元数中的占比。
[0049]
具体完好性测试总体流程如下：
[0050]
首先，为保障应用平台在任意地点均满足完好性性能需求，需开展不同地点下的性能测试，因此开展测试区域选定及测试地点划分工作，该过程在步骤2中具体说明。
[0051]
其次，为测试应用平台在多类型故障场景下的完好性性能，需针对步骤2中已经划分的全部测试地点，依次开展不同故障源hi下的卫星故障及信号异常模拟，该过程在步骤3、4中具体说明。
[0052]
在预设测试时间t内，依次收集每个测试地点在故障场景下(故障场景由步骤2、3、4决确定)所输出完好性响应值，当设定区域内所有测试地点均完成测试后，汇总完好性监测算法输出响应信息用于评估完好性性能，该过程在步骤5中具体说明。
[0053]
测试如附图1所示，为本发明的一种基于北斗全链条故障模拟的完好性测试方法流程图。
[0054]
步骤1，依据应用平台所规定的安全性能要求，设定应用平台的待测试完好性性能指标，待测试完好性性能指标以及待测试故障场景设定；将待测试故障场景初始化；选定待测试的性能指标和操作指标，设置应用平台的分布参数，确定故障场景需求；向模拟器输入北斗卫星历书文件；所述应用平台是指采用北斗定位服务开展导航、定位作业的系统。
[0055]
开展完好性性能测试，需首先实现待测试完好性工作指标以及待测试故障场景初始化。涉及初始化的参数及内容包括：
[0056]
待测试完好性性能测试指标：
[0057]
1、待测试性能指标：待测试完好性性能指标irreq，告警极限l，告警时间t；上述指标用于步骤5中，完好性性能评估；
[0058]
2、操作指标：测试时长t，测试采样间隔ts；上述参数用于规定完好性测试时长与间隔；
[0059]
待测试故障场景：
[0060]
应用平台分布参数：经、纬度lat、lon，分布密度k；上述参数用于步骤2中，测试区域设定与测试地点划分；
[0061]
故障场景需求：待测试故障类型hi(i表示故障类型编号，常见故障类型包括：电离层梯度异常hi、对流层梯度异常h
t
、卫星轨道机动ho、卫星钟跳变hs及周跳hc)；待测试故障对应先验故障概率p
hi
。这里需要说明，故障场景可以设置为多重故障并发情形；上述参数故障类型与先验故障概率参数用于步骤4中，卫星信号异常模拟；卫星钟的跳变，属于卫星导航领域一种常见故障
[0062]
外部文件输入：北斗卫星历书，用于步骤3、4中卫星位置计算。
[0063]
步骤2，开展全链条“卫星-信号-应用平台”中，应用平台测试区域范围划定；
[0064]
为保障应用平台在任意地点均满足完好性性能需求，开展不同地点下的性能测试模拟。首先通过设置经、纬度区间，确定待测试区域大小；其次对测试区域开展网格点划分，
确定具体待测试地点经、纬度坐标。具体待测试区域范围划定的操作流程如下：
[0065]
1、测试区域大小划定：设置待测试区域经纬度lon
°
e～lon
°
w、lat
°
s～lat
°
n；
[0066]
2、测试区域内网格点划分：待测试区域划分为k
°×k°
的网格，选网格线交点为待测试地点点经、纬度参数；
[0067]
3、计算待测地点总数：station＝[(lon
°
e－lon
°
w)+1]
×
[(lat
°
s－lat
°
n)+1]/(k
×
k)；
[0068]
步骤2输入、输出参数归纳：
[0069]
(1)步骤2所需初始化参数：经纬度lon，lat，网格点间隔k；
[0070]
(2)步骤2输出参数：测试地点数目station，测试地点经、纬度。
[0071]
步骤3，开展全链条“卫星-信号-应用平台”中，故障北斗卫星选定
[0072]
完成参数初始化及测试地点位置划分工作后，进行全链条“卫星-信号-应用平台”中，故障卫星选定操作。以步骤2划分的测试地点经、纬度位置以及历书文件计算卫星位置为基础，通过模拟器选定当前历元几何构型下最易发生完好性事件的北斗卫星，支撑完好性压力测试开展。最易发生完好性事件的北斗卫星可通过下述数学模型判定，
[0073][0074]
其中a为量测域投影至定位域转换矩阵，s表示量测域投影至检测统计域转换矩阵，q表示定位误差方向，j表示卫星星号；q在a中表示矩阵的第q行，j表示a的第j列；j在s中表示矩阵的j行j列。设置当前几何构型下的最大slopej北斗卫星为故障卫星。步骤2完成后，可确定故障卫星号k，后续信号异常模拟将在卫星k上开展。
[0075]
步骤3输入、输出参数归纳：
[0076]
步骤3输入参数：测试地点经、纬度；
[0077]
步骤3所需初始化参数：北斗卫星历书；
[0078]
(3)步骤3输出参数：故障卫星号k，用于支撑信号异常模拟。
[0079]
步骤4，开展全链条“卫星-信号-应用平台”中，北斗信号异常模拟
[0080]
以步骤3中选定故障卫星k为基础，开展北斗卫星信号异常模拟。北斗信号异常是指向故障卫星k的测距添加一定误差值，该过程依据为美国航空无线电技术委员会rtca do-245官方文档：卫星钟、轨异常，电离层、对流层梯度异常、周跳多类型故障均可等效转换为测距误差，开展统一测试。
[0081]
故障卫星k的选定操作已在步骤2中完成，接下来需确定卫星发生故障时所产生的测距误差大小。由于测距误差的选取没有明确的范围限制，为避免测距误差选取不确定性，选择注入令应用平台最容易发生完好性事件的误差。上述最易发生完好性事件的测距误差可通过概率模型搜索获取，搜索的目标函数如下所示：
[0082]
max p
nd
·
p
pf
·
p
hi
\*mergeformat(2)
[0083]
式由检测失败概率、定位超限概率、先验故障概率(初始化参数给定)组成，
[0084][0085]
其中η表示检测阈值；n表示可视卫星个数；λ表示非中心化参数；b表示测距误差；σ表示定位误差标准差；q表示标准正态分布尾部描述函数。遍历搜索最易发生完好性事件的测距误差，并将其注入故障卫星开展完好性压力测试；q表示定位误差方向，q在a中表示矩阵的第q行。
[0086]
步骤4输入、输出参数归纳：
[0087]
(1)步骤4输入参数：测试地点经、纬度，上述参数用于计算卫星几何构型(步骤2)；卫星号k，用于确定故障卫星(步骤3)；
[0088]
(2)步骤4所需初始化参数：故障场景hi，先验故障概率p
hi
，故障先验概率，即设置故障发生的频率；两者均是用于计算所需注入测距误差；
[0089]
(3)步骤4输出参数：故障测距误差b。
[0090]
步骤5，评估应用平台的完好性性能。
[0091]
1、定义测试样本数目：根据待测试完好性性能指标ir
req
(初始化参数)，确定待测试样本量n＝(p
hi
/ir
req
)，单个测试地点测试样本数ns＝n/station；
[0092]
2、根据步骤2、3、4输出测试地点位置、故障卫星号k以及测距误差取值b，组成故障场景；
[0093]
3、在故障场景下，激励完好性监测算法，并统计每个测试地点在t时间内输出的完好性参数。所需统计完好性输出参数包括：保护水平p，告警间隔s，待测定位误差e；
[0094]
4、定义完好性丢失：当定位误差e大于告警限制l，但保护水平p未发出告警或告警超时，即{|e|＞l＞p∩s＞t}，则判定完好性丢失，统计完好性丢失总次数n
ir
；
[0095]
5、计算完好性风险：n
ir
/n。
[0096]
步骤5输入、输出参数归纳：
[0097]
(1)步骤5输入参数：测试数目station(步骤2)；
[0098]
(2)步骤5所需初始化参数：测试时长t，测试间隔ts，待测试完好性性能指标ir
req
，先验故障概率p
hi
；卫星钟的跳变，属于卫星导航领域一种常见故障
[0099]
(3)步骤5输出参数：完好性风险。
[0100]
鉴于北斗导航“卫星-信号-应用平台”全链条服务体系潜在面临多类型故障场景，同时为满足北斗导航定位服务严格可信性需求，本文提出一种方法，针对北斗系统全链条开展压力测试，遍历北斗定位服务潜在故障场景，实现高可靠完好性能测试。
[0101]
本发明实施例提供的一种基于北斗全链条故障激励的完好性测试方法，该方法包括：应用平台测试地点规划操作，支撑不同地点下的应用平台完好性性能测试；故障卫星选取操作，支撑不同故障卫星下的应用平台完好性性能测试；故障信号异常模拟操作，支撑不同故障源下的应用平台完好性性能测试。通过应用平台测试地点规划、故障卫星选定以及信号异常模拟操作，严格测试北斗定位服务体系所涉及的卫星-信号传输-应用平台三类关键环节，实现可靠的完好性性能测试。
[0102]
对于本领域技术人员而言，显然本发明实施例不限于上述示范性实施例的细节，
而且在不背离本发明实施例的精神或基本特征的情况下，能够以其他的具体形式实现本发明实施例。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明实施例的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化涵括在本发明实施例内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。此外，显然“包括”一词不排除其他单元或步骤，单数不排除复数。系统、装置或终端权利要求中陈述的多个单元、模块或装置也可以由同一个单元、模块或装置通过软件或者硬件来实现。第一，第二等词语用来表示名称，而并不表示任何特定的顺序。
[0103]
最后应说明的是，以上实施方式仅用以说明本发明实施例的技术方案而非限制，尽管参照以上较佳实施方式对本发明实施例进行了详细说明，本领域的普通技术人员应当理解，可以对本发明实施例的技术方案进行修改或等同替换都不应脱离本发明实施例的技术方案的精神和范围。