Multiple Model-Based Synchronization Approaches for Time Delayed Slaving Data in a Space Launch Vehicle Tracking System

Due to the inherent characteristics of the flight mission of a space launch vehicle (SLV), which is required to fly over very large distances and have very high fault tolerances, in general, SLV tracking systems (TSs) comprise multiple heterogeneous sensors such as radars, GPS, INS, and electrooptical targeting systems installed over widespread areas. To track an SLV without interruption and to hand over the measurement coverage between TSs properly, the mission control system (MCS) transfers slaving data to each TS through mission networks. When serious network delays occur, however, the slaving data from the MCS can lead to the failure of the TS. To address this problem, in this paper, we propose multiple model-based synchronization (MMS) approaches, which take advantage of the multiple motionmodels of an SLV. Cubic spline extrapolation, prediction through an α-β-γ filter, and a single model Kalman filter are presented as benchmark approaches.We demonstrate the synchronization accuracy and effectiveness of the proposed MMS approaches using the Monte Carlo simulation with the nominal trajectory data of Korea Space Launch Vehicle-I.


Introduction
The range safety system (RSS) [1] for the flight mission of a space launch vehicle (SLV) consists of multiple heterogeneous tracking systems (TSs) with mission control systems (MCSs).Critical launch mission details such as time-space-position information (TSPI), launch mission status data, that is, quick look message (QLM), and flight safety information can be acquired from the RSS.Tracking an extensive mission trajectory of an SLV requires widespread multiple TSs so that the RSS covers the entire trajectory of the SLV flight mission.Generally, multiple TSs are spread out over different sites and they automatically hand over observation coverage according to the flight of the SLV.In this circumstance, one of the most important roles of the RSS is to distribute slaving data to each TS for continuous tracking of the SLV.If a critical network delay results in time delayed slaving data to be sent to the TSs, the MCS will not receive accurate SLV tracking data.This problem can lead to significant difficulties for the SLV mission progress and analysis.Therefore, the motivation of this research is to enhance slaving data accuracy by compensating for possible network time delays.The basic solution to this problem is simple (linear or nonlinear) extrapolation of the filtered data.In this case, extrapolation methods simply propagate the TSPI data without regard to the system dynamics of the SLV.On the other hand, a Kalman filter (KF) and its prediction capability [2,3] can reflect the system dynamics of the SLV, which results in better synchronization performance.However, since a KF only utilizes a single dynamic model, in general, the tracking performance of a KF for a maneuvering target is inferior to multiple model estimators [4].In addition, due to stage separation, the flight phase of an SLV is separated into two parts, the propelled flight phase (PFP) and the coasting flight phase (CFP).Hence, the dynamic model of an SLV can be described using multiple models so that the multiple flight phases are properly accounted for.To adaptively select one of the multiple dynamic models according to the flight phase of the SLV, multiple model estimators (MME) such as an interacting multiple model (IMM) [5] and a multiple model adaptive estimator (MMAE) [6] could significantly improve SLV tracking when a network delay results in delayed slaving data transmission.
In the past several decades, considerable research has been undertaken in the field of launch vehicle tracking based on multiple dynamic models; researchers have shown interest in various applications such as the tracking of reentry vehicles, short-range projectiles, and sounding rockets [7][8][9][10][11][12][13]. A reentry vehicle tracking problem known as highly nonlinear dynamics was conducted using a modified IMM, with a different algorithm cycle compared with an ordinary IMM with multiple modes of diverse ballistic coefficients [7].Short-range ballistic munitions or projectiles with multiple models such as spin-stabilized and fin-stabilized models were implemented using an IMM [8,9].Both research alternatives can be applied in the impact point prediction of projectiles.Tactical ballistic missile tracking was carried out using an IMM estimator with three modes [10].The first mode was a constant axial force model for the boosting and reentry phases.The second mode was a ballistic acceleration model that incorporated the gravitational, Coriolis, and centripetal forces for the exoatmospheric phase.The final mode was a standard autocorrelated acceleration Singer's model for malfunction motions of missiles such as reentry tumbling.A multiple IMM algorithm with an unbiased mixing approach for multiple modes of thrusting or for ballistic projectiles was presented [11,12].A sounding rocket with multiple modes of propelled flight or free fall flight was tracked using a multiple model adaptive estimation approach [14].
In this paper, multiple model-based synchronization (MMS) approaches are proposed to synchronize the time delayed slaving data of the RSS.The proposed approaches can be expressed via two distinct multiple models, a nonlinear model and a linear model.The nonlinear model considers comprehensive factors such as thrust, gravity, drag coefficient, Mach number, and air density [10][11][12][13].Although the nonlinear model precisely describes the motion of the SLV, it requires complex information concerning the SLV specifications to be collected in advance.In contrast, in the case of the linear dynamic model, a simple constant velocity (CV) or constant acceleration (CA) model with multiple hypotheses, which takes advantage of Singer's model [16,17], can be utilized [14].To describe the motion of the SLV, the motion modes of both models are separated into two parts, PFP and CFP.We propose a slaving data synchronization approach for the RSS based on MME so that the MCS can adaptively find an appropriate dynamic model at an arbitrary time index, where time delay has occurred.The performance of slaving data synchronization is compared to various benchmark methods such as cubic spline extrapolation, prediction through an -- filter, and a single model KF to demonstrate the effectiveness of the proposed algorithm.
The remainder of the paper is organized as follows.Section 2 presents a statement of the problem for delayed slaving data in RSS.In Section 3, conventional synchronization approaches are illustrated.In Section 4, two proposed MMEbased synchronization approaches are derived.Section 5 presents simulation settings and results; the comparison between different types of synchronization approaches is depicted as an aspect of RMS error of the state vector.Finally, in Section 6, the conclusions of this paper are presented.

