Next Article in Journal
Structural Strength Analysis and Optimization of Commercial Aircraft Nose Landing Gear under Towing Taxi-Out Conditions Using Finite Element Simulation and Modal Testing
Previous Article in Journal
The Impact of Upstream Static Deformation on Flow Past a Cylinder/Flare
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Refined Aircraft Positioning Based on Stochastic Hybrid Estimation with Adaptive Square-Root Unscented Particle Filtering

by
Yangyang Zhang
1,
Zhenxing Gao
2,*,
Kai Qi
3 and
Jiawei Li
1
1
College of Civil Aviation, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
2
College of General Aviation and Flight, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
3
Department of Flight, Shandong Airlines, Jinan 250014, China
*
Author to whom correspondence should be addressed.
Aerospace 2024, 11(5), 413; https://doi.org/10.3390/aerospace11050413
Submission received: 22 March 2024 / Revised: 14 May 2024 / Accepted: 17 May 2024 / Published: 20 May 2024
(This article belongs to the Section Aeronautics)

Abstract

:
The positioning of civil aviation aircraft relative to a geographic reference point on Earth in a Cartesian frame is significant to detect the deviations from the desired path, especially for high-altitude airports or special airports based on performance-based navigation (PBN). To obtain these critical deviations during aircraft approach and landing, it is fundamental to estimate the continuous flight variables and discrete flight modes simultaneously with enough accuracy. With the coordinate conversion between the North, East, and Down (NED) frame and the geographic coordinate system based on World Geodetic System 1984 (WGS-84) considered, this study proposed a non-linear stochastic hybrid estimation algorithm with adaptive square-root unscented particle filtering (ASR-UPF) to estimate the true path. The probabilities of mode transition, represented by the normal cumulative density function of continuous states, determine whether to proceed with mode transitions. In addition, the adaptive update characterized by tracking variable noise and the importance sampling distributions based on the results of square-root unscented Kalman filtering (SR-UKF), as a comparative study of continuous system filtering, were used. The experiments illustrated the ASR-UPF is able to reduce the state estimation error more effectively, and more promptly track the error caused by incorrect mode estimation with adaptability compared to the SR-UKF. A further test with real flight data indicates that the proposed method gives the refined estimation of position and azimuth in NED frame.

1. Introduction

The deviations from the desired path of a civil aviation aircraft are significant for accident analysis, such as collision with obstacles, which needs positioning of the aircraft relative to obstacles, especially for high-altitude airports or special airports based on performance-based navigation (PBN) [1] that are positioned using global positioning systems (GPSs) with waypoints. Meanwhile, its positioning error is smaller than the protection area of traditional instrument flight procedure, and there are many tall obstacles around the airport. Consequently, it is necessary to estimate the position and azimuth in the Cartesian frame. Unfortunately, the position and azimuth of an aircraft relative to a geographic reference point on Earth in a Cartesian frame are not recorded in flight data, which provides support for solving the problem of positioning estimation in a Cartesian frame. When the aircraft is a highly complex nonlinear dynamic system, the continuous motion of the aircraft can make transitions between flight modes. Correspondingly, different modes can also change the patterns of the aircraft’s continuous motion, that is, variable aerodynamic configuration. Meanwhile, the flight data inevitably have uncertain noise due to uncertain factors such as navigation error, sensor measurement error, and external disturbances [2], and when integrated with manual operations, this could cause stochastic transitions between flight modes [3]. Therefore, the model of single-continuous-flight dynamics cannot match the variable aerodynamic configuration caused by mode transition, and a hybrid dynamics model should be established for discrete state and continuous state estimation.
Hybrid dynamics modeling mainly appears in the study of hybrid systems with coupled discrete and continuous states [4]. The earliest hybrid systems modeling formalized the underlying behavior and manual operations of the system using formal method [5]. However, this method, which has poor universality, uses strict formal language to express discrete mode transition, which has no description of continuous behavior. The mixed logical dynamic model (MLDM) contains the logic of continuous motion within discrete state transition logic, describing the transition of discrete states as a linear inequality containing both discrete and continuous variables via Boolean algebra and determines the mode transition based on the true or false logical conditions [6,7]. Nevertheless, MLDM characterized by deterministic probabilities of discrete state transition cannot reflect the uncertainty of flight mode transition. The generalized fuzzy hidden Markov model (GFHMM), extending HMM to the fuzzy domain, obtains the mode transition of automatic flight systems by fuzzy inference with two parameters, transition fuzzy density, and emission fuzzy density, which can be obtained through parameter statistical analysis or clustering algorithms [8,9]. However, the GFHMM does not take into account the mode transition caused by manual operations, and it is essentially a method for continuous time series analysis and prediction.
To demonstrate the stochastic mode transition, interacting multiple model (IMM) and constrained IMM take into account the probability of constant mode transition [10,11,12,13], which is limited to stochastic hybrid systems with discrete mode transitions that are independent of the continuous state variation. The state-dependent transition hybrid estimation (SDTHE) algorithm explicitly calculates the mode transition probabilities that depend on the continuous states [14]. Compared with IMM, SDTHE is more suitable for state estimation for general stochastic hybrid systems. Scholars have used the SDTHE algorithm to infer the intentions of the automatic flight system and compare them with the pilot’s operational intentions to identify abnormal human–machine interactions [15]. By constraining the continuous state, the constrained SDTHE algorithm was developed for the estimation of unmanned aerial vehicle trajectory with turning constraints [16].
In cooperation with state estimation algorithms, this system can carry out state estimation based on the stochastic hybrid dynamic model with coupled discrete and continuous states. The extended Kalman filter (EKF), unscented Kalman filter (UKF), and cubature Kalman filter (CKF) based on KF have limited improvement of estimation accuracy for nonlinear systems [17,18,19], while particle filter could effectively solve this problem using sequential Monte Carlo sampling with the appropriate number of particles [20,21]. However, the importance of sampling distribution commonly uses a particle filter, which is generally a prior distribution, instead of using the posterior probability distribution corrected by measurement data, which would cause issues of low efficiency and sensitivity to singular points and would not effectively utilize measurement data [22]. Meanwhile, taking the posterior probability distribution generated by the square root unscented Kalman filter (SR-UKF) [23] as the importance sampling distribution not only utilizes new measurement data but also maintains computational stability in the form of the square root of covariance. In addition, due to the uncertainty of noise during filtering, adaptive methods need to be adopted to update the noise, which could improve the estimation accuracy of flight data further [24].
Considering that the aircraft is a highly nonlinear system, whose modes are numerous and complex, there is generally a problem of dimension explosion. Even though some solutions for simplifying state equations and measurement equations via aircraft mode abstraction have been put forward [14,25], it would induce oversimplification if only the motion of the center of mass in a Cartesian frame were considered. In other words, calculating the coordinates of an aircraft in the Cartesian frame directly instead of converting coordinates between geographic and Cartesian frame would lead to inaccurate position calculations, which further cause excessive estimation error.
To estimate the position and azimuth of civil aviation aircraft with coupled modes and continuous states is fundamental to flight safety analysis and accident investigation. Focusing on the hybrid estimation of a nonlinear stochastic hybrid system, this study contributes to developing a synthetic filtering system characterized by the adaptive square-root unscented particle filtering (ASR-UPF) algorithm and considering the coordinate transformation via the track projection between a Cartesian frame and a geographic frame based on the World Geodetic System 1984 (WGS-84) to ensure high accuracy of the hybrid model [26]. In a simulation study and experiments on real flight data, the stochastic nonlinear hybrid estimation with ASR-UPF is tested and discussed deeply in this study.
This paper is organized as follows. In Section 2, the methodology of synthetic filtering system integrated with ASR-UPF and track projection is described, which is employed in hybrid system modeling and its estimation algorithm. In Section 3, the proposed algorithm is demonstrated with two examples of simulation tests and experiments on real flight data. The discussion and conclusions are presented in Section 4 and Section 5, respectively.

