Next Article in Journal
Semiconductor Optical Amplifiers with Wide Gain Bandwidth and Enhanced Polarization Insensitivity Based on Tensile-Strained Quantum Wells
Next Article in Special Issue
Deep Learning Soft-Decision GNSS Multipath Detection and Mitigation
Previous Article in Journal
Advancements in Detection and Mitigation Strategies for Petroleum-Derived Contaminants in Aquatic Environments: A Comprehensive Review
Previous Article in Special Issue
LJCD-Net: Cross-Domain Jamming Generalization Diagnostic Network Based on Deep Adversarial Transfer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved RAIM Availability Assessment Method Based on the Characteristic Slope

1
China Transport Telecommunications and Information Center, Beijing 100011, China
2
Graduate School, Beihang University, Beijing 100191, China
3
School of Electronic Information Engineering, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(11), 3283; https://doi.org/10.3390/s24113283
Submission received: 15 April 2024 / Revised: 16 May 2024 / Accepted: 17 May 2024 / Published: 21 May 2024

Abstract

:
The availability assessment is an important step for onboard application in Receiver Autonomous Integrity Monitoring (RAIM)s. It is commonly implemented using the protection level (PL)-based method. This paper analyzes the deficiencies of three kinds of PL-based methods: RAIM availability might be optimistically or conservatively assessed using the classic-PL-base method; might be conservatively assessed using the enhanced-PL-based method, and neither be optimistically nor conservatively assessed using the ideal-PL-based method with the cost of large calculation amount on-board. An improved slope-based RAIM availability assessment method is proposed, in which the characteristic slope is designed as the assessment basis, and its threshold that can exactly match the integrity risk requirement is derived. The slope-based method has the same RAIM availability assessment result as the ideal-PL-based method. Moreover, because the slope threshold can be calculated offline and searched online, the on-board calculation burden can be reduced using the slope-based method. Simulation is presented to verify the theoretical analysis of the RAIM availability assessment performances for the three PL-based and the slope-based methods.

1. Introduction

The integrity of the Global Navigation Satellite System (GNSS) is one of the important factors to ensure civil aviation safety. There are three categories of GNSS integrity augmentation systems: Satellite-Based Augmentation System (SBAS), Ground-Based Augmentation System (GBAS) and Aircraft-Based Augmentation System (ABAS). The first two categories are at the system level, and the latter category is at the user level [1]. ABAS can be implemented with Receiver Autonomous Integrity Monitoring (RAIM), which provides a navigation solution with guaranteed integrity by consistency checking among measurements [2].
For an onboard application, RAIM is executed in two steps: the RAIM availability assessment and the satellite fault detection [1]. The former (RAIM availability assessment) is used to assess in advance whether the navigation solution can meet the integrity risk requirement with the fault detection procedure. For decades, RAIM availability assessment has been achieved by calculating the protection level (PL), the upper bound of the position error corresponding to the integrity risk requirement [3]. The threshold of the PL is the alert limit (AL), the upper bound of the user-allowed position error. If the PL is lower than the AL, RAIM is considered available; otherwise, it is considered unavailable.
Many studies have focused on determining how to obtain a lower PL to improve the availability of RAIM. Some of these studies were devoted to developing new navigation solution calculation methods, for example, the improved Integrity-Optimized RAIM (NIORAIM) [4] and the optimal weighted average solution (OWAS) [5] methods used for the snapshot RAIM algorithm. These methods can obviously decrease the PL with a slight increase in nominal position error. In addition, some studies have committed to accurately modeling the stochastic measurement noise, such as the discrete error-distribution (NavDEN) model proposed by Rife and Pervan [6] and the distribution model considering both elevation angle and orbit type proposed by Fan [7]. These measurement noise models are all helpful for obtaining a tight PL.
However, in the process of pursuing a lower PL, i.e., higher availability of RAIM, there is a key issue that is ignored by most researchers: whether the PL can accurately assess RAIM availability. Milner and Ochieng noted this issue [3]. They qualitatively described that the classic PL, the product of the characteristic slope and the Minimum Detection Bias (MDB) proposed by Brown and Chin [8], was too optimistic for RAIM availability assessment. The slope is a geometric feature-related parameter that qualitatively describes the relationship between positioning error and pseudo-range residual [9]. The reason is “PL < AL” might not mean that the integrity risk satisfies the requirement for the measurement bias less than MDB. Meanwhile, its enhancement, abbreviated as the enhanced PL in the following, which provides an additional term to protect against the variation in position error proposed by Angus [10], is too conservative. The reason is “PL ≥ AL” might not mean that the integrity risk exceeds the requirement. Here, measurement bias means the measurement error caused by the satellite fault, which is different from the measurement noise in the nominal mode. Furthermore, they proposed the ideal PL. It is the minimum PL value that guarantees the integrity risk, satisfying its requirement for arbitrary measurement bias. The ideal PL can prevent RAIM availability assessment from being optimistic or conservative. However, it cannot be solved analytically. A numerical search for the ideal PL begins with an improbably large value [3], which leads to a large amount of calculation, increasing the computational burden of a GNSS receiver or an onboard computer.
In recent years, most researchers have focused on Advanced RAIM (ARAIM), in which PL is calculated after fault detection [11]. ARAIM supports multi-constellation dual-frequency GNSS integrity monitoring. Multiple hypothesis solution separation (MHSS) algorithm is used in ARAIM. How to solve the accurate PL for MHSS is a research hotspot, including the PL calculation method for each fault mode and the optimization strategy for the ultimate PL [12,13,14,15,16]. Jiang and Wang adopted the ideal PL [3] in ARAIM and verified it was more accurate than other PLs for the ARAIM availability assessment [17,18]. ARAIM is still in the theoretical research stage and is not currently being applied in engineering practice.
Compared with ARAIM, RAIM has two deficiencies. The first is that RAIM is designed for a single constellation, monitoring only a single satellite fault [2]. ARAIM is designed for double constellations, monitoring not only the single satellite fault but also the multiple satellite faults and the constellation fault [12]. The second is classic, and the enhanced PLs are not rigorous enough for RAIM availability assessment, while the PL of ARAIM is much more rigorous. However, the on-board calculation of RAIM using the classic or the enhanced PLs is much less than that of ARAIM. For a single constellation, RAIM can still be used, but there is a problem needs to be considered, founding a RAIM availability assessment method both satisfying the rigor and maintaining the low on-board computational burden.
In this paper, a slope-based RAIM availability assessment method is proposed to solve the above problem. The characteristic slope is taken as the assessment basis. Using the ideal slope threshold, this method can achieve a consistent RAIM availability assessment with the ideal-PL-based method. The ideal slope threshold can be calculated offline and searched online because it is only related to one geometric parameter.
The remainder of this paper is organized as follows: Section 2 states some technical backgrounds. Section 3 reviews PL-based RAIM availability assessment methods, including the classic, enhanced and ideal PLs. The deficiencies of the classic and enhanced-PL-based methods can be analyzed quantitatively using the rates of optimistic or conservative assessment. Section 4 proposes the slope-based RAIM availability assessment method after the derivation of the ideal slope threshold. Section 5 gives an overview of the simulation results of the classic-PL-based, enhanced-PL-based, ideal-PL-based and slope-based methods. Section 6 concludes this work with a brief summary. The discussion of this paper takes vertical integrity as an example, using one kind of classic RAIM snapshot algorithm, i.e., the least squares residuals (LSR) algorithm. In this paper, the measurement noise is assumed to be independent white Gaussian noise (WGN).

2. Technical Background

Before discussing the RAIM availability assessment methods, some technical backgrounds need to be stated, including the derivation of the integrity risk requirement for the single-satellite fault mode, the definition of the vertical characteristic slope and the specific meaning of RAIM being available.

