# The Unbiased Characteristic of Doppler Frequency in GNSS Antenna Array Processing

Authors: Yuchen Xie; Zhengrong Li; Feiqiang Chen; Huaming Chen; Feixue Wang
Journal: International Journal of Antennas and Propagation
 (2019)

---

## Abstract

The antenna array technology, especially the spaced-time array processing (STAP), is one of the effective methods used in Global Navigation Satellite System (GNSS) receivers to refrain the power of jamming and enhance the performance of receivers in the circumstance of interference. However, biases induced to the receiver because of many reasons, including characteristic of antennas, front-end channel electronics, and space-time filtering, are extremely harmful to the high precise positioning of receivers. Although plenty of works have been done to calibrate the antenna and to mitigate these biases, achieving a good performance of antijamming, high accuracy, and low complexity at the same time still remains challenging. Different from existing works, this paper leverages the characteristic of GNSS signal’s Doppler frequency in STAP, which is proven to remain unbiased to solve the problem, even when the nonideal antennas are used and the interference circumstance changes. Since the integration of frequency is carrier phase, the unbiased Doppler frequency leads to an accurate estimation of carrier phase which can be used to calibrate the antenna array without extra apparatus or complicating algorithms. Therefore, a simple Doppler-aid strategy may be developed in the future to solve the difficulty of STAP bias mitigation.


---

## Body

## 1. Introduction

Array processing is one of the most effective ways to refrain the jamming aimed at GNSS receiver, and among those array processing methods, the use of controlled reception pattern antenna (CRPA) arrays has attracted more and more attention. CRPA arrays can provide beam forming/null steering in specific direction by adaptively adjusting the weight of each antenna element. Typically, an adaptive finite impulse response (FIR) filter placed behind each element allows a better performance of antijamming, which is known as the space-time array processing [1]. CRPA arrays especially STAP-based antenna arrays are so popular that numerous adaptive algorithms have been researched to provide a high performance of interference suppression [2].Unfortunately, despite of the antijamming benefit that STAP can provide, undesirable and unpredictable biases are also induced to the GNSS receiver [2], which harms the accuracy of navigation solution [3–9]. There are several facts contributing to the STAP bias, including the overall gain and phase response of antennas, mutual coupling of antennas, front-end channel electronics, and space-time filtering. Previous research has focused on the code and carrier phase measurement bias using STAP. For space only processing (SOP), a code phase bias on the order of meters can be observed in simulation [3], and the phase bias of real data is discussed in [4]. The measurements are distorted further in STAP because of the applied FIR filters [5, 7]. In [1], the dependency of bias on interference circumstances is analyzed, which leads to the unpredictable characteristic of STAP-based bias.Therefore, to fulfill the ability of adaptive arrays, researchers have made great efforts to mitigate these biases. One solution is to use an adaptive antenna array that exhibits small biases [8], or to do the precalibration to correct the antenna based bias, but it may be impractical to some low cost receivers [10, 11]. Another way is to apply some specific algorithms that constrain the distortion of signal phase, but this method causes the loss of antijamming freedom degree. Some software based algorithms are also proposed to calibrate the bias [12, 13], but the computational complexity prevents them from the real time application. It seems that high accuracy, low complexity, and real time response are incompatible in the bias mitigation of STAP.Different from previous works, this paper analyzes the performance of Doppler frequency during the processing of antenna arrays. Logically, neither the response of antennas and channels nor the adaptive weights will change the carrier frequency of signal. In this paper, the unbiased characteristic of Doppler frequency of GNSS signal in STAP is proved by using the auto correlation function (ACF), and the simulation results also show that the changing of interference leads to the uncertain variation of phase, but it will not affect the estimation of Doppler frequency of receiver. When the response of channel is not ideal, the unpredictable distortion of phase is more obvious, while the estimation of Doppler frequency is still unbiased. Since the integration of frequency is carrier phase, the unbiased characteristic of Doppler frequency can be critical to high precision GNSS applications. Although the deduction is based on the GNSS signal, the unbiased characteristic is common in array processing. It means that we can correct the phase error induced by STAP or calibrate the antenna instantaneously using the Doppler frequency without any extra complexity of algorithm or hardware.The rest of this paper is organized as follows. First the array model is established and a brief description of the antenna induced bias in phase is given in Section2. Then the unbiased characteristic of the Doppler frequency of GNSS signal is proved in Section 3 with the help of ACF. In Section 4, simulation results are presented to compare the phase error and the frequency error, which further proves the unbiased characteristic of Doppler frequency. Finally, the conclusion is made in Section 5.