2. Methodology

A synthetic filtering system integrated with ASR-UPF and track projection is developed to estimate the relative position and azimuth of an aircraft in a Cartesian frame or local inertial frame, whose orientation is North, East, and Down x n , y n , z n with the defined origin.

2.1. Adaptive Square-Root Unscented Particle Filtering

For continuous states, the state equation at time k 1 and the measurement equation at time k are defined as
x k q   = f   q x k 1 q + ω k 1 q y k = h   q x k q + r k q
where q M is the mode and M is the set consisting of N m modes. ω k 1 q , r k q represent mutually independent Gaussian white noise at time k 1 and k, respectively.
At time k = 0 , we define the initial mode as q 0 , of which the probability is α 0 , q 0 + = 1 . The distribution of the initial continuous state corresponds to N x ^ 0 , q 0 + , Q 0 q 0 , from which sample N particles x ^ 0 , q 0 + , n p , n p = 1 , 2 , 3 , , N , and Q 0 q 0 is the initial system noise covariance matrix obtained from the statistical analysis of historical flight data. The weights and square root of the covariance matrices of each particle are obtained as
w ˜ 0 q 0 , n p = 1 N S 0 q 0 , n p = Q 0 q 0
For the filter of other modes, we define the probabilities, the mean, and the square root of the covariance matrices of the continuous states as
α 0 , q + = 0 x ^ 0 , q + = 0 S 0 q   = 0 q M , q q 0
Thus, the sampled particles and their weights of other modes are
x ^ 0 , q + , n p = 0 S 0 q   , n p = 0 w ˜ 0 q   , n p = 1 N
At time k 1 , the posterior distribution of each particle in each filter is N x ^ k 1 , q k 1 + , n p , S k 1 q k 1 , n p T S k 1 q k 1 , n p . The posterior probability of each discrete mode is α k 1 , q k 1 + , and the probability of mode transitions is γ k q k 1 q k .
As the algorithm basis of ASR-UPF, the SR-UKF is described in the Appendix A1. To improve estimation accuracy, the adaptive filtering scheme is proposed, in which the two noise covariance matrices Q k q k , n p and R k q k , n p of particle n p in mode q k (or the q k t h filter) are able to be updated according to the error statistics. Then, the measurement noise covariance matrix is adjusted based on the a priori estimation in Equation (4) as
R k q k , n p = ( 1 d k ) R k 1 q k 1 , n p + d k [ y ˜ k q k 0 , n p y ˜ k q k 0 , n p T r o w = 0 2 n W r o w q k 0 , n p c ( y ^ k q k 0 , n p ( r o w ) y ^ k q k 0 , n p ) ( y ^ k q k 0 , n p ( r o w ) y ^ k q k 0 , n p ) T ]
where the weights d k [ 0 , 1 ] determine the window length in adaptive filtering. With equal weights, the estimation will gradually become accurate if the noise characteristics remain unchanged. However, the weights with exponential fading are more appropriate to deal with the variation in the noise characteristics:
d k = ( 1 b ) / ( 1 b k + 1 )
where b ( 0 , 1 ) is the forgetting factor. Therefore, the square root of the measurement noise covariance matrix is derived as
R k q k , n p T R k q k , n p = ( 1 d k ) ( R k 1 q k 1 , n p ) T d k y ˜ k q k 0 , n p ( 1 d k ) R k 1 q k 1 , n p T d k y ˜ k q k 0 , n p r o w = 1 2 n d k W r o w q k 0 , n p ( c ) ( y ^ k q k 0 , n p ( r o w ) y ^ k q k 0 , n p ) d k W r o w q k 0 , n p ( c ) ( y ^ k q k 0 , n p ( r o w ) y ^ k q k 0 , n p ) T s g n ( W 0 q k 0 , n p ( c ) ) d k W 0 q k 0 , n p ( c ) 2 ( y ^ k q k 0 , n p ( 0 ) y ^ k q k 0 , n p ) ( y ^ k q k 0 , n p ( 0 ) y ^ k q k 0 , n p ) T
In computer implementation, the QR decomposition is carried out to obtain R k q k , n p , and the Cholesky factor of the rank-1 update is carried out twice as
R k q k , n p = q r ( 1 d k ) R k 1 q k 1 , n p , d k y ˜ k q k 0 , n p T T R k q k , n p = c h o l u p d a t e ( R k q k , n p , d k W 1 : 2 n q k 0 , n p ( c ) y ^ k q k 0 , n p ( 1 : 2 n ) y ^ k q k 0 , n p , ) R k q k , n p = c h o l u p d a t e ( R k q k , n p , d k W 0 q k 0 , n p ( c ) y ^ k q k 0 , n p ( 0 ) y ^ k q k 0 , n p , ( s g n ( W 0 q k 0 , n p ( c ) ) ) ) R k q k , n p = d i a g ( d i a g ( R k q k , n p ) )
where d i a g represents the diagonal matrix. After the adaptive update of R k q k , n p , the measurement update is performed based on traditional SR-UKF, as shown in Equation (A5). To complete the adaptation operation, the process noise covariance matrix Q k q k , n p , updated based on the a priori estimation in Equation (A3) and a posteriori estimation in Equation (A6), is adjusted as follows:
Q k q k , n p = 1 d k Q k 1 q k 1 , n p + d k x ^ k q k 0 , n p + x ^ k q k 0 , n p x ^ k q k 0 , n p + x ^ k q k 0 , n p T U k q k , n p U k q k , n p T
After the QR decomposition and Cholesky factor of rank-1 update,
Q k q k , n p T Q k q k , n p = 1 d k Q k 1 q k 1 , n p d k ( x ^ k q k 0 , n p + x ^ k q k 0 , n p ) T T ( 1 d k ) Q k 1 q k 1 , n p ( d k ( x ^ k q k 0 , n p + x ^ k q k 0 , n p ) ) T d k U k q k , n p ( d k U k q k , n p ) T
In computer implementation, the square root of the process noise covariance matrix is derived by
Q k q k , n p = q r ( 1 d k ) Q k 1 q k 1 , n p , ( d k ( x ^ k q k 0 , n p + x ^ k q k 0 , n p ) ) T T Q k q k , n p = c h o l u p d a t e ( Q k q k , n p , d k U k q k , n p , ) Q k q k , n p = d i a g ( d i a g ( Q k q k , n p ) )
Define N x ^ k q k 0 , n p + , S k q k 0 , n p S k q k 0 , n p T according to Equation (A6) as the importance sampling distribution, from which we obtain x ^ k , q k + , n p particles and calculate the weight of each particle. Then, the weight of each particle is updated and normalized as follows:
w k q k , n p = w ˜ k 1 q k 1 , n p N y k ; x ^ k , q k + , n p , R k q k , n p R k q k , n p T N x ^ k , q k + , n p ; x ^ k q k 0 , n p + , S k q k 0 , n p S k q k 0 , n p T N x ^ k , q k + , n p ; x ^ k 1 , q k 1 + , n p , Q k 1 q k 1 , n p Q k 1 q k 1 , n p T w ˜ k q k , n p = w k q k , n p n p = 1 N w k q k , n p
We determine whether to perform resampling based on the number of valid particles and calculate the number of valid particles as follows [27]:
N ^ e f f = 1 n p = 1 N w k q k , n p 2
If N ^ e f f < N t h r e s h o l d , where N t h r e s h o l d represents the threshold of particle number, we perform random resampling and make w ˜ k q k , n p = 1 N . The algorithm pseudocode is shown in Algorithm 1.    
Algorithm 1: Adaptive Square-Root Unscented Particle Filter
Aerospace 11 00413 i001

