Dual-Polarized Synthetic Antenna Array for GNSS Handheld Applications

Small portable Global Navigation Satellite System (GNSS) receivers have revolutionized personal navigation through providing real-time location information for mobile users. Nonetheless, signal fading due to multipath remains a formidable limitation and compromises the performance of GNSS receivers. Antenna diversity techniques, including spatial and polarization diversity, can be used to mitigate multipath fading; however, the relatively large size of the spatially distributed antenna system required is incompatible with the small physical size constraints of a GNSS handheld receiver. User mobility inevitably results in motion of the handset that can be exploited to achieve diversity gain through forming a spatially distributed synthetic array. Traditionally, such motion has been construed as detrimental as it decorrelates the received signal undermining the coherent integration processing gain generally necessary for acquiring weak faded GNSS signals. In this paper the processing gain enhancement resulting from a dual-polarized synthetic array antenna, compatible with size constraints of a small handset that takes advantage of any user imposed motion, is explored. eoretical analysis and experimental veri�cations attest the effectiveness of the proposed dual-polarized synthetic array technique by demonstrating an improvement in the processing gain of the GNSS signal acquisition operation.


Introduction
e initial GNSS signal acquisition is a necessary step in the process of obtaining a set of pseudorange estimates sufficient for location estimation.e underlying multihypothesis detection problem is typically exasperated by a large search space which limits the processing resources and time that can be appropriated for testing individual hypotheses [1,2].Multipath fading, in the form of spatial signal power �uctuations [3], poses a formidable challenge to GNSS signal acquisition, a problem which has recently attracted signi�cant interest [4,5].Utilizing a higher processing gain through a longer integration time is not an attractive solution as the mean acquisition time is signi�cantly increased [6].In addition, signal decorrelation due to oscillator instability and user mobility effectively limits the coherent integration time [7].
Recently multiple receive antennas have been utilized to enhance signal detection performance in the form of beamforming and diversity systems [8,9].Antenna diversity based on employing multiple antennas with different radiation characteristics has long been in use in diversity systems [10,11].It can be shown that for typical indoor multipath environments resembling a sphere of scatterers model, an antenna separation as small as half a carrier wavelength results in a negligible signal complex correlation coefficient [12,13].Even with this small element spacing, the size of a multielement antenna array is not compatible with the small physical size constraints of the GNSS handheld terminal.
Diversity achievable through approximately collocated antennas of differing polarization responses is an attractive option for the size constrained handset receiver [11,14,15].Most of the work on polarization diversity is based on orthogonally polarized antennas [11,16].In a polarization diversity system, colocated antennas with orthogonal polarizations are considered as uncorrelated diversity branches, provided that mutual coupling is minimal [11,14,16].A twobranch antenna diversity system based on combining outputs of a Right-Hand Circularly Polarized (RHCP) antenna and a Le-Hand Circularly Polarized (LHCP) antenna can be considered for GNSS applications [17].Although polarization diversity based on the circularly polarized GNSS signal provides two orthogonal polarizations, a higher number of diversity branches are typically desirable for further enhancements of the processing gain [5,18].Spatial diversity can provide these additional diversity branches.
As stated earlier, utilizing spatial diversity through an array of spatially distributed antennas is not compatible with the size factor of GNSS handheld terminals.However, an innate property of the mobile receiver is that it is generally undergoing some form of motion as part of the typical usage mode.While such motion is typically regarded as being detrimental as it results in signal decorrelation limiting the processing gain achievable by coherent signal integration, it can be exploited to achieve usable diversity gain.e resulting diversity gain can enhance the performance of GNSS signal acquisition in multipath faded environments.Authors in [5,18] have proposed a technique that utilizes spatial diversity based on motions of a single antenna and forming a spatially distributed synthetic array.Although the concept of synthetic array for signal parameter estimation has been around for several years (see [19]), its application for combating multipath fading was �rst introduced by [18].
Although a unipolarized synthetic array (the previous work, [5], which is herein called UPSA) enhances GNSS signal acquisition, capturing the required number of signal samples in spatially correlated channels requires longer trajectories of several carrier wavelengths.Nondeliberate motion on the part of the user may not naturally provide a trajectory of such length.To utilize shorter trajectories it is necessary to incorporate another form of diversity system in addition to the spatial diversity available.is paper proposes a combination technique based on merging polarization and spatial diversities through utilizing a dual-polarized RHCP-LHCP antenna in a GNSS handheld receiver and exploiting user movements to form a dual-polarized synthetic array, herein called DPSA (see Figure 1).
Note that the use of a dual-polarized antenna in a handheld has previously been reported in the literature (see [11]); nevertheless, the detection performance of a moving dual-polarized antenna in the context of GNSS signal acquisition has not been quanti�ed, before.is paper studies the effect of motion on the detection performance of a GNSS receiver based on a dual-polarized antenna and proposes a practical scheme for producing useful diversity gain.Here, the detection performance of the DPSA technique is compared to the previously proposed con�guration, UPSA, and it is shown that a DPSA con�guration devours a much shorter trajectory compared to UPSA, for given performance requirements.Also, through extensive theoretical analysis and experiments, it is shown that under similar conditions (trajectory length, wireless channel condition, etc.), DPSA outperforms UPSA and signi�cantly improves the initial acquisition performance.e remainder of the paper is organized as follows.
Section 2 de�nes the system model and the assumptions.e dual-polarized synthetic array model and the despreading process are also presented in this section.Section 3 discusses the realizable detection enhancement through comparing the performance of the proposed DPSA with that of the UPSA in a generalized correlated Rayleigh fading channel.Section 4 discusses the measurement results.Conclusions are given in Section 5.

