Robust GNSS Positioning Using Unbiased Finite Impulse Response Filter

: In a typical GNSS receiver, pseudorange and pseudorange rate measurements are gener ‐ ated through the code and carrier tracking loops, respectively. These measurements are then em ‐ ployed to calculate the user’s position and velocity (PV) solutions, which is typically achieved using a Kalman filter (KF) or the least squares (LS) algorithm. However, the LS method only uses the current observation without error analysis. The positioning result is greatly affected by the errors in the observed data. In KF, by using an iterative approach that combines predictions and measure ‐ ments of PV information, more accurate estimates can be obtained because the PV information is time ‐ correlated. Meanwhile, its optimal estimate requires that both the model and noise statistics are exactly known. Otherwise, achieving optimality cannot be guaranteed. To address this issue, this paper proposes and implements a novel GNSS solution method based on an unbiased finite impulse response (UFIR) filter. Two different field tests were conducted. The position results of UFIR are compared with those from the LS and KF methods, and the horizon positioning mean error is improved by 44% and 29%, respectively, which highlights its efficacy. The method offers two primary benefits: it is robust to noise uncertainty, and it leverages historical data within the UFIR framework to provide a more accurate estimate of the current state.


Introduction
Nowadays, location-based services (LBS) have seamlessly integrated into our daily lives and are playing a crucial role in various fields, including civilian and military applications [1][2][3][4].By utilizing a constellation of satellites, the Global Navigation Satellite System (GNSS) has the capability to provide users with global precise positioning, navigation, and timing (PVT) information.Users who possess a GNSS receiver can obtain PVT information if there are sufficient satellites in view [5,6].In a typical GNSS receiver, satellite signals are first down-converted and sampled through an RF front-end.Then, intermediate frequency (IF) signals are utilized for signal acquisition and tracking.From the tracked signals, the receiver calculates pseudorange and pseudorange rate measurements.Finally, based on these measurements, the GNSS receiver can further process these data to calculate its PVT information, which commonly employs either a least squares (LS) method or a Kalman filtering (KF) approach [7,8].
The LS algorithm is a classic mathematical optimization technique.By minimizing the sum of the squared residuals (between observations and estimations), LS can effectively find the best-fitting solution that optimally fits the observed data.Several attempts were made to enhance the LS estimation performance in GNSS.For example, Filić et al.
proposed a weighted LS to take the expected errors into account [9] to reduce the computational burden caused by matrix inversion in the LS.Barbu et al. suggested a novel QR factorization method that effectively improved the precise point positioning performance of high-frequency GNSS [10].Li et al. conducted a study where they integrated Doppler measurements with code pseudorange data to improve LS position estimation [11].María et al. applied an LS solution network over all the recorded data generated by several GNSS devices to achieve sub-meter horizontal accuracy [12].Xu et al. introduced a bias-corrected weighted LS method for accurate location services that auto-corrects model nonlinearity and squared range biases while sequentially estimating unknown location parameters with condition adjustments [13].The LS converges quickly in the initial iteration, but there is no error analysis on the observed data, and it is thus greatly affected by the errors in the current observations.In addition, these LS methods overlook the time correlation inherent in the GNSS measurement data.For the KF approach, the state vector usually consists of the quantities that define the system's state, such as position, velocity, clock bias, and drift.Meanwhile, the observation vector typically includes pseudorange and pseudorange rate increments [14,15].Based on this, smoother results can be obtained by incorporating state transformation constraints.For the KF to achieve optimality, it requires precise knowledge of the noise statistics.However, in practice, it is difficult to obtain exact knowledge of noise statistics for a time-varying system.To address the noise challenges associated with the KF, scholars have proposed numerous methods to address the noise sensitivity issue inherent in the Kalman filter (KF).For instance, Gao et al. introduced an adaptive KF navigation algorithm, termed RL-AKF.The uniqueness of this method lies in its ability to adaptively estimate the process noise covariance matrix by leveraging the principles of reinforcement learning [16].In [17], Ahmed introduced a novel DRKF technique that constructs a statistically robust equivalent weight model using the chi-square test, thereby enhancing placement outcomes.To address the GNSS outliers, Berkay et al. integrated robust KF with variance component estimation to dynamically set weights and variances for multi-GNSS observations in single-frequency positioning [18].
Another approach to approximating noise statistics by using innovation about the process is referred to as the adaptive KF (AKF), but the approximation is limited [15,[19][20][21].In addition, the KF operates in a recursive manner, making it an infinite impulse response (IIR) filter, which allows estimation errors to accumulate over time during the estimation process [22].Instead, a superior filter named the unbiased finite impulse response (UFIR) filter has been emerging as a remarkable replacement.It has been extensively studied and thoroughly analyzed in [23].In the UFIR filter, a fixed-size horizon of past data is collected and processed together to estimate the current state, effectively suppressing the propagation of historical errors to the current epoch [24].Another advantage of UFIR is that it is insensitive to noise interference [25].The UFIR filter stands out for its ability to minimize the mean-square error (MSE) while efficiently utilizing a horizon length of N measurements.Fortunately, the optimal N can be determined in several ways [26].
Because of the above, the primary purpose of this work is to explore the implementation and application of the UFIR method in GNSS solutions to improve the accuracy of position estimation.In addition to its independence from noise statistics, this paper presents several contributions compared to existing LS-and KF-based position estimation methods in GNSS, which are as follows: (1) The proposed method offers a novel approach to estimating position information that effectively eliminates the requirement for calculating error covariances during the iteration process.(2) Because the UFIR operates only with the opt N adjacent points in a batch form, the accumulation of errors is ultimately suppressed.
(3) Two field tests were conducted, and comparative analyses were performed with LS and KF methods.The results show that the UFIR solution can significantly improve GNSS positioning performance.
This paper is structured as follows.In Section 2, some formulas of GNSS positioning are given, and we also review the LS and KF approaches.In Section 3, we present the framework of the UFIR and its implementation in GNSS.In Section 4, two field tests are conducted, and the results of LS, KF, and UFIR are also presented, compared, and analyzed.Section 5 concludes this paper.

Problem Formulation
As mentioned earlier, the user PVT information is derived using the pseudorange and pseudorange rate observations.These observations are generated from two different tracking loops, which include the code tracking loop and the carrier tracking loop.This section sequentially introduces the parameterization model in operation and the estimation methods.Specifically, Section 2.1 presents GNSS formulas for position calculation; Section 2.2 gives the principle of the LS method; and Section 2.3 reviews the KF and its difference with LS.

GNSS Position Calculation
Figure 1 illustrates the geometric structure relating the satellite positions to the receiver's location.Assuming there are i visible satellites, the pseudorange formulas are expressed as follows [6]:

Least Square Method
In the above, each pseudorange measurement equation is a nonlinear function of the user position.After Taylor expansion and ignoring high-order terms, the relation between the offset x and pseudorange values  can be expressed as [6]: where is the state increment of user receiver and  is the pseudorange increment; they represent the provision of the approximate solution around the linearization point and the linearization of the satellite-receiver geometric range, referring to The receiver clock offset u t  is not affected by the linearization; H is the direction cosine matrix, given by: and the [ ] e e e  e denotes the unit vector from the user location to the th i satellite, obtained by follows: ,e ,e r r r where is the approximate distance between the user and the th i satellite.The  can be further expressed as: According to Equation ( 2), the state increment of the user receiver can be written as: Accordingly, one can express the update of the position estimation and the receiver clock bias by adding the position increment x to the current state, as follows: As we can see from the above, the LS method considers that each measurement is independent of each other, and it provides a single estimate of the position that fits the observed measurements.This indicates that the relationship between these measurements has not been investigated.

Kalman Filter
The KF is typically a state-space approach utilized for estimating the state of a dynamic system while considering noisy observations.It involves the use of both the state vector and the measurement vector in the estimation process.In general, the adjacent state vector x can be modeled as [27]: where and c t    present the velocity deviation of the user and clock drift.
Thus, the state propagation equation at epoch k can be given by: in which  denotes the system transition matrix; w is the process noise term; I de- notes an identity matrix; and  is the update interval of the KF.
The process noise covariance matrix Q that related to w is determined based on the rule of thumb, which takes into consideration the user dynamics and oscillator used.The calculation is performed as follows [27]: where v S represents noise power spectral density (PSD) of the receiver velocity, set by the level of user dynamics; and f S and g S are the PSD of receiver clock phase and fre- quency, respectively, calculated as follows: where 0 h and 2 h  represent the clock model coefficients associated with the white frequency modulation noise and random walk frequency modulation noise, respectively, associated with the utilized oscillator.
The covariance matrix k P of the priori sate estimation can be expressed as: The linearized GNSS measurement equation is given as follows: where H is the measurement matrix, given by: in which  and    are the pseudorange increment and pseudorange rate increment, respectively; k v is the measurement noise vector.
Unlike the process noise, the measurement noise depends on the particular phase discriminator implementation, and its covariance matrix is calculated by: diag( ( )) where ( ) x  is the variance of x, diag• represents a diagonal matrix, and is referred to as the measurement residual or innovation sequence.Therefore, it suffices to update the posteriori state estimation based on the innovation k  , using [28]: in which k K denotes the state feedback gain matrix, given by: The updated covariance matrix is defined by: Incorporating the updated state estimation, k x , the final state of the user receiver is obtained by: Unlike the LS, the KF considers that each measurement is time-correlated because state transformation constraints are included.In addition, KF makes an error analysis of the current estimation and aims to minimize the variance of the state error.Compared to the results obtained from the LS method, utilizing the KF approach yields smoother and more authentic results.Nonetheless, because of the sensitivity to noise uncertainty and because the historical information has not been explored, the improvement is limited.A framework that possesses historical information constraints and that is noise independent is anticipated to produce superior estimate results compared to the LS and KF.

Optimization Framework
In this section, the UFIR optimization framework is proposed to obtain superior position estimates.To this end, the extended state-space model is first constructed.Then, the posterior estimation, error covariance, and gain of UFIR are described in detail.Finally, the way to find the optimal horizon length is given.

Extended Model
For UFIR, the extended state vector n,m

X and observation vector n,m
Y are repre- sented in a horizon length of N epochs, within the interval [ 1 ] m n N ,n    .Accordingly, Equations ( 9) and ( 17) can first be rewritten using the time-forward solution by [29]:  , , ) . Generally, matrix G is set as the identity matrix I , and  is time invariant.The noise vectors n w and n v are both white, and the corresponding variance matrices related to the extended model are , , ) The details of the expanded state-space model can be further explored in reference [23].

