一种北斗非机动型广播星历故障实时监测方法技术

本发明专利技术公开了一种北斗非机动型广播星历故障实时监测方法，包括选取广播星历参数外推时间，分析广播星历参数差值的统计分布特性；基于卫星位置对广播星历参数变化的灵敏度分析，得到各个广播星历参数变化对卫星位置的影响大小，对卫星位置影响较大的广播星历参数采用高阶的外推方式，对卫星位置影响较小的广播星历参数采用低阶的外推方式，根据所需导航性能需求获取满足系统要求的阈值，并将其与检验统计量比较，实现对北斗非机动型广播星历故障的实时监测。本发明专利技术对广播星历参数采用高低阶多项式混合的外推方式，在保证较低计算量的同时，可获得更高精度卫星位置偏差，有效监测北斗非机动型广播星历故障，为北斗广播星历故障监测提供参考。监测提供参考。监测提供参考。

【技术实现步骤摘要】
一种北斗非机动型广播星历故障实时监测方法


[0001]本专利技术属于导航
，具体涉及一种北斗非机动型广播星历故障实时监测方法，特别是一种基于星历参数混合外推的北斗非机动型广播星历故障实时监测方法。

技术介绍

[0002]将卫星导航系统应用于飞机精密进近，可有效提高飞行安全和效率。飞机精密进近导航的完好性和连续性是航空用户最关注的性能需求，其中卫星广播星历故障是影响飞机精密进近的主要风险源之一。为保障与生命安全相关用户的定位完好性，对广播星历故障的监测是必不可少的。广播星历故障可分为机动性广播星历故障和非机动型广播星历故障。由于卫星机动发生的概率很小，所以非机动型广播星历故障是影响定位完好性的主要因素。
[0003]全球定位系统监测非机动型广播星历故障的方法是基于卫星运动的轨道重合期，采用零阶保持的方式将上一轨道周期的广播星历参数外推到当前时刻，计算卫星位置偏差，构建监测模型。由于每个广播星历参数对卫星位置的影响大小不同，对所有参数采用同一外推方式是不合理的。除此之外，北斗混合星座的构型及卫星运动的轨道重合周期和全球定位系统存在显著差异，所以传统非机动型广播星历故障监测方法无法直接应用于北斗系统。综上所述，设计一种适用于北斗多类型卫星的非机动型广播星历故障实时监测方法具有重要意义。

技术实现思路

