Detection Performance of Polarization and Spatial Diversities for Indoor GNSS Applications

Multipath fading in the form of signal power fluctuation poses a formidable challenge to GNSS signal detection in harsh multipath environments such as indoors. Antenna diversity techniques such as polarization and spatial diversities can be used to combat multipath fading in wireless propagation channels. This paper studies and compares GPS signal detection performance enhancements arising from the spatial and polarization diversity techniques. Performance enhancements are quantified from a theoretical perspective and later verified based on several test measurements in various indoor environments. Enhancement is quantified based on measuring the correlation coefficient values between diversity branches, SNR levels, and computing the level crossing rate and average fade duration. In addition, the processing gain is quantified and the performance of each individual diversity system is evaluated. Experimental results show that, for a given target detection performance in terms of the probability of false alarm and the probability of detection, the required input SNR level to meet the target detection performance can be significantly reduced utilizing the diversity system.


Introduction
GNSS signal detection in indoor environments is a challenging problem since the signal suffers from insertion loss and fading.The magnitude of the insertion loss depends on the carrier frequency, the obstacle material, and the angle of incidence and can be as high as 30 dB [1,2].In addition to the insertion loss, the electromagnetic waves propagation in indoor environments also suffers from fading arising from the superposition of simultaneously arriving multipath signals with random phases and amplitudes.In reality, a moving handheld GNSS receiver experiences highly dynamic multipath situations resulting in high rate deep fading.In order to overcome these difficulties, various approaches such as High-Sensitivity GPS (HSGPS) [3][4][5][6], Assisted GPS (AGPS) [7,8], modernized GPS signals, and antenna diversity techniques [9,10] have been developed.
Antenna diversity systems are establish based on the fact that multiple antennas with different patterns, polarizations, and locations provide multiple diversity branches such that the moment when the multipath fading affects the output of a diversity branch, another branch may have a reasonable signal power.Therefore, combining the independent copies of transmitted signals received by the diversity branches leads to less fading and higher overall SNR which improves the detection procedure.In antenna diversity systems, the receiver uses multiple antennas with different characteristics to collect statistically independent signals to establish a diversity system.These antennas can be different in terms of polarization or location, resulting in diversity systems known as polarization and spatial diversities, respectively.The signal detection enhancement utilizing a pair of spatially separated antennas is studied in [11].It has been shown that a net processing gain can be achieved utilizing such a system.Colburn et al. [12] evaluate the spatial diversity performance for three different antenna configurations for the indoor 902-928 MHz propagation channel.They show that the indoor propagation channel can be modeled by either a Rician or a Rayleigh distribution.The time and phase diversity techniques are studied in [13] and [14], respectively.
International Journal of Antennas and Propagation Narayanan et al. [15] extensively analyse the polarization diversity and the channel characteristics at 1800 MHz which is close to the L1 frequency.They analysed the effect of the correlation coefficient between two diversity branches on diversity gain thoroughly.In addition, they evaluated the performance of the polarization diversity using different antenna configurations in both Rayleigh and Rician fading environments.In [16], the performance of a dual channel spatial diversity system at 1800 MHz has been considered.Lemieux et al. [17] have experimentally compared the space, frequency, and polarization diversities in the indoor environment for frequency of 900 MHz.Recently, au [18] proposed a new technique based on the motion of a single antenna to form a spatially distributed synthetic array for enhancing the detection performance of the GPS signals in the indoor environments.
The performance of dual-antenna spatial and polarization diversities to mitigate the multipath fading in the indoor environments is studied here.In the polarization diversity system, two structures have been considered.The first configuration consists of two circular polarized antennas, whereas the second one consists of two linear polarized antennas.A comparison is made based on the detection performances of each diversity scheme.Consequently, a detection scheme based on a Rayleigh fading channel is demonstrated and the processing gain realizable through utilizing spatial and polarization diversity techniques is quantified theoretically and experimentally.The second-order statistics of the multipath channel, namely, the Level Crossing Rate (LCR) and Average Fade Duration (AFD) metrics, are experimentally quantified and analyzed.
The remainder of this paper is organized as follows.The GPS signal structure and system model is described in Section 2. GPS signal polarization variations in multipath fading environments are also considered in this section.In Section 3, the signal detection procedure and the diversity approach utilized in this work is discussed and analyzed theoretically.The experiments performed are described in Section 4 and the empirical results are examined as well.Finally, Section 5 concludes the paper.