## 2. Array Model

Although STAP is an effective way to mitigate the interference in GNSS signal [14], biases induced by antennas and algorithms are not negligible for precision applications, especially when a large antenna array with complicated filtering is considered [3]. This section presents a model of array processing for describing errors induced by STAP and proving their dependency on interference circumstances.The antenna model is depicted in Figure1. In this model, K individual antenna element with an M-tap FIR filter following is considered. Ak(f,θ,φ) represents the system response of thekth element in (θ,φ) direction (θandφ stand for azimuth and elevation angles, respectively).Ak is also affected by the frequency of received signals (the effect of mutual coupling is not considered). Fk(f) stands for the electronic which downconverts the signal to baseband and perform analog-to-digital conversion. The output digital signal of each element is then filtered by an M-tap adaptive filter, and the frequency response of each FIR filter is denoted by Wk(f). Then the outputs of the filters are summed for postprocessing.Figure 1
STAP-based adaptive antenna array model.The complex weight in STAP filter can be represented by a stack vector(1)w=w1w2⋮wKwk is an M×1 vector corresponding to thekth filter(2)wk=wk1wk2⋯wkMTwhere wkm is each adaptive weight.The instantaneous digital output snapshot on the taps of thekth front-end channel is denoted by(3)xkn=xknxkn-1⋯xkn-M+1Txk[n] is the received signal on thekth element, which contains desired signal sk[n], undesired interference signal jk[n], and noise ηk[n].(4)xkn=skn+jkn+ηknThe snapshot on each filter is combined into a received signal stack vector(5)xn=x1nx2n⋮xKnThe STAP makes the weighted sum of signal vector by multiplying it with the weight vector(6)yn=wTxnBy applying some criterions [14], such as power minimization [15, 16], multiple constrained minimum variance [17], and minimum mean square error [18], the power of interference signal can be effectively restrained after the adaptive weighted summation.However, characteristics of antennas and channels, as well as STAP algorithms, cause biases to the output signal. To measure them, it is reasonable to analyze their effects on the GNSS receiver cross-correlation.The GNSS receiver correlators perform the cross-correlation by multiplying local C/A code replicasr[n] with the received signal s[n] and then doing the average of the correlation time Nr [19]. (7)Rτ=1Nr∑n=1Nrsn+τrnR(τ), named as the auto correlation function, is the output of the correlator, whereτ is the code delay between the received signal and the local replica.In STAP, the correlation process can be treated as the correlation between the signal of each FIR filter tap and the local C/A code replica, and then all of theseK×M ACFs sum up after multiplying their weight. (8)Rydτ=1Nr∑n=1Nryn+τrn=∑l=1KMalf,θ,φflfwlRτ+τl+j′+η′al(f,θ,φ), fl(f) and wl are the amplitude effects of antenna receiving pattern, channel response, and adaptive weight, respectively, while τl is the total code delay induced by antenna, channel, and weight for each tap. j′ and η′ are the remaining interference and noise after correlation.It can be known from (8) that the output ACF of STAP is a combination of K×M different ACFs whose amplitude and time delay vary from one to another. In fact, the antenna receiving pattern depends on the received signal’s direction and frequency, the channel response changes with time and temperature, and adaptive weights are also affected by interference circumstance. Therefore, the bias induced by STAP in ACF changes with the GNSS signal, the interference, and the environment, which is consequently unpredictable.Simulation results in different interference circumstances support this conclusion, which will be explained in detail in Section4. In this case, even if STAP has successfully refrained the power of interferences, the navigation output which is based on code phase measuring and carrier phase aid may not be accurate. However, it can be proved that Doppler frequency is unbiased after array processing, which can be further used to enhance the measuring accuracy of receiver. The details of deduction will be presented in the next section.