2.2. Stochastic Hybrid Estimation with ASR-UPF

As shown in Figure 1. For discrete states, the set of conditions for mode transition is
D S t ( q k 1 , q k ) = [ x k 1 , q k 1 , ζ ] T | ( L q k 1 q k [ x k 1 T ζ T ] T 0 | q k 1 )
where L q k 1 q k   is a constant matrix, and is ζ the logic judgment variable of mode transition, which refers to the combination of continuous flight parameters that affect mode transition and could be obtained from historical flight data through data statistics. The posterior probability of mode q k 1 is given as follows:
α k 1 , q k 1 + = p q k 1 | Y k 1
where Y k 1 = y 1 , y 2 , , y k 1 represents the time series of measurements. q k 1 represents the mode at time k 1 .
Thus, the conditional posterior probability of a continuous state at time k 1 is N μ q k 1 q k , Σ k 1 q k 1 q k . We define conditional mode transition probability as [15]:
π q k 1 q k S t = p ( q k | q k 1 , x k 1 ) = p ( [ x k 1 , q k 1 , ζ ] T D S t ( q k 1 , q k ) ) = a q k 1 q k   + b q k 1 q k   Φ L q k 1 q k   x k 1   ζ q k 1 q k , Σ k 1 q k 1 q k
where Φ · represents the normal cumulative density function, and ζ q k 1 q k is the judgment variable for mode transition logic, and its distribution can be obtained by statistical analysis of flight data (e.g., we obtain the continuous flight parameters related to mode transition of a certain aircraft model, and obtain the distribution of these parameters in a parametric test). Σ k 1 q k 1 q k is the constant matrix of ζ q k 1 q k . a q k 1 q k   , b q k 1 q k   are scalar constants such that π q k 1 q k S t 0 for all q k 1 , q k M , and q k = 1 N m π q k 1 q k S t = 1 .
Then, the probability of mode transition is derived as follows:
γ k q k 1 q k = p q k | q k 1 , Y k 1 = a q k 1 q k   + b q k 1 q k   Φ L q k 1 q k   x ^ k 1 , q k 1 + μ q k 1 q k ; Σ ˜ k 1 q k 1 q k Σ ˜ k 1 q k 1 q k = Σ k 1 q k 1 q k + L q k 1 q k   P k 1 q k 1 L q k 1 q k  
where μ q k 1 q k is the mean of ζ q k 1 q k and P k 1 q k 1 is the covariance matrix of x ^ k 1 , q k 1 + , which can be obtained as follows:
x ^ k 1 , q k 1 + = n p = 1 N w ˜ k 1 q k 1 , n p x ^ k 1 , q k 1 + , n p Σ ˜ k 1 q k 1 q k = Σ k 1 q k 1 q k + L q k 1 q k   S k 1 q k 1 T S k 1 , q k 1 L q k 1 q k T P k 1 q k 1 = n p = 1 N w ˜ k 1 q k 1 , n p x ^ k 1 , q k 1 + , n p x ^ k 1 , q k 1 + x ^ k 1 , q k 1 + , n p x ^ k 1 , q k 1 + T + n p = 1 N w ˜ k 1 q k 1 , n p S k 1 q k 1 , n p T S k 1 q k 1 , n p S k 1 q k 1 = q r n p = 1 N w ˜ k 1 q k 1 , n p x ^ k 1 , q k 1 + , n p x ^ k 1 , q k 1 +   , n p = 1 N w ˜ k 1 q k 1 , n p S k 1 q k 1 , n p T
Initial value of each filter can be obtained as follows:
x ^ k 1 q k 0 , n p = q k 1 = 1 n m x ^ k 1 , q k 1 + , n p α ¯ k 1 q k q k 1 P k 1 q k 0 , n p = q k 1 = 1 n m S k 1 q k 1 , n p T S k 1 q k 1 , n p + x ^ k 1 , q k 1 + , n p x ^ k 1 q k 0 , n p x ^ k 1 , q k 1 + , n p x ^ k 1 q k 0 , n p T α ¯ k 1 q k q k 1 S k 1 q k 0 , n p = q r α ¯ k 1 q k q k 1 q k 1 = 1 n m S k 1 , | q k 1 T , q k 1 = 1 n m x ^ k 1 , | q k 1 + x ^ k 1 q k 0 T α ¯ k 1 q k q k 1 = p q k 1 | q k , Y k 1 = γ k 1 q k 1 q k α k 1 , q k 1 + q k 1 = 1 n m γ k 1 q k 1 q k α k 1 , q k 1 +
Further, we obtain the filtering values for each mode using the ASR-UPF filter based on Equation (19):
x ^ k q k + = n p = 1 N w ˜ k q k , n p x ^ k , q k + , n p P k q k = n p = 1 N w ˜ k q k , n p x ^ k q k + x ^ k , q k + , n p x ^ k q k + x ^ k , q k + , n p T + S k q k 0 , n p S k q k 0 , n p T
Meanwhile, the residual and its covariance matrix can be obtained as follows:
y ˜ k q k 0 = n p = 1 N w ˜ k q k , n p y k y ^ k q k 0 , n p P y , k q k = n p = 1 N w ˜ k q k , n p S y , k q k , n p S y , k q k , n p T
The mode prior probability is derived as follows:
α k , q k = p q k | Y k = q k 1 = 1 n m γ k 1 q k 1 , q k α k 1 , q k 1 +
We define likelihood function of the q k t h filter as follows:
Λ k , q k   = p y k | q k , Y k 1 = N y ˜ k q k 0 ; 0 , P y , k q k
After that, the posterior probability of mode q k and the mode at time k can be derived as follows:
α k , q k + = p y k | q k , Y k 1 p q k | Y k 1 q k = 1 n m p y k | q k , Y k 1 p q k | Y k 1 = Λ k , q k   α k , q k q k = 1 n n Λ k , q k   α k , q k q k = arg max q k α k , q k +
Thus, the posterior probability of continuous state and its covariance can be obtained as follows:
x ^ k + = j = 1 n n α k , q k + x ^ k q k + P k   = q k = 1 n m α k , q k + x ^ k q k + x ^ k + x ^ k q k + x ^ k + T + P k q k