Signal and Channel Model
2.1.Signal Model.Herein, the demodulated complex baseband antenna output is denoted by r(t).The signal component of r(t) emanating from the kth satellite is denoted as s k (t).This signal can be expressed as s k (t) = A k s ok (t) where s ok (t) is the deterministic complex baseband component of the kth satellite signal that is known to the receiver and A k is the channel response to the incident signal.The received signal is corrupted with additive noise which has an equivalent complex baseband representation denoted by w(t).It is assumed that w(t) is a complex normal random process, independent of the signal with a power spectral density (PSD) that is constant within the bandwidth of the received signal with a level of N o .Therefore, r(t) can then be expressed as Reflector where A k (p(t)) is the complex channel gain which is a function of the antenna position p(t), D k (t) is the navigation data modulation, c k (t) is the PRN code, τ k is the code phase, Δ f k represents the carrier frequency offset, ψ k is the initial phase offset, K is the number of satellites in view, and w(t) is complex additive white Gaussian noise.Herein, the signal detection of individual GNSS satellites will be considered, while the presence of simultaneous GNSS signals is not considered.This is due to the fact that the cumulative selfinterference from other PRNs is negligible compared to other dominant independent noise sources [19].Hence, in the remainder of the text, for the notational convenience, the subscript of k is ignored and the received signal is modeled as where s o (t) = D(t − τ)c(t − τ)e j(2πΔ f t+ψ) is known to the receiver except for the navigation data D(t), the code phase τ, the carrier frequency offset Δ f , and the initial phase offset ψ.In unresolvable Rayleigh-faded multipath environments which is assumed here, A(p) can be modeled by a Complex Normal (CN) random variable, such that A ∼ CN(0, σ 2 A ) where ∼ denotes the PDF of the left-hand side variable and CN(μ, σ 2 ) signifies a complex normal probability density function (PDF) of mean μ and variance σ 2 [20].
The received signal passes through the signal acquisition procedure as an initial operation dealing with the detection and coarse estimation of Doppler and code phase offset.A two-dimensional search is applied on the signal to roughly estimate the Doppler frequency and the code phase offset.

GPS Signal Polarization in Dense Multipath Environments.
The transmitted GPS signal is a right hand circularly polarized (RHCP) electromagnetic (EM) wave which is a special case of the more general elliptically polarized wave.Assume that an RHCP GPS signal reaches the reflector with an incident angle (grazing angle) of θ g as shown in Figure 1.For a circularly polarized signal, the electric field intensity vector E is decomposed into two orthogonal components, E h and E v , called horizontal and vertical components, respectively.Since the horizontal and vertical components of the reflection coefficient of a reflector are not equal [20], the absolute value of the E field components in the reflected signal changes independently.Therefore, the circular polarity of the reflected GPS signal changes to elliptical.In addition to magnitude, the phase of each component changes depending on the Brewster angle, Ψ B (known as the polarization angle as well) as a property of the reflectors [21].For instance, Brewster angle for a metallic reflector is a few degrees [21] while for a concrete surface is around 69 degrees [22].After reflection, the phase of the horizontal component (E h ) always changes by 180 degrees; however, that of the vertical component (E v ) may change either by 0 or 180 degrees depending on whether the grazing angle is below or above the Brewster angle [21].If the grazing angle is higher than the Brewster angle, the vertical component phase changes by 180 degrees and the reflected signal has right hand polarization.On the other hand, if the grazing angle is lower than the Brewster angle, the reflected signal has left-hand polarization.Hence, the reflected wave may have either left-or right-hand elliptical polarization based on the impinging angle.Balanis [23] shows that an elliptical polarized signal can be written as a combination of two RHCP and LHCP waves.Hence, if the grazing angle of the signal is less than the Brewster angle, the RHCP component of the reflected signal is stronger whereas, if the grazing angle of the LOS signal is greater than the Brewster angle, the LHCP component of the reflected signal will be stronger.In indoor environments, since the received signals are highly likely to be subject to reflections, the final received signals can be assumed to be elliptically polarized and, hence, include both RHCP and LHCP components.The arriving RHCP and LHCP waves arise from the superposition of many reflected waves.Since the reflection phenomenon is random, the final RHCP and LHCP waves are uncorrelated and, therefore, can be used for diversity gain.This diversity gain is quantified through theoretical analysis and verified based on actual measurements in this paper.

