CN114280633A - 一种非差非组合精密单点定位完好性监测方法 - Google Patents
一种非差非组合精密单点定位完好性监测方法 Download PDFInfo
- Publication number
- CN114280633A CN114280633A CN202111628146.7A CN202111628146A CN114280633A CN 114280633 A CN114280633 A CN 114280633A CN 202111628146 A CN202111628146 A CN 202111628146A CN 114280633 A CN114280633 A CN 114280633A
- Authority
- CN
- China
- Prior art keywords
- frequency
- satellite
- protection level
- fault
- current epoch
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000012544 monitoring process Methods 0.000 title claims abstract description 20
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 238000000926 separation method Methods 0.000 claims abstract description 13
- 230000001419 dependent effect Effects 0.000 claims description 15
- 230000009977 dual effect Effects 0.000 claims description 15
- 230000010339 dilation Effects 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 11
- 238000012937 correction Methods 0.000 claims description 7
- 230000001934 delay Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000003321 amplification Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 4
- 239000005433 ionosphere Substances 0.000 claims description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 6
- 230000006870 function Effects 0.000 description 29
- 230000003416 augmentation Effects 0.000 description 4
- 230000009191 jumping Effects 0.000 description 4
- 238000004590 computer program Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种非差非组合精密单点定位完好性监测方法。所述方法包括:建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型;基于所述非差非组合PPP函数模型,输入当前历元处的可用观测量,所述可用观测量包含单频点观测量和双频点观测量;基于多重假设解分离理论,列出所有故障子集;针对每个故障子集,进行所述故障子集的时间更新和测量值更新;根据更新后的测量值,计算每个故障子集的保护水平,进而得到所述当前历元的整体保护水平;根据所述当前历元的整体保护水平进行可用性判断。本发明有效抑制了显著多路径效应对PPP函数模型完好性监测算法的影响,对海上常见危险源有很强的适应性,显著提升了PPP位置服务的可靠性。
Description
技术领域
本发明涉及卫星导航定位技术领域,具体涉及一种非差非组合精密单点定位完好性监测方法。
背景技术
精密单点定位(Precise Point Positioning,PPP)技术具有作业灵活、无需基站辅助的优势,可为海洋开发等多种场景提供分米级甚至厘米级的高精度位置服务。PPP的定位性能虽然可以满足大部分场景的要求,但在卫星信号接收受限的海上作业环境下,如钻井平台上的建筑物遮挡,船舶远洋航行受船体及冰山遮挡,以及一些海上的极端天气影响,会导致部分频点丢失以及显著多路径频发的现象,其位置服务的可靠性难以得到保障。
现有技术中,为保障PPP位置服务的可靠性,可以通过将RAIM(ReceiverAutonomous Integrity Monitoring)引入PPP建立PPP-RAIM算法。虽然这一算法进行故障检测时能够一定程度上提高可靠性,但与PPP模型相关性较弱。为进一步提升与PPP模型契合度,基于SBAS校正信息的双频PPP算法,结合ARAIM(Advanced RAIM)完好性算法,可产生10米以内的保护水平,但构建的保护水平(Protection Level,PL)过于保守,超出高精度位置服务应用允许的误差界限范围。基于扩展卡尔曼滤波估计的PPP完好性监测算法可以实现米级至分米级高完好性保证的误差界限,但使用环境受限于开阔收星条件下的航空环境,在卫星受限的海洋作业环境中收到一定的限制。因此,面向卫星受限及显著多径频发的海上作业环境下PPP可靠性保障问题亟待解决。
发明内容
鉴于上述问题,提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的非差非组合精密单点定位完好性监测方法。
根据本发明的一个方面,提供了一种非差非组合精密单点定位完好性监测方法,包括:
建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型;
基于所述非差非组合PPP函数模型,输入当前历元处的可用观测量,所述可用观测量包含单频点观测量和双频点观测量;
基于多重假设解分离理论,列出所有故障子集;
针对每个故障子集,进行所述故障子集的时间更新和测量值更新;
根据更新后的测量值,计算每个故障子集的保护水平,进而得到所述当前历元的整体保护水平;
根据所述当前历元的整体保护水平进行可用性判断。
在一种可选的方式中,所述建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型进一步包括:
根据GNSS观测量与待估参数,构建原始伪距与载波相位基本观测方程;
引入精密卫星轨道和种差改正,对所述原始伪距与载波相位基本观测方程进行参数规整,建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型。
在一种可选的方式中,所述原始伪距与载波相位基本观测方程为:
其中,为原始伪距;为载波相位;上标s和Q分别表示卫星的PRN号和卫星系统;下标r和j分别表示接收机ID和观测值频段号;表示卫星与测站之间的几何距离;表示频率fj对应的载波波长;c表示真空中的光速;dtr和dts,Q分别表示接收机和卫星钟差;表示与卫星高度角有关的湿投影函数;Zw表示测站的天顶湿延迟;表示频率相关的电离层延迟方法因子;和分别表示与频率相关的接收机和卫星端的非校正伪距硬件延迟;和分别表示频率相关的接收机和卫星端的非校正相位硬件延迟;表示载波相位整周模糊度;为伪距噪声;为载波相位的噪声;为f1频段上对应的斜向电离层延迟误差。
在一种可选的方式中,所述非差非组合PPP函数模型为:
其中,下标j表示观测值频段号,取值为1或者2;和分别表示伪距和载波的观测量减去计算量;为方向余弦;x为相对于初始坐标的三维位置增量;为频率相关放大因子;为频率相关的卫星端差分码偏差;为卫星UCD的无电离层组合,包含于外部提供的精密卫星钟差改正;为频率fj对应的载波波长;c表示真空中的光速;dtr和为接收机和卫星钟差;为与卫星高度角有关的湿投影函数;Zw为测站的天顶湿延迟;是频率相关的电离层延迟方法因子和为频率相关的接收机和卫星端的非校正相位硬件延迟;为载波相位整周模糊度;为伪距的噪声;为载波相位的噪声;为f1频段上对应的斜向电离层延迟误差。
在一种可选的方式中,在所述基于多重假设解分离理论,列出所有故障子集之后,所述方法还包括:
进行状态更新;所述状态包含:误差、对流层湿延迟、电离层消息和载波相位模糊状态;
其中,利用对偶膨胀算法对多路径误差进行膨胀处理,以形成多路径误差的高斯边界。
在一种可选的方式中,所述对偶膨胀算法包含一对CDF高斯边界GL,GR,分别在实际误差分布的左侧和右侧,具体描述为:
其中,GL为基于对偶膨胀算法建立的基于多路径误差的高斯左边界;GR为基于对偶膨胀算法建立的基于多路径误差的高斯右边界;Ga为实际多径误差CDF分布;x表示多径误差。
在一种可选的方式中,所述时间更新进一步包括:
对故障子集状态和协方差矩阵分别进行从上一个历元到当前历元的传播。
在一种可选的方式中,所述测量值更新进一步包括:
执行误差处理,以产生预测的伪距和载波相位观测值。
在一种可选的方式中,所述整体保护水平PL为:
其中,PLi和Ti分别为故障子集i的保护水平和检测门限;N为故障子集数量;Pap(Hi)为故障子集i所对应的先验故障概率;Q-1为正态分布的逆;σ(i)为故障子集i的标准差;PHMI表示危险误导信息概率。
在一种可选的方式中,所述根据所述当前历元的整体保护水平进行可用性判断进一步包括:
基于所述整体保护水平PL和所需导航性能确立的告警极限AL进行可用性判断;其中,当PL≥AL时,表示当前历元不可用;当PL<AL时,表示当前历元可用。
根据本发明提供的方案,建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型;基于所述非差非组合PPP函数模型,输入当前历元处的可用观测量,所述可用观测量包含单频点观测量和双频点观测量;基于多重假设解分离理论,列出所有故障子集;针对每个故障子集,进行所述故障子集的时间更新和测量值更新;根据更新后的测量值,计算每个故障子集的保护水平,进而得到所述当前历元的整体保护水平;根据所述当前历元的整体保护水平进行可用性判断。本发明有效抑制了显著多路径效应对PPP函数模型完好性监测算法的影响,对海上常见危险源有很强的适应性,显著提升了PPP位置服务的可靠性。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂,以下特举本发明的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1示出了根据本发明一个实施例的非差非组合精密单点定位完好性监测方法的流程示意图;
图2示出了根据本发明另一个实施例的非差非组合精密单点定位完好性监测方法的流程示意图。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
在实施本发明实施例之前,对下文中涉及的技术术语在此进行统一解释:
GNSS:全球导航卫星系统(Global Navigation Satellite System,GNSS),泛指所有的卫星导航系统,包括全球的、区域的和增强的,如美国的GPS、俄罗斯的Glonass、欧洲的Galileo、中国的北斗卫星导航系统,以及相关的增强系统,如美国的WAAS(广域增强系统)、欧洲的EGNOS(欧洲静地导航重叠系统)和日本的MSAS(多功能运输卫星增强系统)等,还涵盖在建和以后要建设的其他卫星导航系统。
PPP:精密单点定位(Precise Point Positioning,PPP)指的是利用全球若干地面跟踪站的GPS观测数据计算出的精密卫星轨道和卫星钟差,对单台GPS接收机所采集的相位和伪距观测值进行定位解算。
非差非组合模型:非差非组合模型是一种统一的GNSS数据处理模型,任何GNSS系统任何频率的观测值均能按照这种模型建立观测方程。
非差非组合PPP模型:非差非组合PPP模型使用原始伪距和载波相位观测值并参数化倾斜方向电离层延迟,不会放大观测值噪声,利于单频、双频、多频(三频及以上)模型的统一表达。
图1示出了根据本发明一个实施例的非差非组合精密单点定位完好性监测方法的流程示意图。本方法采用非差非组合PPP函数模型,避免了线性组合导致的噪声放大,同时在处理多频观测信息时更具有弹性,并在此基础上,利用多重假设解分离理论建立完好性监测算法,在最大限度的保留原始观测信息的同时,基于该模型的完好性算法可以获得更平滑的协方差和保护水平。以北斗系统为例,具体地,如图1所示,该方法包括如下步骤:
步骤S101,建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型。
本实施例中,构建基于GNSS观测量单双频点混合的非差非组合PPP函数模型,不仅可支持双频点观测量,而且还支持在部分频点缺失后的单频观测量,并可最大化利用观测信息,以提升观测冗余。
建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型,可执行包括但不限于以下处理步骤:
步骤一,根据GNSS观测量与待估参数之间的关系,将原始伪距和载波相位结合而构建原始伪距与载波相位基本观测方程。
步骤二,在构建的原始伪距与载波相位基本观测方程的基础上,引入精密卫星轨道和钟差改正,对原始观测方程进行参数规整。
步骤三,在所述原始伪距与载波相位基本观测方程参数规整的同时,针对部分频点缺失导致观测冗余降低情况,建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型。
具体地,步骤一所述原始伪距与载波相位基本观测方程为:
其中,为原始伪距;为载波相位;上标s和Q分别表示卫星的PRN号和卫星系统;下标r和j分别表示接收机ID和观测值频段号;表示卫星与测站之间的几何距离;表示频率fj对应的载波波长;c表示真空中的光速;dtr和dts,Q分别表示接收机和卫星钟差;表示与卫星高度角有关的湿投影函数;Zw表示测站的天顶湿延迟;表示频率相关的电离层延迟方法因子;和分别表示与频率相关的接收机和卫星端的非校正伪距硬件延迟;和分别表示频率相关的接收机和卫星端的非校正相位硬件延迟;表示载波相位整周模糊度;为伪距噪声;为载波相位的噪声;为f1频段上对应的斜向电离层延迟误差。
具体地,步骤三所述基于GNSS观测量单双频点混合的非差非组合PPP函数模型为:
其中,下标j表示观测值频段号,取值为1或者2;和分别表示伪距和载波的观测量减去计算量;为方向余弦;x为相对于初始坐标的三维位置增量;为频率相关放大因子;为频率相关的卫星端差分码偏差;为卫星UCD的无电离层组合,包含于外部提供的精密卫星钟差改正;为频率fj对应的载波波长;c表示真空中的光速;dtr和为接收机和卫星钟差;为与卫星高度角有关的湿投影函数;Zw为测站的天顶湿延迟;是频率相关的电离层延迟方法因子和为频率相关的接收机和卫星端的非校正相位硬件延迟;为载波相位整周模糊度;为伪距的噪声;为载波相位的噪声;为f1频段上对应的斜向电离层延迟误差。
步骤S102,基于非差非组合PPP函数模型,输入当前历元处的可用观测量,可用观测量包含单频点观测量和双频点观测量。
构建基于单双频点混合的非差非组合PPP函数模型,不仅支持双频点观测量,还支持在部分频点缺失后的单频观测量,以最大化利用GNSS观测信息和提升观测冗余。
为计算可用观测量的整体保护水平PL,并与告警极限AL进行可用性判断,首先需要采集所有的观测历元处的可用观测量。
具体地,基于所构建的单双频点混合的非差非组合PPP函数模型,输入当前历元处的可用观测量,其中,可用观测量包含单频点观测量和双频点观测量,在执行完毕后续的步骤后,并继续依次遍历所有观测历元处的可用观测量。
步骤S103,基于多重假设解分离理论,列出所有故障子集。
在构建所述单双频点混合的非差非组合PPP函数模型的同时,基于多重假设解分离理论,列出当前历元处的所有子集,其中,每个子集对应于一个故障模式。
具体地,包括以下步骤:对每个子集建立索引,其中,索引为0的对应于全视图模式;依次遍历每个子集;如果添加了新的观测值,则更新子集和新状态,新状态包括误差、对流层湿延迟、电离层消息和载波相位模糊状态。
可选地,同时利用对偶膨胀算法对多路径误差进行膨胀处理,以形成多路径误差的高斯边界,以补偿因多径误差畸变带来的不确定性影响,并建立适应复杂海洋作业环境下的威胁模型,以保证在显著多径效应引起的较大测距误差环境下误差界限的完好性。
步骤S104,针对每个故障子集,进行故障子集的时间更新和测量值更新。
所述时间更新是指对故障子集状态和协方差矩阵分别进行从上一个历元到当前历元的传播,可以选择对每个子集状态和协方差分别进行矩阵乘法操作,该操作具有较低的计算复杂度。
可选地,在准备测量值更新时,执行多个误差处理,以产生预测的伪距和载波相位观测值,从而进一步建立适应卫星信号接收受限环境下可信的正态模型及威胁模型,进一步地,为了排查异常值,测量值更新步骤还包括检查测量残差。
具体的,依次遍历当前历元处的每个故障子集,并分别更新每个故障子集的时间和测量值。
步骤S105,根据更新后的测量值,计算每个故障子集的保护水平,进而得到当前历元的整体保护水平。
为对当前历元的整体保护水平进行可用性判断,根据更新后的测量值,计算当前历元的每个故障子集的保护水平,并保留所有故障子集保护水平的最大值,进而得到所述当前历元的整体保护水平。
步骤S106,根据当前历元的整体保护水平进行可用性判断。
基于当前历元的整体保护水平PL和所需导航性能确立的告警极限AL进行可用性判断;其中,当PL≥AL时,表示当前历元不可用,同时,将告警标识设置在当前历元处;当PL<AL时,表示当前历元可用。
本发明上述实施例提供的方案,建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型;基于所述非差非组合PPP函数模型,输入当前历元处的可用观测量,所述可用观测量包含单频点观测量和双频点观测量;基于多重假设解分离理论,列出所有故障子集;针对每个故障子集,进行所述故障子集的时间更新和测量值更新;根据更新后的测量值,计算每个故障子集的保护水平,进而得到所述当前历元的整体保护水平;根据所述当前历元的整体保护水平进行可用性判断。本发明有效抑制了显著多路径效应对PPP函数模型完好性监测算法的影响,对海上常见危险源有很强的适应性,显著提升了PPP位置服务的可靠性。
图2示出了根据本发明另一个实施例的非差非组合精密单点定位完好性监测方法的流程示意图。本方法在基于非差非组合PPP函数模型,依次遍历所有历元处的可用观测量,并对当前历元处的可用观测量,基于多重假设解分离理论,列出所有故障子集,并进行当前历元的每个故障子集进行状态更新,同时利用对偶膨胀算法(Paired Overbounding,PB)对多路径误差进行膨胀处理,以形成多路径误差的高斯边界,根据更新后的测量值,计算每个故障子集的保护水平,进而得到当前历元的整体保护水平,根据当前历元的整体保护水平进行可用性判断,以北斗系统为例,具体地,如图2所示,该方法包括以下步骤:
步骤S201,基于非差非组合PPP函数模型,输入当前历元处的可用观测量。
基于单双频点混合的非差非组合PPP函数模型,输入当前历元处的可用观测量,其中,可用观测量包含单频点观测量和双频点观测量,在执行完毕后续的步骤后,并继续依次遍历所有观测历元处的可用观测量。
步骤S202,基于多重假设解分离理论,列出所有故障子集。
基于多重假设解分离理论,列出当前历元的所有故障子集,总数为N,其中,每个子集对应于一个故障模式。对每个子集建立索引i,其中,子集索引i为0对应于全视图模式。
步骤S203,进行状态更新。
进行状态更新,其中,如果添加了新的观测值,则更新每个子集和状态,同时利用对偶膨胀算法对多路径误差进行膨胀处理,以形成多路径误差的高斯边界。
所述状态包含:误差、对流层湿延迟、电离层消息和载波相位模糊状态。
进一步的,利用对偶膨胀算法对多路径误差进行膨胀处理,以形成多路径误差的高斯边界,以补偿因多径误差畸变带来的不确定性影响。
所述对偶膨胀算法包含一对CDF高斯边界GL,GR,分别在实际误差分布的左侧和右侧,具体描述为:
其中,GL为基于对偶膨胀算法建立的基于多路径误差的高斯左边界;GR为基于对偶膨胀算法建立的基于多路径误差的高斯右边界;Ga为实际多径误差CDF分布;x表示多径误差。
步骤S204,判断当前子集i是否小于总故障子集数N。
如果i小于总故障子集数N,则跳转到步骤S205,更新当前子集的时间和测量值;否则,跳转到步骤S207。
步骤S205,更新当前子集的时间和测量值。
依次遍历当前历元处的每个故障子集,并分别更新每个故障子集的时间和测量值。其中,时间更新是指状态和协方差矩阵从上一个时间历元到当前历元的传播,时间更新是对每个故障子集状态和协方差分别进行。
所述测量值更新进一步包括:
执行误差处理,以产生预测的伪距和载波相位观测值。
步骤S206,根据更新后的测量值,计算当前子集的保护水平。
根据更新后的测量值,计算当前子集的保护水平,并跳转到步骤S204,直至遍历完所有当前历元的每个子集。
步骤S207,计算当前历元的整体保护水平PL。
根据当前历元的每个故障子集的保护水平,并保留所有故障子集保护水平的最大值,进而得到所述当前历元的整体保护水平。
具体地,所述整体保护水平PL为:
其中,PLi和Ti分别为故障子集i的保护水平和检测门限;N为故障子集数量;Pap(Hi)为故障子集i所对应的先验故障概率;Q-1为正态分布的逆;σ(i)为故障子集i的标准差;PHMI表示危险误导信息概率。
步骤S208,判断PL是否超过告警极限AL。
判断当前历元的整体保护水平PL是否小于告警极限AL,如果PL超过了AL的约束,则表示当前历元不可用,并跳转到步骤S209;如果PL未超过AL的约束,则表示当前历元可用,并跳转到步骤S210。
步骤S209,在当前历元处设置告警标志。
在当前历元处设置告警标志,继续执行步骤S210,即进行历元更新并继续遍历下一个历元。
步骤S210,进行当前历元更新并继续遍历下一个历元。
继续遍历下一个历元,检测新的测量值,跳转到步骤S201,直至遍历完所有的历元。
本发明上述实施例提供的方案,利用对偶膨胀方法确立多径误差高斯边界,以补偿因多径误差畸变带来的不确定性影响,建立适应复杂海洋作业环境下的威胁模型,以保证在显著多径效应引起的较大测距误差环境下误差界限的完好性。并在此基础上,利用多重假设解分离理论建立完好性监测算法,在最大限度的保留原始观测信息的同时,基于该模型的完好性算法可以获得更平滑的协方差和保护水平,有效抑制了显著多路径效应对PPP函数模型完好性监测算法的影响,对海上常见危险源有很强的适应性,显著提升了PPP位置服务的可靠性。
在此提供的算法或显示不与任何特定计算机、虚拟系统或者其它设备固有相关。各种通用系统也可以与基于在此的示教一起使用。根据上面的描述,构造这类系统所要求的结构是显而易见的。此外,本发明实施例也不针对任何特定编程语言。应当明白,可以利用各种编程语言实现在此描述的本发明的内容,并且上面对特定语言所做的描述是为了披露本发明的最佳实施方式。
在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
类似地,应当理解,为了精简本发明并帮助理解各个发明方面中的一个或多个,在上面对本发明的示例性实施例的描述中,本发明实施例的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该公开的方法解释成反映如下意图:即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说,如下面的权利要求书所反映的那样,发明方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本发明的单独实施例。
本领域那些技术人员可以理解,可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件,以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外,可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述,本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外,本领域的技术人员能够理解,尽管在此的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如,在下面的权利要求书中,所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
本发明的各个部件实施例可以以硬件实现,或者以在一个或者多个处理器上运行的软件模块实现,或者以它们的组合实现。本领域的技术人员应当理解,可以在实践中使用微处理器或者数字信号处理器(DSP)来实现根据本发明实施例的一些或者全部部件的一些或者全部功能。本发明还可以实现为用于执行这里所描述的方法的一部分或者全部的设备或者装置程序(例如,计算机程序和计算机程序产品)。这样的实现本发明的程序可以存储在计算机可读介质上,或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到,或者在载体信号上提供,或者以任何其他形式提供。
应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制,并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中,不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中,这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。上述实施例中的步骤,除有特殊说明外,不应理解为对执行顺序的限定。
Claims (10)
1.一种非差非组合精密单点定位完好性监测方法,其特征在于,所述方法包括:
建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型;
基于所述非差非组合PPP函数模型,输入当前历元处的可用观测量,所述可用观测量包含单频点观测量和双频点观测量;
基于多重假设解分离理论,列出所有故障子集;
针对每个故障子集,进行所述故障子集的时间更新和测量值更新;
根据更新后的测量值,计算每个故障子集的保护水平,进而得到所述当前历元的整体保护水平;
根据所述当前历元的整体保护水平进行可用性判断。
2.根据权利要求1所述的方法,其特征在于,所述建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型进一步包括:
根据GNSS观测量与待估参数,构建原始伪距与载波相位基本观测方程;
引入精密卫星轨道和种差改正,对所述原始伪距与载波相位基本观测方程进行参数规整,建立基于GNSS观测量单双频点混合的非差非组合PPP函数模型。
3.根据权利要求2所述的方法,其特征在于,所述原始伪距与载波相位基本观测方程为:
4.根据权利要求2所述的方法,其特征在于,所述非差非组合PPP函数模型为:
5.根据权利要求1所述的方法,其特征在于,在所述基于多重假设解分离理论,列出所有故障子集之后,所述方法还包括:
进行状态更新;所述状态包含:误差、对流层湿延迟、电离层消息和载波相位模糊状态;
其中,利用对偶膨胀算法对多路径误差进行膨胀处理,以形成多路径误差的高斯边界。
7.根据权利要求1所述的方法,其特征在于,所述时间更新进一步包括:
对故障子集状态和协方差矩阵分别进行从上一个历元到当前历元的传播。
8.根据权利要求1所述的方法,其特征在于,所述测量值更新进一步包括:
执行误差处理,以产生预测的伪距和载波相位观测值。
10.根据权利要求1所述的方法,其特征在于,所述根据所述当前历元的整体保护水平进行可用性判断进一步包括:
基于所述整体保护水平PL和所需导航性能确立的告警极限AL进行可用性判断;其中,当PL≥AL时,表示当前历元不可用;当PL<AL时,表示当前历元可用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111628146.7A CN114280633B (zh) | 2021-12-28 | 2021-12-28 | 一种非差非组合精密单点定位完好性监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111628146.7A CN114280633B (zh) | 2021-12-28 | 2021-12-28 | 一种非差非组合精密单点定位完好性监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114280633A true CN114280633A (zh) | 2022-04-05 |
CN114280633B CN114280633B (zh) | 2024-05-17 |
Family
ID=80877136
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111628146.7A Active CN114280633B (zh) | 2021-12-28 | 2021-12-28 | 一种非差非组合精密单点定位完好性监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114280633B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114460615A (zh) * | 2022-04-14 | 2022-05-10 | 中国测绘科学研究院 | 一种附加虚拟观测值的北斗三新频点定位方法及系统 |
CN116859417A (zh) * | 2023-07-07 | 2023-10-10 | 哈尔滨工程大学 | 用于北斗ppp-rtk/mems的完好性监测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101843004B1 (ko) * | 2017-11-15 | 2018-03-29 | 한국 천문 연구원 | 다중 위성 항법 시스템의 시스템 바이어스를 이용한 통합 측위 장치 및 그 방법 |
CN108919634A (zh) * | 2018-08-13 | 2018-11-30 | 中国科学院国家授时中心 | 一种北斗三频非差非组合观测值时间传递系统及方法 |
CN111025346A (zh) * | 2019-11-18 | 2020-04-17 | 广州南方卫星导航仪器有限公司 | 一种快速估计gnss精密卫星钟差的方法及存储介质 |
CN111273327A (zh) * | 2020-03-20 | 2020-06-12 | 中国人民解放军61081部队 | 一种基于组合和非组合混合观测模型的精密单点定位方法 |
CN111427068A (zh) * | 2020-03-21 | 2020-07-17 | 哈尔滨工程大学 | 一种动对动平台局域增强卫星a类星历故障完好性监测方法 |
-
2021
- 2021-12-28 CN CN202111628146.7A patent/CN114280633B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101843004B1 (ko) * | 2017-11-15 | 2018-03-29 | 한국 천문 연구원 | 다중 위성 항법 시스템의 시스템 바이어스를 이용한 통합 측위 장치 및 그 방법 |
CN108919634A (zh) * | 2018-08-13 | 2018-11-30 | 中国科学院国家授时中心 | 一种北斗三频非差非组合观测值时间传递系统及方法 |
CN111025346A (zh) * | 2019-11-18 | 2020-04-17 | 广州南方卫星导航仪器有限公司 | 一种快速估计gnss精密卫星钟差的方法及存储介质 |
CN111273327A (zh) * | 2020-03-20 | 2020-06-12 | 中国人民解放军61081部队 | 一种基于组合和非组合混合观测模型的精密单点定位方法 |
CN111427068A (zh) * | 2020-03-21 | 2020-07-17 | 哈尔滨工程大学 | 一种动对动平台局域增强卫星a类星历故障完好性监测方法 |
Non-Patent Citations (4)
Title |
---|
KAZUMA GUNNING 等: "Design and evaluation of integrity algorithms for PPP in kinematic applications", PROCEEDINGS OF THE 31ST INTERNATIONAL TECHNICAL MEETING OF THE SATELLITE DIVISION OF THE INSTITUTE OF NAVIGATION, 31 December 2018 (2018-12-31) * |
伍维甲 等: "GNSS双频导航信号用户完好性监测机制研究", 航空学报, no. 12, 30 August 2012 (2012-08-30) * |
周锋: "多系统GNSS非差非组合精密单点定位相关理论和方法研究", 测绘学报, no. 07, 15 July 2020 (2020-07-15) * |
王式太: "基于M估计的GNSS接收机自主完好性监测研究", 中国博士学位论文全文数据库基础科学辑, no. 09, 15 September 2017 (2017-09-15) * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114460615A (zh) * | 2022-04-14 | 2022-05-10 | 中国测绘科学研究院 | 一种附加虚拟观测值的北斗三新频点定位方法及系统 |
CN114460615B (zh) * | 2022-04-14 | 2022-07-08 | 中国测绘科学研究院 | 一种附加虚拟观测值的北斗三新频点定位方法及系统 |
CN116859417A (zh) * | 2023-07-07 | 2023-10-10 | 哈尔滨工程大学 | 用于北斗ppp-rtk/mems的完好性监测方法 |
CN116859417B (zh) * | 2023-07-07 | 2024-04-30 | 哈尔滨工程大学 | 用于北斗ppp-rtk/mems的完好性监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114280633B (zh) | 2024-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110376621B (zh) | 一种基于北斗三号B2b信号的卫星定位方法及装置 | |
US10012738B2 (en) | Positioning method and positioning apparatus using satellite positioning system | |
EP3792665A1 (en) | Protection level generation methods and systems for applications using navigation satellite system (nss) observations | |
EP3627188A1 (en) | Correction information integrity monitoring in navigation satellite system positioning methods, systems, and devices | |
Paziewski et al. | Assessment of GPS+ Galileo and multi-frequency Galileo single-epoch precise positioning with network corrections | |
US7501981B2 (en) | Methods and apparatus to detect and correct integrity failures in satellite positioning system receivers | |
US10739471B2 (en) | GNSS receiver with a capability to resolve ambiguities using an uncombined formulation | |
CN110007317B (zh) | 一种选星优化的高级接收机自主完好性监测方法 | |
US8742984B2 (en) | Integrity method for differential corrections | |
CN114280633B (zh) | 一种非差非组合精密单点定位完好性监测方法 | |
CN111290005B (zh) | 载波相位的差分定位方法、装置、电子设备及存储介质 | |
CN104483690A (zh) | 一种gnss三频精密单点定位模糊度固定方法 | |
CN111983641B (zh) | 一种用于实时生成北斗星基增强系统完好性参数的方法 | |
EP1862809A1 (en) | GNSS signal integrity indicator | |
US20130090858A1 (en) | System for measuring coseismic movements or vibrations of structures based on global navigation satellite systems-gnss and/or pseudolites | |
CN115826016B (zh) | 一种北斗双频星基增强改正数及完好性参数解算的方法 | |
Gunning et al. | SBAS corrections for PPP integrity with solution separation | |
CN110208835A (zh) | 一种基于消电离层组合的跨系统紧组合差分定位方法 | |
CN104483691A (zh) | 一种gnss组合精密单点定位方法 | |
CN111929705A (zh) | 基于动态基准站的定位方法、装置、电子设备及存储介质 | |
CN115327593B (zh) | 基于无人机的定位方法、系统及存储介质 | |
Hutton et al. | Centimeter-level, robust GNSS-aided inertial post-processing for mobile mapping without local reference stations | |
Kazmierski | Performance of absolute real-time multi-GNSS kinematic positioning | |
Wang et al. | Comparison of three widely used multi‐GNSS real‐time single‐frequency precise point positioning models using the International GNSS Service real‐time service | |
CN114966787A (zh) | 定位方法、装置、设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |