Abstract and Applied Analysis

Volume 2013 (2013), Article ID 972371, 14 pages

http://dx.doi.org/10.1155/2013/972371

Research Article

A Dynamic Cooperative Scheme with Multiple Antennas for Indoor Mobile Robot Localization

Chung-Liang Chang and Bo-Han Wu

Department of Biomechatronics Engineering, National Pingtung University of Science and Technology, No. 1 Shuefu Road, Neipu, Pingtung County 91201, Taiwan

Received 12 June 2013; Accepted 28 August 2013

Academic Editor: Chang-Hua Lien

Copyright © 2013 Chung-Liang Chang and Bo-Han Wu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper proposes a spatial cooperative diversity and decision-making technique to enhance signal detection and indoor mobile robot positioning performance of a global positioning satellite system (GNSS) receiver. Though the adaptive antenna array technique in early research could effectively promote antijamming freedom, overcome time-varying system, and mitigate narrowband and wideband interferences, factors such as the decrease of signal magnitude caused by obstacles (especially in the indoor environment), multipath, and blanking effect caused by the change in antenna direction with the motion of mobile robot can degrade the detection and interference mitigation performance of GNSS receivers. This paper aims to develop a dynamic cooperative scheme to proceed with the switch, selection, combination, and optimization among antennas. In addition, a signal processing experimental platform is also established to receive actually indoor GNSS signals for verification. The proposed scheme is capable of effectively promoting the postcorrelation signal-to-noise ratio (SNR) capability of a GNSS receiver under the indoor environment.

1. Introduction

The global navigation satellite system (GNSS) signal can be considered as a combination of signal power, carrier, spread code, and navigation message, which includes intended and unintended interferences. Once the receiver gets the navigation messages of more than four satellites, then the user position, velocity, and timing can be calculated through internal core processing of the receiver. However, how to overcome the external uncertainty factor (such as interference, and noise) is the problem concerning most of the scholars and, particularly, how to make indoor positioning efficient through the use of the receiver. Generally, the present Global Positioning System (GPS) does not function well in indoor or urban canyon multipath environments, which are regarded as difficult environments. It is hard to copy the attenuated and multipath faded signals caused by transmission around and through concrete and steel building structures.

Klukas et al. analyze the attenuation level of navigation signal through the signal reflected by different materials in indoor environment. The results demonstrate that the attenuation level of signal is between 0.5 and 23 dB [1]. Many scholars provide techniques to improve signal detection probability in indoor environment [2, 3]. The most common method is so called assisted-GPS (A-GPS) [4]. This method provides positioning information from a remote server in the cellular base station through wireless network. Nevertheless, GNSS suppliers expect the indoor positioning system to operate independently without wireless network and also maintain high positioning accuracy. It is difficult to obtain high positioning precision without any assisting device (such as WiFi, UWB, and ZigBee). In the current technique, the use of multiple antennas array combination method seems promising to be able to achieve indoor positioning. Tasi et al. propose Laplacian energy distribution model to derive spatial fading function under circular antenna array [5]. This attenuation signal model is often adopted to simulate and analyze communication channel response in indoor environment.

In early communication systems, the most common method was the use of antenna array with space-time coding or cooperation to increase system capacity under fading environment [6]. This technique was applied to interference suppression in multiinput multioutput (MIMO) system or navigation system later [7, 8]. In recent years, some antenna array processing techniques with intelligent and switching algorithms are proposed against the narrowband/wideband interferences, mutual coupling, and noise [9–11]. However, the application of these methods to GNSS positioning is still at initial stage.

Dietrich Jr. et al. first propose the concept of antenna array in indoor environment. This method utilizes spatial diversity characteristic for antenna array to enhance the signal-to-noise (SNR) ratio [12]. After that, Chiau and Gao not only use this technique in MIMO system to increase the channel capacity but also apply it to mobile device for future Galileo navigation system [13, 14]. The core concept is that each antenna maintains a certain distance from other antennas, and then the independent fading signal is received from each antenna. Thus, the multiple antennas combination method can be applied in this system to increase signal detection probability and SNR.

Broumandan et al. utilize the concept of diversity technique to achieve 7 dB signal processing gain for mobile navigation device in indoor environment [15].Zaheri et al. propose spatial diversity technique using two antennas, which results in 5 dB diversity gain enhancement and effectively suppresses multipath fading effect [16]. This paper proposes dynamic cooperative scheme based on the concept of antenna diversity characteristics. Through the antenna combination and decision logic process, the proposed scheme improves the signal detection probability under indoor environment. Unlike the diversity method proposed earlier [17, 18], this paper provides a dynamic switching and selection process for antenna array system. The experiments are also conducted with a mobile robot in indoor environment, and analysis and evaluation are demonstrated in this paper.

