[发明专利]一种北斗卫星导航系统定位精度门限值计算方法

说明书

本发明提供一种北斗卫星导航系统定位精度门限值计算方法，包括以下步骤：根据星座构型参数、轨道周期，设定观测仰角、仿真步长以及格网点参数，假设2颗中轨道卫星故障，计算星座对所有格网点水平和垂直精度因子；统计得到水平精度因子最大值，统计该点水平精度因子值的概率密度和累计分布函数；找到累计分布函数值为可用性门限值时对应的水平精度因子值，将该水平精度因子值乘以用户等效距离误差，得到全球水平定位精度门限值；统计得到垂直精度因子最大值，统计该点垂直精度因子值的概率密度和累计分布函数；找到累计分布函数值为可用性门限值时对应的垂直精度因子值，将该垂直精度因子值乘以用户等效距离误差，得到全球垂直定位精度门限值。

技术领域

本发明属于北斗卫星导航系统定位精度计算领域，具体涉及一种北斗卫星导航系统定位精度门限值计算方法。

背景技术

北斗卫星导航系统于2012年12月，正式向亚太地区提供服务，目前正在部署全球星座，计划于2020年左右向全球提供服务。全球定位服务一个很关键的指标就是定位精度，包括水平定位精度和垂直定位精度。如何向全球用户承诺一个合理的定位精度，即定位精度门限值，是反映全球卫星导航定位水平的重要因素。目前美国的GPS和俄罗斯的GLONASS系统均已向用户提供导航定位服务，但是都没有公开其定位精度门限值计算方法。

目前，学者们研讨全球卫星导航定位系统服务性能的文献较多，例如《全球卫星导航系统性能标准及指标体系概述》(龚佩佩等，第四届中国卫星导航学术年会电子文集，2013年5月，武汉)、《格洛纳斯服务性能及其初步评估研究》(白羽等，第五届中国卫星导航学术年会电子文集，2014年5月，南京)、《GPS标准定位服务性能规范评估方法研究》(宋晓丽等，第五届中国卫星导航学术年会电子文集，2014年5月，南京)，现有文献直接给出了GPS和GLONASS的定位精度门限，没有给出具体的计算方法。文献《北斗区域系统PDOP可用性状态更新与仿真分析》(张天桥等，无线电工程，2017，47(2):52－56)分析了北斗二号区域系统的位置精度因子可用性状态，文中可用性门限采用的是星座位置精度因子门限值，这在系统设计阶段是可以接受的，但是针对北斗全球系统服务，需要在分析所有点定位精度的基础上，重点关注定位精度最差点，将定位精度最差点满足一定可用性门限值时的定位精度作为全球定位服务的精度门限值。

发明内容

针对现有技术存在的缺陷，本发明提供一种北斗卫星导航系统定位精度门限值计算方法，能够在找到全球定位最差点的基础上，通过统计精度因子值的概率密度和累计分布函数，得到满足可用性门限值时的全球定位精度门限值，可将此定位精度门限值作为向用户公开承诺的定位精度。

本发明采用的技术方案如下：

本发明提供一种北斗卫星导航系统定位精度门限值计算方法，包括以下步骤，如图1所示：

步骤1：假设2颗中地球轨道卫星故障，统计所有地面格网点水平精度因子和垂直精度因子；

北斗卫星导航系统空间星座部分由27颗中地球轨道、5颗地球静止轨道以及3颗倾斜地球同步轨道组成，星座轨道周期取7天，在7天时间内中地球轨道运行13圈，观测仰角取5°，仿真步长取60秒，格网点参数取4°×4°，假设2颗中地球轨道卫星故障，计算北斗卫星导航系统对所有地面格网点水平精度因子和垂直精度因子。27颗中地球轨道卫星中随机假设2颗卫星故障，总共有351种情况。计算所有351种情况对应的所有地面格网点水平精度因子和垂直精度因子。

步骤2：找到水平精度最差点，统计水平精度因子值的概率密度和累计分布函数；

统计得到所有地面格网点水平精度因子的最大值，该最大值对应的格网点即水平精度最差点，计算水平精度最差点在7天内的水平精度因子值，统计水平精度因子值的概率密度和累计分布函数。图2给出了水平精度因子的分布曲线。

步骤3：找到累计分布函数值为可用性门限值时对应的水平精度因子值，乘以用户等效距离误差，得到全球水平定位精度门限值：