Diversity Scheme and Detection Procedure
Having collected GPS signals with the diversity branches, the received signals should be combined using an efficient combiner.Herein, the main focus is to combine signals received from two diversity branches.In order to implement a detector, the likelihood ratio test (LRT) function is utilized.The detection problem is to find a moment in each code and Doppler search where the incoming signal and the locally generated one are synchronized.Two hypotheses H 0 and H 1 are defined as Here, m is the number of diversity branches, x[m] is the correlator output, s[m] is the desired signal in which the code phase and Doppler of the incoming signal and the locally generated one are synchronized, ω[m] is the embedded zero mean white complex Gaussian noise with a covariance matrix of C w = σ 2 ω I 2 where σ 2 ω is the noise variance, and I 2 is 2 × 2 identity matrix.In a Rayleigh-faded channel, the desired signal received at the antenna is a zero mean circular normal distribution with a covariance matrix of C s .Therefore, the signal probability distribution function under H 0 and H 1 hypotheses can be defined as Under these conditions, the final test statistic becomes [24] T In a general case, the received signal from different diversity branches may be correlated.In order to simplify the final test statistic represented in (6), the modal matrix V, a 2 × 2 matrix formed with the eigenvectors of C s , is used to decorrelate the input signals by where Λ s is a matrix of eigenvalues of C s and is formed as Therefore, the covariance matrix of y becomes After some manipulations, the final test statistic is reduced to [24] T y Consequently, the final test statistic can be considered as a decorrelator followed by a weighting function as shown in Figure 2.This detector is known as the Estimator-Correlator (EC).

Detection Performance Evaluation.
The detector performance is evaluated by the probability of false alarm P FA and the probability of detection P D .The probabilities of false alarm and detection for the final test statistic in (9) are quantified as  the P D and P FA .The probability function of y = V T x can be written as Thus, y has a zero mean Gaussian distributions and |y| 2 results in central chi-squared distributions.Hence, the characteristics function of T(y) leads to the summation of two weighted central chi-squared distributions.The probabilities of detection and false alarm for this model are quantified as [25] where To determine a closed-form expression for P FA and P D in the case of distinctive eigenvalues, it can be shown, using a partial fraction expansion, that (13) can be reduced to [24] Similarly, the P FA can be computed as where 3.2.Theoretical Analysis of the Indoor Environment.In the dual-branch diversity system, a general scenario for the received signal covariance matrix can be defined as where x 1 and x 2 are the correlator outputs of the input branches, σ 2 sm = E(s m s * m ) is the variance of s m under H 1 condition for m = {1, 2}, (•) * is a conjugate operator, and . To simplify the general equation for the signal covariance matrix, without loss of generality, assume that where r = σ 2 s2 /σ 2 s1 is the power ratio of the input signals.In addition, the cross-correlation coefficient (ρ) between the received signals is quantified as Therefore, the signal covariance matrix can be rewritten as As shown in (9), the detection test statistics is a function of the eigenvalues of the signal covariance matrix λ s1 and λ s2 .In this case, the eigenvalues of the covariance matrix are [26] The probability of detection for a fixed probability of false alarm (P FA = 0.01) as a function of the signal cross correlation coefficient and the input power difference is shown in Figure 3. Lower performance results from a higher power difference and cross-correlation coefficient between input branches.The interesting observation is that for the correlation coefficient values less than 0.4 the P D for given P FA does not significantly degrade.
In the case of uncorrelated equal power input signals (ρ = 0 and r = 1), the EC detector is simplified to equal gain (EG) combiner as However, combining the correlated unequal power signals by the equal gain (EG) combiner results in the probability of false alarm and probability of detection given by respectively, where λ x1 and λ x2 are the eigenvalues of C x = C s + σ 2 w I 2 and can be shown to be These will be used in the following section to compare the diversity gain performance of EC and EG combiners.

Diversity Gain Evaluation.
The effectiveness of a diversity system is evaluated by a quantity known as the diversity gain.Herein, the diversity gain is defined as the excess in the required input average SNR for a single antenna scheme to achieve the same P D as the combining scheme for a specific P FA .For a given target detection performance in terms of P D and P FA , the diversity gain is defined as where SNR S and SNR D signify the required SNR values for the single branch and the diversity combiner, respectively, to achieve the same detection performance.Figure 4 shows the diversity gain as a function of the cross-correlation coefficient and power difference between the diversity branches for a specific design point of P FA = 0.01 and P D = 0.9 using the estimator correlator (EC) as the detector.The results of Figure 4 show that up to 5.5 dB diversity gain for uncorrelated and equal power input signals can be achieved.As shown by increasing the correlation  coefficient between branches, the diversity gain declines.However, by increasing the power difference between input signals, the diversity gain reduces to 0 dB even for uncorrected input signals, which is due to the signal degradation in one of the branches.
In Figure 5, the diversity gain difference between the estimator correlator (EC) and equal gain (EG) combiner is shown as a function of power difference and crosscorrelation coefficient between the diversity branches.Here, the design point is defined as P FA = 0.01 and P D = 0.9. Figure 5 shows that EC and EG result in the same detection performance for low correlation coefficient between the input branches (up to 0.4).