2. Methodology

The goal of antenna diversity technique is to mitigate the fading and improve the overall quality of the radio link. The basic principle is that the receiver combines more than one copy of a satellite signal where each copy is received through a different branch. As for GNSS applications, the receiver must be able to detect a weak navigation signal with −140 to −160 dBm power strength under outdoor environment [19]. Once the signal is masked or blanked, it is difficult for the receiver to position the user location due to signal attenuation. Thus, some diversity methods, including frequency diversity, time diversity, and antenna diversity, are proposed to solve this problem. In the following section, the dynamic cooperative scheme is demonstrated as a measure of effectiveness for multiple antennas.

2.1. Signal Model

The proposed scheme utilizes different types of antenna array combination to form a specific spatial relationship, which results in the success of navigation signal acquisition. The system model is considered to describe the relationship between each antenna location and reference point in the space. Figure 1 depicts the orientation of the robot, which is described in terms of Euler angles. Let and be the positions vectors in the body-fixed frame (-frame) and earth-centered earth-fixed frame (-frame), respectively. Meanwhile, let and be the orientation positions in the body-fixed frame (-frame) and earth-centered inertial frame (-frame), respectively. The transformation matrix can be utilized to transfer the -frame to the -frame system. Thus, the positions vectors are , where represents the transformation matrix from the -frame to the -frame. As the robot moves, both the antenna position and the rotation matrix are also changed accordingly [20]. The location of each antenna element in the body-fixed frame is assumed to be at , and the pointing angle of the antenna is along , where is the lever arm from the center of the gravity of the vehicle to the center of the antenna, is the angle of the th antenna initial location, and is the radius length. depicts the transpose operator. Thus, assume the received signal from the reference point under indoor environment can be described as follows:

where denotes the reception power of the th satellite signal, are the delay time from the th satellite to the reference point, represents the navigation data bit, is the spread code sequence, is the carrier frequency, is the carrier phase offset, and . Suppose that there are antenna elements employed for the reception and processing of GNSS signals. Then the th received signal can be described as the signal with a propagation time delay:

where depicts the propagation time delay from the th satellite signal to the reference point. The value can be transferred to phase delay. Thus, let and , where depicts the inner product between “” and “”. The unit vector is from the reference point to the th incoming satellite signal, coordinated in the earth-centered earth-fixed frame (-frame), and is the speed of light in m/sec. is the number of the visible satellites in the sky. Under the signal attenuation environment, each satellite signal is received by the receiver antenna from different directions with a certain angular spread in space. The antenna correlation gain and direction gain of weight factor decide the degree of the signal attenuation. The symbols and depict the th and th antennas, respectively. In addition, the phase delay factor caused by multipath effect is also considered. Thus, the signal model in the fading environment is as follows:

where and depicts noise term with Gaussian distribution, such as thermal noise and hardware noise.

Figure 1: Definition of coordinate systems and locations of antenna elements.

2.2. Prediversity Process

The detection performance of weak indoor navigation signal can be enhanced through multiple antennas combination method, such as polarization diversity, spatial diversity, and pattern diversity.

Polarization diversity method combines pairs of antennas with orthogonal polarizations, which can immunize a system from polarization mismatches that would otherwise cause signal fade. The antenna correlation can be depicted as complex and envelope correlations. The envelope correlation coefficient is related to the complex correlation as shown below [21]:

where , , and denotes the electric far field of the antenna from source [22]. is the ratio of the averaged vertical power to time averaged horizontal power in the fading environment in linear form. and are the azimuthal angular and elevation angular density functions of the incoming plane waves. In general, a uniform distribution is a reasonable assumption for the azimuthal angular density function as assumed earlier [23], when user of a mobile terminal moves along a random route. The angular density functions in the elevation direction are not uniformly distributed. Figure 2 shows the electric far field of the antenna.

(a)

(b)

Figure 2: (a) The relation of angular coordinates to Cartesian coordinates. (b) Antenna gain pattern.

In (4), is for ideal case. If the correlation coefficient is greater than zero (), then the diversity gain will be reduced. The analysis shows that, when the correlation is not too close to unity or , the degradation of the diversity gain due to envelope correlation is given by the gain degradation factor (GDF) of the following equation:

It is suitable to adjust the antenna spacing and the number of antennas to achieve access to the independent fading signal by using the diversity scheme. However, smaller antenna spacing will result in mutual coupling effect and degrade the effectiveness of the spatial diversity. Under the optimal antenna spacing arrangement, the use of pattern diversity method achieves different pattern shapes and thus is conducive to access a particular direction of the signal [24, 25]. The antennas serve as parasitic elements to each other and alter their patterns to allow signals to be picked up in different directions. This technique is so called “beamforming.” Assume the complex weight vector with unit norm ; then the gain of antenna system is shown as follows:

