A Holistic Approach to Guarantee the Reliability of Positioning Based on Carrier Phase for Indoor Pseudolite

be guaranteed. Abstract: The integrity monitoring algorithm based on pseudorange observations has been widely used outdoors and plays an important role in ensuring the reliability of positioning. However, pseudorange observations are greatly affected by the error sources such as multipath, clock drift, and noise in indoor pseudolite system, thus the pseudorange observations cannot be applied to high-precision indoor positioning. In general, double differenced (DD) carrier phase observations are used to obtain a high-precision indoor positioning result. What’s more, the carrier phase-based integrity monitoring (CRAIM) algorithm is applied to identify and exclude potential faults of the pseudolites. In this article, a holistic method is proposed to ensure the accuracy and reliability of positioning results. Firstly, if the reference pseudolite operates normally, extended Kalman filter is used for parameter estimation on the premise that the number of common pseudolites meets positioning requirements. Secondly, the innovation sequence in the Kalman filter is applied to construct test statistics and the corresponding threshold is determined from the Chi distribution with a given probability of false alert. The pseudolitehorizontal protection level (HPL) is calculated by the threshold and a prior probability of missed detection. Finally, compared the test statistics with the threshold to exclude the faultypseudolite for the reliability of positioning. The experiment results show that the proposed method improves the accuracy and stability of the results through faults detection and exclusion. This method ensures accuracies at the centimeter level for dynamic experiments and millimeter levels for static ones.