UFIR Estimator
From the extended state dynamic Equation ( 24), state n x with respect to x can be represented as: where G .Combining Equations ( 28) and ( 25), the measurement equation becomes: where Given the extended matrices, the batch UFIR estimate n x at n epoch using the discrete convolution is [30]: Equation ( 30) can be further computed in an iterative form, as: ( ) where [ , ] j m a n   and j G and j x are initialized in short batch form as: where a represents the dimension of the state n x .
Furthermore, the error covariance j P can be calculated using the following method: in which denotes the UFIR filter gain that is used to correct the bias.The initial covariance matrix is given by: It is important to note that matrix  , which is used to correct the estimation bias in UFIR, is determined by H and  , which are already known.However, the gain of KF depends not only on the process noise Q but also on the measurement noise R , making it challenging to set them perfectly.This distinctive feature makes UFIR more robust than KF against the errors of noise statistics.It is also noteworthy that UFIR requires the horizon length N to be optimal, which will be given in the next section.

Turning of UFIR
For the UFIR strategy, the objective of opt N is to minimize the MSE.The memory length opt N is the critical tuning parameter that ensures UFIR is optimal.When opt N < N , the UFIR fits with the model well, but the noise reduction is insufficient.Alternatively, the bias errors become dominant.This is because the UFIR no longer fits with the model, although it better suppresses the noise.Generally, if the ground truth is available in the operation, then opt N can be pro- vided by [31]: where tr( )  denotes the trace operator and  argmin  is the minimization operator.
Otherwise, an alternative approach is given as follows: where is the mean-square value.

