Next Article in Journal
Asymmetric Directional Multicast for Capillary Machine-to-Machine Using mmWave Communications
Next Article in Special Issue
Optimal Design of a Center Support Quadruple Mass Gyroscope (CSQMG)
Previous Article in Journal
Photometric Calibration and Image Stitching for a Large Field of View Multi-Camera System
Previous Article in Special Issue
Mechanical Coupling Error Suppression Technology for an Improved Decoupled Dual-Mass Micro-Gyroscope
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sample-Wise Aiding in GPS/INS Ultra-Tight Integration for High-Dynamic, High-Precision Tracking

School of Electronic and Information Engineering, Beihang University (BUAA), Beijing 100191, China
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(4), 519; https://doi.org/10.3390/s16040519
Submission received: 7 March 2016 / Revised: 3 April 2016 / Accepted: 6 April 2016 / Published: 11 April 2016
(This article belongs to the Special Issue Inertial Sensors and Systems 2016)

Abstract

:
By aiding GPS receiver tracking loops with INS estimates of signal dynamics, GPS/INS ultra-tight coupling can improve the navigation performance in challenging environments. Traditionally the INS data are injected into the loops once every loop update interval, which limits the levels of dynamics accommodated. This paper presents a sample-wise aiding method, which interpolates the aiding Doppler into each digital sample of the local signal to further eliminate the dynamic errors. The relationship between the tracking error and the aiding rate is derived analytically. Moreover, the effects of sample-wise aiding using linear and spline interpolations are simulated and compared with traditional aiding under different INS data update rates. Finally, extensive tests based on a digital IF (intermediate frequency) signal simulator and a software receiver validate the theoretical equations and demonstrate that the dynamic stress error can be significantly reduced by sample-wise aiding.

Graphical Abstract

1. Introduction

GPS (global positioning system) and INS (inertial navigation system) have complementary characteristics and their integration has been widely explored [1]. High dynamics and high precision impose conflicting requirements on noise bandwidths of conventional tracking loops in GPS receivers [2]. GPS/INS ultra-tight coupling or deep integration can relieve this conflict by aiding the tracking loops with inertial measurements and, thus, boost immunity to high dynamics, severe interferences and brief line-of-sight blockages [3,4,5]. The idea of using ultra-tight GPS/INS integration to improve performance in challenging environments has drawn wide attention of the navigation community [6,7,8,9,10,11,12,13]. Yang, Y. & El-Sheimy, N. discussed the implementation and performance analysis of a PLL (phase locked loop) aided by a MEMS (micro electromechanical system) IMU (inertial measurement unit) [10]. Details of the IMU mechanization and state propagation can be found in [8]. Alban, S., et al., used a feed-forward aided PLL model and its transfer function to facilitate the performance analysis of carrier tracking loops with external aiding [11]. Gebre-Egziabher, D., et al., quantified the effects of tracking loop bandwidth, receiver clock jitter, aiding Doppler error and platform vibration on the quality of the carrier phase estimate provided by the aided PLL [12]. Chiou, T.Y. took account of different orders of carrier tracking loops with external aiding and detailed the process of data collection and experimental setup [13]. With the development of vector tracking based receiver techniques [14], the GPS/INS deep coupling architecture (where the tracking loops are implemented as part of a larger navigation filter) has been studied increasingly by a number of authors. Gustafson, D., et al., proposed a vector-based GPS receiver which can be deeply integrated with inertial and alternate sensors [15]. Petovello, M.G. & Lachapelle, G. investigated three different Kalman filter implementation options in both scalar-tracking mode and vector-tracking mode with particular consideration for carrier phase tracking [16]. The analyses and tests in [17,18] showed improved tracking and reacquisition performances of vector-based deep integration at the cost of increased computational loads. Babu, R. & Wang, J. noticed that, typically, the integration filter runs at a rate of 1 to 100 Hz, while the carrier tracking loop normally runs at about 1000 Hz [19]. An interpolator based on a CIC (cascaded integrator comb) filter is proposed to interpolate the low rate Doppler measurements from the navigation filter to 1000 Hz. Babu, R. & Wang, J. concerned the quality of IMU-derived Doppler estimates for integration, and derived the relationship between the baseband I (in-phase)/Q (quadrature) and the INS error for performance analysis [20]. Several studies addressed the algorithms of the pre-filter and the centralized navigation filter [21,22]. Hwang, D.B., et al., provided a unified frame applicable to several different algorithms of GPS/INS deep fusion [23]. Xie, F., et al., proposed an adaptive robust deeply coupled GNSS/INS system based on a federated architecture which can process the GPS L1 and BeiDou B1 signals together [24]. Qin, F., et al., assessed the performance using low-cost IMU and vector tracking-based integration for high dynamic applications [25]. Recent work in [26] presented a federated GPS/INS deep coupling system employing three Kalman filters and providing adaptation ability to improve the navigation performance under various jamming conditions.
All of these reported aiding methods, however, inject the externally-estimated carrier Doppler directly into the tracking loop as a constant during the coherent integration period, which is at least 1 millisecond (ms). This limits the aiding rate not higher than 1000 Hz which is insufficient to the suppression of high dynamic stress errors. For example, the carrier Doppler may change up to 10.5 Hz over a 1 ms interval under a 2000 m/s2 acceleration. Therefore, the effects of the traditional aiding are restricted.
This paper proposes a sample-wise aiding method, which interpolates the aiding data onto each sampling point before controlling local reference signal generation to further eliminate the high-order dynamic errors with only minor modification of the receiver NCO (numerically controlled oscillator). Section 2 goes deep into the digital realization of the Doppler-aided PLL to explore the influence of the aiding data update rate on the scalar tracking performance. Theoretical derivation of the phase error transfer function reveals that the aiding effect has a lot to do with the product of the frequency of signal dynamics and the aiding update interval. Section 3 presents the sample-wise aiding method to break the performance limitation of the traditionally aided PLL, and simulates the effects of linear and spline interpolations in comparison with the traditional aiding. In Section 4 a testbed for the performance limitation of GPS/INS ultra-tight integration is implemented, which employs a digital IF (intermediate frequency) GPS signal simulator to test the tracking errors of a software GPS receiver ultra-tightly coupled with INS data in various high dynamic scenarios. The test results coincide with theoretical expectations and demonstrate that the sample-wise aiding can eliminate the dynamic errors much more effectively than the traditional aiding.