## 3. Doppler Frequency Estimation in STAP

In the single array receiver, the I∖Q orthogonal demodulation is applied to the received signal to move the carrier of it [19]. After orthogonal demodulation and correlation, the output of correlator can be denoted by(9)in=aDnRτcos⁡2πfen+ϕe(10)qn=aDnRτsin⁡2πfen+ϕea is the amplitude of received signal and D[n] is the data bit, both of which can be regarded as 1 for the sake of convenience. R(τ) is the ACF which has been defined in Section 2. fe and ϕe are the frequency and phase discrepancies between local carrier and received signal’s carrier, respectively.The relation among the local carrierflocal, the received signal’s carrier fcarr, the frequency error fe, the Doppler frequency of signal fd, and the standard carrier frequency f0 can be written as(11)fcarr=flocal+fe(12)fd=fcarr-f0=flocal+fe-f0Ideally, fe is zero, the Doppler frequency will exactly be the discrepancy between local generated carrier frequency and the standard frequency, i.e., flocal–f0. But in real situation fe is nonzero and includes the Doppler frequency estimation error and some other noise errors. Although they are difficult to be separated to get the exact Doppler frequency, in simulation test, with setting Doppler frequency and nonsignificant estimation error, the Doppler frequency can be reasonably estimated by calculating fe.Therefore, thefe calculating process of the GNSS receiver in its fine acquisition is firstly introduced. By simplifying (9) and (10) and doing the square, we get(13)i2n=R2τcos22πfen+ϕe(14)q2n=R2τsin22πfen+ϕeFurther, we set(15)zrn=i2n-q2n=R2τcos⁡4πfen+2ϕe(16)zin=2inqn=R2τsin⁡4πfen+2ϕeCombining zr[n] and zi[n] into a complex signal, we get(17)zn=zrn+jzin=R2τcos⁡4πfen+2ϕe+jR2τsin⁡4πfen+2ϕeIt is noticed thatz[n] is a single frequency complex signal with the amplitude of R2(τ), the frequency of 2fe, and the phase of 2ϕe. After doing the Fast Fourier Transform (FFT) toz[n], the maximum in its frequency domain is located at 2fe which is unrelated to its phase error 2ϕe.Therefore,fe can be achieved by(18)fe=12findmaxFFTznwhere findmax means searching for the frequency maximizing FFTzn.In the STAP receiver, the I∖Q orthogonal demodulation is also applied to the received signal, and derived from (8), (9), and (10), it can be rewritten as(19)iarrayn=∑l=1KMblRτ+τlcos⁡2πfen+ϕe+ϕl(20)qarrayn=∑l=1KMblRτ+τlsin⁡2πfen+ϕe+ϕlwhere(21)bl=alf,θ,φflfwlis the coefficient containing all amplitude effects. fe has been defined in (10), and jammer induced Doppler shift at the receiver is assumed to be small compared to satellite Doppler. τl has been defined in (8), and ϕl is the total carrier phase error induced by STAP. As mentioned before, bl, τl, and ϕl vary from one to another.Similarly, the square and multiplication of (19) and (20) are(22)iarray2n=∑l=1KM∑p=1KMblbpRτ+τlRτ+τp·cos⁡2πfen+ϕe+ϕl·cos⁡2πfen+ϕe+ϕp(23)qarray2n=∑l=1KM∑p=1KMblbpRτ+τlRτ+τp·sin⁡2πfen+ϕe+ϕl·sin⁡2πfen+ϕe+ϕp(24)iarrayqarrayn=∑l=1KM∑p=1KMblbpRτ+τlRτ+τp·cos⁡2πfen+ϕe+ϕl·sin⁡2πfen+ϕe+ϕpSo zr_array[n] and zi_array[n] can be denoted as(25)zr_arrayn=iarray2n-qarray2n=∑l=1KM∑p=1KMblbpRτ+τlRτ+τp·cos⁡4πfen+2ϕe+ϕl+ϕp(26)zi_arrayn=2iarraynqarrayn=∑l=1KM∑p=1KMblbpRτ+τlRτ+τp·sin⁡4πfen+2ϕe+ϕl+ϕpand(27)zarrayn=zr_arrayn+jzi_arrayn=∑l=1KM∑p=1KMblbpRτ+τlRτ+τp·cos⁡4πfen+2ϕe+ϕl+ϕp+j∑i=1KM∑j=1KMblbpRτ+τlRτ+τp·sin⁡4πfen+2ϕe+ϕl+ϕpTherefore, zarray[n] is the sum of KM×KM single frequency complex signals. Although they are different in amplitude and carrier phase, as(28)zarray_lpn=blbpRτ+τlRτ+τp(29)ϕarray_lp=2ϕe+ϕl+ϕpthey have the same frequency 2fe. In this case, (18) is still effective in estimating the Doppler frequency.Based on these deductions, it is reasonable to say that unlike code and carrier phase which will be unpredictably shifted because of the changing of interferences, the Doppler frequency of the received signal remains unbiased in array processing.In fact, even in the situation that the bias of phase is too severe to calibrate, or in the situation that the electronic characteristic of analog element changes, which makes the previous calibration ineffective, the Doppler frequency still remains unbiased because interference and STAP only affect the received signal’s phase rather than its frequency.Although the deduction is based on GNSS signal, it is also true to other signals of array processing. As the integration of frequency through time leads to phase, the phase error can also be corrected with the help of Doppler frequency. Therefore, the unbiased and accurate Doppler frequency is especially useful to high precision locating applications. In another way, the integration of the Doppler frequency can also be used to calibrate the bias induced by STAP, as the difference between the integration Doppler frequency and the output carrier phase is the total bias of STAP. In that case, the real time calibration for STAP can be realized to significantly enhance the performance of array processing.Simulation results in Section4 support the conclusion above and further prove the unbiased characteristic of Doppler frequency.