2.3. Track Projection

Taking the flight of an aircraft from the initial approach fix (IAF) to the final approach fix (FAF) as an example, we define the aircraft as having a heading-hold (HH) mode and a turning (TN) mode in terms of lateral, and ot has constant height (CH) mode and descent (DT) mode terms longitudinally. Thus, there are four modes defined in this study, containing m 1 = H H , C H , m 2 = T N , C H , m 3 = H H , D T , and m 4 = T N , D T . The set of modes can be described as M = m 1 , m 2 , m 3 , m 4 .
For simplification, the spherical terrestrial reference was used to calculate the position in a Cartesian frame instead of an ellipsoidal one. Choosing x = x n , y n , z n , η n T as the state vector, according to laws of kinematics, we can easily obtain x n k , y n k , z n k by given x n k 1 , y n k 1 , z n k 1
x n k , y n k , z n k T = x n k 1 , y n k 1 , z n k 1 T + C x ˙ n k 1 , y ˙ n k 1 , z ˙ n k 1 , x ¨ n k 1 , y ¨ n k 1 , z ¨ n k 1 T
where C is a matrix that varies with the modes, which will be discussed in Section 3 in detail. It is believed that the projection only changes the azimuth rather than the length of the arc between two points. Thus, the azimuth η n k by given x n k 1 , y n k 1 , z n k 1 and η n k 1 can be obtained as follows:
η n k = η n k 1 + Δ η n k 1 Δ η n k 1 = arccos cos L 1 cos L 2 cos L 3 / sin L 2 sin L 3
where the second equation is the cosine law with respect to sides of a spherical triangle [28], while L 1 , L 2 , and L 3 represent the length of the three sides, which can be calculated as shown in Figure 2. Δ η n k 1 is the angle corresponding to side L 1 .
Then, the state equation can be built as follows:
x k = f x k 1 , u k 1 + w k 1
where the input matrix u k 1 = x ˙ n k 1 , y ˙ n k 1 , z ˙ n k 1 , x ¨ n k 1 , y ¨ n k 1 , z ¨ n k 1 T .
According to Great Circle projection shown in Figure 2. If the latitudes and longitudes of points A and B are given as ϕ A , λ A and ϕ B , λ B , respectively. The azimuth η from A to B can be obtained as follows:
cos ( A O B ) = cos ( π 2 ϕ B ) cos ( π 2 ϕ A ) + sin ( π 2 ϕ B ) sin ( π 2 ϕ A ) cos ( λ B λ A ) η = arcsin ( sin ( π 2 ϕ B ) × sin ( λ B λ A ) sin ( A O B ) )
where A O B represents the angle between O A and O B . Thus, The Cartesian coordinates of point B relative to point A can be obtained as follows:
x B t o A = L 2 cos ( π 2 η ) y B t o A = L 2 sin ( π 2 η )
where L 2 represents the great circle distance between A and B, which can be calculated as L 2 = R c .
We define the reference point, which is the origin in Cartesian frame, as point A; thus, the coordinate of any point B relative to point A according to Equations (29) and (30) can be described as
x B t o A , y B t o A , η T = h ϕ B , λ B
As mentioned above, we can obtain the latitude and longitude of point B by giving the latitude and longitude of point A, the azimuth η , and the great circle distance L 2 . In a short period of time T s , it is believed that the flight distance in a Cartesian frame could be approximately equal to L 1 . That means that
L 1 = T s x ˙ n k 1 + 1 2 x ¨ n k 1 T s 2 2 + T s y ˙ n k 1 + 1 2 y ¨ n k 1 T s 2 2 1 2
Consequently, the measurement input can be defined as y = ς , γ , h T while for all modes, the measurement equations are all the same as
y k = h 1 x k + r k

3. Results and Discussion

3.1. Scenario and Parameter Settings

We take the instrument approach procedure, the base turn-a-specific pattern of reversal procedure, of a certain aircraft on a certain runway with elevation 3567.5 m and magnetic variation 0.3 W as an example. As shown in Figure 3, during the approach, the aircraft would first hold the bearing 107 from the IAF at 7200 m and descend by 1400 m, and then turn to intercept the inbound track at 274   while descending by 150 m. Finally, it descends to 5350 m and then holds the altitude with the constant bearing until the FAF.
For the leg from point IAF to point WPT1, the aircraft descends at a constant speed. After descending with variant speed of the three-axis NED frame for the leg from WPT1 to WPT2, it then descends at an even speed until the point at nearly 5350 m, where it finally keeps a level flight until the FAF. Thus, we obtain the matrix C of different modes in Equation (26).
C 1 = T s 0 0 0 0 0 0 T s 0 0 0 0 0 0 0 0 0 0 f o r m 1 f o r m 2 C 3 = T s 0 0 0 0 0 0 T s 0 0 0 0 0 0 T s 0 0 0 f o r m 3 C 4 = T s 0 0 T s 2 / 2 0 0 0 T s 0 0 T s 2 / 2 0 0 0 T s 0 0 T s 2 / 2 f o r m 4
where T s = 0.1 s.
We collect data of the altitude H and azimuth η at points WPT1, WPT2, and the point where the aircraft has just been performing a level flight, that is, where the mode transition occurs. It is assumed that these data are independent of each other and follow a normal distribution. Statistical analysis is conducted on these data separately to obtain the means μ H and μ η , as well as the standard deviations σ H and σ η . The results are shown in Table 1 in the International System of Units (SI).
According to the approach procedure, the series of modes is m 3 m 4 m 3 m 1 , so the probability of mode transition is built as shown in Table 1 for simplification. It could indeed be simplified like this with only the path from IAF to FAF considered. The probability of the mode transition γ q k 1 q k is determined by a q k 1 q k   , b q k 1 q k   , L q k 1 q k   , μ q k 1 q k , Σ k 1 q k 1 q k , x ^ k 1 q k 1 and P k 1 q k 1 . Meanwhile, x ^ k 1 q k 1 and P k 1 q k 1 can be obtained from the results by ASR-UPF, and here, a q k 1 q k   = 0 and b q k 1 q k   = 1 were used because of their designability [15]. The variable μ q k 1 q k and its variance Σ k 1 q k 1 q k are obtained from Table 1.
We only consider the mode transitions in Table 2. For the transition from m 3 to m 4 , the azimuth η and altitude H would decrease, and we define L 34 = 0 0 1 0 0 0 0 1 , μ 34 = 2239.767 , 108.127 T and Σ 34 = 3.046 2 0 0 1.862 2 . Thus, L 34 x ^ k 1 , | i   μ 34 0 when the mode transition happens. For the same mode transition from m 4 to m 3 , we define L 43 = 0 0 1 0 0 0 0 1 , μ 43 = 2084.072 , 274.845 T and Σ 43 = 2.536 2 0 0 1.387 2 . For the transition from m 3 to m 1 , it is sufficient to consider the altitude only with L 31 = 0 0 1 0 , μ 31 = 1783.758 and Σ 31 = 1.653 2 .

3.2. Estimation Results

The following tests were conducted to verify the performance of the nonlinear stochastic hybrid estimation of the synthetic filtering system. It is undesirable to test the algorithm by comparing the estimation results with recorded flight data directly because the measurements from in-flight data are not accurate in terms of noise and some unexpected disturbances and cannot be regarded as the true value, and there is even a lack of records fir some data. Instead, a flight dynamics model (FDM) of B737-800 with motion models, aerodynamic models, engine models, and control system models, etc., was used to provide standard flight variables for comparison.
Based on the approach procedure mentioned above, the estimation results of ASR-UPF with 100, 300 and 500 particles and the SR-UKF were compared to each other to illustrate the estimation accuracy. The two noise covariance matrices in the FDM can be set freely because of the robustness in the adaptive scheme [29]. While or fast convergence, the proper values were set with Q 0 = d i a g 1 , 1 , 1 , 0.005 and R 0 = d i a g 0.005 , 0.005 , 1 . The forgetting factor in ASR-UPF was set by b = 0.9 . Figure 4 shows the estimation error of the position and azimuth by SR-UKF (green line) and the ASR-UPF with 100 (black line), 300 (cyan line), and 500 particles (red line). Both estimation errors were almost the same initially. However, the estimation error of ASR-UPF started convergence at about 3.9 s, 3.0 s, and 2.8 s corresponding to ASR-UPF with 100, 300, and 500 particles, respectively, while the error of the SR-UKS starts convergence at about 6.8s and subsequently ramps up in the filtering process. In addition, the x n , y n , z n (or H) and η estimation errors of both algorithms increased significantly at about 399 s and 560 s because there is a delay in mode transition estimation. And the estimation error of z n by both algorithms increased at about 628 s due to the mode transition from m 3 to m 1 with a deceleration descent process. In addition, the estimation errors of ASR-UPF with 300 and 500 particles are close.
The comparison results of mode estimation with the various algorithms are shown in Figure 5, Figure 6 and Figure 7. Obviously, in Figure 5, there is a bit of delay in mode transitions by both algorithms, while the mode probability variation of ASR-UPF in Figure 6 is more stable than that of SR-UKF in Figure 7, which oscillates slightly. Since the performance of ASR-UPF with 300 and 500 particles is close, we just consider ASR-UPF with 300 particles here. The results of ASR-UPF display a better performance in mode estimation regardless of the number of particles and show more accurate mode estimation and more stable mode transitions as well. We also conducted another 50 simulation experiments, and the average number of mode transition error and the root mean square error (RMSE) of the position and azimuth for the respective algorithm are given in Table 3. This shows that the RMSE of ASR-UPF is much smaller than that of SR-UKF, as well as the mode transition error, while the RMSEs by ASR-UPF with 300 particles and 500 particles are close.

3.3. Results of ASR-UPF Based on Real Flight Data

The synthetic filtering system and ASR-UPF were applied to real flight-data-based nonlinear hybrid estimation. There is a problem to be addressed before the estimation based on real flight data. The system update and measurement update may not be performed simultaneously because the sampling rates of related variables are probably different. Under normal circumstances, the system update and measurement update are carried out simultaneously. When there are no measurements, only the system update is executed by unequal interval filtering [30].
Still taking the approach procedure in Figure 3 as an example, the measurements and estimated results according to the flight data are shown in Figure 8. Results indicate that the trend of the estimated height approximates that of the recorded flight data, though there were deviations between the estimated z n and the recorded values. The horizontal coordinates in the NED frame and azimuth without magnetic variation of the aircraft relative to VOR, not the measured in-flight data, could be estimated by the synthetic filtering system as shown in Figure 8a,b. According to the test, the convergence of ASR-UPF was slower than the simulation results by approximately 11s, but the subsequent error slightly increased, as shown in Figure 8c.

4. Discussion