2. Performance of Traditional Doppler Aiding

Both the scalar tracking-based ultra-tight integration and the vector tracking-based deep integration use Doppler-aided tracking loops, which allows the INS or the integrated navigation filter to account for the platform dynamics such that the receiver can narrow its loop bandwidths. To simplify the analysis of the effects of the aiding data update rate, this section will focus on the dynamic process of the scalar tracking based PLL.
According to the fundamental theory of PLL [27], the linear phase model of a PLL can be extended to describe the dynamic process of a PLL aided by external Doppler data, as shown in Figure 1.
In Figure 1, θ e ( t ) = θ i ( t ) θ o ( t ) is the phase error between the input phase θ i ( t ) and the PLL output phase θ o ( t ) ; f a ( t ) is the aiding Doppler provided by the INS estimate; F ( p ) denotes the transfer operator of the loop filter in the time-domain analysis, whose counterpart is F ( j ω ) in the frequency domain, and F ( s ) in the Laplace domain.
For an unaided PLL, the closed-loop phase transfer function is H ( s ) = F ( s ) s + F ( s ) , and the phase error transfer function is H e ( s ) = 1 H ( s ) = s s + F ( s ) . The phase error θ e ( s ) is determined by the input signal phase θ i ( s ) and the phase noise n ( s ) as follows:
θ e ( s ) = H e ( s ) θ i ( s ) H ( s ) n ( s )
Generally the term “high dynamics” is embodied in two aspects: high frequency and large amplitude of the range/velocity/acceleration/jerk variations. For the case study of the behavior of the loop under high dynamics, a typical sinusoidal trajectory in the LOS (line of sight) direction is selected in this paper due to its following characteristics: (1) any complicated LOS dynamics can be decomposed into a number of sinusoidal components with different amplitudes and frequencies; and (2) it can provide severe high-order dynamics with continuous derivatives and limited range of motion, which is easy to control and convenient for testing. Therefore the LOS Doppler is simply modeled as:
f d ( t ) = A sin ( ω t )
The carrier phase of the incoming signal is:
θ i ( t ) = A ω cos ( ω t ) + θ 0
Since second-order and higher-order PLLs are not sensitive to frequency/phase steps, the steady state carrier phase tracking error is only caused by the first item in Equation (3):
θ e ( t ) = A ω | H e ( ω ) | cos [ ω t + Φ e ( ω ) ]
It is seen that the sinusoidal dynamic stress error cannot be completely eliminated, and the amplitude of the error transfer function | H e ( ω ) | represents the ratio of the amplitude of the phase error to that of the incoming phase.
For a Doppler aided PLL, it can be easily derived from Figure 1 such that:
θ e ( s ) = H e ( s ) θ i a ( s ) H ( s ) n ( s )
θ i a ( s ) = θ i ( s ) θ a ( s )
θ a ( s ) = 1 s f a ( s )
where θ i a ( s ) is defined as the equivalent input phase of the aided PLL, and θ a ( s ) is the feedback phase introduced by Doppler aiding.
Ideally, if the INS-estimated Doppler is exactly the derivative of the input phase, there will be θ i a ( t ) = 0 and θ e ( s ) = H ( s ) n ( s ) , which means that the LOS dynamics can be completely eliminated by Doppler aiding. As a matter of fact, the signal frequency shift is also contributed by the receiver clock drift, the satellite clock drift, and even induced by the ionosphere/troposphere propagation, (which are beside the point of the paper). In this case, the design of the loop filter only needs to account for the phase noise instead of tracking the signal dynamics. The analog domain analysis, however, cannot cover all aspects of a digital realization. (1) In reality, a PLL in a modern receiver is implemented in its digital version—the DPLL where, traditionally, the aiding Doppler value remains constant over an external data update interval Ta (which may be equal to, smaller than, or larger than the pre-detection integration time T; nevertheless most reported systems keep T a T as considered in [19]); whereas the Doppler variation within this interval cannot be ignored in high dynamic situations; and (2) the carrier NCO control word actually comes from the Doppler estimate of the previous interval rather than the current interval. Taking these two factors into account, we replace the analog NCO integrator in Figure 1 with a digital accumulator operated at the sampling rate R s , as shown in Figure 2. Note that another input of the adder before the integrator in Figure 1, namely the loop filter output, has been omitted in Figure 2, being dedicated to investigating the feedback phase of the aiding data.
Assuming that the aiding Doppler is provided at an update rate R a = 1 / T a , and the Doppler data update interval is the integral multiple of the digital signal sample interval T s : m = T a T s = R s R a .
At the nth sample time n T s = ( k m + l ) T s = k T a + l T s , l = 0 , 1 , , m 1 , the feedback phase can be expressed as θ a ( n T s ) = p = 0 k 1 f a ( p T a ) T a + f a ( k T a ) l T s . The equivalent input phase of the loop at each aiding Doppler update epoch (n = km, l = 0) is:
θ i a ( k T a ) = θ i ( k T a ) θ a ( k T a ) = θ i ( k T a ) p = 0 k 1 f a ( p T a ) T a = p = 0 k 1 [ θ i ( p T a + T a ) θ i ( p T a ) f a ( p T a ) T a ] + θ i ( 0 )
Assuming that the INS can provide exact Doppler estimates synchronized with the LOS signal Doppler, we have f a ( p T a ) = A sin ( ω p T a ) for the sinusoidal input carrier phase in Equation (3). Substituting them into Equation (8) and neglecting the constant item (because it contributes nothing to the steady state tracking error), we obtain:
θ i a ( k T a ) = p = 0 k 1 [ A ω cos ( ω p T a + ω T a ) + A ω cos ( ω p T a ) A sin ( ω p T a ) T a ] = A p = 0 k 1 [ C a 2 ( ω T a ) sin ( ω p T a ) + C a 1 ( ω T a ) cos ( ω p T a ) ] T a = A C a ( ω T a ) p = 0 k 1 sin [ ω p T a + ϕ a ( ω T a ) ] T a
C a 1 ( ω T a ) = 1 cos ( ω T a ) ω T a = C a ( ω T a ) sin [ ϕ a ( ω T a ) ] ,    C a 2 ( ω T a ) = sin ( ω T a ) ω T a ω T a = C a ( ω T a ) cos [ ϕ a ( ω T a ) ] C a ( ω T a ) = C a 1 2 ( ω T a ) + C a 2 2 ( ω T a ) ,    ϕ a ( ω T a ) = atan 2 ( 1 cos ω T a , sin ω T a ω T a )
When the product of the frequency of signal dynamics and the aiding update interval satisfies ω T a 2 π , the accumulation of discrete sinusoidal signals in Equation (9) can be approximated by the integration of the corresponding analog signal as the following:
θ i a ( k T a ) A C a ( ω T a ) 0 k T a sin [ ω t + ϕ a ( ω T a ) ] d t = A ω C a ( ω T a ) cos [ ω k T a + ϕ a ( ω T a ) ] + A ω C a ( ω T a ) cos [ ϕ a ( ω T a ) ]
Consequently, the phase error of the aided loop will follow the sinusoidal shape with the same angular frequency of the input signal phase as shown in Equation (4), yet with an amplitude of C a ( ω T a ) times that of the unaided PLL. The equivalent phase error transfer function of the aided loops is:
H e a ( ω ) = H e ( ω ) C a ( ω T a ) exp [ ϕ a ( ω T a ) ]
When ω T a > 2 π , we can find a frequency ω b meeting ω b T a < 2 π and ω T a = ω b T a + m 2 π , which yields θ i a ( k T a ) A ω C a ( ω T a ) cos [ ω b k T a + ϕ a ( ω T a ) ] + A ω C a ( ω T a ) cos [ ϕ a ( ω T a ) ] . For this case Equation (12) still suffices. It is accurate enough, provided that the product of ω b and T a satisfies ω b T a 2 π , as can be validated by the simulation in Section 3.2.
Thus, C a ( ω T a ) is the ratio of the phase error amplitudes of aided loop to unaided loop. Figure 3 shows how this ratio varies with ω T a according to Equation (10). Note that the C a ( ω T a ) expression is applicable to any loop orders except for 1.
It can be seen that: (1) If ω T a < 2.3311 , the Doppler aiding can improve tracking performance by C a ( ω T a ) < 1 ; (2) When ω T a 2 π , C a ( ω T a ) 1 , the aiding can hardly help to mitigate the dynamic stress error; (3) When ω T a 2 π , C a ( ω T a ) 1 2 ω T a , lim ω T a 0 C a ( ω T a ) = 0 . Hence, for the traditional Doppler aided loop, the aiding rate is expected to be as high as possible to efficiently reduce the dynamic stress error. However, the practical INS-offered data rate may not satisfy the cases with high frequency component or large amplitude of the LOS dynamics.