## 4. Simulation Results

In this section, simulations of the typical BeiDou receiver’s performance in different interference circumstances are presented to analyze the bias induced by STAP.The simulated receiver is a 4-element circular antenna array receiver, which has one antenna at the origin point with three others surrounded, and the distanced from the origin to each other antenna is half carrier wavelength. The structure of the antenna array is shown in Figure 2.Figure 2
Structure of antenna array.The STAP with an FIR filter back to each antenna is applied, and the power inverse (PI) [20] criterion is chosen to adapt the weight of each tap, which minimizes the output of antenna array processing to mitigate the effect of jammer. Meanwhile, the same signal received in an interference free circumstance and without applying any antijamming method is also processed by the receiver as the reference. The parameters for this simulation are listed in Table 1.Table 1
Simulation parameters.
Parameter Value Signal Type BeiDou Carrier Frequency 1268.52MHz Signal Length 1000 ms Intermediate Frequency 46.52MHz Sampling Frequency 61MHz Code Frequency 10.23MHz Code Length 10230 Signal Noise Ratio (SNR) -15dB Signal Direction (θ=85°, φ=70°) Jamming Noise Ratio (JNR) Jammer1: 50dB, Jammer2: 30dB Antenna Element 4 Array Structure Circular Time Taps 3, 5, 7 Anti Jamming Criterion PIIn the simulation, the circumstance of interferences changes with time and contains different types, directions, and powers of jamming. The simulation includes four steps as follows.Step 1.
Turn on the signal (θ=85°,φ=70°, SNR: -15dB).Step 2.
Turn on the jammer1 (θ=300°,φ=5°, JNR: 50dB) at the 200ms.Step 3.
Turn on the jammer2 (θ=135°,φ=10°, JNR: 30dB) at the 400ms.Step 4.
Change the direction of jammer1 (θ=180°,φ=30°) at the 600ms.Step 5.
Turn off the jammer1 and the jammer2 at the 800ms.
The jammer1 emits a white noise interference with the band of 20.46 MHz centered at 1268.52MHz, while the jammer2 emits a single frequency interference closed to the signal’s frequency. During the whole time of simulation, the signal is combined with white noise.