Introduction
The global navigation satellite system (GNSS) has been widely applied in outdoor positioning. Nevertheless, it is difficult to obtain reliable GNSS signals in sheltered environments or indoors [1], where a pseudolite system can make up for the shortage of GNSS [2,3].
Relative to the GNSS, indoor pseudolite observations have no tropospheric and ionospheric errors. The PL/receiver clock and antenna phase error can also be eliminated by DD observations which greatly reduce the difficulty in processing error sources in indoor PL positioning [4]. However, the multipath effect, noise, and fixed satellite constellations are still the problems indoors. What is worse, the indoor pseudolite system has a great nonlinear effect due to the close distance between the receivers and pseudolites. If initial position estimation has poor accuracy, the nonlinear error caused by linearization will be severe. Therefore, large code measurement errors caused bythe above factors may result in low accuracy and reliability of positioning results.
At present, the research on the indoor pseudolite system mainly focuses on the design of pseudolite constellation or high-precision positioning based on the carrier phase, but few on the reliability of positioning. According to quality control methods of GNSS carrier phase positioning, they can be divided into two categories: priori and posteriori quality control methods. As to priori methods, there are some classical methods, such as cycle slip detection, DD method. Cycle slip detection can effectively eliminate such faults as cycle slip or signal interruption. The double differenced method will eliminate clock difference and weaken troposphere and ionosphere delays. In general, the posteriori quality control method does not distinguish which kind of error sources, such as RAIM, IGG III, and variance component estimation. These methods can effectively improve the positioning accuracy. However, the existing legacy approaches cannot eliminate all faults, such as multipath effect, noise and so on, which is also the current problem. In an indoor environment, the influence of multipath and noise is greater than outdoors. A holistic approach to guarantee the reliability of positioning for indoor pseudolite based on carrier phase is proposed and the feasibility of the method is verified by some experiments.The proposed method creatively considers the fault of indoor positioning source (pseudolite) and other error sources, including complete fault detection, exclusion procedure, andformula derivation of pseudolite protection level. Faults on reference pseudolite are also considered. Although we cannot deal with any specific error source, the data with poor quality can be effectively identified no matter which error source. By excluding the corresponding epochs, the positioning accuracy and reliability will be guaranteed.
The research on indoor pseudolite high-precision positioning algorithms has become a hotspot in recent years and the key is to use carrier phase observations and ambiguity resolution [4][5][6].Reference [4] has successfully introduced AFM (Ambiguity Function Method) into PL-AR (Pseudolite-Ambiguity Resolution). However, due to the multi-peak characteristic of AFM, the reliability of the solutions is sometimes difficult to be guaranteed compared with the LAMBDA method. In addition, the way to ensure the reliability of positioning has not been discussed. References [5,6] proposed a new OTF-AR method for ground-based positioning systems which solely relies on carrier phase measurements and applies to cases where no sufficiently accurate measurements were provided in advance; however, the gross error must be given careful consideration in practical applications and the experiments in the references were conducted under relatively good conditions. Actually, it is the poor conditions that lead to unreliable positioning results. Some research institutes used KPI when LAMBDA and ratio tests are applied to obtain the integer ambiguity resolution [7]. The most mature indoor pseudolite system is produced by LOCATA [8]. ProprietaryTimeLoc technology is utilized for tight time synchronization, allowing the clock difference in the carrier phase measurement equation to be ignored. The horizontal error of the nondifferential positioning is less than 6 cm, and the precision of the altitude direction is about 15 cm [9]. The frequency stability and accuracy of the PL clock are much inferior to the GPS, and PL clock errors are difficult to obtain accurately [2]. Although LOCATA can ensure PL clock synchronization, many pseudolite systems lack this costly technology at present. Reference [10] achieved centimeter-level PPP (Precise point positioning) results with a new prototype pseudolite system without relying on accurate time synchronization. It is believed that PL would be more extensively applicable for integration with GNSS and INS (Inertial Navigation System). Reference [11] proposed a new concept of pseudolite-based navigation to overcome the receiver modification problem of systems. However, only meter-level positioning precision could be achieved in poor geometry tests.References [10,11] did not consider quality control under severe conditions as well. In reference [12], Xu et al. built a new pseudolite-based indoor positioning system and the clock could be synchronized with each other. Nevertheless, only simulated data wereused to obtain high-precision positioning results. References [13,14] believed that the convergence precision of initial coordinate values was important for positioning, but how to ensure centimeter-level positioning precision has not been considered with bad initial coordinate value. Reference [15] used carrier interferometry to conduct analyses on indoor pseudolite positioning. However, this method has great limitations in applications. Reference [16] proposed a new method for indoor pseudolite positioning based on AFM and DD pseudorange observations with good data. However, it is generally known that data quality has a great influence on the AFM.
Current high-precision PL positioning algorithms mainly are applied based on carrier phase observations. However, the indoor environment is more severe than outdoors. Whether the carrier phase observations or the positioning results are reliable remains to be discussed. To the authors' knowledge, there are few methods to guarantee the reliability of the results for indoor pseudolite positioning and no effective strategy to ensure the positioning accuracy whena fault happens. Therefore, the application of receiver autonomous integrity monitoring (RAIM) for pseudolite system is of great significance in practical positioning.
There are two types of RAIM algorithm, snapshot RAIM, and sequential RAIM. The snapshot RAIM includes the parity vector method, the least-square residual method [17][18][19], weighted RAIM algorithm, etc. Sequential RAIM is almost based on the Kalman filter which is an indispensable component of many navigation systems [20]. In view of the above, the carrier phase-based RAIM (CRAIM) algorithm [21] proposed by Feng.was applied tothe high-precision positioning of the indoor pseudolite system. Reference [22] proposed an INS-assisted RAIM algorithm to increase the level of reliability and integrity for an integrated navigation system. Reference [23] explored the application of the snapshot RAIM method for GNSS vessel receivers when space segment geometry was poor. However, the proposed algorithm did not discuss how to determine the protection level with a given probability of missed detection (MD) and false alarm (FA). Reference [24] showed a correlation between traditional Shapiro-Wilk test statistics when multiple faults occur. The paper proposed a criterion to overcome the correlation based on the reliability principle and designed a multi-fault processing method based on an extended Shapiro-Wilk test. Li et al. proposed a complete algorithm to ensure reliability in attitude measurement [25]. In this article, a new integrity monitoring method was proposed for the GNSS-based attitude determination and they constructed the test statistic by using the measurement residuals to detect the potential failure. Furthermore, they derived the attitude protection level based on the worst-case protection concept under the multiple faults hypothesis. Reference [26] indicated that carrier phase observations error does not subject to the standard normal distribution and a carrier phase-based RAIM algorithm using a Gaussian sum filter was proposed. Since pseudorange-based RAIM (PRAIM) cannot meet the requirement of PPP, References [27][28][29] applied CRAIM to GNSS PPP and RTK algorithm for airport surface movement. Similar to GNSS PPP and RTK, carrier phase observations must also be used in high-precision indoor pseudolite positioning. In references [30,31], a specific algorithm for multiple constellations was studied and a method of integrity monitoring-based ambiguity validation for triple-carrier ambiguity resolution was proposed. They considered that ambiguity validation methods include the ratio test, the quantitative evaluation based on the success rate and the equivalent failure rate of ambiguity resolution. However, the ambiguity correctness validation cannot necessarily guarantee the reliability of indoor pseudolite positioning. Sometimes, an incorrectly fixed ambiguity will be obtained because of the severe environment and the poor satellite geometry. Therefore, it is necessary to carry out the integrity monitoring for indoor pseudolite rather than solely using the ambiguity validation method.
The RAIM algorithm is mainly used in safety-critical fields such as aircraft navigation at present and is almost unmanned in indoor pseudolite. In fact, indoor pseudolite positioningis seriously influenced by many environmental factors, such as people's walking, multipath effects from wallsand so on. Sometimes it is difficult to obtain a high-precision position result. In addition, even if a specific positioning result can be obtained, the reliability may not be guaranteed. This paper focuses on the RAIM algorithm of indoor pseudolite system based on EKF, including PRAIM and CRAIM algorithm. Some experiments are operated to verify the feasibility of the indoor pseudolite RAIM algorithm and the reliability of the positioning results, which can expand the availability of pseudolite positioning. According to the actual characteristics of pseudolite, an algorithm to guarantee the accuracy and reliability of indoor positioning for pseudolite based on CRAIM is proposed, which is named PLRAIM (Pseudolite based RAIM). It is especially important to use this method in harsh indoor environments.