2.1. Integrity Risk Requirement for the Single-Satellite Fault Mode

The integrity risk P ( H M I ) [19] is the probability of undetected faults causing unacceptably large errors in the estimated position [20]. HMI is short for hazardous misleading information (HMI). P ( H M I ) can be divided into three fault modes, the nominal mode, the single-satellite fault mode and the multiple-satellite fault mode, expressed as:
P ( H M I ) = i = 0 2 P ( H M I , i F ) = i = 0 2 P ( H M I i F ) P ( i F )
In Equation (1), 0 F , 1 F and 2 F respectively represent the nominal, the single-satellite fault and the multiple-satellite fault modes, where P ( i F ) is the prior probability of the fault mode i F and P ( H M I i F ) is the probability of HMI under the i F fault mode. P ( 0 F ) = ( 1 P s a t ) K and P ( 1 F ) = C K 1 P s a t ( 1 P s a t ) K , where C K 1 means the number of combinations for choosing one element from K elements; P s a t is the prior fault probability of a satellite, and K is the total visible satellite number. Taking vertical plane for example, P ( H M I i F ) is calculated with the following equation [21]:
P ( H M I i F ) = P ( V P E V A L i F ) P ( T s < T d i F )
In Equation (2), V P E is the vertical position error; V A L is the vertical AL; T s and T d are the fault detection test statistic and threshold, respectively.
To ensure the integrity of the navigation system, P ( H M I ) should be less than its requirement, denoted as P r ( H M I ) .
P ( H M I , 0 F ) + P ( H M I , 1 F ) + P ( H M I , 2 F ) < P r ( H M I )
Given a geometry between the user and all-in-view satellites, V A L , α and P s a t , P ( H M I , 0 F ) can be calculated following:
P ( H M I , 0 F ) = ( 1 α ) [ 1 V A L V A L g ( x ; 0 , a v 2 ) d x ] ( 1 P s a t ) K
In Equation (4), α is the allowable false alarm probability under the nominal fault mode, satisfying P ( T s < T d 0 F ) = 1 α ; g ( x ; 0 , a v 2 ) is the probability density function (PDF) of the normal-distributed VPE under the nominal fault mode with mean value 0 and standard deviation value a v , where a v is explained in Appendix A.
According to GNSS Evolutionary Architecture Study (GEAS) report [22], the integrity risk requirement allocated on the multiple satellites fault, denoted as P r ( H M I , 2 F ) , can be set 1.3 × 1 0 8 [22].
To ensure Equation (3) to be true, the integrity risk for single fault mode should satisfy:
P ( H M I , 1 F ) < P r ( H M I 1 F )
where P r ( H M I 1 F ) = [ P r ( H M I ) P r ( H M I , 2 F ) P ( H M I , 0 F ) ] / P ( 1 F ) . Thus, the problem of P ( H M I ) < P r ( H M I ) has evolved into the problem of P ( H M I 1 F ) < P r ( H M I 1 F ) .
Moreover, to ensure that P r ( H M I , 1 F ) is nonnegative, the threshold for a v , denoted as T a v , can be derived from the inequality P r ( H M I ) P r ( H M I , 2 F ) P ( H M I , 0 F ) 0 .
T a ν = V A L / Φ 1 P r ( H M I ) P r ( H M I , 2 F ) 2 ( 1 α ) ( 1 P s a t ) K
where Φ 1 · represents the inverse function of the cumulative distribution function (CDF) for the standard normal distribution. a v T a v indicates that only the integrity risk of the nominal and multiple-satellite fault modes have exceeded the total requirement, i.e., R 0 + P r ( H M I , 2 ) P r ( H M I ) .

2.2. Vertical Characteristic Slope

The vertical characteristic slope is defined according to these two parameters [2],
S l o p e m = a 3 m 2 s m m
where a 3 m and s m m respectively characterizes V P E and T s change caused by the measurement bias of the m-th visible satellite, signed as V S m , is the faulty satellite. The details of a 3 m and s m m can be seen in Appendix A and Appendix B, respectively. For a specific bias value, a faulty satellite with a large slope value will present a high P ( H M I 1 F ) , and a faulty satellite with a small slope value will present a low P ( H M I 1 F ) . The “characteristic slope” will be abbreviated as “slope” hereafter.

2.3. Specific Meaning of RAIM Being Available

RAIM being available refers to P ( H M I ) < P r ( H M I ) , i.e., P ( H M I 1 F ) < P r ( H M I 1 F ) for the arbitrary measurement bias value while RAIM being unavailable refers to P ( H M I 1 F ) P r ( H M I 1 F ) for at least one measurement bias.
A specific example is used to intuitively explain the meaning of RAIM being available. The 32-satellite GPS constellation is used in this example. The pseudorange measurement is assumed to be the dual-frequency ionosphere-free combination of L1 and L5. The standard deviation of the measurement noise for V S i , signed as σ i , is set according to the ARAIM interim report [11]. For the location of 37° N latitude, 117° E longitude and height 0 m and the epoch of UTC 14 March 2019 17:15:00, there are 9 visible satellites with a masking angle of 10°. Their vertical slope values are recorded in Table 1.
Figure 1 presents the base-10 logarithm of P ( H M I 1 F ) for PRN5, PRN6 and PRN13 with the measurement bias ξ b in the interval of ( 0 .   30   m ) . These three satellites have the top three vertical slope values, as shown in Table 1. The ξ b lg P ( H M I 1 F ) curves for these three satellites follow the same order as the slope values, which illustrates that a faulty satellite with a large slope will present a high P ( H M I 1 F ) . The ξ b lg P ( H M I 1 F ) curves for the other six visible satellites must be lower than that of PRN 13 because their vertical slopes are smaller. In this example, RAIM would be unavailable if PRN5 was the faulty satellite because P ( H M I 1 F ) is larger than P r ( H M I 1 F ) for ξ b values in the range from 11 m to 15 m. Therefore, the intersection between the ξ b lg P ( H M I 1 F ) curve of the faulty satellite and the P r ( H M I 1 F ) line means that RAIM is unavailable. RAIM would be available if the faulty satellite was one of the other visible satellites except for PRN5, because its P ( H M I 1 F ) smaller than P r ( H M I 1 F ) at an arbitrary measurement bias. Therefore, the separation between the ξ b lg P ( H M I 1 F ) curve of the faulty satellite and the P r ( H M I 1 F ) line means that RAIM is available.
Because the faulty satellite is unknown in actual situations, RAIM is considered available only if P ( H M I 1 F ) < P r ( H M I 1 F ) is true for the arbitrary measurement bias in the worst case, i.e., the satellite with the maximum slope being faulty.

3. PL-Based RAIM Availability Assessment

The classic PL, enhanced PL, and ideal PL are reviewed in this section.

3.1. Classic PL and Its Enhancement