### 4.1. Ideal Channel Simulation

In the first simulation, we focus on the bias induced by adaptive algorithm; therefore ideal antennas and channels are considered. Figures3 and 4 show the simulation results.Figure 3 (a)
Code phase error in simulation1 (b)
Carrier phase error in simulation1Figure 4 (a)
Tracking loop frequency error in simulation1 (b)
Fine acquisition error in simulation1Figure3(a) is the difference of code phase between the signal after STAP and the reference signal, and Figure 3(b) is its counterpart carrier phase difference. It can be known from the figures that STAP successfully restrains the power of interferences and enables the receiver to keep tracking of the code phase. Nevertheless, even when the jammers are turned off, the code error is not zero because STAP induces bias to the receiver. Taking the 3-tap filter simulation in Figure 3(a) as example, the average error is -0.907m (1ms to 200ms) with the minimum of 0.799m at the 1ms, considering that one chip corresponds to 30m in our simulation. When the jammer1 is turned on and switched, the code phase error varies slightly, as the average error is -0.744m during 200ms to 400ms and -0.904m during 400ms to 600ms. However, when the jammer2 is turned on as well, the error of code phase sees a dramatic jump near the 600ms, after which it fluctuates severely, and the average error is -2.124m during 600ms to 800ms. The situations for 5- and 7-tap filter are similar but the errors are more severe than that of 3-tap filter.As for the carrier phase error, it can be known from Figure3(b) that it strongly depends on the circumstance of interference. In detail, when there is no interference, the carrier phase error fluctuates around zero in the 3-tap filter simulation, with a maximal absolute value of 1.003° at the 131ms. However, when turning on the jammers or switching the direction of jammer1, the carrier phase error jumps dramatically, as can be seen at the 200ms, 400ms, 600ms, and 800ms. Besides, when there exist interferences, the average of error is -16.227° (200ms to 400ms), -20.362° (400ms to 600ms), and 12.110° (600ms to 800ms), respectively, which presents an obvious bias from zero. It can also be noticed that the error fluctuates more drastically during the period from 600ms to 800ms when interferences are more complicated. The number of filter taps also influences the phase error but not significantly as can be seen from 400ms to 600ms in Figure 3(b).The average of code phase error and carrier phase error (denoted aserrcode and errcarr), as well as their standard deviation (denoted as σcode and σcarr), is shown in Tables 2(a) and 2(b) with the maximal value of each row being bold.Table 2  (a)
Average phase error in simulation1
e r r c o d e / e r r c a r r 0 - 200ms 200 – 400ms 400 - 600ms 600 – 800ms 800 – 1000ms total 3 taps -0.907m/-0.936° -0.744m/-16.227° -0.904m/-20.362° -2.124m/12.110° -0.877m/-0.266° -1.114m/-5.195° 5 taps -1.233 m/-0.864° -1.081m/-16.199° -1.210m/-18.508° -2.446m/15.585° -1.211m/-0.094° -1.439m/-4.064° 7 taps -1.247m/-0.773° -1.098m/-16.135° -1.216m/-16.948° -2.264m/17.682° -1.214m/-0.290° -1.410m/-3.329°  (b)
Standard deviation of phase in simulation1
σ c o d e / σ c a r r 0 - 200ms 200 – 400ms 400 - 600ms 600 – 800ms 800 – 1000ms total 3 taps 0.041m/0.382° 0.041m/0.650° 0.088m/0.479° 0.188m/1.080° 0.245m/0.550° 0.532m/12.527° 5 taps 0.036m/0.507° 0.036m/0.589° 0.086m/0.575° 0.200m/1.270° 0.228m/0.714° 0.530m/13.193° 7 taps 0.056m/0.660° 0.040m/0.667° 0.078m/0.610° 0.159m/1.499° 0.200m/0.817° 0.450m/13.543°Figures4(a) and 4(b) present two types of Doppler frequency errors in the simulation using different estimation methods. Data in Figure 4(a) is calculated from the output carrier frequency of the tracking loop while, in Figure 4(b), fine acquisition is applied to each 10ms signal to estimate the accurate Doppler frequency. It can be noticed in Figure 4(a) that as the tracking loop calculates the frequency by using carrier phase, when the interference changes, which causes dramatic jumps to the code phase at 200ms, 400ms, 600ms, and 800ms, the output frequency jumps consequently. However, after this variation, the output frequency returns back to its original value and the error fluctuates around zero; for example, the average error is -0.002Hz for 3-tap filter simulation. This can be further proved in Figure 4(b) where the accurate Doppler frequency is estimated; the errors are exact zero for different taps filter simulations during the whole simulation time.


