Next Article in Journal
The Impact of Distributed Propulsion on the Aerodynamic Characteristics of a Blended-Wing-Body Aircraft
Next Article in Special Issue
Performance Impact Assessment of Reducing Separation Minima for En-Route Operations
Previous Article in Journal
Steering a Swarm of Large-Scale Underactuated Mechanical Systems Using a Generalized Coordinates Transformation
Previous Article in Special Issue
Controlling Aircraft Inter-Arrival Time to Reduce Arrival Traffic Delay via a Queue-Based Integer Programming Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Three-Dimensional ANP Evaluation Method Based on Spatial Position Uncertainty under RNP Operation

1
College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, No. 29 General Avenue, Nanjing 211106, China
2
AVIC the First Aircraft Institute, Xi’an 710000, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(11), 703; https://doi.org/10.3390/aerospace9110703
Submission received: 13 September 2022 / Revised: 30 October 2022 / Accepted: 1 November 2022 / Published: 10 November 2022
(This article belongs to the Special Issue Advances in Air Traffic and Airspace Control and Management)

Abstract

:
Performance-based navigation (PBN) operations based on the required navigation performance (RNP) operations are the trend for civil aviation in the future. In order to further ensure the safety and efficiency of civil aviation operation, RNP will transition from 2D to 3D/4D. Accurately evaluating the actual navigation performance in three-dimensional directions (3D ANP) under civil aircraft RNP operation is important to guarantee the safe flight of civil aircraft. However, the traditional two-dimensional ANP evaluation method mainly focuses on plane navigation performance along the track, lacks the evaluation of the vertical direction. Moreover, the traditional three-dimensional accuracy evaluation method is based on the assumption of three-dimensional independence and uniformity, which can only carry out approximate calculation, and the evaluation result is inaccurate. Therefore, this paper constructs a three-dimensional ellipsoid error probability (EEP) evaluation model for the spatial position uncertainty of the navigation output of the flight management system in three-dimensional directions, and gives the three-dimensional accurate ANP calculation results through iterative numerical integration. The simulation results show that the evaluation method proposed in this paper can accurately evaluate the actual navigation performance of the airborne navigation system in all directions of 3D space, and has higher evaluation accuracy and precision than the traditional ANP evaluation method, which is of great significance to ensure the flight safety of civil aircraft under 3D/4D RNP operation in the future.

1. Introduction