Experimental Results
In the previous sections, the detection performances of the diversity systems in multipath fading environments modeled by the Rayleigh fading were theoretically examined.The experimental measurements described in this section validate the assumptions in the context of the GNSS signal detection in indoor environments.The objective of the experimental measurements conducted herein is to examine

Circular polarized antennas
Linear polarized antennas the dual-antenna diversity systems, namely, the polarization and spatial diversities for a selection of typical indoor locations.

Test Setup.
A set of experiments was performed in order to compare the performance of the polarization and spatial diversity structures.Measurement results in [18,27] have shown that dense indoor GNSS channels decorrelate by almost half a carrier wavelength.Hence, herein spatial diversity was achieved by two spatially RHCP antennas separated by 1.5 wavelength of the GPS L1 signal (30 cm).The spatial diversity main principle is based on the fact that in dense multipath environments the received signals are a random function of antenna location.Herein, a dual circular polarized antenna consisting of RHCP and LHCP components was used to implement circular polarization diversity.More details regarding the utilized antennas can be found in [28].Besides, a dual-polarized antenna consisting of two linear polarized antennas mounted on a horizontal plane and mutually perpendicular, called vertical and horizontal antenna, was utilized to create a linear polarization diversity structure.The principle of this diversity structure is based on the fact that the amplitude and the phase of the vertical and horizontal components of the reflected signals vary independently and can be considered as a diversity structure.
In order to achieve the objectives of this work, two datasets in two different indoor locations were collected.Each dataset consisted of three consecutive data collections with different diversity schemes, namely, the circular polarization (CP), spatial polarization (SP), and linear polarization (LP) structures, as shown in Figure 6.In each case, the diversity branches along with a reference antenna were connected to a synchronized triple port downconverter/digitizer to collect synchronous raw IF samples in the GPS L1 band.The indoor antennas were mounted on a linear motion table to collect independent signal samples in the indoor GPS channels.The table moved by about 2.8 m in each lap with an antenna speed of 2 cm/s.The first data collection set was performed in a laboratory with a high metallic ceiling as shown in Figure 7(a).The second dataset was carried out in the corridor of a building having large window panes in both the south and north directions and a concrete ceiling as shown in Figure 7(b).To aid the indoor receiver process with the navigation data, a reference antenna was located in an open-sky condition within 30 metres of the indoor diversity structures.The reference channel was used to remove the navigation data bit from the received signals to increase the coherent integration time up to 100 ms.The GSNRx-rr software GNSS receiver [29] was used to process the data.

Comparison of the Polarization and Spatial Diversities.
In this section, different metrics have been used to characterize different antenna diversity systems.

Correlation Coefficient.
As shown earlier, the correlation coefficient between diversity branches can be used to characterize the performance of a diversity system.Narayanan et al. in [15] quantify both complex and envelope correlation coefficient and conclude that the value of complex correlation coefficients is constant from run to run and hence it is a reliable metric.Herein, the complex correlation coefficient between the received signals in each diversity branch is quantified as [12] where E[•] is the expected value operator.The results of correlation coefficient measurements for circular polarization (CP), spatial and linear polarization (LP) diversity systems, and different PRNs are shown in Figure 8.According to the results, the correlation coefficient for the majority of the datasets is less than 0.3, which results in a diversity gain.
Considering the results of Figure 5, the estimator-correlator and equal-gain combiner performance does not have any significant difference for the correlation coefficient values less than 0.4.

Average Input SNR.
As mentioned before, the average input SNR plays a crucial role in the efficiency of a diversity system.The measured carrier-to-noise ratio (C/N 0 ) values of the received signals through different diversity branches are quantified and shown in Figure 9. GSNRx-rr was used to measure the C/N 0 values.Considering the results of Figure 9, the level difference of the received signals in the different diversity branches is not considerable and, given the theoretical results (Figure 5), is negligible.According to these findings and considering the results of Figure 5, the equal gain combiner is considered as the combining method herein.The interesting observation is that in all cases the received signal power through the LHCP antenna is slightly more than that of the RHCP antenna.As explained before, this implies that the LHCP component of the reflected signals is stronger than the RHCP component.Thus, considering the LHCP component as an input for the circular polarization, diversity leads to better signal detection performance.

