Adaptive Multi-Sensor Joint Tracking Algorithm with Unknown Noise Characteristics

In this study, to solve the low accuracy of multi-space-based sensor joint tracking in the presence of unknown noise characteristics, an adaptive multi-sensor joint tracking algorithm (AMSJTA) is proposed. First, the coordinate transformation from the target object to the optical sensors is considered, and the observation vector-based measurement model is established. Then, the measurement noise characteristics are assumed to be white Gaussian noise, and the measurement covariance matrix is set as a constant. On this premise, the traditional iterative extended Kalman filter is applied to solve this problem. However, in most actual engineering applications, the measurement noise characteristics are unknown. Thus, a forgetting factor is introduced to adaptively estimate the unknown measurement noise characteristics, and the AMSJTA is designed to improve the tracking accuracy. Furthermore, the lower bound of the proposed algorithm is theoretically proved. Finally, numerical simulations are executed to verify the effectiveness and superiority of the proposed AMSJTA.


Introduction
The monitoring and tracking technology of near-Earth objects is currently one of the research hotspots in the field of aerospace.As the prerequisite and foundation of planetary defense, it has been a topic of concern that has received extensive attention from major institutions and scholars in various related fields [1][2][3].Currently, the increasing prevalence of missile technology as a primary method of aerospace operation has led to numerous countries attaining mastery in this field, establishing it as an independent and effective combat approach.Consequently, the development of anti-missile defense technology and enhancement of missile defense capabilities have become crucial requirements for strategic defense systems.The primary function of a space-based early warning satellite system is to detect incoming threats (mainly ballistic missiles) and provide timely warnings and realtime detection data.A ballistic missile is a type of aircraft that utilizes ballistic control to follow predetermined flight paths.During their weightless free-flight phase, which is also the mid-segment of the flight process, warheads traverse predetermined orbits.Therefore, this segment also serves as a pivotal research area for orbit tracking and determination within early warning systems [4].
At present, the principal detection methods for near-Earth objects (NEOs) mainly include ground-based telescopes and space-based telescopes based on optical sensors and ground-based radars based on electromagnetic waves.As for ground-based radars, owing to the constraints of antenna aperture, power consumption, detection range, and susceptibility to atmospheric conditions and illumination, there might be difficulties, such as limited site selection, an insufficient observation time, a narrow detection area, and more constraints during the NEO process, which makes it difficult to accurately co-locate in complex environments, and a large number of optimization algorithms are required to determine the layout method for collaborative positioning [5][6][7].In comparison, space-based optical sensors achieve real-time searching and tracking for NEOs, with the advantages of low energy consumption, high efficiency, wide area coverage, great measurement accuracy, long detection distances, etc.Moreover, the arrangement of space-based optical sensors is more flexible and less affected by the atmospheric environment [8].These advantages make space-based optical detection a developing trend and the main direction for future NEO observation equipment.A constellation of several satellites also displays great advantages.Compared with a single spacecraft, a constellation has advantages with respect to capability, flexibility, and robustness.Space-based optical detection proves to be a crucial technique to attain space situational awareness.However, optical sensors can only obtain angle measurement information on NEOs and cannot directly obtain distance information.This passive angle positioning method cannot fully reflect the orbital information of the target.In order to obtain the three-dimensional position information of the target, at least two or more optical observation carriers are required to participate in localization.The progression of the cooperative orbit determination without the help of the ground station's aid is one of the features of a satellite constellation.From the perspective of the task level, the basic tasks of space-based detection include various aspects, such as orbit maneuvering, attitude control, imaging models, etc. [9] Usually, comprehensive analytical assessments and numerical simulations of observability are required in the simulation process [10].In this study, the filtering problems in the determination of a target's trajectory are the main focus.
One of the major challenges in multi-space-based sensor joint tracking scenarios is the low accuracy due to the presence of unknown noise characteristics.In tracking applications, the use of the Cartesian coordinate system for the modeling process and the fact that the measurements are usually in a polar or spherical coordinate system, as well as the motion model of the space target, all lead to the issue that the target tracking process is essentially a nonlinear filtering problem [11].Nonlinear filtering algorithms, such as the extended Kalman filter, unscented Kalman filter, and particle filter, can perform well in their respective applicable scenarios to improve measurement accuracy [12][13][14][15].There are also many current studies that have used various types of improved filtering methods for radar measurements or sensor measurements of a target's motion trajectory.For instance, the tracking problem of high-speed gliding target tracking with unknown time-varying maneuvers was investigated in [16], and a robust extended Kalman filter based on the compensation of maneuver observers was proposed to solve the model mismatch problem by using the output of maneuver observers to correct the prediction step of the dynamics.A method for initial trajectory determination was proposed in [17], which can be applied to space-based space target monitoring, solving the effective result and reducing the susceptibility to the initial value of distance in space-based space target monitoring.In [18], a method of quantized extended Kalman filtering under quantized measurement conditions was derived, and the measurement update was based on the conditional mean estimation of the given quantized measurements, which improves the accuracy of trajectory tracking under the condition of quantized noise.There are also methods that utilize artificial intelligence to assist filtering.For example, in [19], an information fusion technology based on artificial neural networks was proposed for multi-sensor integrated autonomous navigation.Most of the current methods have been able to realize high precision under Gaussian noise, but most of them require large amounts of feedback data.We hope to offer an available option to solve the problem of low accuracy in tracking in the presence of unknown noise characteristics when the data are not abundant.Furthermore, the tracking of target trajectories is also of great importance in the task allocation of spacebased systems [20][21][22].
In this study, we propose an adaptive multi-sensor joint tracking algorithm (AMSJTA).First, the orbiting and tracking model of NEOs based on space-based optical observation is investigated.Then, the orbital dynamics model of NEOs and the observation model of optical sensors are established.According to the acquisition mode and characteristics of space-based optical goniometric information, the corresponding tracking process is provided.Noting that the values of the process noise covariance matrix Q and the measurement noise covariance matrix R are highly influential on the values of the extended Kalman filter [23], the adaptive method is introduced here to adaptively update the process noise and the measurement noise at each step of the filtering process, so as to improve the performance of the iterative extended Kalman filter and enhance the accuracy of the filtering.A simple derivation of the Bayesian CRLB for the estimation error is also carried out as an intuitive alternative to the estimation error to be used as a tracking performance metric for the target.
The remainder of this article is organized as follows: Section 2 presents the dynamic observation model of space-based multi-optical sensors for orbital targets.In Section 3, the AMSJTA filtering method is specified, and the related derivation procedure is provided.Section 4 presents a discussion and the simulation results of an application of the proposed algorithm in a specific scenario.Finally, Section 5 presents the conclusions.

Target Motion and Measurement Model
The models in this study are mainly implemented in the Earth-centered, Earth-fixed coordinate system (ECEF), the Earth-centered inertial coordinate system (ECI), and the coordinate system based on the sensor itself, which are described as follows: 1.
The ECEF coordinate system is a typical Cartesian coordinate system, which is fixed relative to the Earth.The origin of the coordinate system O is located in the center of the Earth, the OX axis is in the equatorial plane and points to the meridian where the Greenwich Observatory is located, the OZ axis is perpendicular to the equatorial plane and points toward the North Pole direction when combined with the axis of the rotation of the Earth, and the relationship between the OX, OY, and OZ axes satisfies the right-hand rule.

2.
The ECI coordinate system is also a typical Cartesian coordinate system, which is stationary relative to the fixed star.It is usually considered to be in the inertial space, with the geocentric point O as the coordinate origin.OX is in the equatorial plane pointing to the equinox point.OZ is perpendicular to the equatorial plane and points to the North Pole, and the orientation of the OY axis conforms to the right-hand rule.

3.
In the orbital coordinate system in which the coordinate origin O of the sensor body coordinate system is the geometric center of the sensor, the OX axis is the extension of the line connecting the Earth's center and the coordinate origin, the OZ axis is located in the plane formed by the OX axis and the Earth's rotation axis and is perpendicular to the OX axis, and the orientation of the OY axis satisfies the right-hand rule.
According to the principle of coordinate rotation transformation, the transfer matrix from the orbital coordinate system to the ECI coordinate system is as follows: where ω is the argument of periapsis, i is the inclination, and Ω is the longitude of the ascending node.Firstly, the basic structure and function of the space-based multi-sensor system are analyzed and studied, the implementation process of the space-based multi-sensor early warning and detection mission along with the working mode and detection capability of the space-based multi-sensor are discussed, and the working model in the middle part of the ballistic trajectory tracking of multiple targets is established.
When solving the observation vector, it is necessary to obtain the position information of the observation satellite.The position vector of the observation satellite is first calculated according to the orbital mechanics.Considering that the motion of the LEO satellite around the Earth can be approximated as only being acted upon by the gravitational force of the Earth, if the Earth is equivalent to a mass point, the Earth and the satellite constitute a simple two-body system as shown in Figure 1.The orbits can be constructed by using the orbital elements for orbit reconstruction [24].In this study, the method of dynamical equations is adopted to establish the continuous state equation of the space target.The fourth-order Runge-Kutta method is applied to obtain the nonlinear discrete equation of the space target [25].
satellite around the Earth can be approximated as only being acted upon by the gravitational force of the Earth, if the Earth is equivalent to a mass point, the Earth and the satellite constitute a simple two-body system as shown in Figure 1.The orbits can be constructed by using the orbital elements for orbit reconstruction [24].In this study, the method of dynamical equations is adopted to establish the continuous state equation of the space target.The fourth-order Runge-Kutta method is applied to obtain the nonlinear discrete equation of the space target [25].In the ECEF coordinate system, it is assumed that the position and velocity of space targets are expressed as follows: Considering the effect of the J2 perturbation term, the orbital dynamics equations for the space target are as follows: where . According to these formulas, the continuous equation of the state of the space target is expanded by the fourth-order Runge-Kutta method.The acceleration can be obtained according to Equation (2).Then, the nonlinear discrete equation of the state of the space target can be expressed as follows: The process noise covariance matrix is set as  = (     ).According to [26], it can be described using the Singer model in the following form: In the ECEF coordinate system, it is assumed that the position and velocity of space targets are expressed as follows: Considering the effect of the J2 perturbation term, the orbital dynamics equations for the space target are as follows: where According to these formulas, the continuous equation of the state of the space target is expanded by the fourth-order Runge-Kutta method.The acceleration can be obtained according to Equation (2).Then, the nonlinear discrete equation of the state of the space target can be expressed as follows: The process noise covariance matrix is set as According to [26], it can be described using the Singer model in the following form: For the signals acquired by the sensors, a spatial polar coordinate system is applied for processing as shown in Figure 2. The measurement vectors of the observed target in the coordinate system are defined using the direction angle and the pitch angle as follows: Sensors 2024, 24, x FOR PEER REVIEW 5 of 16  = [  5  20  4  8  3  6  4  8  3  3 For the signals acquired by the sensors, a spatial polar coordinate system is applied for processing as shown in Figure 2. The measurement vectors of the observed target in the coordinate system are defined using the direction angle and the pitch angle as follows: For the optical sensors, they can only acquire the angle information of the target.As a result, in order to determine the spatial location of the target, we need the observation results of at least two sensors and their spatial location data.In the case of dual-satellite observation and orbit determination, measurement errors are considered as follows: Using the observed angle information of at least two sensors, along with the spatial information of the sensors themselves, we have enough information to obtain the spatial information of the observation target.Therefore, the observation can be determined as the spatial position coordinates of the target constructed from the direction and pitch angles obtained from at least two telescope satellites.The observation error is defined as: For optical sensors, since they can only obtain the angle information during observation, we can assume that there is a positioning ray that starts from the center of the optical sensor and points in the direction of the target.Under ideal conditions, the positioning rays of multiple satellites for a single observation target should intersect in space.However, due to the influence of various systematic errors and random errors, it is very likely that these rays from multiple satellites will have no common intersection point in space under practical conditions.Using a system of linear equations to directly solve the tracking problem may often result in the inability to yield the final outcomes.In the case of two For the optical sensors, they can only acquire the angle information of the target.As a result, in order to determine the spatial location of the target, we need the observation results of at least two sensors and their spatial location data.In the case of dual-satellite observation and orbit determination, measurement errors are considered as follows: Using the observed angle information of at least two sensors, along with the spatial information of the sensors themselves, we have enough information to obtain the spatial information of the observation target.Therefore, the observation can be determined as the spatial position coordinates of the target constructed from the direction and pitch angles obtained from at least two telescope satellites.The observation error is defined as: For optical sensors, since they can only obtain the angle information during observation, we can assume that there is a positioning ray that starts from the center of the optical sensor and points in the direction of the target.Under ideal conditions, the positioning rays of multiple satellites for a single observation target should intersect in space.However, due to the influence of various systematic errors and random errors, it is very likely that these rays from multiple satellites will have no common intersection point in space under practical conditions.Using a system of linear equations to directly solve the tracking problem may often result in the inability to yield the final outcomes.In the case of two satellites cooperating, the midpoint of the common perpendicular between two positioning rays in space can be used instead of the point of intersection for the operation as shown in Figure 3.A and B are the intersection points of the positioning rays from two satellites and the common perpendicular.However, in the case of more than two satellites observing simultaneously, Newton's method can be used to iteratively obtain the point in space that minimizes the sum of the distances from each ray [27].
satellites cooperating, the midpoint of the common perpendicular between two positioning rays in space can be used instead of the point of intersection for the operation as shown in Figure 3.A and B are the intersection points of the positioning rays from two satellites and the common perpendicular.However, in the case of more than two satellites observing simultaneously, Newton's method can be used to iteratively obtain the point in space that minimizes the sum of the distances from each ray [27].It is noted that since the final results of the extended Kalman filter are highly influenced by the values of the process noise covariance matrix Q and the measurement noise covariance matrix R, improper values of Q and R can greatly reduce the accuracy of the filter results and even lead to filter dispersion.In most studies of dynamic state prediction, the values of Q and R are set as fixed constants during the estimation process, and it is often difficult to obtain sufficiently good results for dynamic processes by using the exhaustive method.Therefore, the adaptive estimation method used in this study is to adaptively adjust the measurement noise covariance matrix R and the process noise covariance matrix Q.Then, R and Q are updated according to the current state values at each step of the iterative extended Kalman filtering so as to further improve the filtering results.

Adaptive Multi-Sensor Joint Tracking Algorithm
Since the observation equations are nonlinear, extended Kalman filtering is required for the state estimation of the target.According to the established motion equations and observation equations, the traditional extended Kalman filtering process is as follows in Algorithm 1.The status transition matrix F in extended Kalman filtering is different from the case in linear filtering and can be replaced by the Jacobian matrix.The Taylor series expansion can be retained as the third-or fourth-order terms.However, the second-order term is most usually applied for the comprehensive consideration of both computational complexity and performance.It is noted that since the final results of the extended Kalman filter are highly influenced by the values of the process noise covariance matrix Q and the measurement noise covariance matrix R, improper values of Q and R can greatly reduce the accuracy of the filter results and even lead to filter dispersion.In most studies of dynamic state prediction, the values of Q and R are set as fixed constants during the estimation process, and it is often difficult to obtain sufficiently good results for dynamic processes by using the exhaustive method.Therefore, the adaptive estimation method used in this study is to adaptively adjust the measurement noise covariance matrix R and the process noise covariance matrix Q.Then, R and Q are updated according to the current state values at each step of the iterative extended Kalman filtering so as to further improve the filtering results.

Adaptive Multi-Sensor Joint Tracking Algorithm
Since the observation equations are nonlinear, extended Kalman filtering is required for the state estimation of the target.According to the established motion equations and observation equations, the traditional extended Kalman filtering process is as follows in Algorithm 1.The status transition matrix F in extended Kalman filtering is different from the case in linear filtering and can be replaced by the Jacobian matrix.The Taylor series expansion can be retained as the third-or fourth-order terms.However, the second-order term is most usually applied for the comprehensive consideration of both computational complexity and performance.
At the same time, an iterative approach is applied, where the results of the k th moments that have been obtained are processed iteratively again, and the outcome after each iteration is used as the input for the next iteration so that a better performance can be achieved for the filtering results for each state X k .

Residuals of predictions:
v k = Z k+1 − H k+1 Xk+1|k Update process of extended Kalman filtering: Exit the loop when there is no new measured state quantity 4. Output : X k When the data available are limited, we can apply the least squares method for the determination of the initial orbit.The original trail can be determined through secondorder polynomial fitting using the initial data.We can define the residual as the difference between the actual observation and the estimated observation.With the idea of the least squares method, we can minimize the sum of squared residuals of the system within the observation interval.During the iterative process, the covariance can be written as follows: (9)   After each iteration, the covariance proceeds in a decreasing trend.For each iteration, the update of X k is: The iteration is terminated when the reduction in the covariance between two iterations falls below a threshold.
In this study, an adaptive filtering approach is adopted for further improvement in the effect of the EKF to obtain better optimization results, as shown in the following equation: where v k is the error vector at time k, which can be taken as the residual vector or the state error depending on the object to be adapted and can be used to perform the corresponding covariance matrix, R k or Q k .According to [28], the following expressions can be obtained: The traditional extended Kalman filter algorithm has difficulties in accurately determining the system process noise and the measurement noise.Some existing adaptive methods have deficiencies in estimating the covariance of measurement noise while up-dating the state.Therefore, in order to cope with this challenge, this study proposes an improved adaptive method to determine R k and Q k in the EKF process.The covariance matching method can adjust its covariance matrix according to the theoretical value of innovation.In the prediction step of the EKF, the innovation is the difference between the actual measured value and its predicted value, which can be calculated as follows: Compared with some methods that directly use the residuals estimated at the current moment, here, a sliding window is used for the ε k ε k T approximation by taking the average value as follows: Meanwhile, this study introduces the forgetting factor α to adaptively estimate R k , which satisfies the following equation: A smaller α gives more weight to the previous estimate and, therefore, prevents excessive fluctuations in the adaptive process and longer time delays in order to adapt to changes.According to the principle of automated control, d is taken as 0.95.
A brief representation of the AMSJTA process in this study is shown below in Figure 4:  Due to the one-step prediction error of the state, the following expressions can be obtained.
The filtering residuals: The state estimation error: Due to the one-step prediction error of the state, the following expressions can be obtained.
The filtering residuals: The state estimation error: A series of separately defined positive definite matrices exist as follows: According to [29], it is known that such a construction exists; therefore, the initial condition can be satisfied as follows: Then, for the Kalman gain K k+1 at time k − 1, the following expressions can be applied: Next, pseudo-system state transfer matrices and measurement matrices are constructed using statistical linear error propagation methods: Then, the original state transfer equation and measurement equation can be linearized as follows: x where α k and β k are the compensating diagonal matrices used to minimize possible estimation errors.The one-step prediction of the information matrix is as follows: Considering the nonlinear stochastic system presented earlier and the adaptive iterative extended Kalman filtering algorithm described in the previous section, the error is exponentially bounded on the mean square under these assumptions [30].
Finally, a brief discussion of the existence of error bounds for adaptive filtering is presented.Consider the discrete nonlinear system as described below: Since v k and w k are white noise independent of each other here and independent of the initial state density function of the system, assume that the set of measurements up to moment k is denoted as Z k = {z i } k i=1 , the unbiased estimate of X k is denoted as Xk|k , and the covariance matrix is denoted as P k|k .The lower bound of this matrix is denoted as CRLB and satisfies the following equation: where I k is the Fisher information matrix (FIM).
For the above Gaussian model, the FIM can be written as the sum of the following two components: where I k corresponds to the a priori information about the state of the target, and I Z (P k , x k ) is the part generated during the measurement process.
The calculations can be simplified by utilizing predicted values instead of statistical averages, as given by the following equation:

Simulation and Validation
In this section, we present simulation experimental scenarios and results to demonstrate that the improved adaptive iterative extended Kalman filter has better performance in applications.Moreover, in addition to the comparison in terms of performance, we also investigate the possible effects of changes in the initial conditions on the filtering results.A simulation scenario of the joint observation of multiple space-based optical transmitters for a low-orbit target is constructed to validate the performance of the AMSJTA in this study.The parameters used in the simulation are listed in the following Table 1.
For the simulation of the orbit of the near-Earth satellites, both satellites carry spacebased optical sensors.According to Section 2, the method of the four-stage Runge-Kutta expansion under the consideration of the influence of the Earth's J 2 perturbation term is adopted.The initial state of the two satellites and the space target is shown in the following Table 2 in ECEF coordinates.The trajectories of the satellites and a space target in the ECEF coordinate system are shown in the following Figure 5.The sensors of each satellite need to acquire angular degrees in their respective coordinate systems and jointly determine the orbit in the ECI coordinate system, the true trajectory in the ECI coordinate system, the observation trajectory, and the trajectory after adaptive iterative extended Kalman filtering.
The root-mean-square error (RMSE) of the target position can be used to evaluate the performance of the filtering effect.The tracking results of the iterative extended Kalman filter and the filter with the addition of adaptive tuning are compared, as shown in Figure 6.In the results in the figure, it can be observed that after adding the adaptive process, the method improves the tracking accuracy.Therefore, the designed AMSJTA has a better performance compared with the conventional IEKF filter, and the filtered error is basically maintained at around 0.1%.The measurement period of each of the two observation satellites for the target is dt = 1 s.The standard deviation and covariance matrices of the sensor's random measurement noise are as follows: The sensors of each satellite need to acquire angular degrees in their respective coordinate systems and jointly determine the orbit in the ECI coordinate system, the true trajectory in the ECI coordinate system, the observation trajectory, and the trajectory after adaptive iterative extended Kalman filtering.
The root-mean-square error (RMSE) of the target position can be used to evaluate the performance of the filtering effect.The tracking results of the iterative extended Kalman filter and the filter with the addition of adaptive tuning are compared, as shown in Figure 6.In the results in the figure, it can be observed that after adding the adaptive process, the method improves the tracking accuracy.Therefore, the designed AMSJTA has a better performance compared with the conventional IEKF filter, and the filtered error is basically maintained at around 0.1%.In order to better demonstrate the advantages of the adaptive method, the initial values are artificially biased to observe the performance of the traditional filtering meth and the filtering method with adaptive improvement.The results of 1000 Monte Ca simulations are taken in the experiment.
First, take Q = 10Q 0 .In this case, according to Figure 7, it can be clearly obser that the conventional IEKF process estimation bias is significantly improved, while use of adaptive methods can still maintain a relatively acceptable filtering performan greatly enhancing the stability of the original algorithm.In another case, instead of Gaussian white noise, the inserted system noise Q is m to be a random noise that gradually decreases with time, oriented as shown in Figur In this example, the adaptive processing of the filtering can still converge to a suita accuracy relatively quickly, but the error of the conventional filtering method becom larger and even appears to have a possible tendency toward divergence.In order to better demonstrate the advantages of the adaptive method, the initial Q k values are artificially biased to observe the performance of the traditional filtering method and the filtering method with adaptive improvement.The results of 1000 Monte Carlo simulations are taken in the experiment.
First, take Q = 10Q 0 .In this case, according to Figure 7, it can be clearly observed that the conventional IEKF process estimation bias is significantly improved, while the use of adaptive methods can still maintain a relatively acceptable filtering performance, greatly enhancing the stability of the original algorithm.In order to better demonstrate the advantages of the adaptive method, the initial   values are artificially biased to observe the performance of the traditional filtering method and the filtering method with adaptive improvement.The results of 1000 Monte Carlo simulations are taken in the experiment.
First, take Q = 10Q 0 .In this case, according to Figure 7, it can be clearly observed that the conventional IEKF process estimation bias is significantly improved, while the use of adaptive methods can still maintain a relatively acceptable filtering performance, greatly enhancing the stability of the original algorithm.In another case, instead of Gaussian white noise, the inserted system noise Q is made to be a random noise that gradually decreases with time, oriented as shown in Figure 8.In this example, the adaptive processing of the filtering can still converge to a suitable accuracy relatively quickly, but the error of the conventional filtering method becomes larger and even appears to have a possible tendency toward divergence.In another case, instead of Gaussian white noise, the inserted system noise Q is made to be a random noise that gradually decreases with time, oriented as shown in Figure 8.In this example, the adaptive processing of the filtering can still converge to a suitable accuracy relatively quickly, but the error of the conventional filtering method becomes larger and even appears to have a possible tendency toward divergence.The results obtained from the simulation experiments are basically consistent with the previous conclusions.The results of this study demonstrate that the proposed adaptive algorithm can obtain a better filtering performance and improve the robustness of the extended Kalman filtering algorithm.

Conclusions
As the most concise and direct method for nonlinear problems such as orbital target tracking, the extended Kalman filter is deficient in stability.The proposed AMSJTA method compensates for the shortcomings of the EKF's instability and significantly improves the accuracy of tracking orbital targets.The improved algorithm combines the least-squares criterion and adaptive iterative extended Kalman filtering for orbit determination.Based on the covariance matching theory, this study established and utilized the known residual quantities of the filtering process to update the measurement state noise and the process noise covariance.The proposed adaptive algorithm improves the accuracy of the estimation of target orbits and achieves better performance than traditional extended Kalman filtering with unknown measurement noise characteristics.The results obtained from the simulation experiments are basically consistent with the previous conclusions.The results of this study demonstrate that the proposed adaptive algorithm can obtain a better filtering performance and improve the robustness of the extended Kalman filtering algorithm.

Conclusions
As the most concise and direct method for nonlinear problems such as orbital target tracking, the extended Kalman filter is deficient in stability.The proposed AMSJTA method compensates for the shortcomings of the EKF's instability and significantly improves the accuracy of tracking orbital targets.The improved algorithm combines the least-squares criterion and adaptive iterative extended Kalman filtering for orbit determination.Based on the covariance matching theory, this study established and utilized the known residual quantities of the filtering process to update the measurement state noise and the process noise covariance.The proposed adaptive algorithm improves the accuracy of the estimation of target orbits and achieves better performance than traditional extended Kalman filtering with unknown measurement noise characteristics.

Figure 1 .
Figure 1.A geometric diagram of multiple space-based optical sensors for orbital target observations.

Figure 1 .
Figure 1.A geometric diagram of multiple space-based optical sensors for orbital target observations.

Figure 2 .
Figure 2. Diagram of polar coordinate observation of a target in the sensor body coordinate system.

Figure 2 .
Figure 2. Diagram of polar coordinate observation of a target in the sensor body coordinate system.

Figure 4 .
Figure 4. Diagram of the adaptive multi-sensor joint tracking algorithm process.

Figure 4 .
Figure 4. Diagram of the adaptive multi-sensor joint tracking algorithm process.

Figure 5 .
Figure 5. Satellite and space target trajectories.The measurement period of each of the two observation satellites for the target is dt = 1 s.The standard deviation and covariance matrices of the sensor's random measurement noise are as follows: σ r = 1000m, σ a = 0.1°, σ e = 0.1° R 0 = diag([σ r , σ a , σ e ])

Figure 5 .
Figure 5. Satellite and space target trajectories.

Figure 8 .
Figure 8.(a) Diagram of the inserted system noise Q change over time.(b) Comparison of RMSE after IEKF and AMSJTA filtering when Q decreases over time.

Figure 8 .
Figure 8.(a) Diagram of the inserted system noise Q change over time.(b) Comparison of RMSE after IEKF and AMSJTA filtering when Q decreases over time.

Table 2 .
Initial state vectors for space units.