The increasing diversification of air transport and the continuous growth of air transport volume have put forward new demands. The air transport system urgently needs to be upgraded and innovated upon. The International Civil Aviation Organization (ICAO), the Federal Aviation Administration (FAA), the European Aviation Safety Agency (EASA) and Civil Aviation Administration of China (CAAC) actively promote the performance-based navigation (PBN) mode based on the required navigation performance (RNP) operations to enhance the management capability of the future air traffic control four-dimensional track (4DT) operation [1,2,3]. The development of RNP operations and the innovation of airborne navigation systems have freed civil aircraft from the dependence on ground-based navigation facilities, made flight procedures more flexible, maximized the use of airspace, further refined the design of flight safety interval standards and further improved the safety, efficiency, capacity, access, flexibility, predictability and resilience of the National Airspace System (NAS). At the same time, these have reduced the impact of air transport on the environment.
The required navigation performance is developed from the rules for implementation of area navigation (RNAV), which is a RNAV flight mode with airborne performance monitoring and warning functions. When implementing RNP procedures, performance monitoring and alarm functions must be provided in different forms according to the different installations, structures and layouts of airborne navigation systems [4]:
(1)
Real-time actual navigation performance evaluation and prediction;
(2)
Monitor the airborne navigation performance, provide an alarm when it fails to meet the RNP requirements of the current airspace;
RNP uses total system error that includes the navigation system error (NSE), flight technical error (FTE) and, path definition error (PDE) as the measurement standards for navigation capability. In the designated airspace, as long as the aircraft meets the required navigation performance requirements, it is allowed to fly in the airspace. The limit value of RNP is expected to be achieved for at least 95% of the flight time.
ICAO has specified the minimum navigation performance standard for civil aircraft for different stages of RNP operation [5]. As shown in Figure 1, in the RNP specification, the concept of a safety tunnel is used to quantitatively describe the specific requirements of all flight phases of the aircraft from take-off, climb, cruise, descent, approach and landing, that is, each flight phase has a virtual tunnel, which specifies the vertical and horizonal boundary of the flight path perimeter of the flight phase, and the aircraft navigation performance must meet the RNP requirements of the flight phase.
The flight management system takes the actual navigation performance (ANP) as an index to evaluate the accuracy of the navigation system, so as to judge whether the aircraft is in a safe tunnel [6]. At this stage, the ANP model uses the probability circle radius of estimated position uncertainty (EPU) to represent the actual navigation performance [7], and uses the 95% position error circle radius in the horizontal direction as the actual navigation performance of civil aircraft during RNP flight. This method reduces the three-dimensional position information of civil aircraft to two-dimensional [8,9,10,11,12], which has the problem of lack of vertical information and cannot characterize the three-dimensional actual navigation performance of civil aircraft.
At the same time, at present, the spherical error probability (SEP) model is commonly used to evaluate the spatial position accuracy, that is, 50% of the error sphere radius is used as the evaluation result of the three-dimensional position accuracy. However, during the RNP operation, FMS uses different sensors to provide fusion navigation information. The GNSS outputs three-dimensional position information, but presents different performance in the horizontal and vertical directions due to constellation distribution [13]. The land-based radio navigation system (DME/VOR) only provides horizontal position. The barometric altimeter only provides height information [14]. This results in the unequal variance of the three-dimensional position error of FMS output; the SEP model cannot obtain an accurate solution, so only approximate calculations can be carried out [15]. In addition, the spatial position uncertainty of civil aircraft during RNP flight is also reflected in the correlation of navigation results in the three-dimensional direction. At this time, using the SEP model for ANP evaluation will further reduce the evaluation accuracy and effectiveness.

2. Related Work

2.1. Traditional Two-Dimensional ANP Model and Calculation Method

Traditional ANP calculation only considers the position uncertainty caused by a single sensor and mainly focuses on the navigation and positioning performance of the aircraft along the track direction.
The horizontal ANP evaluation is regarded as the calculation process of the horizontal position error probability circle radius with the estimated position of the main navigation sensor as the center of the circle.
E represents the mean value of the position-estimated error of the FMS along the horizontal plane.
E = [ μ x μ y ] T
V is the covariance matrix of horizontal position estimation error for flight management system.
V = σ x 2 σ x y σ x y σ y 2 = σ x 2 ρ σ x σ y ρ σ x σ y σ y 2
Then the joint probability density function of horizontal position error can be expressed as:
f x , y = 1 2 π V 1 / 2 exp 1 2 E T V 1 E = 1 2 π σ x σ y 1 ρ 2 exp ( x μ x ) 2 σ x 2 + ( y μ y ) 2 σ y 2 2 ρ ( x μ x ) ( y μ y ) σ x σ y 2 ( 1 ρ 2 )
The correlation coefficient is:
ρ = σ x y σ x σ y
The major and minor axes of the error ellipse after rotation transformation coincide with the coordinate axis. At this time, the joint probability density function can be expressed as:
f m x , y = 1 2 π det V m 1 / 2 exp 1 2 E T V m 1 E = 1 2 π σ m x σ m y exp x 2 2 σ m x 2 y 2 2 σ m y 2
The covariance matrix becomes:
V m = σ m x 2 σ m y 2 = Λ
Within:
σ m x = l m a j o r σ m y = l m a j o r
The probability P of the real position falling into the two-dimensional normal distribution integral circle is:
P = S 1 2 π det V m 1 / 2 exp 1 2 E T V m 1 E d S
As shown in Figure 2, the ANP value estimated by the position of the flight management system at this time can be obtained by calculating the radius of the integral circle.
Let us consider:
x = r cos θ y = r sin θ
Then Equation (8) becomes:
P = 0 2 π 0 k σ m y 1 2 π σ m x σ m y exp 1 2 ( r 2 cos 2 θ σ m x 2 + r 2 sin 2 θ σ m y 2 ) r d r d θ
For specific probability (P = 95%) and given covariance matrix V, the solution of ANP is converted to the solution of scale coefficient k, but Equation (10) can only be obtained through numerical calculation. Therefore, the traditional method approximates k (P = 95%), as follows:
k = 0.4852 r a t i o 3 + 1.9625
Within:
r a t i o = σ m y / σ m x
ANP is expressed as:
A N P = k * l m a j o r

