CN103592657B - 一种基于钟差辅助的低可见星下单模raim实现方法 - Google Patents
一种基于钟差辅助的低可见星下单模raim实现方法 Download PDFInfo
- Publication number
- CN103592657B CN103592657B CN201310455355.5A CN201310455355A CN103592657B CN 103592657 B CN103592657 B CN 103592657B CN 201310455355 A CN201310455355 A CN 201310455355A CN 103592657 B CN103592657 B CN 103592657B
- Authority
- CN
- China
- Prior art keywords
- receiver
- raim
- satellite
- sigma
- clock
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000001514 detection method Methods 0.000 claims abstract description 33
- 239000011159 matrix material Substances 0.000 claims abstract description 18
- 239000013598 vector Substances 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 3
- 238000000556 factor analysis Methods 0.000 claims 2
- 239000012895 dilution Substances 0.000 abstract description 2
- 238000010790 dilution Methods 0.000 abstract description 2
- 230000008030 elimination Effects 0.000 description 6
- 238000003379 elimination reaction Methods 0.000 description 6
- 239000000243 solution Substances 0.000 description 6
- 101150116940 AGPS gene Proteins 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 5
- 125000000769 L-threonyl group Chemical group [H]N([H])[C@]([H])(C(=O)[*])[C@](O[H])(C([H])([H])[H])[H] 0.000 description 2
- 230000003416 augmentation Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000005433 ionosphere Substances 0.000 description 1
- 230000035484 reaction time Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/23—Testing, monitoring, correcting or calibrating of receiver elements
-
- 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/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/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于钟差辅助的低可见星下单模RAIM实现方法。它是在环路跟踪可见卫星较少情况下实现单模RAIM的一种方法。在接收机时钟频率稳定的前提下,通过牛顿插值模型进行钟差预测,在五颗可见卫星下能够进行故障卫星检测与排除,从而提高接收机定位精确性、连续性和稳定性。该方法首先通过卫星信号的载噪比等参数确定加权最小二乘中的权重因子,其次通过接收机的虚警率、漏警率、观测方程和权重矩阵进行RAIM可用性检测,利用非精密条件下能够容忍的虚警率进行检测门限值确定,通过牛顿插值模型进行接收机钟差预测,以实现故障卫星的检测与排除,最后对添加钟差辅助方程后的接收机精度因子进行分析比较。本发明在具有一颗冗余卫星条件下能够实现单颗故障卫星的检测与排除,通过钟差辅助模型降低了卫星的空间精度因子,提高了接收机定位准确性和***鲁棒性。
Description
技术领域
本发明属于卫星导航领域,具体的说是一种基于牛顿插值法钟差模型预测的低可见星下接收机自主完好性监测方法。
背景技术
全球导航卫星***(Global Navigation Satellite System,GNSS)在军事和民用领域都扮演着极其重要的地位,GNSS可为海、陆、空各个层面的物体进行定位、导航和授时,包括船舶远洋航行和进港引水,汽车自主导航以及飞机航路引导和进场降落,电力、邮电和通信网络***授时和校频等。接收机自主完好性监测(ReceiverAutonomous IntegrityMonitoring,RAIM)确保了用户定位、授时结果的精确性、连续性和鲁棒性,且为GNSS***的正常运行提供了保证。
在某些场景中,比如建筑物遮挡、城市峡谷或者室内条件下,接收机环路能够跟踪到的可见卫星数量很少,难以接收六颗卫星信号进行故障卫星检测与排除,故需要借助辅助信息提高***完好性监测性能。通常可采用卫星导航***差分定位服务中的广域增强***(Wide Area Augmentation System,WAAS)和局域增强***(Local Area AugmentationSystem,LAAS)进行完好性辅助监测,还可以采用辅助GPS(AGPS)和GPS-INS组合方式提高***性能,然而综合考虑外界环境、成本和接收机性能要求等因素,在本地接收机时钟源频率稳定的前提下,可不借助外界辅助信息进行完好性监测,即借助本地接收机钟差模型预测进行五颗可见卫星下的故障星检测与排除。本发明针对低可见卫星下(五颗可见卫星)单模GNSS定位,且提出了加权最小二乘法中权重因子的确定方法,最后给出定位结果的同时对精度因子(Dilution ofPrecision,DOP)进行了分析说明。
在卫星导航***领域的快照算法中,广泛采用最小二乘法、奇偶矢量法和加权最小二乘法,考虑到最小二乘法与奇偶矢量法定位结果的等价性,本发明采用加权最小二乘法(Weighted Least Squares,WLS)进行定位解算。
卫星导航***中校正后伪距测量值可表示为:
其中ρn代表伪距测量值,代表第n颗卫星钟差项,代表电离层延时,代表对流层延时。
RAIM算法中采用的加权最小二乘法可表示成:
b=G·X+ε (2)
其中b∈Rn×1,n为单模卫星导航***中参与解算的卫星数目,b代表校正后伪距与接收机至卫星几何距离之差,R表示实数域,G∈Rn×4,G代表接收机至卫星的方向余弦矩阵,x∈Rn×1,代表接收机位置坐标和钟差值,ε∈Rn×1,代表接收机测量误差矢量,正常情况下服从均值为0,方差为σ2的高斯分布。
本发明假定单模导航***中最多只存在1颗故障卫星,且假定第i颗卫星出错,加权最小二乘法中的第n颗卫星的测量误差标准差为σn,则权重矩阵W可表示成一个Rn×n的对角阵:
其中令C=WTW,则加权最小二乘解为:
x=(GTCG)-1GTCb (4)
定位后残余矢量b为:
其中矩阵S定义为:S=I-G(GTCG)-1GTC。
定义伪距残余平方加权和(Weighted Sum ofSquares for Error,WSSE)εWSSE:
由上式知,εWSSE表示加权后的残余矢量的长度平方。因为由加权最小二乘法得到的定位解可使加权残余各分量的平方和εWSSE最小,所以εWSSE值大小体现着各个测量值之间的一致性程度,且服从自由度(Degree ofFreedom,DOF)为N-4的χ2分布,其中N为参与解算的卫星个数。
漏警率Pmd和虚警率Pfa可作如下定义:
其中Tthr代表检测门限值,H0代表不存在故障卫星,H1代表存在故障卫星。ε=0时,εWSSE服从中心χ2分布;ε≠0时,εWSSE服从非中心χ2分布,则虚警率PFA和漏警率PMD分别表示为:
其中Г(α)代表伽玛函数,代表第一阶修正贝塞尔函数。在实际应用中,加权最小二乘法主要针对单星进行故障检测与排除,传统模式下,利用五颗卫星可进行故障检测,利用六颗卫星可进行故障排除,但在建筑物遮挡、城市峡谷或室内时,有时不能满足六颗可见卫星的情形,因此,有必要探索出新的方法使接收机只存在一颗冗余卫星的情况下(即五颗可见卫星)就能够进行故障卫星检测与排除,进而实现接收机自主完好性监测。
发明内容
本发明的目的是克服上述背景中的不足之处,提出一种适用于单模卫星导航***中只存在一颗冗余卫星就能够实现单星故障检测和排除的RAIM方法。
该方法通过卫星信号的载噪比确定权重矩阵因子,通过接收机的虚警率、漏警率进行RAIM的水平保护门限值(Horizontal Protection Level,HPL)确定,并通过虚警率和自由度进行RAIM检测门限值的确定,只存在一颗故障卫星的前提下,我们采用牛顿插值法进行接收机时钟预测,进而实现低可见星下单星故障检测与排除,同时给出了添加辅助方程后精度因子(Dilution ofPrecision,DOP)的变化趋势。
为了实现上述目的,本发明提出了一种基于钟差辅助的低可见星下单模RAIM实现方法,该方法包括以下六个步骤:
步骤一:权重因子确定。
根据接收机单模导航***类别,鉴相器因子、码环带宽、预检积分时间和载噪比等进行权重因子估算。本发明采用如下公式进行式(3)中权重因子确定:
其中σtDLL是接收机延时锁定环(DLL)1σ热噪声,σatm是1σ的大气延时,且热噪声可以看作载噪比C/N0的函数:
其中F1和F2代表鉴相器因子,d代表码片相关间隔,Bn代表码环噪声带宽,T代表预检积分时间,λc代表每一码片的波长,在GPS L1频点中λc=293.05[m],Compass B1I频点中λc=146.53[m],本发明假定电离层和对流程传播误差1σ值分别为5m和1.5m,且σatm为5.25m。
1σ的锁相环(PLL)测量误差是由多种噪声构成的,锁相环的热噪声估计公式如下:
GPS L1波段载波长度为λL1=0.1903[cycles/s],Compass B1I波段载波长度为λL1=0.1920[cycles/s],伪距观测误差通常包括接收机低频分量影响和大气延时,故我们的1σ锁相环(PLL)测量误差协方差为:
其中和用于计算每个接收机通道方差估计值,且这些参数用于加权最小二乘中自适应权重因子的估计方程中,用于计算接收机位置的加权因子,用于计算接收机速度的加权因子。
步骤二:RAIM可用性判定。
根据漏警率Pmd和虚警率Pfa确定最小检测误差Pblas,水平保护门限值的确定由测量误差到定位误差的几何因子SLOPEmax决定。
其中 ei∈Rn×1,且第i行值为1,其余全为0,C为权重矩阵。
HPL=SLOPEmax×Pblas (15)
水平定位误差PEH计算公式如下:
其中Ai,j代表A矩阵中的第i行和第j列,且A=(GTCG)-1GTC。
假定单模卫星导航***中只存在一颗故障卫星,且发生在第i颗,则第i颗卫星发生故障后有:
其中S=I-G(GTCG)-1GTC,Sn代表矩阵S中对角线上的第i个元素。SLOPEmax最终结果如下:
故通过SLOPEmax和Pblas可确定水平保护门限值,进而判定RAIM的可用性。
步骤三:确定RAIM检测门限。
加权残余向量的各分量加权平方和εWSSE呈自由度为N-4的χ2分布,本发明中N=5。根据虚警率Pfa确定检测门限值Threshold,若超过门限值,则执行步骤四,否则直接输出经过完好性监测得出的结果。
步骤四:确定钟差辅助方程。
采用具有承袭性的牛顿插值法进行接收机时钟模型预测,通过差商确定插值因子,在接收机时钟频率稳定的情形下,根据历史定位结果的接收机钟差值,可预测新的钟差值,等价于人为添加了钟差辅助方程。
f(t)在f[tk]处的零阶差商为f[tk]=f(tk),k=0,1,2...n,f(t)在ti,tj处的一阶差商为0≤i,j≤n,f(t)的k阶差商表示为:
由此构造出的牛顿插值多项式为:
f(t)=f(t0)+f[t0,t1](t-t0)+f[t0,t1,t2](t-t0)(t-t1)+...
f[t0,t1,...tn](t-t0)(t-t1)...(t-tn-1)+f[t,t0,t1,...tn](t-t0)(t-t1)...(t-tn)
本发明中滑动窗时间区间设置为ΔT=4s,每一解算历元间隔为1s,故在前ΔT的正确结果下,我们可以预测第(ΔT+1)次的接收机钟差值,接收机利用加权最小二乘法解算过程中等价于人为添加了如下钟差辅助方程:
其中ξ为接收机钟差预测误差,故将上式带入式(2)中,有
其中
步骤五:故障卫星检测与排除。
添加钟差辅助方程后,在五颗卫星中选择四颗并结合钟差辅助方程式(20),中会有一种组合使得加权平方和εWSSE小于门限值,从而实现故障卫星的检测与排除。
步骤六:WLS定位解算。
排除掉故障卫星或者在没有故障卫星的情况下,根据步骤一确定出的权重矩阵,利用公式(4)解算出接收机位置、钟差信息。
步骤七:调整滑动窗区间,并进行DOP性能分析。
成功解算出接收机位置及钟差后,调整钟差模型滑动窗区间。人为添加辅助方程后,DOP方程中的权系数阵修改为: 且为对称矩阵。故与未添加钟差辅助方程时H=(GTG)-1比较,GDOP、HDOP、PDOP、TDOP值都会降低,等价于间接的提高了接收机定位精度。
本发明的优点在于以下四个方面:
(1)本方法利用卫星信号的载噪比等信号参数作为权重矩阵中的加权因子,提高了***精确性。
(2)本方法利用牛顿插值法进行接收机钟差模型估计,预测准确、占用内存空间小、易于硬件实现。
(3)本方法较传统方法相比,最大的优势在于能够利用五颗卫星进行故障卫星检测与排除,特别适用于城市峡谷或室内接收机等信号不好的情形。
(4)本方法结构简单,故障排除率高,通过钟差辅助方程降低了DOP数值,进一步提高了接收机定位精度,且满足单模卫星导航***中RAIM要求。
附图说明
图1是本发明在低可见星下单模RAIM方法流程图;
图2是本发明接收机虚警率、漏警率关系图;
图3是本发明接收机定位结果与载噪比关系比较图;
图4是本发明在低可见星下单模RAIM接收机位置与误差关系图;
图5是添加辅助方程与原始方程DOP比较图;
具体实施方式
下面将结合附图和实例对本发明作进一步的详细说明。
如图1所示,按照发明内容步骤一至步骤七实现的一种基于钟差辅助方程的低可见星下单模RAIM方法流程图。接收机根据码环、载波环跟踪的卫星信号载噪比信息、信号传播过程的大气延时确定权重矩阵;根据非精密用户容许的虚警率、漏警率确定最小检测偏差,并通过误差几何因子SLOPEmax确定RAIM可用性;根据卡方分布的自由度及虚警率,确定RAIM检测门限值;根据牛顿插值模型添加接收机钟差辅助方程,并将定位后加权残余向量平方和与门限值比较,若存在故障卫星则进行排除,最后利用加权最小二乘法进行接收机定位解算,同时对卫星空间分布的几何因子进行分析说明。
如图2所示,当不存在故障卫星时,定位后加权残余向量平方和呈中心卡方分布,如图中的蓝色线所示,存在故障卫星时,呈非中心卡方分布,如图中的红色线所示,根据接收机能够容忍的虚警率,确定RAIM检测门限值,其中绿色区域围成的面积表示虚警率,洋红色区域围城的面积表示漏警率,虚警率降低的同时,漏警率会增加,故在定位解算时门限值的确定需要在二者之间进行平衡,本发明的虚警率设置为1×10-5,漏警率设置为0.2。
如图3所示,横坐标代表解算时间历元,纵坐标代表接收机在WGS-84坐标系下的位置偏差,其中数据来自IGS(international GNSS service)中的WAMC观测站在2013年3月14日10:01发布的星历数据,载噪比越高,信号质量越好,其中造成的定位误差越低。
如图4(a)所示,本发明通过选星算法选择五颗跟踪卫星进行定位解算,在没有故障卫星情况下,由加权最小二乘解算出的WAMC观测站三维误差基本在3m之内;如图4(b)所示,在解算历元80s之后,对用于解算的1号卫星人为以5m误差递增,最大误差递增到100m,之后再正常定位解算,观测站位置在存在故障卫星时间段内存在较大偏差;如图4(c)所示,在同样条件下添加基于钟差辅助方程的RAIM方法,初始条件下由于伪距偏差较少,不能检测到故障星的存在,故存在7m以内的偏差,当偏差逐渐增大时,接收机能够成功检测到故障卫星并进行隔离,故之后接收机的位置坐标恢复正常。
如图5所示,添加接收机钟差辅助方程后,卫星几何精度因子GDOP、空间位置精度因子PDOP和钟差精度因子都相应降低,其中CA代表ClockAuxiliary。
本发明提出了一种基于钟差辅助的低可见星下单模RAIM实现方法,该方法具体包括以下六个步骤:
步骤一:权重因子确定。
根据接收机单模导航***类别,鉴相器因子、码环带宽、预检积分时间和载噪比等进行权重因子估算。本发明以GPS单模导航***为例,σtDLL是接收机延时锁定环(DLL)1σ热噪声,σatm是1σ的大气延时,其中F1、F2取值都为1,码片相关间隔d=1/2[chips],码环噪声带宽Bn=2[Hz],预检积分时间T=2[ms],GPSL1频点每一码片的波长λc=293.05[m],电离层和对流程传播误差1σ值分别为5m和1.5m,且σatm为5.25m。
步骤二:RAIM可用性判定。
根据漏警率Pmd=0.2、虚警率Pfa=1×10-5、权重因子和观测矩阵,确定最小检测误差Pblas,水平保护门限值的确定由测量误差到定位误差的几何因子SLOPEmax决定。
通过SLOPEmax和Pblas可确定水平保护门限值,进而判定RAIM的可用性。
步骤三:确定RAIM检测门限。
加权残余向量的各分量加权平方和εWSSE呈自由度为N-4的χ2分布,本发明中N=5,虚警率Pfa=1×10-5,故确定检测门限值Threshold=19.5114。
步骤四:确定钟差辅助方程。
采用具有承袭性的牛顿插值法进行接收机时钟模型预测,通过差商确定插值因子,在接收机时钟频率稳定的情形下,根据历史定位结果的接收机钟差值,可预测新的钟差值,等价于人为添加了钟差辅助方程。
f(t)在f[t0]处的零阶差商为f[t0]=f(t0),f(t)在t0,t1处的一阶差商为 f(t)的三阶差商表示为:
故f(t)的牛顿均差型三次插值钟差辅助方程为:
f(t)=f(t0)+f[t0,t1](t-t0)+f[t0,t1,t2](t-t0)(t-t1)+f[t0,t1,t2,t3](t-t0)(t-t1)(t-t2)
接收机钟差滑动窗区间设置为ΔT=4s,每一解算历元间隔为1s,故在前ΔT的正确结果下,我们可以预测第(ΔT+1)次的接收机钟差值,接收机利用加权最小二乘法解算过程中等价于人为添加了如下钟差辅助方程:
故修正后的定位方程组为:
其中
步骤五:故障卫星检测与排除。
添加钟差辅助方程后,在五颗卫星中选择四颗并结合钟差辅助方程,中会有一种组合使得加权平方和εWSSE小于门限值,从而实现故障卫星的检测与排除。
步骤六:WLS定位解算。
排除掉故障卫星或者在没有故障卫星的情况下,根据步骤一确定出的权重矩阵,牛顿迭代法x=(GTCG)-1GTCb解算出接收机位置等信息。
步骤七:调整滑动窗区间,并进行DOP性能分析。
成功解算出接收机位置及钟差后,调整钟差模型滑动窗区间。人为添加辅助方程后,DOP方程中的权系数阵修改为: 故与未添加钟差辅助方程时对比H=(GTG)-1,GDOP、PDOP、TDOP值都会降低,等价于间接提高了接收机定位精度。未添加钟差辅助方程时,GDOP平均值为2.496,PDOP均值为2.196,TDOP均值为1.184,添加钟差辅助方程后,GDOP平均值为1.956,PDOP均值为1.801,TDOP均值为0.764。
Claims (7)
1.一种基于钟差辅助的低可见星下单模RAIM实现方法,其特征在于,包括以下七个步骤:
步骤A:根据卫星信号载噪比、鉴相器因子、码环噪声带宽、预检积分时间、码片波长、码片相关间隔参数、卫星信号大气延时确定权重矩阵,如公式1所示:
其中,σtDLL是接收机延时锁定环热噪声,σatm是大气延时,且热噪声可以看作载噪比C/N0的函数,如公式2所示:
其中,F1和F2代表鉴相器因子,d代表码片相关间隔,Bn代表码环噪声带宽,T代表预检积分时间,λc代表每一码片的波长;
步骤B:根据虚警率、漏警率进行RAIM可用性判定;
步骤C:确定RAIM检测门限;
步骤D:利用牛顿插值模型确定钟差辅助方程;采用具有承袭性的牛顿插值法进行接收机时钟模型预测,通过差商确定插值因子,在接收机时钟频率稳定的情形下,根据历史定位结果的接收机钟差值,进行新的钟差值预测,等价于人为添加了钟差辅助方程,采用牛顿均差型三次插值方法,钟差预测值如公式3所示:
f(t)=f(t0)+f[t0,t1](t-t0)+f[t0,t1,t2](t-t0)(t-t1)+f[t0,t1,t2,t3](t-t0)(t-t1)(t-t2) (3)
其中,f(t0)代表f(t)在f[t0]处的零阶差商,f[t0,t1]代表一阶差商,f[t0,t1,t2]代表二阶差商,f[t0,t1,t2,t3]代表三阶差商;
接收机利用加权最小二乘法解算过程中等价于人为添加了公式4钟差辅助方程:
其中,ξ为接收机钟差预测误差,σ代表加权最小二乘法中卫星的测量误差标准差;
步骤E:故障卫星检测与排除;
步骤F:利用加权最小二乘法进行接收机定位解算;
步骤G:调整钟差滑动窗区间,并进行精度因子分析。
2.根据权利要求1所述的低可见星下单模RAIM实现方法,其特征在于:步骤B通过虚警率、漏警率确定最小检测偏差,并通过误差几何因子SLOPEmax确定水平保护门限值,进而判断RAIM可用性。
3.根据权利要求1所述的低可见星下单模RAIM实现方法,其特征在于:步骤C根据卡方分布的自由度及虚警率要求,确定RAIM检测门限值。
4.根据权利要求1所述的低可见星下单模RAIM实现方法,其特征在于:步骤D构造均差型三次牛顿插值多项式,调整接收机至卫星的方向余弦矩阵和定位后残余矢量,并将所述步骤D中的钟差辅助方程添加到公式5中:
b=G·x+ε (5)
其中b∈Rn×1,n为单模卫星导航***中参与解算的卫星数目,b代表校正后伪距与接收机至卫星几何距离之差,R表示实数域,G∈Rn×4,G代表接收机至卫星的方向余弦矩阵,x∈Rn ×1,代表接收机位置坐标和钟差值,ε∈Rn×1,代表接收机测量误差矢量。
5.根据权利要求1所述的低可见星下单模RAIM实现方法,其特征在于:步骤E通过钟差辅助方程与原来颗卫星组合,对伪距残余平方加权和与步骤C中的门限值进行比较,实现故障卫星的检测与排除,所述伪距残余平方加权和公式如下所示:
其中,εWSSE代表伪距残余平方加权和,W代表权重矩阵,代表定位解算后的残余矢量,C=WTW。
6.根据权利要求1所述的低可见星下单模RAIM实现方法,其特征在于:步骤F根据步骤A中的权重矩阵进行接收机加权最小二乘法定位解算。
7.根据权利要求1所述的低可见星下单模RAIM实现方法,其特征在于:步骤G调整牛顿插值钟差滑动窗区间,并进行添加钟差辅助方程后的精度因子分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310455355.5A CN103592657B (zh) | 2013-09-30 | 2013-09-30 | 一种基于钟差辅助的低可见星下单模raim实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310455355.5A CN103592657B (zh) | 2013-09-30 | 2013-09-30 | 一种基于钟差辅助的低可见星下单模raim实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103592657A CN103592657A (zh) | 2014-02-19 |
CN103592657B true CN103592657B (zh) | 2016-08-31 |
Family
ID=50082867
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310455355.5A Expired - Fee Related CN103592657B (zh) | 2013-09-30 | 2013-09-30 | 一种基于钟差辅助的低可见星下单模raim实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103592657B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054216A (zh) * | 2016-05-24 | 2016-10-26 | 中国人民解放军信息工程大学 | 基于gdop和uere的多模gnss选星方法 |
CN109061683B (zh) * | 2018-06-29 | 2021-05-25 | 北京航空航天大学 | 一种优化水平保护级的h-araim*** |
CN109946722B (zh) * | 2019-04-01 | 2023-04-11 | 成都新橙北斗智联有限公司 | 一种多***多频段定位方法及*** |
CN110231636B (zh) * | 2019-04-29 | 2021-03-26 | 深圳市联和安业科技有限公司 | Gps与bds双模卫星导航***的自适应无迹卡尔曼滤波方法 |
CN110907953B (zh) * | 2019-10-18 | 2022-04-29 | 湖北三江航天险峰电子信息有限公司 | 一种卫星故障识别方法、装置及软件接收机 |
CN112505728B (zh) * | 2020-12-10 | 2024-03-22 | 中电科西北集团有限公司 | 时间完好性监测方法 |
CN113093248A (zh) * | 2021-04-11 | 2021-07-09 | 湘潭大学 | 一种范围一致性raim方法 |
CN113341915A (zh) * | 2021-05-27 | 2021-09-03 | 中国船舶重工集团公司第七0七研究所 | 一种四陀螺谐振惯导***故障诊断隔离方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1846147A (zh) * | 2003-09-05 | 2006-10-11 | 纳夫科姆技术公司 | 用于接收机自主的完整性监视和故障检测及消除的方法 |
CN101089651A (zh) * | 2007-07-12 | 2007-12-19 | 北京航空航天大学 | 接收机自主完好性监测的可用性预测方法 |
-
2013
- 2013-09-30 CN CN201310455355.5A patent/CN103592657B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1846147A (zh) * | 2003-09-05 | 2006-10-11 | 纳夫科姆技术公司 | 用于接收机自主的完整性监视和故障检测及消除的方法 |
CN101089651A (zh) * | 2007-07-12 | 2007-12-19 | 北京航空航天大学 | 接收机自主完好性监测的可用性预测方法 |
Non-Patent Citations (4)
Title |
---|
基于灰色理论的钟差预测模型辅助RAIM算法研究;滕云龙 等;《电子测量与仪器学报》;20100331;第24卷(第3期);257-262 * |
接收机时钟辅助RAIM算法研究;陈小平 等;《宇航学报》;20090131;第30卷(第1期);271-275 * |
时钟偏差辅助的GPS完整性监测算法;黄继勋 等;《北京航空航天大学学报》;20010430;第27卷(第2期);161-163 * |
时钟改进模型辅助RAIM算法研究;王永超 等;《电子学报》;20070630;第35卷(第6期);1084-1088 * |
Also Published As
Publication number | Publication date |
---|---|
CN103592657A (zh) | 2014-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103592657B (zh) | 一种基于钟差辅助的低可见星下单模raim实现方法 | |
US10018729B2 (en) | Selected aspects of advanced receiver autonomous integrity monitoring application to kalman filter based navigation filter | |
Groves et al. | Height aiding, C/N0 weighting and consistency checking for GNSS NLOS and multipath mitigation in urban areas | |
US7501981B2 (en) | Methods and apparatus to detect and correct integrity failures in satellite positioning system receivers | |
Sairo et al. | Weighted dilution of precision as quality measure in satellite positioning | |
KR101574819B1 (ko) | 비정상적 의사-거리 측정에 대해 무선 내비게이션 수신기 사용자를 보호하는 방법 | |
JP5232994B2 (ja) | Gps受信装置およびgps測位補正方法 | |
CN103592658A (zh) | 多模卫星导航***中基于选星算法的raim新方法 | |
US20160097859A1 (en) | Monitor based ambiguity verification for enhanced guidance quality | |
Mahmoud et al. | VANETs positioning in urban environments: A novel cooperative approach | |
KR20180088365A (ko) | 분리된 공식을 사용하여 모호성을 해결할 수 있는 능력을 가진 gnss수신기 | |
CN109085619B (zh) | 多模gnss***的定位方法及装置、存储介质、接收机 | |
Kuusniemi et al. | GNSS signal reliability testing in urban and indoor environments | |
GB2462926A (en) | GPS multipath reduction | |
CN109031356B (zh) | 特征斜率加权的最小二乘残差接收机自主完好性监测方法 | |
CN105487088A (zh) | 一种卫星导航***中基于卡尔曼滤波的raim算法 | |
Jia et al. | Ground vehicle navigation integrity monitoring for multi-constellation GNSS fused with cellular signals of opportunity | |
CN105158778A (zh) | 多***联合实施载波相位差分故障卫星剔除方法及其*** | |
Gioia | GNSS Navigation in difficult environments: Hybridization and Reliability | |
Gaglione et al. | Robust estimation methods applied to GPS in harsh environments | |
Shallberg et al. | WAAS measurement processing; current design and potential improvements | |
Leppakoski et al. | RAIM and complementary Kalman filtering for GNSS reliability enhancement | |
CN113721267B (zh) | 一种基于双接收机载波相位差分的gnss欺骗干扰检测方法 | |
KR101290085B1 (ko) | 다중 기준국 환경에서 대류층 지연 변칙현상 모니터링 방법 및 이를 이용한 시스템 | |
Mun et al. | Robust positioning algorithm for a yard transporter using GPS signals with a modified FDI and HDOP |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160831 Termination date: 20170930 |
|
CF01 | Termination of patent right due to non-payment of annual fee |