Next Article in Journal
Objective Video Quality Assessment Based on Machine Learning for Underwater Scientific Applications
Previous Article in Journal
Wideband Spectrum Sensing Based on Riemannian Distance for Cognitive Radio Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Analytical Model for BDS B1 Spreading Code Self-Interference Evaluation Considering NH Code Effects

1
School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai 200240, China
2
Centre for Transport Studies, Imperial College London, London SW7 2BU, UK
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 13 February 2017 / Revised: 19 March 2017 / Accepted: 21 March 2017 / Published: 23 March 2017
(This article belongs to the Section Remote Sensors)

Abstract

:
The short spreading code used by the BeiDou Navigation Satellite System (BDS) B1-I or GPS Coarse/Acquistiion (C/A) can cause aggregately undesirable cross-correlation between signals within each single constellation. This GPS-to-GPS or BDS-to-BDS correlation is referred to as self-interference. A GPS C/A code self-interference model is extended to propose a self-interference model for BDS B1, taking into account the unique feature of the B1-I signal transmitted by BDS medium Earth orbit (MEO) and inclined geosynchronous orbit (IGSO) satellites—an extra Neumann-Hoffmann (NH) code. Currently there is no analytical model for BDS self-interference and a simple three parameter analytical model is proposed. The model is developed by calculating the spectral separation coefficient (SSC), converting SSC to equivalent white noise power level, and then using this to calculate effective carrier-to-noise density ratio. Cyclostationarity embedded in the signal offers the proposed model additional accuracy in predicting B1-I self-interference. Hardware simulator data are used to validate the model. Software simulator data are used to show the impact of self-interference on a typical BDS receiver including the finding that self-interference effect is most significant when the differential Doppler between desired and undesired signal is zero. Simulation results show the aggregate noise caused by just two undesirable spreading codes on a single desirable signal could lift the receiver noise floor by 3.83 dB under extreme C/N0 (carrier to noise density ratio) conditions (around 20 dB-Hz). This aggregate noise has the potential to increase code tracking standard deviation by 11.65 m under low C/N0 (15–19 dB-Hz) conditions and should therefore, be avoided for high-sensitivity applications. Although the findings refer to Beidou system, the principle weakness of the short codes illuminated here are valid for other satellite navigation systems.

1. Introduction

The BeiDou Navigation Satellite System (BDS) is a relatively new member of the Global Navigation Satellite System (GNSS) family. Its global role has recently been recognized by formal acceptance by the International Maritime Organization (IMO) as a key part of its World-Wide Radio-Navigation System (WWRNS). This is in addition to the adoption of 1575.42 MHz for the future B1c component (instead of 1561.098 MHz currently used by B1 signal) for interoperability with other GNSS’s. For both mass-market and professional GNSS users, BDS will provide improved accuracy, integrity, continuity and availability. There are several factors that determine a GNSS receiver’s tracking sensitivity, one of which is self-interference. This is due to cross-correlation between codes from other satellites, within the same GNSS constellation. Self-interference problem is present in GNSS with short (1023 or 2046 chips in one code period) spreading codes. The fact that this problem has little influence on existing GNSS such as GPS indicates it is not a major problem, but in some exceptional cases when signal is extremely weak (down to 20 dB-Hz or even weaker), this could cause trouble for receiver correlators, which requires close investigation. Take acquisition for example, the effects could be two-fold. First, the self-interference between a strong and a weak signal may surpass the weak signal autocorrelation peak. Second, the very same self-interference contributes to the pre-detection noise floor, effectively reducing the processing gain of a weak signal; this could amount to a 9 dB processing gain degradation in weak signal acquisition [1]. In addition, in reseach such as multiple access capacity of GNSS, self-interference is a non-negligible factor when considering the upper limit of navigation satellites in the sky [2]. All these examples necessitates evolution of existing self-interference evaluation methods.
Since GPS achieved full operational capability (FOC), research and development activities on the understanding and mitigation of self-interference were historically focused on GPS L1 C/A code. Typically, there are two methods to study this effect: simulation or an analytical way. Simulation can provide highly accurate results but at a cost in terms of computation time. A single location on earth surface takes about 11 CPU hours to gather the data for a whole day: 30-s data sets spaced every 15 min [3], and this translates into 13.752 CPU seconds for each second of simulated data. The accuracy of the simulation approach is 0.5 degrees (1 σ) for carrier tracking standard deviation, with coherent integration time equal to 0.02 s [4]. The analytical method is expected to achieve a lower but acceptable accuracy with a much reduced computation overhead. This paper intends to introduce a method for BDS self-interference evaluation on top of existing GPS self-interference models.
The analytical model of GPS L1 C/A self-interference has long been investigated. Self-interference when introduced as extra noise injected into a receiver, like cross correlation between different GNSS’s, is evaluated using spectrum separation coefficient (SSC) [5,6]. It is common practice that analytical models treat the short code as a simple wide sense stationary (WSS) random process using a single parameter auto-correlation function (ACF) [7]. However, this simple model for C/A code self-interference ignores the fact that the spreading code repeats every 1023 chips and treats C/A code as an infinite long sequence of random plus/minus ‘1’ s. Efforts to improve the accuracy of the analytical model were made by Dierendonck and Hegarty [8], Hegarty and Dierendonck [9], Shibata and Maeda [10] and Dierendonck et al. [11]. More complex WSS models were derived to predict C/A code self-interference by considering the real L1 C/A signal comprised of 20 repetitions (which was previously ignored) of spreading code of 1ms long in one navigation data bit. In fact, the short spreading codes including the GPS L1 C/A code (and also BDS B1-I code), are more accurately modelled as a cyclostationary random process using two-parameter ACF [12]. Although simple, the existing analytical methods generate results worse than (typically with accuracy of 0.2 degrees (1 σ) for carrier tracking standard deviation over a 10-s sliding window) or even different from those from pure simulation methods using a constellation simulator, a software receiver, and a software signal generator [13]. In addition, to date all analytical models on self-interference are for GPS L1 C/A code, and there are no models for self-interference within the BDS constellation.
To evaluate BDS self-interference, a simple analytical model is developed from an existing model for GPS self-interference [12] to predict the impact of B1-I code self-interference on BDS receiver functions that are dependent primarily on the sum of the correlations (e.g., carrier phase tracking and data demodulation). The focus is on self-interference between BDS MEO and IGSO satellites. Self-interference effects between BDS GEOs are not studied because GEO signals do not interfere with MEO signals due to the high symbol rate [14]. The effects of the NH code [15], only present in the BDS MEO and IGSO B1-I signal, on self-interference evaluation, is addressed during the development of the analytical model.
The proposed model is based on the perspective that self-interference can be modelled by the introduction of an extra white noise in receiver operations. The noise has the same effects on code and carrier tracking as self-interference. The accuracy of the model is ensured by using a two-parameter ACF to characterize a B1-I code either with or without navigation data. This type of cyclostationary analytical model will be validated as an accurate method to evaluate self-interference for short codes, producing excellent agreement with actual receiver observables, using hardware simulator intermediate frequency (IF) data and a software receiver.
The Section “Signal and Correlator Model” provides basic assumptions and signal models to be used in the development of the proposed analytical model. Section “SSC and Equivalent White Noise Level” derives analytical models for the variance at the correlator output for three situations: the data bits of the desired and undesired signals are (1) aligned or (2) misaligned by an integer number of code chips or (3) misaligned by a fraction of a code chip. These three models are the basis of the proposed analytical model. Section “Proposed Model for BDS MEO/IGSO B1-I Self-Interference” presents the proposed model based on SSCs derived in the previous section. Section “Model Validation” demonstrates the validity of the model using live BDS data record and a modified software receiver. Section “Results” analyzes typical behavior of B1-I self-interference and its impact on receiver range observables using the developed model. Section “Conclusions” summarizes the expected contributions and presents a framework for further efforts.

