Next Article in Journal
Understanding Customers’ Insights Using Attribution Theory: A Text Mining and Rule-Based Machine Learning Two-Step Multifaceted Method
Next Article in Special Issue
A Single-Launch Deployment Strategy for Lunar Constellations
Previous Article in Journal
Spherical Atomic Norm-Inspired Approach for Direction-of-Arrival Estimation of EM Waves Impinging on Spherical Antenna Array with Undefined Mutual Coupling
Previous Article in Special Issue
Analysis of Resonance Transition Periodic Orbits in the Circular Restricted Three-Body Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Gaussian Mixture Model for Uncertainty Propagation Using Virtual Sample Generation

1
School of Astronautics, Harbin Institute of Technology, Harbin 150001, China
2
Department of Mathematics and Theories, Peng Cheng Laboratory, No. 2, Xingke 1st Street, Nanshan, Shenzhen 518000, China
3
Deep Space Exploration Labortory, Beijing 100089, China
4
School of Aerospace Engineering, Beijing Institute of Technology, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(5), 3069; https://doi.org/10.3390/app13053069
Submission received: 18 January 2023 / Revised: 17 February 2023 / Accepted: 20 February 2023 / Published: 27 February 2023
(This article belongs to the Special Issue Astrodynamics and Celestial Mechanics)

Abstract

:
Orbit uncertainty propagation plays an important role in the analysis of a space mission. The accuracy and computation expense are two critical essences of uncertainty propagation. Repeated evaluations of the objective model are required to improve the preciseness of prediction, especially for long-term propagation. To balance the computational complexity and accuracy, an adaptive Gaussian mixture model using virtual sample generation (AGMM-VSG) is proposed. First, an unscented transformation and Cubature rule (UT-CR) based splitting method is employed to adaptive update the weights of Gaussian components for nonlinear dynamics. The Gaussian mixture model (GMM) approximation is applied to better approximate the original probability density function. Second, instead of the pure expensive evaluations by conventional GMM methods, virtual samples are generated using a new active-sampling-based Kriging (AS-KRG) method to improve the propagation efficiency. Three cases of uncertain orbital dynamical systems are used to verify the accuracy and efficiency of the proposed manuscript. The likelihood agreement measure (LAM) criterion and the number of expense evaluations prove the performance.

1. Introduction

The analysis of orbit uncertainty propagation is a key step in space mission design [1,2]. Because of the existence of uncertainty, e.g., injection errors and measurement errors, the spacecraft will deviate from their nominal orbits [3,4], which may change the orbital stability and impact the performance of exploration mission [5,6]. Therefore, it is necessary to investigate the evolution of orbital uncertainty under multi-source errors [7], which provides feedback coupling inputs for control and navigation subsystem design [8,9].
Generally, uncertainty propagation is used to characterize the spacecraft’s state moment (usually represented by mean and covariance matrix) or its probability density function (PDF) [10,11]. Theoretically, the uncertainty description of the orbital state can only be accurately obtained by solving Fokker-Plank Equation (FPE) [12,13] or carrying out Monte Carlo (MC) simulations [9]. Though the FPE provides accurate solutions for the propagation of the PDF, it is usually difficult to solve the equation as it is a high-dimensionally defined and nonlinear perturbed Keplerian-driven partial differential equation. The MC simulation is always considered a reliable method whose result approaches the true probability distribution when the number of MC samples tends to be infinity [14]. As a result, either FPE or MC simulations suffer from drawbacks related to a heavy computation burden [9]. To improve efficiency, some approximations are made, which can be categorized into gradient-enhanced methods and gradient-free methods [15,16].
Gradient-enhanced methods, such as state transfer matrix (STM) [17,18] and state transition tensors (STT) [19,20,21], analyze the orbital uncertainty by simplifying the relative nonlinear motion with the Taylor series terms [22,23]. For gradient-enhanced methods, the gradients or higher-order derivatives are incorporated for approximation [20], and analytical or semi-analytical solutions are obtained [24,25]. However, gradient-enhanced methods are only suitable for cases whose gradient information is easy to obtain. For some of the practical problems, the dynamic motion is complex, discontinuous, and nondifferentiable, and the derivatives are unavailable [16]. The gradient-free methods (also called non-intrusive methods [26]) allow for the usage of already existing orbit propagation tools as the black-box and avoid any simplifications on the fully perturbed orbital dynamics [27,28], making them conveniently extend into orbit uncertainty propagation problems [29,30]. Recently, some sample-based methods, for example, unscented transformation (UT) [31,32], Conjugate Unscented Transform (CUT) [33], and Cubature rules (CRs) [34,35], have been well investigated. These methods approximate the probability distribution at a future time by nonlinearly integrating a few sigma samples, which are deterministically chosen according to the given initial distribution [9]. The drawbacks of the above sample-based methods are that they can only deliver a second-order approximation of the distribution, which is inadequate for long-term propagation whose ultimate uncertainty distribution inevitably turns out to be highly non-Gaussian [36,37].
In order to fully characterize the non-Gaussian orbital uncertainty, the Gaussian mixture model (GMM) is commonly employed [38,39]. The basic idea of GMM is that an arbitrary PDF can be approximated, as closely as desired, by a finite sum of weighted Gaussian components (also called Gaussian PDFs). For each Gaussian component of GMM, only the mean and covariance need to be propagated because a Gaussian PDF can be completely determined with its first two moments [40,41]. Therefore, the above methods, such as QRs, UT, CUT, and CRs, can be used to propagate each Gaussian component. The GMM can better represent non-Gaussian PDFs and results in a more accurate prediction [42]. However, the GMM suffers from a heavy computation burden. Taking the CMM combined with a fifth-order CR as an example, in this case, 2190–14,600 sigma-samples are to be integrated, which is hundreds of times larger than the conventional single-Gaussian propagation using fifth-order CR [16]. Thus, it is necessary and promising to investigate approaches to enhance the efficiency of GMM.
In recent, the use of virtual sample generation (VSG) has been well discussed. The VSG is a kind of technique that aims at providing artificial samples to enlarge the current data set, which was first introduced in the field of machine learning [43,44]. Apart from the field of machine learning, the VSG technique has also been successfully applicated to engineering optimization problems. Tang et al. proposed a VSG-based optimization method to reduce the required expensive evaluation [45]. Motivated by the concept of the VSG, it is possible to improve the efficiency of the GMM by incorporating visual samples in propagation.
In this paper, an adaptive Gaussian mixture model using virtual sample generation (AGMM-VSG) is proposed. The virtual sample generation (VSG) technique [45,46] is employed to reduce the required computation for propagating each Gaussian component. Furthermore, a novel Gaussian component splitting strategy based on the virtual sample is employed for efficiently updating weights.
The paper is organized as follows. In Section 2, the orbital uncertainty propagation problem with a single Gaussian model and a GMM is described. Section 3 presents details of the proposed AGMM-VSG, followed by detailed discussions of the adaptively splitting strategy for Gaussian components and the VSG-based propagation method. The proposed AGMM-VSG is applied to three cases of uncertainty propagation for uncertain orbital dynamical systems in Section 4. Finally, Section 5 gives the conclusions.