Algorithm
The algorithm includes two aspects, the calculation of the protection level and fault detection. All algorithms in this paper are implemented in horizontal positioning.

Observation Model in Pseudolite System
Due to the complex indoor environment and the poor stability of commercial pseudolite clocks, the pseudolite high-precision positioning algorithms are generally realized based on DD carrierphase positioning mode. The error of clock synchronization and antenna phase can be eliminated or weakened by theDD model [32]. The DD observations do not include tropospheric and ionospheric delay because the signal from the pseudolite itself does not pass through the troposphere and ionosphere. Furthermore, since the pseudolite constellation is fixed, there are no ephemeris errors and relativistic effects. Therefore, the indoor single-frequency pseudolite DD observation equations can be formulated as [33][34][35] where ∇Δ is a DD operator; r and b represent the rover and base stations, respectively. k and s represent pseudolites. The pseudorange and carrier phase observations are represented by P and φ. ρ represents the geometric distance between pseudolites and the receivers. e and ε denote the noise of the pseudorange and carrier phase observations, respectively. M and m represent the multipath effects of pseudorange and carrier phase observations. N indicates the integer ambiguities of carrier phase observations. λ indicates the wavelength.

Parameter Estimation
Extended Kalman filter (EKF) is applied to estimate the parameters ina high-precision positioning solution. EKF is mainly to solve the linearization problem in KF. It can be considered that the coordinates remain the same because the rover station is a state of low dynamic or static in indoor positioning. The transition matrix is I and the covariance matrix of the system process noise is O. The measurement noise should be considered. The observation and the state model are as follows, where k is the current epoch. X represents the state vector. Ф is the transition matrix I. ω is the dynamic process noise of the system. Z is a measurement vector after linearization. H is a design matrix [Ir λ] and v indicates the measurement noise. If the dynamic process noise and the measurement noise of the system are Gaussian white noise, then the parameters have the following statistical characteristics, After the initial coordinate being obtained, the corresponding initial state parameter X0 and its variance-covariance matrix P0 can be given. The basic formula of EKF follows, where subscript k denotes the current epoch. X is state vector; P is covariance matrix. K is Kalman gain matrix and V is innovation vector. W is weight matrix. R is measurement noise matrix. (+) and (−) denote predicted and updated values. When positioning, the iterative EKF is usually used to weaken the linearization error. In order to avoid the problem of divergence caused by the imprecise of a given dynamic model during filtering, the weight can be adjusted according to the residual matrix to achieve the robust Kalman filter [36].