3. Sample-Wise Doppler-Aided PLL

3.1. Ideal Sample-Wise Doppler Aiding

In order to reduce the high dynamic stress error, it is effective to increase the aiding rate. Actually, the external aiding data can be interpolated to each sample point of the digital signal, and the control words of the carrier NCO for generating the local reference carrier can be updated at a rate up to the sampling rate R s without major modification of the receiver architecture, as shown in Figure 4. This is the main idea of the sample-wise Doppler aiding method proposed in this section. Compared with traditionally-aided DPLL in Figure 2, the only additional cost of sample-wisely aided DPLL is the frequency interpolation unit before the NCO phase accumulator.
In Figure 4 the control word of the carrier NCO is f s a ( 0 : T s : T p ) = interp1 [ f a ( 0 : T a : T p ) ] , where interp1[·] is the one-dimensional data interpolation function and T p is the time length of the data. In contrast to the traditionally aided DPLL, the equivalent input carrier phase in a sample-wisely Doppler aided DPLL becomes:
θ i a ( n T s ) = θ i ( n T s ) θ a ( n T s ) = θ i ( n T s ) q = 0 n 1 f s a ( q T s ) T s
If there are still no aiding Doppler error and no time alignment error, ignoring the interpolation error (which will be investigated in the next subsection), we can substitute f s a ( q T s ) = A sin ( q T s ) and θ i ( t ) = A ω cos ( ω t ) into Equation (13) to obtain:
θ i a ( n T s ) A ω C a ( ω T s ) cos [ ω n T s + ϕ a ( ω T s ) ] + A ω C a ( ω T s ) cos [ ϕ a ( ω T s ) ]
It is easy to meet that:
C a ( ω T s ) 1 2 ω T s
Thereby the equivalent phase error transfer functions of the sample-wisely aided loop is:
H e s a ( ω ) = H e ( ω ) C a ( ω T s ) exp [ ϕ a ( ω T s ) ]
The equations indicate that the aiding performance can be significantly improved by increasing the aiding rate from R a to R s . Smaller dynamic stress error can be achieved with a higher sampling frequency.
As stated in [14], a second-order PLL is superior to a third-order PLL when aided. Taking typical third-order and second-order PLL designs in [2] as examples, Figure 5 plots the amplitudes of the phase error transfer functions of the unaided second-order PLL, the unaided third-order PLL, the traditionally-aided second-order PLL, and the sample-wisely-aided second-order PLL, respectively, with the PLL bandwidth BPLL = 15 Hz, the Doppler aiding rate R a = 1 / T a = 1 KHz , and the signal sampling rate R s = 1 / T s = 10 MHz .
Figure 5 shows that:
  • When the frequency of the LOS dynamic vibration goes much higher than the natural frequency of the loop, the unaided carrier tracking loop has an amplitude magnification of near 1 with hardly any suppression on the dynamic stress error.
  • Except for the frequency range of 6.7 < ω < 28.6 (rad/s), the unaided third-order PLL shows better dynamic performance than the unaided second-order PLL.
  • For the dynamic vibrations with ω > 2331 rad / s , the traditional Doppler aiding cannot improve tracking performance efficiently.
  • Since T s is normally several orders of magnitude smaller than T a , the phase error can be significantly reduced by sample-wise Doppler aiding.

3.2. Effects of Interpolation Techniques

In practice, we cannot obtain the exact value of the signal Doppler at each sample point by interpolation. This section explores the effects of different interpolation techniques on the dynamic stress error, still ignoring the errors in the INS estimates. Two typical interpolation approaches are investigated: (1) the linear interpolation; and (2) the cubic spline interpolation.
A linear interpolation is generally the simplest that yields a good result when up-sampling. Figure 6 illustrates the block diagram of a second-order PLL sample-wisely-aided by the Doppler data linearly-interpolated from the INS data rate R a to the sampling rate R s .
Figure 6 follows the phase model of Figure 1 in the time domain and omits the input noise for simplicity. The components shown in purple are operated at the sampling rate R s and often implemented in hardware for real-time processing; whereas the components shown in white and green are operated at the loop filter update rate R (usually the reciprocal of the pre-detection integration time T ) and the INS data update rate R a respectively. One exception is the carrier phase discriminator, which is usually composed of the carrier lookup table, the multipliers for carrier and code wipe-offs, and the integration and dump accumulator. The accumulator performs low-pass filtering and down-samples the data from R s to R = 1 / T . The cost of replacing traditional aiding with sample-wise aiding is the linear interpolation unit shown in the top block of Figure 6, which can also be regarded as an extra stage of NCO accumulator.
In high dynamic situations, the simple linear interpolation may not meet the requirements of curve smoothness and approximation. More complicated up-sampling should be used to construct a smoother velocity and approximate complex trajectory. Among the interpolation techniques including polynomials, splines, Gaussian/low-pass processes, wavelets, etc., the cubic spline is very popular for its tradeoff between accuracy and computational cost. The trajectory constructed by cubic splines makes the velocity, acceleration, and jerk (VAJ) continuous, and the VAJ of the initial and final positions can be flexibly configured [28,29].
Simulations are carried out for the two interpolation techniques from a wide range of INS data rates to three sampling rates. Figure 7 compared their effects on the ratios of phase error amplitudes of aided PLLs to unaided loops, with (a) ω = 1 rad / s , (b) ω = 50 π rad / s , respectively. The simulation results of traditional aiding are also included and compared with the theoretical prediction in Equation (10).
It can be seen that:
  • Equation (10) is accurate to quantify the effectiveness of traditional Doppler aiding, regardless of the aiding data error and the time alignment error.
  • When ω T a > 2.33 rad / s , any aiding technique cannot suppress the dynamic stress error efficiently and, thus, is not worth the bother.
  • When Ta begins to decrease from its upper threshold T u ( ω ) 2.33 / ω , sample-wise aiding with linear interpolation shows better tracking performance than traditional aiding, yet its improvement cannot catch up with the spline interpolation, until that T a reaches a certain lower threshold, denoted by T l ( ω , R s ) . When T a < T l ( ω , R s ) , the ratio of the phase error amplitudes of sample-wisely aided PLL to unaided PLL will be only determined by the theoretical limit C a ( ω T s ) 0.5 ω T s , independent of T a , and cannot be improved by employing other more advanced interpolation techniques. In addition, before reaching the lower threshold, there is a middle threshold T m ( ω , R s ) of Ta, such that for the T m ( ω , R s ) < T a < T u ( ω ) region, the improvement effect of sample-wisely aided PLL to unaided PLL is only determined by T a and the interpolation method, independent of R s , and can be improved by employing other more advanced interpolation techniques.
It is interesting to divide each sample-wise aiding curve into the following four regions by using three thresholds of Ta as shown in Table 1.
By contrast, there are only two regions for the traditional aiding: (1) T a < T u ( ω ) : C a ( ω T a ) 0.5 ω T a ; (2) T a T u ( ω ) : C a ( ω T a ) 1 .
It is also seen that the spline interpolation extends both the lower threshold T l ( ω , R s ) and the middle threshold T m ( ω , R s ) towards longer T a than the linear interpolation. Therefore, the pressure of high INS data updating rate can be relieved by employing advanced interpolation techniques to achieve the same performance as the linear interpolation. Additionally, higher dynamic frequency ω demands shorter T a thresholds. Yet the quantitative relationships of the T a thresholds T l ( ω , R s ) and T m ( ω , R s ) with the signal dynamic frequency ω and the sampling rate R s with different interpolation techniques are quite complicated and still under investigation.