where channel gain is the same as in (3) and is the convariance matrix of the channel response vector . Assuming independent fading environment, , then array gain is equal to 1 (), which depicts the array providing no gain over a single antenna. If the fading is fully correlated between the antennas, the matrix and , where denotes expected operator. Thus, the gain of antenna system is ().

2.3. Postdiversity Process

After obtaining independent fading paths, the next criterion in diversity technique is antenna combinations. The signals received at each branch can be combined in different ways to mitigate the effects of fading. If the two signals are uncorrelated, it is rare that the two fading signals will be in a deep null at the same time. Thus, the combined signal creates a higher postcorrelation SNR (post-SNR) at the output compared to a single branch resulting in a diversity gain.

2.3.1. Multiple Antennas Combination Methods

It is necessary to synchronize the clock timings for each branch before using the antenna diversity technique. The decision module relies on the difference in the signal detection probability of each branch and performs the antenna selection, switch, and combination. Figure 3 depicts the system diagram including preprocessing module, postprocessing module, switched module, and logic decision module. Among them, the mode selection is built in logic decision module, which can perform different antenna combinations.

Figure 3: Cooperative antenna array system.

This proposed scheme can combine the received signal of each antenna through different methods, such as switched combination (SWC), selection combination (SC), equal gain combination (EGC), and maximum ratio combination (MRC) [26, 27].

The concept of SC method is that the rear of each antenna is equipped with a radio frequency (RF) front-end module. Each module calculates the SNR of its corresponding received signal and then selects the highest SNR as system output [27]. The SWC method quite resembles the selection method, but it is different in the priori setting of signal to noise ratio threshold employed in this method. If the SNR of selected antenna is below threshold value, then the system will switch to another antenna. The drawback of SWC is that the SNR of switched antenna may not be the highest one among them. The SC and SWC methods can only have one antenna for system output; and remaining antennas are not used. To effectively make use of all antennas and enhance signal detection performance, this method sums the received data of all the antennas to output. However, the summed signal may appear incoherent due to the inconsistent signal phase of each antenna. Thus, the rear of each antenna is multiplied by a phase compensator to render the received signal phase coherent. EGC is a special case of the maximum ratio combination. The antenna weights are set as one. Consequently, there is no need to adjust the weight of the individual signals before adding them, which indicates that all the antenna signals are weighed equally in spite of their SNRs. Thus, the total SNR applied to the detector after the diversity reception of selected channels is given by

where represents the total noise power expressed as the weighted sum of the noise powers of all received signals in the diversity reception. is the average noise power at each branch.

As for EGC, , the average SNR is thus given by

Assume the noise in each of the branches is independent of noise power . Under Rayleigh fading channels, when ,

where is the variance operator. Under , assuming that each of the diversity branches receives uncorrelated fading signal, one gets

Thus, the average output SNR in (8) becomes

where depicting the average SNR per branch of diversity. In fact, EGC is an exception of maximum ratio combination. The received signal of each antenna module is processed through weighting, phase compensation, and sum. The reason for weighting is because the received SNR of each antenna is not the same. Low SNR of certain antennas will degrade the signal detection performance of the whole system. Thus, each signal is processed through weighting prior to phase compensation in order to maximize the SNR of the whole system. A diversity antenna with a lower SNR is given a lower weighting, and a diversity antenna with a higher SNR is given a higher weighting. In EGC, the weighting for all is set as 1. Thus, (7) is shown as

where is the envelope that is applied to the detector. is the desired signal at the th antenna, and is the weight of the th antenna. Rappaport showed that the maximum SNR of the detector is expressed as [7]

where and indicate the SNR and the average noise power at the th antenna, respectively. Thus, it can be said from (13) that the output SNR applied to the detector utilizing selected branches equals the sum of SNR of various antenna signals. The probability density function of the output SNR can be shown as a chi-square random variable with degree of freedom with average SNR being equal to , and it can be demonstrated as [28]

The cumulative distribution function can be given by

where is the detection threshold.

2.3.2. Decision-Making Logic