The Procedure of PLRAIM
The reliability of the positioning algorithm to obtain the rover coordinates is greatly affected by multipath and ambiguity resolution. Therefore, this paper focuses on studying the reliability of positioning results. Integrity monitoring mainly includes two aspects: fault detection and alert and availability of the RAIM algorithm. The main idea is to detect and excludefaultsbased on the consistency between redundant observations [37,38]. The construction of test statistics is one of the key factors, consequently. Considering the actual situation of our indoor pseudolite positioning system, the steps of PLRAIM are as follows: (a) Obtain data of rover and base station and ensure there are more than three common visible pseudolites in order for the horizontal positioning to be performed; (b) Judge fault on reference pseudolite (rpl). After the reference pseudolite is determined by the elevation angle, the priori residuals of DD pseudorange observations can be calculated. If the residuals of all DD observations perform largely, it can be considered that the reference pseudolite is faulty and it should be replaced by another one; (c) When the number of common pseudolites meets requirements and reference pseudolite is normal, parameter estimation and ambiguity resolution will be performed. The classical LAMBDA method is used for ambiguity resolution (AR) [39,40]; (d) Judge the availability of RAIM. Compare HPL calculated by the probability of missed detection (Pmd) with HAL. We set Pmd = 1 × 10 −3 [21]. As to HAL, it is only used to set a maximum error range for comparison with HPL. Referring to GNSS positioning, the accuracy of high-precision positioning needs at least a centimeter-level. Therefore, in high-precision indoor pseudolite positioning, the positioning accuracy withcentimeter-level would be reliable and available. The HAL of indoor pseudolite high-precision positioning is set to 0.1m; (e) Fault detection of pseudolite positioning system. Construct the test statistic T and compare it with threshold Th determined from the probability of false alert (Pfa). We set Pfa = 3.33 × 10 −9 [21]. If T is larger than the threshold, the pseudolite may be faulty; (f) Fault exclusion. The most intuitive method is to remove pseudolites in turn and compare the test statistic. Then exclude this pseudolite or delete corresponding epochs when a fault occurs; (g) Positioning. After the above steps, the reliability of results is effectively improved. The flow chart follows, In Figure 1, HPL is calculated before fault detection. If HPL is greater than HAL, the current observation condition (pseudolite constellation) would not be good enough and the unavailable situation of the system shown in Figure 2 may appear. In Figure 2b, HMI means hazardously misleading information and MI means misleading information. There are two important concepts in fault detection, false alert (FA) and missed detection (MD). Pfa refers to the probability of the fault detected by the system and alerts on the condition that the pseudolite is not faulty or HPE is less than HPL. FA occurs when the navigation is normal in Figure 2. Pmd is the probability that the system fails to detect it while in fact,the fault happens and HPE is greater than HPL. The MD shown as MI and HMI in Figure 2 mainly occurs when the navigation is not normal. HPL guarantees that Pmd will not exceed demand under the current pseudolite constellation and the threshold ensures that there is no FA in fault detection.
The determination of HPL and test statistic threshold will be presented in the next section in conjunction with the pseudolite system. In summary, under the premise that RAIM availability is guaranteed (HPL<HAL), the following four situations in Figure 3 may occur during fault detection.

Fault Detection
In fault detection, the most critical step is the construction of test statistic T and the calculation of threshold Th. The test statistic is a value calculated from a set of sample data used to determine whether the null hypothesis should be rejected in the hypothesis test. The construction of T depends on the probability model (distribution) and the assumptions of the problem. As long as T is constructed, Th can be determined from a probability model with a given confidence level.
According to (2), the residual vector of the DD observations is obtained as, The simultaneous (6) and (7) can be written as, The residual squared SSE is a scalar, The posterior unit weight mean error and test statistic can be constructed as, where t represents the necessary number of observations. In experiments, parameters to be estimated are the horizontal coordinates x, y, and ambiguity, thereby t = 3. The innovation sequence V provides the corresponding source of information for PLRAIM. Nevertheless, when DD observations are applied, the noise is no longer independent. The dependence is a mathematical correlation rather than a physical correlation. To de-correlate the dependence, T can be constructed as [21], where residual vectorr = [rL1rP1] T , thus, the statistic of each observation can be expressed as, In (12), TL1 is a test statistic based on carrier phase L1 observations and TP1 is constructed based on pseudorange observations P1. These two statistics form the core of the CRAIM and PRAIM algorithm respectively. Some faults can be detected when using PRAIM, for example, indoor multipath. However, there are also some faults that cannot be reflected from the test statistic based on pseudorange (e.g.cycle slips) and such faults cannot be ignored in high-precision indoor pseudolite positioning. Consequently, the proposed indoor PLRAIM can also achieve the goal of data pre-processing such as cycle slip detection and gross error detection. Moreover, since the results of indoor pseudolite positioning are greatly affected by data quality, the PLRAIM algorithm is also a guarantee for the availability and reliability of positioning.
The innovation sequence in a Kalman filter is assumed to subject to Gaussian distribution (Mehra 1970), N(0, σ 2 ). According to the statistical distribution theory, test statistic T 2 /σ 2 subjects the chi-square distribution with a degree of freedom of m. When there is no deviation in the observations, the system is under a normal state. If an alert is performed, it would be a false alert. The following probability equation can be established from a given Pfa, The threshold Th can be determined by (13) and T/σ is used to compare with Th. If T/σ is greater than Th, it can be considered that the fault exists in observations that should be excluded. If T/σ is less than the threshold, the positioning result is reliable and the observations can be applied for positioning.
In addition, the selection principle of reference pseudolite is also based on the elevation angle in indoor pseudolite positioning. However, elevation angles of all pseudolites are relatively close owing to the small scene indoors. Currently, it seems not reliable to select a reference pseudolite according to the highest elevation angle. Therefore, the possibility of reference pseudolite fault is also considered in PLRAIM. The fault of reference pseudolite is easier to detect. Once the fault in reference pseudolite occurs, the prior residuals of DD pseudorange observations of all pseudolites will be large and the reference pseudolite should be replaced in time.