2. Problem Statement and Methodology

2.1. Problem Statement

Consider a nonlinear dynamical system governed by the following differential equation [42]:
x ˙ ( t ) = f ( x ( t ) , t ) x t 0 = x 0
where x ( t ) R n is the state of the system, f ( x ( t ) , t ) R n represents the sufficiently differentiable nonlinear dynamics of the system, and x 0 is the initial condition. For black-box systems, the differential equation (or gradient information) is unavailable; thus the corresponding discrete-time nonlinear system of Equation (1) is as follows:
x t k + 1 = F x t k , d t k x t 0 = x 0
In practice, it is usually assumed that the unknown uncertainties as Gaussian distribution (or Gaussian PDF), which can be represented by [10]:
p ( x ; m , P ) = N ( x ; m , P ) = 1 | 2 π P | exp 1 2 ( x m ) T P 1 ( x m )
where N ( x ; m , P ) denotes the normal distribution with x , m and P being the input variables, mean and covariance, respectively.

2.2. Single Gaussian Model-Based Uncertainty Propagation

Generally, uncertainty propagation is used to determine the specific statistical moments of dynamics (1) or (2), with a given initial probability distribution and propagation period. The state moment, i.e., mean, and covariance matrix or probability density function (PDF) are specific quantitative measures of statistical moments. The mean and covariance matrix could be defined as follows with given random vector x R n and PDF p ( t ) :
m = E [ x ] = ξ p ( ξ ) d ξ P = E ( x m ) ( x m ) T = ( ξ m ) ( ξ m ) T p ( ξ ) d ξ
where ξ is a random variable.
The Gaussian weighted integral (4) can be approximated efficiently with different quadrature rules, such as the unscented transformation (UT), Conjugate Unscented Transform (CUT), and Cubature rules (CRs). In this paper, the fifth-order CR is utilized for propagating the single Gaussian distribution. The fifth-order CR is detailed in Ref. [47], which will not be repeated here.

2.3. Gaussian Mixture Model Approach

GMM is a direct extension of the Gaussian PDF, which is given by a sum of weighted Gaussian PDFs:
p ( x ) p ˜ ( x ) = i = 1 N α i p g x ; m i , P i
where N represents the number of components of the GMM, α i , m i , and P i denote the weights, mean, and the covariances associated with component i . Each Gaussian PDF p g x ; m i , P i for x is defined per Equation (3). With given GMM per Equation (5), the final mean and covariance matrix is computed by:
m = i = 1 N α i m i P = i = 1 N α i P i + m i m i T m m T

3. Methodology of Adaptive Gaussian Mixture Model Using Virtual Sample Generation

3.1. Overall Procedure of AGMM-VSG

Before introducing the methodology of the proposed AGMM-VSG, some concepts in AGMM-VSG are explained as follows:
Definition 1.
(point, response, and sample): As is shown in Table 1, the point is a kind of pure coordinate. In AGMM-VSG, point refers to any vector of state x t k at time t k . The response is taken as the output of the black-box model (2), i.e., the state vector x t k + 1 , which corresponds to the point x t k . A sample is a combination of a point and its response.
Definition 2.
(real sample and virtual sample): The real sample represents a combination of a point and its real response generated by invoking expensive black-box simulations while the virtual sample refers to the combination of a point and its virtual response generated without calling expensive simulations. The virtual responses are only generated using the active KRG model. Note that the number of real samples indicates the computational expense of the uncertainty propagation. Moreover, in this paper, the virtual response is generated by the active sampling-based KRG model, and the expense to generate virtual samples is negligible.
Definition 3.
(cheap point and superior cheap point): Cheap point is a point generated according to the given PDF. Note that the actual response of a cheap point is not calculated, and therefore, the computational cost is negligible. The cheap points are drawn as a candidate pool for constructing the AS-KRG model. The superior cheap point has a promising performance for improving the global accuracy of the KRG model, which is classified on basis of the variance predicted using the KRG model. The process of using both the cheap point and superior cheap point will be detailed in Section 3.3.
The flowchart of AGMM-VSG is shown in Figure 1, and the procedure is explained in detail as follows:
Step 1: The parameters of the problem are determined, including the black-box discrete dynamics, the initial PDF (or initial mean and covariance matrix), the propagating period t f , and discrete-time d t . Other necessary items, such as the maximal number of samples for constructing the KRG model and splitting library, will be discussed in the next section. The current time epoch is set to be zero, i.e., t = 0 .
Step 2: In the first iteration, the single Gaussian model is used, with its mean and covariance matrix equaling the initial PDF. The number of Gaussian mixture components is one.
Step 3: The KRG model is first constructed using initial real samples and then is gradually refined with newly added real samples. The constructed KRG model is utilized for generating virtual samples.
Step 4: Propagate each component using fifth-order CR. Note that only virtual samples are used when propagating; thus, the expense of this step is negligible. Renew the current time epoch, i.e., t = t + d t .
Step 5: For each component, determine whether the splitting process should be executed. If it should, then apply a univariate splitting library in a specified direction, i.e., the principal (eigenvector) direction of the covariance matrix.
Step 6: If t < t f , jumps to Step 3. Otherwise, it outputs the results.
The algorithms of Gaussian components splitting and virtual sample generation-based propagation are detailed in Section 3.2 and Section 3.3, respectively.
Figure 1. Flowchart of AGMM-VGS.
Figure 1. Flowchart of AGMM-VGS.
Applsci 13 03069 g001

3.2. Gaussian Component Nonlinearity Detecting

