1. Introductions
The United States was the first country to start the new generation of space satellite navigation and position system. Subsequently, Russia and China have studied global navigation positioning systems and developed them rapidly. The research and development of GPS was to serve military needs and, as early as 1972, the United States Navy transit project and the Air Force 621B project realized the application of GPS. With the development of economic globalization, smart phones, vehicles and ships, which are based on the combination of GPS and modern technology, have gradually penetrated into people’s daily life [
1]. A country’s important national infrastructure needs a perfect and stable Global Positioning System, and the precision is an important parameter in the development of GPS [
2].
At present, there are two methods to improve the positioning precision of GPS. The first is to use the Differential Global Positioning System (DGPS) [
3]. However, this method has several limitations: (1) the equipment limitations: the receiver of a differential signal must be used in practical applications. (2) Restrictions of work area: the work area is limited by the differential network of a wide area. (3) Non-autonomous: the transmitting source is required to improve precision [
4]. The other is to use the positioning data which are collected by the receiver for error correction. This method has poor real-time performance and is vulnerable to external interference [
5,
6]. Therefore, it is particularly important to improve the positioning precision of GPS without changing the hardware device.
Scholars have built the error model of GPS to improve the positioning precision of GPS. In Ref. [
7], Zhao Shan et al. used the ARMA (3,2) model to obtain the single error model of a user by the combination of position operation and clock displacement filtering, which has a certain stability. In Ref. [
8], Liu Di et al. made a long-term observation of a static point to obtain the same characteristics of error sequence and elevation error sequence between longitude and latitude, and finally established the AR (n) error model of GPS. In Ref. [
9], Wang Rong et al. analyzed the error of GPS and established the error model of GPS. In Ref. [
10], Zhiqiang Liu et al. modeled the time-selective channel as an AR process and used a Kalman filter to track the time change. In Ref. [
11], Tan-Jan Ho et al. proposed a framework modeling which is based on multiple AR models, and developed the channel predictor MAR. In Ref. [
12], Christos Komnanakis et al. used a Kalman filter to track a low-order autoregressive model, which is similar to the change in multiple input and multiple output channels.
The models used AR models based on time series by the scholars. If the output sequence of a model is
, there is a non-linear problem of parameter estimation, and the algorithm is complex and difficult. In order to solve the non-linear problems and effectively improve the positioning precision of GPS, we use the TAS and Kalman filter to analyze the characteristics of positioning error data, set up the Autoregressive (AR) model of GPS, complete the parameter estimation of the AR model by the least square method, and estimate the state of the system by a Kalman filter, so as to make the positioning precision more reliable and accurate [
13].
2. The Error Model of GPS
The error data of GPS are a discrete random variable of time series, which is different from the analyzed dynamic data because the time series is a realization of random processes and has a different physical background [
14,
15]. The processing method of the corresponding model is an approximate description according to the data characteristics, to determine the type of model suitable for the time series [
16].
There are three important models of finite order linear in TSA, Moving Average (MA) model, Autoregressive (AR) model, and Autoregressive Moving Average (ARMA) model [
17]. The type of model can be determined by analyzing the autocorrelation and partial correlation of the error signal, and the identification methods of the three models are shown in
Table 1.
As shown in
Table 1, the autocorrelation function of MA is truncated, while the partial autocorrelation function is trailing. The autocorrelation function and partial autocorrelation function of ARMA are trailing. The autocorrelation function of AR has the trailing property, while the partial autocorrelation function has the truncation property, where truncation means that the time-order autocorrelation function or partial autocorrelation function is 0 when the order is greater than a constant K. The trailing property means that the autocorrelation function or partial autocorrelation function fluctuates near zero after a certain order.
If a time series is generated by a certain type of model, it should theoretically have corresponding statistical characteristics [
18,
19]. Therefore, the sample autocorrelation function and partial autocorrelation function of time series can be calculated, and the characteristics can be compared with the characteristics of a theoretical autocorrelation function and partial autocorrelation function of different types of series, and then the model type suitable for the series can be judged [
20].
The estimated value of the autocorrelation function is a measurement to describe the dependence between values of random processes at different times.
After obtaining the error observation data
, the sample autocorrelation function
of the error sequence is
where the estimate value of autocovariance is
,
. The estimate data of mathematical expectation are
.
Using a numerical method and MATLAB programming, the response simulations of the autocorrelation function and partial autocorrelation function of longitude and latitude are obtained.
By using
instead of
, the estimation of the partial autocorrelation function can be obtained recursively. If the autocorrelation function
of the sample is truncated in step, it can be determined as an AR (M) sequence. If
is not truncated, it is an ARMA sequence. From
Figure 1a and
Figure 2a, the error autocorrelation sequences of longitude and latitude have coordinates fluctuations near zero after order 4, the red line represents the autocorrelation function data of latitude and longitude error. From
Figure 1b and
Figure 2b, the autocorrelation function of longitude and latitude shows trailing, and the partial autocorrelation function shows truncation, the blue line represents the partial correlation function error data of latitude and longitude error, * represents the data point of latitude and longitude error per second. Therefore, the error of longitude and latitude can be expressed by the AR (M) model, that is,
The recursive formula of can be used to find by using instead of .
4. The Application of Kalman Filter in the Error Model
4.1. Discretization of State Equations of Continuous Systems
The actual physical system is generally continuous, and the dynamic characteristics are described by continuous differential equations. Therefore, the discretizations of the system equation and the observation equation are needed [
28,
29].
The system state equation describing the dynamic characteristics of the physical system is:
where the driving source
of the system is the white noise process, which is
where q is v(t) variance intensity matrix.
σ(t-
τ) is a function of Dirac
σ.
According to the linear system theory, the discretization of the system state equation is:
where the one-step transfer matrix a
satisfies the equation:
where
. The equation is the real-time calculation formula of a one-step transfer matrix.
The discretization state equation of continuous system also includes the equivalent discretization of the excited white noise process
.
Equation (11) can be abbreviated as:
where
. Then, for
defined in Equation (13), it is:
where,
is Kronecker
function. The variance matrix
of
satisfies the following equation:
where,
, the Equation (16) is the real-time calculation formula of
.
4.2. The Basic Equation of Discrete Kalman Filter
Kalman is linear minimum variance estimation [
30]. For Kalman model, the state equation and observation equation of discrete linear system are respectively:
where,
is one step transition matrix of time, which is
to
.
is the measurement matrix.
is the system noise driving.
is the observation noise sequence.
is the excitation noise sequence of the system. Both
and
are satisfied
where,
is the variance matrix of the system noise sequence, it is a non-negative matrix.
is the variance matrix of noise sequence on both sides, it is positive definite matrix.
is Kronecker δ function.
If is conformed to Equation (19), the measured value is conformed to Equation (20), the system noise and the measurement noise are conformed to Equation (21), and the system noise variance matrix is non-negative definite, the measurement noise variance matrix is positive definite, and the measurement of time K is , the equation can be solved for estimation .
Then
where,
is the one-step prediction equation of state.
is the one-step prediction of mean square error.
is the gain equation of filter.
is the estimation equation of filter.
is the optimal mean square error at time K.
As long as the initial value and is given, according to the measurement at time k, the estimated state at time k can be deduced.
4.3. Kalman Filter Based on AR Model
The state vector of GPS is
The model is based on , the state space model of Kalman filter model is Equations (9) and (10), where and have statistical characteristics of Equation (21).
The process noise turn the difference equation of AR (10) model into state Equation (24), it is
Let
, Equation (24) can also be abbreviated as
The observation equation is
where,
.
Because the statistical characteristics of
and
are consistent with Kalman filter, the recursive expression of Kalman filter based on AR model is show as:
where,
is the one-step prediction equation of filter state.
is the one-step prediction of mean square error.
is the gain equation of filter. If the Kalman filter gain value
is very small, the filtering result is closer to the recursive result of the system state estimation value. If
is large, the filtering result is closer to the state variable calculated of the observed value.
is the estimation equation of filter.
is the optimal mean square error at time K.
There are two ways to deal with Q: one is that Q is a certain value. The second is that Q is an uncertain random variable. Therefore, the Q value in this paper is a definite value, which is the process noise variance of the system, and its value is . When the state transition process has been determined, the smaller of Q is the better. When Q gradually increases, the convergence of filter slows down and the disturbance of the state variable becomes larger. The value of R is related to the characteristics of the device and is the input value of the filter.
6. Conclusions
In this study, the error of GPS is researched, and the statistical characteristics of GPS are analyzed and simulated. The least square method is applied to estimate the parameters, and obtain the basic equation of a discrete Kalman filter by the continuous Kalman filter. To eliminate the random error of GPS dynamic positioning data, the error model of GPS is combined with a Kalman filter and the experimental results show that the smaller the mean square error, the better the precision of the prediction model on the experimental data. The dynamic positioning precision after correction has been significantly improved, and this method can effectively improve the positioning precision of GPS and support for the application of GPS.