### 4.2. Imperfect Channel Simulation

In the second simulation, imperfect antennas and channels are taken into consideration. The characteristics of channels are presented in Figures5(a) and 5(b), whose amplitude response waves randomly rang from -0.5dB to 0.5dB and phase response is nonlinear with a maximal shift of 30°. The other parameters and steps are exactly the same as those in the first simulation, and the results are shown in Figures 5–7.Figure 5 (a)
Channel amplitude response (b)
Channel phase responseFigure 6 (a)
Code phase error in simulation2 (b)
Carrier phase error in simulation2Figure 7 (a)
Tracking loop frequency error in simulation2 (b)
Fine acquisition error in simulation2Comparing Figure6(a) with Figure 3(a), it is obvious that the nonideal response of channels worsens the error of phase to vary more randomly; the gap between the maximum and minimum errors is about 5.039m. Similarly, the comparison between Figures 6(b) and 3(b) also shows a more drastic and random variation of the carrier phase, and the stable states of two pictures are different as well, which suggests new biases are induced because of channel characteristics.e r r c o d e, errcarr, σcode, and σcarr of the second simulation are shown in Table 3.Table 3  (a)
Average phase error in simulation2
e r r c o d e / e r r c a r r 0 - 200ms 200 – 400ms 400 - 600ms 600 – 800ms 800 – 1000ms total 3 taps 0.242m/96.960° 1.803m/-159.035° 1.645m/-156.356° 3.210m/ -83.518° 1.029m/1.545° 1.593m/-60.829° 5 taps 0.108m/68.139° 0.901m/-160.019° 0.817m/-158.376° 3.016m/ -91.127° 2.247m/16.234° 1.408m/-66.017° 7 taps 0.126m/32.475° 0.756m/-159.938° 0.706m/-158.293° 2.541m/ -94.550° 1.805m/46.786° 1.179m/-68.082°  (b)
Standard deviation of phase in simulation2
e r r c o d e / e r r c a r r 0 - 200ms 200 – 400ms 400 - 600ms 600 – 800ms 800 – 1000ms total 3 taps 1.062m/146.941° 0.407m/0.609° 0.096m/14.355° 0.455m/19.584° 2.183m/ 2.151° 1.593m/119.025° 5 taps 1.170m/162.459° 0.260m/0.723° 0.080m/13.630° 0.629m/134.521° 1.908m/ 178.308° 1.478m/142.145° 7 taps 1.193m/173.186° 0.227m/0.711° 0.079m/13.246° 0.552m/34.942° 2.121m/ 172.617° 1.403m/141.533°On the contrary, it can be figured out in Figures7(a) and 7(b) that the frequency error remains unbiased even with the consideration of channel effect. Although the deviation in Figure 7(a) is larger than that in Figure 4(a), the error returns to fluctuate around zero very soon, and the fine acquisition result in Figure 7(b) is the same zero as that in Figure 4(b).Based on the simulation results above, it can be concluded that STAP induces unpredictable bias into receivers, which causes errors in the estimation of code and carrier phase, and the situation is even worse when antennas and channels are nonideal. However, thanks to the unbiased characteristic of Doppler frequency, the estimation of frequency in our simulation remains stable no matter how the interference circumstance changes.


## 4.1. Ideal Channel Simulation