It is noted that only the linear propagation can maintain the Gaussian property, and any deviation in the nonlinear dynamics from the linearized solution may indicate nonlinear effects. This deviation can be handled by specifying a threshold and monitoring the nonlinearity [48]. That is to say, when the nonlinearity causes a more significant non-Gaussian property, a component of the GMM can be split into finer Gaussian components. A split GMM component is considered “finer” because Gaussians with smaller covariances would remain more Gaussian compared with ones with larger covariances, which will make the uncertainty propagation more representative. Thus, an adaptive GMM with splitting technologies is necessary, especially for long-term propagation. In this section, a novel unscented transformation and Cubature rule-based (UT-CR) method for monitoring (or detecting) the nonlinearity is introduced, and the Gaussian component splitting will be detailed in the next section.
Usually, the term entropy is used for detecting nonlinearity [42]. For a Gaussian distribution with the determinant of a scaled form of the covariance matrix, the corresponding differential entropy (DE), in terms of the logarithm, is given as:
H ( x ) = 1 2 log | 2 π e P |
where the temporal derivative of the DE is given by:
H ˙ ( x ) = 1 2 trace P 1 P ˙
where P ˙ is the temporal derivative of the covariance matrix P .
In case of a linearized representation, the time rate of the DE is written as:
H ˙ ( x ) = 1 2 trace f x
where f x = f / x represents the linearized dynamics Jacobian.
It is worth noting that for linearized dynamics, H ˙ ( x ) = 0 , meaning that the differential entropy is constant for the linearized system, making the DE unable to be used for splitting the Gaussian components for the linearized system. In this paper, a UT-CR-based splitting strategy is utilized to detect nonlinearity.
The novel dynamic nonlinearity index measures the departure from linearity as [49]:
v t , t 0 = Φ C R 5 t , t 0 Φ U T t , t 0 Φ U T t , t 0
where v t , t 0 is the dynamic nonlinearity index and · represents the norm. The Φ C R 5 t , t 0 and Φ U T t , t 0 are the STM obtained by fifth-order CR (CR5) and UT, respectively. It is noted that the UT could capture the solutions accurately to the third order while the CR5 could capture accurately to the fifth order; thus, the nonlinearity with fourth and fifth orders are introduced with the index v t , t 0 .
The STM Φ C R 5 t , t 0 and Φ U T t , t 0 are obtained as follows. First, the posterior covariances (transformed distribution) P C R 5 t , t 0 and P U T t , t 0 are calculated using fifth-order CR and UT, given the same prior probability P t 0 = P 0 . Consider the following first-order approximation of the covariance of the transformed distribution [50]:
P C R 5 ( t ) = Φ C R 5 t , t 0 P 0 Φ C R 5 T t , t 0 P U T ( t ) = Φ U T t , t 0 P 0 Φ U T T t , t 0
If the Cholesky factorization P 0 = L 0 L 0 T , P C R 5 ( t ) = L C R 5 ( t ) L C R 5 T ( t ) and P U T ( t ) = L U T ( t ) L U T T ( t ) , where L 0 , L C R 5 ( t ) and L U T ( t ) are the corresponding unique lower triangular “square root” matrixes. Therefore, Equation (11) is rewritten as:
L C R 5 ( t ) L C R 5 T ( t ) = Φ C R 5 t , t 0 L 0 L 0 T Φ C R 5 T t , t 0 L U T ( t ) L U T T ( t ) = Φ U T t , t 0 L 0 L 0 T Φ U T T t , t 0
A transformation equation can be written as:
L C R 5 ( t ) = Φ C R 5 t , t 0 L 0 L U T ( t ) = Φ U T t , t 0 L 0
The solution for Equation (13) is as:
Φ C R 5 t , t 0 = L C R 5 ( t ) L 0 1 Φ U T t , t 0 = L U T ( t ) L 0 1
Substitute Equation (14) into Equation (10), we arrive at the following expression of dynamic nonlinearity index:
v t , t 0 = L C R 5 ( t ) L 0 1 L U T ( t ) L 0 1 L U T ( t ) L 0 1 = L C R 5 ( t ) L U T ( t ) L 0 1 L U T ( t ) L 0 1 = L C R 5 ( t ) L U T ( t ) / L 0 L U T ( t ) / L 0 = L C R 5 ( t ) L U T ( t ) L U T ( t )
The fewer the values the nonlinearity index indicates, the less significant the nonlinearity. Note that the virtual samples, which will be discussed in the next subsection, can be used for propagating the P C R 5 ( t ) and P U T ( t ) to reduce computation cost.

3.3. Gaussian Component Splitting

Once the nonlinearity is detected, then the splitting is to be executed. Current splitting algorithms can be summarized into two categories: (1) Split a Gaussian PDF in all directions or (2) Split a Gaussian PDF by applying a univariate splitting library in a specified direction [48]. Solving an L2-minimization optimization problem is usually required when using the first kind of method, which incurs additional computational consumption. In this paper, splitting in one direction is employed. The multivariate Gaussian PDF is split by applying a univariate splitting library in the direction of maximum eigenvalue, which is given per the following steps:
First a standard normal distribution p ( x ; 0 , 1 ) is split into N homoscedastic Gaussian distributed components as N p ( x ; 0 , 1 ) i = 1 N α ˜ i p g x ; m ˜ i , σ ˜ i 2 . This step is independent from the concerned problem; therefore, parameters α ˜ i , m ˜ i , σ ˜ i can be computed offline and saved as a splitting library. The computation of the splitting parameters is well-provided by [42], and one case with is illustrated in Table 2.
Then, consider the spectral factorization of the covariance matrix P i as:
P i = V Λ i V T
where Λ i = diag λ 1 , , λ k , , λ n is the eigenvalues matrix with λ k being the k-th largest eigenvalue of the covariance matrix P i . The i-th component is approximated by p g x ; m i , P i j = 1 N α ˜ j p g , j x ; m i , j , P i , j with mean and covariance matrix given as:
m i , j = m i + λ k m ˜ j v k
P i , j = V Λ ˜ i V T
where v k is the k-th eigenvector of P i and Λ i ˜ = diag λ 1 , , λ k , , λ n is the set of eigenvalues of the i-th new component.

3.4. Virtual Sample Generation-Based Propagation Strategy