To improve the estimation accuracy, this study provides a hybrid estimation of positioning, in a Cartesian frame, of an aircraft based on flight data with different modes or patterns of the aircraft’s continuous motion considered. The mode estimation, characterized by state-dependent transitions, ensures the accurate modeling for the continuous system. Meanwhile, the ASR-UPF introducing the SR-UKF to generate the importance sampling distribution fully utilizes measurements and guarantees the convergence while maintaining the computational stability of the filter. Otherwise, the adaptive update by tracking noise for reducing errors caused by modeling and measuring improves the robustness to the stochastic noise and enables aircraft positioning with better accuracy.
Though the proposed algorithm is validated based on real flight data, the estimation accuracy needs to be improved further. Also, the resolution and accuracy of inertial measurements, used as the control variable, depend on different flight-data-acquisition systems. Therefore, the estimation performance based on the inertial measurements of a specific aircraft needs to be explored in depth. In addition, this study only considered the transition between two modes because of the specific approach procedure, although there were four modes that were artificially defined. More complex circumstances with transitions between more than two modes, like changing the flight schedule due to some unexpected events (such as encountering windshear), should be studied in depth.

5. Conclusions

This study proposes a nonlinear stochastic hybrid filtering system integrated with the Great Circle projection and adaptive square root unscented particle filter. Compared with other research, more accurate state and measurement equations are provided through the projection from longitudes and latitudes in the WGS-84 frame to coordinates in the NED frame, although existing errors are caused by using a spherical reference system for coordinate conversion. The ASR-UPF is recommended for use within position and azimuth estimation applications based on flight data. Since the inertial measurements are insusceptible to external disturbance, without loss of generality, the proposed algorithm could be applied for the refinement of other variables that are not recorded in flight data as well. As a data-driven method, it is a suitable candidate for discrete and continuous state estimation of civil aviation aircraft in specific scenarios.

Author Contributions

Data curation, K.Q.; Formal analysis, K.Q. and J.L.; Funding acquisition, K.Q.; Investigation, Z.G. and Y.Z.; Methodology, Y.Z. and Z.G.; Project administration, Z.G.; Resources, K.Q.; Software, Y.Z.; Supervision, Z.G.; Validation, Y.Z. and J.L.; Visualization, Y.Z. and J.L.; Writing—original draft, Y.Z. and Z.G.; Writing—review and editing, Y.Z. and Z.G. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Natural Science Foundation of China (U2333202, 52272351), the Aeronautical Science Foundation of China (2022Z066052002), and the Civil Aviation Administration of China (ASSA2023/22).

Data Availability Statement

Restrictions apply to the availability of these data. Data were obtained from Shandong Airline and are available from Kai, Qi with the permission of Shandong Airline.

Acknowledgments

The authors acknowledge the flight data provided by Shandong Airlines.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PBNPerformance Based Navigation
GPSGlobal Positioning System
NEDNorth, East, and Down
WGS-84World Geodetic System 1984
ASR-UPFAdaptive square-root unscented particle filter
SR-UKFSquare-root unscented Kalman filter
MLDMMixed logical dynamic model
GFHMMGeneralized fuzzy hidden Markov model
IMMInteracting multiple model
SDTHEState-dependent transition hybrid estimation
EKFExtended Kalman filter
UKFUnscented Kalman filter
CKFCubature Kalman filter

Appendix A

Appendix A.1

The fundamental square root unscented Kalman filtering algorithm, from which the ASR-UPF is derived, is described as follows. The SR-UKF is initialized according to Equation (19).
At each time step t k 1 , a group of 2 n + 1 sigma points are selected as
Ø k 1 q k 0 0 + , n p = x ^ k 1 q k 0 , n p Ø k 1 q k 0 r o w + , n p = x ^ k 1 q k 0 , n p + n + λ S k 1 q k 0 , n p , r o w , r o w = 1 , . . . , n Ø k 1 q k 0 n + r o w + , n p = x ^ k 1 q k 0 , n p n + λ S k 1 q k 0 , n p , r o w , r o w = 1 , . . . , n
where S k 1 q k 0 , n p , r o w represents the r o w t h -column of square root matrix S k 1 q k 0 , n p . A group of weights related to each sigma point is set by
W 0 q k 0 , n p ( m ) = λ n + λ W 0 q k 0 , n p ( c ) = λ n + λ + 1 p g 2 + p h W r o w q k 0 , n p ( m ) = W r o w q k 0 , n p ( c ) = 1 2 ( n + λ ) r o w = 1 , 2 , , 2 n
where W r o w q k 0 , n p ( m ) and W r o w q k 0 , n p ( c ) represent the weights related to the mean and covariance of x ^ each. The constant p a [ 0.0001 , 1 ] determines the spread of the sigma points around x ^ . The constant p b , used to incorporate prior knowledge of the distribution of x ^ , is set as p b = 2 for Gaussian distribution. λ = n ( p a 2 1 ) is the scaling parameter. The update process of SR-UKF is described by
x ^ k q k 0 , n p ( r o w ) = f ( Ø k 1 q k 0 ( r o w ) , n p + ) x ^ k q k 0 , n p = r o w = 0 2 n W r o w q k 0 , n p m x ^ k q k 0 , n p ( r o w ) S k q k 0 , n p = q r r o w = 0 2 n W r o w q k 0 , n p c ( x ^ k q k 0 , n p ( r o w ) x ^ k q k 0 , n p ) T , Q k 1 q k 1 , n p T S k q k 0 , n p = c h o l u p d a t e ( S k q k 0 , n p W 0 q k 0 , n p ( c ) ( x ^ k q k 0 , n p ( r o w ) x ^ k q k 0 , n p ) T , s g n ( W 0 q k 0 , n p ( c ) ) )
where cholupdate represents the Cholesky factor update and sgn represents the sign of W 0 q k 0 , n p ( c ) . As the square root of the error covariance matrix, S k q k 0 , n p is computed by QR decomposition and Cholesky factor update successively. Note that if S is the original Cholesky factor of P = A A T , then the Cholesky factor of the rank-1 update (or downdate) P ± m m T is denoted as S = c h o l u p d a t e ( S , m , ± ) . If m is a matrix with M columns, the result is M consecutive updates of the Cholesky factor by using the M columns of m .
The next step is to carry out the measurement update. Another group of 2 n + 1 sigma points is selected according to the priori estimation x ^ k q k 0 , n p as
Ø k q k 0 , n p ( 0 ) = x ^ k q k 0 , n p Ø k q k 0 , n p ( r o w ) = x ^ k q k 0 , n p + n + λ S k q k 0 , n p r o w , r o w = 1 , , n Ø k q k 0 , n p ( n + r o w ) = x ^ k q k 0 , n p n + λ S k q k 0 , n p r o w , r o w = 1 , , n
The priori estimation of measurements is
y ^ k q k 0 , n p ( r o w ) = h ( Ø k q k 0 , n p ( r o w ) ) y ^ k q k 0 , n p = r o w = 0 2 n W r o w q k 0 , n p m y ^ k q k 0 , n p ( r o w )
The measurement update is performed as
S y , k q k , n p = q r r o w = 0 2 n W r o w q k 0 , n p ( c ) y ^ k q k 0 , n p ( r o w ) y ^ k q k 0 , n p R k q k , n p S y , k q k , n p = c h o l u p d a t e ( S y , k q k , n p , W 0 q k 0 , n p ( c ) ( y ^ k q k 0 , n p ( 0 ) y ^ k q k 0 , n p ) , s g n ( W 0 q k 0 , n p ( c ) ) ) P x y , k q k , n p = r o w = 0 2 n W r o w q k 0 , n p ( c ) ( x ^ k q k 0 , n p ( r o w ) x ^ k q k 0 , n p ) ( y ^ k q k 0 , n p ( r o w ) y ^ k q k 0 , n p ) T K k q k , n p = P x y , k q k , n p S y , k q k , n p ( S y , k q k , n p ) T 1 x ^ k q k 0 , n p + = x ^ k q k 0 , n p + K k q k , n p ( y k y ^ k q k 0 , n p ) U k q k , n p = K k q k , n p S y , k q k , n p S k q k 0 , n p = c h o l u p d a t e ( S k q k 0 , n p , U k q k , n p , )
where the square root of the measurement noise covariance matrix S y , k q k , n p is obtained by QR decomposition and Cholesky factor of rank-1 as well. P x y , k q k , n p is the cross-covariance between x ^ k q k 0 , n p and y ^ k q k 0 , n p . K k q k , n p is the filtering gain. With the above operations, the square root of error covariance matrix S k q k 0 , n p is obtained.