In indoor environment, the signal is susceptible to multipath, blanking, and narrowband interference. A decision of all the interferences in spatial cooperative diversity technique depends on signal probability distribution model [11]. Related decision logic is a technique depicted in the literature review [29]. Effective decision logic and the selection of proper spatial cooperative diversity technique can enhance signal detection performance. In antenna diversity combination, input signals, after despreading denoted as , received via diversity branches, are considered. In the detection process for GNSS signals, the performance can be evaluated through adoption of an individual binary hypothesis testing, defined as and . corresponds to the scenario, when the incoming signal does not correspond to the hypothesis related to local replica , whereas, under state, it is supposed that the frequency and code delay of the incoming signal are synchronized with the local replica. The detection structure for such binary hypothesis test results in the Likelihood Ratio Test (LRT), when the additive noise is zero mean Gaussian distribution [14]. The LRT for an antenna diversity structure consisting of branches can be considered as received signal estimator followed by a correlator. The different types of antenna combination methods can be performed under condition where correlation coefficient is equal to zero [30]. Assume the incoming signal includes desired signal component, and then the signal flag is set as one. Otherwise, the signal flag is set as zero. denotes the total times of , which decides the type of antenna combination method. The definition of is the same as in (7). Table 1 demonstrates the decision-making logic process for the proposed scheme. The design criterion of decision-making logic depends on value that corresponds to the four operation modes including SC, SWC, MRC, and EGC.

Table 1: Decision-making logic.

Combination mode

Mode 1: “SC”
,

Mode 2: “SWC”
,

Mode 3: “MRC”
,

Mode 4: “EGC”
,

2.3.3. Correlation Coefficient Evaluation

Each of the multiple antennas is separated spatially to achieve independent fading paths. The relative phase of the fading signal is obviously different at any two antennas because of the spatial separation, that is, antenna spacing between both the antennas. Thus, the correlation coefficient can be obtained through (2), which can be calculated by using the two signal correlations. Assume the received signal envelopes for the two antennas are and , respectively, where and are the received complex signals on the first and second antenna, respectively, shown in (2). The correlation between the two antennas is given by [1]

where describes time delay from antenna 1 to antenna 2. represents the autocorrelation function of the received envelope . The detailed deduction of (16) is elaborated by Klukas et al. [1]. In (16), the is given by

where and are the real and imaginary parts of , respectively. Supposeing imaginary parts of correlation are equal to zero, then ; meanwhile , substituting above value into (16). Thus, the correlation in (16) can be calculated as

Utilizing the Bessel function expression, the can be expressed as

where denotes the spacing between first antenna and second antenna, is the wavelength of the incoming plane wave, and represents zero-order Bessel function of first kind [31, 32]. indicates the total received envelope power. Using (19), the autocovariance function is given by

After normalization, the antenna correlation between the first and second antenna elements is expressed as

The envelope correlation between the two receiver channels is calculated and compared with theoretical values found by Stüber [31]. The resulting output is coherently in-phase and quadrature integrated prior to temporal or spatial processing after code and carrier removal. Assume a coherent integration output of seconds for the first and second antennas is and respectively; in that case their envelopes are and , respectively; then the envelope correlation between the two antennas is depicted in terms of mean-removed versions of and as

The depiction in (22) holds true for a uniform angle of arrival distribution in azimuth and omnidirectional receiving antennas matched to the polarization of the incoming radio wave. The symbol depicts the variance calculation operator with deviation . The envelope correlation of signals in (22) received by any two antennas suffers from several factors such as the mutual coupling, radiation pattern, direction of arrival of the radio frequency (RF) wave, and antenna gain.

3. Performance Criteria

The proposed scheme combines the polarization diversity and pattern diversity to form the desired array gain pattern. Because the right-hand circular polarization (RHCP) or linear polarization is used in common GNSS antenna, the antenna spacing is only an arrangement factor used to avoid the effect of mutual coupling and cross-polarization [14]. Therefore, this section illustrates the performance evaluation criteria of the proposed system.

3.1. Signal Detection Probability

As detection performance of the spatial combination technique is a function of only the probability of false alarm and the probability of detection , thus, assume that, under hypothesis, the frequency and code phase of the incident signal are synchronized with locally generated signal. Otherwise, the signal is not presented (under hypothesis). For single antenna case, the detection statistic is a chi-square two-degrees-of-freedom central probability density function random variable for both and conditions. As a result, the false alarm probability and detection probability are shown as follows [33]:

where is probability function. As detection performance of a multiple antenna system, the EGC is explained with uncorrelated channel condition: that is, ,

Figure 4 demonstrates the distribution of signal detection probability using different switch methods under . The diagram shows that the signal detection probability of proposed scheme can be regarded as upper and lower boundary range. That is because the proposed decision logic is based on Table 1. Thus, different combinations of the antennas are selected and adopted.

Figure 4: The probability of detection for different spatial combination method.

3.1.1. Diversity Gain

The performance analysis of spatial cooperative diversity technique can be assessed through calculated diversity gain. The evaluation standard lies in the difference between the acquired SNR of single antenna system and that of the spatial cooperative diversity technique. For example, in multipath environment, the channel is Rayleigh. Then, the cumulative distribution function (CDF) can be calculated as suggested earlier [34]:

