Next Article in Journal
Multi-Mode Electromagnetic Ultrasonic Lamb Wave Tomography Imaging for Variable-Depth Defects in Metal Plates
Previous Article in Journal
Implementation Strategies for a Universal Acquisition and Tracking Channel Applied to Real GNSS Signals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tracking Architecture Based on Dual-Filter with State Feedback and Its Application in Ultra-Tight GPS/INS Integration

School of Automation, Beijing Institute of Technology (BIT), Beijing 100081, China
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(5), 627; https://doi.org/10.3390/s16050627
Submission received: 28 January 2016 / Revised: 24 April 2016 / Accepted: 26 April 2016 / Published: 2 May 2016
(This article belongs to the Section Physical Sensors)

Abstract

:
If a Kalman Filter (KF) is applied to Global Positioning System (GPS) baseband signal preprocessing, the estimates of signal phase and frequency can have low variance, even in highly dynamic situations. This paper presents a novel preprocessing scheme based on a dual-filter structure. Compared with the traditional model utilizing a single KF, this structure avoids carrier tracking being subjected to code tracking errors. Meanwhile, as the loop filters are completely removed, state feedback values are adopted to generate local carrier and code. Although local carrier frequency has a wide fluctuation, the accuracy of Doppler shift estimation is improved. In the ultra-tight GPS/Inertial Navigation System (INS) integration, the carrier frequency derived from the external navigation information is not viewed as the local carrier frequency directly. That facilitates retaining the design principle of state feedback. However, under harsh conditions, the GPS outputs may still bear large errors which can destroy the estimation of INS errors. Thus, an innovative integrated navigation filter is constructed by modeling the non-negligible errors in the estimated Doppler shifts, to ensure INS is properly calibrated. Finally, field test and semi-physical simulation based on telemetered missile trajectory validate the effectiveness of methods proposed in this paper.

1. Introduction

The navigation errors of GPS are not divergent with time, but a receiver usually performs less well in highly dynamic situations because there are external signals to be dealt with. On the other hand, though INS has the problem of error accumulation, it realizes autonomous navigation and is robust to vehicle dynamics. Since GPS and INS have complementary advantages, their integration can overcome each other's defects and provide a better navigation performance [1]. Unlike loose or tight GPS/INS integration, in which GPS is just an aid to INS, the ultra-tight GPS/INS integration uses INS to couple with GPS tracking channels. In this sense, the ultra-tight form is especially applicable to high dynamic situations [2], and has gradually become the mainstream approach for designing GPS/INS integrated systems.
The ultra-tight GPS/INS integration originates from vector tracking [3]. Its core concept is using the known ephemeris and the corrected inertial navigation information to calculate tracking parameters. Namely, satellite signal locking is assisted by the integrated navigation filter, and GPS tracking channels no longer remain independent of each other. According to the structural design of the filter models, the ultra-tight GPS/INS integration methods can be roughly divided into two types. In one case, a single high-dimensional filter is used to estimate the GPS tracking parameters of each channel and navigation parameters together [4,5], and in the other case, each channel has a baseband signal pre-filter which estimates the tracking parameters for itself, and the integrated navigation filter just needs to estimate the navigation parameters. The second category is actually a form of federated filtering, and usually has higher computational efficiency [6,7,8].
The baseband signal preprocessing realizes optimal estimation of signal characteristic quantities by KF. Thus, under the application background of the ultra-tight GPS/INS integration, it can restrain phase noise more effectively than the loop filter of the classical GPS tracking channel [9,10]. For each channel, the baseband signal preprocessing is generally completed by a single pre-filter [11]. Both carrier tracking errors and code tracking errors are included in the state vector, namely, these two kinds of tracking errors are tightly coupled through the state and observation equations. However, the precision of carrier and code loops do not have the same order of magnitude. The code loop has great tolerance for tracking errors due to the long wavelength of CA code, while the carrier loop is much more sensitive [12,13]. As preprocessing continues, the code tracking errors would be delivered to the carrier loop and cause the degradation of carrier tracking precision or even loss of lock. In addition, even if the baseband signal preprocessing is applied in a GPS receiver, the Doppler shifts estimated under highly dynamic conditions may still contain significant errors. In the ultra-tight GPS/INS integration, these incorrect GPS outputs eventually lead the estimation of INS errors to be contaminated [14]. That generally means the whole system suffers performance degradation, or even completely crashes.
To solve these problems, this paper makes creative contributions in both the tracking and navigation domain. Two independent pre-filters with state feedback are constructed to replace the classical loop filters. The advantages of decoupling between carrier tracking and code tracking are analyzed in detail. Since the single-filter and dual-filter preprocessing models are constructed based on the same principle, their comparison can be free from the influence of other factors. In fact, compared with a popular error-state pre-filter, even the single-filter model has better tracking performance due to the benefits of state feedback. In the ultra-tight GPS/INS integration, local carrier generation is still controlled by state feedback, and the carrier frequency derived from the external navigation information is just used to correct the state feedback value. Meanwhile, a specific integrated navigation filter is presented by taking the pseudorange rate estimation error in each channel as the additional state variable. This scheme ensures the estimation of original state variables, especially the INS errors, is immune to these GPS tracking errors, and hence enables the integrated system to maintain high accuracy even when GPS receiver operates under harsh conditions.
The rest of this paper is organized as follows: first, a baseband signal preprocessing model is designed based on a dual-filter structure, while the state feedback control values of carrier and code loops are derived in detail. Then, the specific integrated navigation filter is proposed, and a concrete approach to aid the tracking channel is elaborated. Furthermore, the results of field tests and semi-physical simulation are given and discussed. Finally, the work in this paper is summarized.

2. Construction of Baseband Signal Preprocessing Model Based on Dual-filter Structure