Average Fade Duration and Level Crossing Rate.
In indoor GPS multipath channels modeled by Rayleigh fading, the received signal amplitude is subject to rapid fluctuations.
As the GPS receiver moves, the rate of change of amplitude 2 changes.The Rayleigh fading distribution just describes the first-order statistics of the signal amplitude for short distances where the mean level of the received signal can be assumed constant.The first-order statistics does not provide any information regarding the fading time or distance.It just provides an overall percentage of time or location where the signal amplitude lies below a specific threshold.Hence, it is of interest to characterize the rate at which fades of any depth happen and the average fade duration below any given depth.
In indoor GPS multipath channels, the fading envelope varies and the fading rate and envelope amplitude are a function of time.In order to obtain a quantitative measure of the fading occurring in the channel, the level crossing rate (LCR) and the average fade duration (AFD) are quantified as two statistical parameters of the channel model [30].LCR is quantified as the expected rate at which the received signal envelope crosses a specific level, in either the positive or negative direction.The AFD is the average fade duration below any given depth.In other words, AFD is quantified as the average period of time for which the received signal envelope is below a specific level.
2  The threshold for the LCR and AFD quantification procedure can be defined based on the level (signal envelope) or signal SNR.Since in GPS applications the SNR values has been extensively used, in this work the signal post-SNR is utilized for the LCR and AFD quantification.In statistical detection theory, the post-SNR can be quantified by the deflection coefficient as [24] where E( signals in all diversity structures show less fading compared to the single antenna cases.Figure 11 shows the measured AFD values and the significant improvement in the combined signals in all diversity structures.According to these results, overall all diversity structures result in the same improvement in the fading mitigation.

Diversity Gain.
The diversity gain for a specific design point in terms of P FA and P D is quantified as the reduction in the input average SNR of the diversity branches to provide the same detection performance utilizing a single channel receiver.The measured diversity gains for different diversity systems are shown in Figure 12 where the diversity gain is measured for P FA = 0.01.The average diversity gain of the collected data for circular polarization, spatial, and linear polarization diversity structures is 3.7 dB, 2.8 dB, and 2.2 dB, respectively.As discussed previously, since in indoor environments the GNSS signal is reshaped to the elliptically polarized waves which can be observed by the RHCP and LHCP antennas, the polarization diversity employing the RHCP and LHCP antennas improves the GPS signal detectability and leads to higher performance.However, spatial diversity systems take only the power of the RHCP component of the received signals, which leads to lower processing gain.

Conclusions
Spatial and polarization diversity systems were employed to combat fading in dense multipath environments.Performance enhancements arising from utilizing the spatial and polarization diversity techniques were first studied theoretically.It was shown that the detection performance of the EG combiner is similar to that of the EC for low correlation coefficients values (less than 0.4) and low power differences between diversity branches.Experiments were performed to compare polarization and spatial diversity structures.It was observed that both spatial and polarization diversity result in sufficiently low correlation coefficients between the received signals with a low power difference between input branches.However, the received signal power through the LHCP antenna outweighs the signal power received by the RHCP antenna.According to the results for the target detection design point in term of P FA and P D , the polarization diversity leads to higher diversity gain in indoor environments compared to the spatial diversity schemes.The diversity gain can reach up to 5 dB for the practical range of detection design point.LCR and AFD values, representing the secondorder statistics of the fading channel, were also measured for different diversity branches.

Figure 2 :
Figure 2: LRT detector scheme for a two-antenna diversity system.

Figure 3 :
Figure 3: P D for P FA = 0.01 for the EC combiner versus correlation coefficient and power difference between diversity branches.

6 CFigure 4 : 1 C
Figure 4: Diversity gain versus correlation coefficient and power difference between diversity branches.

Figure 5 :
Figure 5: EC and EG diversity gain difference versus correlation coefficient and power difference between diversity branches.

Figure 6 :
Figure 6: Three antenna diversity structures mounted on a moving linear table.

Figure 7 :
Figure 7: Data collection sites in indoor locations (a) Energy High Bay (b) corridor.

Figure 8 :
Figure 8: Correlation coefficients between two diversity branches for different diversity structures.

Figure 9 :
Figure 9: Average input SNR for different diversity structures.

Figure 12 :
Figure 12: Diversity gain for various diversity schemes.