4. Test Results and Discussion

4.1. Test Setups and Scenarios

The scalar tracking based GPS/INS ultra-tight integration has been implemented in our software receiver. Both the traditional Doppler aiding and the sample-wise aiding with spline interpolation are incorporated into the receiver PLL and DLL (delay locked loop), and their tracking performance in high dynamic scenarios can be tested by using our digital IF GPS signal simulator. The accurate INS Doppler estimates needed in the tests can simply be obtained from the user trajectory provided by the simulator and the satellite ephemerides.
Three factors have been considered in the design of dynamic scenarios for the tests: (1) to ensure that the receiver achieves steady state tracking before the aiding begins, no VAJ jumps are included in the trajectory, because carrier tracking loops are sensitive to dynamic jumps. For example, third-order loops are sensitive to jerk jumps; (2) the trajectory should be controlled within a limited spatial range; and (3) severe high-order dynamics should be included.
As a result the following sinusoidal user trajectory in the height direction is constructed:
h = h 0 + D [ 1 cos ( ω t ) ]
where, h is the height of the receiver antenna phase center, h0 is the initial height, D is the amplitude of the sinusoidal trajectory, and ω is its angular frequency. In the following tests, we chose ω = 1 rad / s to simplify the control and calculation of dynamic parameters of different orders. For example, if D = 500 m, then the maximum velocity, acceleration, and jerk are 500 m/s, 500 m/s2, and 500 m/s3, respectively. With a satellite elevation angle of α = 28.67 ° , the LOS velocity is:
V L O S = ( D sin α ) ω sin ( ω t ) 239.88 sin ( t )
The receiver first tracks the high dynamic signals using a third-order PLL, and then changes to a second-order PLL aided by a priori known carrier Doppler after decoding the navigation message. Meanwhile the reception time (in GPS time) of each individual sample point can be obtained by using the single-SV (satellite vehicle) timing method proposed in [30]. For high dynamic scenarios, the known user trajectory rather than a static user location is utilized in the timing process. In this way the time alignment necessary for the aiding can be achieved.
The loop parameters recommended in Tables 5.6 in [2] are adopted in the tests. In addition, the PLL is configured with a loop bandwidth of 15 Hz and a loop update interval of 1ms, while the DLL taking a 1 Hz bandwidth and a 20 ms loop update interval. For traditional aiding, the external Doppler data are injected to the PLL every 1ms, and to the DLL every 20 ms whereas, for the sample-wise aiding, the receiver first interpolates the aiding data using cubic splines from 1000 Hz to the sampling rate, and then injects the interpolated carrier Doppler and code Doppler values into the carrier NCO and the code NCO, respectively.
Finally the code phase errors and carrier phase errors can be obtained by subtracting the simulator recorded data from the receiver measurements.

4.2. Aiding Effects on Carrier Phase Errors with no Thermal Noise in the Signal

To exclude the effects of thermal noises and cross-correlation interferences, a single-SV scenario is necessary to achieve high C/N0 (carrier-to-noise ratio) and, thus, high precision dedicated for the dynamic stress error test [30]. Although no noise is added to the output signal of the digital IF signal simulator, the C/N0 estimated by the receiver is around 128 dB-Hz owing to both the inevitable quantization noise and the limited range of the wideband-narrowband power ratio method (PRM) [31]. This is near the best PRM-estimated C/N0 achieved by our test platform so far. The carrier phase tracking errors in a 5 g (≈50 m/s2) dynamic condition are shown in Figure 8.
In the first 30 s, the receiver is in unaided state as the epoch time has not been solved. The top plot in Figure 8 shows that the carrier phase error is greatly reduced after external Doppler aiding being turned on. The bottom plot discards the first 30s data and examines the aided period. It can be seen that the carrier phase error remains a sinusoidal wave with the same angular frequency of 1 rad/s, as predicted by Section 2.
The tests have been conducted using six sets of signals with the dynamics of 5 g, 20 g, 50 g, 100 g, 200 g, and 500 g, respectively. The carrier phase tracking errors using traditional Doppler aiding are shown in Figure 9a. Figure 9b calculates their amplitudes by combining Equations (4), (10) and (12), and the experimental results shown in the red triangle symbols are well in accordance with the theoretical expectation in the blue line.
Table 2 summarizes the test results of the error amplitudes for the unaided third-order PLL and the aided second-order PLL and their ratios. The theoretical calculation of this ratio follows:
r = | H e 3 ( j ω ) | C a ( ω T a ) | H e 2 ( j ω ) | 230.04
The test results demonstrate the correctness of the equations of carrier tracking errors in Section 2. The small discrepancies between the test results and the theoretical value of the ratio are due to the impacts of the time alignment errors and the quantization errors of the test platform, as well as the round-off errors of the error amplitudes.
Corresponding to Figure 9a, Figure 10 shows the tracking results using the sample-wise aided PLL with 5 g, 50 g, and 500 g dynamics. It is seen that the dynamic errors have been further eliminated. It should be noted that the tracking error achieved by the test platform even in a static user scenario is at the 0.8 × 10−6 m level in terms of the standard deviation of the carrier phase error. The errors in Figure 10 are also contributed by the alignment errors between the simulator and the receiver (related to the signal dynamics) and the quantization errors.

4.3. Aiding Effects on Code/Carrier Phase Errors with a C/N0 of 80 dB-Hz