2. Signal and Correlator Model

Figure 1 shows a generic correlator model, whose input is the undesired (interfering) signal, x1(k)(t − Δk)ej2πf, and the desired (victim) signal is x2(t). There is a differential Doppler frequency, i.e., difference between Doppler frequencies, denoted as f, between the undesired and desired signals. These two signals are correlated by mixing and integration over an arbitrary navigation data bit period, [kTb, (k + 1)Tb), where Tb = 20 ms and k = 0, 1, 2, ....
To simplify the problem without loss of generality, we focus on an individual undesired signal including a specific spreading code identified by PRN (Pseudo Random Noise) number, x1(t). This undesired signal can be any one of the undesired signals x1(k)(t) in the input to the correlator in Figure 1. Therefore, the undesired signal without Doppler or time delay is modelled as:
x 1 ( t ) = k = d k u = 0 19 α u v = 0 2045 c v p T c ( t [ 40 , 920 k + 2046 u + v ] T c )
where pTc(t) is a rectangular pulse of width Tc:
p T c ( t ) = { 1 0 t < T c 0 e l s e
and cv ∈ {−1, 1} is spreading code chips with a chipping rate 1/Tc = 2.046 MHz. The chip width is Tc; dk is navigation data bits with bit rate 1/Tb = 50 bps, and the bit width is Tb. αu is the NH code adopted for B1-I signal and takes a fixed pattern of (0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 1, 0, 0, 1, 1, 1, 0), which spans one navigation data bit. The spreading code here is treated as a random periodic sequence of length 2046 chips, with chips independent from each other. The navigation data bit stream is also treated as a random periodic sequence, but with an infinite length.
In order to facilitate the development in the appendices, spreading code period, T, is introduced, and there is an obvious relationship among Tc, T, and Tb:Tb = 20T = 40,920Tc. This relationship is illustrated in Figure 2. The receiver local code replica or desired signal is modelled similarly as:
x 2 ( t ) = k = l = 0 2045 c k p T c ( t [ 2046 k + l ] T c )
using the same notation as Equation (1). The only difference between the desired and undesired code is that the undesired one has a time delay of ∆ seconds, which is normally within the range of ±20 ms for BDS users on or near the horizon of the Earth. Both signals are cyclostationary, not WSS, and therefore, their characteristics are accurately represented using a two-parameter ACF [16]. Since the period of cyclostationarity is 20 ms, the model derived in the following sections only takes into account time delays ranging from 0 to 20 ms, i.e., 0 ≤ ∆ ≤ 20 ms, which is equal to the time delay modulo 20 ms.
In a real receiver, both the incoming signal and local replica have a Doppler shift. However, for simplicity and without loss of generality, we assume here that only the incoming signal has its Doppler shift, f, while the local replica has zero Doppler. This is simply because the receiver carrier tracking loop is after all, tracking the differential Doppler between the incoming and the local replica, not the absolute Doppler of each signal.

3. SSC and Equivalent White Noise Level

The extent of the undesired signal’s impact on the desired signal is characterized by an equivalent white noise level, which is derived from SSC, whose analytical forms are derived in this section. The k-th correlator output, yk, is a random variable. Its mean value is:
E { y k } = k T b ( k + 1 ) T b E { x 1 ( t Δ ) } E { x 2 ( t ) } e j 2 π f t d t = 0
and its variance is:
E { | y k | 2 } = E { y k y k } = E { k T b ( k + 1 ) T b k T b ( k + 1 ) T b x 1 ( s Δ ) x 1 ( t Δ ) x 2 ( s ) x 2 ( t ) e j 2 π f s e j 2 π f t d s d t } = k T b ( k + 1 ) T b k T b ( k + 1 ) T b E { x 1 ( s Δ ) x 1 ( t Δ ) } E { x 2 ( s ) x 2 ( t ) } e j 2 π f s e j 2 π f t d s d t = k T b ( k + 1 ) T b k T b ( k + 1 ) T b R 11 ( s Δ , t Δ ) R 22 ( s , t ) e j 2 π f s e j 2 π f t d s d t
where the dot within the curly brace means “dot product”. The development of Equation (5) uses the definition of auto-correlation function of two complex signals. The two-parameter ACF’s of both the undesired signal and receiver local replica are derived in Appendix A. The result is as follows:
R 22 ( s , t ) = p = p T c ( t 2046 p T c v T c )
where R22(s, t) is the two-parameter ACF of the local replica, the integer v is selected as:
v = mod ( floor ( s / T c ) , 2046 )
and:
R 11 ( s Δ , t Δ ) = p = 0 19 p T c ( t Δ [ 40 , 920 k + 2046 p + v ] T c )
where R11(s, t) is the two-parameter ACF of the undesired signal, ∆ is the time delay of the undesired signal and the integer k and v are selected as:
k = floor ( ( s Δ ) / T b )
v = mod ( floor ( ( s Δ ) / T c ) , 2046 )
where function mod(a, m) returns the remainder after division of a by m where a is the dividend and m is the divisor. Function floor(x) rounds x to the nearest integer less than or equal to x.
The development of Equations (8)–(10) has taken into account the special structure of the B1-I signal which has an NH code modulated on it: in subsection “Desired signal ACF” of Appendix A it can be seen that since the NH code is a fixed pattern finite-length code, it is treated as a deterministic variable here.
Figure 3 illustrates an example of ACF for the desired signal R22(s, t) with s ∈ [0, Tc). For every fixed value of the first parameter, s, this two parameter ACF is composed of 20 pulses of width Tc with unity amplitude, over arbitrary 20-ms correlation interval, which is equal to one data bit width. For example, if s falls within the time duration of the repetition of the first chip of the spreading code sequence, the ACF takes on a value of unity during each repetition of this chip and zero elsewhere. The same ACF results if s fell within repetition of the first chip, e.g., s ∈ [7T, 7T+Tc). The ACF for the undesired signal can be obtained by shifting the desired signal’s ACF by ∆.
The rest of this section is divided into three subsections. Analytical model for variance at the correlator output under three circumstances are derived: (i) data bits of the undesired and desired signal are aligned, (ii) data bits of these two signals are not aligned by an integer multiples of a code chip, and (iii) data bits of these two signals are misaligned by an arbitrary delay. This variance is then converted to the equivalent white noise level using SSC, which is the metric to measure how the self-interference affects receiver performance [17].

3.1. Data Bits Aligned

When the data bits modulated on the undesired and desired signal are aligned, i.e., ∆ = 0, the situation is illustrated in Figure 4.
In this case, ACF’s of the undesired and desired signals become identical over any arbitrary integration interval. The variance of the correlator output could be represented as:
E { | y k | 2 } = k T b ( k + 1 ) T b k T b ( k + 1 ) T b R 11 ( s , t ) R 22 ( s , t ) e j 2 π f s e j 2 π f t d s d t = k T b ( k + 1 ) T b k T b ( k + 1 ) T b R 22 ( s , t ) e j 2 π f s e j 2 π f t d s d t
the last line of which comes from the fact that the ACF’s of the incoming and local code are identical: both of them are unity pulse trains. As derived in Appendix B, substituting Equation (6) into Equation (11), we have:
E { | y k | 2 } = T T c [ sin ( π f T c ) ( π f T c ) sin ( π f T b ) ( π f T ) ] 2 = ( T b ) 2 T c T [ sin c ( f T c ) sin c ( f T b ) ] 2
where:
sin c ( t ) { 1 t = 0 sin ( π t ) π t t 0
is the sampling function. On the other hand, assume that the correlator shown in Figure 1 was driven by additive Gaussian white noise with power spectral density I0, then its output should have zero mean and a variance of:
E { | y k | 2 } I 0 T b
which is the product of the power spectral density and equivalent bandwidth 1/Tb. Comparing Equations (12) and (14), and an equivalent white noise with power spectral density I0 is obtained defined as:
I 0 = E { | y k | 2 } T b = T c T T b [ sin c ( f T c ) sin c ( f T b ) ] 2
because this equivalent white noise generates the same variance values as a true interference does at the output of the correlator. Based on this definition, we can further define the spectral separation coefficient [17] for B1-I as SSCB1-I:
I 0 P R T c T T b [ sin c ( f T c ) sin c ( f T b ) ] 2 P R S S C B 1 I
where, PR is inserted for the more general case where an undesired signal has an arbitrary power level in watts.

3.2. Data Bits Misaligned by an Integer Number of Code Chips

When the data bits modulated on the undesired and desired signals are aligned, i.e., ∆ is not zero, the situation is illustrated in Figure 5.
In this case, the ACF of the incoming/undesired signal is:
R 11 ( s Δ , t Δ ) R 22 ( s , t ) = { p = 0 19 p T c ( t Δ ( 40 , 920 k + 2046 p + v ) T c ) 0 s < Δ p = 0 19 p T c ( t Δ ( 2046 p + v ) T c ) Δ s < T b
where k and v is set up using Equations (7), (9) and (10). As derived in Appendix B, assuming that the data bit misalignment is an integer multiple of the spreading code period T, i.e., ∆ = KT for 0 K < 20, where K is an integer, results in:
S S C B 1 I = T T c T b [ sin c ( f T c ) ] 2 [ sin 2 ( π f K T ) + sin 2 ( π f ( 20 K ) T ) sin 2 ( π f T ) ]
which, for users near or at the Earth’s horizon, could be approximated as:
S S C B 1 I T T c T b [ sin 2 ( π f K T ) + sin 2 ( π f ( 20 K ) T ) sin 2 ( π f T ) ]
More generally, assuming that the data bit misalignment is an integer multiple of the spreading code chip width Tc, i.e., ∆ = KT + CTc, where 0 K < 20, 0 < C < 2046, K and C is an integer, yields:
S S C B 1 I = 2046 C 2046 S S C B 1 I , K + C 2046 S S C B 1 I , K + 1
where SSCB1-I,K is the SSC in Equation (18) for any users or SSC in Equation (19) for users near or at the Earth’s horizon, when ∆ = KT.

3.3. Data Bits Misaligned by a Fraction of a Code Chip

Section 3.2Data Bits Misaligned by an Integer Number of Code Chips” summarizes SSC expressions that is valid when the differential time delay between undesired and desired signals, ∆ is multiples of a single code chip Tc. With this constraint, R11(s − ∆, t − ∆)R22(s,t) is a pulse train of width Tc, given by Equation (17).
Now we take a step further to examine what the SSC should look like when ∆ is only a fraction of a code chip. In this case, R11(s − ∆, t − ∆)R22(s,t) is still a pulse train of uniformed width, but the pulse width is no longer equal to Tc. As illustrated in Figure 6 and Figure 7, for fixed s, R11(s − ∆, t − ∆)R22(s,t) is either a train of pulses of width ∆ or Tc − ∆.
When the differential time delay takes on arbitrary value, we proceed by using τ = mod(∆, Tc), 0 τ < Tc. Therefore, R11(s − ∆, t − ∆)R22(s,t) will be a train of pulses of width τ for τ/Tc percent of Tc, or a train of pulses of width Tcτ for (Tcτ)/Tc percent of Tc. The total effect of such arbitrary time delay, for the differential Doppler, is that a factor of Tc in the SSCB1-I expressions in the Section 3.2Data Bits Misaligned by an Integer Number of Code Chips” is replaced by [12]:
τ T c τ + ( T c τ ) T c ( T c τ )
which implies the expectation of receiving such a time difference. It can be seen from Equation (21) that spreading code self-interference is maximized when τ = 0, i.e., when the misalignment between the desired and undesired signals happens to be multiples of a code chip; it can be seen also that self-interference is minimized then τ = Tc/2. Moreover, the average of Equation (21) can be computed as:
1 T c 0 T c [ τ 2 T c + ( T c τ ) 2 T c ] d τ = 2 3 T c
which, when substituted for Tc in Equation (19) yields:
S S C B 1 I 2 3 T T c T b [ sin 2 ( π f K T ) + sin 2 ( π f ( 20 K ) T ) sin 2 ( π f T ) ]
which is used to compute SSC in the following sections.

4. Proposed Model for BDS MEO/IGSO B1-I Self-Interference

Based on Equations (20) and (23), the proposed model for BDS MEO/IGSO B1-I self-interference, or the aggregate equivalent white noise due to self-interference can be expressed as:
I 0 , m = p m + i = 1 i m n ( S S C B 1 I ) m i = p m + 2 T c 2 3 T b sin 2 ( π f m i T ) i = 1 i m n { ( 2046 C m i ) [ sin 2 ( π f m i K m i T ) + sin 2 ( π f m i [ 20 K m i ] T ) ] + C m i [ sin 2 ( π f m i [ K m i + 1 ] T ) + sin 2 ( π f m i [ 19 K m i ] T ) ] }
where I0,m is the aggregate equivalent white noise imposed on the m-th signal by self-interference, pm is the received power at the antenna, (SSCB1-I)im is the SSC between the i-th and m-th signals, which is determined by differential Doppler fim and differential time delay Kim T + Cim Tc.
Table 1 shows an example of how to use the developed analytical model, i.e., Equations (20) and (23) to predict SSCB1-I and equivalent white noise I0.
Note that the last line of Table 1 comes from the sum of I0 of PRN 8 and 9. At a specific time, there are 8 satellites visible to the receiver. In Table 1, the first column is the PRN number. The second column is the relative received power, relative to −160 dBW. The third column is the true (geometric) range from satellite to the receiver. The fourth column is the true Doppler. These four columns are the parameters used to compute the fifth and sixth columns, i.e., SSCB1-I and equivalent white noise I0, respectively. Since here it is assumed that PRN 6 is the desired signal, it is unnecessary to compute its own SSC and I0.
Table 2 presents the procedure to compute SSC and I0. The procedure described in Table 2 is simply a re-iteration of the development of the model Equation (24).

5. Model Validation

Live BDS signal records from a hardware signal simulator are collected to test if the proposed model is correct within a certain boundary, set by established benchmark work.
The validity of the model is not directly proved because the proposed model includes the NH code effect and there is no current work covering this unique signal component on B1 self-interference evaluation, though there does exist effective ways to overcome NH code impact on various stages of signal processing [15]. Since no other formula is currently available, an indirect method is applied to justify the assumptions and the proposed model. Self-interference effect on receiver noise is intangible or indirect to see or to feel but can appropriately be shown and visualized by its impact in the measurement domain (i.e., DLL output or a step further, code range), and in this regard, we choose to apply Betz’s well-established expressions linking C/N0 with discriminator output statistics. Since no validation of Betz’s model was provided in Betz’s paper [18] or his later publications, and therefore in order to use Betz’s formula, we had first to prove it (or at least, show its effectiveness using real measurements), which is presented in this section and does form an integral part of the completeness and soundness of the validation logic.

5.1. Criterion

Model validity is checked by comparing the predicted and measured tracking standard deviation, using the expression [18]:
var { τ k u | τ k s } = var { e ( ε ) | τ k s } C s 4 K 2
where K is code tracking loop gain and Cs the measured signal power. var{e(ε) | τsk} is measured discriminator output variance. var{τuk | τsk} is variance of the unsmoothed code delay estimate τuk, which is predicted by using Equation (26) [19]:
var { τ k u | τ k s } = B L B r / 2 B r / 2 G s ( f ) G w ( f ) sin 2 ( π f Δ ) d f ( 2 π ) 2 C s ( B r / 2 B r / 2 f G s ( f ) sin ( π f Δ ) d f ) 2 [ 1 + 1 T ( C s N 0 ) e f f B r / 2 B r / 2 G s ( f ) d f ]
which is the analytical form of code tracking variance. The term ‘(Cs/N0)eff’ is effective carrier to noise ratio, and is equal to CsI0, where Cs is the measured signal power and I0 the equivalent white noise as defined in Equation (16) and transformed into dB. Here I0 includes self-interference effects, to be evaluated using the proposed model (Equations (16) and (20)). Both Cs and I0 are in dB and (Cs/N0)eff should be converted back to non-logarithmic unit when it is used in Equation (26). If the model is correct, the squared root of the two sides of Equation (25) should agree within a small error tolerance for a specific channel in a software receiver:
| var { τ k u | τ k s } var { e ( ε ) | τ k s } C s 4 K 2 | < λ
where λ is a tolerance value and we set it as follows. We compare our prediction accuracy with accuracy of most recent analytical methods for self-interference, such as the work by Cerruti et al. [4], where the accuracy of prediction, i.e., the difference between predicted and measured carrier tracking standard deviation is 0.5 degrees. This translates into a prediction accuracy of 0.0264 cm (0.5/360 × 19) in range by phase measurements, which is in turn, translated into code tracking prediction accuracy of 6.6 cm by a scale factor of 250, since code tracking loop is 100 (with multipath) to 250 (without multipath; this is the case with current validation settings described in Table 3, and is also the case with the benchmark work by Cerruti et al. [3]) times noisier than carrier tracking loop in terms of tracking standard deviation [20]. Therefore, we set λ = 0.066 m.

5.2. Test Setup

The test setup is illustrated in Figure 8. Since we want to compare predicted and measured tracking standard deviations, we have to collect live signal records to measure the code tracking standard deviation using a software receiver. Without loss of generality, we choose two IGSO satellites, PRN 13 and 14 instead of the whole current constellation to test the developed model. This is because: (i) In validation, only two satellites (desired and undesired signal) are necessary; (ii) PRN 13 & 14 are the two of the four currently operational IGSO satellites in BDS, and since MEO and IGSO in BDS use the same signal modulation scheme; IGSO signal will suffice to represent a combination of MEO and IGSO signals. Therefore, first we collect IF data containing both PRN 13 and PRN 14 signals and then we collect IF data containing only the PRN 13 signal. Finally, we obtain the measured code tracking standard deviation of PRN 13 using the first dataset, which includes the interference effect imposed by PRN 14. This interference effect is predicted using the estimated power (Cs) from the second dataset (free of interfering signal PRN 14) and our proposed model. As a result, a hardware signal simulator and not live signals from satellites, must be used as the signal source.
We use a HWA-RNSS-7400, which is a multi-constellation GNSS signal simulator by HWA Create Co., Ltd. (Beijing, China) [21]. It features BDS B1, B2, and B3 signal simulation and is widely used within BDS industries. We also use a SAS6862A, which is a wide-band, dual channel (GPS L1+BDS B1 and BDS B3) radio frequency front end (RFFE) with a highly stable (better than 0.005 ppm) oven controlled crystal oscillator (OCXO) [22]. We use a multi-constellation software receiver called UNIversal COmmunications and Radio Navigation receiver (UNICORN), developed by the authors at the Centre for Transport Studies (CTS), Imperial College London and used to obtain the ESA’s certificate for the “First 50 Successful Galileo Position Fixes”. The receiver has been adapted to process the B1 signal.
Parameters used to configure UNICORN and drive the analytical model are the same and are provided in Table 3. The simulation starts at UTC 00:10:00, 27 Feb 27, 2014 and two 37-s IF data was recorded using SAS6862A. The receiver is assumed to be located at 31° N, 121.5° E, at an altitude of 5 m in Shanghai, China.
UNICORN has to produce two measurements for the analytical model: C/N0 and code loop gain, which is indispensable for successful evaluation of tracking error variance through the model. Within UNICORN, several lines of code were inserted to calculate tracking error standard deviation.

5.3. Results

Finally, the results from UNICORN and the analytical model are compared to see if the predicted values agree with measured values according to Equation (27). Figure 9 shows the results.
In Figure 9, it can be observed that the discrepancy between measured and predicted standard deviations are bounded by 0.0002 B1 primary code chips, which translates into a prediction accuracy of 0.0293 m, less than half of the threshold value of 0.66 m. This proves the validity of the analytical model.

6. Results

Based on the validity of the model, we now apply this model to show some insights into the typical behaviour of self-interference in BDS B1-I signals. Firstly, we try to find when the equivalent noise introduced by the self-interference approaches its maximum value.
The SSC predicted from Equation (23) is illustrated in Figure 10 and Figure 11. Here the differential time delay between the desired and undesired signals is an integer multiple of spreading code periods. Since the SSC values for K > 10 is equal to SSC values with 20 − K, only cases with K 10 are shown here. To facilitate viewing, the results are shown separately when K is odd or even. From Figure 10 and Figure 11 we can conclude that when the differential Doppler between desired and undesired signals is bounded by 50 Hz, self-interference will introduce the most serious noise into receiver processing chain.
Secondly, in order to analyze the self-interference effects in current BDS configuration, a simulation case is setup. The orbit propagation, receiver acquisition and tracking are performed by the Spirent SimGEN simulator software suite. SimGEN is capable of full and versatile case generation. Users are allowed to control multiple GNSS and regional satellite constellations including BDS, signal propagation, multipath and obstruction effects, antenna patterns, vehicle trajectory and various error models, thus vastly facilitating our simulation. For the current simulation, three inputs are needed from SimGEN, namely true range, received power, and carrier Doppler shift. However, due to the version of SimGEN, it is not able to simulate BDS constellation. The problem is bypassed by using Two Line Element (TLE) file retrieved from NORAD [23] and feeding it to the Motion control panel for GPS in SimGEN. The simulation length is 20 minutes, starting from UTCG1927, 26 April 2015, with a measurement rate of 1 Pulse Per Second (PPS). The receiver is placed in Shanghai, China. The result is shown in Figure 12. The constellation information is shown in Appendix C.
As illustrated in Figure 12, we can obtain such observations:
(1)
Within most of the simulation interval, the individual equivalent noise and aggregate noise is below −215 dBW/Hz, which is much lower than a typical receiver noise floor of about −201.5 dBW/Hz, and therefore would not cause much trouble for most applications. Here trouble means worsening of such receiver capabilities as acquisition sensitivity and tracking jitter.
(2)
During 200 and 240 s, the equivalent noise of PRN 8 and the aggregate noise reaches −200 dBW/Hz. For a typical value of receiver noise floor (−201.5 dBW/Hz), this noise floor will be transferred to −197.67 dBW/Hz. In this case, PRN 6’s effective carrier to noise density ratio (Cs/N0)eff is 3.83 dB lower than the original carrier to noise density ratio, should the self-interference caused by PRN 8 be absent. This could pose a serious problem to high-sensitivity processing during acquisition, tracking or navigation bit modulation. The potential hazard caused by this (C/N0)eff decrease could be best exemplified by Figure 13, which shows the code tracking standard deviation versus (Cs/N0)eff , using established analytical method by Zhang and Zhan [19]. The B1-I receiver used in Figure 13 is assumed to have a precorrelation bandwidth of 4 MHz, a sampling frequency of 16 MHz, and 2-bit quantization. The coherent integration interval is 20 ms and a Dot Product discriminator is used. When (C/N0)eff decreases from 19 to 15 dB-Hz, tracking standard deviation increases from 11.07 m to 22.72 m. In this regard, B1-I self-interference could pose a potential hazard on pseudorange measurement for high-sensitivity receivers. However, it must be underlined that for most receivers (from mass market to medium-grade receiver), a (C/N0)eff value under 20 dB-Hz, though may be possible, is an exceptional case. When (C/N0)eff decreases from 24 to 20 dB-Hz, tracking standard deviation increases from 5 m to 9 m, which does not impose significant impact on final positioning accuracy for most low- to medium-grade receivers.
(3)
Around 240 s for PRN 8 and around 1200 s for PRN 9, the self-interference effect on the desired signal PRN 6 reaches its maximal value for each undesired signal respectively. This corresponds to the time when these two PRN’s Doppler shift (modulo 1 kHz) are aligned (zero differential Doppler shift) with the desired signal, PRN 6. The same situation has also been found in GPS C/A-to-C/A self-interference [4].
In this regard, the proposed model could be used in B1 receiver as a code-tracking performance indicator, since code tracking standard deviation is directly determined by (Cs/N0)eff., the effective carrier-to-noise power density ratio, which is in turn, directly affected by self-interference effects, as indicated in Figure 13. In the simplest case, the common satellite selection step (e.g., depending on DOP, Dilution of Precision) should be complemented by a new step to narrow down the list of satellites participating in navigation solution, based on the result of self-interference calculation according to Equations (20) and (23), as suggested in Figure 14.
In Figure 14, a typical receiver baseband (correlator + measurements) plus navigation module is shown, aided by a satellite selection module, whose action depends on a Signal Quality Monitoring (SQM) module [24] on top of its normal operations, monitoring the equivalent noise spectral density imposed by self-interference,. This may be an indication to receiver manufacturers for their design of strategies for satellite selection in Position, Velocity and Time (PVT) computation and signal acquisition in case of very low C/N0, which could be estimated stably and accurately through novel estimators [25].

7. Conclusions

A simple analytical model is developed to predict the impact of B1-I code self-interference on BDS receiver functions that are dependent primarily on the sum of the correlations (e.g., carrier phase tracking and data demodulation). The model is based on the fact that the code is actually a cyclostationary process rather than WSS, as previously assumed in SSC evaluation. A two-parameter ACF, in contrast to the usual one-paramter ACF, is applied to characterize B1-I code either with or without navigation data. The assumption of cyclostationarity is essential in evaluating self-interference effects for short codes, producing excellent agreement with actual receiver observables. The model is found to be able to accurately predict self-interference effects on receiver code tracking standard deviation, when the difference of time delay and Doppler between desired and undesired signals varies.
Even during the days when BDS is evolving into its Phase III stage, there is a large user group still using the legacy BDS signal. The proposed model could be used in B1 receivers as a code-tracking performance indicator and an aid to satellite selection during navigation solution.

Acknowledgments

This work is jointly supported by National Natural Science Foundation of China (No. 61403253) and the International Postdoctoral Exchange Fellowship from China Postdoctoral Council (No. 53, 2013), which supported the primary author’s two-year postdoctoral research at Centre for Transport Studies, Imperial College London.

Author Contributions

Xin Zhang conceived and designed the experiments, performed the experiments, analyzed the data, and wrote the paper. Shaojun Feng helped conceive the experiments and restructure several parts after Xin Zhang completed the first edition of the manuscript. Xingqun Zhan helped improve the abstract, introduction, and conclusions. Both Feng and Zhan provided multiple proof readings of the latest version of the manuscript. Shaojun Feng and Washington Ochieng provided the software simulation tool and Xingqun Zhan provided the hardware tools. Washington Ochieng offered final proof reading of the manuscript. All revisions were committed by Xin Zhang.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Autocorrelation Correlation Function Derivations

First, we will derive ACF for the undesired signal. Since the signal is assumed cyclostationary, a two-parameter ACF instead of one-parameter (a single time difference between a stationary signal and its lagged version) ACF is developed here.

Appendix A.1. Undesired Signal ACF

We start from the definition of two-parameter ACF:
R 11 ( s , t ) E { x 1 ( s ) x 2 ( t ) }
Substitute Equation (1) into Equation (A1) and we will arrive at:
R 11 ( s , t ) = E { k = d k u = 0 19 α u v = 0 2045 c v p T c ( s [ 40 , 920 k + 2046 u + v ] T c ) n = d n p = 0 19 α p q = 0 2045 c q p T c ( t [ 40 , 920 n + 2046 p + q ] T c ) } = k = n = u = 0 19 p = 0 19 v = 0 2045 q = 0 2045 E { c v c q } E { d k d n } α u p T c ( s [ 40 , 920 k + 2046 u + v ] T c ) α p p T c ( t [ 40 , 920 n + 2046 p + q ] T c ) = k = u = 0 19 p = 0 19 v = 0 2045 α u p T c ( s [ 40 , 920 k + 2046 u + v ] T c ) α p p T c ( t [ 40 , 920 k + 2046 p + v ] T c ) = p = 0 19 p T c ( t [ 40 , 920 k + 2046 p + v ] T c )
where k and v are selected as:
k = floor ( s / T b ) v = mod ( floor ( s / T c ) , 2046 )
such that:
0 s [ 40 , 920 k + 2046 u + v ] T c < T c
We can see from the last line of Equation (A2) that the effects of NH code are cancelled out by the 20 ms integration time. Next, we will derive ACF for the desired signal in a similar way.

Appendix A.2. Desired Signal ACF

Similar to the process revealed in Equation (A2) through Equation (A4), the receiver local replica (desired signal) ACF could be derived as:
R 22 ( s , t ) = E { k = u = 0 2045 c k p T c ( s [ 2046 k + u ] T c ) n = v = 0 2045 c n p T c ( t [ 2046 n + v ] T c ) } = k = p T c ( t [ 2046 k + v ] T c )
where:
v = mod ( floor ( s / T c ) , 2046 )
This ACF for the desired signal will be used in Appendix B where correlator output variance will be derived. This correlator output variance will be transferred to SSC and finally, the analytical form of equivalent white noise, i.e., model for the effects of self-interference will be developed.

Appendix B. Correlator Output Variance Derivations

We will develop analytical forms of correlator output variance in two situations: (i) data bits are aligned or (ii) data bits are misaligned by an integer number of primary code chips. The third situation where data bits are misaligned by an arbitrary fraction of a primary code chips is dealt with in Section 3.3Data Bits Misaligned by a Fraction of a Code Chip” of Section 3 “SSC and Equivalent White Noise Level”. For simplicity, we use notation a to denote function floor (a).

Appendix B.1. Data Bits Aligned

In this case, the variance at the output of correlator will be:
E { | y k | 2 } = k T b ( k + 1 ) T b k T b ( k + 1 ) T b p = p T c ( t 2046 k T c v T c ) e j 2 π f s e j 2 π f t d s d t = k T b ( k + 1 ) T b e j 2 π f s k T b ( k + 1 ) T b p = p T c ( t 2046 k T c v T c ) e j 2 π f t d t d s
the inner part of which is a Fourier transform:
F [ p T c ( t ) ] = T c sin c ( f T c ) e j π f T c
and therefore we have:
E { | y k | 2 } = T c sin c ( f T c ) e j π f T c k = 0 19 e j 2 π f k T k T b ( k + 1 ) T b e j 2 π f ( s v T c ) d t d s
where:
k = 0 19 e j 2 π f k T = 1 e j 2 π f T b 1 e j 2 π f T
where exponential sum formula is applied. Without loss of generality, we focus on the first data bit:
0 T b e j 2 π f ( s v T c ) d s = k = 0 19 k T ( k + 1 ) T e j 2 π f ( s ( s / T c 2046 s / T ) T c ) d s = k = 0 19 e j 2 π f k T m = 0 2045 m T c ( m + 1 ) T c e j 2 π f ( s s / T c T c ) d s = k = 0 19 e j 2 π f k T m = 0 2045 0 T c e j 2 π f s d s = 2046 T c sin c ( f T c ) e j π f T c k = 0 19 e j 2 π f k T
which assumes k = 0 in Equation (A9), and finally we have:
E { | y k | 2 } = T c [ sin c ( f T c ) ] 2 1 e j 2 π f T b 1 e j 2 π f T 1 e j 2 π f T b 1 e j 2 π f T = T T c sin 2 ( π f T c ) ( π f T c ) 2 sin 2 ( π f T b ) sin 2 ( π f T )
which is the variance at the output of correlator. Next, we will derive analytical form of correlator output variance when the data bits are misaligned by an integer number of code chips.

Appendix B.2. Data Bits Misaligned by an Integer Number of Code Chips

In this case, the variance at the output of correlators will be:
E { | y k | 2 } = s = 0 Δ t = 0 T b p = 0 19 p T c ( t Δ p T + T b v T c ) e j 2 π f ( s t ) d t d s + s = Δ T t = 0 T b p = 0 19 p T c ( t Δ p T v T c ) e j 2 π f ( s t ) d t d s
whose integration area will be divided into four smaller regions and be dealt with separately to facilitate derivation:
D 1 = { s :   0 s < Δ ,   mod ( s , T c ) < C T c } D 2 = { s :   0 s < Δ ,   mod ( s , T c ) C T c } D 3 = { s :   Δ s < T b ,   mod ( s , T c ) < C T c } D 4 = { s :   Δ s < T b ,   mod ( s , T c ) C T c }
This leads to a four-part integration:
E { | y k | 2 } = s = 0 s D 1 Δ t = 0 T b p = 0 19 p T c ( t Δ p T + T b v T c ) e j 2 π f ( s t ) d t d s + s = 0 s D 2 Δ t = 0 T b p = 0 19 p T c ( t Δ p T + T b v T c ) e j 2 π f ( s t ) d t d s + s = Δ s D 3 T t = 0 T b p = 0 19 p T c ( t Δ p T v T c ) e j 2 π f ( s t ) d t d s + s = Δ s D 4 T t = 0 T b p = 0 19 p T c ( t Δ p T v T c ) e j 2 π f ( s t ) d t d s
which becomes:
E { | y k | 2 } = s = 0 s D 1 Δ e j 2 π f s p = 19 K 19 e j 2 π f ( Δ + p T T b + v T c ) d s + s = 0 s D 2 Δ e j 2 π f s p = 20 K 19 e j 2 π f ( Δ + p T T b + v T c ) d s + s = Δ s D 3 T e j 2 π f s p = 0 18 K e j 2 π f ( Δ + p T + v T c ) d s + s = Δ s D 4 T e j 2 π f s p = 0 19 K e j 2 π f ( Δ + p T + v T c ) d s
The first integration of Equation (A16) is:
s = 0 s D 1 Δ e j 2 π f s p = 19 K 19 e j 2 π f ( Δ + p T T b + v T c ) d s = e j 2 π f T b p = 19 K 19 e j 2 π f p T s = 0 s D 1 Δ e j 2 π f ( s Δ v T c ) d s = e j 2 π f T b p = 19 K 19 e j 2 π f p T n = 0 K n T C T c + n T e j 2 π f ( s Δ v T c ) d s = e j 2 π f T b p = 19 K 19 e j 2 π f p T n = 0 K e j 2 π f ( n K 1 ) T m = 0 C 1 ( m C ) T c ( m C + 1 ) T c e j 2 π f ( α α / T c T c ) d α = C T c sin ( π f T c ) ( π f T c ) e j π f T c sin 2 ( π f ( K + 1 ) T ) sin 2 ( π f T )
The second integration of Equation (A16) is:
s = 0 s D 2 Δ e j 2 π f s p = 20 K 19 e j 2 π f ( Δ + p T T b + v T c ) d s = e j 2 π f T b p = 20 K 19 e j 2 π f p T s = 0 s D 2 Δ e j 2 π f ( s Δ v T c ) d s = e j 2 π f T b p = 20 K 19 e j 2 π f p T n = 0 K 1 C T c + n T ( n + 1 ) T e j 2 π f ( s Δ v T c ) d s = e j 2 π f T b p = 20 K 19 e j 2 π f p T n = 0 K 1 e j 2 π f ( n K ) T ( n K ) T ( n K + 1 ) T C T c e j 2 π f ( α α / T c T c ) d α = e j 2 π f T b p = 20 K 19 e j 2 π f p T n = 0 K 1 e j 2 π f ( n K 1 ) T ( 2046 C ) 0 T c e j 2 π f α d α = ( 2046 C ) T c sin ( π f T c ) ( π f T c ) e j π f T c sin 2 ( π f K T ) sin 2 ( π f T )
The third integration of Equation (A16) is:
s = Δ s D 3 T b e j 2 π f s p = 0 18 K e j 2 π f ( Δ + p T + v T c ) d s = p = 0 18 K e j 2 π f p T s = Δ s D 3 T b e j 2 π f ( s Δ v T c ) d s = p = 0 18 K e j 2 π f p T n = 0 18 K e j 2 π f ( n K ) T ( n K ) T ( n K + 1 ) T C T c e j 2 π f ( α α / T c T c ) d α = p = 0 18 K e j 2 π f p T n = 0 18 K e j 2 π f n T m = 0 C 1 m T c ( m + 1 ) T c e j 2 π f ( α m T c ) d α = C T c sin ( π f T c ) ( π f T c ) e j π f T c sin 2 ( π f ( 19 K ) T ) sin 2 ( π f T )
The fourth integration of Equation (A16) is:
s = Δ s D 4 T b e j 2 π f s p = 0 19 K e j 2 π f ( Δ + p T + v T c ) d s = ( 2046 C ) T c sin ( π f T c ) ( π f T c ) e j π f T c sin 2 ( π f ( 20 K ) T ) sin 2 ( π f T )
Finally, we combine Equation (A17) through Equation (A20) and will arrive at:
E { | y k | 2 } = T c 2 [ sin c ( f T c ) ] 2 ( ( 1023 C ) sin 2 ( π f K T ) + sin 2 ( π f ( 20 K ) T ) sin 2 ( π f T ) + C sin 2 ( π f ( K + 1 ) T ) + sin 2 ( π f ( 19 K ) T ) sin 2 ( π f T ) )
This completes the development of the variance of correlator output when the data bits are misaligned by an integer multiple of primary code chips. Specifically, when C = 0, i.e., the differential time delay between undesired and desired signals are integer multiples of a code period, we arrive at:
E { | y k | 2 } = T T c [ sin c ( f T c ) ] 2 sin 2 ( π f K T ) + sin 2 ( π f ( 20 K ) T ) sin 2 ( π f T )
which is a special case and we put it here for ease of reference.

Appendix C. Simulation Case For Section “Results”

The receiver is placed at 31°13.3′ N, 121°27.12′ E, Shanghai, China. The TLE retrieved for the simulation results is provided in Table A1.
Table A1. TLE for the simulation results in Table 1.
Table A1. TLE for the simulation results in Table 1.
SatelliteNo.Inclination (degrees)Right Ascension of the Ascending Node (degrees)EccentricityArgument of Perigee (degrees)Mean Anomaly (degrees)Mean Motion (revol. per day)
362871.622410.93090.0003751185.323528.15981.00265967
365901.508529.10930.0004961300.4656312.76271.00277041
3682854.3192201.27970.0039362203.8556127.66771.00249186
372100.928251.94350.0005461156.242831.27821.00271739
3725654.0451320.34270.0030637199.8493312.39011.00256291
3738456.880581.18490.0025638195.5632159.19831.00273265
3776354.6352203.51920.0035745204.2305243.69251.00288677
3794854.1427319.78620.0029394199.0165307.88741.00268580
380911.091426.44220.0000830112.2111164.48661.00268964
3825055.722380.52010.0023395212.374012.31011.86233381
3825155.658580.00010.0025721205.481017.64551.86234185
3877454.7564200.24060.0034012171.3479188.75711.86251530
3877554.8606199.76560.0019260217.0437222.96541.86233085
389530.137420.01790.0001222337.9805294.81311.00272848
For ease of reference, TLE elements used in simulation are extracted and titled.

References

  1. Morton, Y.; Miller, M.; Tsui, J.; Lin, D.; Zhou, Q. GPS Civil Signal Self-Interference Mitigation During Weak Signal Acquisition. IEEE Trans. Signal Process. 2007, 55, 5859–5863. [Google Scholar] [CrossRef]
  2. Gao, G.X.; Enge, P. How Many GNSS Satellites are Too Many? IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 2865–2874. [Google Scholar]
  3. Cerruti, A.P.; Rushanan, J.J.; Winters, D.W. Modeling C/A on C/A interference. In Proceedings of the 2009 International Technical Meeting of The Institute of Navigation, Anaheim, CA, USA, 26–28 January 2009; pp. 142–156. [Google Scholar]
  4. Cerruti, A.P.; Betz, J.W.; Rushanan, J.J. Further investigations into C/A-to-C/A interference. In Proceedings of the 2014 International Technical Meeting of The Institute of Navigation, San Diego, CA, USA, 27–29 January 2014; pp. 349–361. [Google Scholar]
  5. Liu, L.; Zhan, X.; Liu, W.; Niu, M. Assessment of Radio Frequency Compatibility between Compass Phase II and Other GNSSs. J. Navig. 2011, 64, 55–72. [Google Scholar]
  6. Liu, W.; Du, G.; Zhan, X.; Zhai, C. Assessment of Radio Frequency Compatibility Relevant to the Galileo E1/E6 and Compass B1/B3 Bands. J. Navig. 2010, 63, 419–434. [Google Scholar] [CrossRef]
  7. Hegarty, C.; Tran, M.; Lee, Y. Simplified techniques for analyzing the effects of non-white interference on GPS receivers. In Proceedings of the 15th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GPS 2002), Portland, OR, USA, 24–27 September 2002; pp. 620–629. [Google Scholar]
  8. Dierendonck, A.J.V.; Hegarty, C. Methodologies for assessing intrasystem and intersystem interference to satellite navigation systems. In Proceedings of the 15th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS 2002), Portland, OR, USA, 24–27 September 2002; pp. 1241–1250. [Google Scholar]
  9. Hegarty, C.; Dierendonck, A.J.V. Recommendations on digital distortion requirements for the civil GPS signals. Proceeding of IEEE/ION-PLANS-2008: Position, Location and Navigation Symposium, Monterey, CA, USA, 6–8 May 2008; pp. 1090–1099. [Google Scholar]
  10. Shibata, T.; Maeda, H. Extended theory of spectral separation coefficient for GNSS signal interference. In Proceedings of the 2010 International Technical Meeting of the Institute of Navigation, San Diego, CA, USA, 25–27 January 2010; pp. 930–940. [Google Scholar]
  11. Dierendonck, A.J.V.; Erlandson, R.J.; Shallberg, K.; Ericson, S. The inadequacy of the spectral separation coefficient and aggregate gain factor for quantifying the effects of GPS C/A code self interference. Proceedings pf the 26th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS+ 2013), Nashville, TN, USA, 16–20 September 2013; pp. 1435–1444. [Google Scholar]
  12. Hegarty, C. A simple model for C/A-code self-interference. In Proceedings of the 27th International Technical Meeting of The Satellite Division of the Institute of Navigation (ION GNSS+ 2014), Tampa, FL, USA, 8–12 September 2014; pp. 3484–3494. [Google Scholar]
  13. Golshan, R.; Fan, T.; Arrendondo, A.; Stansell, T. Implications of C/A-on-CA interference on carrier tracking loop performance. In Proceedings of the 27th International Technical Meeting of The Satellite Division of the Institute of Navigation (ION GNSS+ 2014), Tampa, FL, USA, 8–12 September 2014; pp. 3510–3525. [Google Scholar]
  14. Dierendonck, A.J.V.; Erlandson, R.; McGraw, G.; Coker, R. Determination of C/A code self-interference using cross-correlation simulations and receiver bench tests. In Proceedings of the 15th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GPS 2002), Portland, OR, USA, 24–27 September 2002; pp. 630–642. [Google Scholar]
  15. Bhuiyan, M.Z.H.; Söderholm, S.; Thombre, S.; Ruotsalainen, L.; Kuusniemi, H. Overcoming the Challenges of BeiDou Receiver Implementation. Sensors 2014, 14, 22082–22098. [Google Scholar] [CrossRef] [PubMed]
  16. Gardner, W.A. Introduction to Random Processes: With Applications to Signals and Systems, 2nd ed.; Mcgraw-Hill Inc.: New York, NY, USA, 1990. [Google Scholar]
  17. Betz, J.W.; Goldstein, D.B. Candidate designs for an additional civil signal in GPS spectral bands. In Proceedings of the 15th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GPS 2002), Portland, OR, USA, 24–27 September 2002; pp. 622–631. [Google Scholar]
  18. Betz, J.W.; Kolodziejski, K.R. Generalized theory of code tracking with an early-late discriminator part II: Noncoherent processing and numerical results. IEEE Trans. Aerosp. Electron. Syst. 2009, 45, 1557–1564. [Google Scholar] [CrossRef]
  19. Zhang, X.; Zhan, X.Q. Joint determination of precorrelation bandwidth, sampling frequency and quantization in wideband compass receivers. Navigation 2012, 59, 195–208. [Google Scholar] [CrossRef]
  20. Misra, P.; Enge, P. GPS Measurements and Error Sources. In Global Positioning System-Signals, Measurements, and Performance, 2nd ed.; Ganga-Jamuna Press: Lincoln, MA, USA, 2006. [Google Scholar]
  21. Hwa Create Co. HWA-RNSS-7400 Technical Specifications; Hwa Create Co., Ltd.: Beijing, China, 2015. [Google Scholar]
  22. Caspl Co. Multiconst Ellation Multifrequency Data Logging and Replay System (SAS6862A) Technical Specifications; CASPL Co., Ltd.: Shanghai, China, 2015. [Google Scholar]
  23. NORAD Two-Line Element Sets Current Data. North American Aerospace Defense Command. Available online: http://www.celestrak.com/NORAD/elements/ (accessed on 25 April 2015).
  24. Bhuiyan, M.Z.H.; Söderholm, S.; Kuusniemi, H.; Thombre, S.; Ruotsalainen, L. Utilization of a Novel Channel Quality Index for Improved Multi-GNSS Positioning in GNSS-denied Environments. In Proceedings of the 5th International Galileo Science Colloquium, Braunschweig, Germany, 27–29 October 2015. [Google Scholar]
  25. Muthuraman, K.; Borio, D. C/N0 estimation for modernized GNSS signals: Theoretical bounds and a novel iterative estimator. Navigation 2010, 57, 309–323. [Google Scholar] [CrossRef]