2.2. Traditional Three-Dimensional ANP Model and Calculation Method

For the RNP operation, FMS needs to monitor the navigation performance in both the horizontal and vertical direction. The spherical error probability (SEP) model is often used to describe the accuracy of the falling point of an object in three-dimensional space, which can be used to evaluate the 3D ANP.
The position estimation error (x, y, z) of the FMS in the X-Y-Z three-dimensional space follows normal distribution x N ( μ x , σ x 2 ) , y N ( μ y , σ y 2 ) , z N ( μ z , σ z 2 ) .
The spatial three-dimensional joint probability density function is:
f x , y , z = 1 ( 2 π ) 3 / 2 σ x σ y σ z exp 1 2 ( x μ x ) 2 σ x 2 + ( y μ y ) 2 σ y 2 + ( z μ z ) 2 σ z 2
The traditional SEP model considers that the variances of the three dimensions are equal and independent of each other, and the mean value is zero, which is μ x = μ y = μ z = 0 , σ x = σ y = σ z = σ .
Let:
x = r sin α cos β y = r sin α sin β z = r cos α
Then the probability is:
P = 2 R ( 2 π ) 1 / 2 σ exp ( R 2 2 σ 2 ) + 2 ( 2 π ) 1 / 2 σ 0 R exp ( r 2 2 σ 2 ) d r
The approximate solution of the specified probability can be obtained by numerical integration.
As shown in Figure 3, the approximate solution of SEP (95%) obtained by Grubbs [16] method is:
R ( 95 % ) = 1.65 × σ m ( v / 9 + ( 1 v / 9 ) ) 3 / 2
Within:
σ m = σ x 2 + σ y 2 + σ z 2
v = 2 ( σ x 4 + σ y 4 + σ z 4 ) σ m 4

3. Method

3.1. Three-Dimensional Ellipsoid Error Probability EEP Model

Considering the spatial position uncertainty of civil aircraft flying, this paper establishes a three-dimensional ellipsoidal EEP model and performs three-dimensional ANP calculation on this basis.
Use E to represent the mean value of the FMS estimation error in three-dimensional space:
E = [ μ x μ y μ z ]
V is the covariance matrix of the 3D position estimation error of the FMS:
V = σ x 2 σ x y 2 σ x z 2 σ x y 2 σ y 2 σ y z 2 σ x z 2 σ y z 2 σ z 2 = σ x 2 ρ x y σ x σ y ρ x z σ x σ z ρ x y σ x σ y σ y 2 ρ y z σ y σ z ρ x z σ x σ z ρ y z σ y σ z σ z 2
where ρ x y , ρ x z , ρ y z are correlation coefficients:
ρ x y = σ x y σ x σ y ρ x z = σ x z σ x σ z ρ y z = σ y z σ y σ z
The joint probability density function of the three-dimensional position error (x, y, z) is:
f x , y , z = 1 ( 2 π ) 3 / 2 det V 1 / 2 exp 1 2 E T V 1 E
Similarity transformation is performed on the third-order covariance matrix:
V = A Λ A 1
A is the similar transformation matrix V, the diagonal matrix can be obtained:
Λ = λ 1 λ 2 λ 3
Within λ 1   λ 2   λ 3 are the eigenvalues of the covariance matrix, and the covariance matrix becomes uncorrelated:
V m = σ m x 2 σ m y 2 σ m z 2 = Λ
The joint probability density function can be expressed as:
f x , y , z = 1 ( 2 π ) 3 / 2 σ m x σ m y σ m z exp 1 2 x 2 σ m x 2 + y 2 σ m y 2 + z 2 σ m z 2
The probability P that the real position falls into the three-dimensional normal distribution integral ellipsoid domain is:
P = Ω 1 ( 2 π ) 3 / 2 det V m 1 / 2 exp 1 2 E T V m 1 E d Ω
As shown in Figure 4, for a given probability P, the three-dimensional ANP calculation process can be converted into the calculation of the three-axis length of the three-dimensional normal distribution integral ellipsoid Ω.