The carrier tracking loop is usually the weak link in a receiver and the carrier phase error characterizes the receiver dynamic performance. For the traditionally-aided receiver tracking in high dynamic conditions, however, the code tracking loop may also become the dynamically weak part since the DLL has a narrower bandwidth and a longer integration time and sometimes a narrow correlator spacing for the sake of noise and interference suppression. Both the carrier phase error and the code phase error are concerned in this subsection. The signal C/N0 is configured to be 80 dB-Hz in the test for the following two reasons: (1) the 1σ thermal noise code tracking jitter is about 0.0059 m under 80 dB-Hz and a narrow correlator spacing of 0.1chip. If the C/N0 is set too low, all the test results of the aided receiver will be buried in the thermal noise jitter; and (2) one job in the future is to experiment based on a RF (radio frequency) signal simulator and a RF signal collection device, where the 80 dB-Hz C/N0 has already been achieved at RF level using a single-SV scenario [31].
Figure 11, Figure 12 and Figure 13 show the test results of the traditional aiding method and the sample-wise method with user dynamics of 5 g, 50 g and 500 g, respectively. Figure 14 summarizes the standard deviations of carrier tracking errors and code tracking errors of the two aiding methods in the test results. For comparison, Figure 14 also shows the theoretical values of the thermal noise floors at 80 dB-Hz and 45 dB-Hz (in normal signal conditions) as well as the theoretical values of the total 1σ errors at 80 dB-Hz (including contributions from both the dynamic stress and the thermal noise) for the traditional aiding method.
It can be seen that: (1) under a high C/N0 of 80 dB-Hz the dynamic stress error remains the main error source in all the traditional aiding results, and the total carrier phase errors becomes much closer to the test results in last subsection with the dynamics increasing; whereas all the sample-wise aiding results are dominated by the thermal noise jitters; and (2) although the dynamic stress errors can be significantly reduced by traditional Doppler aiding, they can still become dominant in carrier phase errors when the dynamics go higher than 3400 m/s2, and dominant in code phase errors when the dynamics go higher than 390 m/s2 under a normal C/N0 of 45 dB-Hz. As a contrast, the tracking errors remain almost the same values near the noise floor with the increase of user dynamics when using sample-wise aiding, which indicates that the dynamic stress errors have been almost completely eliminated by the sample-wise Doppler-aided receiver.

5. Conclusions

  • By examining the digital implementation of the Doppler-aided PLL in scalar tracking-based GPS/INS ultra-tight integration, a theoretical expression of the ratio of phase error amplitudes of an aided loop to an unaided loop is derived, which indicates how the aiding effects are impacted by the product of the frequency of signal dynamics and the INS aiding rate.
  • A new sample-wise aiding method in GPS/INS ultra-tight integration for high-dynamic, high-precision tracking has been presented. Its effectiveness and advantage over traditional aiding have not only been analytically expressed and numerically simulated, but also physically tested using a digital IF signal simulator and a software receiver.
  • The simulation with respect to the effects of two interpolation techniques on the sample-wise aiding performance helps to find the boundary conditions of the aiding methods in terms of the INS aiding data update rate.
Although the research in this paper is based on the scalar tracking architecture in GPS/INS integration, the sample-wise aiding method can be introduced into the vector tracking-based deep coupling to satisfy high-dynamic, high-precision requirements. In the vector tracking case, the navigation filter output can be interpolated to the sampling rate before controlling the local signal generation of each satellite channel. Similarly, the method can also be extended to the fields of map aiding and integration of GNSS with other wireless positioning. Related future work also includes the analysis of the performance limitation due to the INS data error and the misalignment between GPS and INS.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No. 61271197).

Author Contributions