Problem Statement for Delayed Slaving Data in RSS
The transmission of slaving data from the MCS to multiple TSs facilitates seamless tracking of the SLV in a sparsely located multiple TS environment.If a data transmission delay problem occurs, it can cause an SLV tracking failure.As depicted in Figure 1, the antennas of TSs are pointing at the SLV by controlling their attitude according to slaving data from the MCS.In this situation, a slight time delay in slaving data can give rise to large differences between the antenna beam and the SLV due to the fast movement of the SLV.
where   is a linear or nonlinear propagation function and  = [  ]  is a slaving state vector that is composed of x-axis, y-axis, and z-axis positions but is not limited to the position components.
Let us assume a synchronized slaving state vector to be an unknown function of the delayed slaving state vector whose values are known only until time .We then define a cubic spline extrapolation function  ep  ( + ), where  = , ,  at a specific synchronization time  +  that is extrapolated based on the known function   () that has a real value, with  + 1 points, where  0 ≤  ≤   .

Synchronization
Using an -- Filter.When the state estimation covariance for a time invariant system converges under suitable conditions to a steady-state value, explicit expressions of the steady-state covariance and filter gain can be obtained.The resulting steady-state filters for noisy kinematic models are known as - and -- filters [20].In this paper, a combined -- filter using both an expanding memory polynomial filter (EMF) and a fading memory polynomial filter (FMF) was used [21,22].At first, the EMF is represented as follows: where Δ is the sampling time and   , , V  , , and   , are the position, velocity, and acceleration estimates of the EMF, respectively.In addition, the variance reduction factor (VRF) [22] of the EMF, that is, VRF  , can be represented as On the other hand, the FMF and its VRF (VRF  ) can be written as follows: where 0 <  < 1.Both filters are conducted in parallel but only one of them is selected as a final estimate by comparing the VRFs.During the early tracking phase, the EMF is selected as the final estimate; however, at a certain time, where VRF  is larger than VRF  , FMF is selected as the final estimate: where Finally, synchronization using an -- filter is completed by linear propagation using the system matrix   such that

Kalman Predictor.
The motion of the SLV is simply depicted as a discretized Wiener process acceleration model [20]: where the state vector   ∈ R 9 consists of the position, velocity, and acceleration components along axis, -axis, and -axis, respectively; that is, .The system matrix and covariance matrix of the system noise   are represented as (10) and (11), respectively: A radar measurement for the SLV gives the spherical coordinate observations such that where V , ∼ N(0,  , ) and  , = diag( 2  ,  2  ,  2  ).Using a 3D debiased converted measurement [23], we can transform the original nonlinear equations (12) into their linear form as Here, V , is the converted measurement noise expressed in terms of Cartesian coordinates; that is, V  ∼ N(0,   ): Prediction or fixed-lead prediction in mean square means the synchronization of the slaving data is the estimation of the state at a future time +, where  > 0 beyond the observation interval; that is, based on data up to an earlier time [20,24], The optimal predictor or synchronized state X+| ≜ XKP +| and its error covariance  +| ≜  KP +| are given by the Kalman predictor equations [2,24]: where XKF | and  KF | are the KF estimate and covariance, respectively.