The classic PL, denoted as P L c , is defined as the product of the maximum slope S l o p e m a x and the minimum detectable bias λ a as follows:
P L c = S l o p e m a x λ a
where λ a is the noncentral parameter of the fault detection test statistic T s under single-satellite fault mode [23].
The enhanced PL, denoted as P L e , has an additional term that protects against variation in the random error of the position solution on the basis of P L c as follows [10]:
P L e = S l o p e m a x λ a + α ( P M D ) σ
where α ( P M D ) = Φ 1 ( 1 P M D ) and σ is the standard deviation of the position error distribution.
Slope reflects the relationship between the position error and the pseudorange residual. PL is the projection of the pseudorange residual on the position error as shown in Equations (8) and (9). According to λ = ξ b σ m 2 2 s m m derived in Appendix B, the measurement bias ξ b for a faulty satellite V S m , which causes T s to obey χ 2 ( K 4 , λ a ) , is ξ b = σ m 2 λ a s m m . Based on the position error derivation in Appendix A, this ξ b makes the VPE under single-satellite fault mode obey:
V P E N ( S l o p e m λ a , a v 2 )
Contrasting Equation (8) with Equation (10), the classic vertical PL, denoted as V P L c , is the expectation of V P E at ξ b = σ m 2 λ a s m m with the maximum slope. Therefore, P ( V P E V P L c 1 F ) = 0.5 at ξ b = σ m 2 λ a s m m in the worst case. Similarly, the enhanced vertical PL, denoted as V P L e , is the upper quantile of the P M D of V P E distribution at ξ b = σ m 2 λ a s m m with the maximum slope. Therefore, P ( V P E V P L e 1 F ) = P M D at ξ b = σ m 2 λ a s m m in the worst case.

3.2. Deficiency of Classic- and Enhanced-PL-Based RAIM Availability Assessment

Here defines H0 as “ ξ b , P ( H M I 1 F ) < P r ( H M I 1 F ) ” and H1 as “ ξ b , P ( H M I 1 F ) P r ( H M I 1 F ) ” in the worst case. Therefore “ V P L < V A L H 1 ” or “ V P L V A L H 0 ” respectively mean the RAIM availability assessment is optimistic or conservative.
For ξ b
V P L < V A L P V P E V P L 1 F P T s < T d 1 F > P ( V P E V A L 1 F ) P ( T s < T d 1 F ) P a ( H M I 1 F ) > P ( H M I 1 F )
and
V P L V A L P a ( H M I 1 F ) P ( H M I 1 F )
where P ( V P E V P L 1 F ) P ( T s < T d 1 F ) = P a ( H M I 1 F ) ,
Similar to P ( H M I 1 F ) , P a ( H M I 1 F ) changes with the ξ b value. Moreover, the value of P a ( H M I 1 F ) is known at ξ b = σ m 2 λ a s m m , P a ( H M I 1 F ) = 0.5 P M D for V P L c and P a ( H M I 1 F ) = P M D 2 for V P L e . A ξ b lg P a ( H M I 1 F ) curve can be generated, which is higher than the ξ b lg P ( H M I 1 F ) curve with V P L < V A L and lower than the ξ b lg P ( H M I 1 F ) curve with V P L V A L . The positional relationship among the ξ b lg P ( H M I 1 F ) curve, the ξ b lg P a ( H M I 1 F ) curve, and the lg P r ( H M I 1 F ) line jointly determine whether the PL-based RAIM availability assessment is conservative or optimistic. The following is a RAIM availability assessment discussion according to the positional relationship between the ξ b lg P a ( H M I 1 F ) curve and the lg P r ( H M I 1 F ) line.

3.2.1. Separation

In this situation, there are three kinds of positional relationships among the ξ b lg P ( H M I 1 F ) curve, the ξ b lg P a ( H M I 1 F ) curve and the lg P r ( H M I 1 F ) line, as shown in Figure 2.
The ξ b lg P ( H M I 1 F ) curve intersects with the lg P r ( H M I 1 F ) line, and is higher than the blue dashed ξ b lg P a ( H M I 1 F ) curve, shown as the brown P 1 ( H M I 1 F ) curve in Figure 2. In this situation, RAIM is unavailable and the VPL is larger than the VAL, i.e., V P L V A L H 1 , meaning a successful detection of “RAIM being unavailable”.
The ξ b lg P ( H M I 1 F ) curve is separated from the lg P r ( H M I 1 F ) line, and is higher than the blue dashed curve ξ b lg P a ( H M I 1 F ) curve, shown as the black curve P 2 ( H M I 1 F ) in Figure 2. In this situation, RAIM is available but the VPL is larger than the VAL, i.e., V P L V A L H 0 , meaning a conservative assessment.
The ξ b lg P ( H M I 1 F ) curve is separated from the lg P r ( H M I 1 F ) line, and is lower than the blue dashed ξ b lg P a ( H M I 1 F ) curve, shown as the blue curve in Figure 2. In this situation, RAIM is available, and the VPL is smaller than the VAL, i.e., V P L < V A L H 0 , meaning a successful detection of “RAIM being available”.
Therefore, for the condition that the ξ b lg P a ( H M I 1 F ) curve is separated from the lg P r ( H M I 1 F ) line, the RAIM availability might be conservatively assessed.

3.2.2. Intersection

In this situation, there are also three kinds of positional relationships among the ξ b lg P ( H M I 1 F ) curve, the ξ b lg P a ( H M I 1 F ) curve and the lg P r ( H M I 1 F ) line, as shown in Figure 3.
The ξ b lg P ( H M I 1 F ) curve intersects with the lg P r ( H M I 1 F ) line, and is higher than the blue dashed ξ b lg P a ( H M I 1 F ) curve, shown as the brown P 1 ( H M I 1 F ) curve in Figure 3 In this situation, RAIM is unavailable and the VPL is larger than the VAL, i.e., V P L V A L H 1 , meaning a successful detection of “RAIM being unavailable”.
The ξ b lg P ( H M I 1 F ) curve intersects with the lg P r ( H M I 1 F ) line, and is lower than the blue dashed ξ b lg P a ( H M I 1 F ) curve, shown as the black P 2 ( H M I 1 F ) curve in Figure 3. In this situation, RAIM is unavailable, but the VPL is smaller than the VAL, i.e., V P L < V A L H 1 , meaning an optimistic assessment.
The ξ b lg P ( H M I 1 F ) curve is separated from the lg P r ( H M I 1 F ) line, and is lower than the blue dashed ξ b lg P a ( H M I 1 F ) curve, shown as the blue curve P 3 ( H M I 1 F ) in Figure 3. In this situation, RAIM is available, and the VPL is smaller than the VAL, i.e., V P L < V A L H 0 , meaning the successful detection of “RAIM being available”.
Therefore, for the condition that the ξ b lg P a ( H M I 1 F ) curve intersects with the lg P r ( H M I 1 F ) line, the RAIM availability might be optimistically assessed.

3.2.3. Tangency