where is the normalized SNR or threshold, is the average of SNR, depicts the calculated instantaneous SNR, and represents the probability of SNR below threshold . The choice of antenna indicates that the signal detection probability of whole antenna system is times. Assuming that the antenna branches have independent signals and equal mean SNRs, the probability of all branches having a SNR below equals the probability for a single antenna branch increased to the power as

Usually is set as 1%. At low instantaneous SNR (), (26) approximates when is given; diversity gain can be calculated as

where represents the SNR using cooperative diversity technique, whereas stands for the SNR of the first antenna without using cooperative diversity technique. The result of diversity gain can evaluate the signal detection performance of the whole system. In GNSS application, the diversity gain measurements initially involve the estimation of SNR at the output of the coherent integration. Thus, here the post-SNR for the single and multiple antenna channels combination method is regarded as the deflection coefficient. Thus, (27) is rewritten as follows [33]:

where represents the test statistic. Because the detection statistics are no longer modeled as Gaussian random variables, the deflection coefficient offers only a rough measure of diversity gain.

Figure 5 shows the reduction of the probability of fading below a given threshold or normalized SNR when increasing the number of antennas. This figure also shows that, for % reliability level, there are diversity gains of 9.2 dB and 7.3 dB for the two-antenna and four-antenna selection combiners, respectively.

Figure 5: Probability of Rayleigh fading signals for different numbers of diversity antennas.

3.1.2. Correlation Coefficient

Furthermore, the spatial relationship between any two antennas is used to calculate the spatial correlation. Ideally, the spatial correlation function between any two antennas is low, and the difference of calculated SNR among each antenna model is small. The square value of complex correlation coefficient is quantified in the following equation:

where and are zero mean complex values of the correlator output of each received signal channel. “*” is the complex conjugate operator. Higher correlation coefficient in spatial diversity may result in lower diversity gain.

4. Test Hardware

This chapter describes hardware setup that is the use of multiple antennas for different spatial combing including switching, selection, equal gain, and maximum gain method. The establishment of antenna array system is also demonstrated in this chapter.

4.1. Hardware Specification

Figure 6 shows the mobile robot appearance. The robot is equipped with motor module, and the use of motor and gear is adopted to control the rotation speed of the antenna array. The RF front-end module and NI CB-2162 are mounted on computer side, and ends of the antennas are connected to the GPS501 hardware box. Four commercial GPS front-ends are adopted to filter, down-convert, and digitize incoming RF signal to obtain digital samples. A master clock is utilized to control the synchronization among four different channels so that the samples can be acquired simultaneously. The digitized samples are acquired at the rate of 5.714 MHz. The digital samples are obtained using a PCI-based data acquisition board (PCI-6536) through a data acquisition software (NI LabView). The data samples from each antenna element are stored and postprocessed for interference mitigation analysis. Figure 7 shows the structure of array system hardware. Firstly, switched board (NI CB-2162) acquires the “SGN” and “MAG” binary digital signal from the four RF front-end modules. After that, the 8 bit signal, which consists of four groups of “SGN” and “MAG” signal, is sent to an NI PCIe-6536 high-speed digital I/O cards and then stored in the personal computer. The NovAtel reference antenna is adopted to estimate the Doppler frequency, code phase, and navigation data bits, and these estimates are utilized to demodulate the GPS digital intermediate frequency (IF) data with respect to the indoor antenna.

(a)

(b)

Figure 6: Platform of the moving vehicle.

Figure 7: Schematic of the experiment setup.

4.2. Experimental Setup

The two selected indoor locations are situated at biomechatronics engineering building, namely, the BIOME building of Pingtung University of Science and Technology (NPUST). Figure 8 illustrates the schematic plan view of the second floor of the building and the indoor test locations. The first location is situated in a corridor, and the other location is inside a room of the 2nd floor. The indoor test location number 1 inside the corridor has its southern view blocked with the building wall, but northern of the antenna is open through the glass windows. For verification purposes, a second location on the same floor of the same building is also investigated. The patch antenna indoors are facing the glass window panes towards the western direction. Figure 9 depicts the sky-view plot at the time as recorded by the reference receiver (NovAtel Inc.) located at the rooftop. The roof reference antenna is stationary; the indoor patch antennas are mounted on a rotation platform. The antenna spacing is set as 1 λ. The experiments are divided into static and dynamic scenarios, respectively. After that, the received data is stored and processed using MATLAB software for signal acquisition and tracking.

Figure 8: Experimental test location in BIOME building, NPUST.

Figure 9: Sky-view plot of GNSS.

5. Experimental Results