Proposed MME-Based Synchronization Approaches
4.1.Synchronization Using an IMM with Singer's Linear Model.Singer [16] described 2D manned maneuvering targets in range-bearing coordinates.This model can be adapted to the SLV kinematics in 3D Cartesian coordinates with multiple flight phases: where   ∈ R 3×1 is a white noise process along the Cartesian axis  = , , .The parameter  = 1/  is the maneuver correlation time constant, and  2  is the acceleration variance describing maneuver intensity.In a steady state, To describe SLV kinematics, the model must cope with crucial nonzero mean acceleration maneuvers during the propelled phase.In addition, after each stage's engine burns out, a multiple model approach is applied to the coasting flight; that is, one model describes PFP, whereas the other depicts the CFP.Empirically tuned, independent probability density functions (PDFs) represented by TUM describe the accelerations of the SLV in the local coordinate frame.Figure 2 shows the means and variances of the acceleration processes of the SLV in this paper.The discrete-time model with state transition matrix Ψ  (  , Δ) is as follows: In the case of the PFP,  , is a nonzero mean white noise sequence caused by the nonzero mean acceleration   seen in Figure 2. A nonzero mean white noise sequence for the PFP should be considered in the target kinematics when implementing the state propagation stage in the KF.Thus, the deterministic input  , caused by  , along -axis, -axis, and -axis is derived as follows: -axis and -axis can be derived in the same manner as shown in (20).The maneuver excitation covariance [10], which represents the uncertainty of the SLV kinematics model, is ] .
In addition, the measurement matrix   for Singer's model can be depicted as where ]  is measurement noise (as shown in ( 14)) with error covariance.

IMM with Singer's Linear
Model.From ( 19)-( 22) in Section 4.1, we can rewrite the Markov jump linear systems, where the th model of the finite multiple model set M = { (1) , . . .,  () } obeys the following equations: Here,  2  for the SLV can be represented by TUM as in Figure 2. The PDFs of Figures 2(a)-2(d) are experimentally sampled from the nominal acceleration profile of the SLV.The superscript (i) denotes quantities pertinent to model  () in M, and the jumps of the system mode are assumed to have transition probabilities: where  ()  denotes the event in which model  () matches the system mode in effect at time .In our application,  = 2, and  (1) and  (2) are the propelled and coasting modes, respectively.

Synchronization Using an IMM with a Nonlinear Ballistic
Model.For a nonlinear ballistic model, the state vector for the propelled mode is denoted as where   is the drag coefficient and   is the thrust.Generally, the drag coefficient varies significantly with the Mach number regime: subsonic, transonic, and supersonic.Therefore, we take advantage of the dynamic model considering a Mach number-dependent multiplier [11,12] such that (31) The first term on the right side of (31) represents the thrust (m 2 /s) of the SLV in , , and  directions.Two distinct multiple modes of a nonlinear SLV model can be separated by the existence of thrust.In other words, the existence of thrust signifies propelled flight mode, whereas zero thrust signifies coasting flight mode. is the magnitude of the velocity V = [ ẋ ẏ ż ]  , that is, the SLV speed (m/s).The second part of (31) is the drag term, which is related to velocity and altitude; that is,  = −()/2, where () =  0  − is the air density (kg/m 3 ) at an altitude  (m) and  is the air density constant (m −1 ) [9]. is the drag coefficient and   is the Mach numberdependent drag coefficient multiplier, which is approximated by the cubic spline curve shown in Figure 3.In this paper, for the drag characteristics of the SLV, which are applied in the subsequent simulation section, drag coefficients of the Saturn  launch vehicle [15] are used.The third part of (31) is a gravity term.Gravity  is the standard acceleration due to gravity at sea level, which is assumed to be constant throughout the trajectory, with a value of 9.812 m/s 2 .ω1 , ω2 , and ω3 are assumed to be continuous time zero mean white Gaussian noises.The drag coefficient and thrust acceleration are represented as Wiener processes with a slow variation [11,12].We can modify the dynamic equations ( 30) and (31) as a compact form such that where The state vector equation ( 33) is discretized by a secondorder Taylor expansion [26].Then, (33) can be written as a discretized continuous time system with white process noise such that where   is the Jacobian of (33) evaluated at x  [26] and   is the discretized continuous time process noise for the sampling interval Δ.The corresponding covariance matrix of the discretized process noise is where  3 is the 3 × 3 identity matrix and the continuous time process noise intensities  V ,   , and   are the corresponding power spectral densities.The measurement matrix  NL for the nonlinear ballistic multiple model can be depicted as where ]  is measurement noise with error covariance that is, ( 14).An IMM algorithm for nonlinear dynamics with different sizes of the mode state vector is summarized as follows [25]: (v) Model-conditioned reinitialization (for  = 1, 2, . . ., ): (a) predicted mode probability: (vi) Model-conditioned filtering (for  = 1, 2, . . ., ): (a) predicted estimate and covariance: The nonlinear MMS of time delayed slaving data is completed by propagating the combined estimate based on the current mode's estimated vector such that where  is a lead-time for synchronization ( =  ⋅ Δ).If the SLV is in the PFP ( (1)  NL, >  (2)  NL, ), the current state estimate includes the thrust term; that is, the state vector x NL, = [      ẋ  ẏ  ż      ] .If the SLV is not in the PFP, the state vector does not include the thrust term (  = 0).