3.2. Three-Dimensional ANP Calculation Method

FMS provides positioning function based on navigation information fusion. This paper mainly considers the ANP model and calculation method after obtaining the position covariance matrix from the FMS. The influence of time-varying sensor performance on covariance matrix is beyond the scope of this paper.
For a given probability P, the covariance matrix, which considers the contribution of each sensor, can be obtained by the integrated navigation calculation module of the FMS. It can be seen from the Formula (27), when P = 0.95, the ANP is determined by σ m x , σ m y and σ m z :
x = r sin α cos β y = r σ y σ x sin α sin β z = r σ z σ x cos α
Within α 0 , π , β 0 , 2 π , the joint probability density function is transformed into spherical coordinates:
P = σ y σ z σ x 2 0 R 0 π 0 2 π f ( r sin α cos β , r sin α sin β , r cos β ) r 2 sin α d r d α d β
The probability integration becomes:
P = 1 ( 2 π ) 3 / 2 σ m x 3 0 R r 2 exp ( r 2 2 σ m x 2 ) d r 0 π sin α d α 0 2 π d β = 2 R ( 2 π ) 1 / 2 σ m x 3 0 R r 2 exp ( r 2 2 σ m x 2 ) d r = 2 R ( 2 π ) 1 / 2 σ m x exp ( R 2 2 σ m x 2 ) + 2 ( 2 π ) 1 / 2 σ m x 0 R exp ( r 2 2 σ m x 2 ) d r
Let w = R / σ m x , then:
P + 1 2 = Φ ( w ) w ( 2 π ) 1 / 2 exp ( w 2 2 )
where Φ ( w ) is the standard normal distribution function. For a given P, Equation (32) cannot obtain an analytical solution. In this paper, an approximate solution is obtained by iterative numerical integration, as shown in Figure 5:
.
Take the value to calculate the allowable error ξ = 10 6 ( P 0.95 < ξ ),when R / σ m x = 2.79548 , the probability is 95%. It can be obtained that the three-axis radius lengths of the integral ellipsoid domain are: R x = 2.79548 σ m x , R y = 2.79548 σ m y , R z = 2.79548 σ m z . The ANP of the FMS in three-dimensional space can be represented by R x , R y and R z .
As shown in Figure 6, the volume of the EEP model is smaller than that of the SEP model, and the disparity increases rapidly with the increase of the ratio of three-axis standard deviation. That is, under the same probability, the EEP model is closer to the real position error and has higher evaluation accuracy.
At the same time, the projection of the EEP model and SEP model in the XOY, XOZ, and YOZ planes are compared, respectively. As shown in Figure 7, Figure 8 and Figure 9, the EEP model more accurately represents the difference of position error variance of FMS position estimation in every direction, and realizes the accurate evaluation of the three-dimensional direction of position uncertainty.
The calculation process of three-dimensional ANP evaluation method based on the EEP model is shown in Figure 10 (Algorithm 1).
Algorithm 1: Calculation process of three-dimensional ANP evaluation method
Step 1: Get the position error covariance matrix V from the FMS.
Step 2: Similarity transform the position error covariance V to obtain the positive definite covariance matrix Vm.
Step 3: Use numerical integration iterative to obtain the inclusion probability P of the EEP model.
Step 4: If P 0.95 < ξ , proceed to the next step; otherwise, return to step 3.
Step 5: Take the three-axis radius Rx, Ry, Rz of the integrated ellipsoid domain are the ANP of the FMS in three-dimensional directions.

4. Experimental Setup