This chapter analyzes the performance of proposed scheme through experiments. The initial step is to have a mobile robot to move to the designated position. Upon arrival, the mobile robot stops moving and starts its experiment of static signal reception. The duration of static data is half an hour. The dynamic rotation speeds are 1, 4, and 20 rpm. The reception time in terms of dynamic data is half an hour. The above data are utilized for the performance analysis using the following three steps. First of all, the indoor fading scenario for all of the antennas is demonstrated, where the distribution characteristics indicate the indoor fading scenario with respect to envelope correlation and Rician factor. The envelope correlation is compared with the theoretical correlation coefficient as the antenna spacing varies. Secondly, the effects of cooperative combination are depicted. The effects of EGC and the proposed scheme are particularly compared. Besides, the comparison with single antenna case is also conducted, and post-SNR improvements for the proposed decision scheme are also described. At last, analyses with receiver operating characteristics (ROC) curves with regard to improvement in detection probability in contrast to that of the single antenna processing are demonstrated.

5.1. Fading Environment

It is necessary to extend coherent integration time to obtain high-sensitivity operation due to weak and severe nonline of sight (NLOS) signal conditions in indoors. Coherent and noncoherent integration times are 8 ms and 1 sec, respectively. As an example of indoor fading characteristics, over the time for PRN 3 corresponding to the dataset with 1 λ antenna spacing is shown in Figure 10. The average for all of the antennas is roughly 25.1 dB-Hz. Figure 11 demonstrates the acquisition result of a weak GPS signal (PRN 3) from each antenna with 1 λ antenna spacing, where the correlation peak is also zoomed in plots of Doppler frequency and code offset.

Figure 10: plots for the two antennas at the indoor fading environment.

(a)

(b)

(c)

(d)

Figure 11: Weak GNSS signal acquisition results from each antenna branch.

In Figure 12, the average antenna correlation acquired through experimental measurements supports the theoretical envelope correlation prediction demonstrated in (12). This figure illustrates that the measured antenna envelope correlation is lower than the theoretical correlation curve, especially in terms of smaller antenna spacing. This is because the indoor patch antenna employed in this paper is not an ideal one, and the radiation pattern of the indoor patch antenna is not omnidirectional. Figure 12 depicts that the measured squared magnitude envelopes maintain a similar shape for all of the available pseudorandom noise (PRNs), but they differ with respect to the actual correlation value. This is because the actual antenna envelope correlation counts on the angular multipath spread in addition to the antenna characteristics.

Figure 12: Squared envelope correlation as a function of antenna spacing, indoor location.

Figure 13 depicts the deflection coefficient measures for the single and proposed decision schemes corresponding to the indoor data set of 0.5 λ antenna spacing for PRN 19. The cumulative distribution function is also plotted for the same dataset. The deflection coefficient for the spatial combination techniques has gone through a lower number of fades in contrast to the single antenna processing. Correspondingly, the deflection coefficient values increase drastically for a smaller value of cumulative probability. The performance gap between single antenna and diversity antennas and the performance improvement results from the diversity gain with cooperative signal combination are also shown in Figure 13.

(a)

(b)

Figure 13: (a) Deflection coefficient as a function of time. (b) Probability of detection versus deflection coefficient, PRN 19,  λ.

Figure 14 shows the results of a research on the achievable post-SNR for the different antenna spacings averaged over all PRNs. A coherent integration period of 5 ms and a total observation period of 1 s for noncoherent summation are taken into account. Both EGC and proposed scheme show better performance than single antenna performance. Nevertheless, proposed scheme presents a slightly superior performance in contrast to that of EGC. This is because the performance of proposed scheme is affected by the correlation coefficients between the antennas. The correlation coefficient is reduced by increasing the antenna spacing, and both EGC and proposed scheme present similar performance.

Figure 14: Averaged post-SNR as a function of antenna spacing (indoor location;  ms, noncoherent integration time is 1 sec, PRN 19).

5.2. Post-SNR

Table 2 illustrates the average post-SNR results that use different types of cooperative diversity method under the worst case. In Table 2, the proposed scheme with 1 λ antenna spacing provides 2.3 dB post-SNR performance enhancement in contrast to EGC method under the static experiment (without antenna array rotation). The performance of proposed scheme is the same as that of MRC method. In addition, only the use of proposed scheme with 1 λ antenna spacing and 20 rpm can acquire the desired signal under the dynamic experiment (with antenna array rotation).

Table 2: Results of static and dynamic experiment in worst case ( ms, noncoherent integration time is 0.5 sec).

Average post-SNR (Th: 17 dB)
Single antenna case
SWC
MRC
EGC
Proposed scheme

1st antenna

(North)
2nd antenna