Simulation Results
To demonstrate the performance of the proposed synchronization approaches for delayed slaving data, we simulated the SLV tracking problem based on the nominal flight trajectory of the Korea Space Launch Vehicle-I (KSLV-I).In the simulation, the radar measurement noise intensities of ( 12  small motion.A single KF prediction model cannot exactly describe the motion of both PFP and CFP, whereas the proposed multiple model-based approaches work well.As shown in Figure 4, regardless of the delay , the proposed multiple model-based approaches find synchronized position vectors.In the comparison of synchronization performance between linear IMM and nonlinear IMM synchronization, the nonlinear IMM-based synchronization approach shows the best performance, except for -axis, where CA motion is dominant.In addition, the difference between the proposed multiple model-based approaches is very small as shown in Figure 4, but the complexities of the algorithms for real-time applications are quite dissimilar.Hence, the operator may adaptively select one of the proposed approaches according to one's environment.

Conclusions
In this paper, we investigated the time synchronization approaches of delayed slaving data in the RSS for SLV tracking.One of the most important roles of the RSS is to distribute slaving data to each TS for continuous tracking of the SLV.If there is a critical network delay resulting in time delayed slaving data being sent to each TS, the MCS will not receive accurate SLV tracking data.This problem can give rise to significant difficulties for the SLV mission progress and analysis.To overcome this problem, we proposed MMS approaches which take advantage of the multiple motion models of an SLV.The linear IMM-based synchronization approach was developed using Singer's model with ternary uniform mixtures and the nonlinear IMM-based synchronization approach was derived from a nonlinear ballistic model with a drag coefficient.For verification of the proposed algorithms, SLV tracking simulations using KSLV-I and the radar measurement data generated from nominal trajectory were conducted.To demonstrate the superiority of time synchronization performance in these simulations, we compared the proposed algorithm with benchmark approaches for absolute error between the nominal trajectory data and the synchronized slaving data; the simulation results demonstrated that the proposed MMS approaches performed competitively.

Figure 1 :
Figure 1: Illustration of problem statement for delayed slaving data in RSS.

Figure 2 :
Figure 2: Nonzero mean acceleration PDF in the PFP model (a) along -axis, (b) along -axis, and (c) along -axis.Zero mean acceleration PDF in the CFP Model (d) along -axis, -axis, and -axis.

Figure 4 :
Figure 4: Synchronization errors at PFP and CFP.
, ,  , , and  , are unknown coefficients.To find the unknown coefficients, we make the assumption that   () should be continuous in [ 0 ,   ].Therefore, the following equation containing unknown coefficients  , ,  , ,  , , and  , is obtained:  , (  ) =  ,+1 (  ) ) are selected as   = 15 m,   = 0.01 deg, and   = 0.01 deg.The nominal flight sequence of KSLV-I is as follows.First, the payload fairing separates during the first stage flight at 215.4 s.After the first stage engine shutdown at 228.7 s, the upper stage separates from the first stage and enters the CFP.The second stage continues in the CFP until the kick motor ignition at 395 s, and the vehicle then enters the PFP.At the end of kick motor combustion (452.7 s), the upper stage enters the target orbit in CFP as in the previous separation.Finally, the satellite is inserted into the target orbit after it separates from the upper stage during CFP at 540 s.Since delayed slaving data from the MCS to the TSs can occur in both the PFP and CFP, synchronization simulations are conducted at arbitrary points of the flight phases.Figure4shows the synchronization error of the delayed slaving data with respect to the benchmark and the proposed approaches for delays of  = our application.This is why the performance of the Kalman prediction-based synchronization is worse than the proposed approaches.Nonetheless, we can see that the -axis synchronization errors shown in Figures4(b) and 4(e) are small.This result is observed because the -axis motion of KSLV-I in a local coordinate is relatively smaller than the other axis motion, and the CA model approximately expresses this Mathematical Problems in Engineering