In the first simulation, we focus on the bias induced by adaptive algorithm; therefore ideal antennas and channels are considered. Figures3 and 4 show the simulation results.Figure 3 (a)
Code phase error in simulation1 (b)
Carrier phase error in simulation1Figure 4 (a)
Tracking loop frequency error in simulation1 (b)
Fine acquisition error in simulation1Figure3(a) is the difference of code phase between the signal after STAP and the reference signal, and Figure 3(b) is its counterpart carrier phase difference. It can be known from the figures that STAP successfully restrains the power of interferences and enables the receiver to keep tracking of the code phase. Nevertheless, even when the jammers are turned off, the code error is not zero because STAP induces bias to the receiver. Taking the 3-tap filter simulation in Figure 3(a) as example, the average error is -0.907m (1ms to 200ms) with the minimum of 0.799m at the 1ms, considering that one chip corresponds to 30m in our simulation. When the jammer1 is turned on and switched, the code phase error varies slightly, as the average error is -0.744m during 200ms to 400ms and -0.904m during 400ms to 600ms. However, when the jammer2 is turned on as well, the error of code phase sees a dramatic jump near the 600ms, after which it fluctuates severely, and the average error is -2.124m during 600ms to 800ms. The situations for 5- and 7-tap filter are similar but the errors are more severe than that of 3-tap filter.As for the carrier phase error, it can be known from Figure3(b) that it strongly depends on the circumstance of interference. In detail, when there is no interference, the carrier phase error fluctuates around zero in the 3-tap filter simulation, with a maximal absolute value of 1.003° at the 131ms. However, when turning on the jammers or switching the direction of jammer1, the carrier phase error jumps dramatically, as can be seen at the 200ms, 400ms, 600ms, and 800ms. Besides, when there exist interferences, the average of error is -16.227° (200ms to 400ms), -20.362° (400ms to 600ms), and 12.110° (600ms to 800ms), respectively, which presents an obvious bias from zero. It can also be noticed that the error fluctuates more drastically during the period from 600ms to 800ms when interferences are more complicated. The number of filter taps also influences the phase error but not significantly as can be seen from 400ms to 600ms in Figure 3(b).The average of code phase error and carrier phase error (denoted aserrcode and errcarr), as well as their standard deviation (denoted as σcode and σcarr), is shown in Tables 2(a) and 2(b) with the maximal value of each row being bold.Table 2  (a)
Average phase error in simulation1
e r r c o d e / e r r c a r r 0 - 200ms 200 – 400ms 400 - 600ms 600 – 800ms 800 – 1000ms total 3 taps -0.907m/-0.936° -0.744m/-16.227° -0.904m/-20.362° -2.124m/12.110° -0.877m/-0.266° -1.114m/-5.195° 5 taps -1.233 m/-0.864° -1.081m/-16.199° -1.210m/-18.508° -2.446m/15.585° -1.211m/-0.094° -1.439m/-4.064° 7 taps -1.247m/-0.773° -1.098m/-16.135° -1.216m/-16.948° -2.264m/17.682° -1.214m/-0.290° -1.410m/-3.329°  (b)
Standard deviation of phase in simulation1
σ c o d e / σ c a r r 0 - 200ms 200 – 400ms 400 - 600ms 600 – 800ms 800 – 1000ms total 3 taps 0.041m/0.382° 0.041m/0.650° 0.088m/0.479° 0.188m/1.080° 0.245m/0.550° 0.532m/12.527° 5 taps 0.036m/0.507° 0.036m/0.589° 0.086m/0.575° 0.200m/1.270° 0.228m/0.714° 0.530m/13.193° 7 taps 0.056m/0.660° 0.040m/0.667° 0.078m/0.610° 0.159m/1.499° 0.200m/0.817° 0.450m/13.543°Figures4(a) and 4(b) present two types of Doppler frequency errors in the simulation using different estimation methods. Data in Figure 4(a) is calculated from the output carrier frequency of the tracking loop while, in Figure 4(b), fine acquisition is applied to each 10ms signal to estimate the accurate Doppler frequency. It can be noticed in Figure 4(a) that as the tracking loop calculates the frequency by using carrier phase, when the interference changes, which causes dramatic jumps to the code phase at 200ms, 400ms, 600ms, and 800ms, the output frequency jumps consequently. However, after this variation, the output frequency returns back to its original value and the error fluctuates around zero; for example, the average error is -0.002Hz for 3-tap filter simulation. This can be further proved in Figure 4(b) where the accurate Doppler frequency is estimated; the errors are exact zero for different taps filter simulations during the whole simulation time.