Yanhong Kou conceived the sample-wise aiding idea. Yanhong Kou and Han Zhang elaborated the problem formulation, designed the simulations and the experiments. Han Zhang performed the tests. Yanhong Kou analyzed the data and wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Farrell, J.A.; Barth, M. The Global Positioning System & Inertial Navigation; McGraw-Hill: New York, NY, USA, 1999. [Google Scholar]
  2. Ward, P.W.; Betz, J.W.; Hegarty, C.J. Satellite Signal Acquisition, Tracking, and Data Demodulation. In Understanding GPS Principles and Applications, 2nd ed.; Kaplan, E.D., Hegarty, C.J., Eds.; Artech House: Norwood, MA, USA, 2006; pp. 153–241. [Google Scholar]
  3. Jwo, D. GPS receiver performance enhancement via inertial velocity aiding. J. Navig. 2001, 54, 105–117. [Google Scholar] [CrossRef]
  4. Abbott, S.; Lillo, W.E. Global Positioning Systems and Inertial Measuring Unit Ultratight Coupling Method. US Patent No. 6516021, 4 February 2003. [Google Scholar]
  5. 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, San Diego, CA, USA, 25–27 April 2006; pp. 44–53.
  6. Alban, S. Design and Performance of a Robust GPS/INS Attitude System for Automobile Applications. Ph.D. Thesis, Stanford University, Stanford, CA, USA, 2004. [Google Scholar]
  7. Lashley, M.; Bevly, D.M. A Comparison of the Performance of a Non-Coherent Deeply Integrated Navigation Algorithm and a Tightly Coupled Navigation Algorithm. In Proceedings of the 21st International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS 2008), Savannah, GA, USA, 16–19 September 2008; pp. 2123–2129.
  8. Groves, P.D. Principles of GNSS, Inertial, and Multisensor Integrated Navigation Systems; Artech House: Norwood, MA, USA, 2008. [Google Scholar]
  9. Qin, F.; Zhan, X.; Du, G. Performance Improvement of Receivers Based on Ultra-Tight Integration in GNSS-Challenged Environments. Sensors 2013, 13, 16406–16423. [Google Scholar] [CrossRef]
  10. Yang, Y.; El-Sheimy, N. Improving GPS Receiver Tracking Performance of PLL by MEMS IMU Aiding. In Proceedings of the 19th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS 2006), Fort Worth, TX, USA, 26–29 September 2006; pp. 2192–2201.
  11. Alban, S.; Akos, D.M.; Rock, S.M.; Egziabher, D.G. Performance Analysis and Architectures for INS-Aided GPS Tracking Loops. In Proceedings of the Institute of Navigation’s National Technical Meeting, Anaheim, CA, USA, 22–24 January 2003.
  12. Gebre-Egziabher, D.; Razavi, A.; Enge, P.; Gautier, J.; Pullen, S.; Pervan, B.S.; Akos, D. Sensitivity and Performance Analysis of Doppler-Aided GPS Carrier-Tracking Loops. Navigation 2005, 52, 49–60. [Google Scholar] [CrossRef]
  13. Chiou, T.Y. GPS Receiver Performance Using Inertial-Aided Carrier Tracking Loop. In Proceedings of the ION GNSS 18th International Technical Meeting of the Satellite Division, Long Beach, CA, USA, 13–16 September 2005.
  14. Spilker, J.J., Jr. Fundamentals of Signal Tracking Theory. In Global Positioning System: Theory and Applications; Parkinson, B.W., Spilker, J.J., Jr., Eds.; American Institute of Aeronautics and Astronautics Inc.: Reston, VA, USA, 1994; Volume I, pp. 290–305. [Google Scholar]
  15. Gustafson, D.; Dowdle, J.R.; Elwell, J.M. Deeply-Integrated Adaptive GPS-Based Navigator with Extended-Range Code Tracking. U.S. Patent US6331835B1, 18 December 2001. [Google Scholar]
  16. Petovello, M.G.; Lachapelle, G. Comparison of vector-based software receiver implementations with application to ultra-tight GPS/INS integration. In Proceedings of the 18th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS '06), Fort Worth, TX, USA, 26–29 September 2006; pp. 1790–1799.
  17. Lashley, M. Modeling and Performance Analysis of GPS Vector Tracking Algorithms. Ph.D. Thesis, Auburn University, Auburn, AL, USA, 18 December 2009. [Google Scholar]
  18. Lashely, M.; Bevly, D.M. Performance Comparison of Deep Integration and Tight Coupling. Navigation 2013, 60, 159–178. [Google Scholar] [CrossRef]
  19. Babu, R.; Wang, J. Comparative study of interpolation techniques for ultra-tight integration of GPS/INS/PL sensors. J. Glob. Posit. Syst. 2005, 1, 192–200. [Google Scholar] [CrossRef]
  20. Babu, R.; Wang, J. Ultra-Tight GPS/INS/PL Integration: A system Concept and Performance Analysis. GPS Solut. 2009, 13, 75–82. [Google Scholar] [CrossRef]
  21. Sivananthan, S.; Weitzen, J. Sub-Optimality of the Prefilter in Deeply Coupled Integration of GPS and INS. Navigation 2010, 57, 175–184. [Google Scholar] [CrossRef]
  22. Won, J.H.; Dotterbock, D.; Eissfeller, B. Performance Comparison of Different Forms of Kalman Filter Approaches for Vector-Based GNSS Signal Tracking Loop. Navigation 2010, 57, 185–199. [Google Scholar] [CrossRef]
  23. Hwang, D.B.; Lim, D.W.; Cho, S.L.; Lee, S.J. Unified approach to ultra-tightly-coupled GPS/INS integrated navigation system. IEEE Aerosp. Electron. Syst. Mag. 2011, 26, 30–38. [Google Scholar] [CrossRef]
  24. Xie, F.; Liu, J.; Li, R.; Hang, Y. Adaptive robust ultra-tightly coupled global navigation satellite system/inertial navigation system based on global positioning system/BeiDou vector tracking loops. IET Radar Sonar Navig. 2014, 8, 815–827. [Google Scholar] [CrossRef]
  25. Qin, F.; Zhan, X.; Zhan, L. Performance assessment of a low-cost inertial measurement unit based ultra-tight global navigation satellite system/inertial navigation system integration for high dynamic applications. IET Radar Sonar Navig. 2014, 8, 828–836. [Google Scholar] [CrossRef]
  26. Li, K.; Zhao, J.; Wang, X.; Wang, L. Federated ultra-tightly coupled GPS/INS integrated navigation system based on vector tracking for severe jamming environment. IET Radar Sonar Navig. 2016. [Google Scholar] [CrossRef]
  27. Gardner, F.M. Phaselock Techniques, 3rd ed.; John Willy & Sons: New York, NY, USA, 2005. [Google Scholar]
  28. Stoer, J.; Bulirsch, R. An Introduction of Numerical Analysis, 3rd ed.; Springer: New York, NY, USA, 2002. [Google Scholar]
  29. Bartels, R.H.; Beatty, J.C.; Barsky, B.A. An Introduction to Splines for Use in Computer Graphics and Geometric Modeling; Morgan Kaufmann Publishers, Inc.: San Francisco, CA, USA, 1987. [Google Scholar]
  30. Kou, Y.; Sui, J.; Chen, Y.; Zhang, Z.; Zhang, H. Test of Pseudorange Accuracy in GNSS RF Constellation Simulator. In Proceedings of the 25th International Technical Meeting of The Satellite Division of the Institute of Navigation (ION GNSS 2012), Nashville, TN, USA, 17–21 September 2012; pp. 161–173.
  31. Sharawi, M.S.; Rochester, M.I.; Akos, D.M.; Aloi, D.N. GPS C/N0 estimation in the presence of interference and limited quantization levels. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 227–238. [Google Scholar] [CrossRef]