References

  1. ICAO. ICAO 9613 Performance-Based Navigation (PBN) Manual; ICAO: Montreal, QC, Canada, 2013; Available online: https://standards.globalspec.com/std/1687641/icao-9613 (accessed on 16 May 2024).
  2. Pretto, M.; Giannattasio, P.; De Gennaro, M. Mixed analysis-synthesis approach for estimating airport noise from civil air traffic. Transp. Res. Part D Transp. Environ. 2022, 106, 103248. [Google Scholar] [CrossRef]
  3. Nandiganahalli, J.S.; Lee, S.; Hwang, I. Formal Verification for Mode Confusion in the Flight Deck Using Intent-Based Abstraction. J. Aerosp. Inf. Syst. 2016, 13, 343–356. [Google Scholar] [CrossRef]
  4. Menci, M.; Papi, M.; Smarrazzo, F. Existence of solutions for coupled hybrid systems of differential equations for mi-croscopic dynamics and local concentrations. Commun. Pure Appl. Anal. 2023, 22, 2146–2168. [Google Scholar] [CrossRef]
  5. Degani, A.; Shafto, M.; Kirlik, A. Modes in Human-Machine Systems: Constructs, representation and classification. Int. J. Aviat. Psychol. 1999, 9, 125–138. [Google Scholar] [CrossRef]
  6. Sirmatel, I.I.; Geroliminis, N. Mixed logical dynamical modeling and hybrid model predictive control of public transport operations. Transp. Res. Part B Methodol. 2018, 114, 325–345. [Google Scholar] [CrossRef]
  7. Pataro, I.M.L.; Gil, J.D.; da Americano, C.M.V.; Guzmán, J.L.; Berenguel, M. A nonlinear control approach for hybrid solar thermal plants based on operational conditions. Renew. Energy 2022, 183, 114–129. [Google Scholar] [CrossRef]
  8. Lyu, H.; Nandiganahalli, J.S.; Hwang, I. Human Automation Interaction Issue Detection Using a Generalized Fuzzy Hidden Markov Model. In Proceedings of the AIAA Information Systems-AIAA Infotech at Aerospace, Grapevine, TX, USA, 9–13 January 2017. [Google Scholar] [CrossRef]
  9. Lyu, H.; Nandiganahalli, J.S.; Hwang, I. Automation intent inference using the GFHMM for flight deck mode confusion detection. J. Aerosp. Inf. Syst. 2018, 15, 172–177. [Google Scholar] [CrossRef]
  10. Arroyo, C.A.; Asensio, V.M. Adaptive IMM-UKF for Airborne Tracking. Aerospace 2023, 10, 698. [Google Scholar] [CrossRef]
  11. Yin, J.; Yang, Z.; Luo, Y. Adaptive tracking method for non-cooperative continuously thrusting spacecraft. Aerospace 2021, 8, 244. [Google Scholar] [CrossRef]
  12. Ebrahimi, M.; Ardeshiri, M.; Khanghah, S.A. Bearing-only 2D maneuvering target tracking using smart interacting multiple model filter. Digit. Signal Process. 2022, 126, 103497. [Google Scholar] [CrossRef]
  13. Xu, W.; Xiao, J.; Xu, D.; Wang, H.; Cao, J. An Adaptive IMM Algorithm for a PD Radar with Improved Maneuvering Target Tracking Performance. Remote Sens. 2024, 16, 1051. [Google Scholar] [CrossRef]
  14. Seah, C.E.; Hwang, I. Stochastic Linear Hybrid Systems: Modeling, Estimation, and Application in Air Traffic Control. IEEE Trans. Control Syst. Technol. 2009, 17, 563–575. [Google Scholar] [CrossRef]
  15. Lee, S.; Hwang, I.; Leiden, K. Intent Inference-Based Flight-Deck Human-Automation Mode-Confusion Detection. J. Aerosp. Inf. Syst. 2015, 12, 503–518. [Google Scholar] [CrossRef]
  16. Lee, J.; Hwang, I.; Shim, D.H. UAS surveillance in low-altitude airspace with geofencing: Constrained stochastic linear hybrid systems approach. In Proceedings of the 2018 AIAA Information Systems-AIAA Infotech @ Aerospace, Kissimmee, FL, USA, 8–12 January 2018. [Google Scholar] [CrossRef]
  17. Shan, D.; Wang, D.; He, D.; Zhang, P. Position Sensorless Vector Control System for Lawnmower Permanent Magnet Synchronous Motor Based on Extended Kalman Filter. Energies 2024, 17, 1230. [Google Scholar] [CrossRef]
  18. Khan, A.B.; Akram, A.S.; Choi, W. State of Charge Estimation of Flooded Lead Acid Battery Using Adaptive Unscented Kalman Filter. Energies 2024, 17, 1275. [Google Scholar] [CrossRef]
  19. Lv, Y.; Liu, S.; Gao, Y.; Dai, J.; Ren, Z.; Liu, Y. An Ultra-Wideband Indoor Localization Algorithm with Improved Cubature Kalman Filtering Based on Sigmoid Function. Appl. Sci. 2024, 14, 2239. [Google Scholar] [CrossRef]
  20. Tivay, A.; Hahn, J. A Population-Informed Particle Filter for Robust Physiological Monitoring Using Low-Information Time-Series Measurements. IEEE Trans. Biomed. Eng. 2023, 70, 2298–2309. [Google Scholar] [CrossRef] [PubMed]
  21. Zare, H.; Weber, T.K.D.; Ingwersen, J.; Nowak, W.; Gayler, S.; Streck, T. Combining Crop Modeling with Remote Sensing Data Using a Particle Filtering Technique to Produce Real-Time Forecasts of Winter Wheat Yields under Uncertain Boundary Conditions. Remote Sens. 2024, 14, 1360. [Google Scholar] [CrossRef]
  22. Kanthalakshmi, S.; Raghappriya, M. Active fault diagnosis of 2 DoF helicopter using particle filter-based log-likelihood ratio. Int. J. Control 2022, 95, 3148–3165. [Google Scholar] [CrossRef]
  23. Peng, S.; Zhang, A.; Liu, D.; Cheng, M.; Kan, J.; Pecht, M. State-of-Charge Estimation of Lithium-Ion Batteries Based on Dual-Coefficient Tracking Improved Square-Root Unscented Kalman Filter. Batteries 2023, 9, 392. [Google Scholar] [CrossRef]
  24. Janjanam, L.; Saha, S.K.; Kar, R.; Mandal, D. Hammerstein-Wiener nonlinear system identification by using honey badger algorithm hybridized Sage-Husa adaptive Kalman filter with real-time applications. Int. J. Electron. Commun. 2022, 151, 154218. [Google Scholar] [CrossRef]
  25. Seah, C.E.; Hwang, I. State Estimation for Stochastic Linear Hybrid Systems with Continuous-State-Dependent Transitions: An IMM Approach. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 376–392. [Google Scholar] [CrossRef]
  26. Fazilova, D. Uzbekistan’s coordinate system transformation from CS42 to WGS84 using distortion grid model. Geod. Geodyn. 2022, 13, 24–30. [Google Scholar] [CrossRef]
  27. Wang, P.; Wang, M.; Yang, Y. Particle Filter with improved MSV Resample for Manoeuvring Target Tracking. In Proceedings of the 2020 Chinese Automation Congress (CAC), Shanghai, China, 6–8 November 2020; pp. 6522–6527. [Google Scholar] [CrossRef]
  28. Petrera, M.; Suris, Y.B. Spherical geometry and integrable systems. Geom. Dedicata 2014, 169, 93–98. [Google Scholar] [CrossRef]
  29. Gao, Z.; Zhang, Y.; Ochieng, W.Y.; Xiang, Z. Adaptive Air-Data Smoothing Estimation with Customized Wind Model Based on Flight Data. J. Aerosp. Eng. 2022, 35, 04022054. [Google Scholar] [CrossRef]
  30. Gao, Z.; Wang, H.; Xiang, Z. Flight Data-Based Wind Disturbance and Air Data Estimation. Atmosphere 2021, 12, 470. [Google Scholar] [CrossRef]