According to the objectives and mission characteristics of the ANP verification, the civil aircraft RNP flight multi-source information fusion simulation and verification platform is constructed, as shown in Figure 11.
The simulation and verification platform was composed of a civil aircraft RNP flight simulation module, multi-sensor navigation data generation module, multi-source navigation information fusion, actual navigation performance evaluation module, and flight display module. It can simulate the output of civil aircraft airborne multi-source redundant configuration sensors under RNP flight, simulate the information processing of FMS, fuse and process the multi-source redundant navigation information, and evaluate the ANP of the FMS so as to verify the method proposed in this paper.
Simulation with actual RNP flight operation [17], which can meet the simulation function requirements of typical civil aviation flight environment.
The RNP flight procedure from Beijing Capital Airport (ZBAA) to Shanghai Hongqiao International Airport (ZSSS) is executed in the simulation. The initial position is 40.073333° N, 116.598333° E at an altitude of 30.4 m, the initial heading angle is 175°, and the end point is 31.196667° N, 121.335000° E at an altitude is 8 m. The simulation duration is 7081 s. The corresponding waypoint positions along the route are shown in Table 1, and the fly track simulation output is shown in Figure 12.
The parameters of airborne navigation sensors are set according to the minimum aviation system performance standards for RNP operation [18]. The bias of gyros is set as 0.01°/h and the bias of accelerometers is set as 10−4 g. The GNSS (GPS, Beidou) uses the real ephemeris for simulation. The pseudo range measurement error is set to 10.5 m [19]. The positions of the land-based navigation stations are shown in Table 2. The barometric altimeter is also used to aid vertical navigation [20].

5. Results

The three-dimensional position estimated error of FMS is compared with the traditional ANP method, SEP method, and the EEP method that is proposed in this paper.
Figure 13, Figure 14 and Figure 15 show the comparison results of the position-estimated error of the FMS in the longitude, latitude, and vertical direction with the ANP evaluation values.
The comparison of actual navigation performance evaluation algorithm results are shown in Table 3.
In order to compare the evaluation accuracy and precision of the traditional ANP model, SEP model, and EEP model proposed in this paper, the recall and mean absolute percentage error (MAPE) statistics are shown in Table 4 and Table 5.
Statistical analysis shows that the EEP model proposed in this paper has no false evaluation compared with the SEP model, and the evaluation accuracy in longitude direction and latitude direction is effectively improved compared with the traditional ANP evaluation method and SEP evaluation method.
In order to verify the effectiveness of the proposed algorithm, 50 groups of Monte Carlo simulations were carried out under the same simulation conditions. The simulation results are shown in Figure 16.
For the 50 groups of Monte Carlo simulation results, the evaluation accuracy and precision of the traditional ANP model, SEP model, and EEP model proposed in this paper are statistically analyzed in Table 6 and Table 7.
Statistical analysis shows that the EEP model proposed in this paper has no false evaluation, and the evaluation accuracy in the longitude direction is 43.62% higher than that of the traditional ANP evaluation method and 34.91% higher than that of the SEP evaluation method. In the latitude direction, it is 34.55% higher than that of the traditional ANP evaluation method and 13.28% higher than that of the SEP evaluation method.

6. Conclusions

During RNP operation, it is important to monitor the actual navigation performance of FMS accurately. Traditional ANP calculation mainly focuses on the navigation and positioning performance of the aircraft along the track direction and considers the position uncertainty caused by a single sensor. The lack of comprehensive consideration of spatial position uncertainty will lead to the loss of dimension perception of aircraft navigation performance and inaccurate evaluation of ANP.
In this paper, to meet the needs of future development of civil aircraft RNP to 3D/4D, aiming at the lack of evaluation dimension of the traditional two-dimensional actual navigation performance evaluation method and the inaccuracy of the SEP model, a three-dimensional EEP model is proposed based on the uncertainty of the three-dimensional position of civil aircraft with unequal variances and relevance, and the corresponding three-dimensional ANP calculation method is given.
(1)
Compared with the traditional two-dimensional ANP model and SEP model, it fully considers the spatial uncertainty of civil aircraft, has smaller model inclusion volume, and improves the accuracy of the evaluation model.
(2)
The simulation results show that the EEP model proposed in this paper can effectively improve the evaluation accuracy on the basis of ensuring the accuracy of the evaluation of the actual navigation performance of the FMS in three-dimensional directions and more accurately represent the uncertainty of the position of the FMS in the flight space–time of the civil aircraft RNP to ensure flight safety.