HPL and the Availability of PLRAIM
This article focuses on horizontal positioning, so only the horizontal protection level (HPL) is analyzed. HPL refers to the horizontal positioning error that can be guaranteed when a fault occurs in current pseudolite constellation before fault detection. If HPL exceeds a given HAL, the current pseudolite constellation is not available. In other words, HPL is to ensure that the maximum Pmd under the current pseudolite constellation does not exceed the demand. There are usually two methods to calculate: The first way is to use the geometric factor SLOPE between horizontal position error (HPE) and test statistic [41]. SLOPE represents the mapping of the test statistic to HPE, reflecting the relationship between them. If there is a deviation in the test statistic, it will be reflected in the SLOPE.
The horizontal position errors are the first two elements of KkVk inan extended Kalman filter. Let S = (I -HkKk) T (I -HkKk), assuming that the deviation occurs only in the i th pseudolite and the other pseudolites are normal. HSLOPE(i) can be used to indicate the sensitivity of horizontal position error to the deviation of the i th observations, where Sii represents the i th element in diagonal matrix S. SLOPE is purely a geometric factor. Corresponding to the same test statistic, Pmd is the highest one when a fault occurs on apseudolite with the largest SLOPE. Therefore, the larger the HSLOPE, the more difficult it is to detect. In other words, when SLOPE is the largest, its statistic is the smallest in the case of the same HPE, which means that its residuals are the smallest, too. At this time, the minimum detectable deviation bias (MDB) is set to pbias and it can be calculated according to the threshold and a given Pmd.
If a fault occurs, T/σ should theoretically be greater than Th. If T/σ is smaller than Th, it indicates that a missed detection happens. The following probability equation can be established for a given Pmd, where λ is the noncentral parameter, pbias 2 = σ 2 λ/(n−t). σ represents the standard deviation of carrier phase measurements. The horizontal protection level HPL can be calculated as,

HPL HSLOPE pbias
The other way is to use the position estimate uncertainty, in which the first three elements of covariance matrix (P) is applied, The positioning error can be regarded as Gaussian noise and kH(Pr) can be obtained based on a given confidence probability Pr (usually 1 -Pmd). In (17),σH represents the horizontal position uncertainty, 11 22 = H P P σ + where P represents covariance matrix in expression (4). According to (16) and (17), HPL can be expressed as,

Experiment and Analyses
Five pseudolites have been set up in the laboratory. The ranging signals with C/A codes are transmitted on the GPS L1 frequency. The RF front end of the receiver is a USRP and a full series of GNSS signals can be captured with a suitable daughter board. Taking a certain mark as the origin point, the two horizontal directions of the floor are X and Y axis, respectively. The Z-axis is perpendicular to the floor and a right hand indoor independent coordinate system is established. All the experiments are carried out on the floor. In addition, the scope of the laboratory is about 10m×7m×4m. The pseudolites constellation is as shown in Figure 4a,b. The HDOP (Horizontal Dilution of Precision) on the floor under this condition of pseudolite layout is shown in Figure 4c. HDOP is less than 6 in 92.8% of the laboratory area, and the mean value of HDOP is 3.1055.
Two sets of experiments were performed. It should be noted that this pseudolite system is in good condition in most cases. In most of the other experiments, high precision positioning results can be obtained in our laboratory. The faults in the following two groups of experiments are accidental. However, these kinds of accidental faults should also be considered in high-precision positioning.The first group of experiments with a zero baseline is named Experiment 1. The purpose is to verify the availability of an indoor positioning system by eliminating all the spatial correlation errors (for example, multipath).Through conducting a zero-baseline experiment, the receiver bias is estimated.A total of 20 epochs are sampled with a sampling interval of 1s. The second group is a short baseline dynamic experiment, including the static stage before moving, named Experiment 2.The sampling interval is set to 1 s, and a total of 123epochs are sampled. Firstly, the static measurement of the rover station at an initial position is maintained and then the rover station is controlled to move back and forth on the fixed rail for two periods. One period is divided into two datasets and a total of four datasets are obtained.