Experiments and Results
In this section, we conducted land-vehicle experiments and collected data from real scenarios to assess the efficiency of the proposed approach.To be specific, two separate experiments were performed in different areas: one in an urban environment and the other in a suburban environment, utilizing authentic GNSS signal reception.These scenarios were selected based on frequent occurrence challenges in the field to the majority of GNSS applications.By examining these common and representative situations, we aim to ensure that our findings are not limited to niche cases but have broader applicability and relevance across a range of scenarios.Experimental results derived from these scenarios can be generalized to many similar situations in GNSS solutions.Furthermore, the existing LS and KF positioning methods are also employed as benchmarks for a fair comparison and analysis.

Data Description
To verify the proposed method, we utilized datasets collected from two different scenarios.Figure 2 illustrates the experimental setup, in which two GNSS antennas were mounted on the top of the car.A commercial U-blox NEO-M8T GNSS receiver, sourced from u-blox, Thalwil, Switzerland, was employed for providing a ground truth of positioning, at a frequency of 2 Hz.The raw signal was first collected by a GNSS IF front-end sampler and then stored on a laptop.Finally, the raw IF data were processed by a GNSS software-defined receiver (SDR), version 1.0, which is developed based on MATLAB.The SDR processing was conducted using a laptop equipped with an Intel Core i5-12400H processor, running at 4.4 GHz and 16 GB of RAM.The processor is a product of Intel Corporation, headquartered in Santa Clara, CA, USA.Detailed configurations of the equipment are listed in Table 1.

