CN102901971B - Parity vector method-based double-satellite failure recognition method - Google Patents
Parity vector method-based double-satellite failure recognition method Download PDFInfo
- Publication number
- CN102901971B CN102901971B CN201210365524.1A CN201210365524A CN102901971B CN 102901971 B CN102901971 B CN 102901971B CN 201210365524 A CN201210365524 A CN 201210365524A CN 102901971 B CN102901971 B CN 102901971B
- Authority
- CN
- China
- Prior art keywords
- fault
- satellite
- satellites
- parity vector
- detect
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention relates to global satellite navigation system receiver autonomous integrity monitoring technology and discloses a parity vector method-based double-satellite failure recognition method, aiming at the problems of false positives and false negatives caused by fault deviation offsetting when the parity vector method is used for recognizing two fault satellites. According to the technical scheme, the parity vector method is used for recognizing one fault satellite; with the fault satellite as the basis, four fault-free satellites are found out, and the information of the fault-free satellites is used for roughly locating, so that the fault satellites can be recognized; the recognized fault satellites are removed, and then the position resolution is carried out again, so that the locating accuracy is improved; therefore, the problems of false positives or false negatives caused by fault deviation offsetting can be avoided. The method solves the problem of the fault deviation offsetting caused by parity vector residual error and realizes the detection and the recognition for a plurality of fault satellites. After the method is used for detecting and recognizing satellite failure, the locating accuracy is improved. The method is mainly used for monitoring the autonomous integrity of a global satellite navigation system receiver.
Description
Technical field
The present invention relates to satellite navigation and location technology, be particularly related to GPS (Global Position System) (Global Navigation Satellite System, GNSS) receiver autonomous integrity monitoring (Receiver Autonomous Integrity Monitoring, RAIM) method.
Background technology
GNSS receiver, in obtaining navigational solution, must be taken the impact of pseudorange deviation on navigational solution into account.The principal element that causes pseudorange deviation has that larger satellite clock floats, the fault of the incorrect and satellite ingredient of navigation message data, is referred to as satellite failure here.Because system cannot guarantee the reaction time to satellite failure, the fast monitored of satellite failure is only carried out at user side so, thereby has occurred receiver autonomous integrity monitoring method.
RAIM utilizes redundancy distance measuring signal to detect to cause the satellite mistake of larger deviations.Conventionally RAIM comprises two functions: fault detect (Fault Detection) and Fault Identification (Fault Identification).Comparatively conventional RAIM method comprises: Agonists by Distance Comparison Method, least square method and parity vector method.But these methods are all under the supposed premise based on a fault satellites.Along with overhead satellites number constantly increases, the probability that two or multi-satellite break down simultaneously can not be left in the basket again.
Traditional fault satellites recognition methods can not be directly used in many fault satellites identifications, and not perfect about the recognition methods of many fault satellites.For example, parity vector method there will be fault deviation to offset while identification for two fault satellites, and brings thus the situation of failing to judge and judging by accident.
The principle of a fault satellites of parity vector method identification is as follows:
The basic observation equation of receiver positioning calculation is:
y=Hx+e (1)
In formula (1), y is the residual vector of Pseudo-range Observations, and H is observing matrix, and x is receiver location residual error and receiver clock error to be calculated, the measuring error that e is Pseudo-range Observations.
In the time that non-fault satellite exists, think that measuring error e obeys variance and is
normal distribution; In the time having fault satellites to exist, there is deviation in the pseudorange of fault satellites, and now, e is Normal Distribution no longer.
When application parity vector method is carried out fault satellites detection with identification, first observing matrix H is decomposed into the long-pending of an orthogonal matrix Q and a upper triangular matrix R, even:
H=QR (2)
By the transposed matrix Q of matrix Q
tpiecemeal:
In formula (3), Q
p∈ R
(n-4) × nfor odd even space matrix, wherein n is number of satellite, Q
xfor Q
tfront four lines:
p=Q
py=Q
pe (4)
In formula (4), p is parity vector.Pseudorange residual sum of squares (RSS) SSE can construct with parity vector p:
SSE=p
Tp (5)
While not there is not fault satellites,
obeying degree of freedom is the centralization card side distribution of (n-4),
while there is fault satellites,
the card side that obeys decentralization distributes,
λ is decentralization parameter.Therefore, can construct fault detect judgement amount
As overall false-alarm probability α
0after given, can calculate fault detect thresholding σ
t, nfor:
In formula (7), icdf is (1-α) fractile that computer card side distributes.
time think have fault satellites to exist, and then identification fault satellites, otherwise, think non-fault satellite, the receiver information of current calculating is accurately.The fault detect judgement amount here
with fault detect thresholding σ
t, n, for carrying out fault detect.
From formula (4), known parity vector p is through Q by observational error e
pprojection gets, and can carry out fault satellites identification with the geometric properties between them.Structure Fault Identification judgement amount:
When all satellites are during without pseudorange deviation,
therefore work as the false-alarm probability α of single satellite
1give timing, can calculate Fault Identification thresholding:
In formula (9), erf (x) is Gauss error function, erf
-1(x) be its inverse function.Work as r
ibe greater than T
rtime, can judge that i satellite is fault satellites, there is deviation in its pseudorange.After i satellite rejected, utilize all the other satellites again to position and resolve to improve positioning precision.The Fault Identification judgement amount r of formula (8) and formula (9)
iwith Fault Identification thresholding T
r, n, for carrying out Fault Identification.
Summary of the invention
Technical matters to be solved by this invention, there will be fault deviation to offset while identification exactly, and bring thus the shortcoming of failing to judge and judging by accident for parity vector method for two fault satellites, and a kind of double star fault recognition method based on parity vector method is provided.
The present invention solve the technical problem, and the technical scheme of employing is that the double star fault recognition method based on parity vector method, comprises step:
(1), from an existing n satellite, reject satellite i, n is receiver visible satellite quantity in global position system, n >=7;
(2), to remaining (n-1) satellite according to fault detect judgement amount
with fault detect thresholding σ
t, (n-1)carry out fault detect, if
think that remaining non-fault satellite in satellite exists, and detects and finishes; Otherwise, carry out fault satellites identification by parity vector method, according to Fault Identification judgement amount r
imaximal value r
j=max (r
1, r
2..., r
n), satellite j is rejected;
(3), to remaining (n-2) satellite according to fault detect judgement amount
with fault detect thresholding σ
t, (n-2)carry out fault detect, if
confirm satellite i non-fault, make non-fault number of satellite add one, enter step (4); Otherwise, the combinations of satellites (i that record is rejected
k, j
k) and calculate position and the clock correction x of receiver by least square method
k, forward step (5) to;
(4), judge whether the non-fault number of satellite of current gained reaches 4, if reach 4 rough position and clock correction that calculate receiver
and identify fault satellites from residue (n-4) satellite, and again resolve receiver location and clock correction, detect and finish; Otherwise, enter step (5);
(5), make i=i+1, if i > n, non-fault number of satellite is less than 4, now, from reject combinations of satellites with pseudorange residual error statistic λ
kminimum value λ
min=min (λ
1, λ
2..., λ
m) corresponding combinations of satellites is as fault satellites, and calculate the receiver location corresponding with this combinations of satellites and clock correction as positioning result, detect and finish; Otherwise return to step (1).
Wherein, SSE is pseudorange residual sum of squares (RSS);
Fault detect thresholding σ
t, (n-1)obtained by following formula:
Wherein, icdf is (1-α) fractile that computer card side distributes;
Fault Identification judgement amount r
iobtained by following formula:
Wherein, SSE is pseudorange residual sum of squares (RSS);
Fault detect thresholding σ
t, (n-2)obtained by following formula:
Wherein, icdf is (1-α) fractile that computer card side distributes.
Wherein, H
0for 4 observing matrixes that non-fault satellite forms of current gained,
for H
0transposed matrix, y
0for the Pseudo-range Observations vector of 4 non-fault satellites of current gained.
Concrete, receiver location and clock correction x
kobtained by following formula:
Wherein, H
kfor having rejected two satellite (i
k, j
k) after the observing matrix that forms of remaining (n-2) satellite, y
kfor having rejected two satellite (i
k, j
k) after the Pseudo-range Observations vector of remaining (n-2) satellite;
Pseudorange residual error statistic λ
kobtained by following formula:
The invention has the beneficial effects as follows, solved the fault deviation cancellation problem of parity vector residual error, realized the detection identification of many fault satellites.Detect identification satellite fault by the present invention, can improve positioning precision.
Accompanying drawing explanation
Fig. 1 is process flow diagram of the present invention.
Embodiment
Below in conjunction with drawings and Examples, describe technical scheme of the present invention in detail.
The present invention is take fault satellites of parity vector method identification as basis, find out 4 non-fault satellites, then carry out coarse localization by the information of non-fault satellite, and identify thus fault satellites, after being rejected, re-start again location compute to improve positioning precision, thereby avoided being offset and failing to judge or erroneous judgement problem of bringing by fault deviation.
Of the present inventionly mainly contain two gordian techniquies, the one, the identification of non-fault satellite take fault satellites of parity vector method identification as basis, is found out 4 non-fault satellites in order to coarse localization from an existing n satellite; The 2nd, calculated rough position and the clock correction of receiver by the information of 4 non-fault satellites
wherein
the rough coordinates of receiver,
receiver clock correction.And by
calculate the pseudorange residual delta of all the other (n-4) satellites
i, and determine detection threshold according to the false-alarm probability of fault satellites identification, identify fault satellites.
Elaborate with regard to two gordian techniquies below:
1, non-fault satellite recognition principle
If there are two fault satellites to exist, after successively having rejected two satellites by parity vector method in n satellite, if be still greater than fault detect thresholding according to the fault detect judgement amount of remaining (n-2) satellite calculating gained, that is:
first satellite rejecting must be non-fault satellite, do not have deviation in pseudorange.The principle of identification non-fault satellite can be divided into following two kinds of situations and discuss:
(1) if be still less than fault detect thresholding according to the fault detect judgement amount of remaining (n-2) satellite calculating gained, that is:
time, in remaining (n-2) satellite, may not there is not fault satellites or have two fault satellites, but under the effect of offsetting in fault deviation, make
cannot judge in this case whether first rejected satellite is fault satellites.
(2) if be still greater than fault detect thresholding according to the fault detect judgement amount of remaining (n-2) satellite calculating gained, that is:
time, in remaining (n-2) satellite, may there is a fault satellites; Also now fault deviation counteracting is not obvious, has two fault satellites.No matter be which kind of situation, first satellite of rejecting must be all non-fault satellite, otherwise, in remaining (n-1) satellite, only have a fault satellites to exist, now adopt parity vector method another fault satellites can be identified, thereby
If it should be noted that note:
for event A; First satellite of rejecting is that non-fault satellite is event B.
Discuss for following three propositions:
Proposition 1: if B is A.This is a wrong proposition, when second satellite of rejecting is during also for satellite without reason, has two fault satellites, if the deviation that now breaks down is offset in remaining (n-2) satellite
therefore, even if first satellite is non-fault satellite,
also be likely less than σ
t (n-2).
Proposition 2: if
this proposition is the converse negative proposition of proposition 1, is also therefore wrong, and also just explanation can not be because of for it
just judge that two satellites of rejecting are fault satellites.
Proposition 3: if
even reject first satellite be fault satellites
this proposition is set up, in the time that first satellite of rejecting is fault satellites, in remaining (n-1) satellite, only there is a fault satellites, use parity vector method this satellite can be identified, so non-fault satellite exists in remaining (n-2) satellite, therefore
Method and proposition 3 that the present invention identifies non-fault satellite are mutually converse negative proposition, thereby the correctness of identification non-fault satellite method has been described.
2, fault satellites recognition principle
Can detect 4 non-fault satellites by method above, the information of 4 non-fault satellites just can calculate by the method for least square rough position and the receiver clock correction of receiver thus
Work as acquisition
after can compute pseudo-ranges residual delta
i:
In the time there is not pseudorange deviation in satellite i, pseudorange observational error
otherwise, e
ibe not 0 Gaussian distribution by obeying average.Therefore, as single satellite false-alarm probability α
1have to timing:
Through type (11) can calculate fault satellites identification thresholding T
Δ.Δ
iwith T
Δall identify for fault satellites.
Work as Δ
i> T
Δtime can judge that satellite i is as fault satellites, can realize thus the identification of two fault satellites.
In addition, when fault satellites is during more than one, due to the cause that fault deviation is offset, the Fault Identification judgement amount r of parity vector method
i(i=1,2 ..., maximal value r n)
j=max (r
1, r
2..., r
n), r
maxcorresponding satellite j may not be fault satellites.And from reducing the angle of calculated amount, in the time realizing, the rejecting of first satellite can adopt the method for rejecting one by one, in the time of second fault satellites of identification, re-use parity vector method.
From non-fault satellite recognition principle, the present invention is mainly by rejecting the fault detect judgement amount after two satellites
decide whether first satellite of rejecting is fault satellites.Calculating
time, need to meet n-2>=5, therefore, only have and in the time that number of satellite is more than or equal to 7 (n>=7), just can carry out double star fault detect.But when number of satellite is less and fault deviation is offset when serious, even if all satellites are all searched one time, the non-fault number of satellite obtaining may still be less than 4.Now can not calculate the coarse information of receiver, thereby cannot carry out fault satellites identification.
It is right that non-fault number of satellite is less than 4 o'clock needs
time all satellites of rejecting combination (i
k, j
k) test.Corresponding every kind of fault satellites combination, can be calculated receiver coordinate and the clock correction x of one group with residue (n-2) satellite by least square method
k=(x
k, y
k, z
k, b
k) k=1,2 ..., m, m is
corresponding all combinations of satellites numbers.
Equally according to x
kcan calculate the pseudorange residual error statistic λ of residue (n-2) satellite
k, as decision statistics:
If λ
l=min (λ
1, λ
2..., λ
m), (i
l, j
l) for fault satellites combination, after rejecting, position and resolve with residue satellite again, obtain receiver exact position.
Embodiment
Fig. 1 has provided flow process of the present invention.Its concrete grammar is as follows:
(1) when, reject satellite i from an existing n satellite, initial, i is 1;
Calculate the fault detect thresholding σ of (n-1) satellite
t, (n-1):
If
think that remaining non-fault satellite in satellite exists, and detects and finishes; Otherwise, carry out fault satellites identification by parity vector method, calculate Fault Identification judgement amount r
i:
Try to achieve r
i(i=1,2 ..., maximal value r n)
j=max (r
1, r
2..., r
n), after being rejected, satellite j enters step (3);
(3), remaining (n-2) satellite is carried out to fault detect, the fault detect judgement amount of calculating (n-2) satellite:
Wherein, SSE is pseudorange residual sum of squares (RSS);
Barrier detection threshold σ
t, (n-2)obtained by following formula:
Wherein, icdf is (1-α) fractile that computer card side distributes.
If
confirm satellite i non-fault, make non-fault number of satellite add one, enter step (4); Otherwise, the combinations of satellites (i that record is rejected
k, j
k) and calculate position and the clock correction x of receiver by least square method
k,
k=1,2 ..., N, wherein, H
kfor having rejected two satellite (i
k, j
k) after the observing matrix that forms of remaining (n-2) satellite, y
kfor having rejected two satellite (i
k, j
k) after the Pseudo-range Observations vector of remaining (n-2) satellite.Forward step (5) to;
(4), judge whether the non-fault number of satellite of current gained reaches 4, if reach 4 rough position and clock correction that calculate receiver
In formula (18), H
0for 4 observing matrixes that non-fault satellite forms of current gained,
for H
0transposed matrix, y
0for the Pseudo-range Observations vector of 4 non-fault satellites of current gained.
Obtain
after, can carry out compute pseudo-ranges residual delta according to formula (10) and formula (11)
iwith fault satellites identification thresholding T
Δ, and identify fault satellites from remaining (n-4) satellite, and again resolve receiver location and clock correction information, detect and finish; Otherwise enter step (5);
(5), make i=i+1, if i > n, non-fault number of satellite is less than 4, now, according to x
kobtain pseudorange residual error statistic λ
k, from reject combinations of satellites with pseudorange residual error statistic λ
kminimum value λ
min=min (λ
1, λ
2..., λ
m) corresponding combinations of satellites is as fault satellites, and calculate the receiver location corresponding with this combinations of satellites and clock correction as positioning result, detect and finish; Otherwise return to step (1).
Claims (5)
1. the double star fault recognition method based on parity vector method, comprises step:
(1), from an existing n satellite, reject satellite i, n is receiver visible satellite quantity in global position system, n >=7;
(2), to remaining (n-1) satellite according to fault detect judgement amount
with fault detect thresholding σ
t, (n-1)carry out fault detect, if
think that remaining non-fault satellite in satellite exists, and detects and finishes; Otherwise, carry out fault satellites identification by parity vector method, according to Fault Identification judgement amount r
imaximal value r
j=max (r
1, r
2, r
n), satellite j is rejected;
(3), to remaining (n-2) satellite according to fault detect judgement amount
with fault detect thresholding σ
t, (n-2)carry out fault detect, if
confirm satellite i non-fault, make non-fault number of satellite add one, enter step (4); Otherwise, the combinations of satellites (i that record is rejected
k, j
k) and calculate position and the clock correction x of receiver by least square method
k, forward step (5) to;
(4), judge whether the non-fault number of satellite of current gained reaches 4, if reach 4 rough position and clock correction that calculate receiver
and identify fault satellites from residue (n-4) satellite, and again resolve receiver location and clock correction, detect and finish; Otherwise, enter step (5);
(5), make i=i+1, if i > n, non-fault number of satellite is less than 4, now, from reject combinations of satellites with pseudorange residual error statistic λ
kminimum value λ
min=min (λ
1, λ
2..., λ
m) corresponding combinations of satellites is as fault satellites, and calculate the receiver location corresponding with this combinations of satellites and clock correction as positioning result, detect and finish; Otherwise return to step (1).
2. the double star fault recognition method based on parity vector method according to claim 1, is characterized in that fault detect judgement amount
obtained by following formula:
Wherein, SSE is pseudorange residual sum of squares (RSS);
Fault detect thresholding σ
t, (n-1)obtained by following formula:
Wherein, icdf is (1-α) fractile that computer card side distributes;
Fault Identification judgement amount r
iobtained by following formula:
3. the double star fault recognition method based on parity vector method according to claim 1, is characterized in that, detects judgement amount
obtained by following formula:
Wherein, SSE is pseudorange residual sum of squares (RSS);
Fault detect thresholding σ
t, (n-2)obtained by following formula:
Wherein, icdf is (1-α) fractile that computer card side distributes.
4. the double star fault recognition method based on parity vector method according to claim 1, is characterized in that rough position and clock correction
obtained by following formula:
5. the double star fault recognition method based on parity vector method according to claim 1, is characterized in that, receiver location and clock correction x
kobtained by following formula:
Wherein, H
kfor having rejected two satellite (i
k, j
k) after the observing matrix that forms of remaining (n-2) satellite, y
kfor having rejected two satellite (i
k, j
k) after the Pseudo-range Observations vector of remaining (n-2) satellite;
Pseudorange residual error statistic λ
kobtained by following formula:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210365524.1A CN102901971B (en) | 2012-09-27 | 2012-09-27 | Parity vector method-based double-satellite failure recognition method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210365524.1A CN102901971B (en) | 2012-09-27 | 2012-09-27 | Parity vector method-based double-satellite failure recognition method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102901971A CN102901971A (en) | 2013-01-30 |
CN102901971B true CN102901971B (en) | 2014-05-14 |
Family
ID=47574311
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210365524.1A Expired - Fee Related CN102901971B (en) | 2012-09-27 | 2012-09-27 | Parity vector method-based double-satellite failure recognition method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102901971B (en) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454650B (en) * | 2013-08-20 | 2015-06-24 | 北京航空航天大学 | Method for monitoring satellite integrity with vision as auxiliary |
CN104749587B (en) * | 2013-12-31 | 2017-03-29 | 清华大学 | Receiver pseudorange fault monitoring method and receiver |
CN106646526B (en) * | 2017-02-09 | 2019-05-31 | 南京航空航天大学 | A kind of receiver-autonomous integrity detection method that can detect identification various faults simultaneously |
CN107783154B (en) * | 2017-09-22 | 2019-07-23 | 北京时代民芯科技有限公司 | A kind of receiver-autonomous integrity fault detection and method for removing |
CN108507590B (en) * | 2018-03-20 | 2023-03-07 | 千寻位置网络(浙江)有限公司 | Constant speed evaluation method and system and vehicle-mounted terminal |
CN109709583B (en) * | 2018-11-16 | 2022-08-05 | 南京航空航天大学 | Composite fault self-adaptive detection method for multiple-loop iteration sliding window accumulation |
CN110596736B (en) * | 2019-10-15 | 2021-04-02 | 中国电子科技集团公司第五十四研究所 | GNSS observation abnormal value detection and isolation method |
CN110907953B (en) * | 2019-10-18 | 2022-04-29 | 湖北三江航天险峰电子信息有限公司 | Satellite fault identification method and device and software receiver |
CN111123304B (en) * | 2019-11-28 | 2021-12-24 | 北京航空航天大学 | Visual navigation integrity monitoring and calculating method |
CN113777630B (en) * | 2021-02-10 | 2022-05-27 | 北京航空航天大学 | Fault monitoring method and system for single receiver of ground-based augmentation system |
CN113341438B (en) * | 2021-06-02 | 2023-03-24 | 成都天奥信息科技有限公司 | Multi-satellite fault identification method and system based on gross error inverse solution |
CN115390099B (en) * | 2022-10-28 | 2023-02-03 | 北京航空航天大学 | Maximum value joint chi-square fault elimination method based on odd-even vector projection |
CN115561782B (en) * | 2022-11-18 | 2023-02-28 | 北京航空航天大学 | Satellite fault detection method in integrated navigation based on odd-even vector projection |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1996041A (en) * | 2006-12-20 | 2007-07-11 | 北京航空航天大学 | Method for monitoring GNSS receiver autonomous integrity based on multi-satellite failure recognition |
CN101520503A (en) * | 2009-03-19 | 2009-09-02 | 北京航空航天大学 | Method for detecting fault satellite of satellite navigation system |
-
2012
- 2012-09-27 CN CN201210365524.1A patent/CN102901971B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1996041A (en) * | 2006-12-20 | 2007-07-11 | 北京航空航天大学 | Method for monitoring GNSS receiver autonomous integrity based on multi-satellite failure recognition |
CN101520503A (en) * | 2009-03-19 | 2009-09-02 | 北京航空航天大学 | Method for detecting fault satellite of satellite navigation system |
Non-Patent Citations (6)
Title |
---|
GNSS 中基于双星故障的RAIM 中的数学关系;陈灿辉 等;《遥测遥控》;20110331;第32卷(第2期);第33-38页 * |
卫星故障诊断的最优奇偶向量法;杨静 等;《航空学报》;20020331;第23卷(第2期);第183-186页 * |
杨静 等.卫星故障诊断的最优奇偶向量法.《航空学报》.2002,第23卷(第2期),第183-186页. |
用于检测双星故障的RA IM算法分析与研究;郭承军 等;《计算机应用研究》;20110331;第28卷(第3期);第925-926,944页 * |
郭承军 等.用于检测双星故障的RA IM算法分析与研究.《计算机应用研究》.2011,第28卷(第3期),第925-926,944页. |
陈灿辉 等.GNSS 中基于双星故障的RAIM 中的数学关系.《遥测遥控》.2011,第32卷(第2期),第33-38页. |
Also Published As
Publication number | Publication date |
---|---|
CN102901971A (en) | 2013-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102901971B (en) | Parity vector method-based double-satellite failure recognition method | |
CN100582811C (en) | Method for monitoring GNSS receiver autonomous integrity based on multi-satellite failure recognition | |
CN102654407B (en) | Multiple-fault detecting device and detecting method for tightly-integrated inertial satellite navigation system | |
US10670729B2 (en) | System and method to provide an ASIL qualifier for GNSS position and related values | |
CN101281248B (en) | Multi-fault recognizing method applied to combined satellite navigation system | |
US20110118979A1 (en) | Automotive location data integrity | |
CN101799524B (en) | Method for autonomously monitoring receiver integrity of global navigation satellite system | |
CN106707304A (en) | Satellite navigation receiver fault satellite detection method | |
CN107783154B (en) | A kind of receiver-autonomous integrity fault detection and method for removing | |
CN102135621B (en) | Fault recognition method for multi-constellation integrated navigation system | |
CN105487088A (en) | RAIM algorithm in satellite navigation system based on Kalman filtering | |
CN105738925A (en) | Method for monitoring satellite receiver autonomous integrity special for train positioning | |
CN104199051B (en) | Method for detecting and identifying satellite navigation RAIM (Receiver Autonomous Integrity Monitoring) multi-satellite faults | |
CN101806903A (en) | Receiver autonomous integrity monitoring (RAIM) method used for satellite navigation system | |
CN104267410A (en) | Method and device for excluding multiple faults in airborne integrity monitoring | |
CN105549033A (en) | Integrity processing method based on least square residual error edge detection | |
CN115265594B (en) | Multi-source PNT information elastic fusion navigation multi-level autonomous integrity monitoring method and system | |
CN105116391A (en) | Marine radar error calibration-oriented valid target sequence set combined analysis method | |
CN103983986A (en) | Particle filter-based improved RAIM (Receiver Autonomous Integrity Monitoring) anti-deception jamming method | |
CN115420284A (en) | Fault detection and identification method for combined navigation system | |
CN105511481A (en) | Satellite borne orbit determination optimization method | |
CN116859415A (en) | Quick, stable and high-precision multi-fault satellite identification and positioning method | |
CN113341438B (en) | Multi-satellite fault identification method and system based on gross error inverse solution | |
CN115062097A (en) | Positioning integrity monitoring method, terminal and computer storage medium | |
CN110907953B (en) | Satellite fault identification method and device and software receiver |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20140514 Termination date: 20150927 |
|
EXPY | Termination of patent right or utility model |