Figure 1. Generic correlator model.
Figure 1. Generic correlator model.
Sensors 17 00663 g001
Figure 2. Spreading code, navigation data bit, and NH code relationship for B1-I.
Figure 2. Spreading code, navigation data bit, and NH code relationship for B1-I.
Sensors 17 00663 g002
Figure 3. ACF for the undesired signal, which is a pulse train of 20 pulses, each of which has a width of Tc. Note that this figure is not to scale.
Figure 3. ACF for the undesired signal, which is a pulse train of 20 pulses, each of which has a width of Tc. Note that this figure is not to scale.
Sensors 17 00663 g003
Figure 4. Undesired and desired signals with navigation data bits perfectly aligned.
Figure 4. Undesired and desired signals with navigation data bits perfectly aligned.
Sensors 17 00663 g004
Figure 5. Undesired and desired signals misaligned by ∆.
Figure 5. Undesired and desired signals misaligned by ∆.
Sensors 17 00663 g005
Figure 6. One repetition of R11(s − ∆, t − ∆)R22(s,t) as the result of data bit misalignment, with 0 ≤ s < ∆.
Figure 6. One repetition of R11(s − ∆, t − ∆)R22(s,t) as the result of data bit misalignment, with 0 ≤ s < ∆.
Sensors 17 00663 g006
Figure 7. One repetition of R11(s − ∆, t − ∆)R22(s,t) as the result of data bit misalignment, with ∆ ≤ s < Tc.
Figure 7. One repetition of R11(s − ∆, t − ∆)R22(s,t) as the result of data bit misalignment, with ∆ ≤ s < Tc.
Sensors 17 00663 g007
Figure 8. Test setup for model validation.
Figure 8. Test setup for model validation.
Sensors 17 00663 g008
Figure 9. Measured vs. predicted tracking standard deviation.
Figure 9. Measured vs. predicted tracking standard deviation.
Sensors 17 00663 g009
Figure 10. SSCB1-I variation as a function of differential Doppler (between desired and undesired signals) and differential time delay. K is even.
Figure 10. SSCB1-I variation as a function of differential Doppler (between desired and undesired signals) and differential time delay. K is even.
Sensors 17 00663 g010
Figure 11. SSCB1-I variation as a function of differential Doppler (between desired and undesired signals) and differential time delay. K is odd.
Figure 11. SSCB1-I variation as a function of differential Doppler (between desired and undesired signals) and differential time delay. K is odd.
Sensors 17 00663 g011
Figure 12. Time evolution of equivalent noise level for PRN 6, 8, and 9 in current BDS constellation.
Figure 12. Time evolution of equivalent noise level for PRN 6, 8, and 9 in current BDS constellation.
Sensors 17 00663 g012
Figure 13. Code tracking standard deviation for different C/N0 values.
Figure 13. Code tracking standard deviation for different C/N0 values.
Sensors 17 00663 g013
Figure 14. Proposed Signal Quality Monitoring (SQM) based on self-interference evaluation when selecting satellites for navigation.
Figure 14. Proposed Signal Quality Monitoring (SQM) based on self-interference evaluation when selecting satellites for navigation.
Sensors 17 00663 g014
Table 1. Model utilization example.
Table 1. Model utilization example.
PRNRelative Received Power (dB)True Range (m)Doppler (Hz)SSCB1-I (dB/Hz)I0 (dBW/Hz)
6611.06736,868,693−563.659--
8811.25536,071,721−525.663−62.7−211.4
9910.45939,531,328578.098−76.6−226.2
Total−211.3
Table 2. Model utilization procedure.
Table 2. Model utilization procedure.
1:  while GPS measurement do
2:  Use true range r to obtain transit time t by t = r/c,
      where c is the speed of light in vacuum.