Positioning Performance Evaluation
Figure 3 shows two test trajectories and satellite views plotted on Google Earth.The test scenarios were classified into two types: bumpy road sections (Scene 1) and smooth road sections (Scene 2).During Scene 1, only five satellites were visible due to blocked GNSS signals, and the experiment lasted for 165 s.For the second test, an open area with no obstructions was chosen, and the experiment lasted for 245 s with a total of seven satellites available. . . .
The process and measurement noise covariance matrices, expressed in squared meters, are initialized as follows [27]: where diag [ ]  represents a diagonal matrix; I is an identity matrix of dimension m; and m is the number of the tracking channel.The observation noise covariance matrix is then adaptively computed by Equation (19).The horizon lengths of these two tests are calculated using Equation (39)-specifically,

Test I: Urban Area, Bumpy Road
The positioning results in the horizontal direction for the LS, KF, and UFIR correction are illustrated in Figure 4.The horizon positioning error and elevation positioning error of these three methods are presented in Figures 5 and 6, respectively.The entire driving route is divided into three periods.The car travels relatively smoothly in the first and third periods.As we can see, in the first-and third-time intervals, the errors obtained by the LS algorithm were the most significant, with horizontal positioning errors of approximately 11 m and vertical positioning errors of 24 m.The positioning results of the KF algorithm outperformed those of the LS, reducing the horizontal positioning error to around 8 m and the vertical error to 20 m.When employing the UFIR method, the horizontal and vertical positioning errors were further improved to 5.79 m and 15 m, respectively.The results indicate that the UFIR method is markedly superior to the other two algorithms.However, in the second period, the car encounters a bumpy road.We consider this situation as process interference.The reason is that the tilt of the antenna causes a decrease in GNSS signal gain.In this case, as a result, the positioning errors of all three methods have diverged.Specifically, the maximum horizontal positioning error for LS reached 165 m, while for KF it was 147 m, and for UFIR, it was reduced to 124 m.It is also observed that the LS algorithm yields the largest positioning error.Although the KF algorithm shows some improvement in positioning accuracy, the enhancement is limited, while the positioning results of the UFIR-based GNSS solution method are relatively smooth.The comparative results illustrate that, in terms of position accuracy, the relationship between these filters is LS < KF < UFIR.
Accordingly, the quantitated horizontal position results of LS, KF, and UFIR have been listed in Table 2, including the median, mean values, and standard deviation (STD).It can be seen that, compared with the KF, the mean values of horizon positioning errors from UFIR decreased by 28.33%, and the median values decreased by 30.48%.Additionally, both the KF and UFIR have a common objective of minimizing the discrepancy between the predicted states and the observed states over time, thereby ultimately reducing the error in the estimated states.The residual error, also known as the estimation error or innovation, is a fundamental concept used to evaluate the quality and accuracy of the estimated states.That is to say, analyzing the trajectory of the state elements 0 x offers a straightforward indicator of the filter's capability to precisely anticipate the actual states.Hence, a residual error analysis is conducted as part of the stability analysis.
To achieve this, we introduced the relationship between the residual error vector k e and states k x and its 1-sigma bound, which can be found by [32]:

H x
is the measurement matrix that relates the predicted state Figure 7 presents the trends of the state errors in the x, y, and z directions, along with their corresponding 1-sigma confidence intervals.These state errors were computed along the trajectory represented in Figure 4.In Figure 7, only 1-sigma bound is shown.It is observed that both the KF and UFIR residual errors are within the ±1-sigma bounds.This indicates that the residual error will also satisfy the conventional ±3-sigma criterion.Even so, the UFIR demonstrates a higher degree of stability and a reduced level of noise variance in its residual errors when compared to KF.