Figure 1. Structure of the nonlinear SHE algorithm.
Figure 1. Structure of the nonlinear SHE algorithm.
Aerospace 11 00413 g001
Figure 2. Great Circle projection.
Figure 2. Great Circle projection.
Aerospace 11 00413 g002
Figure 3. Base turn chart.
Figure 3. Base turn chart.
Aerospace 11 00413 g003
Figure 4. Error comparison of position and azimuth accuracy.
Figure 4. Error comparison of position and azimuth accuracy.
Aerospace 11 00413 g004
Figure 5. Comparison of the accuracy of mode transitions.
Figure 5. Comparison of the accuracy of mode transitions.
Aerospace 11 00413 g005
Figure 6. Mode probability estimation by ASR-UPF with 300 particles.
Figure 6. Mode probability estimation by ASR-UPF with 300 particles.
Aerospace 11 00413 g006
Figure 7. Mode probability estimation by SR-UKF.
Figure 7. Mode probability estimation by SR-UKF.
Aerospace 11 00413 g007
Figure 8. Estimation with respect to real flight data.
Figure 8. Estimation with respect to real flight data.
Aerospace 11 00413 g008
Table 1. Means and variances of parameters.
Table 1. Means and variances of parameters.
ParametersWPT1WPT2Level Flight
μ H / μ η σ H / σ η μ H / μ η σ H / σ η μ H σ H
H m 2241.7673.0462086.0722.5361786.7581.653
η 108.1271.862274.8451.387--
Table 2. The probabilities of mode transition.
Table 2. The probabilities of mode transition.
Modes m 1 m 2 m 3 m 4
m 1 1000
m 2 0000
m 3 γ 31 0 1 γ 34 γ 31 γ 34
m 4 00 γ 43 1 γ 43
Table 3. RMSE of 50 simulations by SR-UKF and ASRUPF.
Table 3. RMSE of 50 simulations by SR-UKF and ASRUPF.
AlgorithmRMSE
x n (m) y n (m) z n (m) η (°)
SR-UKF43.58629.59463.8604.862
ASR-UPF with 100 particles20.65821.39525.6744.206
ASR-UPF with 300 particles14.37712.59320.8053.551
ASR-UPF with 500 particles14.47612.20319.2313.704
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

Zhang, Y.; Gao, Z.; Qi, K.; Li, J. Refined Aircraft Positioning Based on Stochastic Hybrid Estimation with Adaptive Square-Root Unscented Particle Filtering. Aerospace 2024, 11, 413. https://doi.org/10.3390/aerospace11050413

AMA Style

Zhang Y, Gao Z, Qi K, Li J. Refined Aircraft Positioning Based on Stochastic Hybrid Estimation with Adaptive Square-Root Unscented Particle Filtering. Aerospace. 2024; 11(5):413. https://doi.org/10.3390/aerospace11050413

Chicago/Turabian Style

Zhang, Yangyang, Zhenxing Gao, Kai Qi, and Jiawei Li. 2024. "Refined Aircraft Positioning Based on Stochastic Hybrid Estimation with Adaptive Square-Root Unscented Particle Filtering" Aerospace 11, no. 5: 413. https://doi.org/10.3390/aerospace11050413

APA Style

Zhang, Y., Gao, Z., Qi, K., & Li, J. (2024). Refined Aircraft Positioning Based on Stochastic Hybrid Estimation with Adaptive Square-Root Unscented Particle Filtering. Aerospace, 11(5), 413. https://doi.org/10.3390/aerospace11050413

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