Assumptions and System Model
e received complex GNSS baseband signal is denoted here by .e signal component of  emanating from the th satellite is represented by     which is a function of time, , as well as antenna position . denotes the position vector of the phase center of the receiver antenna.    can be decomposed to            where     is the complex baseband component of the kth satellite signal and is known to the receiver.   is the channel response to the incident signal at the antenna position of .e received signal, , is corrupted with additive noise which has an equivalent complex baseband representation denoted by .It is assumed that  is a complex normal random process, independent of the signal, and has a Power Spectral Density (PSD) that is constant within the bandwidth of the received signal.Herein, the signal detection of individual GNSS signals will be considered, while ignoring the presence of other GNSS signals.is is reasonable since the GNSS signal level is well below the noise �oor, thus cumulative interference from other GNSS satellites is negligible.Hence, in the remainder of this paper, for the notational convenience, the superscript "" associated with individual satellites is omitted and the received signal is given by          + where is the navigation data modulation,  is the Pseudorandom Noise (PRN) code,  is the code phase, Δ represents carrier frequency offset (due to the Doppler of the GNSS signal as well as any frequency offset of the receivers local oscillator), and  is the initial phase offset.   is known to the receiver except for the navigation data, the code phase, the carrier frequency offset, and the initial phase offset .
For indoor GNSS applications it can be shown that    is sufficiently narrowband such that a �at fading model can be assumed [3,5].It is also assumed that during the data collection, hereaer denoted by the signal snapshot period, , the channel is temporally stationary, and therefore coherent integration results in positive processing gain.Consequently, during the snapshot interval, , the channel gain, , is assumed to be a function of position only, not time.
where   and   are the signal and the noise components of the th collected sample,   .To simplify the expressions to follow and without any loss of generality, the noise covariance is assumed to be normalized such that   ∼ CN 1, where CN  2  is the circular normal distribution with   2 representing the complex mean and the variance, respectively.e variance of the signal component,   , aer Δ seconds of coherent integration is  2  .Note that  2  is the variance of the complex signal component aer  seconds of coherent integration.Consequently, a signal-to-noise ratio (SNR) can be de�ned as Note that noise is assumed to be solely intrinsic and independent of the signal.Consequently,   ,   , and   can be represented as a vector by where " †" denotes a matrix transpose, the superscripts "" stands for the principle polarization, and    denotes the th signal sample collected based on an RHCP antenna.
Similarly, the baseband signal samples collected based on a DPSA con�guration can be represented by where "" stands for the orthogonal polarizations and    denotes the th signal sample collected based on an LHCP antenna.
It is further assumed that the baseband output signals of the orthogonally polarized antennas are uncorrelated and identically distributed [20].In other words      *     and var     var    for all ,   1 …  .Here,  represents an expectation operator and " * " is a complex conjugate.

