Open Access is an initiative that aims to make scientific research freely available to all. To date our community has made over 100 million downloads. It’s based on principles of collaboration, unobstructed discovery, and, most importantly, scientific progression. As PhD students, we found it difficult to access the research we needed, so we decided to create a new Open Access publisher that levels the playing field for scientists across the world. How? By making research easy to access, and puts the academic needs of the researchers before the business interests of publishers.
We are a community of more than 103,000 authors and editors from 3,291 institutions spanning 160 countries, including Nobel Prize winners and some of the world’s mostcited researchers. Publishing on IntechOpen allows authors to earn citations and find new collaborators, meaning more people see your work not only from your own field of study, but from other related fields too.
To purchase hard copies of this book, please contact the representative in India:
CBS Publishers & Distributors Pvt. Ltd.
www.cbspd.com

customercare@cbspd.com
EE Department, King Fahd University of Petroleum & Minerals, Saudi Arabia
A. H. Muqaibel
EE Department, King Fahd University of Petroleum & Minerals, Saudi Arabia
A. S. AlAhmari
EE Department, King Fahd University of Petroleum & Minerals, Saudi Arabia
H.R. Khan
EE Department, King Fahd University of Petroleum & Minerals, Saudi Arabia
R. A. AlNimnim
EE Department, King Fahd University of Petroleum & Minerals, Saudi Arabia
*Address all correspondence to:
1. Introduction
Recently, wireless mobile communication systems have experienced a tremendous growth and became an integral part of people’s daily life worldwide. This global predominance of wireless communications has been ever more pronounced with the success of new generations of wireless communication standards that support a rich set of valueadded features in addition to basic phone services. Among these features is the possibility to offer radiolocation services whereby the mobile system (MS) position is determined by combining relevant information (such as signal time of arrival or angle of arrival) from different base stations (BSs) having radio links with the intended mobile. This positioning capability can support many services ranging from medical emergency help, security and law enforcement, ontheroad assistance, locationdependent commercial advertisement, etc. As such, mobile radiolocation has been mandated by several of the recently introduced wireless standards, and is being widely deployed by cellular network operators worldwide (Rappaport et al., 1996).
With the deployment of the 3^{rd} generation wireless cellular standards such as the Universal Mobile Telecommunication System (UMTS) (Dahlman et al., 2000), the use of wideband code division multiple access (CDMA) signals is poised to offer highly accurate positioning capabilities with timeofarrival (TOA) information owing to the fine timing resolution of the high chip rate wideband spreadspectrum waveforms used. There are however some impairments such as multipath fading, multiuser interference and noise, that can affect the performance of mobile positioning schemes. In particular, because of the need to use signal detection at several base stations in mobile positioning, the problem of nearfar interference at remote base stations (whereby a faraway mobile weak signal can be overwhelmed by strong signals from closein mobiles) can constitute a major limiting factor. In CDMA networks, this is further exacerbated by power control (Viterbi, 1995). Indeed, power control loops operate in such a way to maintain the received power from different users at the same level at their respective serving base station. However, at other nonserving base stations (not actively involved in a call with the intended mobile), the mobile received power level can be extremely low, thereby giving rise to a problem of signal hearability, which, in turn, will affect the accuracy of the mobile positioning algorithms (Gosh & Love, 1998). A thorough analysis of these different aspects under realistic channel modeling and network traffic loading conditions is therefore necessary in order to obtain an accurate assessment of the achievable positioning performance.
In this work, we study the accuracy of mobile radiolocation under nearfar interference, and show that it can vary considerably depending on the mobile link quality with the base stations involved in position determination. We first present a detailed performance analysis of mobile TOA estimation in broadband CDMA wireless cellular networks at the different base stations surrounding the mobile terminal. In most previous works, the proposed radiolocation algorithms have simplistically assumed that the timeofarrival measurements are Gaussiandistributed, with a fixed known variance that is commonly set irrespective of the actual RF channel and interference conditions or other system parameters. Moreover, the same TOA variance is typically assigned for all base stations involved in positioning, which is an inaccurate and overly optimistic assumption, as our results will show. Instead, our approach is the use of complete statistics for the mobile timing estimation errors, derived by taking into account realistic system parameters at the different base stations of interest.
Since precise TOA estimation in CDMA receivers is typically implemented by delaylocked loop (DLL) tracking systems (Viterbi, 1995), a detailed study of the DLL TOA tracking is introduced under fading and multiuser interference conditions assuming a cellular model with multiple tiers of base stations, and it is shown that the DLL tracking performance can vary widely depending on the level of multipleaccess interference and RF propagation conditions of the links between the mobile and base stations. Using the TOA data collected at the base stations involved in mobile positioning, a numericallyefficient, quasioptimum algorithm, based on Approximate Maximum Likelihood (AML) estimation (Chan et al., 2006), is presented, and a generalization of this algorithm is also derived under the realistic assumption of unequal TOA estimation errors at the different base stations.
Based on our analysis, the impact of mobile link condition and its relative position with respect to the base stations involved in its positioning is fully quantified for a number of scenarios depending on the nearfar interference environment and the level of Soft Handoff (SHO) connectivity of the mobile with the base stations. Our results show that positioning accuracy is greatly improved when the mobile station is in 2way, and particularly in 3way SHO (i.e., with two and three base stations, respectively), compared to single connectivity with only the home serving base station.
The rest of the chapter is organized as follows. First, in Section 2, we present an overview of radiolocation techniques, focusing in particular on networkbased approaches with TOA processing suitable for cellular systems. The system modeling and analysis of nearfar interference in CDMA networks is presented in Section 3, followed by a detailed analysis of the performance of delaylocked loops for timeofarrival tracking in Section 4. Then, in Section 5, the approximate maximum likelihood positioning algorithm is presented, and various illustrative examples and numerical results are discussed in Section 6, followed by summary and final conclusions in Section 7.
The concept of wireless radiolocation refers to the determination of the geographic position information of a mobile user in terms of its geographic coordinates with respect to a reference point. Wireless location is also commonly referred to as mobile positioning, radiolocation, or localization. Position location techniques can be classified into two main categories: handsetbased and networkbased. A well known example of handsetbased radiolocation is the Global Positioning System (GPS) and other similar systems (Kaplan, 1996). The other category consists of networkbased techniques that utilize the existing wireless cellular networks to obtain location information (Sayed, et al., 2005). In this work, we mainly focus on the networkbased approach which integrates seamlessly with the widely deployed mobile cellular infrastructures throughout the world. The basic architecture of such systems is illustrated in Figure 1.
The networkbased positioning schemes rely on data collected by several base stations surrounding the mobile station of interest, and can be based on the mobile’s signal strength (SS), angle of arrival (AOA), time of arrivals (TOA) or time difference of arrival (TDOA) measurements. Using these measurements, specific geometric and/or statistical signal processing algorithms are used to determine the mobile location. Hybrid methods involving more than one of type of measurements can also be used (Sayed, et al., 2005). In general, locating a mobile in twodimensions requires a minimum of three sets of measurements from corresponding base stations, although for AOA methods, two base stations are sufficient. However, in the presence of noisy measurements, statistical signal processing algorithms using data collected from multiple base stations are preferred in order to resolve the ambiguities arising from multiple crossings of the lines of position, and to improve the positioning accuracy.
2.2. TimeofArrival Mobile Radiolocation
TOAbased techniques offer several advantages compared to the other methods, including low cost and ease of use. The TOA data is readily available from timing synchronization mechanisms at the different base stations, without requiring complex hardware as with the AngleofArrival methods. In particular, with the widespread deployment of the latest 3G CDMAbased wireless cellular networks, the spreadspectrum signaling waveforms offer high time resolution and good robustness visàvis the radio channel impairments (fading, shadowing, and nearfar interference), and are therefore wellsuited to aid in achieving the required accuracy in position location services.
In the TOAbased mobile radiolocation approach, the distance between an MS and a BS is measured by finding the oneway propagation delay under direct lineofsight (LOS) propagation conditions. The TOA measurements at the different base stations are therefore directly proportional to the mobilebase distance separation. The involved base stations are assumed to have a common timing reference, with a known mobile transmission time. Geometrically, the mobile position will trace a circle centered at the base stations. By using three base stations to resolve ambiguities, the mobile position is given by the intersection of these circles as illustrated in Figure 2. It should be pointed out that, in the presence of noisy TOA data due to interference and synchronization errors, the three circles may not intersect at a single point. Therefore, the geometric approach will not produce a single intersection point, and “statistical” techniques are typically adopted to process the noisy data, as will be further discussed subsequently when we introduce the AML algorithm.
Consider a cellular CDMA network employing signaling schemes with quaternary phaseshift keying (QPSK) modulation and complex spreading (conforming to 3G/UTMS standards (Dahlman et al., 2000). At the transmitter, the complex baseband signal for a given mobile user is given by
s˜(t)=∑m=−∞∞smh(t−mTc)E1
where {s_{m}} are the complex spreading chip symbols, and T_{c} the chip duration. The information data bits are omitted here for simplicity as this will be applicable for a pure pilot signal that can be used for positioning purposes. The impulse response h(t) of the pulseshaping filter is assumed to be a rootraised cosine filter with rolloff factor 22%, as recommended in the UMTS standard. As will be seen next, the Fourier transform H(f) of the filter impulse response h(t) will have a major impact on the otheruser interference in cellular CDMA networks.
For wireless fading channels, the received signal at the output of the radio frequency (RF) receiver filter is written as
r(t)=PRe{∑i=1Ms˜(t−τi(t))ai(t)ej2πfct}+n(t)E2
where P is the received signal power, M is the number of resolvable multipath signals, with a_{i} and τ_{i} denoting the ith path complex Gaussian tap factor (with a Rayleighfading magnitude) and its propagation delay, respectively, and f_{c} is the carrier frequency. The noise signal n(t) models the total noiseplusinterference terms and is assumed to be zeromean Gaussian random process. For the different resolvable multipath signal epochs, we mainly focus on the first arriving signal which will be tracked to estimate its TOA. The combined multipleaccess interference (MAI) terms from both samecell and othercell users can be modeled as being Gaussian distributed, which is a valid assumption for a large number of users. In this case, it can be shown that the total composite power spectral density (PSD), I_{0}, that captures the effect of both thermal noise and MAI terms (assumed to be statistically independent) is given by (Viterbi, 1995)
where N_{o} is due to the thermal noise component and the factor ρ_{ο} represents both samecell and othercell MAI terms. The function H(f) denotes the Fourier transfer of the chip shaping filter impulse response h(t). The MAI PSD term is typically the dominant factor in CDMA systems, and will depend on network loading, fraction of intercell to intracell MAI, channel pathloss and shadowing models. For typical system parameters of interest, it is found that the ρ_{o} MAI factor is approximately 1.6(K_{u}1), where K_{u} is the number of users per cell (Viterbi, 1995).
3.2. NearFar Interference and Soft Handoff Impact
We consider a cellular network with a central cell and two tiers of surrounding cells as illustrated in Figure 3, where mobiles are assumed to be uniformly distributed across the coverage area.
The mobile received power at a given base station, BS_{i}, is multiplied by an attenuation factor, α, that reflects distance path loss, p(d), and lognormal shadowing factor, ξ, according to (Viterbi, 1995):
α(dBSi,ξBSi)=p(dBSi)10ξBSi/10E4
The pathloss factor is assumed to follow the model:
p(d)=10nlog10(d)E5
where n is the path loss exponent and d the mobile’s relative position with respect to the base station. For a given mobile location, shadowing visàvis the different base stations is typically assumed to be partially correlated, lognormallydistributed, and given by (Viterbi, 1995):
ξBSi=aξc+bξiE6
where ξ_{c} and ξ_{i} denote the common and independent terms, respectively, and a^{2}+b^{2}=1.
With closedloop CDMA power control, a given mobile will have the same received power at its home (serving) base station compared with other intracell mobiles. On the other hand, at the neighboring base stations, the mobile will be received at much lower power due to the socalled nearfar problem (which cannot be mitigated due to the lack of “intercell” power control). However, when the mobile is in softhandoff with other neighboring base stations, its received power is relatively close to that at the home cell, and this greatly improve radiolocation accuracy, as will be discussed next.
Since TOA estimation accuracy depends on the timing synchronization mechanism, which is in turn affected by the received interference levels at the different base stations, the nearfar interference at the nonserving base stations will have a major impact on the final mobile positioning accuracy. To illustrate this point, we consider a system where the mobile is served by the center base station BS_{1} and is radiolocated using TOA data from three or more base stations BS_{1}, BS_{2},…, BS_{7} (sorted in a descending order from BS_{1} that receives the highest average received power). The nearfar interference impact can be conveniently assessed by defining the ratio of its average received power at BS_{i} compared to BS_{1} (for which the power will be normalized to 1, and used as a reference value). We then define the following:
βi=Pi/P1E7
where P_{i} is the received power at BS_{i} and β_{1}=1≥ β_{2} ≥ β_{3} ≥ ... ≥ β_{7}. A wide fluctuation in the β factors is expected depending on the mobile position relative to the base stations of interest. It should be noted that, due to power control, all mobiles will be received at equal power (β_{}=1) at their respective home serving base stations, but much lower values of β are expected at faraway base stations (because of the nearfar problem). This however will depend on the relative position of the mobile with respect to the other base stations (i.e., its proximity to the cell border). Soft Handoff (SHO) is one of the salient features of CDMA cellular systems which allow the mobile to be simultaneously connected to more than one serving base station. In fact, the possibility of SHO calls enables a stronger signal reception at multiple base stations, and this will in turn improve positioning accuracy.
To further investigate this point, we consider different scenarios denoted by Cases 1, 2, and 3, respectively. Case 1 refers to a mobile in close proximity (within half the cell radius, R) of its serving BS_{1}, with a signal at least 10dB above that at the other nearest two base stations. Case 2 represents a twoway soft handover scenario, with the mobile power at BS_{2} within 3dB from that at BS_{1}. Finally, Case 3 corresponds to a 3way soft handover situation where the mobile signal at both BS_{2} and BS_{3} is within 3dB compared to BS_{1}. The numbers chosen here merely serve to illustrate the variability in received signal power (and its subsequent impact on TOA estimation accuracy), but do not attempt to model the specific thresholds and soft handover mechanisms used in CDMA standards.
Table 1 and Table 2 give the different β factors for the three cases of interest, and for different values for the pathloss and shadowing models, where a typical 50% correlation factor is assumed (a=b). For the numerical results, we assumed typical parameters for the radio propagation channel model. A twosegment pathloss model with breakpoint at distance d_{o}=200m and exponents n=2 and 4, respectively, was used, with a load of 20 users per cell, and a cell radius of 2Km. The relative powerfactors given in Table 1 and Table 2 clearly show that large variations occur across the base stations depending on the mobile soft handoff conditions and its proximity to one or more base station, and this will impact the accuracy of signal estimation at the different base stations as will be discussed next.
β_{1}
β_{2}
β_{3}
β_{4}
β_{5}
β_{6}
β_{7}
Case 1 (mobile close to home BS)
1
0.0216
0.0113
0.0069
0.0045
0.0031
0.0021
Case 2 (mobile in twoway SHO)
1
0.6982
0.2215
0.1202
0.0735
0.0485
0.0331
Case 3 (mobile in threeway SHO)
1
0.7922
0.6353
0.2993
0.1701
0.1065
0.0706
Table 1.
Relative received power factors for various mobile softhandoff link conditions. Shadowing s.t.d σ=8dB.
β1
β2
β3
β4
β5
β6
β7
Case 1 (mobile close to home BS)
1
0.0248
0.0125
0.0072
0.0045
0.0030
0.0020
Case 2 (mobile in twoway SHO)
1
0.7000
0.2281
0.1258
0.0761
0.0486
0.0322
Case 3 (mobile in threeway SHO)
1
0.7985
0.6443
0.3403
0.1953
0.1252
0.0808
Table 2.
Relative received power factors for various mobile softhandoff link conditions. Shadowing s.t.d σ=12dB.
Timing synchronization for CDMA signals is typically implemented in two steps consisting of an initial coarse timing acquisition (within an uncertainty range on the order of onechip interval), followed by fine time tracking achieved by a delaylocked loop (DLL) mechanism (Viterbi, 1995). In this work, we assume that the initial timing acquisition has been achieved, and focus on the more accurate DLL tracking as the main signal timing estimation mechanism. For each base station involved in mobile positioning, a DLL device continuously attempts to bring the local code timing estimate in perfect alignment with the incoming mobile signal. However, this timing estimation will be subject to error due to noise, fading and multipleaccess interference. In the following, we consider TOA estimation based on a noncoherent DLL (NCDLL) scheme. In practice, the NCDLL is preferred over a coherent structure because of its insensitivity to carrier phase error and data modulation. Figure 4 depicts a block diagram showing the different processing stages of a NCDLL code tracking loop. Because QPSK spreading is used, the NCDLL employs both I & Q branches where the I & Q received signals, after down conversion and chip matched filtering, are fed to two early & late branches which correlate the spreadspectrum waveforms with advanced and delayed code sequence replicas. The outputs obtained at the I & Q channels of the early & late branches can be obtained as (Viterbi, 1995)
YI±=ANPTCR(τ±δ)cosθE8
YQ±=ANPTCR(τ±δ)cosθE9
where A is the fading signal envelope modeled as a Rayleigh random variable, and θ is its uniform phase. δ is the early & late timing offset (typically set to T_{c}/2), P is the signal power, and N is the number of accumulated chips. The function R(τ) is a correlation function given by the convolution of the impulse responses of the pulseshaping filter and its matched filter (Viterbi, 1995):
R(τ)=h(τ)∗h(−τ)=∫−∞∞H(f)2cos(2πfcτ)dfE10
The DLL I & Q correlator outputs are combined as shown in Figure 4, and a resulting discriminator metric ΔZ is lowpass filtered to form an error signal used to drive a numericallycontrolled oscillator (NCO) that controls the code timing adjustment. The DLL discriminator output can be obtained as (Viterbi, 1995):
ΔZ=βAS(ε)+ηE11
where β is the power reduction factor reflecting the impact of DLL operation at different base stations (compared to the main serving base station, as discussed previously). The term η represents the combined Gaussian noise term with variance I_{0}. The function S(ε) is known as the normalized DLL Scurve, and is given by
S(ε)=R2(ε−δ)−R2(ε+δ)E12
where ε is the normalized timing error given by
ε=(τ−τ^)/TcE13
with τ denoting the correct TOA and τ^ the estimated one.
A common figure of merit for assessing the DLL performance is based on the tracking jitter variance. For additive white Gaussian (AWGN) channels, a simple upperbound approximation for the tracking jitter variance, valid for linearized firstorder DLL models, is derived in (Viterbi, 1995) as
σε≈2Io2+4NEcIoN2Ec2κ2E14
where E_{c} is the received chip energy, κ is the slope of the Scurve at the origin, and I_{0} is the variance of the thermal noise & MAI terms. However, there is no simple equivalent result for the case of frequencyselective multipath fading channels considered in this work, and the above approximation is only valid for high mobile received SignaltoNoise Ratio (SNR) conditions, which is not necessarily the case at the remote base stations (other than the home serving one), as discussed previously. In this work, we resort to a more accurate performance analysis valid for all SNR conditions and based on the computation of the full statistics (probability density function) of the DLL tracking loop, as discussed next.
4.2. DLL Tracking Error Statistics
For the purpose of analyzing the impact of TOA estimation error on mobile positioning accuracy, we need to obtain accurate statistics for the TOA measurements at each of the base stations involved in tracking the given mobile TOA. Hypothetical distributions (Gaussian models) are commonly assumed for the TOA timing estimation error. The same variance is usually assigned to the TOA error at all base stations involved in mobile positioning, and this variance is sometimes set rather arbitrarily. In this work, we use results based on a more rigorous analysis with full derivation of the TOA error probability density function (PDF). The results are obtained following the approach presented in (Su & Yen, 1997) and extended to quadraturespread CDMA signals in (Khan, 2009). Assuming a discretetime model, the analysis is based on a discretetime Markov model for the residual DLL error, according to the following equation:
εk=εk−1−KNCO[βAk−12S(εk−1)+ηk−1]E15
where K_{NCO} is the NCO gain, and η_{k} is the additive Gaussian noise term. It can be seen that the residual tracking error follows a discretetime Markov process for which the first order probability distribution can be obtained using the KolmogorovChapman equation (Su & Yen, 1997):
pk(ε)=∫−∞∞fk−1(εx)pk−1(xε0)dxE16
where ԑ_{0} is the initial timing error, p_{k1}(xԑ_{0}) is the PDF of x given ԑ_{0}, and f_{k1}(ԑx) is the transition pdf of ԑ_{k} given x. Through detailed analysis, one can calculate the exact expressions f_{k1}(ԑx) under assumption of a Rayleigh fading channel. Using numerical integration, it is then possible to iterate the KolmogorovChapman equation to get the PDF of the TOA estimation error. The lengthy details of this derivation are not included here, but can be found in (Khan, 2009).
The procedure outlined above can be done for different scenarios reflecting mobile TOA estimation at a given base station of interest. For our purpose, we analyze the performance of TOA DLL tracking at the different base stations for each one of the three cases described in Section 3. The results are shown in Figures 5 through 7, where it is clearly seen that the TOA residual error behavior can vary widely depending on the mobile position visàvis the tracking base station. For example, in Figure 5 which corresponds to a mobile very close to its home serving base station, the residual error at the home base station is wellconfined and nearly Gaussiandistributed (with a standard deviation found to be on the order of 0.15T_{c}), whereas for the other two base stations, the timing errors remain nearly uniformlydistributed with a standard deviation of 0.29T_{c} (which shows that the DLL loops at those base stations are not effectively tracking the mobile signal TOA). On the other hand, in Figure 6 and Figure 7 corresponding to 2way and 3way SHO, respectively, the mobile time tracking performance at BS2 and BS3 is markedly better, with distributions approaching that of the home BS1. It is to be noted that the number of users per cell (assumed the same for all cells, for simplicity) can also have a major impact on the timing estimation accuracy, regardless of the mobile position scenario (i.e., for all different cases discussed previously). Indeed, as shown in Figure 8, when the number of users increases, the tracking error PDFs are more widespread, and have an increasing error s.t.d (given with a normalization factor of 1/T_{c}).
In timeofarrival radiolocation techniques, the distance is calculated as the propagation time multiplied by the speed of light c. Lineofsight (LOS) propagation is assumed, whereby the mobile signal travels on a direct path at the constant speed of light in free space. It is further assumed that the base stations are synchronized and the mobile transmission time is known (set to zero for simplicity). The TOA measurements, produced at each base station, are therefore directly proportional to the mobilebase distance separation. Geometrically, circles centered at the base stations can be drawn with the calculated distance as the radius. With the help of three base stations, the mobile location can be found geometrically as the intersection of the corresponding circles. However, in the presence of noise and interference, the three circles may not intersect at a single point. Therefore, the geometric approach is not suitable, and several other “statistical” techniques have been proposed (Caffrey, 1999), (Sayed et al., 2005) to process the noisy data. Many of these are based on iterative algorithms using leastsquares or gradient search minimization. On the other hand, since the positioning equations involved are typically nonlinear, some traditional approaches based on linearization followed by a gradient search were proposed (Niezgoda & Ho, 1994). However, these approaches suffer from sensitivity to initialization errors and convergence problems. More recently, some researchers introduced new closedform linear techniques (Caffrey & Stuber, 2000), (Chan & Ho, 1994), but the drawback of these methods is that optimum location estimates can only be found at high SNR values, which may not always be the case in practice.
Another interesting solution was proposed by (Chan et al., 2006) and shown to have nearoptimum performance, with the added advantage of reduced complexity. Their method uses an Approximate Maximum Likelihood (AML) algorithm, which we adopt in this study. In the following, we give a general overview of the AML algorithm processing steps, and introduce a modification tailored towards the case of unequal TOA variances at the different base stations, and relevant to mobile radiolocation under nearfar interference conditions (which is the focus of this work), as will be outlined subsequently. Based on the assumed cellular geometry, the true distances between the BSs and MS are given by
Ri=(x−xi)2+(y−yi)2E17
i=1,2,……,N_{BS}
where N_{BS} is the number of involved BSs. The measured distances, l_{i}, is given by
li=Ri+εiE18
i=1,2,……,N_{BS}
where ε_{i} is the DLL timing error. In matrix form, this is written as
l=R+εE19
Dividing by c to get the measured TOA vector T, we obtain
T=Rc+εc=T0+eE20
where
R=[R1…RN]t=R(Θ)E21
and
e=[e1....eN]tE22
is the vector of additive measurement noise, and T^{0} is the vector of true TOAs. The original AML algorithm proposed by (Chan et al., 2006) assumes that all BSs have an equal error variance. Hence, the elements of e are assumed to be independent, zeromean Gaussian random variables with covariance matrix
Q=E{eet}=σ2IE23
The AML algorithm can be generalized to account for different values of the error variance at different BSs. If the variance of the error can be estimated at each BS, one can use this information to improve the localization performance by giving more trust to BSs with lower error variance. In this case the elements of e are assumed to be independent, zero mean Gaussian random variables with covariance matrix:
Q=E{eet}=diag[σ12…σN2]E24
The conditional probability density function (PDF) of T given Θ is given by
The ML estimate of the MS position (x, y) is the Θ that minimizes J (Trees, 1968). Minimizing J is done by setting its gradient with respect to Θ to zero. Considering first the derivative with respect to the position variable x, we have
∂J∂x=∑i=1NBS1σi2∂ei2∂x=∑i=1NBS1σi22ei⋅∂ei∂xE27
By expressing the timing error in terms of the difference between the true TOA and measured one, this gives
The above steps starting with Equation (27) can be repeated for ∂J/∂y in a straightforward manner. Finally, by setting the gradient of J with respect to Θ to zero, we get the two ML equations
In a more compact form, Equation (34) can be rewritten as
AΘ = b
(38)
with the matrix A and vector b being functions of Θ. A suboptimal solution based on a linear model (Chan & Ho, 1994) can be used as a first initial estimate of Θ, which will in turn give starting values of A and b. Then, solving Equation (33) produces a new value of Θ to update A and b, and subsequently Θ. This iterative procedure first gives an approximate maximum likelihood (AML) estimator, which can then be iterated a number of times to obtain a final solution. The final solution takes the Θ that gives the smallest J in Equation (26). This ensures that the AML will not diverge. In fact, simulation results presented in (Chan et al., 2006) show that the AML can nearly achieve the Cramer–Rao lower bound (CRLB) with a small number of iterations (typically on the order of five updates are found to be sufficient). It should be noted that, for the special case of equal measurement variance assumption (as in the original AML), the common σ term drops from Equation (32). However, with the modified AML in Equation (33), the quantities g_{i} and h_{i} will be different for different BSs, and this is found to yield some improvement in performance as will be discussed next.
5.2. Positioning Accuracy with the Modified AML Algorithm
To illustrate the impact of the modified AML that takes into account the unequal error variances, we consider an example of three BSs with different measurement error statistics for mobiles located in three different regions, as described in the different cases of Section 3. Figure 9 shows the geometry of the layout used for this purpose, where distances are shown in meters. The dense scatter points represents a total of 10^{5} noisy mobile locations generated according to the TOA statistics for the different cases of interest which are referred to as MS1, MS2, MS3, corresponding to mobiles classified in Section 3 as Case 1, 2, and 3, respectively. For each mobile location, the AML and modified AML algorithms are executed to obtain the estimated mobile coordinates, and by comparing with the known mobile position, the resulting positioning error can be computed. The different parameters for the radio channel and relative received power factors used are based on the results given in Table 1 of Section 3.
Various aspects of mobile positioning performance are illustrated in Figure 10 based on the cumulative distribution function (CDF) of the radiolocation error. First, it can be seen that the modified AML which takes into account the differences in TOA noise statistics at the various base stations (as proposed in this work) outperforms, albeit slightly, the conventional AML which assumes equal TOA noise variance at all base stations. This is more pronounced for mobiles categorized as MS1 (i.e., which are in close proximity to their home base stations). The other major observation from Figure 10 is that a large difference in radiolocation accuracy is present depending on the mobile relative position with respect to the base stations. It is clearly seen that MS2 and particularly MS3 mobiles, corresponding to mobiles in 2way and 3way soft handoff, respectively, achieve much better performance as opposed to MS1. This is due to the poor TOA accuracy at BS2 and BS3 for the latter case, owing to the overwhelming nearfar interference experienced by the mobile at the remote bases stations, as discussed in Section 4. As an example, the probability of the residual positioning error being less than 20m is almost one for MS3 mobiles, while it is on the order of 70% for MS2, and only 50% for MS1–type mobiles. Similar observations also hold for other distances. Therefore, as highlighted throughout this study, the accuracy of mobile positioning is best when the mobile is in close proximity to a border cell region where soft handoff connectivity is established with one or two more cells involved in its positioning in addition to its home serving cell.
This study dealt with the performance analysis of timeofarrival (TOA) techniques for mobile positioning in CDMA wireless cellular networks. Since several base stations are typically needed for mobile radiolocation, the problem of weak signal hearability at remote base stations is a major challenge, and a detailed analysis of this issue was presented by taking into account the nearfar interference usually present in CDMA cellular networks.
TOAbased positioning methods are wellsuited for wide deployment of radiolocation services since the synchronization circuits that can extract timing information are readily available at the base stations receiving the mobile signal. For CDMA signals, the delaylocked loop (DLL) is commonly used as a TOA estimation device, and a detailed analysis of DLLbased TOA tracking was presented taking into account the cellular network layout, cell loading and otheruser interference, as well as RF channel shadowing and pathloss that affect received signal strength at the different base stations. In particular, it was shown that the TOA tracking error statistics can vary considerably depending on the mobile link conditions with respect to the base stations involved in its positioning, and soft handoff (SHO) links with two or three base stations for mobiles close to bordercell regions were found to improve the precision of mobile TOA tracking.
Using TOA information, a mobile radiolocation positioning algorithm based on a computationallyefficient, nearoptimum approximate maximum likelihood (AML) processing was presented, and a generalization to the case of unequal timing error variances at the different base stations was also derived. Numerical results were presented to quantify the achievable positioning accuracy of the AML algorithms, and it was also shown that SHO radio links with two or particularly three base stations have a clear impact on the precision of mobile radiolocation. Finally, it should be noted that the problem of nonlineofsight (NLOS) propagation, which constitutes another major challenge to mobile positioning accuracy, was not included in this study, and would be addressed in future work.
This work was gratefully supported by King Abdulaziz City for Science & Technology and King Fahd University of Petroleum & Minerals, Saudi Arabia.
References
1.CaffreyJ. J.1999 Wireless Location in CDMA Cellular Radio Systems. Kluwer Academic Publishers.
2.CaffreyJ. J.StuberG. L.2000 Effects of multiple access interference on the noncoherent delay lock loop.
3.ChanTY.HangH.ChingCP.2006Exact and approximate maximum likelihood localization algorithms.
4.ChanY. T.HoK. C.1994A simple and efficient estimator for hyperbolic location.
5.ChanY. T.YauC. H.ChingP. C.2006Exact and approximate maximum likelihood localization algorithms.
6.DahlmanE.GudmundsonB.NilssonM.SkoldA.2000 UMTS/IMT2000 based on wideband CDMA. ber 2000) 70 80.
7.GoshA.LoveR.1998 Mobile station location in a DSCDMA system,
8.KaplanE.1996Understanding GPS: Principles and Applications. Norwood, MA: Artech House.
9.KhanH. R.2009 DLL Code Tracking for CDMA Signals under Fading and Multiple Access Interference, M.S. Thesis, King Fahd University of Petroleum & Minerals, February 2009.
10.NiezgodaG. H.HoK. C.1994 Geolocation by combined range difference and range rate measurements,
11.RappaportT.ReedJ.WornerB.1996Position location using wireless communications on highways of the future. IEEE Communications Magazine
, (October 1996) 33 41 .
13.SuS. L.YenN. Y.1997Performance of digital code tracking loops for directsequence spreadspectrum signals in mobile radio channels. IEEE Transactions on Communications
, 45 (May 1997) 596 604 .
14.TreesH. L. V.1968Detection, Estimation, and Modulation Theory: Part 1, Wiley, New York.
15.ViterbiA. J.1995CDMA: Principles of Spread Spectrum Communication, AddisonWesley.
Written By
M. A.Landolsi, A. H. Muqaibel, A. S. AlAhmari, H.R. Khan and R. A. AlNimnim