[0004]针对上述现有技术，本专利技术要解决的技术问题是提供一种基于星历参数混合外推的北斗非机动型广播星历故障实时监测方法，在保证较低计算量的同时得到更高精度的广播星历参数预测值，进一步获得更高精度的卫星位置偏差。
[0005]为解决上述技术问题，本专利技术的一种北斗非机动型广播星历故障实时监测方法，包括以下步骤：
[0006]步骤1：基于卫星的轨道运动重合周期及位置偏差选取广播星历参数外推时间；
[0007]步骤2：计算广播星历参数外推时间间隔的星历参数差值，分析其统计分布特性；
[0008]步骤3：卫星位置偏差对广播星历参数变化的灵敏度分析；
[0009]步骤4：根据灵敏度分析中广播星历参数对卫星位置的影响大小，对广播星历参数进行影响等级分类；
[0010]步骤5：对不同影响等级中的广播星历参数选择不同阶数的多项式外推；
[0011]步骤6：基于卫星位置偏差的统计分布特性，构建检验统计量和最小可检测误差；
[0012]步骤7：根据系统的连续性要求和完好性风险要求获取检验统计量阈值，通过阈值和检验统计量比较的方式，实现北斗非机动型广播星历故障实时监测，具体为：若计算得到的检验统计量小于检验统计量阈值，则满足系统的连续性要求，否则发生误警；在满足连续性要求的前提下，若计算的最小可检测误差小于其给定检测门限，则满足系统的完好性风
险要求，否则发生漏检。
[0013]本专利技术还包括：
[0014]1.步骤1中广播星历参数外推时间为1天。
[0015]2.步骤2中计算广播星历参数外推时间间隔的星历参数差值，分析其统计分布特性具体为：
[0016]设当前k时刻，相对于星历参考时间t
oe
，广播星历参数为P
k
，相邻日时刻的历史广播星历参数为P
k
‑1，计算广播星历参数偏差可表示为：
[0017]ΔP＝P
k
‑
P
k
‑1[0018]其中，P代表由15个广播星历参数组成的矩阵；
[0019]通过计算相邻两天的广播星历参数差值得到其统计分布特性。
[0020]3.步骤3中卫星位置偏差对广播星历参数变化的灵敏度分析具体为：
[0021]通过使用北斗接口控制文件
‑
2013中定义的卫星位置非线性函数计算出卫星位置，卫星位置计算模型具体为：
[0022][0023]其中，x、y、z为卫星三轴位置的计算函数，t为当前时刻，相对于t
oe
，p
i
表示不同广播星历参数，i＝1，2，
…
，15；
[0024]卫星位置对广播星历参数变化的灵敏度系数通过卫星位置函数分别对每个广播星历参数的微分计算，具体为：
[0025][0026]其中，δx(t)、δy(t)、δz(t)分别为卫星三轴位置误差，δp表示星历参数的变化量；
[0027]简写为：
[0028]δr(t)＝AδP
[0029]其中，δr(t)表示t时刻卫星位置误差；
[0030]则广播星历参数变化对卫星位置误差的影响为：
[0031][0032]其中，σ(δr
i
)是卫星位置误差标准差，σ(δp
i
)是相邻日广播星历参数差值的标准差；
[0033]根据上述公式可以得到每个广播星历参数变化分别对北斗地球静止轨道卫星、倾斜地球同步轨道卫星和中圆地球轨道卫星三种轨道卫星位置的影响大小。
[0034]4.步骤4中根据灵敏度分析中广播星历参数对卫星位置的影响大小，对广播星历参数进行影响等级分类具体为：根据步骤3中的灵敏度分析结果，将对北斗地球静止轨道卫星、倾斜地球同步轨道卫星和中圆地球轨道卫星三种轨道卫星位置影响程度不同的广播星历参数分为三类：将对卫星位置影响误差超过1000m的广播星历参数归类为影响等级一，将影响卫星位置误差低于100m的广播星历参数归类为影响等级三，将其余广播星历参数归类为影响等级二。
[0035]5.步骤5中对不同影响等级中的广播星历参数选择不同阶数的多项式外推具体为：
[0036]根据步骤4中的广播星历参数影响等级分类，对影响等级三范围内的星历参数采用一阶的外推方式，具体为：
[0037]P
k
＝P
k
‑1+(P
k
‑1‑
P
k
‑2)
[0038]其中，P
k
表示当前广播星历的k时刻，P
k
‑
i
表示第i天前广播星历的k时刻；
[0039]对影响等级一中内的参数采用三阶的外推方式，具体为：
[0040]P
k
＝P
k
‑1+(P
k
‑1‑
P
k
‑2)+[(P
k
‑1‑
P
k
‑2)
‑
(P
k
‑2‑
P
k
‑3)]+{[(P
k
‑1‑
P
k
‑2)
‑
(P
k
‑2‑
P
k
‑3)]‑
[(P
k
‑2‑
P
k
‑3)
‑
(P
k
‑3‑
P
k
‑4)]}
[0041]对影响等级二中的参数采用二阶的外本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种北斗非机动型广播星历故障实时监测方法，其特征在于，包括以下步骤：步骤1：基于卫星的轨道运动重合周期及位置偏差选取广播星历参数外推时间；步骤2：计算广播星历参数外推时间间隔的星历参数差值，分析其统计分布特性；步骤3：卫星位置偏差对广播星历参数变化的灵敏度分析；步骤4：根据灵敏度分析中广播星历参数对卫星位置的影响大小，对广播星历参数进行影响等级分类；步骤5：对不同影响等级中的广播星历参数选择不同阶数的多项式外推；步骤6：基于卫星位置偏差的统计分布特性，构建检验统计量和最小可检测误差；步骤7：根据系统的连续性要求和完好性风险要求获取检验统计量阈值，通过阈值和检验统计量比较的方式，实现北斗非机动型广播星历故障实时监测，具体为：若计算得到的检验统计量小于检验统计量阈值，则满足系统的连续性要求，否则发生误警；在满足连续性要求的前提下，若计算的最小可检测误差小于其给定检测门限，则满足系统的完好性风险要求，否则发生漏检。2.根据权利要求1所述的一种北斗非机动型广播星历故障实时监测方法，其特征在于：步骤1所述广播星历参数外推时间为1天。3.根据权利要求1所述的一种北斗非机动型广播星历故障实时监测方法，其特征在于：步骤2所述计算广播星历参数外推时间间隔的星历参数差值，分析其统计分布特性具体为：设当前k时刻，相对于星历参考时间t
oe
，广播星历参数为P
k
，相邻日时刻的历史广播星历参数为P
k
‑1，计算广播星历参数偏差可表示为：ΔP＝P
k
‑
P
k
‑1其中，P代表由15个广播星历参数组成的矩阵；通过计算相邻两天的广播星历参数差值得到其统计分布特性。4.根据权利要求1所述的一种北斗非机动型广播星历故障实时监测方法，其特征在于：步骤3所述卫星位置偏差对广播星历参数变化的灵敏度分析具体为：通过使用北斗接口控制文件
‑
2013中定义的卫星位置非线性函数计算出卫星位置，卫星位置计算模型具体为：其中，x、y、z为卫星三轴位置的计算函数，t为当前时刻，相对于t
oe
，p
i
表示不同广播星历参数，i＝1，2，
…
，15；卫星位置对广播星历参数变化的灵敏度系数通过卫星位置函数分别对每个广播星历参数的微分计算，具体为：
其中，δx(t)、δy(t)、δz(t)分别为卫星三轴位置误差，δp表示星历参数的变化量；简写为：δr(t)＝AδP其中，δr(t)表示t时刻卫星位置误差；则广播星历参数变化对卫星位置误差的影响为：其中，σ(δr
i
)是卫星位置误差标准差，σ(δp
i
)是相邻日广播星历参数差值的标准差；根据上述公式可以得到每个广播星历参数变化分别对北斗地球静止轨道卫星、倾斜地球同步轨道卫星和中圆地球轨道卫星三种轨道卫星位置的影响大小。5.根据权利要求1所述的一种北斗非机动型广播星历故障实时监测方法，其特征在于：步骤4所述根据灵敏度分析中广播星历参数对卫星位置的影响大小，对广播星历参数进行影响等级分类具体为：根据步骤3中的灵敏度分析结果...

【专利技术属性】
技术研发人员：李亮，刘翼，翁铭泽，杨福鑫，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：