In this situation, there are two kinds of positional relationships among the ξ b lg P ( H M I 1 F ) curve, the ξ b lg P a ( H M I 1 F ) curve and the lg P r ( H M I 1 F ) line, as shown in Figure 4.
The ξ b lg P ( H M I 1 F ) curve intersects with the lg P r ( H M I 1 F ) line, and is higher than the blue dashed ξ b lg P a ( H M I 1 F ) curve, shown as the brown P 1 ( H M I 1 F ) curve in Figure 4. In this situation, RAIM is unavailable, and the VPL is larger than the VAL, i.e., V P L V A L H 1 , meaning a successful detection of “RAIM being unavailable”.
The ξ b lg P ( H M I 1 F ) curve is separated from the lg P r ( H M I 1 F ) line, and is lower than the blue dashed ξ b lg P a ( H M I 1 F ) curve, shown as the black P 2 ( H M I 1 F ) curve in Figure 4. In this situation, RAIM is available and the VPL is smaller than the VAL, i.e., V P L < V A L H 0 , meaning a successful detection of “RAIM being available”.
Therefore, for the condition that the ξ b lg P a ( H M I 1 F ) curve is tangent to the lg P r ( H M I 1 F ) line, both optimistic and conservative assessments can be prevented.
Based on the above analysis, the accuracy of PL-based RAIM availability assessment depends on the positional relationship between the ξ b lg P a ( H M I 1 F ) curve and the lg P r ( H M I 1 F ) line. There is a risk of conservative assessment when the ξ b lg P a ( H M I 1 F ) curve is separated from the lg P r ( H M I 1 F ) line and risk of optimistic assessment when the ξ b lg P a ( H M I 1 F ) curve intersects with the lg P r ( H M I 1 F ) line. PL-based RAIM availability assessment is accurate only if the ξ b lg P a ( H M I 1 F ) curve is tangent to the lg P r ( H M I 1 F ) line. Because only the P a ( H M I 1 F ) value at a specific ξ b , i.e., ξ b = σ m 2 λ a s m m , is determined for V P L c or V P L e , the position relationship between the entire ξ b lg P a ( H M I 1 F ) curve and the lg P r ( H M I 1 F ) line is uncertain. Consequently, both optimistic and conservative assessments might happen when using V P L c or V P L e to assess whether vertical RAIM is available.
Because V P L c < V P L e , the ξ b lg P a ( H M I 1 F ) curve for V P L c is much higher than that for V P L e . Thus the possibility of intersection between the ξ b lg P a ( H M I 1 F ) curve and the lg P r ( H M I 1 F ) line for V P L c is much higher than that for V P L e , which may lead to an optimistic assessment, while the possibility of separation between the ξ b lg P a ( H M I 1 F ) curve and the lg P r ( H M I 1 F ) line for V P L e is much higher than that for P L c , which may lead to conservative assessment. Consequently, the optimistic assessment risk of using V P L c is higher than that of using V P L e for vertical RAIM availability; in contrast, the conservative assessment risk of using V P L e is higher than that of using V P L c .

3.3. Ideal Protection Level

According to the above analysis, the ideal positional relationship between the ξ b lg P a ( H M I 1 F ) curve and the lg P r ( H M I 1 F ) line is tangency, which can prevent both optimistic and conservative RAIM availability assessments. The ideal VPL, denoted as V P L d , proposed by Milner and Ochieng, satisfies this condition. It matches the exact required integrity risk for the worst-case bias (WCB), the measurement bias presenting the highest integrity risk. Thus V P L d forms a ξ b lg P a ( H M I 1 F ) curve tangent to the lg P r ( H M I 1 F ) line. If V P L d < V A L , the ξ b lg P ( H M I 1 F ) curve must be separated from the lg P r ( H M I 1 F ) line, which means that vertical RAIM is available; otherwise the ξ b lg P ( H M I 1 F ) curve must be tangent to or intersect with the lg P r ( H M I 1 F ) line, which means that vertical RAIM is unavailable.
V P L d is the solution of:
max [ P a ( H M I 1 F ) ] = P r ( H M I 1 F )

4. Slope-Based RAIM Availability Assessment

In addition to the ideal PL, there is another ideal test statistic for RAIM availability assessment: the slope. Both optimistic and conservative RAIM availability assessments can be prevented using the slope once an ideal threshold is found. The following is a deviation of this ideal threshold.

4.1. Derivation of the Ideal Threshold for the Slope

The ideal slope threshold derivation begins with searching for a condition satisfying P ( H M I 1 F ) < P r ( H M I 1 F ) for all possible measurement bias values of an arbitrary faulty satellite. To ensure P ( H M I 1 F ) < P r ( H M I 1 F ) constantly true, the maximum value of P ( H M I 1 F ) should be less than P r ( H M I 1 F ) ,
max [ P ( H M I 1 F ) ] < P r ( H M I 1 F )
Where:
P ( H M I 1 F = P ( V P E V A L 1 F P ( T s < T d 1 F = [ 1 V A L V A L g ( x ; a 3 m ξ b σ m 2 , a v 2 ) d x ] 0 T d f ( y ; ξ b σ m 2 2 s m m , K 4 ) d y
Substituting Equation (7) into Equation (15),
P ( H M I 1 F ) = [ 1 V A L V A L g ( x ; a 3 m ξ b σ m 2 , a v 2 ) d x ] 0 T d f ( y ; ξ b σ m 2 2 a 3 m S l o p e m 2 , K 4 ) d y )
Setting t = x a 3 m ξ b / σ m 2 a v , Equation (16) can be transformed into:
P ( H M I 1 F ) = [ 1 V A L a 3 m ξ b σ m 2 a v V A L a 3 m ξ b σ m 2 a v g ( t ; 0,1 ) d t ] 0 T d f ( y ; a 3 m ξ b σ m 2 2 1 S l o p e m 2 , K 4 ) d y )
Setting μ s = a 3 m ξ b σ m 2 , P ( H M I 1 F ) can be taken as a function of μ s as follows:
Φ ( μ s ) = P ( H M I 1 F ) = [ 1 V A L μ s a v V A L μ s a v g ( t ; 0,1 ) d t ] 0 T d f ( y ; μ s S l o p e m 2 , K 4 ) d y )
When substituting Equation (18) into Equation (14), an ideal threshold must exist for S l o p e m with the given P r ( H M I 1 F ) , V A L , a v , T d and K , denoted as T S l o p e , which satisfies the limit situation
max Φ ( μ s ) S l o p e m = T S l o p e = P r ( H M I 1 F )
According to Equations (18) and (19), once P r ( H M I 1 F ) , V A L , T d and K are given, T S l o p e is determined by only one parameter a v .
Figure 5 presents the numerically solved T S l o p e for a v values with a step of 0.001 at P r ( H M I ) = 1 × 1 0 7 , P r ( H M I , 2 ) = 1.3 × 1 0 8 , V A L = 50   m , P s a t = 1 × 1 0 5 , α = 1 × 1 0 6 and K = 7 , , 11 . As shown in Figure 5, T S l o p e decreases with the increase of a v . Each a v T S l o p e curve exhibits a nearly constant segment at the beginning and a sharply decreasing segment at the end. The a v T S l o p e curve ends when a v reaches T a v , and the spacing between two adjacent curves obviously decreases as K increases.

4.2. Practical Meaning of the Ideal Slope Threshold

Here, the specific example in Section 2.3 is used to intuitively explain the practical meaning of T S l o p e . The a v value for this example is 4.944. Given P r ( H M I ) = 1 × 1 0 7 , P r ( H M I , 2 F ) = 1.3 × 1 0 8 , V A L = 50   m , P s a t = 1 × 1 0 5 and α = 1 × 1 0 6 , the ideal slope threshold can be numerically solved according to Equation (19), and T S l o p e = 5.390 . Figure 6 presents the μ s lg P ( H M I | 1 F ) curves for S l o p e = 6.060 (PRN5), S l o p e = 3.010 (PRN6) and T S l o p e = 5.390 . P ( H M I | 1 F ) for T S l o p e is calculated according to Equation (19). The μ s lg P ( H M I | 1 F ) curve for T S l o p e is tangent to the lg P r ( H M I 1 F ) line. The μ s lg P ( H M I | 1 F ) curve for S l o p e = 6.060 , which is larger than T S l o p e , is intersected with the lg P r ( H M I 1 F ) line. In contrast, the μ s lg P ( H M I | 1 F ) curve for S l o p e = 3.010 , which is smaller than T S l o p e , is separated from the lg P r ( H M I 1 F ) line.
It can be deduced that the slope value determines the positional relationship between the μ s lg P ( H M I | 1 F ) curve and the lg P r ( H M I 1 F ) line. T S l o p e forms a tangent μ s lg P ( H M I | 1 F ) curve to the lg P r ( H M I 1 F ) line. If the slope value of a faulty satellite is larger than T S l o p e , its μ s lg P ( H M I | 1 F ) curve would be intersected with the lg P r ( H M I 1 F ) line, meaning that P ( H M I 1 F ) P r ( H M I 1 F ) can be satisfied for at least one possible measurement bias value. If the slope value of a faulty satellite was smaller than T S l o p e , its μ s lg P ( H M I | 1 F ) curve would be separated from the lg P r ( H M I 1 F ) line, meaning that P ( H M I 1 F ) < P r ( H M I 1 F ) can be ensured at an arbitrary measurement bias value. From the analysis of this specific example, the slope is an ideal test statistic for RAIM availability assessment, with the ideal threshold T s l o p e calculated according to Equation (19). Therefore, for each observation epoch, the slope of a visible satellite can be considered as “large slope” if it was larger than T s l o p e and considered as “small slope” if it was smaller than T s l o p e .