It is noted that for a black-box model, only sample-based gradient-free methods can be used for single Gaussian component propagation, and fifth-order CR is selected in this paper, as is shown in Equation (19).
λ F ( x ) p ( x ; m , P ) d x = F ( Ax + m ) p ( x ; m , P ) d x i = 1 N p W i F γ i
where γ i i = 1 , , N p are Cubature points according to the fifth-order CR.
Although the Gaussian component splitting technology could reduce unnecessary components, the computation burden, if real samples are involved when propagating, is still intractable. However, the computation burden can be released by using virtual samples with the information hidden in the raw samples (some already obtained real samples):
F ( x ) p ( x ; m , P ) d x i = 1 N p W i F γ i i = 1 N p W i F ˜ γ i
where F ˜ γ i is the response of the virtual samples.
Generally, virtual samples can be constructed using prior knowledge, the addition of noise, and existing datasets. In this paper, an active sampling-based KRG model based on existing datasets is proposed for generating virtual samples for the black-box uncertainty model.
KRG [16,51] is an interpolation metamodeling technique as shown in Equation (20).
f K R G ( x ) = G ( x ) + Z ( x )
where f K R G ( x ) is the output function of the KRG model, G ( x ) captures the major trend of input–output mapping, and Z ( x ) is a zero-mean, unit-variance stationary Gaussian process that determines the approximation ability of the KRG model. The covariance matrix of Z ( x ) is given as follows:
Cov Z x i , Z x j = σ 2 R R x i , x j
where σ 2 is the process variance of Z ( x ) , R is the correlation matrix, and R ( · , · ) is a Gaussian correlation function that controls the smoothness of the KRG model. Note that in addition to the output predicted using the KRG model, the variance, which can be taken as the degree of confidence of the prediction, can also be obtained.
The performance of the KRG model will influence the accuracy of the virtual samples, which, in turn, will influence the accuracy of propagation. Although investigations have been conducted to improve the global accuracy of the KRG model, it is still difficult to use a static metamodel to accurately approximate a sophisticated problem (e.g., high-dimensional nonlinear problem) with sparsely scattered samples. However, to improve the convergence and efficiency, the adaptive sampling is actually realized by the management of the online KRG models [46]. In this paper, an adaptive sampling called the “active sampling method” is adopted, with which the superior cheap points are adaptively selected.
The least-confidence strategy is adopted in this paper to perform the active sampling. The superior sample is selected according to the following metric:
x * = arg max V a r R B F ( x )
where V a r R B F ( x ) is the variance of the output prediction f K R G ( x ) . Comparison studies in Ref. [46] have proven that the KRG generally has better approximation accuracy than other technologies, for example, RBF and Gaussian Process (GP) [52], in the condition of relatively sparse samples and lack of prior knowledge. In contrast, the construction expense of KRG is more expensive than that of the Radial basis function (RBF) [53]. Therefore, to further reduce the expense, in this paper, up to n a d d samples, with the larger variances, are added into the sample set at each iteration, for constructing the KRG model, where n a d d is the user-defined parameter. In this paper, the KRG model is constructed using MATLAB toolbox SURROGATES.
The procedure of this strategy is exhibited in detail as follows:
Step 1: The parameters are determined, including the total number of cheap points n c , the number of superior cheap points added at each iteration n a d d , the maximal number of real samples n r , and the number of initial real samples n i .
Step 2: Generate n c points according to the current PDF p ( x , t ) as the candidate points set; denote as S c 0 = x 1 , x 2 , , x n c . Select n i points from the candidate points set S c 0 , note as S s = x i 1 , x i 2 , , x i n i . Evaluate the points in S s using black-box dynamical model (2), and then responses are obtained. All of these n i points and their corresponding responses make up the sample set. Construct an RBF model based on the existing samples. Renew the candidate points set as S c 1 = x x S c 0 and x S s . The iteration counter k is set to be one, i.e., k 1 . The current number of real samples n s is set to be n s = n i .
Step 3: Evaluate all the points in candidate points set S c k using the KRG model. Sort these points according to their predicted variances. The top n a d d superior points are selected, noted as S a d d k . Evaluate these points and then add them to the sample set. Construct the RBF model. Remove the selected superior points from the candidate points set as S c k + 1 = x x S c k and x S s k , k k + 1 , n s n s + n a d d .
Step 4: If n s < n r , return to Step 3; elsewise jump to Step 5.
Step 5: Use the virtual sample evaluated by the KRG model to propagate.
To further demonstrate the process of active sampling based KRG model, a one-dimensional numerical benchmark is tested. The benchmark describes the mapping relation between the final position error along the x-axis and the initial velocity error along the y-axis, which is given as:
g Δ v y t 0 = Δ r x t f
The spacecraft moves on the equatorial plane and the dynamics considered the Earth’s gravity and atmospheric drag, which is governed by [42]:
x = r T , v T T = r x , r y , v x , v y T x ˙ = f ( x ) = v x , v y , μ r x r 3 ρ ( h ) β v rel v rel , x / 2 , μ r y r 3 ρ ( h ) β v rel v rel , y / 2 T
where r = r x 2 + r y 2 , V r e l , x = V x ω r y , V r e l , y = V y ω r x , and ω is the angular velocity of the Earth. ρ ( h ) is the atmospheric density as a function of the altitude of the object. Additionally, it is assumed that the atmospheric density is described by an exponential atmosphere model via:
ρ ( h ) = ρ 0 exp h h 0 / h s
with ρ 0 = 3.614 × 10 13 kg / m 3 , h 0 = 700 km and h s = 88.667 km . β is the ballistic coefficient of the object, and it is taken to be 1.4 in this paper.
The elements of the nominal orbit are illustrated in Table 3, and the propagation time is set to be one nominal orbit period, which is about 5331.288 s. The points are set to be evenly distributed in range. The cheap points, real samples, and newly added real samples (superior cheap points) are illustrated in Figure 2. As is shown in Figure 2, the newly added real samples (blue triangle) are selected where the KRG model is the most unconfident. Note that the global accuracy of the KRG model is improved when samples are added sequentially.

4. Numerical Simulations

In this section, the developed AGMM-VSG method for uncertainty propagation is applicated to the propagation of uncertainty for two test problems, i.e., a high-Earth orbit case and a low-Earth orbit case. The performance of the proposed AGMM-VSG method is discussed by comparing it with the conventional fifth-order Cubature rule and adaptive Gaussian mixture model method.

4.1. High-Earth Orbit Case