The classical GPS tracking channel, which typically adopts second-order or third-order phase control technology, has fixed loop bandwidth. Although a narrower bandwidth is helpful to reduce the thermal noise, in highly dynamic situations, the tracking channel may not retain useful high frequency information and suffers from significant dynamic stress error. These usually lead to signal distortion and tracking loss. Thus, a trade-off between noise resistance and response to dynamics is required while setting the loop bandwidth. In view of the abovementioned fact, the classical control approach cannot cope well with high dynamics, and hence KF, which belongs to the modern control category, is taken into account. KF is a kind of optimal estimator, which allows precise modeling of the tracking loop. With its help, baseband signal preprocessing can achieve accurate control of the local signal generation, even in highly dynamic situations [15].
The baseband signal preprocessing model constructed in this paper contains two pre-filters. In this model, the estimation of carrier and code tracking errors are relatively separate, the same as the manner adopted by classical GPS receiver. Figure 1 just takes carrier loop for example to show the details of this model. It is observed that the loop filter in the classical tracking channel has been replaced by a pre-filter equipped with state feedback. Obviously, designing a suitable KF is the key to accurate carrier tracking.
This paper sets X k c a r r = [ Δ θ k c a r r   f k c a r r   f ˙ k c a r r ] T as the state vector of pre-filter 1, where k is the epoch, Δ θ k c a r r is the carrier phase error (unit: rad), f k c a r r is the input carrier frequency, and f ˙ k c a r r is the changing rate of the input carrier frequency. These state variables have the following recurrence relations:
Δ θ k , k 1 c a r r = 1 T c o h 0 T c o h ( Δ θ k 1 c a r r + 2 π f k 1 c a r r t + π f ˙ k 1 c a r r t 2 ) d t 1 T c o h 0 T c o h 2 π u k 1 c a r r t d t = Δ θ k 1 c a r r + π f k 1 c a r r T c o h + 1 3 π f ˙ k 1 c a r r T c o h 2 π u k 1 c a r r T c o h
where T c o h is the coherent integration time, u k 1 c a r r is the local carrier frequency at epoch k 1 .
The output of carrier phase discriminator at epoch k, notated as Δ θ k , k 1 c a r r , is the average carrier phase error over an integration interval, due to the presence of Doppler shift. Thus, this output can be expanded as:
Δ θ k , k 1 c a r r = 1 T c o h 0 T c o h ( Δ θ k 1 c a r r + 2 π f k 1 c a r r t + π f ˙ k 1 c a r r t 2 ) d t 1 T c o h 0 T c o h 2 π u k 1 c a r r t d t = Δ θ k 1 c a r r + π f k 1 c a r r T c o h + 1 3 π f ˙ k 1 c a r r T c o h 2 π u k 1 c a r r T c o h
Unfortunately, if Equation (2) is directly used to generate the observation model for pre-filter 1, there will be a time lag between the state vector and observation [16]. Hence, by moving Δ θ k 1 c a r r , f k 1 c a r r and f ˙ k 1 c a r r to the left side of the equal sign, Equation (1) is rewritten as follows (more details are available in Appendix A):
{ Δ θ k 1 c a r r = Δ θ k c a r r 2 π f k c a r r T c o h + π f ˙ k c a r r T c o h 2 + 2 π u k 1 c a r r T c o h f k 1 c a r r = f k c a r r f ˙ k c a r r T c o h f ˙ k 1 c a r r = f ˙ k c a r r
Substitution of Equation (3) into Equation (2) results in:
Δ θ k , k 1 c a r r = Δ θ k c a r r π f k c a r r T c o h + 1 3 π f ˙ k c a r r T c o h 2 + π u k 1 c a r r T c o h
Then, in terms of Equations (1) and (4), the state equation and observation equation of pre-filter 1 can be expressed as:
{ X k c a r r = Φ k , k 1 c a r r X k 1 c a r r + M k 1 c a r r + W k 1 c a r r Z k c a r r = H k c a r r X k c a r r + N k c a r r + V k c a r r
where:
Φ k , k 1 c a r r = [ 1 2 π T c o h π T c o h 2 0 1 T c o h 0 0 1 ]
M k 1 c a r r = [ 2 π T c o h u k 1 c a r r 0 0 ] T
H k c a r r = [ 1 π T c o h 1 3 π T c o h 2 ]
N k c a r r = π T c o h u k 1 c a r r
Z k c a r r is the observation, and can be directly obtained from carrier phase discriminator, namely:
Z k c a r r = Δ θ k , k 1 c a r r = arctan Q p I p
W k 1 c a r r , V k c a r r represent the system noise and observation noise, respectively. Assume these two noises have the following characteristics:
{ E [ W k ] = 0 , E [ W k W j T ] = Q k c a r r δ k j E [ V k ] = 0 , E [ V k V j T ] = R k c a r r δ k j E [ W k V j T ] = 0
where δ k j is the Kronecker delta. Then, the corresponding iterative algorithm for pre-filter 1 is given below:
X ^ k , k 1 c a r r = Φ k , k 1 c a r r X ^ k 1 c a r r + M ^ k 1 c a r r
P k , k 1 c a r r = Φ k , k 1 c a r r P k 1 c a r r Φ k , k 1 c a r r , T + Q k 1 c a r r
K k c a r r = P k , k 1 c a r r H k c a r r , T ( H k c a r r P k , k 1 c a r r H k c a r r , T + R k c a r r ) 1
P k c a r r = ( I K k c a r r H k c a r r ) P k , k 1 c a r r
X ^ k c a r r = X ^ k , k 1 c a r r + K k c a r r ( Z k c a r r N k c a r r H k c a r r X ^ k , k 1 c a r r )
The carrier loop is in essence a Phase Lock Loop (PLL), and aims to achieve the elimination of the carrier phase error. Since Δ θ k + 1 c a r r = Δ θ k c a r r + 2 π ( f k c a r r u k c a r r ) T c o h + π f ˙ k c a r r T c o h 2 , by assuming Δ θ k + 1 c a r r = 0 , the feedback at epoch k can be obtained as:
u k c a r r = Δ θ ^ k c a r r 2 π T c o h + f ^ k c a r r + 1 2 f ˙ ^ k c a r r T c o h
It can be seen that u k c a r r is not exactly equivalent to the estimated input carrier frequency f ^ k c a r r . And further discussions about u k c a r r and f ^ k c a r r can be found in Section 3 and Section 4.1.
Similarly, the state vector of pre-filter 2 is set to X k c ode = [ Δ θ k c o d e   f k c o d e   f ˙ k c o d e ] T , where Δ θ k c o d e   is the code phase error (unit: chip), f k c o d e is the input code frequency, and f ˙ k c o d e is the changing rate of the input code frequency. Meanwhile, in consideration of the following formula:
f d = f c a r r f I F = f L 1 f C A ( f c o d e f C A )
where f d is the Doppler shift, f I F is the Intermediate Frequency (IF), f L 1 (1575.42 MHz) is the L1 carrier frequency, f C A (1.023 MHz) is the nominal code frequency, the state equation and observation equation of pre-filter 2 can be written as:
{ X k c o d e = Φ k , k 1 c o d e X k 1 c o d e + M k 1 c o d e + W k 1 c o d e Z k c o d e = H k c o d e X k c o d e + N k c o d e + V k c o d e
where:
Φ k , k 1 c o d e = [ 1 T c o h 1 2 T c o h 2 0 1 T c o h 0 0 1 ]
M k 1 c o d e = [ T c o h u k 1 c o d e 0 0 ] T
H k c o d e = [ 1 1 2 T c o h 1 6 T c o h 2 0 f L 1 f C A 0 ]
N k c o d e = [ 1 2 T c o h u k 1 c o d e 0 ] T
Z k c o d e = [ Δ θ k , k 1 c o d e f k c a r r f I F + f L 1 ]
W k 1 c ode , V k c ode are respectively the system and observation noises of pre-filter 2, while the observation Δ θ k , k 1 c o d e obtained from code phase discriminator can be denoted by:
Δ θ k , k 1 c o d e = 1 2 I E 2 + Q E 2 I L 2 + Q L 2 I E 2 + Q E 2 + I L 2 + Q L 2
u k 1 c o d e is the local code frequency at epoch k 1 . Moreover, at epoch k , the local code frequency, namely the feedback for code tracking, can be expressed as:
u k c o d e = Δ θ ^ k c o d e T c o h + f ^ k c o d e + 1 2 f ˙ ^ k c o d e T c o h
Since code estimation is processed after carrier stripping, its errors have no effect on carrier tracking. At the same time, once these two pre-filters are skillfully tuned, they can accurately track GPS signals in high dynamic situations. Taking pre-filter 1 for example, its equivalent bandwidth can be calculated from the elements of the gain matrix K k c a r r and is proportional to Q k c a r r ( 1 , 1 ) , Q k c a r r ( 2 , 2 ) and Q k c a r r ( 3 , 3 ) [17]. Herein, Q k c a r r ( 1 , 1 ) , Q k c a r r ( 2 , 2 ) and Q k c a r r ( 3 , 3 ) are the diagonal elements of Q k c a r r . Q c a r r ( 1 , 1 ) and Q c a r r ( 2 , 2 ) (the epoch k is omitted to reveal these two values are constant) can be obtained from the receiver clock model [11], while Q k c a r r ( 3 , 3 ) is related to the vehicle dynamics. If Q k c a r r ( 3 , 3 ) is estimated in real-time by any efficient adaptive algorithm, the loop bandwidth will vary in a time-varying optimal manner. This paper uses another more practical method: Q c a r r ( 3 , 3 ) is set to a large constant value based on the general knowledge of vehicle dynamics. Since Q c a r r ( 1 , 1 ) and Q c a r r ( 2 , 2 ) are unchanged, the low variances of Δ θ ^ k c a r r and f ^ k c a r r can be maintained [3]. In other words, although the variance of f ˙ ^ c a r r increases, phase and frequency estimations which normally deserve concern can hardly be affected. Additionally, R c a r r is set according to the knowledge of the phase noise variance, while P 0 c a r r is usually larger than Q c a r r to leave space for state corrections. In this case, pre-filter 1 gradually settles down to a wide bandwidth, and hence is robust to high dynamics. Benefiting from the above mentioned advantages, this baseband signal preprocessing model can be appropriate for highly dynamic situations and make GPS have better tracking and navigation performance.

3. Realization of the Specific Ultra-tight GPS/INS Integration

The carrier tracking loop employing pre-filter 1 is in essence equivalent to a third-order loop. However, a N-th order phase locked loop cannot correctly track the signal whose phase varies by time to the N th or higher-order power (the estimates of signal parameters can still have low variance, but not low bias) [18]. That indicates, even improved by baseband signal processing, an independent GPS receiver is vulnerable to extreme high dynamics, and may give inaccurate tracking and navigation results.
On the contrary, INS is robust to vehicle dynamics. Based on this fact, INS can be used to assist GPS, which usually means baseband signal processing is connected closely with the estimation of navigation parameters. In the ultra-tight GPS/INS integration shown in Figure 2, the integrated navigation filter is responsible for obtaining corrected inertial navigation information, and delivering them back to GPS tracking channels. Because the feedback from the corrected inertial navigation information has the advantage of high accuracy in a short time, GPS tracking channel can perform better under highly dynamic conditions. The detailed work procedure of this ultra-tight GPS/INS integration is as follows: in each GPS tracking channel, in-phase and quadrature signals are processed by phase discriminators and pre-filters to estimate tracking parameters, such as carrier frequency and code phase; the observations of the integrated navigation filter can be converted from tracking parameters as well as deduced from inertial navigation information and satellite ephemeris, and then the corrected inertial navigation information can be obtained; finally, auxiliary tracking parameters, carrier frequency, for example, can be deduced from these external navigation information and satellite ephemeris, and then be applied to assist signal tracking. Through the above analysis, it can be seen that identifying conversion relations between tracking parameters and navigation parameters as well as finding a rational way to assist tracking are the keys to designing a good ultra-tight GPS/INS integration.
In the ideal situation, pseudorange rate and variation of pseudorange output by GPS tracking channel can be respectively derived from:
ρ ˙ G = f d c f L 1
Δ ρ G = Δ θ c o d e c f C A
where c is the speed of light. But in fact, under highly dynamic conditions, the carrier phase of satellite signal may be proportional to time to the power of N or higher value, and hence the estimated carrier frequency may contain a large error δ f c a r r . Since the incorrect GPS outputs could destroy the estimation of INS errors, it is necessary to analyze this error in detail. Therefore Equation (27) is rewritten as:
ρ ˙ G = ( f ^ c a r r f I F ) c f L 1 = ( f d + δ f c a r r ) c f L 1
On the other hand, Equation (28) can remain unchanged due to the fact that the code loop dynamic stress error is usually small. ρ G , which represents the current pseudorange output by the GPS tracking channel, can be calculated from Δ ρ G and the previous pseudorange. Then, in Earth Centered Earth Fixed (ECEF) coordinate system, the observation equation of integrated navigation filter can be expressed as:
{ δ ρ = ρ I ρ G = e 1 δ x + e 2 δ y + e 3 δ z δ t u + v ρ δ ρ ˙ = ρ ˙ I ρ ˙ G = e 1 δ x ˙ + e 2 δ y ˙ + e 3 δ z ˙ δ t r u + δ v G + v ρ ˙
where ρ I , ρ ˙ I respectively represent the pseudorange and pseudorange rate deduced from inertial navigation information and satellite ephemeris, δ x , δ y , δ z , δ x ˙ , δ y ˙ , δ z ˙ are vehicle position errors and velocity errors in ECEF coordinate system, e 1 , e 2 , e 3 are unit vectors in the x , y , and z direction, δ t u , δ t r u represent clock bias and drift respectively, δ v G is the error of pseudorange rate output by the GPS tracking channel and equals to δ f c a r r c / f L 1 , v ρ , v ρ ˙ are observation noises. It is noted that Equation (30) should usually be transformed into the geodetic coordinate system to facilitate INS error correction. Then, taking all the available tracking channels into consideration, the observation matrix H ˜ can be obtained naturally, where the superscript " ~ " denotes the navigation filter designed in this paper.
If there are m satellites being locked by tracking channels, the state vector of the integrated navigation filter can be expressed as:
X ˜ = [ X δ v G , 1 δ v G , 2 δ v G , m ] T
where X is the transpose of the state vector adopted by a traditional integrated navigation filter, namely
X = [ φ E φ N φ U δ v E δ v N δ v U δ L δ λ δ h ε b x ε b y ε b z b x b y b z δ t u δ t r u ]
and the first 15 variables of X represent three INS error states each in attitude, velocity, position, gyro bias, accelerometer bias. Accordingly, the state transition matrix Φ and the noise matrix G of the traditional integrated navigation filter can be expanded to:
Φ ˜ = [ Φ 17 × 17 0 17 × m 0 m × 17 I m × m ]
G ˜ = [ G 17 × 8 0 17 × m 0 m × 8 I m × m ]
In addition, H ˜ can also be written as:
H ˜ = [ H 2 m × 17 [ 0 m × m I m × m ] ]
where H is the traditional observation matrix. The specific information of Φ , G and H can be easily found in many references [19,20,21,22] and will not be discussed in detail here.
According to the conversion relations indicated by Equations (27) and (28), the corrected inertial navigation information can in turn feed back into GPS tracking channels. The corresponding carrier frequency estimate is usually accurate enough, and can be used to directly control local carrier generation. However, in the new baseband signal preprocessing model, it is hoped that the phase error can return to 0 in the shortest amount of time, so the external information assists in local carrier generation through correcting the state feedback value, and the input carrier frequency estimated output by pre-filter 1 can be expected to remain relatively stable and accurate. In other words, instead of being given by Equation (17), the local carrier frequency u k c a r r is calculated from:
u k c a r r = Δ θ ^ k c a r r 2 π T c o h + { f L 1 c [ e 1 ( v x s v x ) + e 2 ( v y s v y ) + e 3 ( v z s v z ) ] + f I F } | k + 1 2 f ˙ ^ k c a r r T c o h
where v x , v y , v z are the corrected vehicle velocity components in ECEF coordinate system, v x s , v y s , v z s are the satellite velocity components. Meanwhile, because vehicle dynamics have less effect on the code loop, generation of local code is still based upon the state feedback value u c o d e given by Equation (26), namely, the structure of the code loop is unchanged.

4. Results and Discussion

This paper carried out static field tests and semi-physical simulations based on a telemetered missile trajectory to verify the validity of the proposed methods. In the static field test, the integration of GPS and INS is not involved, and only GPS is utilized to compare the performance of the baseband signal preprocessing models based on signal-filter and dual-filter structures. In the simulation, the navigation results of the ultra-tight GPS/INS integration are analyzed, and special attention is given to the further improvement of carrier frequency tracking accuracy.

4.1. Static Field Test

Although the baseband signal preprocessing model and the ultra-tight GPS/INS integration can both restrain noise better, the field test is carried out in an open-sky area ( N : 39 ° 57 ' 36.16 " , E : 116 ° 19 ' 1.58 " , Hgt : 32.9 m ) and does not involve the weak signal problem. After collecting signals with a radio frequency front-end, seven satellites are chosen according to their high Carrier-to-Noise Ratio ( C / N 0 ) listed in Table 1. The sky plot of these satellites is shown in Figure 3, and the elevation of each satellite is always greater than 10 ° during the test.
Since weak signals are beyond the scope of this contribution, there is no need to employ INS to aid GPS tracking in static field tests. Hence this part just compares the performance of GPS receivers equipped with single pre-filter (pre-filter-S) and double pre-filters (pre-filter-D). First of all, in the tracking domain, the reference Doppler shifts can be provided by cubic spline fit interpolation. Then, taking PRN25 for example, Figure 4 shows the difference between Doppler shift tracking errors obtained by these two aforesaid methods. According to the statistical analysis, the Root Mean Square (RMS) value of the Doppler shift tracking errors corresponding to pre-filter-S is 0.487 Hz. For pre-filter-D, this value is 0.198 Hz. Obviously, pre-filter-S has a bad Doppler shift estimation outcome due to the influence of code phase error, while pre-filter-D achieves a significantly higher degree of precision. The decoupling between carrier tracking and code tracking proposes a 59.3 percent reduction in Doppler shift tracking error. Staying with the PRN25 example, Figure 5 shows the difference between the code phase tracking errors obtained by these two methods. It is observed that the new model proposed in this paper has slightly better performance. Specifically, the RMS values of code phase tracking errors corresponding to pre-filter-S and pre-filter-D are respectively 1.55   × 10 2 chip and 1.43   × 10 2 chip. The decoupling reduces the code phase tracking error by 7.7 percent.
To analyze the impact of tracking parameters on position and velocity detection, receiver implements navigation through Least Squares (LS). The position estimation errors in the ECEF coordinate system are shown in Figure 6, while their corresponding RMS values are listed in Table 2. It can be seen that pre-filter-D makes positioning more precise. Then, Figure 7 illustrates the velocity estimation errors in ECEF frame, and Table 3 lists their RMS values. Obviously, there are multiple glitch impulses in the velocity estimation result output by the receiver employing pre-filter-S, which is quite similar to the situation in Doppler shift tracking. Sometimes the velocity errors in the y and z direction are even close to 0.3 m/s. On the contrary, the other receiver achieves comparatively stable and accurate outcomes. The absolute maximum values of the velocity errors are smaller than 0.1, 0.13 and 0.12 m/s, respectively.
It can be drawn from the analysis above that, using the new preprocessing model based on dual-filter makes GPS receiver behave better in the static field tests. Furthermore, it is worth mentioning that the frequency difference between input and local carrier is shown in Figure 8. Although the GPS tracking channel generates a local carrier with a high fluctuation from the input carrier, it can be interpreted as resulting from the design principle that Δ θ k + 1 c a r r must return to 0 in the shortest amount of time. The analysis of Figure 4 indicates the precision of Doppler shift tracking is not reduced, that is, pre-filter 1 is insensitive to the wide fluctuation of local carrier frequency.

4.2. Semi-Physical Simulation Based on Telemetered Missile Trajectory

Highly dynamic conditions, which are needed for validating the effectiveness of the ultra-tight GPS/INS integration, are extremely difficult to achieve in practice. To simulate this condition, in this paper, a missile trajectory lasting 180.23 s is derived from Fiber Optic Gyroscope (FOG) and accelerometer data, which are telemetered during a missile test, and is treated as the reference trajectory. The position of the missile relative to the launch site in geodetic coordinate system, and the velocity of the missile in ECEF coordinate system are shown in Figure 9. Figure 10 shows the (strapdown) INS used in the missile test. The statistical characteristics of the FOG and accelerometer are listed in Table 4. If the same errors in turn are added to the original FOG and accelerometer data, a set of new data can be generated. The trajectory corresponding to these new inertial measurements can be used as the INS navigation result, which needs to be calibrated by the ultra-tight GPS/INS integration. In addition, based on the reference trajectory as well as the satellite configuration at the missile launch site shown in Figure 11, the IF signal can be simulated.
Figure 12 provides a close look at the absolute acceleration of the missile following the reference trajectory. It is observed that this value does not remain unchanged during the missile ascent phase. Indeed, this value has abrupt changes at about 4 s and 9 s. That usually means that even baseband signal preprocessing cannot ensure the carrier frequency estimation accuracy. Taking PRN5 as example, the Doppler shift tracking error obtained from the independent GPS receiver equipped with pre-filter-D is shown in Figure 13. Obviously, the rapid variation of acceleration introduces a significant error into the estimated Doppler shift after 4 s, and the sudden reversal of jerk at 9 s causes the sign of this error to change. As a comparison, the same dynamic simulation is also performed with the ultra-tight GPS/INS integration. The local carrier generation can be directly controlled by the frequency derived from the external navigation information (GPS/INS-f), but even so, the carrier tracking loop is still equivalent to a third-order phase locked loop. Hence, in the integrated system, the Doppler shift tracking error follows the same trend as the absolute acceleration. However, the absolute maximum value of this error drops to 5 Hz from 15 Hz, and the whole tracking process goes more smoothly. According to Equation (36), the external navigation information can also be used to correct the state feedback value. In this case, local carrier generation is controlled by state feedback (GPS/INS-u), and the absolute maximum value of the Doppler shift estimation error is 3 Hz. The following discussions focus on GPS/INS-u. It can be seen from Equation (29) that δ v G is proportional to the estimation error of input carrier frequency, so Figure 13 also shows that the estimation of δ v G , 1 (PRN5 is locked by channel #1) is accurate. In this sense, the proposed navigation filter can ensure INS is properly calibrated. If the traditional 17-dimensional filter is applied, in fact, the significant Doppler shift tracking error caused by missile ascent will lead the navigation results to gradually deviate from the reference trajectory until the whole integrated system cannot work.
Since there is no loss of lock in dynamic simulation, the independent GPS receiver equipped with the new preprocessing model can implement navigation. However, in navigation domain, if KF is applied, the large error contained in the estimated Doppler shift may lead to inaccurate positioning results. Hence the independent GPS receiver still uses LS to obtain navigation results. Figure 14 and Figure 15 respectively show its position and velocity estimation errors in ECEF coordinate system. And Figure 16 and Figure 17 show the INS navigation errors in the same coordinate system. Meanwhile, navigation errors of the ultra-tight GPS/INS integration, namely the difference between the corrected inertial navigation information and reference values are shown in Figure 18 and Figure 19. Obviously, the velocity components estimated by the independent GPS receiver have large errors, which are on the order of 3 m/s, during the missile ascent phase. This is consistent with the carrier tracking result. For the case of INS, there is no sudden change in navigation errors. That is because INS is robust to vehicle dynamics. However, INS has the problem of error accumulation. At 180 s, the position errors are 76.3, −21.1 and −60.4 m, while the velocity errors are 0.6, 0.2 and −0.5 m/s. On the other hand, employing the navigation filter with expanded dimension, the ultra-tight GPS/INS integration can provide more accurate navigation results. The absolute maximum values of position errors and velocity errors are smaller than 0.5 m and 0.1 m/s, respectively.
As mentioned before, the traditional navigation filter cannot maintain the estimation accuracy of INS errors when the GPS receiver operates under harsh conditions. Once the incorrect navigation information is utilized to generate local carrier frequency, the Doppler shift tracking error will further increase, so it can be concluded that the proposed navigation filter elegantly avoids this vicious circle. Here, a tight GPS/INS integration based on the traditional navigation filter is realized to show the harmful effect of carrier tracking errors in detail. The corresponding navigation errors are shown in Figure 20 and Figure 21. It can be observed that the tight integration gradually achieves acceptable navigation accuracy as time goes on, but during the missile ascent phase, the velocity errors of this integration are on the order of 2 m/s, and are much larger than those of INS.

5. Conclusions

This paper proposes a baseband signal preprocessing model, which has two pre-filters and uses state feedback values to control the generation of local carrier and code. Compared with the conventional approach based on a single-filter structure, this model limits the impact of code phase error on carrier tracking, and hence enables the receiver to achieve better navigation precision. Furthermore, this model can ensure the stationarity of the estimation of input carrier frequency, even in the case of local carrier frequencies with wide fluctuation. In fact, compared with a popular error-state pre-filter, even pre-filter-S proposed in this paper performs better in the tracking domain due to the benefits of state feedback.
To further promote the tracking and navigation performance, a specific ultra-tight GPS/INS integration is designed on the basis of the new preprocessing model. Since the pseudorange rate estimation error in each channel has been treated as the additional state variable, the proposed navigation filter ensures the estimation of INS errors will not be contaminated. Moreover, in this integration, although the corrected inertial navigation information is used to estimate the carrier frequency, the tracking channel still aims to make the phase error return to 0 in the shortest amount of time. The results of the semi-physical simulation based on a telemetered missile trajectory indicate that navigation solutions with higher accuracy can be achieved. Meanwhile, because tracking is assisted by external navigation information, the error contained in the estimated Doppler shift can be reduced substantially.

Acknowledgments

The authors wish to acknowledge the financial support of the National Natural Science Foundation of China (Grant No. 61153002).

Author Contributions

Xi Zhang provided insights in formulating the ideas and wrote the paper. He also performed the main part of the experiments. Lingjuan Miao offered the experimental data telemetered in a missile test, and critically reviewed the paper. Haijun Shao contributed greatly to experimental data evaluation. All authors read, revised and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

First, the items in Equation (1) are numbered as follows:
Δ θ k c a r r = Δ θ k 1 c a r r + 2 π f k 1 c a r r T c o h + π f ˙ k 1 c a r r T c o h 2 2 π u k 1 c a r r T c o h
f k c a r r = f k 1 c a r r + f ˙ k 1 c a r r T c o h
f ˙ k c a r r = f ˙ k 1 c a r r
Obviously, Equation (A3) can be rewritten as:
f ˙ k 1 c a r r = f ˙ k c a r r
Substitution of Equation (A4) into Equation (A2) results in:
f k c a r r = f k 1 c a r r + f ˙ k c a r r T c o h
then:
f k 1 c a r r = f k c a r r f ˙ k c a r r T c o h
and substitution of Equations (A4) and (A5) into Equation (A1) results in:
Δ θ k c a r r = Δ θ k 1 c a r r + 2 π ( f k c a r r f ˙ k c a r r T c o h ) T c o h + π f ˙ k c a r r T c o h 2 2 π u k 1 c a r r T c o h = Δ θ k 1 c a r r + 2 π f k c a r r T c o h π f ˙ k c a r r T c o h 2 2 π u k 1 c a r r T c o h
then:
Δ θ k 1 c a r r = Δ θ k c a r r 2 π f k c a r r T c o h + π f ˙ k c a r r T c o h 2 + 2 π u k 1 c a r r T c o h
Finally, Equation (3) can be obtained by combining Equations (A4)–(A6) together.

Appendix B

The filter, whose state vector is [ Δ θ k c a r r   Δ f k c a r r   Δ f ˙ k c a r r   Δ θ k c o d e ] T , is a popular baseband signal preprocessing model. Herein, Δ f k c a r r = f k c a r r u k c a r r , Δ f ˙ k c a r r = f ˙ k c a r r 0 . To focus on analyzing the benefits of decoupling between carrier tracking and code tracking, this filer should be further improved. First, consider the equations derived from Equation (18):
f c o d e = β f c a r r + β ( f L 1 f I F )
f ˙ c o d e = β f ˙ c a r r
where β = f C A f L 1 . The new state variables Δ θ k c a r r , f k c a r r , f ˙ k c a r r and Δ θ k c o d e have the following recurrence relations:
{ Δ θ k c a r r = Δ θ k 1 c a r r + 2 π f k 1 c a r r T c o h + π f ˙ k 1 c a r r T c o h 2 2 π u k 1 c a r r T c o h f k c a r r = f k 1 c a r r + f ˙ k 1 c a r r T c o h f ˙ k c a r r = f ˙ k 1 c a r r Δ θ k c o d e = Δ θ k 1 c o d e + β f k 1 c a r r T c o h + 1 2 β f ˙ k 1 c a r r T c o h 2 u k 1 c o d e T c o h + β ( f L 1 f I F ) T c o h
The relation between the observation Δ θ k , k 1 c a r r and state variables can still be represented by Equation (4), while Δ θ k , k 1 c o d e is written as:
Δ θ k , k 1 c o d e = Δ θ k c o d e 1 2 β f k c a r r T c o h + 1 6 β f ˙ k c a r r T c o h 2 + 1 2 u k 1 c o d e T c o h 1 2 β ( f L 1 f I F ) T c o h
Then, the state equation and observation equation of pre-filter-S can be expressed as:
{ X k c a c o = Φ k , k 1 c a c o X k 1 c a c o + M k 1 c a c o + W k 1 c a c o Z k c a c o = H k c a c o X k c a c o + N k c a c o + V k c a c o
where:
Φ k , k 1 c a c o = [ 1 2 π T c o h π T c o h 2 0 0 1 T c o h 0 0 0 1 0 0 β T c o h 1 2 β T c o h 2 1 ]
M k 1 c a c o = [ 2 π T c o h u k 1 c a r r 0 0 β ( f L 1 f I F ) T c o h T c o h u k 1 c o d e ] T
H k c a c o = [ 1 π T c o h 1 3 π T c o h 2 0 0 1 2 β T c o h 1 6 β T c o h 2 1 ]
N k c a c o = [ π T c o h u k 1 c a r r 1 2 β ( f L 1 f I F ) T c o h + 1 2 T c o h u k 1 c o d e ] T
Z k c a c o = [ Δ θ k , k 1 c a r r Δ θ k , k 1 c o d e ]
while the feedback for carrier tracking is still represented by Equation (17), the feedback for code tracking is expressed as:
u k c o d e = 1 T c o h Δ θ ^ k c o d e + β f ^ k c a r r + 1 2 β f ˙ ^ k c a r r T c o h + β ( f L 1 f I F )
Furthermore, it is noted that Δ θ k , k 1 c o d e is uncorrelated with Δ θ k , k 1 c a r r [15]. In fact, based on statistical analysis, it can be drawn that the correlation coefficient between Δ θ k , k 1 c a r r and Δ θ k , k 1 c o d e is usually smaller than 0.01. The measurement covariance matrix R k c a c o , therefore, is set as:
R k c a c o = [ σ 1 2 0 0 σ 2 2 ]
where σ 1 2 and σ 2 2 are the variances at the outputs of the carrier and code discriminators respectively. In general, σ 1 2 and σ 2 2 can be calculated as follows [10]:
σ 1 2 = 1 2 C / N 0 T c o h ( 1 + 1 2 C / N 0 T c o h )
σ 2 2 = d 4 C / N 0 T c o h [ 1 + 2 ( 2 d ) C / N 0 T c o h ]
where d is the early-late correlator spacing.
Pre-filter-S usually has better performance than the popular error-state filter (pre-filter-E). Taking PRN25 for example, Figure B1 shows the differences between tracking errors obtained by pre-filter-E and pre-filter-S. According to statistic analysis, the RMS values of Doppler shift and code phase tracking errors corresponding to pre-filter-E are respectively 1.028 Hz and 1.62   × 10 2 chip. For pre-filter-S, these two values are 0.487 Hz and 1.55   × 10 2 chip.
Figure B1. Tracking errors of pre-filter-E and pre-filter-S.
Figure B1. Tracking errors of pre-filter-E and pre-filter-S.
Sensors 16 00627 g022

References

  1. Gao, Y.B.; Liu, S.F.; Atia, M.M.; Noureldin, A. INS/GPS/LiDAR Integrated Navigation System for Urban and Indoor Environments Using Hybrid Scan Matching Algorithm. Sensors 2015, 15, 23286–23302. [Google Scholar] [CrossRef] [PubMed]
  2. Jwo, D.J.; Yang, C.F.; Chuang, C.H.; Lin, K.C. A Novel Design for the Ultra-Tightly Coupled GPS/INS Navigation System. J. Navig. 2012, 65, 717–747. [Google Scholar] [CrossRef]
  3. Tang, X.H.; Falco, G.; Falletti, E.; Presti, L.L. Theoretical analysis and tuning criteria of the Kalman filter-based tracking loop. GPS Solut. 2014, 19, 489–503. [Google Scholar] [CrossRef]
  4. Buck, T.M.; Wilmot, J.; Cook, M.J. A High G, MEMS Based, Deeply Integrated, INS/GPS, Guidance, Navigation and Control Flight Management Unit. In Proceedings of the 2006 IEEE/ION Position Location and Navigation Symposium, Coronado, CA, USA, 25–27 April 2006; pp. 772–794.
  5. Landis, D.; Thorvaldsen, T.; Fink, B.; Sherman, P.; Holmes, S. A Deep Integration Estimator for Urban Ground Navigation. In Proceedings of the 2006 IEEE/ION Position Location and Navigation Symposium, Coronado, CA, USA, 25–27 April 2006; pp. 927–932.
  6. Ohlmeyer, E.J. Analysis of an Ultra-Tightly Coupled GPS/INS System in Jamming. In Proceedings of the 2006 IEEE/ION Position Location and Navigation Symposium, Coronado, CA, USA, 25–27 April 2006; pp. 44–53.
  7. Dyrad, L.; Woessner, B.; Jovancevic, A.; Ganguly, S. Ultra Tightly Coupled GPS/INS Receiver for TSPI Applications. In Proceedings of the 20th International Technical Meeting of the Satellite Division of the Institute of Navigation, Fort Worth, TX, USA, 25–28 September 2007; pp. 2563–2572.
  8. Lashley, M.; Bevly, D.M.; Hung, J.Y. Analysis of Deeply Integrated and Tightly Coupled Architectures. In Proceedings of the 2010 IEEE/ION Position, Location and Navigation Symposium, Indian Wells, CA, USA, 3–6 May 2010; pp. 382–396.
  9. Patapoutian, A. On Phase-Locked Loops and Kalman Filters. IEEE Trans. Commun. 1999, 47, 670–672. [Google Scholar] [CrossRef]
  10. O'Driscoll, C.; Petovello, M.G.; Lachapelle, G. Choosing the Coherent Integration Time for Kalman Filter-Based Carrier-Phase Tracking of GNSS Signals. GPS Solut. 2011, 15, 345–356. [Google Scholar] [CrossRef]
  11. Petovello, M.G.; Lachapelle, G. Comparison of Vector-Based Software Receiver Implementations with Application to Ultra-Tight GPS/INS Integration. In Proceedings of the 19th International Technical Meeting of the Satellite Division of the Institute-of-Navigation, Fort Worth, TX, USA, 26–29 September 2006; pp. 1790–1799.
  12. Crane, R.N. A Simplified Method for Deep Coupling of GPS and Inertial Data. In Proceedings of the 2007 Institute of Navigation, National Technical Meeting, San Diego, CA, USA, 22–24 January 2007; pp. 311–319.
  13. Luo, Y.; Babu, R.; Wu, W.Q. Double-Filter Model with Modified Kalman Filter for Baseband Signal Pre-processing with Application to Ultra-Tight GPS/INS Integration. GPS Solut. 2012, 16, 463–476. [Google Scholar] [CrossRef]
  14. Chiang, K.W.; Duong, T.T.; Liao, J.K. The Performance Analysis of a Real-Time Integrated INS/GPS Vehicle Navigation System with Abnormal GPS Measurement Elimination. Sensors 2013, 13, 10599–10622. [Google Scholar] [CrossRef] [PubMed]
  15. O’Driscoll, C.; Lachapelle, G. Comparison of Traditional and Kalman Filter Based Tracking Architectures. In Proceedings of the 13th European Navigation Conference, Naples, Italy, 3–6 May 2009.
  16. Miao, J.F. Research on Key Technologies for Robust and Adaptive GPS Software Receiver. Ph.D. Thesis, Nanjing University of Aeronautics and Astronautics, Nanjing, China, December 2009. [Google Scholar]
  17. Won, J.H.; Pany, T.; Eissfeller, B. Characteristics of Kalman filters for GNSS signal tracking loop. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 3671–3681. [Google Scholar] [CrossRef]
  18. Xie, G. Principles of GPS and Receiver Design; Publishing House of Electronics Industry: Beijing, China; pp. 276–277.
  19. Titterton, D.; Weston, J. Strapdown Inertial Navigation Technology; The Institution of Engineering and Technology: Reston, VA, USA, 2004; pp. 377–418. [Google Scholar]
  20. Gao, G.J.; Lachapelle, G. INS-Assited High Sensitivity GPS Receivers for Degraded Signal Navigation. In Proceedings of the 19th International Technical Meeting of the Satellite Division of The Institute of Navigation, Fort Worth, TX, USA, 26–29 September 2006.
  21. Kim, K.H.; Lee, J.G.; Park, C.G. Adaptive Two-Stage Extended Kalman Filter for a Fault-Tolerant INS-GPS Loosely Coupled System. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 125–137. [Google Scholar]
  22. Hong, S.P.; Lee, M.H.; Chun, H.H. Observability of error states in GPS/INS integration. IEEE Trans. Veh. Technol. 2005, 54, 731–743. [Google Scholar] [CrossRef]
Figure 1. The carrier loop using pre-filter.
Figure 1. The carrier loop using pre-filter.
Sensors 16 00627 g001
Figure 2. The structure of the ultra-tight GPS/INS integration.
Figure 2. The structure of the ultra-tight GPS/INS integration.
Sensors 16 00627 g002
Figure 3. Satellites sky plot of the field test.
Figure 3. Satellites sky plot of the field test.
Sensors 16 00627 g003
Figure 4. Doppler shift tracking error for PRN25 in field test.
Figure 4. Doppler shift tracking error for PRN25 in field test.
Sensors 16 00627 g004
Figure 5. Code phase tracking error for PRN25 in field test.
Figure 5. Code phase tracking error for PRN25 in field test.
Sensors 16 00627 g005
Figure 6. Position estimation errors with field data.
Figure 6. Position estimation errors with field data.
Sensors 16 00627 g006
Figure 7. Velocity estimation errors with field data.
Figure 7. Velocity estimation errors with field data.
Sensors 16 00627 g007
Figure 8. Frequency difference between local and input carriers in field tests.
Figure 8. Frequency difference between local and input carriers in field tests.
Sensors 16 00627 g008
Figure 9. The reference trajectory for simulation.
Figure 9. The reference trajectory for simulation.
Sensors 16 00627 g009
Figure 10. The (strapdown) INS used in the missile test.
Figure 10. The (strapdown) INS used in the missile test.
Sensors 16 00627 g010
Figure 11. Satellites sky plot of the simulation.
Figure 11. Satellites sky plot of the simulation.
Sensors 16 00627 g011
Figure 12. The absolute acceleration of the missile.
Figure 12. The absolute acceleration of the missile.
Sensors 16 00627 g012
Figure 13. Doppler shift tracking error for PRN5 and the estimation of δ v G , 1 .
Figure 13. Doppler shift tracking error for PRN5 and the estimation of δ v G , 1 .
Sensors 16 00627 g013
Figure 14. GPS position estimation errors with simulation data.
Figure 14. GPS position estimation errors with simulation data.
Sensors 16 00627 g014
Figure 15. GPS velocity estimation errors with simulation data.
Figure 15. GPS velocity estimation errors with simulation data.
Sensors 16 00627 g015
Figure 16. INS position estimation errors.
Figure 16. INS position estimation errors.
Sensors 16 00627 g016
Figure 17. INS velocity estimation errors.
Figure 17. INS velocity estimation errors.
Sensors 16 00627 g017
Figure 18. GPS/INS position estimation errors.
Figure 18. GPS/INS position estimation errors.
Sensors 16 00627 g018
Figure 19. GPS/INS velocity estimation errors.
Figure 19. GPS/INS velocity estimation errors.
Sensors 16 00627 g019
Figure 20. GPS/INS (tight integration, traditional navigation filter) position estimation errors.
Figure 20. GPS/INS (tight integration, traditional navigation filter) position estimation errors.
Sensors 16 00627 g020
Figure 21. GPS/INS (tight integration, traditional navigation filter) velocity estimation errors.
Figure 21. GPS/INS (tight integration, traditional navigation filter) velocity estimation errors.
Sensors 16 00627 g021
Table 1. C / N 0 of selected satellites.
Table 1. C / N 0 of selected satellites.
Satellite C / N 0 (dB · Hz)Satellite C / N 0 (dB Hz)
PRN2248.881PRN1445.612
PRN3148.025PRN3245.693
PRN2545.987PRN1245.202
PRN1848.086
Table 2. RMS of position estimation errors.
Table 2. RMS of position estimation errors.
MethodRMS of Position Error (m)
X Y Z
Pre-filter-S0.3590.5850.519
Pre-filter-D0.1710.3670.315
Table 3. RMS of velocity estimation errors.
Table 3. RMS of velocity estimation errors.
MethodRMS of Velocity Error (m/s)
X Y Z
Pre-filter-S0.0440.0850.072
Pre-filter-D0.0260.0360.036
Table 4. Gyro and accelerometer specifications.
Table 4. Gyro and accelerometer specifications.
GyrosAccelerometers
In-Run Bias Error0.5 ° / h 2   × 10 4   g
Noise (1 σ , smoothing window: 10 ms)15 ° / h 2   × 10 2   g

Share and Cite

MDPI and ACS Style

Zhang, X.; Miao, L.; Shao, H. Tracking Architecture Based on Dual-Filter with State Feedback and Its Application in Ultra-Tight GPS/INS Integration. Sensors 2016, 16, 627. https://doi.org/10.3390/s16050627

AMA Style

Zhang X, Miao L, Shao H. Tracking Architecture Based on Dual-Filter with State Feedback and Its Application in Ultra-Tight GPS/INS Integration. Sensors. 2016; 16(5):627. https://doi.org/10.3390/s16050627

Chicago/Turabian Style

Zhang, Xi, Lingjuan Miao, and Haijun Shao. 2016. "Tracking Architecture Based on Dual-Filter with State Feedback and Its Application in Ultra-Tight GPS/INS Integration" Sensors 16, no. 5: 627. https://doi.org/10.3390/s16050627

APA Style

Zhang, X., Miao, L., & Shao, H. (2016). Tracking Architecture Based on Dual-Filter with State Feedback and Its Application in Ultra-Tight GPS/INS Integration. Sensors, 16(5), 627. https://doi.org/10.3390/s16050627

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop