CN107219537A - 一种融合选星与完好性检测的多系统兼容定位方法 - Google Patents
一种融合选星与完好性检测的多系统兼容定位方法 Download PDFInfo
- Publication number
- CN107219537A CN107219537A CN201710379986.1A CN201710379986A CN107219537A CN 107219537 A CN107219537 A CN 107219537A CN 201710379986 A CN201710379986 A CN 201710379986A CN 107219537 A CN107219537 A CN 107219537A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- satellite
- msubsup
- satellites
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 239000011159 matrix material Substances 0.000 claims abstract description 52
- 239000013598 vector Substances 0.000 claims abstract description 47
- 238000001514 detection method Methods 0.000 claims abstract description 32
- 238000012544 monitoring process Methods 0.000 claims abstract description 10
- 239000011378 shotcrete Substances 0.000 claims description 15
- 230000009466 transformation Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 230000010365 information processing Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000012536 packaging technology Methods 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 230000035945 sensitivity Effects 0.000 claims description 4
- 230000009286 beneficial effect Effects 0.000 description 5
- 238000004806 packaging method and process Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000013024 troubleshooting Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/20—Integrity monitoring, fault detection or fault isolation of space segment
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/23—Testing, monitoring, correcting or calibrating of receiver elements
- G01S19/235—Calibration of receiver components
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/28—Satellite selection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Security & Cryptography (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开一种融合选星与完好性检测的多系统兼容定位方法,包括:对多卫星导航系统中的所有卫星进行卫星分组并构造各组的凸包集合;利用奇偶矢量算法得到含卫星测距误差信息的距离残差矢量及残差敏感矩阵;依据距离残差矢量判断有误故障卫星,依据残差敏感矩阵剔除故障卫星;使用礼品包装技术计算整个系列的凸包,实现同时完成多个卫星导航系统的选星和完好性监测;进行时间系统统一和坐标系统统一;进行多系统兼容定位。本发明的多系统兼容定位方法,利用分组的思想将系统分开,屏蔽了系统间差异,有助于提高剔除多系统多故障卫星的概率,加快了选星的速度,提高了故障的检测概率和可检测个数,有助于定位性能的提高。
Description
技术领域
本发明涉及融合选星与完好性检测的多系统兼容定位方法。更具体地,涉及一种融合选星与完好性检测的多卫星导航系统兼容定位方法。
背景技术
随着多卫星导航系统的逐步建立,对于接收机的信息处理也带来了挑战,除了保证有良好的卫星子集组合外,还需要确保参与定位的卫星健康。目前选星和完好性检测(RAIM)算法一般都是独立分开的,在多卫星导航系统时代,这种先选星再做完好性检测的处理方式,对于低成本接收机以及消费电子将是很大的挑战。
随着卫星导航星座的增加,同时多个故障卫星的可能性正在增加,尤其是不能忽视两个故障的概率。而传统的完好性检测方法都集中在单一的故障,其中,经典检测单一故障的RAIM算法包括伪距比较法、最小二乘残差法和奇偶矢量法。现有改进的RAIM算法大多是基于遍历的想法,进行复杂运算并多次验证,以确定多个卫星故障,运算量大、时间长且故障星越多越容易出错。这些算法的检出率相对较高,但所消耗的时间也很长。还有一些方法如GS和OWAS,这些方法可以检测同一个星座的多个故障,而两个星座的多故障就无效了。现有完好性兼容定位方法存在以下缺点:
a)故障卫星检测数量有限:大多数方法针对单一卫星系统或单一故障卫星进行检测;
b)运算量大且时间长:多个卫星故障的算法大多是基于遍历的想法,运算量大、时间长,满足不了卫星接收机实时定位的需求;
c)多系统检测概率低:可以检测同一个星座的多个故障,而两个星座的多故障就无效了。
因此,需要提供一种融合选星与完好性检测的多系统兼容定位方法,解决现有方法中故障卫星检测数量有限、运算量大且时间长和多系统检测概率低的问题。
发明内容
本发明的目的在于提供一种融合选星与完好性检测的多系统兼容定位方法。
为达到上述目的,本发明采用下述技术方案:
一种融合选星与完好性检测的多系统兼容定位方法,包括:
S1、对多卫星导航系统中的所有卫星进行卫星分组并构造各组的凸包集合;
S2:利用奇偶矢量算法得到含卫星测距误差信息的距离残差矢量及残差敏感矩阵;
S3:依据距离残差矢量判断有误故障卫星,依据残差敏感矩阵剔除故障卫星;
S4:使用礼品包装技术计算整个系列的凸包,实现同时完成多个卫星导航系统的选星和完好性监测;
S5:进行时间系统统一和坐标系统统一;
S6:进行多系统兼容定位。
优选地,步骤S1进一步包括如下子步骤:
S101、设A={b1,…,bk}为所有卫星的集合,若k/m为整数,根据卫星分组将集合A分为个组,每个组中卫星的数量为m,为对k/m的值进行向上取整操作,设j=1;
S102、计算凸包的划痕,找到布局平面图中第j组的左顶点对第j组中的其他点以关于左顶点的方位角大小顺序标记为删除靠近的点;
S103、设i=2、判断是否是逆时针方向,若是则将放入第j组的凸包集合Bj;
S104、i++后转入步骤S13,遍历完第j组中所有的点后转入步骤S15;
S105、j++后判断是否成立,若是则转入步骤S12,若否则步骤S1流程结束。
优选地,步骤S1进一步包括如下子步骤:
S106、设A={b1,…,bk}为所有卫星的集合,若k/m不为整数,根据卫星分组将集合A分为个组,每个组中卫星的数量最大为m,为对k/m的值进行向上取整操作,设j=1;
S107、计算凸包的划痕,找到布局平面图中第j组的左顶点对第j组中的其他点以关于左顶点的方位角大小顺序标记为删除靠近的点;
S108、设i=2、判断是否是逆时针方向,若是则将放入第j组的凸包集合Bj;
S109、i++后转入步骤S108,遍历完第j组中所有的点后转入步骤S110;
S110、j++后判断是否成立,若是则转入步骤S107,若否则转入步骤S111;
S111、计算第组凸包的划痕,第组卫星数量为找到布局平面图中第组的左顶点对第j组中的其他点以关于左顶点的方位角大小顺序标记为删除靠近的点,执行步骤S108、S109,步骤S1流程结束。
优选地,步骤S2进一步包括以下步骤:
S201:建立卫星导航系统观测方程式:
Δρ=GΔx+ε,
其中,Δρ为真实伪距与近似伪距偏离的差值向量,G为卫星与用户之间的几何观测矩阵,Δx为用户位置相对于真实位置的偏差向量,ε为伪距测量误差向量;
S202:确定随机误差项εunite各分量相互独立且服从相同分布,σ=σ1=σ2…=σn,协方差阵为σ2In,其中,σ、σ1、…σn为随机误差项各分量的标准差,In为n维单位矩阵;
S203:根据测量方程,得出Δxunite的最小二乘估值:
其中,为Δxunite的最小二乘估值,Gunite为多系统卫星与用户之间的几何观测矩阵,为Gunite的转置,Δρunite为多系统卫星真实伪距与近似伪距偏离的差值向量,Δxunite为多系统定位用户位置相对于真实位置的偏差向量;
S204:计算Δρunite的最小二乘估值:
S205:计算距离残差矢量w:
其中,S为残差敏感矩阵。
优选地,步骤S3进一步包括以下步骤:
S301:对几何观测矩阵Gunite进行QR分解:Gunite=QR;
S302:取Q的最后(n-5)列组成奇偶变换矩阵V,且有
p=VTw=VTΔρunite=VTεunite,
S=VVT,
pTp=wTw,
其中,p为奇偶矢量;
S303:无故障卫星时,p的协方差矩阵为:
其中,n为参与定位解算的卫星数,且n>7;
S304:判断pTp与检验门限T的大小,若pTp>T,则检测到有卫星故障,否则无卫星故障;
S305:若检测到有卫星故障,利用下式对卫星进行识别,以剔除故障卫星:
其中,vi为奇偶变换矩阵V的第i行,比较各Fi的大小,最大者判断为故障卫星,进行剔除。
进一步优选地,基于允许的虚警率和漏警率确定检验统计量得检验门限T。
优选地,步骤S4进一步包括以下步骤:
S401:在凸包集合B中选出最小x坐标b0;
S402:画通过b0的一条垂直线,对于任一Bi极值点,沿着这条线绕着b0逆时针旋转时最先遇到的点为下一点;
S403:给定bi和bi+1,下一点bi+2点应属于B且使∠bibi+1bi+2最大;
S404:重复上述步骤S401~S403,直至回到出发点b0。
优选地,步骤S5进一步包括以下步骤:
S501:通过接收机的基带信息处理系统跟踪到卫星后,通过测量卫星信号的到达时间可获得初步时间,进行接收机的时延自校准;
S502:当跟踪到多颗卫星后,通过定位处理计算出卫星时间,完成授时以北斗时为基准;
S503:在面向用户显示时,根据计算出来的时间差,进行时间转换输出各个卫星系统时间;
S504:在兼容定位算法中,将坐标系统统一到北斗的CGCS2000坐标上,当用户位置、卫星位置或者通过界面输入坐标需要转换时,通过布尔沙模型或者莫洛金斯基模型进行转换。
优选地,多系统兼容定位时,通过在最小二乘法的基础上增加权矩阵,对不同的卫星系统赋予不同的权重值来权衡各个系统对定位结果的影响,获得卫星数量和整体精度上的平衡,完成多系统兼容定位。
本发明的有益效果如下:
本发明中融合选星与完好性检测的多系统兼容定位方法,利用分组的思想将系统分开,屏蔽了系统间差异,有助于提高剔除多系统多故障卫星的概率。本发明实现了在选出最优星座结构的基础上,减少了卫星接收机待处理信号的通道数,降低了计算量,加快了选星的速度,提高了故障的检测概率和可检测个数,有助于定位性能的提高,解决了选星后若所选用卫星有故障会严重影响定位的问题。
本发明对新一代接收机硬件要求降低,成本减少,可推广应用于已老一代接收机的算法升级,避免了老一代接收机的硬件淘汰。对于导航的消费类电子应用以及对于通道少的接收机完成兼容多系统有重大参考意义。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明。
图1示出融合选星与完好性检测的多系统兼容定位方法步骤图。
图2示出时间系统统一过程示意图。
图3示出坐标系统统一过程示意图。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
选星与完好性检测都是为了提高卫星导航定位性能而提出的,目前由于两者单独使用,而造成计算量大和性能受影响。目前各个系统的可用性门限和故障监测门限情况不一,所以在考虑多系统多故障的时候,如果能将系统进行分组监测,那么可以根据不同系统的不同状况,确定相对应合理门限。
本发明将包装技术选星方法融入到奇偶矢量RAIM算法中,即在多系统定位前,利用包装技术选星进行卫星分组,分组后进行奇偶矢量RAIM算法检测,然后将筛选后的卫星分组进行组合选星,从而达到多系统快速正确的定位。
如图1所示,一种融合选星与完好性检测的多系统兼容定位方法,包括:
S1、对多卫星导航系统中的所有卫星进行卫星分组并构造各组的凸包集合。根据k/m是否为整数,可以判断是否所有分组中卫星数量均为m,具体分为以下两种情况:
情况一若k/m为整数:
S101、设A={b1,…,bk}为所有卫星的集合,根据卫星分组将集合A分为个组,每个组中卫星的数量为m,为对k/m的值进行向上取整操作,设j=1;
S102、计算凸包的划痕,找到布局平面图中第j组的左顶点对第j组中的其他点以关于左顶点的方位角大小顺序标记为删除靠近的点;
S103、设i=2、判断是否是逆时针方向,若是则将放入第j组的凸包集合Bj;
S104、i++后转入步骤S13,遍历完第j组中所有的点后转入步骤S15;
S105、j++后判断是否成立,若是则转入步骤S12,若否则步骤S1流程结束。
情况二若k/m不为整数:
S106、设A={b1,…,bk}为所有卫星的集合,根据卫星分组将集合A分为个组,每个组中卫星的数量最大为m,为对k/m的值进行向上取整操作,设j=1;
S107、计算凸包的划痕,找到布局平面图中第j组的左顶点对第j组中的其他点以关于左顶点的方位角大小顺序标记为删除靠近的点;
S108、设i=2、判断是否是逆时针方向,若是则将放入第j组的凸包集合Bj;
S109、i++后转入步骤S108,遍历完第j组中所有的点后转入步骤S110;
S110、j++后判断是否成立,若是则转入步骤S107,若否则转入步骤S111;
S111、计算第组凸包的划痕,第组卫星数量为找到布局平面图中第组的左顶点对第j组中的其他点以关于左顶点的方位角大小顺序标记为删除靠近的点,执行步骤S108、S109,步骤S1流程结束。
S2:利用奇偶矢量算法得到含卫星测距误差信息的距离残差矢量及残差敏感矩阵。
S201:建立卫星导航系统观测方程式:
Δρ=GΔx+ε,
其中,Δρ为真实伪距与近似伪距偏离的差值向量,G为卫星与用户之间的几何观测矩阵,Δx为用户位置相对于真实位置的偏差向量,ε为伪距测量误差向量;
S202:确定随机误差项εunite各分量相互独立且服从相同分布,
σ=σ1=σ2…=σn,协方差阵为σ2In,其中,σ、σ1、…σn为随机误差项各分量的标准差,In为n维单位矩阵;
S203:根据测量方程,得出Δxunite的最小二乘估值:
其中,为Δxunite的最小二乘估值,Gunite为多系统卫星与用户之间的几何观测矩阵,为Gunite的转置,Δρunite为多系统卫星真实伪距与近似伪距偏离的差值向量,Δxunite为多系统定位用户位置相对于真实位置的偏差向量;
S204:计算Δρunite的最小二乘估值:
S205:计算距离残差矢量w:
其中,S为残差敏感矩阵。
S3:依据距离残差矢量判断有误故障卫星,依据残差敏感矩阵剔除故障卫星。
S301:对几何观测矩阵Gunite进行QR分解:Gunite=QR;
S302:取Q的最后(n-5)列组成奇偶变换矩阵V,且有
p=VTw=VTΔρunite=VTεunite,
S=VVT,
pTp=wTw,
其中,p为奇偶矢量;
S303:无故障卫星时,p的协方差矩阵为:
其中,n为参与定位解算的卫星数,且n>7;
S304:判断pTp与检验门限T的大小,若pTp>T,则检测到有卫星故障,否则无卫星故障;
S305:若检测到有卫星故障,利用下式对卫星进行识别,以剔除故障卫星:
其中,vi为奇偶变换矩阵V的第i行,比较各Fi的大小,最大者判断为故障卫星,进行剔除。
S4:使用礼品包装技术计算整个系列的凸包,实现同时完成多个卫星导航系统的选星和完好性监测。
S401:在凸包集合B中选出最小x坐标b0;
S402:画通过b0的一条垂直线,对于任一Bi极值点,沿着这条线绕着b0逆时针旋转时最先遇到的点为下一点;
S403:给定bi和bi+1,下一点bi+2点应属于B且使∠bibi+1bi+2最大;
S404:重复上述步骤S401~S403,直至回到出发点b0。
S5:进行时间系统统一和坐标系统统一。
S501:通过接收机的基带信息处理系统跟踪到卫星后,通过测量卫星信号的到达时间可获得初步时间,进行接收机的时延自校准;
S502:当跟踪到多颗卫星后,通过定位处理计算出卫星时间,完成授时以北斗时为基准;
S503:在面向用户显示时,根据计算出来的时间差,进行时间转换输出各个卫星系统时间;
S504:在兼容定位算法中,将坐标系统统一到北斗的CGCS2000坐标上,当用户位置、卫星位置或者通过界面输入坐标需要转换时,通过布尔沙模型或者莫洛金斯基模型进行转换。
S6:进行多系统兼容定位。
多系统兼容定位时,通过在最小二乘法的基础上增加权矩阵,对不同的卫星系统赋予不同的权重值来权衡各个系统对定位结果的影响,获得卫星数量和整体精度上的平衡,完成多系统兼容定位。
本发明中选择一种计算量小且能选出较优星座的凸包选星算法,选星过程中其需要进行分组,此时改进的算法引入经典的奇偶矢量法进行RAIM监测排除故障星,之后进行最后融合并给出较优且健康的卫星组合。本发明中将多系统选星和完好性检测结合的方法,不仅降低了计算量,同时也尽可能的选出最优的星座结构,有助于定位性能的提高。该方法利用分组的思想将系统分开屏蔽了系统间差异,有助于提高剔除多系统多故障卫星的概率。
本发明的定位方法原理如下:
首先,B={b1,...,bn}是所有卫星集合,根据卫星分组卫星集合B被区分为组,每个大小至多为m,计算凸包的划痕,找到最左边的最顶点b0,对所有其他的点以关于b0的角度分别标记为{b1,...,bm},为了防止约束,删除靠近b0的点。然后,d1←b1,d2←b2,i←2这时判断di,di-1,bi是逆时针方向,则bi属于集合D,i++。最后,依次类推得到凸包集合D。其次,利用奇偶矢量算法得到距离残差矢量含卫星测距误差信息,可用作判断有无故障星的依据和残差敏感矩阵剔除故障卫星。最后,使用的礼品包装技术计算整个系列的凸包,即同时完成了多个卫星导航系统的选星和完好性监测。对于通道较少的接收机使用此算法,可以采用分时复用的方式工作,大大节约了硬件成本。
利用奇偶矢量算法得到距离残差矢量含卫星测距误差信息,可用作判断有无故障星的依据和残差敏感矩阵剔除故障卫星,具体过程如下:
根据最小二乘原理,卫星导航系统观测方程式:
Δρ=GΔx+ε,
其中,Δρ为真实伪距与近似伪距偏离的差值向量;G为卫星与用户之间的几何观测矩阵;Δx为用户位置相对于真实位置的偏差向量;ε为伪距测量误差向量。
奇偶矢量法基于以下假设:假设随机误差项εunite各分量相互独立且服从相同分布,即σ=σ1=σ2…=σn,协方差阵为σ2In。根据测量方程,得出Δxunite的最小二乘估值:
于是有
得到距离残差矢量w
其中,S为残差敏感矩阵。
对几何观测矩阵Gunite进行QR分解,即Gunite=QR,取Q的最后(n-5)列组成奇偶变换矩阵V。显然,V中各列相互正交,大小归一,并且与Gunite的各列相互正交。因此,存在下式:
p=VTw=VTΔρunite=VTεunite,
可以证明:
S=VVT,
因此,下式成立:
pTp=wTw,
无故障卫星时,p的协方差矩阵为:
当参与定位解算的卫星数大于7颗时,各距离残差的平方和pTp服从自由度为(n-5)的χ2分布。当有一颗卫星出现故障时,各距离残差的平方和wTw服从自由度为(n-5)的非中心χ2分布。当系统为三个系统时,自由度为(n-6),单系统情形自由度为(n-4)。
RAIM算法用各距离残差的平法和作为检验统计量,基于允许的虚警率和漏警率可以确定检验统计量得检验门限T。若pTp>T,则检测到有卫星故障,否则无卫星故障。
若检测到有卫星故障,利用下式对卫星进行识别,以剔除故障卫星。
其中,vi为奇偶变换矩阵V的第i行,比较各Fi的大小,最大者判断为故障卫星,应当剔除。
利用礼品包装技术计算整个系列的凸包时,全新集合里选出最小x坐标b0,也就是集合里的一个顶点。画通过b0的一条垂直线以找到下一个点。对于任一Bi极值点,沿着这条线绕着b0逆时针旋转时最先遇到的点就是要找的点。给定b0和b1,下一点b2点应属于B且使∠b0b1li最大。继续这些包装步骤,直到又回到出发点b0。每个包装的步骤可以通过在每个子集Bi上执行极值查询。每次查询扫描所有的m个点。如此循环反复,凸包的构造将短时间内完成。
如图2和图3所示,时间系统统一和坐标系统的统一过程如下:
多系统兼容接收机具有已知的三维位置(X,Y,Z),通过接收到任何一颗可用的北斗二代卫星、GPS卫星、GALILEO卫星或者GLONASS卫星的信号,进而测量卫星信号的到达时间则可与该颗卫星建立一个一维方程,解算这个方程,即可得出用户接收机的时间T。解算出的T是锁定于卫星时基的,用于控制用户接收机的定时输出。若接收到一颗以上的卫星,则可以具备有一定的容错性并改善授时的精度。多个卫星导航系统的时间差可以通过多方程组求解或者利用采用电文里播发的时间差信息来补偿修正。
北斗卫星导航系统参考坐标系为CGCS2000,GPS卫星导航系统参考坐标系为WGS84,GLONASS卫星导航系统参考坐标系为PZ-90,GALILEO卫星导航系统参考坐标系为CTRF。GALILEO地面参考框架CTRF是国际地球参考时ITRS的一个独立实现,与WGS84一样,CTRF通过一系列伽利略地面监控站的坐标来实现。两个系统之间的差异是厘米级的,可以忽略。北斗与GPS的轨道参数只有非常细小的区别,两者坐标系的差异也很小,唯一与其他三个卫星导航系统参考坐标系差别最大的是GLONASS。以PZ-90转换成WGS84坐标系统为例,它们间存在一定的平移与旋转。确定坐标系转换参数的数学模型中,一般用布尔沙模型(Bursa Model)或者莫洛金斯基模型(Molodensky Model),可以认为两种模型的转换结果是等价的,以下参数采用的是Bursa模型。则坐标转换的七参数布尔沙模型为:
其中:为坐标系A到坐标系B的三个旋转参数;为坐标系A到坐标系B的三个平移参数;为坐标系A到坐标系B之间的尺度参数。
兼容定位包含了信息处理策略和兼容加权定位。多个卫星导航系统的兼容接收机最大可见星数从单系统的8、9颗提高到接近40颗,频点数由原来的1~3个频点增加到11个频点。应用模式的多样需要信息处理策略的多样,这给信息处理的各个方面都提出了挑战。多系统兼容接收机用户端算法首先从单一系统来看(以北斗导航系统为例),当通过信号前端接收处理后给出了B1、B2、B3各个频点的观测数据,数据经过电文解析和伪距生成,即可解算卫星位置以及处理伪距(载波相位平滑伪距和多频组合伪距),再进行定位和完好性处理。然后,在进行多系统兼容算法前实现时空统一。最后在兼容定位中,由于多个卫星导航系统观测值分别来自于不同系统,其精度取决于码速率以及环路算法实现和环境的影响,测距精度各不相同。往往对确定权矩阵所基于的观测值噪声特性了解不多,很难适当地确定权矩阵,通常赋予一个简单的权阵比如单位阵,这样可能会对复杂的数据噪声过于简单化。因此,将使得最小二乘解失去它的最小方差性质,导致不合理的单点定位结果。为提高定位精度,采用加权最小二乘法进行解算,通过在最小二乘法的基础上增加权矩阵,对不同的卫星系统赋予不同的权重值来权衡各个系统对定位结果的影响,获得卫星数量和整体精度上的平衡,完成多系统兼容定位。
本发明对于多系统采用分组处理的方法不仅能够加快选星的速度,而且也能提高故障的检测概率和可检测个数。该算法不仅将多星座选星和RAIM监测融合,通过新的选星算法加快RAIM检测算法且提高检测概率,而且解决了选星后若所选用卫星有故障会严重影响定位的缺点。本发明不仅减少了卫星接收机待处理信号的通道数,且接收机运算量降低,对新一代接收机硬件要求降低,成本减少。同时,也可推广应用于已老一代接收机的算法升级,避免了老一代接收机的硬件淘汰。对于导航的消费类电子应用以及对于通道少的接收机完成兼容多系统有重大参考意义。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (9)
1.一种融合选星与完好性检测的多系统兼容定位方法,其特征在于,包括:
S1、对多卫星导航系统中的所有卫星进行卫星分组并构造各组的凸包集合;
S2:利用奇偶矢量算法得到含卫星测距误差信息的距离残差矢量及残差敏感矩阵;
S3:依据所述距离残差矢量判断有误故障卫星,依据残差敏感矩阵剔除故障卫星;
S4:使用礼品包装技术计算整个系列的凸包,实现同时完成多个卫星导航系统的选星和完好性监测;
S5:进行时间系统统一和坐标系统统一;
S6:进行多系统兼容定位。
2.根据权利要求1所述的融合选星与完好性检测的多系统兼容定位方法,其特征在于,所述步骤S1进一步包括如下子步骤:
S101、设A={b1,…,bk}为所有卫星的集合,若k/m为整数,根据卫星分组将集合A分为个组,每个组中卫星的数量为m,为对k/m的值进行向上取整操作,设j=1;
S102、计算凸包的划痕,找到布局平面图中第j组的左顶点对第j组中的其他点以关于左顶点的方位角大小顺序标记为删除靠近的点;
S103、设i=2、判断是否是逆时针方向,若是则将放入第j组的凸包集合Bj;
S104、i++后转入步骤S13,遍历完第j组中所有的点后转入步骤S15;
S105、j++后判断是否成立,若是则转入步骤S12,若否则步骤S1流程结束。
3.根据权利要求1所述的融合选星与完好性检测的多系统兼容定位方法,其特征在于,所述步骤S1进一步包括如下子步骤:
S106、设A={b1,…,bk}为所有卫星的集合,若k/m不为整数,根据卫星分组将集合A分为个组,每个组中卫星的数量最大为m,为对k/m的值进行向上取整操作,设j=1;
S107、计算凸包的划痕,找到布局平面图中第j组的左顶点对第j组中的其他点以关于左顶点的方位角大小顺序标记为删除靠近的点;
S108、设i=2、判断是否是逆时针方向,若是则将放入第j组的凸包集合Bj;
S109、i++后转入步骤S108,遍历完第j组中所有的点后转入步骤S110;
S110、j++后判断是否成立,若是则转入步骤S107,若否则转入步骤S111;
S111、计算第组凸包的划痕,第组卫星数量为找到布局平面图中第组的左顶点对第j组中的其他点以关于左顶点的方位角大小顺序标记为删除靠近的点,执行步骤S108、S109,步骤S1流程结束。
4.根据权利要求1所述的融合选星与完好性检测的多系统兼容定位方法,其特征在于,所述步骤S2进一步包括以下步骤:
S201:建立卫星导航系统观测方程式:
Δρ=GΔx+ε,
其中,Δρ为真实伪距与近似伪距偏离的差值向量,G为卫星与用户之间的几何观测矩阵,Δx为用户位置相对于真实位置的偏差向量,ε为伪距测量误差向量;
S202:确定随机误差项εunite各分量相互独立且服从相同分布,σ=σ1=σ2…=σn,协方差阵为σ2In,其中,σ、σ1、…σn为随机误差项各分量的标准差,In为n维单位矩阵;
S203:根据测量方程,得出Δxunite的最小二乘估值:
<mrow>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
<mi>T</mi>
</msubsup>
<msub>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
<mi>T</mi>
</msubsup>
<msub>
<mi>&Delta;&rho;</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
其中,为Δxunite的最小二乘估值,Gunite为多系统卫星与用户之间的几何观测矩阵,为Gunite的转置,Δρunite为多系统卫星真实伪距与近似伪距偏离的差值向量,Δxunite为多系统定位用户位置相对于真实位置的偏差向量;
S204:计算Δρunite的最小二乘估值:
<mrow>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>&rho;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
<mi>T</mi>
</msubsup>
<msub>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
<mi>T</mi>
</msubsup>
<msub>
<mi>&Delta;&rho;</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>;</mo>
</mrow>
S205:计算距离残差矢量w:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>w</mi>
<mo>=</mo>
<msub>
<mi>&Delta;&rho;</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>&rho;</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
<mi>T</mi>
</msubsup>
<msub>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
<mi>T</mi>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>&Delta;&rho;</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>S&epsiv;</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>S</mi>
<mo>=</mo>
<msub>
<mi>I</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
<mi>T</mi>
</msubsup>
<msub>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msubsup>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
<mi>T</mi>
</msubsup>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,S为残差敏感矩阵。
5.根据权利要求4所述的融合选星与完好性检测的多系统兼容定位方法,其特征在于,所述步骤S3进一步包括以下步骤:
S301:对几何观测矩阵Gunite进行QR分解:Gunite=QR;
S302:取Q的最后(n-5)列组成奇偶变换矩阵V,且有
p=VTw=VTΔρunite=VTεunite,
S=VVT,
pTp=wTw,
其中,p为奇偶矢量;
S303:无故障卫星时,p的协方差矩阵为:
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>pp</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>E</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>V</mi>
<mi>T</mi>
</msup>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
</msub>
<msubsup>
<mi>&epsiv;</mi>
<mrow>
<mi>u</mi>
<mi>n</mi>
<mi>i</mi>
<mi>t</mi>
<mi>e</mi>
</mrow>
<mi>T</mi>
</msubsup>
<mi>V</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</msub>
<mo>,</mo>
</mrow>
其中,n为参与定位解算的卫星数,且n>7;
S304:判断pTp与检验门限T的大小,若pTp>T,则检测到有卫星故障,否则无卫星故障;
S305:若检测到有卫星故障,利用下式对卫星进行识别,以剔除故障卫星:
<mrow>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msup>
<mi>p</mi>
<mi>T</mi>
</msup>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<mo>|</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
</mrow>
</msqrt>
</mfrac>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>n</mi>
</mrow>
其中,vi为奇偶变换矩阵V的第i行,比较各Fi的大小,最大者判断为故障卫星,进行剔除。
6.根据权利要求5所述的融合选星与完好性检测的多系统兼容定位方法,其特征在于,基于允许的虚警率和漏警率确定检验统计量得检验门限T。
7.根据权利要求2或3所述的融合选星与完好性检测的多系统兼容定位方法,其特征在于,所述步骤S4进一步包括以下步骤:
S401:在凸包集合B中选出最小x坐标b0;
S402:画通过b0的一条垂直线,对于任一Bi极值点,沿着这条线绕着b0逆时针旋转时最先遇到的点为下一点;
S403:给定bi和bi+1,下一点bi+2点应属于B且使∠bibi+1bi+2最大;
S404:重复上述步骤S401~S403,直至回到出发点b0。
8.根据权利要求1所述的融合选星与完好性检测的多系统兼容定位方法,其特征在于,所述步骤S5进一步包括以下步骤:
S501:通过接收机的基带信息处理系统跟踪到卫星后,通过测量卫星信号的到达时间可获得初步时间,进行接收机的时延自校准;
S502:当跟踪到多颗卫星后,通过定位处理计算出卫星时间,完成授时以北斗时为基准;
S503:在面向用户显示时,根据计算出来的时间差,进行时间转换输出各个卫星系统时间;
S504:在兼容定位算法中,将坐标系统统一到北斗的CGCS2000坐标上,当用户位置、卫星位置或者通过界面输入坐标需要转换时,通过布尔沙模型或者莫洛金斯基模型进行转换。
9.根据权利要求1所述的融合选星与完好性检测的多系统兼容定位方法,其特征在于,多系统兼容定位时,通过在最小二乘法的基础上增加权矩阵,对不同的卫星系统赋予不同的权重值来权衡各个系统对定位结果的影响,获得卫星数量和整体精度上的平衡,完成多系统兼容定位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710379986.1A CN107219537A (zh) | 2017-05-25 | 2017-05-25 | 一种融合选星与完好性检测的多系统兼容定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710379986.1A CN107219537A (zh) | 2017-05-25 | 2017-05-25 | 一种融合选星与完好性检测的多系统兼容定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107219537A true CN107219537A (zh) | 2017-09-29 |
Family
ID=59944750
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710379986.1A Pending CN107219537A (zh) | 2017-05-25 | 2017-05-25 | 一种融合选星与完好性检测的多系统兼容定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107219537A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109728868A (zh) * | 2018-11-27 | 2019-05-07 | 中国科学院光电研究院 | 一种基于多重完好性检验的gnss基站组网时间同步方法 |
CN109861739A (zh) * | 2019-01-22 | 2019-06-07 | 北京电子工程总体研究所 | 一种通信卫星指定时空范围内的效能值排序方法和系统 |
WO2019203678A1 (en) * | 2018-04-18 | 2019-10-24 | Limited Liability Company "Topcon Positioning Systems" | Using sdp relaxation for optimization of the satellites set chosen for positioning |
CN111965668A (zh) * | 2020-07-14 | 2020-11-20 | 南京航空航天大学 | 一种面向卫星多故障的raim方法 |
CN113654573A (zh) * | 2021-10-15 | 2021-11-16 | 成都云智北斗科技有限公司 | 地面机动载体组合导航系统粗对准方法及其设备 |
CN113866800A (zh) * | 2021-12-03 | 2021-12-31 | 深圳市麦微智能电子有限公司 | 一种gnss高精度导航信息汇总系统 |
CN113933876A (zh) * | 2021-11-16 | 2022-01-14 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 多星通讯时差定位数据融合处理方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105549033A (zh) * | 2015-12-07 | 2016-05-04 | 郑州威科姆科技股份有限公司 | 一种基于最小二乘残差边缘检验的完好性处理方法 |
-
2017
- 2017-05-25 CN CN201710379986.1A patent/CN107219537A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105549033A (zh) * | 2015-12-07 | 2016-05-04 | 郑州威科姆科技股份有限公司 | 一种基于最小二乘残差边缘检验的完好性处理方法 |
Non-Patent Citations (3)
Title |
---|
LEI ZHANG ET AL.: "Enhancements of the satellite selection method for multi-constellation GNSS using convex geometry", 《IET INTERNATIONAL RADAR CONFERENCE 2013》 * |
NURIA BLANCO-DELGADO ET AL.: "A Convex Geometry Approach to Dynamic GNSS Satellite Selection for a Multi-Constellation System", 《PROCEEDINGS OF THE 22ND INTERNATIONAL TECHNICAL MEETING OF THE SATELLITE DIVISION OF THE INSTITUTE OF NAVIGATION (ION GNSS 2009)》 * |
NURIA BLANCO-DELGADO ET AL.: "Satellite Selection Method for Multi-Constellation GNSS Using Convex Geometry", 《IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019203678A1 (en) * | 2018-04-18 | 2019-10-24 | Limited Liability Company "Topcon Positioning Systems" | Using sdp relaxation for optimization of the satellites set chosen for positioning |
CN109728868A (zh) * | 2018-11-27 | 2019-05-07 | 中国科学院光电研究院 | 一种基于多重完好性检验的gnss基站组网时间同步方法 |
CN109728868B (zh) * | 2018-11-27 | 2020-10-13 | 中国科学院光电研究院 | 一种基于多重完好性检验的gnss基站组网时间同步方法 |
CN109861739A (zh) * | 2019-01-22 | 2019-06-07 | 北京电子工程总体研究所 | 一种通信卫星指定时空范围内的效能值排序方法和系统 |
CN109861739B (zh) * | 2019-01-22 | 2021-09-03 | 北京电子工程总体研究所 | 一种通信卫星指定时空范围内的效能值排序方法和系统 |
CN111965668A (zh) * | 2020-07-14 | 2020-11-20 | 南京航空航天大学 | 一种面向卫星多故障的raim方法 |
CN111965668B (zh) * | 2020-07-14 | 2023-10-24 | 南京航空航天大学 | 一种面向卫星多故障的raim方法 |
CN113654573A (zh) * | 2021-10-15 | 2021-11-16 | 成都云智北斗科技有限公司 | 地面机动载体组合导航系统粗对准方法及其设备 |
CN113654573B (zh) * | 2021-10-15 | 2022-03-08 | 成都云智北斗科技有限公司 | 地面机动载体组合导航系统粗对准方法及其设备 |
CN113933876A (zh) * | 2021-11-16 | 2022-01-14 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 多星通讯时差定位数据融合处理方法 |
CN113933876B (zh) * | 2021-11-16 | 2023-05-23 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 多星通讯时差定位数据融合处理方法 |
CN113866800A (zh) * | 2021-12-03 | 2021-12-31 | 深圳市麦微智能电子有限公司 | 一种gnss高精度导航信息汇总系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107219537A (zh) | 一种融合选星与完好性检测的多系统兼容定位方法 | |
CN106168672B (zh) | 一种gnss多模单频rtk周跳探测方法及装置 | |
CN104714244B (zh) | 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法 | |
CN107132558A (zh) | 惯性辅助的多频多模gnss周跳修复方法及系统 | |
CN110007317B (zh) | 一种选星优化的高级接收机自主完好性监测方法 | |
CN105487088B (zh) | 一种卫星导航系统中基于卡尔曼滤波的raim算法 | |
JP3875276B2 (ja) | Gps受信機の障害検出の方法およびシステム | |
CN108873029B (zh) | 一种实现导航接收机钟差建模的方法 | |
CN106054216A (zh) | 基于gdop和uere的多模gnss选星方法 | |
CN106646537B (zh) | 一种抗多径的gnss快速选星方法及装置 | |
CN103630914B (zh) | 一种gnss基线解算参考卫星选择方法 | |
CN104267410B (zh) | 机载完好性监测中多故障的排除方法和装置 | |
CN109061683B (zh) | 一种优化水平保护级的h-araim系统 | |
CN107505636A (zh) | 海基jpals的定位域mrcc方法及装置 | |
CN105549033A (zh) | 一种基于最小二乘残差边缘检验的完好性处理方法 | |
CN104459722B (zh) | 一种基于多余观测分量的整周模糊度可靠性检验方法 | |
GB2487256A (en) | Location fix from an unknown position | |
CN109752739B (zh) | 一种观测数据处理方法、装置、终端设备及存储介质 | |
CN108089210A (zh) | 一种惯性信息辅助raim检测方法 | |
CN111538056B (zh) | 动态精密单点定位解算方法 | |
CN106093991A (zh) | 一种用于gnss定位的模糊度快速恢复方法及系统 | |
CN104407368A (zh) | 基于北斗卫星导航系统的高精度差分信息生成方法 | |
CN107422342A (zh) | Gnss卫星钟差实时估计质量控制方法 | |
CN107966722A (zh) | 一种gnss钟差解算方法 | |
CN106932793A (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170929 |