First, consider a spacecraft orbiting on a high-Earth orbit (HEO). Assume that the motion of the spacecraft is confined to the equatorial plane, which allows the position r to be described by two scalar values r x and r y , and the velocity v to be described by two scalar values v x and v y . Thus, the state vector and equations of motion that describe the nonlinear dynamical system are written as follows [42]:
x = r T , v T T = r x , r y , v x , v y T x ˙ = f ( x ) = v x , v y , μ r x r 3 , μ r y r 3 T
where μ = 398 , 600 is the gravitational constant of the central body (i.e., Earth in this paper).
Take notice that although dynamics (27) are of the analytical form, it is still considered a black-box model by invoking the kinetic integral F ( t ) = 0 t f ( x ) d x .
The elements of the nominal orbit are listed in Table 4, given by a semimajor axis of 35,000 km and an eccentricity of 0.2.
The initial state distribution is taken to be Gaussian, with the initial mean represented by the nominal orbit elements. Additionally, the initial covariance is set to be diagonal with standard deviations in the position elements of 1 km and the velocity elements of 1 m/s:
P 0 = diag 1 km 2 , 1 km 2 , 10 6 km 2 / s 2 , 10 6 km 2 / s 2
To characterize the performance of the proposed AGMM-VSG method, both the AGMM-VSG and fifth-order CR methods are applied to the problem of predicting the evolution of uncertainty. The propagation time is one orbit period (around 65,164.833 s). The values of parameters to be used in AGMM-VSG are set as t f = 65,164.833 s, d t = t f / 100 . The criterion for splitting the Gaussian components are v > 10 4 , where v is the nonlinearity index (Equation (15)) of the corresponding component. Parameters for active sampling are given as: n c = 200 , n r = 100 , n i = 50 , and n a d d = 10 .
In addition, a Monte Carlo simulation using 1000 samples drawn from the initial distribution is performed. The position and velocity marginal PDF contours using AGMM-VSG and fifth-order CR are illustrated in Figure 3 and Figure 4, respectively. The outputs of the Monte Carlo simulations are shown represented with a grey x. Depart from the scatter plot, the PDFs of each state variable are also illustrated. The PDFs of each state variable of AGMM-VSG and fifth-order CR are obtained by calculating the integral of the multivariate PDFs, while PDFs of each state variable of 1000 Monte Carlo samples are obtained using MATLAB function ksdensity. It can be seen from Figure 3 and Figure 4 that, compared with fifth-order CR, the proposed method is better at representing the curvature exhibited by the Monte Carlo samples. For AGMM-VSG, only 65 Gaussian components and 100 real samples are requested. The computation cost of AGMM-VSG is a little higher than that of fifth-order CR, which uses a total of 33 evaluations ( 2 n 2 + 1 Cubature-points and n = 4 ), while significantly outperforming fifth-order CR in terms of accuracy. Moreover, if use real samples for propagating, a total of 2145 real samples (65 components and 33 evaluations for each component) are needed, which indicates the signification of the utilization of virtual samples for reducing computational expense.
To further provide a relative comparison of the performance of the two solutions obtained (from the AGMM-VSG and fifth-order CR), a criterion called the likelihood agreement measure (LAM) is utilized [42]. The LAM is used to evaluate the likelihood agreement between two distributions. Given two PDFs, noted as p ( x ) and q ( x ) , the LAM is defined as:
LAM ( p , q ) = p ( x ) q ( x ) d x
where Equation (29) represents the density of overlap between the two PDFs. The greater agreement for one PDF with another PDF, the larger their densities.
However, the analytical distribution is intractable to obtain. In a practical problem, the true distribution is approximated using the Monte Carlo samples, and the LAM is rewritten as
LAM ( p , q ) = i = 1 N j = 1 O α i γ j p g x j ; m i , P i
where Q is the total number of samples within the Monte Carlo simulation.
The accuracy of different methods can be compared using the same set of Monte Carlo samples. The LAM of the proposed AGMM-VSG and fifth-order CR are drawn in Figure 5. Moreover, the LAM of the AGMM method, where only real samples are employed, is also tested as the standard results. To be clear, the likelihood measure of the AGMM-VSG is normalized by the value for the AGMM method. It can be observed that the LAM of AGMM-VSG is close to that of AGMM, while clearly outperforming the fifth-order CR method.

4.2. Low-Earth Orbit Case

The dynamic system for the low-Earth-orbit propagation is given per Equation (25). The elements of the nominal orbit are illustrated in Table 3, and the propagation time is set to be one nominal orbit period. The initial covariance is taken to be diagonal with standard deviations of 1.3 km in position r x , 0.5 km in position r y , 2.5 m/s in velocity v x , and 5 m/s in velocity v y , that is:
P 0 = diag 0 . 5 2 km 2 , 2 . 5 2 km 2 , 2 . 5 2 × 10 6 km 2 / s 2 , 5 2 × 10 6 km 2 / s 2
Similarly to the previous example, the performance of AGMM-VSG is verified by comparing it with fifth-order CR. The parameters involved in AGMM-VSG are set to be the same as the previous HEO case, and the results are shown in Figure 6 and Figure 7. From Figure 6 and Figure 7, it is seen that both implementations of the AGMM-VSG approach enable a better representation of the curvature exhibited by the Monte Carlo samples.
Additionally, 25 components are used to represent the final PDF, and only 100 evaluations are requested. The LAMs of AGMM, AGMM-VSG, and fifth-order CR are also normalized to provide a clear relative measure. The results are illustrated in Figure 8. It is easy to observe that the fifth-order CR is outperformed by both AGMM and AGMM-VSG.
With the validation of the two numerical simulations validation, the advantages of the AGMM-VSG method in accuracy and efficiency can be recapitulated. The proposed AGMM-VSG could represent the distribution as well as AGMM, and is better than the fifth-order CR method, while it is much more efficient than the AGMM.

4.3. TianQin Orbit Uncertainty Propagation

To demonstrate the effectiveness and practicality of the proposed AGMM-VSG in solving uncertainty propagation of computationally expensive, black-box systems, as well as the application in real-world space-based GW observatories, a case of TianQin Orbit (TQO) is investigated.
The state of the detector is defined as x = r x , r y , r z , v x , v y , v z T . The high-fidelity dynamics model includes an 8 × 8 spherical-harmonic Earth gravity field (JGM-3), the point-mass gravity field from the Moon, and the Sun (the ephemeris JPL405). The spacecraft of TianQin each have a disturbance reduction system (DRS) needed to reduce the effects of the non-gravitational forces [54]; thus, the solar radiation pressure and the atmospheric drag are neglected.
The elements of the nominal orbit of one detector of TianQin are shown in Table 5. The orbit is propagated in one period, which is about 314,710.3177 s.
Both the AGMM-VSG and the fifth-order CR are applicated to the above TQO problem of predicting the evolution of uncertainty. In employing the AGMM-VSG, the criterion for splitting the Gaussian components is set to be v > 10 3 , and the number of evaluations for constructing the surrogates in each iteration is set to 300. The final state uncertainty and their corresponding means and 3 σ ellipsoids, using proposed AGMM-VSG and MC, are projected in Figure 9. It is witnessed that the accuracy of AGMM-VSG and MC are almost the same when considering the first two moments. Moreover, the projection of the predicted PDF from the AGMM-VSG, together with the 1000 Monte Carlo samples, are shown in Figure 10. It is clearly observed that both the implementations of the position and velocity of the AGMM-VSG approach characterize the non-linearity of the MC samples-approximated distribution.