## 4.2. Imperfect Channel Simulation

In the second simulation, imperfect antennas and channels are taken into consideration. The characteristics of channels are presented in Figures5(a) and 5(b), whose amplitude response waves randomly rang from -0.5dB to 0.5dB and phase response is nonlinear with a maximal shift of 30°. The other parameters and steps are exactly the same as those in the first simulation, and the results are shown in Figures 5–7.Figure 5 (a)
Channel amplitude response (b)
Channel phase responseFigure 6 (a)
Code phase error in simulation2 (b)
Carrier phase error in simulation2Figure 7 (a)
Tracking loop frequency error in simulation2 (b)
Fine acquisition error in simulation2Comparing Figure6(a) with Figure 3(a), it is obvious that the nonideal response of channels worsens the error of phase to vary more randomly; the gap between the maximum and minimum errors is about 5.039m. Similarly, the comparison between Figures 6(b) and 3(b) also shows a more drastic and random variation of the carrier phase, and the stable states of two pictures are different as well, which suggests new biases are induced because of channel characteristics.e r r c o d e, errcarr, σcode, and σcarr of the second simulation are shown in Table 3.Table 3  (a)
Average phase error in simulation2
e r r c o d e / e r r c a r r 0 - 200ms 200 – 400ms 400 - 600ms 600 – 800ms 800 – 1000ms total 3 taps 0.242m/96.960° 1.803m/-159.035° 1.645m/-156.356° 3.210m/ -83.518° 1.029m/1.545° 1.593m/-60.829° 5 taps 0.108m/68.139° 0.901m/-160.019° 0.817m/-158.376° 3.016m/ -91.127° 2.247m/16.234° 1.408m/-66.017° 7 taps 0.126m/32.475° 0.756m/-159.938° 0.706m/-158.293° 2.541m/ -94.550° 1.805m/46.786° 1.179m/-68.082°  (b)
Standard deviation of phase in simulation2
e r r c o d e / e r r c a r r 0 - 200ms 200 – 400ms 400 - 600ms 600 – 800ms 800 – 1000ms total 3 taps 1.062m/146.941° 0.407m/0.609° 0.096m/14.355° 0.455m/19.584° 2.183m/ 2.151° 1.593m/119.025° 5 taps 1.170m/162.459° 0.260m/0.723° 0.080m/13.630° 0.629m/134.521° 1.908m/ 178.308° 1.478m/142.145° 7 taps 1.193m/173.186° 0.227m/0.711° 0.079m/13.246° 0.552m/34.942° 2.121m/ 172.617° 1.403m/141.533°On the contrary, it can be figured out in Figures7(a) and 7(b) that the frequency error remains unbiased even with the consideration of channel effect. Although the deviation in Figure 7(a) is larger than that in Figure 4(a), the error returns to fluctuate around zero very soon, and the fine acquisition result in Figure 7(b) is the same zero as that in Figure 4(b).Based on the simulation results above, it can be concluded that STAP induces unpredictable bias into receivers, which causes errors in the estimation of code and carrier phase, and the situation is even worse when antennas and channels are nonideal. However, thanks to the unbiased characteristic of Doppler frequency, the estimation of frequency in our simulation remains stable no matter how the interference circumstance changes.


## 5. Conclusion

This paper analyzes the bias induced by STAP of the phase of the GNSS antenna array receiver and proves the unbiased characteristic of Doppler frequency of it. Simulation results show that the distortion of phase is unpredictable, and it will be even worse when the nonideal antennas are used or the interference circumstance changes. On the contrary, the Doppler frequency remains unbiased in these situations, which can be used to estimate an unbiased carrier phase to enhance the accuracy of positioning. Since a good-performance, low-complexity, and real-time bias mitigation is difficult to be realized by traditional methods, the Doppler-aid carrier phase correction may be a simple and effective way to achieve this goal.



---

*Source: 5302401-2019-04-24.xml*