Figure 1. Linear phase model of Doppler-aided PLL.
Figure 1. Linear phase model of Doppler-aided PLL.
Sensors 16 00519 g001
Figure 2. Feedback phase of traditionally aided DPLL.
Figure 2. Feedback phase of traditionally aided DPLL.
Sensors 16 00519 g002
Figure 3. Ratio of phase error amplitudes of aided loop to unaided loop.
Figure 3. Ratio of phase error amplitudes of aided loop to unaided loop.
Sensors 16 00519 g003
Figure 4. Feedback phase of sample-wisely aided DPLL.
Figure 4. Feedback phase of sample-wisely aided DPLL.
Sensors 16 00519 g004
Figure 5. Amplitudes of phase error transfer functions.
Figure 5. Amplitudes of phase error transfer functions.
Sensors 16 00519 g005
Figure 6. Block diagram of second-order PLL aided by linearly-interpolated Doppler.
Figure 6. Block diagram of second-order PLL aided by linearly-interpolated Doppler.
Sensors 16 00519 g006
Figure 7. Ratio of phase error of aided PLL to unaided PLL, with (a) ω = 1 rad / s ; (b) ω = 50 π rad / s .
Figure 7. Ratio of phase error of aided PLL to unaided PLL, with (a) ω = 1 rad / s ; (b) ω = 50 π rad / s .
Sensors 16 00519 g007
Figure 8. Carrier phase error using traditional aiding with 5 g dynamics and no thermal noise. (a) From unaided (0–30s) to aided (30–200s) status; (b) Aided period (30–200s).
Figure 8. Carrier phase error using traditional aiding with 5 g dynamics and no thermal noise. (a) From unaided (0–30s) to aided (30–200s) status; (b) Aided period (30–200s).
Sensors 16 00519 g008
Figure 9. Carrier phase errors using traditional aiding with different dynamics and no thermal noise. (a) Test results; and (b) Comparison with theoretical error amplitude.
Figure 9. Carrier phase errors using traditional aiding with different dynamics and no thermal noise. (a) Test results; and (b) Comparison with theoretical error amplitude.
Sensors 16 00519 g009
Figure 10. Carrier phase errors using sample-wise aiding with different dynamics and no thermal noise.
Figure 10. Carrier phase errors using sample-wise aiding with different dynamics and no thermal noise.
Sensors 16 00519 g010
Figure 11. Tracking errors with 5 g dynamics and 80 dB-Hz C/N0. (a) Carrier phase error; (b) Code phase error.
Figure 11. Tracking errors with 5 g dynamics and 80 dB-Hz C/N0. (a) Carrier phase error; (b) Code phase error.
Sensors 16 00519 g011
Figure 12. Tracking errors with 50 g dynamics and 80 dB-Hz C/N0. (a) Carrier phase error; (b) Code phase error.
Figure 12. Tracking errors with 50 g dynamics and 80 dB-Hz C/N0. (a) Carrier phase error; (b) Code phase error.
Sensors 16 00519 g012
Figure 13. Tracking errors with 500 g dynamics and 80 dB-Hz C/N0. (a) Carrier phase error; (b) Code phase error.
Figure 13. Tracking errors with 500 g dynamics and 80 dB-Hz C/N0. (a) Carrier phase error; (b) Code phase error.
Sensors 16 00519 g013
Figure 14. Comparison of tracking errors using traditional and sample-wise aiding.
Figure 14. Comparison of tracking errors using traditional and sample-wise aiding.
Sensors 16 00519 g014
Table 1. Four regions of sample-wise aiding in terms of Ta.
Table 1. Four regions of sample-wise aiding in terms of Ta.
Partitions of TaFeatures of Aiding Effect
T a < T l ( ω , R s ) C a ( ω T s ) 0.5 ω T s , ω T s determined, independent of T a and interpolation techniques
T l ( ω , R s ) T a < T m ( ω , R s ) decreasing as either T a or T s decreases
T m ( ω , R s ) T a < T u ( ω ) ω T a determined, independent of T s , improved by advanced interpolation technique significantly (at kind of exponential rate)
T a T u ( ω ) No much improvement even with smaller T s and advanced interpolation
Table 2. Test results of carrier phase errors using traditional Doppler aiding.
Table 2. Test results of carrier phase errors using traditional Doppler aiding.
Dynamics (Acceleration)Error Amplitude (m) for Unaided 3rd-Order PLL, B = 15 HzError Amplitude (m) for Aided 2nd-Order PLL, B = 15 HzRatio of Unaided to Aided Error Amplitudes
5 g0.003451.50 × 10−5230.00
20 g0.01386.03 × 10−5228.86
50 g0.03451.49 × 10−4231.54
100 gN/A3.00 × 10−4N/A
200 gN/A6.00 × 10−4N/A
500 gN/A1.50 × 10−3N/A

Share and Cite

MDPI and ACS Style

Kou, Y.; Zhang, H. Sample-Wise Aiding in GPS/INS Ultra-Tight Integration for High-Dynamic, High-Precision Tracking. Sensors 2016, 16, 519. https://doi.org/10.3390/s16040519

AMA Style

Kou Y, Zhang H. Sample-Wise Aiding in GPS/INS Ultra-Tight Integration for High-Dynamic, High-Precision Tracking. Sensors. 2016; 16(4):519. https://doi.org/10.3390/s16040519

Chicago/Turabian Style

Kou, Yanhong, and Han Zhang. 2016. "Sample-Wise Aiding in GPS/INS Ultra-Tight Integration for High-Dynamic, High-Precision Tracking" Sensors 16, no. 4: 519. https://doi.org/10.3390/s16040519

APA Style

Kou, Y., & Zhang, H. (2016). Sample-Wise Aiding in GPS/INS Ultra-Tight Integration for High-Dynamic, High-Precision Tracking. Sensors, 16(4), 519. https://doi.org/10.3390/s16040519

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