5. Conclusions

In this paper, to reduce the computational burden while maintaining prediction accuracy in orbital uncertainty propagation, an adaptive Gaussian mixture model using virtual sample generation—noted as AGMM-VSG—is proposed. The splitting algorithm is applied to the dominant eigenvalue of the covariance matrix of the Gaussian components. An active sampling-based KRG algorithm is employed for generating virtual samples, and then propagations are executed based on these expense-less virtual samples.
The performance of the active sampling-based KRG algorithm is validated using an LEO case considering the atmospheric drag. Numerical simulation shows that the proposed active sampling-based KRG algorithm can accurately approximate the orbit state using only 10 samples.
Moreover, the proposed AGMM-VSG is successfully applied to solve the orbit uncertainty propagation problem in LEO, HEO, and Tianqin orbit cases. It was demonstrated that the obtained PDF contours using AGMM-VSG are more representative of the curvature of the true distribution than the fifth-order CR method. Moreover, likelihood agreement measure tests indicate that the new approach shows close performance in approximating the propagation of uncertainty with the AGMM method, while only 100 evaluations are required. The computational cost of the proposed AGMM-VSG method is only 1/20 that of the conventional GMM method. All in all, the newly presented approach has advantages in accuracy and efficiency over competitive methods.

Author Contributions