Test II: Suburban Area, Smooth Road
Because the first test was carried out for about 165 s and in a very specific scenario, this does not usually happen in real applications.In this test, we consider a generalized one, in which the experiment was conducted on a smooth road and lasted for about 204 s.The automobile kept static for about 35 s and then moved to the southeast with a moderate dynamic along a dyke-dam.The trajectories and positioning errors in the horizon direction are shown in Figures 8 and 9, respectively.From Figure 9, 11.36 m of horizon mean error is obtained by the LS method, with the STD of 7.36 m.The error decreases to 8.95 and 5.30 m, respectively, after using the KF method.KF can achieve slightly better position results than those of LS.This is mainly because KF combines the current measurement data with the previous estimate and optimizes the estimation by predicting the system state and its covariance.In other words, the LS method typically provides a single estimation without considering future state predictions.However, when utilizing the UFIR approach, the mean positioning error is significantly reduced in most epochs, as demonstrated in Figure 10.As expected, the horizon mean error decreases to 6.35 m after using the UFIR method.The STD also decreases significantly to 2.65 m compared with that of the LS and KF methods.These results are consistent with the results of test I. Similarly to the first test, the results show that the UFIR outperforms the LS and KF for both the static and dynamic periods.
Additionally, as mentioned earlier, the primary distinction among the LS, KF, and UFIR approaches lies in the amount of historical information utilized during the iteration.The horizon length means how much the previous measurements are taken into consideration.This also means that the more optimal the opt N , the closer the estimated results are to the optimal value.Therefore, we also presented the positioning error of UFIR with different horizon lengths, as shown in Figure 11.Note that as the horizon length increases, the positioning error continuously decreases.The optimal estimation result in UFIR is achieved with a horizon length of 250.However, if N exceeds 250, the positioning error will start to increase continuously.This is because the trace of error covariance is concave on the horizon length.Naturally, a larger horizon length implies a greater amount of historical data to take into account, resulting in longer computation times.For simplicity, we extracted the first 100 s of GNSS data to evaluate the computational load under different horizon lengths.The computational time of UFIR with different values of N is listed in Table 3.The optimal estimated computation time of UFIR (1808.87s) is 1.18 times that of KF (1531.23 s), while the accuracy has been improved by 34%.This is acceptable because of the abundant onchip computing resources.Additionally, when it comes to memory usage, as illustrated in Table 3, the memory consumption for UFIR escalates with an increase in the horizon length.For a comparative analysis, the KF consumed 95.00 Mb of memory in the computational experiments.When the optimal horizon length (N = 250) is employed, UFIR requires a memory capacity of 162.85 Mb, which is approximately 1.7-fold higher than that of KF.Considering the marked improvement in localization accuracy attained with UFIR, the increased memory requirement is deemed a reasonable compromise.To further elucidate the performance of the UFIR in varying noise conditions, an additional set of experiments was conducted.This extended analysis focuses specifically on the robustness of the UFIR algorithm under dynamic noise environments, thereby offering a more comprehensive understanding of its applicability in real-world scenarios.Specifically, we artificially introduced Student's t-distributed noise with five degrees of freedom into the observed pseudoranges.To simulate a situation where noise is growing over time, we applied varying scale factors ranging from 1 to 100.The growing noise was introduced into the last 30 s of the data.Figure 12 shows the horizontal positioning accuracy of UFIR and KF under increasing noise conditions.As we can see, in the first 70 s, the horizontal positioning error for KF is 8.08 m, while for UFIR it is 5.1 m.However, there is a noticeable increase in the horizontal positioning error when using the KF as the level of noise increased.In contrast, the UFIR method exhibited a relatively stable performance under the same conditions.This behavior can be attributed to the sensitivity of the KF to noise, particularly when it is non-Gaussian.UFIR, on the other hand, is better suited for handling zero-mean noise, displaying a greater ability to suppress it.

Discussion
To validate the efficacy of UFIR, the above two targeted experiments were designed to address different application scenarios.Figures 5, 6, 9, and 10 and Table 2 illustrate the positioning results of LS, KF, and UFIR.Taking test II as an example, the LS method (mean error: 11.36 m) in this study has the worst position results.The reason is that LS relies solely on spatial geometric structures to estimate position and assumes that there is no temporal correlation between states.By incorporating the state transformation constraints and updating the state prediction with current measurements, smoother positioning results are obtained for KF (mean error: 8.95 m).Nevertheless, KF is essentially a Markov process, and the current state estimate is only achieved between two nearest neighbors.This means that historical information is excluded when estimating the current state.Thus, KF produces better estimate results.But, the improvement is limited by its recursive property.For GNSS, a typical time-varying system, the state and measurement are highly correlated in time.Unlike KF, UFIR is iterative, and the amount of historical data used during iteration and invariance to noise uncertainty make it more accurate and robust.As a result, the mean horizon positioning error of UFIR is reduced to 6.35 m, which is significantly better than the other two algorithms.Meanwhile, incorporating historical data will result in an increase in the computational load.Thus, we evaluate the computational time under different horizon lengths N. The results in Table 3 indicate that as N increases, the computational load and memory used also increase, but compared to the improvement in positioning accuracy, it is still acceptable.Further analysis on noise sensitivity also corroborates that the UFIR-based method exhibits superior robustness to zero-mean noise compared to the KF method.These results suggest that by introducing historical information and eliminating reliance on noise, higher accuracy and more robust results can be achieved.
Consequently, the UFIR-based GNSS positioning estimator offers improved accuracy and robustness compared to LS-and KF-based estimators.While it does come with the drawback of increased computational load, this is no longer a significant concern for modern computers and microprocessors.

Conclusions
Currently, GNSS positioning is the primary source for providing a globally referenced position.However, accurate positioning estimation is a great challenge for GNSS receivers due to the harsh application environment.LS and KF are currently the main methods for GNSS positioning estimation.However, in LS, the observations between different epochs are independent of each other, neglecting the time correlation between states.In KF, although the adjacent states are correlated with each other through the state transition matrix, the rich historical data have not been effectively utilized.Given this, the

Horizontal position errors (m)
UFIR is proposed to estimate the GNSS position in this paper.UFIR inherits the advantages of KF and takes historical data into account by expanding the state-space model, which will improve positioning performance.Furthermore, it exhibits superior robustness because it does not rely on noise statistics.Two comparative experiments were performed to validate UFIR's positioning performance.The results show that the proposal demonstrates a significant improvement in positioning accuracy.When compared to the LS and KF methods, the mean positioning error was reduced by 44% and 29%, respectively.UFIR overall outperforms these two methods in both positioning accuracy and robustness, only sacrificing some acceptable computing resources.
the receiver's position and the th i satellite position vectors in an earth-centered and earth-fixed (ECEF) frame, respectively; u t  is the receiver clock offset; i  denotes the pseudorange between the th i satellite and the GNSS receiver; i I and i T are ionospheric delay and tropospheric delay errors; i  represents other error factors, such as pseudorange noises and unmodeled random error; and c denotes the speed of the light.

Figure 3 .
Figure 3. Trajectories and satellites sky-plot of two tests scene.

Figure 4 .
Figure 4. Positioning results of test one in Google map.

Figure 7 .
Figure 7. Residual error (in x, y, z directions) and the corresponding 1-sigma bounds of test 1.

Figure 8 .
Figure 8. Positioning results of test two in Google map.

Figure 11 .
Figure 11.Horizon positioning errors under different N.

Figure 12 .
Figure 12.Horizon positioning accuracy of KF and UFIR under growing noise conditions.

Author Contributions:
Conceptualization, J.D. and B.X.; methodology, J.D.; software, B.X.; formal analysis, J.D.; investigation, B.X.; resources, B.X.; data curation, J.D.; writing-original draft preparation, J.D.; writing-review and editing, B.X.; visualization, J.D.; supervision, B.X.; project administration, J.D.; funding acquisition, J.D. and L.D. All authors have read and agreed to the published version of the manuscript.Funding: This work was supported by the National Natural Science Foundation of China (grant number 60904085); the Jiangsu Funding Program for Excellent Postdoctoral Talent; the Foundation of National Key Laboratory of Transient Physics; and the Foundation of Defence Technology Innovation Special Filed.Data Availability Statement: Not applicable.
the extended vectors.

Table 1 .
Parameter settings of the experimental equipment.

Table 2 .
Statistical analysis of the horizontal position errors.

Table 3 .
Computational load with different horizon lengths N.