3:  Use relative received power p_rel to obtain actual
      received power p by p = p_rel + p_ref, where
      p_ref is the reference power set in the software simulator.
4:  Obtain differential Doppler f by differencing the first satellite with the rest.
5:  Obtain differential transit time ∆ by differencing the first satellite with the rest.
6:  Obtain the number of 1ms in ∆, K by K = floor (absolute value of ∆ × 103).
7:  Obtain the number of code chips in ∆, C by C = round( (∆ − K) × chip_rate × 10−3).
8:  if C is equal to code period (2046) then
9:     reset C.
10:     increment K.
11:  end if
12:  Use Equations (20) and (23) to obtain SSCB1-I.
13:I0 = p + SSCB1-I.
14:  end while
Table 3. Parameters used in model validation.
Table 3. Parameters used in model validation.
SectionsParametersValues
Satellite Transmitter (TX)TX bandwidth BT30.69 MHz
Target signalBDS B1
Atmospheric Effects, multipath & RF interferencenoneNot Applicable
RFFEPre-correlation Bandwidth Br4 MHz
Sampling frequency fs62 Msps
Number of quantization levels NQ4
Noise floor−203.5 dBW/Hz
Power level at antenna−160 dBW
Software ReceiverDLL one-sided noise bandwidth BL1 Hz
Coherent Integration Time (CIT) CT20 ms
Non-coherent Integration Time (NIT) NT1
Discriminator typeDot product
Early-late spacing ∆0.5 chips

Share and Cite

MDPI and ACS Style

Zhang, X.; Zhan, X.; Feng, S.; Ochieng, W. An Analytical Model for BDS B1 Spreading Code Self-Interference Evaluation Considering NH Code Effects. Sensors 2017, 17, 663. https://doi.org/10.3390/s17040663

AMA Style

Zhang X, Zhan X, Feng S, Ochieng W. An Analytical Model for BDS B1 Spreading Code Self-Interference Evaluation Considering NH Code Effects. Sensors. 2017; 17(4):663. https://doi.org/10.3390/s17040663

Chicago/Turabian Style

Zhang, Xin, Xingqun Zhan, Shaojun Feng, and Washington Ochieng. 2017. "An Analytical Model for BDS B1 Spreading Code Self-Interference Evaluation Considering NH Code Effects" Sensors 17, no. 4: 663. https://doi.org/10.3390/s17040663

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