Conceptualization, T.X. and Z.Z.; methodology and software, H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. U20B2001), the National Key R&D Program of China (Grant No. 2020YFC2201200) and the Basic Scientific Research Project (Grant No. JCKY2020903B002).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wen, C.; Peng, C.; Gao, Y. Reachable domain for spacecraft with ellipsoidal Delta-V distribution. Astrodynamics 2018, 2, 265–288. [Google Scholar] [CrossRef]
  2. Qiao, D.; Zhou, X.; Li, X. Analytical configuration uncertainty propagation of geocentric interferometric detection constellation. Astrodynamics 2023. [Google Scholar] [CrossRef]
  3. Oguri, K.; McMahon, J.W. Stochastic primer vector for robust low-thrust trajectory design under uncertainty. J. Guid. Control Dyn. 2022, 45, 84–102. [Google Scholar] [CrossRef]
  4. Shu, P.; Yang, Z.; Luo, Y.Z. Higher-Order Lambert Problem Solution Based on Differential Algebra. J. Guid. Control Dyn. 2022, 45, 1913–1926. [Google Scholar] [CrossRef]
  5. Folkner, W.; Hechler, F.; Sweetser, T.; Vincent, M.; Bender, P. LISA orbit selection and stability. Class. Quantum Gravity 1997, 14, 1405. [Google Scholar] [CrossRef] [Green Version]
  6. Tan, Z.; Ye, B.; Zhang, X. Impact of orbital orientations and radii on TianQin constellation stability. Int. J. Mod. Phys. D 2020, 29, 2050056. [Google Scholar] [CrossRef]
  7. Zhou, X.; Qin, T.; Meng, L. Maneuvering Spacecraft Orbit Determination Using Polynomial Representation. Aerospace 2022, 9, 257. [Google Scholar] [CrossRef]
  8. Ye, B.B.; Zhang, X.; Zhou, M.Y.; Wang, Y.; Yuan, H.M.; Gu, D.; Ding, Y.; Zhang, J.; Mei, J.; Luo, J. Optimizing orbits for TianQin. Int. J. Mod. Phys. D 2019, 28, 1950121. [Google Scholar] [CrossRef]
  9. Luo, Y.Z.; Yang, Z. A review of uncertainty propagation in orbital mechanics. Prog. Aerosp. Sci. 2017, 89, 23–39. [Google Scholar] [CrossRef]
  10. Terejanu, G.; Singla, P.; Singh, T.; Scott, P.D. Uncertainty propagation for nonlinear dynamic systems using Gaussian mixture models. J. Guid. Control Dyn. 2008, 31, 1623–1633. [Google Scholar] [CrossRef]
  11. Han, H.; Qiao, D.; Chen, H.; Li, X. Rapid planning for aerocapture trajectory via convex optimization. Aerosp. Sci. Technol. 2019, 84, 763–775. [Google Scholar] [CrossRef]
  12. Fuller, A. Analysis of nonlinear stochastic systems by means of the Fokker–Planck equation. Int. J. Control. 1969, 9, 603–655. [Google Scholar] [CrossRef]
  13. Chen, N.; Majda, A.J.; Tong, X.T. Rigorous Analysis for Efficient Statistically Accurate Algorithms for Solving Fokker–Planck Equations in Large Dimensions. SIAM/ASA J. Uncertain. Quantif. 2018, 6, 1198–1223. [Google Scholar] [CrossRef] [Green Version]
  14. Junkins, J.L.; Akella, M.R.; Alfrined, K.T. Non-Gaussian error propagation in orbital mechanics. Guid. Control. 1996, 1996, 283–298. [Google Scholar]
  15. Ghoreishi, S.; Allaire, D. Adaptive uncertainty propagation for coupled multidisciplinary systems. AIAA J. 2017, 55, 3940–3950. [Google Scholar] [CrossRef]
  16. Jia, B.; Xin, M. Active Sampling Based Polynomial-Chaos–Kriging Model for Orbital Uncertainty Propagation. J. Guid. Control Dyn. 2021, 44, 905–922. [Google Scholar] [CrossRef]
  17. Park, R.S.; Scheeres, D.J. Nonlinear mapping of Gaussian statistics: Theory and applications to spacecraft trajectory design. J. Guid. Control Dyn. 2006, 29, 1367–1375. [Google Scholar] [CrossRef] [Green Version]
  18. Chen, Q.; Qiao, D.; Wen, C. Reachable domain of spacecraft after a gravity-assist flyby. J. Guid. Control Dyn. 2019, 42, 931–940. [Google Scholar] [CrossRef]
  19. Chen, Q.; Qiao, D.; Wen, C. Orbital Element Reachable Set After Gravity Assists of Planets in Elliptical Orbits. J. Guid. Control Dyn. 2020, 43, 989–997. [Google Scholar] [CrossRef]
  20. Yang, Z.; Luo, Y.Z.; Zhang, J. Nonlinear semi-analytical uncertainty propagation of trajectory under impulsive maneuvers. Astrodynamics 2019, 3, 61–77. [Google Scholar] [CrossRef]
  21. Shu, P.; Yang, Z.; Luo, Y.Z.; Sun, Z.J. Collision Probability of Debris Clouds Based on Higher-Order Boundary Value Problems. J. Guid. Control Dyn. 2022, 45, 1512–1522. [Google Scholar] [CrossRef]
  22. Park, R.S.; Scheeres, D.J. Nonlinear semi-analytic methods for trajectory estimation. J. Guid. Control Dyn. 2007, 30, 1668–1676. [Google Scholar] [CrossRef] [Green Version]
  23. Majji, M.; Junkins, J.L.; Turner, J.D. A high order method for estimation of dynamic systems. J. Astronaut. Sci. 2008, 56, 401–440. [Google Scholar] [CrossRef]
  24. Zhou, X.; Qin, T.; Ji, M.; Qiao, D. A LSTM assisted orbit determination algorithm for spacecraft executing continuous maneuver. Acta Astronaut. 2022. [Google Scholar] [CrossRef]
  25. Boone, S.; McMahon, J. Orbital guidance using higher-order state transition tensors. J. Guid. Control Dyn. 2021, 44, 493–504. [Google Scholar] [CrossRef]
  26. Qiao, D.; Zhou, X.; Zhao, Z.; Qin, T. Asteroid Approaching Orbit Optimization Considering Optical Navigation Observability. IEEE Trans. Aerosp. Electron. Syst. 2022, 58, 5165–5179. [Google Scholar] [CrossRef]
  27. Wang, Z.; Wang, X.; Liang, Y.; Yang, F. A new sigma-point filter–uniform random sampling Kalman filter. IFAC-PapersOnLine 2017, 50, 3853–3858. [Google Scholar] [CrossRef]
  28. Li, X.; Qiao, D.; Li, P. Bounded trajectory design and self-adaptive maintenance control near non-synchronized binary systems comprised of small irregular bodies. Acta Astronaut. 2018, 152, 768–781. [Google Scholar] [CrossRef]
  29. Han, H.; Qiao, D.; Chen, H. Optimization of aeroassisted rendezvous and interception trajectories between non-coplanar elliptical orbits. Acta Astronaut. 2019, 163, 190–200. [Google Scholar] [CrossRef]
  30. Jia, B.; Xin, M. Short-arc orbital uncertainty propagation with arbitrary polynomial chaos and admissible region. J. Guid. Control Dyn. 2020, 43, 715–728. [Google Scholar] [CrossRef]
  31. Julier, S.J.; Uhlmann, J.K.; Durrant-Whyte, H.F. A new approach for filtering nonlinear systems. In Proceedings of the American Control Conference-ACC’95, Seattle, WA, USA, 21–23 June 1995; Volume 3, pp. 1628–1632. [Google Scholar]
  32. Van Der Merwe, R. Sigma-point Kalman Filters for Probabilistic Inference in Dynamic State-Space Models; Oregon Health & Science University: Portland, OR, USA, 2004. [Google Scholar]
  33. Meng, D.; Miao, L.; Shao, H.; Shen, J. A seventh-degree cubature Kalman filter. Asian J. Control. 2018, 20, 250–262. [Google Scholar] [CrossRef]
  34. Cui, X.; Jing, Z.; Luo, M.; Guo, Y.; Qiao, H. A new method for state of charge estimation of lithium-ion batteries using square root cubature Kalman filter. Energies 2018, 11, 209. [Google Scholar] [CrossRef] [Green Version]
  35. Zhang, X.; Yan, Z.; Chen, Y. High-degree cubature Kalman filter for nonlinear state estimation with missing measurements. Asian J. Control. 2022, 24, 1261–1272. [Google Scholar] [CrossRef]
  36. Chen, P.; Zabaras, N.; Bilionis, I. Uncertainty propagation using infinite mixture of gaussian processes and variational bayesian inference. J. Comput. Phys. 2015, 284, 291–333. [Google Scholar] [CrossRef] [Green Version]
  37. Li, X.; Warier, R.R.; Sanyal, A.K.; Qiao, D. Trajectory tracking near small bodies using only attitude control. J. Guid. Control Dyn. 2019, 42, 109–122. [Google Scholar] [CrossRef]
  38. Jones, B.A.; Weisman, R. Multi-fidelity orbit uncertainty propagation. Acta Astronaut. 2019, 155, 406–417. [Google Scholar] [CrossRef]
  39. Sun, P.; Colombo, C.; Trisolini, M.; Li, S. Hybrid Gaussian Mixture Splitting Techniques for Uncertainty Propagation in Nonlinear Dynamics. J. Guid. Control Dyn. 2022, 1–11. [Google Scholar] [CrossRef]
  40. Vittaldev, V.; Russell, R.P.; Linares, R. Spacecraft uncertainty propagation using gaussian mixture models and polynomial chaos expansions. J. Guid. Control Dyn. 2016, 39, 2615–2626. [Google Scholar] [CrossRef] [Green Version]
  41. Vittaldev, V.; Russell, R.P. Space object collision probability using multidirectional gaussian mixture models. J. Guid. Control Dyn. 2016, 39, 2163–2169. [Google Scholar] [CrossRef]
  42. DeMars, K.J.; Bishop, R.H.; Jah, M.K. Entropy-based approach for uncertainty propagation of nonlinear dynamical systems. J. Guid. Control Dyn. 2013, 36, 1047–1057. [Google Scholar] [CrossRef]
  43. Zhou, X.; Cheng, Y.; Qiao, D.; Huo, Z. An adaptive surrogate model-based fast planning for swarm safe migration along halo orbit. Acta Astronaut. 2022, 194, 309–322. [Google Scholar] [CrossRef]
  44. Zhou, X.; Wang, S.; Qin, T. Multi-Spacecraft Tracking and Data Association Based on Uncertainty Propagation. Appl. Sci. 2022, 12, 7660. [Google Scholar] [CrossRef]
  45. Tang, Y.; Long, T.; Shi, R.; Wu, Y.; Gary Wang, G. Sequential radial basis function-based optimization method using virtual sample generation. J. Mech. Des. 2020, 142, 111701. [Google Scholar] [CrossRef]
  46. Shi, R.; Long, T.; Baoyin, H. Multi-fidelity and multi-objective optimization of low-thrust transfers with control strategy for all-electric geostationary satellites. Acta Astronaut. 2020, 177, 577–587. [Google Scholar] [CrossRef]
  47. Zhao-Ming, L.; Wen-Ge, Y.; Dan, D.; Yu-Rong, L. A novel algorithm of fifth-degree cubature Kalman filter for orbit determination at the lower bound approaching to the number of cubature points. Acta Phys. Sin. 2017, 66. [Google Scholar] [CrossRef]
  48. Vishwajeet, K.; Singla, P. Adaptive split/merge-based Gaussian mixture model approach for uncertainty propagation. J. Guid. Control Dyn. 2018, 41, 603–617. [Google Scholar] [CrossRef]
  49. Omran, A.; Newman, B. Nonlinearity Index Theory for Flight Mechanics Applications. In Proceedings of the AIAA Atmospheric Flight Mechanics Conference, Toronto, ON, Canada, 2–5 August 2010; p. 7624. [Google Scholar]
  50. Manchester, Z.; Kuindersma, S. Derivative-free trajectory optimization with unscented dynamic programming. In Proceedings of the IEEE 55th Conference on Decision and Control (CDC), Las Vegas, NV, USA, 12–14 December 2016; pp. 3642–3647. [Google Scholar]
  51. Bouhlel, M.A.; Martins, J.R. Gradient-enhanced kriging for high-dimensional problems. Eng. Comput. 2019, 35, 157–173. [Google Scholar] [CrossRef] [Green Version]
  52. Gao, A.; Liao, W. Efficient gravity field modeling method for small bodies based on Gaussian process regression. Acta Astronaut. 2019, 157, 73–91. [Google Scholar] [CrossRef]
  53. Sun, F.; Su, W.Y.; Wang, M.Y.; Wang, R.J. RBF-POD reduced-order modeling of flow field in the curved shock compression inlet. Acta Astronaut. 2021, 185, 25–36. [Google Scholar] [CrossRef]
  54. Jia, F.; Li, X.; Qiao, D.; Zhou, X. Semi-analytical configuration optimization of geocentric gravitational wave observatory. Acta Astronaut. 2023, 202, 522–534. [Google Scholar] [CrossRef]