Multipath Model.
Multipath fading at the output of the receiver antenna in severe and isotropic multipath environments such as indoors can be modeled according to the "sphere of scatterers model" [20][21][22][23][24].According to this model, the magnitude of the complex correlation coefficient, |  |, between complex baseband voltages at the output of two copolarized isotropic antennas with spatial separation of "Δ" follows from [23,24] as where  is the carrier wavelength.Hence, according to (8), the received signal decorrelates by a half a carrier wavelength antenna displacement.erefore, by moving an antenna along a "  2" long smooth trajectory in a multipath fading environment, one can collect up to  spatially uncorrelated signal samples [5,18,20].Consider that these M signal samples resemble  diversity branches experiencing �at Rayleigh fading.(ere is some loss in generality by

ISRN Communications and Networking
ignoring the resolvable multipath and assuming �at fading.Nevertheless, the frequency selective fading results in an alternate source of diversity which can be considered in processing gain evaluations [25]).

Detection Performance and Processing Gain Analysis
is paper focuses on GNSS signal acquisition in a multipath faded environment.is generally implies a multihypothesis detection process as part of the general acquisition over a given search space.Generally as the hypotheses are tested with orthogonal despreading functions, the performance can be assessed by considering an individual hypothesis testing.
In this paper   ( is de�ned to correspond to the despreading function for a speci�c individual hypothesis test such that the detection process is binary [1,2].erefore ℋ1 represents the case when the incoming signal corresponds to the hypothesis associated with     , and ℋ0 represents the case when the incoming signal corresponds to the hypothesis associated with   .e approximation in the analysis is that it is assumed that for the ℋ1 state,   ( corresponds exactly to the spreading signal of the incoming signal except for an arbitrary complex amplitude coefficient.Hence, under ℋ1 it is assumed that the frequency and code phase of the incident signal is synchronized with   (.A more detailed analysis would include the incremental loss associated with a slight mismatch of the spreading signal of the incoming signal and   (.However, herein the objective is to compare the signal acquisition performances of the DPSA and UPSA techniques.Hence, this incremental additional loss which affects all cases in approximately the same way will not be considered.e design of the detection algorithm is based on target values of the probability of detection (  ) and the probability of false alarm (  ) associated with the search hypothesis.Assuming that the statistics of the signal and noise components are known, for a given spatial extent of the synthetic array, the target detection performance requirements of {  ,   } map into a speci�c required average signal-to-noise ratio (SNR).Regardless of further algorithm details, for a given utilization detection scenario in terms of   and   , the goal is to optimize the receiver processing such that the required average SNR is minimized.Consequently, the problem is modeled as one of choosing between ℋ0, the noise-only hypothesis, and ℋ1, the signal-present hypothesis as where, , , and  were de�ned earlier (see ( 5) and ( 7)).Under generalized correlated Rayleigh fading conditions, the distribution of  is according to a zero-mean jointly circular normal PDF with a covariance that is denoted by   .In other words,   CN(,   .e signal vector , can be mapped into a single detection variable , through a diversity combiner.Consequently, a decision is made based on comparing  to a given threshold, , which is typically selected based on a constraint of   [26]. As the carrier phase of the received signal is unknown, a square law combiner (SLC) is utilized to map  to  as where   0 corresponds to UPSA and   1 to DPSA scenarios [21].Accordingly, the measured   and   are related to the cumulative distribution function(CDF) of  and to the threshold, , based on where ℱ denotes the CDF of  under ℋ0 or ℋ1 [26].
Consequently, the corresponding values of average SNRs determined for a target {  ,   } will be compared in order to provide a measure of performance enhancement as where G is the realized gain and  DPSA ,  UPSA , and  Static are the required SNRs for achieving {  ,   } based on a DPSA, UPSA, or static antenna scenarios, respectively.

Performance Evaluation in Independent Rayleigh Fading
Channel.In this subsection, the performances of the UPSA, the DPSA, and those of a static antenna will be compared.
Observing the statistics of SLC under uncorrelated branch variables is insightful due to the simplicity of resulting mathematical expressions.Speci�cally, the advantage of the DPSA model over the UPSA model can be readily realized by noting the achievable extra degrees of freedom at the SLC output chi-squared PDF which will be discussed shortly.Assume a linearly polarized antenna is translated along a smooth /2-long trajectory during the signal snapshot interval of duration  such that  spatially uncorrelated samples can be collected (UPSA scenario).Consequently, it can be shown that ℋ0  CN(,    under the noise-only hypothesis and ℋ1  CN(, ( UPSA /      under the signal-present hypothesis.Here,   denotes an  by  identity matrix.ese samples can be combined using the SLC combiner of (10).Consequently the statistics of the detection variable  is where  2  ( 2 ,  2  denotes a chi-squared distribution with  degrees of freedom ( DOF) with noncentrality parameter  2 and the common variance of the corresponding Gaussian components  2 .
In a DPSA scenario (Figure 1), a doublet of orthogonally polarized antennas is translated over the same trajectory such that 2 orthogonal sets of  spatially uncorrelated samples can be measured.Consequently, it can be shown that   CN(  2 ) under the noise-only hypothesis and   CN( ( DPSA /) 2 +  2 ) under the signalpresent hypothesis.erefore, the SLC output statistics, , for DPSA is e extra DOF realized for the DPSA setup is a simple evidence of performance improvement.Figure 2 demonstrates the realizable processing gain (12) of a DPSA in comparison with a static antenna for various values of   and .As shown in this �gure, for any target detection requirement of {     }, there is a speci�c  that results in a minimum required SNR and therefore a highest processing gain, .erefore, there is an optimum  for which  is maximum which here is denoted by  opt .e existence of an optimum  is due to the increase of diversity gain as well as coherence losses as M increases [5].As shown there is about 20 dB gain utilizing the DPSA algorithm for target detection   = .99 and   = .. opt for this scenario is seven, which maps into a reasonable trajectory size for the GNSS handheld applications with the trajectory of six half carrier wavelengths corresponding to a total trajectory length of 54-75 cm for typical GNSS frequencies.
Figure 3 demonstrates the optimum  for different values of   and   .As shown, the required number of spatial samples,  opt , is lower for DPSA as compared with UPSA.In other words, DPSA devours a shorter trajectory than UPSA does for a given {     }. Figure 4 demonstrates the realizable processing gain of DPSA as compared with UPSA computed from (12).As can be seen from this �gure, for lower  and higher detection requirements, DPSA's performance is signi�cantly higher than UPSA's.By increasing  the performance of the DPSA converges to that of the UPSA.is is due to the fact that both techniques are taking advantage of diversity gain regardless of the source of the diversity (polarization or spatial diversity).

Performance Evaluation in Correlated Rayleigh
Fading Channel.In this subsection, DPSA performance under correlated Rayleigh fading channel will be discussed and compared to those of the UPSA and the static antenna scenarios.As stated earlier, the received signal in a generalized correlated Rayleigh fading channel is distributed according to a zero-mean jointly circularly normal random variable.e covariance of  can be given by   =   +  2 under  and   =  2 under , where   is the covariance of the signal component .Consequently, it can be shown that the CDF of  (see (10) for DPSA,  = ) scenario can be given by [26] ℱ DPSA  =  shown in the Appendix, if the nonzero eigenvalues of   are distinct, a closed form expression can be derived for the CDF of SLC output for DPSA scenario using the partial fraction expansion as [21] ℱ In a UPSA con�guration, the covariance of  can be given by   =   +   under  and   =   under , where   is the covariance of the signal component .Consequently, it can be shown that the CDF of  (see (10) for UPSA,  = ) is given by [26]  where   are the eigenvalues of   .Similar to the DPSA scenario, if the nonzero eigenvalues   are distinct, a closed form expression can be derived for the CDF of the SLC output as [26] ℱ where As shown in (8) the correlation between two spatial diversity branches is a function of their spatial separation.Hence, the branch correlation can be controlled by selecting the synthetic array element spacing.Figure 5 demonstrates the achievable processing gain of DPSA compared to a static antenna for the speci�c target performance requirements of   = 0.01,   = 0.99 for different values of  and antenna trajectory length, .As shown in Figure 5 for a given , increasing  initially enhances the processing gain which is due to a decrease in branch correlation arising from longer array element spacing.However since a half-wavelength antenna spacing is enough to provide uncorrelated samples and to achieve full diversity gain, a longer trajectory,   /, does not improve the detection performance.
In addition, it can be seen from Figure 5 that increasing  for a speci�c  can result in performance losses.is �rstly arises from the fact that a larger  necessitates a smaller integration subinterval, Δ = /, resulting in smaller processing gain.Secondly, increasing  for a speci�c  increases the correlation among spatial samples which reduces the potential diversity gain.As a result and similar to the uncorrelated scenario, for any , there is an optimum  for which the processing gain of DPSA is maximum.Figure 6 shows the realizable processing gain of the DPSA algorithm compared to UPSA for various values of  and .As shown in this �gure, for smaller values of  and , DPSA signi�cantly outperforms UPSA which is due to the signi�cance of polarization diversity where the potential for spatial diversity is rather limited.Note that the performance of DPSA converges towards that of UPSA where the potential for spatial diversity is more available (i.e., larger  and ).
Figure 7 is a plot of optimum  as a function of trajectory length, , for two different practical detection performance  requirements, {  = 0.05,   = 0.95} and {  = 0.01,   = 0.99}.Similar to the spatially uncorrelated scenario of Figure 3, a smaller  opt is realized by utilizing DPSA versus UPSA.

Measurements and Analysis
Experimental measurements represented in this section are based on capturing live GPS L1 C/A signals in an indoor environment.An extensive set of data was captured under controlled conditions, that is, constant antenna motion speed, in a large two and half story workshop as shown in Figure 8(a).Note that in practice a controlled motion on a known trajectory is not required and any form of user motion  can be utilized for achieving diversity.Figure 8(b) shows a sky plot of available satellites during the course of measurements.A layout of data collection site is also shown in Figure 9.
A dual polarized RHCP-LHCP active GPS antenna with a semispherical power pattern was mounted on a linear moving table to create the synthetic antenna array.e linear table was set to move the antenna with a constant speed of 2 cm/s.Since GPS signal is very weak inside buildings mainly due to the attenuation caused by walls and ceilings, a long coherent integration is required for initial acquisition.erefore, the removal of navigation data is essential in order to facilitate this long coherent integration.As a result, a reference static antenna with a clear sky view of the GPS satellites was located outdoor to assist the navigation data removal of the indoor antennas.In other words, the estimated navigation data bits from the reference antenna were used to demodulate the indoor signals.e antennas were connected to a synchronous two-channel RF frontend which downconverts the GPS signal to baseband which is eventually sampled at a 10 MHz rate aer passing through commensurate antialiasing �lters.e indoor receiver aiding process is shown in Figure 10.
Note that the complex baseband signal,   , is composed of fast and slow fading components as   =         for  =     (21) where  = √ −,   is the signal phase,   is the signal envelope aer local mean removal (the fast fading part of the signal), and   represents the local-envelope mean (the slow-fading part of the signal [27]).  is assumed to be constant over the course of measurements.In realistic cases,   ≠   for  ≠  and varies slowly due to path loss and shadowing [28].Nonetheless, the received signal must be normalized to its envelope local mean in order to �lter-out the slow-fading component of the signal [11,14,27].is normalization is essential since it maps the received signal samples into samples of an ergodic process and consequently validates the use of time average instead of expected value operator [21].Hence, herein the received signal was normalized to its amplitude local mean in order to mitigate shadowing and path loss.In addition, a three-point moving average �lter was applied to the processed signal to mitigate high frequency noise and digitizing errors [11].e magnitude of complex correlation coefficient between the RHCP and LHCP branches was then calculated given the standard correlation coefficient estimation as where ⟨⋅⟩ denotes a time average and "" is a matrix Hermitian.e results are demonstrated in Table 1.As can be seen from this table, there is negligible correlation between RHCP and LHCP signals.Note that the RHCP and LHCP signals are both coherently demodulated by the same clock such that the measured correlation is mainly a function of antenna polarizations.
Figure 11 shows the magnitude of complex correlation between spatially separated samples as a function of their spatial distance plotted for PRN 30.As can be seen from this �gure, signal�s spatial decorrelation is in good agreement with the theoretical model of (8) which is due to the nearly spherical distribution of scatterers in the measurement site.
Figure 12 shows a histogram of signal magnitude based on measurements under ℋ1 state and its corresponding theoretical Rayleigh PDF for PRNs 4, 12, and 30.As can be seen from this �gure, the measured histogram is in good agreement with the theoretical Rayleigh fading model assumed in this paper.
An SLC combiner was employed to combine a set of  =  signal samples.Two different sample spacings of Δ =  and Δ =  were considered, and the signal snapshot period was selected to be  =  ms.e Receiver Operating Characteristics (ROC) curves of Figure 13 provide a comparison between the performances of DPSA and UPSA and those of a static antenna.As can be seen from the ROC plots of Figure 13, a signi�cant improvement in signal detection is realizable for DPSA as compared with UPSA and static antenna scenarios.As expected, DPSA always outperforms UPSA (the method which is used in [5]) due to the availability of an extra diversity branch at each spatial point.Also as can be seen from Figure 13, a smaller element spacing degrades the performances of both DPSA and UPSA due to a higher correlation between the available diversity branches.
Table 2 summarized the measured processing gain (DPSA versus UPSA) for  = 1, ,  based on a target performance requirement of   = .1,  = .99 for  = .5,.5, .75 (.eoretical values of processing gains are also computed based on (12).Both RHCP and LHCP signals (denoted by R and L in the table) were considered for the UPSA scenario.It is evident from the measurement results that DPSA signi�cantly outperforms UPSA.Note that  = 1 corresponds to a static antenna scenario.

Conclusions
In this paper the potential diversity gain resulting from the combination of using a small dual-polarized antenna that undergoes arbitrary spatial motion was considered for multipath mitigation in GNSS handheld receivers.Spatial movement of the antenna during signal collection, which is normally regarded as a detriment, was exploited to achieve additional processing gain through the added diversity.e arbitrary motion assumed was within the range of a typical usage case of the mobile device.e result was equivalent to a dual-polarized spatially distributed synthetic array that provides signi�cant diversity gain even in a spatially correlated channel where the previously proposed UPSA technique is less effective.As shown, the required trajectory length for achieving a speci�c target detection requirement is much shorter for DPSA as compared to UPSA.is makes DPSA more compatible with realistic usage mode of a handset receiver.eoretical analysis show that for a practical target detection requirement (e.g.,   = 0.99,   = 0.01), a signi�cant processing gain of about 10 d� can be achieved through utilizing DPSA, compared to UPSA.is processing gain is due to the extra degrees of freedom made available through polarization diversity in a DPSA con�guration.e detection performance of DPSA was also compared to that of the UPSA under different channel conditions, and it was shown that the detection performance of DPSA signi�cantly outperforms that of the UPSA especially when the channel is spatially correlated.Consequently, an extensive set of measurements utilizing indoor GPS signals were performed to verify the assumption of the fading channel and the estimated processing gains.e theoretical evaluations and the experimental analysis proved to be in good agreement.
It is evident from these results that DPSA outperforms the UPSA by exploiting the additional polarization diversity available through the dual polarized antenna.is additional processing gain maps directly into an enhancement in GNSS signal acquisition. Appendix

15 F 2 : 15 F 3 :
=   are the eigenvalues of the covariance matrix,   .Note that    =    arises from the assumption of identical distribution of the signal in either polarizations.As Processing gain of DPSA as compared with a static antenna for   = .. opt for different values of   and   .

F 4 :
Processing gain () of DPSA compared to UPSA for   = 0.01 and different values of  and  0 .

F 8 :
Data collection location and test setup and sky plot of available satellites.

F 9 :F 10 :
Floor plan of the measurement site.Navigation bit aiding process.

T 1 :
Measured magnitude of complex correlation, |  |, between the received RHCP and LHCP signals aer 100 ms of coherent integration.sphere of scatterers model Measured-RHC polarized signal Measured-LHC polarized signal F 11: e magnitudes of complex spatial correlation of the received RHCP/LHCP signals (PRN 30).

F 12 :
Measured and theoretical Rayleigh PDF of signal magnitude.
De�ne the th subinterval that extends over the time interval of        + Δ where   is the starting instance of the th subinterval.In the UPSA scenario, the signal captured in each subinterval is correlated with    resulting in a set of  spatial array samples.e th synthetic array sample is denoted by   and given by  | 2   1.It is assumed that the bandwidth of    is signi�cantly larger than the inverse of the channel coherence time such that (3) can be approximated by | appendix derives the CDF of summation of correlated central chi-squared random variables with four degrees of freedom based on the partial fraction expansion as (A.3)