Fault Identification
Since the antenna is shared between the rover and the base stations, the error source such as multipath has been eliminated by DD and the factor affecting positioning only remains noise. Therefore, the DD pseudorange residuals should be close to zero in theory. In Experiment 1, the receiver antenna is placed near (0.64, 0.08) and the positioning results based on carrier phase observations are shown in Figure 5. The color bar represents epochs. According to Figure 5, the ambiguity fixed rate is 25%. However, even if the ambiguities can be fixed in some epochs, it is fixed incorrectly. There is a large deviation in positioning results, and they are not reliable. In Experiments 1, the reference pseudolite is PL5 and DD pseudorange observations are shown in Figure 6. In theory, DD pseudorange residuals in Figure 6 are close to 0, but it is obvious that there are abnormalities in DD observations of all pseudolites. The same circumstance is also found in DD carrier phase observations. It returned to normal around the 10th epoch. Therefore, it can be judged that the reference pseudolite PL5 is faulty. The reason for the bad pseudorangequality of PL5 may be an accidental fault frompseudolite hardware itself.

Fault Exclusion
Two strategies are considered for fault exclusion to ensure the reliability of positioning. The first one called Strategy 1 is to replace the reference pseudolite from PL5 to PL6 and retain the observations of original reference pseudolite. The other one named Strategy 2 is eliminating PL5 directly. In theory, Strategy 1 is more rigorous because subsequent observations return to normal. After the 10th epoch, the STD of PL6-5 for DD pseudorange observations is 0.27 m and 0.04 cycles forthe DD carrier phase. Therefore, after observations return to normal, the signals of PL5 should continue to be received. Figure 7 shows the DD pseudorange and carrier phase observations after using both two strategies. Figure 8 shows the test statistics after applying two strategies. The construction of test statistics and the determination of thresholdare shown in expression (12) and (13).

Test statistic Threshold
Test statistic Threshold According to Figure 7, the DD observations of all pseudolites return to normal after the reference pseudolite is replaced. When using Strategy 1, the test statistics exceed threshold during the previous 10 epochs because the data quality of PL5 is poor, but test statistics remain normal in the last 10 epochs. When using Strategy 2, test statistics are below the threshold since PL5 is removed. It can be considered that the fault has been excluded and high-precision positioning can be performed.
In theory, as long as PL5 returns to normal, the test statistics are less than the set threshold due to a large reduction in residuals of DD observations. In Experiment 1, the faulty pseudolite PL5 is the reference pseudolite, so the method of replacing reference PL is used in Strategy 1. However, even if the reference PL is replaced, PL5 is still in a faulty state in the first 10 epochs and test statistics are higher than the threshold. However, PL5 is eliminated in Strategy 2 and test statistics are less than the threshold, which can ensure that there is no obvious deviation in the whole positioning process. After PL5 returns to normal, more redundant observations will be received in Strategy 1 and the theoretical positioning precision will be higher than Strategy 2. Figure 9 shows the results of horizontal positioning with two strategies and the color bar represents epochs. According to Figure 9, it can be found that the fixed rate has been increased to 50% and the results of the fixed solution are close to true value when using Strategy 1. The STD of horizontal results is 5.4 mm (X direction is 2.7 mm and Y direction 4.7mm). When applying Strategy 2, the fixed rate has been increased to 100% and the STD of horizontal results is 9.8 mm (X direction is 3.7 mm and Y direction 9.1mm).
This phenomenon is also mutually confirmed with the previous analysis. In Strategy 1, when PL5 is faulty, the positioning result shows a large deviation and correct ambiguity fixed solutions cannot be obtained. At present, the positioning results of Strategy 2 are stabilized near true value. After PL5 returns to normal, the positioning result of Strategy 1 is also close to a true value and STD is better than Strategy 2. Nevertheless, we cannot decide which strategy is better without in combination with specific scenarios. The applicable scenarios of the two strategies and the best methods in the dynamic scenario will be discussed in this article.
Combined with Figure 6, the following conclusions can be obtained: if observations return to normal, PL5 can be resumed for positioning. Considering Strategies 1 and 2 comprehensively, Strategy 2 is better from the perspective of a fixed rate. However, considering the precision of the positioning results, Strategy 1 is better. No matter which strategy is adopted, it will be more accurate and reliable than before.
In practical applications, both strategies are acceptable for users according to the results (the precision can reach less than 1cm) and the positioning results are more reliable than before. The availability can be also improved. However, for dynamic positioning (or real-time navigation), Strategy 2 is more suitable than Strategy 1, because Strategy 2 can ensure that a correct positioning result is obtained in each epoch. The continuity and integrity of dynamic positioning will be guaranteed. From the perspective of redundant observations and final positioning precision, the operating state of PL5 (or faulty pseudolite) should be monitored in real-time. Once it returns to normal, the stations should continue to receive its signal, which will make the constellation configuration better. The number of redundant observations can also be added.
The following experiment is a dynamic one and the most reasonable strategy will be determined according to the real-time state of faulty pseudolite.