Figure 2. Active sampling for one-dimensional numerical benchmark.
Figure 2. Active sampling for one-dimensional numerical benchmark.
Applsci 13 03069 g002
Figure 3. PDF contours using AGMM-VSG with Monte Carlo samples of HEO case.
Figure 3. PDF contours using AGMM-VSG with Monte Carlo samples of HEO case.
Applsci 13 03069 g003
Figure 4. PDF contours using fifth-order CR with Monte Carlo samples of HEO case.
Figure 4. PDF contours using fifth-order CR with Monte Carlo samples of HEO case.
Applsci 13 03069 g004
Figure 5. Likelihood agreement measure for the AGMM, AGMM-VSG, and fifth-order CR of HEO case.
Figure 5. Likelihood agreement measure for the AGMM, AGMM-VSG, and fifth-order CR of HEO case.
Applsci 13 03069 g005
Figure 6. PDF contours using AGMM-VSG with Monte Carlo samples of LEO case.
Figure 6. PDF contours using AGMM-VSG with Monte Carlo samples of LEO case.
Applsci 13 03069 g006
Figure 7. PDF contours using fifth-order CR with Monte Carlo samples of LEO case.
Figure 7. PDF contours using fifth-order CR with Monte Carlo samples of LEO case.
Applsci 13 03069 g007
Figure 8. Likelihood agreement measure for the AGMM, AGMM-VSG, and fifth-order CR of LEO case.
Figure 8. Likelihood agreement measure for the AGMM, AGMM-VSG, and fifth-order CR of LEO case.
Applsci 13 03069 g008
Figure 9. Final uncertainty using AGMM-VSG and MC for TQO case.
Figure 9. Final uncertainty using AGMM-VSG and MC for TQO case.
Applsci 13 03069 g009
Figure 10. PDF contours using AGMM-VSG with Monte Carlo samples for TQO case.
Figure 10. PDF contours using AGMM-VSG with Monte Carlo samples for TQO case.
Applsci 13 03069 g010
Table 1. Description of point, response and sample.
Table 1. Description of point, response and sample.
ItemSymbolDescriptionOther Statement
P o i n t x t k Input of the black-box modelPure coordinate
R e s p o n s e x t k + 1 Output of the black-box model-
S a m p l e x t k , x t k + 1 Combination of a point and its response-
Table 2. Five-component splitting library with β = 2 and λ = 0.0025 [42].
Table 2. Five-component splitting library with β = 2 and λ = 0.0025 [42].
Component i α i m ˜ i σ ˜ i
10.0763216491−1.68997291110.4422555386
20.2474417860−0.80092838340.4422555386
30.352473130000.4422555386
40.24744178600.80092838340.4422555386
50.07632164911.68997291110.4422555386
Table 3. Elements of the nominal orbit of LEO case.
Table 3. Elements of the nominal orbit of LEO case.
a/kmei/deg Ω /deg ω /degn/deg
659600000
Table 4. Elements of the nominal orbit of HEO case.
Table 4. Elements of the nominal orbit of HEO case.
a/kmei/deg Ω /deg ω /degn/deg
35,0000.20000
Table 5. Elements of the nominal orbit of TQO case.
Table 5. Elements of the nominal orbit of TQO case.
a/kmei/deg Ω /deg ω /degn/deg
99,995.5723230.00043094.697997210.445892358.62446361.329603
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

Xu, T.; Zhang, Z.; Han, H. Adaptive Gaussian Mixture Model for Uncertainty Propagation Using Virtual Sample Generation. Appl. Sci. 2023, 13, 3069. https://doi.org/10.3390/app13053069

AMA Style

Xu T, Zhang Z, Han H. Adaptive Gaussian Mixture Model for Uncertainty Propagation Using Virtual Sample Generation. Applied Sciences. 2023; 13(5):3069. https://doi.org/10.3390/app13053069

Chicago/Turabian Style

Xu, Tianlai, Zhe Zhang, and Hongwei Han. 2023. "Adaptive Gaussian Mixture Model for Uncertainty Propagation Using Virtual Sample Generation" Applied Sciences 13, no. 5: 3069. https://doi.org/10.3390/app13053069

APA Style

Xu, T., Zhang, Z., & Han, H. (2023). Adaptive Gaussian Mixture Model for Uncertainty Propagation Using Virtual Sample Generation. Applied Sciences, 13(5), 3069. https://doi.org/10.3390/app13053069

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