Author Contributions

Conceptualization, Y.D. and J.L.; methodology, Y.D. and J.L.; validation, Y.D., Z.L. and J.L.; formal analysis, Y.D. and X.S.; investigation, Y.D., B.Z. and J.L.; data curation, Y.D., P.L. and J.L.; writing—original draft preparation, Y.D.; writing—review and editing, Y.D., Z.L. and J.L.; visualization, J.L.; supervision, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundamental Research Funds for the Central Universities (NS2022039, NF2022001), and National Natural Science Foundation of China (61973160).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

PBNPerformance Based Navigation
RNPRequired Navigation Performance
RNAVArea Navigation
ANPActual Navigation Performance
EEPEllipsoid Error Probability
SEPSpherical Error Probability
EPUEstimated Position Uncertainty
FMSFlight Management System
DMEDistance Measuring Equipment
VOR VHFOmnidirectional Radio
MAPEMean Absolute Percentage Error

References

  1. Joint Planning and Development Office (JPDO). Concept of Operations for the Next Generation Air Transportation System; JPDO: Washington, DC, USA, 2013. [Google Scholar]
  2. SESAR Consortium. The ATM Target of Operations; SESAR Consortium: Brussels, Belgium, 2007. [Google Scholar]
  3. Civil Aviation Administration of China. Statistical Bulletin on the Development of Civil Aviation Industry of 2019; Development Planning Division: Beijing, China, 2020.
  4. International Civil Aviation Organization (ICAO). Doc 9613: Performance Based Navigation (PBN) Manual; ICAO: Montréal, QC, Canada, 2013. [Google Scholar]
  5. Radio Technical Commission for Aeronautics (RTCA). RTCA DO-283A, Minimum Operational Performance Standards for Required Navigation Performance for Area Navigation; RTCA Inc.: Washington, DC, USA, 2013. [Google Scholar]
  6. FAA Aeronautical Center. Assessing the Validity of Using Actual Navigation Performance (ANP) Information for Supporting Designated Flight Inspection Operations; FAA Aeronautical Center: Oklahoma City, OK, USA, 2008.
  7. Yoon, M.G.; Kim, J.K. Evaluation methodology for safety maturity in air navigation safety. J. Air Transp. Manag. 2022, 98, 102159. [Google Scholar] [CrossRef]
  8. Fu, L.; Zhang, J.; Li, R. Real-time total system error estimation: Modeling and application in required navigation performance. Chin. J. Aeronaut. 2014, 27, 1544–1553. [Google Scholar] [CrossRef] [Green Version]
  9. Zhao, H.; Xu, X.; Zhang, J.; Zhu, Y.; Yang, C.; Hong, S. Lateral Flight Technical Error Estimation Model for Performance Based Navigation. Chin. J. Aeronaut. 2011, 24, 329–336. [Google Scholar] [CrossRef] [Green Version]
  10. Chamlou, R. TIS-B: Calculation of navigation accuracy category for position and velocity parameters. In Proceedings of the 23rd Digital Avionics Systems Conference, Salt Lake City, UT, USA, 28 October 2004; IEEE: Piscataway, NJ, USA Computer Science. ; pp. D.3-1–D.3-13. [Google Scholar]
  11. Carl, E.J. Analysis of Fatigue-Crack Propagation and Fracture Data: AIAA-2009-1363; AIAA: Reston, VA, USA, 1973. [Google Scholar]
  12. Zhang, Q.; Zhang, J.; Zhu, Y.; Liu, W. Navigation Accuracy Category-Position Models and Estimate Position Uncertainty Calculations for TIS-B System. Chin. J. Aeronaut. 2009, 22, 419–425. [Google Scholar] [CrossRef] [Green Version]
  13. Lu, D.; Schnieder, E. Performance Evaluation of GNSS for Train Localization. IEEE Trans. Intell. Transp. Syst. 2015, 16, 1054–1059. [Google Scholar] [CrossRef]
  14. Liu, M.; Lai, J.; Li, Z.; Liu, J. An adaptive cubature Kalman filter algorithm for inertial and land-based navigation system. Aerosp. Sci. Technol. 2016, 51, 52–60. [Google Scholar] [CrossRef]
  15. Przemieniecki, J.S. Mathematical Methods in Defense Analyses; American Institute of Aeronautics and Astronautics Inc.: Washington, DC, USA, 2000. [Google Scholar]
  16. Gura, A.; Gersten, R.H. On Analysis of n-Dimensional Normal Probabilities; Space and missile system organization of air force system command Los Angeles air force station: Los Angeles, CA, USA, 1970. [Google Scholar]
  17. Pamplona, D.A.; de Barros, A.G.; Alves, C.J. Performance-Based Navigation Flight Path Analysis Using Fast-Time Simulation. Energies 2021, 14, 7800. [Google Scholar] [CrossRef]
  18. Radio Technical Commission for Aeronautics (RTCA). RTCA DO-236C, Minimum Aviation System Performance Standards: Re-quired Navigation Performance for Area Navigation; Radio Technical Commission for Aeronautics (RTCA): Washington, DC, USA, 2014. [Google Scholar]
  19. Radio Technical Commission for Aeronautics (RTCA). RTCA DO-229C, Minimum Operational Performance Standards for Global Positioning System/Wide Area Augmentation System Airborne Equipment; Radio Technical Commission for Aeronautics (RTCA): Washington, DC, USA, 2001. [Google Scholar]
  20. Aeronautical Radio Inc. (ARINC). ARINC 712-7, Airborne ADF System; ARINC: Annapolis, MD, USA, 1997. [Google Scholar]