(East)
3rd antenna

(South)
4th antenna

(West)

Static experiment (PRN 19)

0.5 λ
×
17.3
17.3
17.7
17.3~17.7
19.5
18.3
19.5

1 λ
17.9
17.8
17.9
18.3
17.8~18.3
20.9
18.6
20.9

2.5 λ
×
×
×
×
×
19.1
18.0
19.1

Dynamic experiment (PRN 19; antenna spacing: 1 λ)

1 rpm
17.7
17.5
17.7
17.9
17.5~17.9
18.8
18.5
18.8

4 rpm
17.7
17.3
17.1
17.0
17.1~17.7
17.5
×
17.5

20 rpm
×
×
17.4
17.1
×
18.7
×
18.7

“×”: unlocked.

5.3. Detection Performance

Performance evaluation with respect to the deflection coefficient is only a rough account of the improvements, not an accurate measure. Thus, the detection performance has been conducted based on the indoor test data in the form of ROC. The coherent integration period of 100 ms results in good probability of the detection. So as to acquire test statistics for , that is, under the presence of a signal, the correct code phase estimate is adopted and the corresponding detection output is computed. Incorrect code phase is employed to obtain the test statistics under condition. The average post-SNR values for the PRN processed in the dataset for the static experiment are in the range from 17 to 20 dB-Hz. Figure 15 depicts the CDF of the detection output for single antenna and cooperative antennas at 1 λ antenna spacing. The detection output was based on 100 ms coherent integration time and 6 sec noncoherent integration time. Under and hypotheses, single antenna performance and cooperative combination with EG and proposed scheme are demonstrated. Figure 15 readily presents the advantages of adopting multiple antennas for signal detection purposes. For any given percentage of time, the proposed scheme presents slightly superior detection performance than that of EGC.

Figure 15: CDF of detection output in the case of four antennas,  ms, and noncoherent integration time is 6 sec, in the indoor location.

6. Conclusion

This research demonstrates the dynamic cooperative scheme to enhance the signal detection performance with regard to multiple antennas under indoor environment. Simulation result illustrates that the use of multiple antennas with decision logic scheme can effectively improve post-SNR performance in contrast to single antenna case. Meanwhile, real-time data reception is conducted under indoor environment. It is successfully verified that the proposed scheme indeed effectively enhances signal detection probability and improve post-SNR. In the future, this technique can combine beamforming algorithm to be applied to the multinavigation satellite signal reception and the interference mitigation system.

Acknowledgments

Many thanks are due to the editor and reviewers for their valuable comments to refine this paper. Besides, the authors would like to thank the National Science Council in Taiwan for supporting this research under Grant NSC 102-2221-E-020-007.

References

R. Klukas, O. Julien, L. Dong, M. E. Cannon, and G. Lachapelle, “Effects of building materials of UHF ranging signals,” GPS Solutions, vol. 8, no. 1, pp. 1–18, 2004.

B. Peterson, D. Bruckner, and S. Heye, “Measuring GPS signals indoors,” in Proceedings of the 10th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS '97), pp. 615–624, September 1997.

F. van Digglen and C. Abraham, “Indoor GPS technology,” in Proceedings of CTIA Wireless-Agenda, pp. 1–10, 2001.

F. van Diggelen, “Indoor GPS theory & implementation,” in Proceedings of the IEEE Position Location and Navigation Symposium, pp. 240–247, April 2002.

J.-A. Tsai, R. M. Buehrer, and B. D. Woerner, “Spatial fading correlation function of circular antenna arrays with Laplacian energy distribution,” IEEE Communications Letters, vol. 6, no. 5, pp. 178–180, 2002.

G. J. Foschini, “Layered space-time architecture for wireless communication in a fading environment when using multi-element antennas,” Bell Labs Technical Journal, vol. 1, no. 2, pp. 41–59, 1996.

H. Bolcskei, D. Gesbert, C. B. Papadias, and A. van der Veen, Space-Time Wireless Systems—From Array Processing to MIMO Communications, Cambridge University Press, Cambridge, UK, 2006.

R. L. Fante and J. J. Vaccaro, “Wideband cancellation of interference in a GPS receive array,” IEEE Transactions on Aerospace and Electronic Systems, vol. 36, no. 2, pp. 549–564, 2000.

S.-J. Kim and R. A. Iltis, “STAP for GPS receiver synchronization,” IEEE Transactions on Aerospace and Electronic Systems, vol. 40, no. 1, pp. 132–144, 2004.

C.-L. Chang, “Self-tuning synthesis filter against mutual coupling and interferences for GNSS and its implementation on embedded board,” EURASIP Journal on Advances in Signal Processing, vol. 2010, Article ID 123625, 13 pages, 2010.