Fault Identification
A dynamic experiment is carried out in this chapter (named Experiment 2). Firstly, the static measurement of the rover station at the initial position is maintained and then the rover station is controlled to move back and forth on the fixed rail for two periods. Each period possesses two datasets. The horizontal coordinate of the start point is about (1.2, 0.3) and the endpoint coincides with the start point. The positioning results are shown in Figures 10 and 11. The color bar represents the epochs. The red line represents the reference trajectory.
According to Figure 10, it can be found that the positioning results are basically inconsistent with the reference line, which indicates a large deviation comparing with the actual situation. This conclusion can also be drawn from Figure 11. Although the ambiguity fixed rate is 31.7%, it is fixed incorrectly. High-precision positioning cannot be guaranteed, and the results are not reliable.
Similar to Experiment 1, judging whether reference pseudolite performs normally is the first step and reference pseudolite in Experiment 2 is PL6. Figure 12 presents the DD carrier-phase observations of four pseudolites groups. Since Experiment 2 takes "stationary-motionstationary-motion-stationary" as a period, the observations are periodic and contain two cycles in theory.   Figure 12 shows that DD carrier-phase observations of PL6-1, PL6-4, and PL6-5 remain normal, which confirms the regular pattern of two motion periods and observations are stable. The reason for such regular observations is that the pseudolite constellation is fixed. The observations will maintain this periodicity as long as the rover station moves back and forth along the fixed rail. Frequent cycle slips and data interruptions only occur in DD carrier-phase observations of PL6-8, indicating that the reference pseudolite PL6 is normal in Experiment 2 and the signals of PL8 may be faulty. It is possible that some environmental factors (such as multipath, shelters) lead to frequent signals interruption and cycle slips of PL8. However, further fault identification is required to determine faulty pseudolite or epochs theoretically. The test statistics and threshold for fault identification in Experiment 2 are shown in Figure 13. The construction of test statistics and the determination of thresholdare shown in (12) and (13). The CRAIM test statistics in Figure 13a obviously exceeds the threshold. It is less than or close to the threshold in some epochs, where the signal of PL8 is interrupted instead of returning to normal. The number of visible pseudolites (# npl) has also changed from five to four. When the number of visible pseudolites is four, the test statistics wereless than the threshold. At this moment, the receivers lose lock to PL8 signals.It indicates that other pseudolites' observations are normal except PL8. Besides, the incorrect ambiguity fixed solutions in Figures 10 and 11 are also verified, because a wrong ambiguity value can cause slow growth effects in both test statistics and position error. The PRAIM test statistics in Figure 11b also exceed the threshold in very few epochs and it is not as significant as test statistics of CRAIM. According to the introduction, pseudorange is greatly affected by indoor environments, such as multipath, clock, noise and so on. Sometimes the pseudorangebased test statistics cannot truly reflect the fault because test statistics may become very large due to other factors even if there is no fault. More importantly, the PRAIM test statistics can only reflect noise and multipath but cannot recognize cycle slips. Therefore, pseudorange-based RAIM is not

Epoch
Test Statistic Test statistic Threshold applicable in high-precision indoor PL positioning, which illustrates the significance and purpose of this study again.
Since the cycle slip is the most threatening fault for the carrier phase-based positioning system and it can frequently occur in reality, the analysis under cycle slip fault cases is also included. The simulated cycle slips (3 cycles) on PL1-L1, PL4-L1, PL5-L1 and PL6-L1 during the last 20 epochs are added in the dynamic data collection for further analysis. In order to avoid the influence of PL8, PL8 is excluded first, and the simulated cycle slips are added respectively. The corresponding test statistics are shown in Figure 14. In this case, the error source is only simulated cycle slips. According to Figure 14, the data quality of the other four pseudolites except PL8 is good. What's more, the test statistics based on carrier phase observations are sensitive to cycle slips and this kind of fault can be also effectively detected.

Fault Exclusion
In combination with Figures 12 and 13, in order to ensure the accuracy and reliability of positioning, the best way to exclude the fault is to eliminate PL8. Figure 15 shows the test statistics based on pseudorange and carrier phase after PL8 is eliminated.
According to Figure 15, it can be found that the test statistics based on pseudorange and carrier phase return to normal after PL8 areexcluded and the carrier phase-based statistics are much smaller than the threshold. Since the normal DD observations can be received after fault exclusion (PL8 is eliminated). The observations are only affected by general factors such as noise and the test statistics based on posterior residuals are less than the threshold. At the same time, the correctness of AR is verified. Although the pseudorange-based test statistics are below the threshold, they are relatively large, which further confirms the unavailability of pseudorange observations in high-precision positioning for indoor pseudolites. In theory, the best way is to continue to receive the signals of PL8 when it returns to normal. However, it is revealed that PL8 has been in a faulty state during the dynamic process according to Figure 13. Therefore, the signals of PL8 are not received during the whole experiment. In addition, PL8 is normal in the first 18 epochs, but the rover station is in a static state, so whether the PL8 is used at this stage does not affect the positioning results. The result without PL8 is shown in Figure 16. The red line represents the reference trajectory.  According to Figure 16 a,b, the motion trajectory conforms to the actual situation basically that the carrier is controlled to move back and forth along a straight rail. The positioning precision is at centimeter-level in ambiguity fixed solution and the fixed-rate reaches 99.2%. Therefore, it can be judged that PL8 is faulty and the normal observations cannot be obtained during the whole experiment according to the test statistics and threshold. The signals of PL8 should not be received in the whole process. After fault exclusion, all DD observations with good quality can be obtained and indoor positioning results with high-precision can be achieved, which ensures the reliability and availability of indoor pseudolite positioning.

The Horizontal Protection Level
In Experiment 2, HPL is calculated according to position estimate uncertainty, which is obtained by the first two elements of the covariance matrix P. In practical applications, HPL being less than HAL is a prerequisite for RAIM. When the system is operating, the circle centered on the motion position with the HPL as the radius will contain the true position in this epoch. Figure 17 shows the HPL calculated in Experiment 2. The calculation of HPL is obtained by (19). The color bar represents the epochs.  Figure 17a shows the motion trajectory of the first dataset in the second period and the corresponding HPL. In this dataset, all ambiguities are fixed, and the constellation is stable without the replacement of the reference pseudolite. Figure 17b shows the HPL of the full dataset in Experiment 2. The HPL changes within 5 cm and is relatively stable throughout the movement except for the first epoch. The reason is that the ambiguity has not been fixed and the variance is relatively large, which is why the peak appears in Figure 17b. In addition, according to the definition of HPL and the theory of RAIM, it can be considered that HPE is smaller than HPL, that is, the positioning precision of the dynamic experiment is below 5cm. Since there is no true value in positioning, this method also provides a reference index for the accuracy assessment of positioning. Therefore, according to Experiments 1 and 2, it can be considered that the proposed indoor PLRAIM algorithm has the abilitytoensuring the accuracy and integrity of pseudolite system. All test statistics and HPL in Figures 15 and 17 are relatively stable, which also proves the reliability of AR.

Conclusions and Discussion
Due to the influence of various error factors, pseudorange observations are not reliable and the accuracy and reliability of positioning cannot be guaranteed in indoor high-precision pseudolite positioning. Considering these problems, a holistic indoor pseudolite integrity monitoring algorithm PLRAIM is proposed and the specific steps are shown in Section 2.2. The feasibility of the whole algorithm is analyzed by static and dynamic experiments. The ambiguity fixed rate of both experiments is greatly improved when using PLRAIM and the reliability of AR is also guaranteed due to the stability of test statistics and HPL. In terms of positioning precision, the RMS of the zerobaseline static experiment reaches mm level when PLRAIM is applied and the maximum deviation in the dynamic experiment is within 5 cm, which ensures the positioning precision well.
The proposed method can be further improved. In fact, it is not reliable to choose the reference pseudolite of an indoor pseudolite system based on the elevation angle because of the fixed pseudolite constellation and the small scope of the indoor environment. Therefore, the elevation angles of all pseudolites are basically very close. Besides, whether the reference pseudolite is faulty is analyzed based on the stability of DD observations.A more suitable fault identification method for reference pseudolite is the next research focus. In addition, the threshold of test statistics is determined by a given Pfa, but it may not be the best way for an indoor pseudolite system. How to adjust the threshold dynamically according to the actual situation is also the next important issue.
Furthermore, this set of pseudolite hardware and software devices will be applied to some actual engineering situations or more complex indoor environments in the future. More actual data will be collected for a longer period. The available state of pseudolites will be monitored in real-time to verify the robustness of our method.