Figure 1. RNP for each leg.
Figure 1. RNP for each leg.
Aerospace 09 00703 g001
Figure 2. Two-dimensional ANP model.
Figure 2. Two-dimensional ANP model.
Aerospace 09 00703 g002
Figure 3. Spherical Error Probable.
Figure 3. Spherical Error Probable.
Aerospace 09 00703 g003
Figure 4. Three-dimensional position error ellipsoid.
Figure 4. Three-dimensional position error ellipsoid.
Aerospace 09 00703 g004
Figure 5. Numerical integration results of P and R / σ m x .
Figure 5. Numerical integration results of P and R / σ m x .
Aerospace 09 00703 g005
Figure 6. Comparison between SEP volume and EEP model volume.
Figure 6. Comparison between SEP volume and EEP model volume.
Aerospace 09 00703 g006
Figure 7. Comparison between EEP model and SEP model in XOY plane.
Figure 7. Comparison between EEP model and SEP model in XOY plane.
Aerospace 09 00703 g007
Figure 8. Comparison between EEP model and SEP model in XOZ plane.
Figure 8. Comparison between EEP model and SEP model in XOZ plane.
Aerospace 09 00703 g008
Figure 9. Comparison between EEP model and SEP model in YOZ plane.
Figure 9. Comparison between EEP model and SEP model in YOZ plane.
Aerospace 09 00703 g009
Figure 10. Three-dimensional ANP calculation process.
Figure 10. Three-dimensional ANP calculation process.
Aerospace 09 00703 g010
Figure 11. RNP flight simulation and verification platform.
Figure 11. RNP flight simulation and verification platform.
Aerospace 09 00703 g011
Figure 12. Flight simulation trajectory.
Figure 12. Flight simulation trajectory.
Aerospace 09 00703 g012
Figure 13. Comparison between longitude position error and ANP.
Figure 13. Comparison between longitude position error and ANP.
Aerospace 09 00703 g013
Figure 14. Comparison between latitude position error and ANP.
Figure 14. Comparison between latitude position error and ANP.
Aerospace 09 00703 g014
Figure 15. Comparison between altitude position error and ANP.
Figure 15. Comparison between altitude position error and ANP.
Aerospace 09 00703 g015
Figure 16. Statistical comparison chart of Monte Carlo simulation.
Figure 16. Statistical comparison chart of Monte Carlo simulation.
Aerospace 09 00703 g016
Table 1. Simulation track waypoint table.
Table 1. Simulation track waypoint table.
Way PointLongitude/°Latitude/°
ZBAA116.59833340.073333
ELKUR116.66500038.639722
DOXAB116.62500038.291111
OVNUG116.69222238.110278
PANKI116.99972237.295833
GUSIR117.04166737.206667
YQG117.21500036.833333
DALIM117.21305636.417778
ABTUB117.36777836.000556
UDINO117.80333334.822500
DPX117.99833334.278333
OMUDI118.26666733.970556
LAGAL118.70555633.467222
ATVAD118.85416733.295833
SUBKU118.94027833.197222
NIXEM119.16055632.941944
XUTGU119.57555632.463889
PIMOL119.76166732.245833
VMB120.19166731.743333
SASAN120.31944431.589444
ZSSS121.33500031.196667
Table 2. Location information of ground stations.
Table 2. Location information of ground stations.
Station NumberLongitude/°Latitude/°Altitude/m
1116.139.637
2117.239.1328117
3118.28566135.22202549
4118.88205031.73621945
5121.32593131.14360565
Table 3. Comparison of actual navigation performance evaluation algorithm results.
Table 3. Comparison of actual navigation performance evaluation algorithm results.
Statistical Items (RMSE)Longitude/mLatitude/mAltitude/m
Position error0.5570.59431.6269
Traditional ANP3.01123.0112-
SEP2.49912.49912.4991
EEP1.65272.27993.4975
Table 4. Comparison of actual navigation performance evaluation algorithm results.
Table 4. Comparison of actual navigation performance evaluation algorithm results.
RecallLongitudeLatitudeAltitude
Traditional ANP100%100%-
SEP100%100%97.14%
EEP100%100%100%
Table 5. Comparison of actual navigation performance evaluation algorithm results.
Table 5. Comparison of actual navigation performance evaluation algorithm results.
MAPELongitudeLatitudeAltitude
Traditional ANP50.728337.9040-
SEP46.423034.013845.4766
EEP31.613130.303755.5729
Table 6. Precision comparison of actual navigation performance evaluation algorithm results in Monte Carlo simulation.
Table 6. Precision comparison of actual navigation performance evaluation algorithm results in Monte Carlo simulation.
RecallLongitudeLatitudeAltitude
Traditional ANP100%100%-
SEP100%100%97.14%
EEP100%100%100%
Table 7. Accuracy comparison of actual navigation performance evaluation algorithm results in Monte Carlo simulation.
Table 7. Accuracy comparison of actual navigation performance evaluation algorithm results in Monte Carlo simulation.
MAPELongitudeLatitudeAltitude
Traditional ANP50.728337.9040-
SEP46.423034.013845.4766
EEP31.613130.303755.5729
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dai, Y.; Lai, J.; Zhu, B.; Li, Z.; Sun, X.; Lv, P. Three-Dimensional ANP Evaluation Method Based on Spatial Position Uncertainty under RNP Operation. Aerospace 2022, 9, 703. https://doi.org/10.3390/aerospace9110703

AMA Style

Dai Y, Lai J, Zhu B, Li Z, Sun X, Lv P. Three-Dimensional ANP Evaluation Method Based on Spatial Position Uncertainty under RNP Operation. Aerospace. 2022; 9(11):703. https://doi.org/10.3390/aerospace9110703

Chicago/Turabian Style

Dai, Yuting, Jizhou Lai, Bin Zhu, Zhimin Li, Xiao Sun, and Pin Lv. 2022. "Three-Dimensional ANP Evaluation Method Based on Spatial Position Uncertainty under RNP Operation" Aerospace 9, no. 11: 703. https://doi.org/10.3390/aerospace9110703

APA Style

Dai, Y., Lai, J., Zhu, B., Li, Z., Sun, X., & Lv, P. (2022). Three-Dimensional ANP Evaluation Method Based on Spatial Position Uncertainty under RNP Operation. Aerospace, 9(11), 703. https://doi.org/10.3390/aerospace9110703

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