C. L. Chang, “Multiplexing scheme for anti-jamming global navigation satellite system receivers,” IET Radar, Sonar & Navigation, vol. 6, no. 6, pp. 443–457, 2012.

C. B. Dietrich Jr., K. Dietze, J. R. Nealy, and W. L. Stutzman, “Spatial, polarization, and pattern diversity for wireless handheld terminals,” IEEE Transactions on Antennas and Propagation, vol. 49, no. 9, pp. 1271–1281, 2001.

C. C. Chiau, Study of the diversity antenna array for the MIMO wireless communication systems [Ph.D. thesis], Queen Mary University of London, London, UK, 2006.

Y. Gao, C. C. Chiau, X. Chen, and C. G. Parini, “Design of diversity antenna array for Galileo/GPS receivers,” in Proceedings of the European Conference on Antennas and Propagation (EuCap '06), pp. 1–5, Nice, France, November 2006.

A. Broumandan, J. Nielsen, and G. Lachapelle, “Indoor GNSS signal acquisition performance using a synthetic antenna array,” IEEE Transactions on Aerospace and Electronic Systems, vol. 47, no. 2, pp. 1337–1350, 2011.

M. Zaheri, A. Broumandan, V. Dehghanian, and G. Lachapelle, “Detection performance of polarization and spatial diversities for indoor GNSS applications,” International Journal of Antennas and Propagation, vol. 2012, Article ID 879142, 10 pages, 2012.

E. D. Kaplan and C. J. Hegarty, Understanding GPS: Principles and Applications, Artech House, Boston, Mass, USA, 2006.

E. L. Houghton and P. W. Carpenter, Aerodynamics for Engineering Students, Edward Arnold, London, UK, 4th edition, 1993.

R. Clarke, “A statistical theory of mobile radio reception,” Bell System Technical Journal, vol. 47, pp. 957–1000, 1968.

M. A. Jensen and Y. Rahmat-Samii, “Performance analysis of antennas for hand-held transceivers using FDTD,” IEEE Transactions on Antennas and Propagation, vol. 42, no. 8, pp. 1106–1112, 1994.

T. Taga, “Analysis for mean effective gain of mobile antennas in land mobile radio environments,” IEEE Transactions on Vehicular Technology, vol. 39, no. 2, pp. 117–131, 1990.

P. Mattheijssen, M. H. A. J. Herben, G. Dolmans, and L. Leyten, “Antenna-pattern diversity versus space diversity for use at handhelds,” IEEE Transactions on Vehicular Technology, vol. 53, no. 4, pp. 1035–1042, 2004.

R. Vaughan, “Switched parasitic elements for antenna diversity,” IEEE Transactions on Antennas and Propagation, vol. 47, no. 2, pp. 399–405, 1999.

S. R. Saunders, Antennas and Propagation for Wireless Communication Systems, John Wiley & Sons, New York, NY, USA, 1999.

C.-L. Chang, H.-N. Shou, and J.-C. Juang, “Spatial cooperative diversity technique in performance enhancement of GNSS receivers,” in Proceedings of the 8th Asian Control Conference (ASCC '11), pp. 1239–1244, May 2011.

M. Tarkiainen and T. Westman, “Predictive switched diversity for slow speed mobile terminals,” in Proceedings of the 47th IEEE Vehicular Technology Conference, vol. 3, pp. 2042–2044, May 1997.

A. Sendonaris, E. Erkip, and B. Aazhang, “User cooperation diversity—part I and part II,” IEEE Transactions on Communications, vol. 51, no. 11, pp. 1927–1948, 2003.

T. S. Rappaport, Wireless Communications Principles and Practice, Pearson Education International, Upper Saddle River, NJ, USA, 2nd edition, 2002.

C.-K. Park and K.-S. Min, “A study on spatial correlation characteristic of array antenna for multi antenna system,” in Proceedings of the Asia-Pacific Microwave Conference (APMC '05), vol. 3, December 2005.

W. C. Jakes, Microwave Mobile Communications, IEEE Press, Piscataway, NJ, USA, 1994.

G. L. Stüber, Principles of Mobile Communications, Springer, New York, NY, USA, 2nd edition, 2000.

S. M. Kay, Fundamentals of Statistical Signal Processing, Volume II: Detection Theory, Prentice Hall PTR, Upper Saddle River, NJ, USA, 1998.

S. R. Saunders, Antennas and Propagation for Wireless Communication Systems, John Wiley & Sons, New York, NY, USA, 2nd edition, 2007.

I. E. Telatar, “Capacity of multi-antenna Gaussian channels,” European Transactions on Telecommunications, vol. 10, no. 6, pp. 585–595, 1999.