4.3. Comparison of the Ideal Slope Threshold and the Ideal Protection Level

The functions of T s l o p e and V P L d are identical, forming a tangent P ( H M I 1 F ) curve to the lg P r ( H M I 1 F ) line. Referring to the equation for T s l o p e , i.e., Equation (18), the equation for V P L d , can be formulated as follows:
P a ( H M I 1 F = P ( V P E V P L d 1 F P ( T s < T d 1 F = [ 1 V P L d μ s a υ V P L d μ s a υ g ( t ; 0,1 ) d t ] 0 T d f ( y ; μ s S l o p e m a x 2 , K 4 ) d y )
Analyzing this equation, V P L d is determined by a v and S l o p e m a x with the given P r ( H M I 1 F ) , T d and K . Therefore, T s l o p e is only related to a v while V P L d is related to both a v and S l o p e m a x .
Considering that the PDF of the noncentral χ 2 distribution is too complicated, both T S l o p e and V P L d should be solved numerically, which will sharply increase the computation burden of a GNSS receiver or an onboard computer. Thanks to the one-to-one correspondence between a v and T S l o p e , T S l o p e can be calculated offline for discrete a v values in the range from 1 to T a v and saved in a receiver. It can be on-board searched from the presaved data according to the specific a v value. However, V P L d needs to be calculated online after both a v and S l o p e m a x obtained.

4.4. Slope-Based RAIM Availability Assessment Method

Slope-based RAIM availability assessment should be implemented in the worst case to fully prevent the integrity risk. Using slope-based method, RAIM is considered available if the maximum slope is less than the ideal threshold T S l o p e ; otherwise, it is considered unavailable. The specific execution process for slope-based RAIM availability assessment is presented in Figure 7.
As shown in Figure 7, the inputs are the observation matrix H and the weighted matrix P. The first step is an assessment based on the total number of visible satellites K. RAIM is considered unavailable if K 4 because fault detection cannot be executed with less than 5 visible satellites. For K > 4 , the second step is an assessment based on the a v value. RAIM is considered to be unavailable if a v T a v because only the sum of P ( H M I , 0 F ) and P ( H M I , 2 F ) has exceed P r ( H M I ) when a v exceeds T a v . For a v < T a v , the last step is an assessment based on the maximum slope S l o p e m a x . RAIM is finally considered available if S l o p e S l o p e m a x .
T a v can be calculated online according to Equation (6) with given V A L , P s a t , α , P r ( H M I ) and P r ( H M I 2 F ) . T S l o p e should be calculated offline with discrete a v with a small step size for different numbers of visible satellites and saved in the GNSS receiver or the onboard computer. Assuming that the total number of visible satellites is K and a v is in the interval of ( a v 1 , a v 2 ] , where a v 1 and a v 2 are the indexes of two adjacent discrete points [ a v 1 , T S l o p e 1 ] and [ a v 2 , T S l o p e 2 ] presaved for K visible satellites, T S l o p e corresponding to a v should be assigned as:
T S l o p e = T S l o p e 1 T S l o p e 2 a v 1 a v 2 ( a v a v 1 ) + T S l o p e 1
This value is calculated according to a linear fit for T S l o p e in the interval of ( a v 1 , a v 2 ] . Because the a v T S l o p e curve is convex as shown in Figure 5, the assigned T S l o p e value is smaller than the real T S l o p e value, which may cause a small conservative assessment risk. However, if the a v step size is small enough, the assigned T S l o p e value would be nearly equal to the true value, thereby preventing the small conservative assessment risk.
A simulation is designed to find the desirable range a v step. 1.2 × 1 0 6 times of a v and S lop e m a x for different location and epochs are collected. As shown in Table 2, the times of conservative assessment increases with step size widen. Taking into account both the amount of calculation and conservative assessment rate, 0.01~0.02 is the desirable range a v step because it is the maximum step size with 0 time of conservative assessment.
It should be mentioned that the calculation amount of the ideal-PL and T S l o p e is exactly the same. The ideal PL procedure begins with an improbably large VPL of 2000 m and halves the search step by checking if the corresponding integrity risk exceeds the requirement [3]. For each step of the iteration, the integrity risk needs to be calculated for different bias values with a fixed step. This process is computationally intensive. Similarly, the T S l o p e procedure begins with an improbably large slope value of 15 and halves the search step. Compared with the ideal-PL-based method, the slope-based method separates the process of numerical iteration from on-board RAIM availability assessment, reducing the burden of on-board computing.

5. Simulation

To compare the performance of the PL-based and slope-based methods, the vertical RAIM availability assessment for a 32-satellite GPS constellation is simulated in worldwide (latitude 60° S~60° N and longitude 180° W~180° E) for a whole day (13 March 2019 0:00:00~24:00:00). The simulation area is meshed as the grid of 1 ° × 1 ° and the simulation time step is 300 s. The masking angle is set to 10°. The dual-frequency ionosphere-free combination of L1 and L5 is assumed to be the pseudorange measurement. The standard derivation of the measurement noise is set according to the ARAIM interim report [11].
The vertical RAIM availability assessment is executed for all grid points, i.e., 121 × 361 = 43,681 grid points, using the classic-PL-based ( V P L c ), the enhanced-PL-based ( V P L e ), the ideal-PL-based ( V P L d ) and the slope-based methods respectively at each simulation epoch. In the simulation, the parameters related to vertical RAIM availability are P r ( H M I ) = 1 × 1 0 7 , V A L = 50 m [21], P r ( H M I , 2 ) = 1.3 × 1 0 8 [22], P s a t = 1 × 1 0 5 [11], α = 1 × 1 0 6 , and P M D = 0.1 % . The T S l o p e values are pre-calculated before the simulation with an a v step size of 0.01.

5.1. Specific Example Analysis

Three specific examples for a single grid point and a simulation epoch are chosen to show the RAIM availability assessment using the four methods in detail. Table 3 records the VPLs, S l o p e m a x and T S l o p e values for these specific examples. Figure 8, Figure 9, Figure 10 and Figure 11 present the real μ s lg P ( H M I 1 F ) curve for S l o p e m a x , the μ s lg P ( H M I 1 F ) curve for T S l o p e , and the μ s lg P a ( H M I 1 F ) curves for V P L c , V P L e and V P L d in these examples, respectively. In Figure 8, Figure 9, Figure 10 and Figure 11, the μ s lg P ( H M I 1 F ) curve for T S l o p e and the μ s lg P a ( H M I 1 F ) curve for V P L d are tangent to the lg P r ( H M I 1 F ) line, while the ξ b lg a P ( H M I 1 F ) curve for V P L e is separated from the lg P r ( H M I 1 F ) line. The positional relationship between the μ s lg a P ( H M I 1 F ) curve for V P L c and the lg P r ( H M I 1 F ) line could be either intersection or separation.
As shown in Figure 8, the real μ s lg P ( H M I 1 F ) curve for S l o p e m a x intersects with the lg P r ( H M I 1 F ) line, which means that RAIM is actually unavailable in the first example. The μ s lg P a ( H M I 1 F ) curve for V P L c is higher than the real μ s lg P ( H M I 1 F ) curve, i.e., V P L c < V A L , which means that RAIM is assessed to be available when using V P L c . All of the μ s lg P ( H M I 1 F ) curve for T S l o p e , the μ s lg P a ( H M I 1 F ) curves for V P L e and V P L d are lower than the real μ s lg P ( H M I 1 F ) curve, i.e., S l o p e S l o p e m a x , V P L e V A L and V P L d V A L , which means that RAIM is assessed to be unavailable when using S l o p e m a x , V P L e or V P L d . Therefore, the RAIM availability is optimistically assessed when using V P L c and successfully assessed using when using S l o p e m a x , V P L e or V P L d .
As shown in Figure 9, the real μ s lg P ( H M I 1 F ) curve for S l o p e m a x is separated from the lg P r ( H M I 1 F ) line, which means that RAIM is actually available in the second example. The μ s lg P ( H M I 1 F ) curve for T S l o p e , and the μ s lg P a ( H M I 1 F ) curves for V P L c and V P L d are higher than the real μ s lg P ( H M I 1 F ) curve, i.e., S l o p e < S l o p e m a x , V P L c < V A L and V P L d < V A L , which means that RAIM is assessed to be available when using S l o p e m a x , V P L c or V P L d . The μ s lg a P ( H M I 1 F ) curve for V P L e is lower than the real μ s lg P ( H M I 1 F ) curve, i.e., V P L e V A L , which means that RAIM is assessed to be unavailable when using V P L e . Therefore, RAIM availability is conservatively assessed when using V P L e and successfully assessed when using S l o p e m a x , V P L c or V P L d .
As shown in Figure 10, the real μ s lg P ( H M I 1 F ) curve for S l o p e m a x is separated from the lg P r ( H M I 1 F ) line, which means that RAIM is actually available in the third example. Both the μ s lg P ( H M I 1 F ) curve for T S l o p e and the μ s lg a P ( H M I 1 F ) curve for V P L d are higher than the real μ s lg P ( H M I 1 F ) curve, i.e., S l o p e < S l o p e m a x and V P L d < V A L , which means that RAIM is assessed to be available when using S l o p e m a x or V P L d . Additionally, both the μ s lg a P ( H M I 1 F ) curves for V P L c and V P L e are lower than the real μ s lg P ( H M I 1 F ) curve, i.e., V P L c V A L and V P L c V A L , which means that RAIM is assessed to be unavailable when using V P L c or V P L e . Therefore, RAIM availability is conservatively assessed when using V P L c or V P L e , and it is successfully assessed when using S l o p e m a x or V P L d in this example.
These three specific examples intuitively illustrate that both optimistic and conservative RAIM availability assessments might happen when using the classic-PL-based method; only conservative assessment might happen when using the enhanced-PL method; and both optimistic and conservative assessments can be prevented when using the ideal-PL-based or the slope-base methods.

5.2. Simulation Results Statistical Analysis

Both the rates of optimistic and conservative assessments for RAIM availability are calculated for each grid point using the classic-PL-based, enhanced-PL-based, ideal-PL-based and slope-based methods. The optimistic assessment rate is indicated by the ratio between the count of epochs in which RAIM is assessed to be available but is actually unavailable and the count of epochs in which RAIM is actually unavailable. The conservative assessment rate is indicated by the ratio between the count of epochs in which RAIM is assessed to be unavailable but is actually available and the count of epochs in which RAIM is actually available.
Table 4 records the specific optimistic and conservative assessment data for all 4 kinds of RAIM availability assessment methods. As shown in Table 4, both the rates of optimistic and conservative assessments are 0 for each grid point using the slope-based and ideal-PL-based methods. These prove that both optimistic and conservative assessments can be prevented. The reason for the 0 conservative assessment rate of the slope-based method is that the a v step size is small enough (0.01) for offline calculation of T S l o p e , leading to the assigned T S l o p e value being nearly equal to its true value.
For the classic-PL-based method, there are both grid points with nonzero optimistic assessment rate and grid points with nonzero conservative assessment rate. Figure 11 and Figure 12 present the optimistic and conservative assessment rates for each grid point using the classic-PL-based method, respectively. As shown in Table 4 and Figure 11, there are 2032 grid points with nonzero optimistic assessment rate, 4.65% of the total grid points, and there are 41 grid points with a 100% optimistic assessment rate, meaning that optimistic assessment always happens when RAIM is unavailable at these grid points. As shown in Table 4 and Figure 12, there are 16,678 grid points with nonzero conservative assessment rate, 38.18% of the total grid points, and the maximum conservative assessment rate is 2.59%. Comparing Figure 12 with Figure 11, the coverage area of conservative assessment is much larger than that of optimistic assessment, but the maximum optimistic assessment rate is much higher than the maximum conservative assessment rate for a single grid point using the classic-PL-based method.
For the enhanced-PL-based method, there are no grid points with a nonzero optimistic assessment rate. This finding illustrates that the value of P a ( H M I 1 F ) is small enough at ξ b = σ m 2 λ a s m m for the enhanced PL, which leads to the ξ b lg P a ( H M I 1 F ) curve always being separated from the P r ( H M I 1 F ) line, preventing optimistic assessment. However, there are 43,038 grid points with nonzero conservative assessment rate, as shown in Table 4 and Figure 13, representing 98.53% of the total grid points. The maximum conservative assessment rate is 8.15%. Comparing Figure 13 with Figure 12, the conservative assessment risk using the enhanced PL is much higher than that using the classic PL, which is represented by the much larger conservative assessment coverage area and the higher conservative assessment rate for a single grid point.
According to the above simulation results, the performance of the classic PL is the worst on RAIM availability assessment for both the risks of optimistic and conservative assessments. In particular, an optimistic assessment might cause HMI, which is intolerable. The performance of the enhanced PL is better than that of classic PL because the optimistic assessment is prevented. However, the risk of conservative assessment is significant, reducing RAIM continuity. The performances of the slope and the ideal PL are optimal, preventing both optimistic assessment and conservative assessment.

6. Conclusions

According to the theory and simulation analysis, the RAIM availability might be optimistically or conservatively assessed using the classic-PL-based method and might only be conservatively assessed using the enhanced-PL-based method. Using the ideal-PL-based method, both optimistic and conservative assessment can be prevented. However, the calculation of the ideal PL brings a heavy computational burden to the GNSS receiver or the onboard computer. The slope-based method has the same RAIM availability assessment result as the ideal-PL-based method. Because the ideal slope threshold is only related to one geometric parameter, it can be calculated offline and searched online. Thus, the on-board calculation burden can be reduced using the slope-based method. This improved method can be used in RAIM for single GNSS constellation. Further, a semi-physical simulation experiment will be implemented using the common on-board processor to verify the real-time performance of the slope-based method. Moreover, the RAIM availability assessment performance will be verified using massive actual measurement data.

Author Contributions

Conceptualization, D.S.; Methodology, J.Z.; Data curation, J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded jointly by the National Key R&D Program, China (No. 2022YFB3904400).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A. Distribution of Position Error

According to the pseudorange position principle, the linearized observation equation between the receiver and the visible satellites is as follows:
Δ Z = H Δ X + ε
In (A1), Δ Z is the K-dimensional pseudorange offset vector from the real receiver position to the linearization point, where K is the total number of visible satellites. Δ X = [ Δ x , Δ y , Δ z , Δ t 1 , , Δ t n ] T , where [ Δ x , Δ y , Δ z ] T and [ Δ t 1 , , Δ t n ] T are respectively the offset vectors for the position and clock errors from the real receiver position to the linearization point; Δ t n is the clock error for the n-th GNSS. H is a K × ( 3 + n ) linearized observation matrix in the local Cartesian coordinate (ENU). ε is the observation error vector. The observation error includes both noise and bias. The measurement noise for the i-th visible satellite, denoted as V S i , is assumed to be the independent WGN with mean zero and variance σ i 2 ( i = 1 , , K ). σ i is set according to the ARAIM interim report [21]. It should be mentioned that the measurement noise is modeled as non-zero mean Gaussian distribution with a nominal bias in ARAIM interim report [21]. This nominal bias is not considered here to simplify the formula derivation. It really needs to be considered in practical application.
σ i 2 = σ U R E , i 2 + σ t r o p o , i 2 + σ u s e r , i 2
where σ U R E , i is the standard deviation of a Gaussian error distribution modelling the range component of the signal-in-space error, containing the ephemeris and satellite clock errors. σ t r o p o , i is the standard deviation of the tropospheric error. σ u s e r , i is the standard deviation of the user error, including the multipath error and the receiver noise. σ U R E , i , σ t r o p o , i and σ u s e r , i are calculated using the equations from the ARAIM interim report [21].
The weighted least squares solution for X in (A1) is as follows:
X ^ = ( H T PH ) 1 H T PZ
In (A3), P = d i a g [ ( σ i 2 ) 1 ] is the weighted matrix, where d i a g [ ( σ i 2 ) 1 ] presents a diagonal matrix with the diagonal element ( σ i 2 ) 1 in the i-th row.
According to (A1) and (A3), the solution error for X , represented by μ , can be expressed as follows:
μ = X ^ X = ( H T PH ) 1 H T P ε
Set A = ( H T PH ) 1 H T , μ = AP ε . Thus, the VPE can be expressed as follows:
μ 3 = ( a 31 ε 1 σ 1 2 + a 32 ε 1 σ 2 2 + + a 3 K ε K σ K 2 )
where a 3 i is the element in the 3-th row of A .
Under the nominal mode:
μ 3 N ( 0 , a v 2 )
where N ( μ , σ 2 ) represents a normal distribution with mean μ and variance σ 2 ;   a v = i = 1 K a 3 i 2 σ i 2 .
Under the single-satellite fault mode:
μ 3 N ( a 3 m ξ b σ m 2 , a v 2 )
where V S m presents the faulty satellite and the measurement bias for V S m is ξ b .
Therefore, the PDF of μ 3 under the single-satellite fault mode is g ( x ; a 3 m ξ b σ m 2 , a v 2 ) , where g ( x , m , σ 2 ) represents the PDF of a normal distribution with the mean value m and the variance σ 2 . Therefore:
P ( V P E V A L 1 F ) = 1 V A L V A L g ( x , a 3 m ξ b σ m 2 b , a v 2 ) d x

Appendix B. Distributions of Test Statistic for LSR Algorithm

The pseudorange residual vector, recorded as ω , is defined by combining (A1) and (A3) as follows:
ω = Z Z ^ = Z H X ^ = [ I H ( H T P H ) 1 H T P ] ε
Set S = P 1 H ( H T PH ) 1 H T ,
ω = SP ε = S 1 ε 1 σ 1 2 + S 2 ε 2 σ 2 2 + + S K ε K σ K 2
where S i is the i-th column of S .
The weighted sum of squares for ω , recorded as S S E , can be calculated with (A9)
S S E = ω T P ω
The test statistic T s for faulty satellite detection in the LSR algorithm is S S E :
T s = ω T P ω
Under the nominal mode:
T s χ 2 ( K 4 )
where χ 2 ( v ) represents a standard χ 2 distribution with the freedom degree of v .
Under the single-satellite fault mode:
T s χ 2 ( K 4 , λ )
where χ 2 ( v , λ ) represents non-central χ 2 distribution with the freedom degree of v and the decentralized parameter of λ . λ can be considered as the difference between the mean values of test statistics for nominal and single-satellite fault modes, calculated as follows.
With the assumption that the single faulty satellite is V S m and the measurement bias for V S m is ξ b , ω can be rewritten as follows:
ω = ω n o + S m ξ b σ m 2
In (A15), ω n o represents the pseudorange residual vector with no faulty satellite:
ω n o = S 1 ε 1 σ 1 2 + + S m 1 ε m 1 σ m 1 2 + S m ξ m σ m 2 + S m + 1 ε m + 1 σ m + 1 2 + + S K ε K σ K 2
where ξ m and ξ b are respectively the measurement noise and bias for V S m , i.e., ε m = ξ m + ξ b . Therefore, the mean value of T s under the single satellite faulty mode is:
E [ T s ] = E [ ω T P ω ] = E [ ( ω n o + S m ξ b σ m 2 ) T P ( ω n o + S m ξ b σ m 2 ) ] = E [ ω n o T P ω n o + ω n o T P ( S m ξ b σ m 2 ) + ( S m ξ b σ m 2 ) T P ω n o + ( S m ξ b σ m 2 ) T P ( S m ξ b σ m 2 ξ b ) ]
In (A17), E [ ω n o T P ω n o ] is the mean value of T s under the nominal mode and E [ ω n o T P ( S m ξ b σ m 2 ) ] = E [ ( S m ξ b σ m 2 ) T P ω n o ] = 0 . Therefore:
λ = E [ ξ b σ m 2 2 S m T P S m ]
S m T P S m is the m-th diagonal element for S T PS . Because S T PS = S :
S m T P S m = s m m
where s m m is the diagonal element in the m-th row of S . Substituting (A19) into (A18):
λ = ξ b σ m 2 2 s m m
Therefore, the PDF of T s under the single-satellite fault mode is f ( y ; ξ b σ m 2 2 s m m , K 4 ) , where f ( y ; λ , K ) represents the PDF of a noncentral χ 2 distribution with the noncentral parameter λ and the degrees of freedom K . P ( T s < T d 1 F ) can be calculated as follows:
P ( T s < T d 1 F ) = 0 T d f ( y , K 4 , λ ) d y

References

  1. Feng, S.; Ochieng, W.Y.; Walsh, D.; Ioannides, R. A Measurement Domain Receiver Autonomous Integrity Monitoring Algorithm. GPS Solut. 2006, 10, 85–96. [Google Scholar] [CrossRef]
  2. Brown, R.G. A Baseline GPS RAIM Scheme and a Note on the Equivalence of Three RAIM Methods. Navigation 1992, 39, 301–316. [Google Scholar] [CrossRef]
  3. Milner, C.D.; Ochieng, W.Y. Weighted RAIM for APV: The Ideal Protection Level. J. Navig. 2011, 64, 61–73. [Google Scholar] [CrossRef]
  4. Madonna, P.; Viola, S.; Sfarzo, L. NIORAIM Algorithm Applied to a Multiconstellation GNSS: Analysis of Integrity Monitoring Performances in Various Phases of Flight. In Proceedings of the IEEE/ION Position, Location and Navigation Symposium, Indian Wells, CA, USA, 4–6 May 2010; IEEE: Piscataway, NJ, USA, 2010; pp. 1258–1263. [Google Scholar]
  5. Lee, Y.C. Two New RAIM Methods Based on the Optimally Weighted Average Solution (OWAS) Concept. Navigation 2007, 54, 333–345. [Google Scholar] [CrossRef]
  6. Rife, J.; Pervan, B. Overbounding Revisited: Discrete Error-Distribution Modeling for Safety-Critical GPS Navigation. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 1537–1551. [Google Scholar] [CrossRef]
  7. Fan, G.; Xu, C.; Zhao, J.; Zheng, X. A Distribution Model of the GNSS Code Noise and Multipath Error Considering Both Elevation Angle and Orbit Type. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2019, 233, 1900–1915. [Google Scholar] [CrossRef]
  8. Brown, R.G.; Chin, G.Y. GPS RAIM: Calculation of the Threshold and Protection Radius Using Chi-Square Methods—A Geometric Approach; Radio Technical Commission for Aeronautics: Washington, DC, USA, 1994. [Google Scholar]
  9. Kaplan, E.D.; Hegarty, C. (Eds.) Understanding GPS/GNSS: Principles and Applications; Artech House: New York, NY, USA, 2017. [Google Scholar]
  10. Angus, J.E. RAIM with Multiple Faults. Navigation 2006, 53, 249–257. [Google Scholar] [CrossRef]
  11. GPS-Galileo Working Group C ARAIM Technical Subgroup Interim Report. Available online: https://www.gps.gov/policy/cooperation/europe/2013/working-group-c/ (accessed on 14 April 2024).
  12. Blanch, J.; Walker, T.; Enge, P.; Lee, Y.; Pervan, B.; Rippl, M.; Spletter, A.; Kropp, V. Baseline Advanced RAIM User Algorithm and Possible Improvements. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 713–732. [Google Scholar] [CrossRef]
  13. Blanch, J.; Walter, T.; Enge, P. Theoretical Results on the Optimal Detection Statistics for Autonomous Integrity Monitoring. Navig. J. Inst. Navig. 2017, 64, 123–137. [Google Scholar] [CrossRef]
  14. Joerger, M.; Langel, S.; Pervan, B. Integrity Risk Minimisation in RAIM Part 2: Optimal Estimator Design. J. Navig. 2016, 69, 709–728. [Google Scholar] [CrossRef]
  15. Meng, Q.; Liu, J.; Zeng, Q.; Feng, S.; Xu, R. Improved ARAIM Fault Modes Determination Scheme Based on Feedback Structure with Probability Accumulation. GPS Solut. 2019, 23, 16. [Google Scholar] [CrossRef]
  16. Zhai, Y.; Joerger, M.; Pervan, B. Fault Exclusion in Multi-Constellation Global Navigation Satellite Systems. J. Navig. 2018, 71, 1281–1298. [Google Scholar] [CrossRef]
  17. Jiang, Y.; Wang, J. A New Approach to Calculate the Vertical Protection Level in A-RAIM. J. Navig. 2014, 67, 711–725. [Google Scholar] [CrossRef]
  18. Jiang, Y.; Wang, J. A New Approach to Calculate the Horizontal Protection Level. J. Navig. 2016, 69, 57–74. [Google Scholar] [CrossRef]
  19. Walter, T.; Enge, P.; Blanch, J.; Pervan, B. Worldwide Vertical Guidance of Aircraft Based on Modernized GPS and New Integrity Augmentations. Proc. IEEE 2008, 96, 1918–1935. [Google Scholar] [CrossRef]
  20. RTCA DO-245A; Minimum Aviation System Performance Standard for the Local Area Augmentation System (LAAS). RTCA: Washington, DC, USA, 2004.
  21. Bhattacharyya, S.; Gebre-Egziabher, D. Kalman Filter–Based RAIM for GNSS Receivers. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 2444–2459. [Google Scholar] [CrossRef]
  22. Panel, G. GNSS Evolutionary Architecture Study (GEAS) Phase II Panel Report; FAA: Washington, DC, USA, 2010; pp. 10–12. [Google Scholar]
  23. Dan, S.; Chuang, S.; Zhipeng, W.; Cheng, W.; Guifei, J. Correlation-Weighted Least Squares Residual Algorithm for RAIM. Chin. J. Aeronaut. 2020, 33, 1505–1516. [Google Scholar]
Figure 1. P ( H M I 1 F ) for different satellites.
Figure 1. P ( H M I 1 F ) for different satellites.
Sensors 24 03283 g001
Figure 2. Separation.
Figure 2. Separation.
Sensors 24 03283 g002
Figure 3. Intersection.
Figure 3. Intersection.
Sensors 24 03283 g003
Figure 4. Tangency.
Figure 4. Tangency.
Sensors 24 03283 g004
Figure 5. T S l o p e for different a v values.
Figure 5. T S l o p e for different a v values.
Sensors 24 03283 g005
Figure 6. Practical meaning of T S l o p e .
Figure 6. Practical meaning of T S l o p e .
Sensors 24 03283 g006
Figure 7. Slope-based RAIM availability assessment process.
Figure 7. Slope-based RAIM availability assessment process.
Sensors 24 03283 g007
Figure 8. RAIM availability judgment for the first specific example.
Figure 8. RAIM availability judgment for the first specific example.
Sensors 24 03283 g008
Figure 9. RAIM availability assessment for the second specific example.
Figure 9. RAIM availability assessment for the second specific example.
Sensors 24 03283 g009
Figure 10. RAIM availability assessment for the third specific example.
Figure 10. RAIM availability assessment for the third specific example.
Sensors 24 03283 g010
Figure 11. Optimistic assessment rate using the classic-PL-based method.
Figure 11. Optimistic assessment rate using the classic-PL-based method.
Sensors 24 03283 g011
Figure 12. Conservative assessment rate using the classic-PL-based method.
Figure 12. Conservative assessment rate using the classic-PL-based method.
Sensors 24 03283 g012
Figure 13. Conservative assessment rate using the enhanced-PL-based method.
Figure 13. Conservative assessment rate using the enhanced-PL-based method.
Sensors 24 03283 g013
Table 1. Vertical slope values for different visible satellites.
Table 1. Vertical slope values for different visible satellites.
PRN256121315192529
S l o p e 0.1626.0603.0101.5762.4251.6172.0360.9251.636
Table 2. Step and times of conservative assessment.
Table 2. Step and times of conservative assessment.
Step0.0010.010.020.030.040.050.10.20.5
Times of conservative assessments000222224
Table 3. Specific examples for RAIM availability assessment.
Table 3. Specific examples for RAIM availability assessment.
NumLocationUTC TimeKavSlopemaxVPLcVPLeVPLdTSlope
18° N, 82° W03:55:0086.3465.23245.28766.16850.4245.172
245° S, 132° W00:00:0082.7045.33046.13955.03745.4755.887
353° N, 122° E08:35:00101.3585.66650.65655.30549.6835.702
Table 4. Performance comparison for the four kinds of RAIM availability assessment methods.
Table 4. Performance comparison for the four kinds of RAIM availability assessment methods.
ItemGrid Points of Optimistic AssessmentGrid Points of Conservative AssessmentMaximum Optimistic Assessment Rate on Single Grid PointMaximum Conservative Assessment Rate on Single Grid Point
Method
classic-PL-based203216,678100%2.59%
enhanced-PL-based043,03808.15%
Ideal-PL-based0000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, J.; Song, D.; Wang, J. An Improved RAIM Availability Assessment Method Based on the Characteristic Slope. Sensors 2024, 24, 3283. https://doi.org/10.3390/s24113283

AMA Style

Zhao J, Song D, Wang J. An Improved RAIM Availability Assessment Method Based on the Characteristic Slope. Sensors. 2024; 24(11):3283. https://doi.org/10.3390/s24113283

Chicago/Turabian Style

Zhao, Jing, Dan Song, and Jitao Wang. 2024. "An Improved RAIM Availability Assessment Method Based on the Characteristic Slope" Sensors 24, no. 11: 3283. https://doi.org/10.3390/s24113283

APA Style

Zhao, J., Song, D., & Wang, J. (2024). An Improved RAIM Availability Assessment Method Based on the Characteristic Slope. Sensors, 24(11), 3283. https://doi.org/10.3390/s24113283

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop