4.1.1. Parameter Setting
To facilitate the display and comparison of experimental results, the initial position of the offshore target T1 is used as a reference point, and some parameters are set in the local geographic coordinate system centered on the reference point. The data rate is set to
, which means every 0.05 s, the platform sensor observes the sea target to obtain observation information, and the inertial guide updates its own attitude information and low-precision geographic coordinate information. The coordinates of UAVs A and B in the local geographic coordinate system of the reference point are
,
; the standard deviations of sensor errors (distance, azimuth, azimuth) and attitude angle errors (yaw, pitch, roll) of UAVs A and B are shown in
Table 3 and
Table 4, respectively.
Where
and
are the geodesic coordinate systematic error scaling factor and attitude angle systematic error scaling factor, respectively, which are used to facilitate the adjustment of error variation in the simulation experiment. The flow of simulation data generation is shown in
Figure 3.
The motion posture is set as follows: the initial position of the sea target T2 with respect to the reference point is 198,000 m in geodetic length and 0 in geodetic azimuth. Both sea targets T1 and T2 move 2000 m in the direction of geodetic azimuth of 0 in 200 s.
The root mean square error (RMSE) is used for the error evaluation index. The root mean square error of the systematic error of the geodesic coordinate estimation and the position estimation error are calculated as
where
can be taken as
,
,
. T is the total duration.
,
denotes the latitude systematic error and longitude systematic error of the UAV at moment
, and
is the position estimation error.
In order to visualize the degree of error reduction, the error reduction rate
is calculated by Equation (53)
4.1.2. Simulation Experimental Result
In the above scenario taken from
,
, the simulation is carried out in the set experimental scenario. The estimated effect of the geodesic coordinate systematic error of UAVs A and B is shown in
Figure 4 and
Figure 5.
Figure 3 and
Figure 4 show the relationship between estimated value of the geodesic coordinate systematic error and the true value of the geodesic coordinate systematic error of UAVs A and B. It is obvious that the estimated value of systematic error all converges to the true value.
Table 5 shows the estimated error of the geodetic coordinate systematic error, which is already a small amount compared to the set true value. The systematic error reduction rate exceeds 70% in different platforms. It shows that the proposed geodetic coordinate systematic error estimation method is effective. The proposed geodetic coordinate systematic error estimation method estimated the systematic error by using EKF in essence. Therefore, the initial value of the system error should be set reasonably.
Figure 6 shows the tracking effect using different algorithms in scenario one, and
Table 6 shows the position estimation error of different algorithms.
Figure 6 shows that it is obvious that IRT has the best tracking result compared to comparing algorithms. The reason for the error in EC1-PLKF and EC1-EKF is that these error estimation methods do not take the effect of attitude angle error and random error in geodesic coordinates into account. EC2-PLKF and EC2-EKF do not take the effect of geodesic coordinate error into account. RT ignores the effect of attitude angle error and geodesic error. In addition, the estimation of the nonlinear filter is slightly better than that of the linear filter.
Figure 7 shows the tracking position error versus time, from which the conclusion of the tracking performance of various algorithms is approximately the same as above. In addition, it can be obtained that the tracking methods using EC1 and EC2 have a better performance in terms of convergence speed, converging in about 5 s, while RT and IRT have a larger initial error and converge slightly slower than the other algorithms because they use the observed value of the target position as the initial tracking condition, which has a larger error compared to other algorithms that have gone through initial compensation.
Table 7 shows the average single-step time of different algorithms, where RT has the shortest time because it does not need to estimate any systematic error value but directly estimates the target state. IRT has the longest time because it adds a new proposed geodesic systematic error estimation algorithm compared to RT and considers the influence of attitude angle. Although the time consumption is increased, it remains in the millisecond range and does not affect the performance of the algorithm for practical application.
In order to fully represent the performance of each algorithm under different systematic errors, simulation experiments are performed for , varying from 1 to 100.
Figure 8 shows the variation of the position estimation error of different algorithms as the geodetic coordinate systematic error increases, where
varies from 1 to 100 and
, focusing on the influence of geodetic coordinate systematic error on the tracking algorithm.
Figure 9 and
Table 8 reflect the position estimation error data under typical geodetic coordinate systematic error.
Figure 8 shows that IRT, EC1-PLKF, and EC1-EKF can effectively suppress the effect brought by the error as the geodesic coordinate systematic error increases, while the error of EC2-PLKF, EC2-EKF, and RT keeps increasing. When
is larger than 20, the tracking effect of EC2-PLKF, EC2-EKF, and RT starts to be worse than that of IRT because the influence of the geodesic coordinate systematic error on the target position estimation starts to dominate at this time.
When is greater than 5, the tracking effect of EC1-PLKF and EC1-EKF is gradually superior to EC2-PLKF and EC2-EKF. When the geodesic systematic error is small, random error affects the estimation effect of the systematic error of method EC1, which affects the tracking effect. When the systematic error of the geodesic coordinate becomes larger, the effect of random error gradually becomes smaller, and the tracking effect of EC1-PLKF and EC1-EKF is always better than that of EC2-PLKF and EC2-EKF.
When the systematic error increases, the position estimation error of IRT becomes slightly larger because IRT assumes that the geodesic coordinate is true when constructing the coordinate transformation of the observed data by measurement equations. When the geodesic coordinate systematic error keeps increasing, the error from this assumption will keep increasing.
Figure 10 represents the variation of position estimation error for different algorithms as the attitude angle systematic error increases, where
varies from 1 to 100 and
, focusing on the effect of attitude angle systematic error on the tracking algorithm.
Figure 11 and
Table 9 reflect the position estimation error under typical attitude angle systematic error.
The EC1-PLKF and EC1-EKF have the fastest tracking error growth rate because their error estimation algorithms do not estimate the attitude error, and the effect of the increasing attitude error is directly applied to the observation data. The EC2-EKF has a certain suppression effect on the attitude systematic error, but the estimation capability of the algorithm is limited and cannot completely suppress the effect caused by the attitude systematic error. RT and IRT are least affected by the attitude systematic error, and the tracking effect of IRT improves by about 25% compared with RT when the attitude systematic error is 0.1° and about 13% compared with RT when the attitude systematic error is 1°, indicating that the tracking effect of IRT improves with the decreasing attitude error. The reason for this is that when considering the attitude systematic error, it is assumed that the attitude angle is true during the derivation process, and when the attitude angle systematic error is increasing, the error brought by this assumption will increase continuously. IRT The method should be used in the presence of geodetic coordinate system errors that are not excessive.
In addition, the tracking performance of the EKF filter is slightly better than that of the PLKF filter, as seen in
Table 8 and
Table 9, which can improve the tracking accuracy by about 20 m to 30 m.