# A Hybrid RSS/TOA Method for 3D Positioning in an Indoor Environment

Authors: Smita Tiwari; Donglin Wang; Michel Fattouche; Fadhel Ghannouchi
Journal: ISRN Signal Processing
 (2012)

---

## Abstract

This paper investigates 3D positioning in an indoor line of sight (LOS) and nonline of sight (NLOS) combined environment. It is a known fact that time-of-arrival-(TOA-) based positioning outperforms other techniques in LOS environments; however, multipath in an indoor environment, especially NLOS multipath, significantly decreases the accuracy of TOA positioning. On the other hand, received-signal-strength-(RSS-) based positioning is not affected so much by NLOS multipath as long as the propagation attenuation can be correctly estimated and the multipath effects have been compensated for. Based on this fact, a hybrid weighted least square (HWLS) RSS/TOA method is proposed for target positioning in an indoor LOS/NLOS environment. The identification of LOS/NLOS path is implemented by using Nakagami distribution. An experiment is conducted in the iRadio lab, in the ICT building at the University of Calgary, in order to (i) demonstrate the availability of Nakagami distribution for the identification of LOS and NLOS path, (ii) estimate the pass loss exponent for RSS technique, and (iii) verify our proposed scheme.


---

## Body

## 1. Introduction

In the last few years, the interest in indoor location has significantly increased due to the importance of the 3D indoor localization with a high level of accuracy. Wireless location estimation [1] is usually based on measuring one or more radio signal propagation attributes. Different attributes of a received radio signal such as time of arrival (TOA), time difference of arrival (TDOA), and angle of arrival (AOA) have been widely applied in location estimation of a target, such as a mobile station (MS) [2–6]. Unfortunately, all of the aforementioned techniques require a line of sight (LOS) link between MS and a number of base stations (BSs) of known positions. In a nonline of sight (NLOS) environment, their performances can degrade significantly.In most indoor environments, NLOS, one of the dominant factors significantly affecting positioning accuracy, is inevitable. A variety of existing tests can identify whether a measurement fits in an LOS or a NLOS profile. Historical range measurements have been used for NLOS identification purposes in [7–9]. By applying a range scaling algorithm, the true range can be estimated by constrained minimization [7]. In [8], NLOS identification is implemented via a running variance of range estimates. Experimental results conducted in [9] indicate that the standard deviation of the range measurements in an NLOS environment is much larger than its LOS counterpart, which can also be used for identification. Another category in LOS/NLOS identification is by examining the distribution of the received signal. In a LOS channel, the TOA measurements are typically Gaussian distributed, whereas in a NLOS channel, they can be exponentially distributed [10]. Accordingly, LOS propagation can be identified by testing whether the TOA measurements have a Gaussian distribution. LOS/NLOS is identified in [11] by performing a likelihood test with prior knowledge such as kurtosis, root mean square (RMS) delay, and excess delay. The distribution function of the received power envelope is used in [12] to identify LOS/NLOS. The authors in [13] take advantage of the fact that the received power envelope in a LOS environment has a Ricean distribution, while it has a Rayleigh distribution in a NLOS environment. Using average fading duration and the level crossing rate of the power envelopes is another way to deal with LOS/NLOS identification [14]. A Nakagami-m distribution is proposed in [15] for LOS/NLOS identification, which is proved to be the optimal radio channel model in an indoor environment.When identifying a NLOS channel, it is necessary to mitigate its effects. A method is proposed in [16] to correct NLOS range error based on the prior knowledge of the range statistics. In [17], the authors try to weight the NLOS measurements by the residue of the measured range. A hard weight selection in weighted least square (WLS) is introduced for a NLOS correction [11]. In this paper, a hybrid weighted least square (HWLS) RSS/TOA method is proposed for target positioning in an indoor LOS/NLOS environment based on the fact that (i) TOA-based positioning outperforms other techniques in LOS environments; however, multipath in an indoor environment, especially NLOS multipath, significantly decreases its accuracy. (ii) On the other hand, the received-signal-strength-(RSS-) based positioning is not affected as much by NLOS multipath as long as the propagation attenuation can be correctly estimated and the multipath effects have been compensated for. So, TOA can be applied for LOS cases while RSS can be used for NLOS cases. By conducting an indoor experiment, three objectives have been reached in this paper: (i) to demonstrate the availability of Nakagami distribution for the identification of LOS and NLOS, (ii) to estimate the pass loss exponent based on RSS data, and (iii) to verify the proposed HWLS RSS/TOA.The paper is organized as follows. Section2 proposes the novel 3D positioning algorithm, HWLS RSS/TOA, for LOS/NLOS environments, and derives its Cramer-Rao lower bound (CRLB). Section 3 describes our experimental setup. In Section 4, the identification of LOS/NLOS is implemented, the path loss model is estimated, and the experimental results are reported, followed by a conclusion in Section 5.


## 2. Proposed 3D Positioning Algorithm for LOS/NLOS Environments

### 2.1. Positioning Scheme in an Indoor Environment

It is a known fact that TOA-based positioning [18] outperforms other techniques such as RSS-based positioning in a LOS environment. However, multipath in a dense urban area or in an indoor environment, especially NLOS multipath, significantly decreases its accuracy. On the other hand, RSS is not affected as much by NLOS multipath as long as the propagation attenuation can be correctly estimated and the multipath effects have been compensated for. Based on these facts, we are proposing a hybrid weighted least square (HWLS) RSS/TOA method for positioning in LOS/NLOS environments, which contains the following procedure by assuming N BSs (Txs) of known positions in a wireless location network which attempts to estimate the position of MS (Rx).(a)
Initialize by settingi=1.(b)
Consider the received signal from theith BS and identify whether the direct path is LOS or NLOS using Nakagami-m distribution. If it is LOS, go to (c); otherwise, jump to (d);(c)
Use a WLS TOA for LOS range estimation. Ifi=N, jump to (f); otherwise, let i=i+1 and go back to (b).(d)
Estimate the propagation exponentα as well as its shadowing effect, thus obtaining a path loss model.(e)
Based on the WLS RSS proposed in Subsection2.2 below, estimate the NLOS range. If i=N, go to (f); otherwise, let i=i+1 and go back to (b).(f)
Implement target positioning based on allN range measurements.


### 2.2. Proposed 3D WLS

#### 2.2.1. RSS

Letxi=[xi,yi,zi]T, i=1,2,…,N, be the known coordinates for the ith BS and let x=[x,y,z]T be the unknown coordinates of MS. The estimated range by means of RSS, R̂iRSS, between the ith BS and MS, can be modeled as(1)R̂iRSS=FiPitPir+ωi,
where Pir is the received power, Pit is the transmitted power, Fi represents an unknown but deterministic coefficient, and ωi denotes the range error. The received power Pir can be further expressed as(2)Pir=FiPitdiα,
where α is the propagation attenuation, and di is the range between the ith BS and MS. Substituting (2) into (1), we have(3)R̂iRSS=[(x-xi)2+(y-yi)2+(z-zi)2]α/2+ωi.
First, let us consider a noise-free environment, where (3) can be expressed as(4)R̂iRSS=[(x-xi)2+(y-yi)2+(z-zi)2]α/2.
Or equivalently,(5)R̂iRSS2/α=x2+xi2-2xxi+y2+yi2-2yyi+z2+zi2-2zzi.
Furthermore,(6)xxi+yyi+zzi-0.5Q=12(xi2+yi2+zi2-R̂iRSS2/α),
where Q is defined as(7)Q=x2+y2+z2.
In matrix notation, (6) can be expressed as(8)GΘ=H,
where(9)G=Δ[x1y1z1-0.5⋮⋮⋮⋮xNyNzN-0.5],H=Δ12[x12+y12+z12-R̂1RSS2/α⋮xN2+yN2+zN2-R̂NRSS2/α],andtheunknownvariableΘ=Δ[xyzQ].
Now considering a noisy environment as in (3), (8) is replaced by(10)ω=GΘ-H,
where the range error vector is ω=Δ[ω1,ω2,…,ωN]T. By minimizing the range error ω, the least square (LS) solution gives an estimate of Θ as(11)Θ̂LSRSS=argmin⁡Θ(GΘ-H)T(GΘ-H)=(GTG)-1GTH.
To enhance its estimation accuracy, a weighted matrix W is used in (11) so that Θ̂LSRSS satisfies the fundamental constraint given in (7), and (11) is replaced by(12)Θ̂WLSRSS=argmin⁡Θ(GΘ-H)TWRSS(GΘ-H)=(GTWRSSG)-1GTWRSSH,
where WRSS=SJS,S=Δdiag[d1RSS2-α/α,d2RSS2-α/α,…,dNRSS2-α/α],J=Δdiag[σ1RSS2,σ2RSS2,…,σNRSS2], and σiRSS2, i=1,2,…,N, is the estimation variance using RSS.


#### 2.2.2. TOA

The 2D LS TOA is derived in [19]. Appendix 5 extends it to the 3D WLS TOA case. In Appendix 5, WLS TOA is used instead of LS TOA in order to enhance its estimation accuracy in a LOS environment.


### 2.3. CRLB for the Proposed Scheme

The CRLB for 2D TOA positioning is derived in [19]. By using a similar derivation, the CRLB for 3D TOA positioning can be obtained as(13a)CRLBWLSTOA(x̂)={E[∑i=1N(x-xi)22σi2di2∑i=1N(x-xi)(y-yi)2σi2di2∑i=1N(x-xi)(z-zi)2σi2di2∑i=1N(x-xi)(y-yi)2σi2di2∑i=1N(y-yi)22σi2di2∑i=1N(y-yi)(z-zi)2σi2di2∑i=1N(x-xi)(z-zi)2σi2di2∑i=1N(y-yi)(z-zi)2σi2di2∑i=1N(z-zi)22σi2di2]}-1,and for each dimension, we have(13b)CRLBWLSTOA(x)=[CRLBWLSTOA(x̂)]11,(13c)CRLBWLSTOA(y)=[CRLBWLSTOA(x̂)]22,(13d)CRLBWLSTOA(z)=[CRLBWLSTOA(x̂)]33,
while the CRLB for 3D RSS positioning can be obtained as(14)CRLBWLSRSS(x̂)={E[∑i=1N(x-xi)22σi2di2(2-α)∑i=1N(x-xi)(y-yi)2σi2di2(2-α)∑i=1N(x-xi)(z-zi)2σi2di2(2-α)∑i=1N(x-xi)(y-yi)2σi2di2(2-α)∑i=1N(y-yi)22σi2di2(2-α)∑i=1N(y-yi)(z-zi)2σi2di2(2-α)∑i=1N(x-xi)(z-zi)2σi2di2(2-α)∑i=1N(y-yi)(z-zi)2σi2di2(2-α)∑i=1N(z-zi)22σi2di2(2-α)]}-1.Assume in our indoor experiment that there are M LOS BSs out of N BSs, for example, the ith BS, 1≤i≤M, corresponds to a LOS link between itself and MS, while the jth BS, M+1≤j≤N, corresponds to a NLOS link. In other words, our proposed method has M TOA measurements and N-M RSS measurements. Therefore, using a similar derivation as [19], the CRLB for the proposed HWLS RSS/TOA can be expressed as(15a)CRLBWLSRSS(x̂)={E[∂2ln⁡(p(r∣x))∂x2∂2ln⁡(p(r∣x))∂x∂y∂2ln⁡(p(r∣x))∂x∂z∂2ln⁡(p(r∣x))∂x∂y∂ln⁡(p(r∣x))∂y2∂2ln⁡(p(r∣x))∂y∂z∂2ln⁡(p(r∣x))∂x∂z∂2ln⁡(p(r∣x))∂y∂z∂ln⁡(p(r∣x))∂z2]}-1,where(15b)∂2ln⁡(p(r∣x))∂x2=∑i=1M(x-xi)22σi2di2+∑i=M+1N(x-xi)22σi2di2(2-α),(15c)∂2ln⁡(p(r∣x))∂y2=∑i=1M(y-yi)22σi2di2+∑i=M+1N(y-yi)22σi2di2(2-α),(15d)∂2ln⁡(p(r∣x))∂z2=∑i=1M(z-zi)22σi2di2+∑i=M+1N(z-zi)22σi2di2(2-α),(15e)∂2ln⁡(p(r∣x))∂x∂y=∑i=1M(x-xi)(y-yi)2σi2di2+∑i=M+1N(x-xi)(y-yi)2σi2di2(2-α),(15f)∂2ln⁡(p(r∣x))∂x∂z=∑i=1M(x-xi)(z-zi)2σi2di2+∑i=M+1N(x-xi)(z-zi)2σi2di2(2-α),(15g)∂2ln⁡(p(r∣x))∂y∂z=∑i=1M(y-yi)(z-zi)2σi2di2+∑i=M+1N(y-yi)(z-zi)2σi2di2(2-α).


## 2.1. Positioning Scheme in an Indoor Environment

It is a known fact that TOA-based positioning [18] outperforms other techniques such as RSS-based positioning in a LOS environment. However, multipath in a dense urban area or in an indoor environment, especially NLOS multipath, significantly decreases its accuracy. On the other hand, RSS is not affected as much by NLOS multipath as long as the propagation attenuation can be correctly estimated and the multipath effects have been compensated for. Based on these facts, we are proposing a hybrid weighted least square (HWLS) RSS/TOA method for positioning in LOS/NLOS environments, which contains the following procedure by assuming N BSs (Txs) of known positions in a wireless location network which attempts to estimate the position of MS (Rx).(a)
Initialize by settingi=1.(b)
Consider the received signal from theith BS and identify whether the direct path is LOS or NLOS using Nakagami-m distribution. If it is LOS, go to (c); otherwise, jump to (d);(c)
Use a WLS TOA for LOS range estimation. Ifi=N, jump to (f); otherwise, let i=i+1 and go back to (b).(d)
Estimate the propagation exponentα as well as its shadowing effect, thus obtaining a path loss model.(e)
Based on the WLS RSS proposed in Subsection2.2 below, estimate the NLOS range. If i=N, go to (f); otherwise, let i=i+1 and go back to (b).(f)
Implement target positioning based on allN range measurements.


## 2.2. Proposed 3D WLS

### 2.2.1. RSS

Letxi=[xi,yi,zi]T, i=1,2,…,N, be the known coordinates for the ith BS and let x=[x,y,z]T be the unknown coordinates of MS. The estimated range by means of RSS, R̂iRSS, between the ith BS and MS, can be modeled as(1)R̂iRSS=FiPitPir+ωi,
where Pir is the received power, Pit is the transmitted power, Fi represents an unknown but deterministic coefficient, and ωi denotes the range error. The received power Pir can be further expressed as(2)Pir=FiPitdiα,
where α is the propagation attenuation, and di is the range between the ith BS and MS. Substituting (2) into (1), we have(3)R̂iRSS=[(x-xi)2+(y-yi)2+(z-zi)2]α/2+ωi.
First, let us consider a noise-free environment, where (3) can be expressed as(4)R̂iRSS=[(x-xi)2+(y-yi)2+(z-zi)2]α/2.
Or equivalently,(5)R̂iRSS2/α=x2+xi2-2xxi+y2+yi2-2yyi+z2+zi2-2zzi.
Furthermore,(6)xxi+yyi+zzi-0.5Q=12(xi2+yi2+zi2-R̂iRSS2/α),
where Q is defined as(7)Q=x2+y2+z2.
In matrix notation, (6) can be expressed as(8)GΘ=H,
where(9)G=Δ[x1y1z1-0.5⋮⋮⋮⋮xNyNzN-0.5],H=Δ12[x12+y12+z12-R̂1RSS2/α⋮xN2+yN2+zN2-R̂NRSS2/α],andtheunknownvariableΘ=Δ[xyzQ].
Now considering a noisy environment as in (3), (8) is replaced by(10)ω=GΘ-H,
where the range error vector is ω=Δ[ω1,ω2,…,ωN]T. By minimizing the range error ω, the least square (LS) solution gives an estimate of Θ as(11)Θ̂LSRSS=argmin⁡Θ(GΘ-H)T(GΘ-H)=(GTG)-1GTH.
To enhance its estimation accuracy, a weighted matrix W is used in (11) so that Θ̂LSRSS satisfies the fundamental constraint given in (7), and (11) is replaced by(12)Θ̂WLSRSS=argmin⁡Θ(GΘ-H)TWRSS(GΘ-H)=(GTWRSSG)-1GTWRSSH,
where WRSS=SJS,S=Δdiag[d1RSS2-α/α,d2RSS2-α/α,…,dNRSS2-α/α],J=Δdiag[σ1RSS2,σ2RSS2,…,σNRSS2], and σiRSS2, i=1,2,…,N, is the estimation variance using RSS.


### 2.2.2. TOA

The 2D LS TOA is derived in [19]. Appendix 5 extends it to the 3D WLS TOA case. In Appendix 5, WLS TOA is used instead of LS TOA in order to enhance its estimation accuracy in a LOS environment.


## 2.2.1. RSS

Letxi=[xi,yi,zi]T, i=1,2,…,N, be the known coordinates for the ith BS and let x=[x,y,z]T be the unknown coordinates of MS. The estimated range by means of RSS, R̂iRSS, between the ith BS and MS, can be modeled as(1)R̂iRSS=FiPitPir+ωi,
where Pir is the received power, Pit is the transmitted power, Fi represents an unknown but deterministic coefficient, and ωi denotes the range error. The received power Pir can be further expressed as(2)Pir=FiPitdiα,
where α is the propagation attenuation, and di is the range between the ith BS and MS. Substituting (2) into (1), we have(3)R̂iRSS=[(x-xi)2+(y-yi)2+(z-zi)2]α/2+ωi.
First, let us consider a noise-free environment, where (3) can be expressed as(4)R̂iRSS=[(x-xi)2+(y-yi)2+(z-zi)2]α/2.
Or equivalently,(5)R̂iRSS2/α=x2+xi2-2xxi+y2+yi2-2yyi+z2+zi2-2zzi.
Furthermore,(6)xxi+yyi+zzi-0.5Q=12(xi2+yi2+zi2-R̂iRSS2/α),
where Q is defined as(7)Q=x2+y2+z2.
In matrix notation, (6) can be expressed as(8)GΘ=H,
where(9)G=Δ[x1y1z1-0.5⋮⋮⋮⋮xNyNzN-0.5],H=Δ12[x12+y12+z12-R̂1RSS2/α⋮xN2+yN2+zN2-R̂NRSS2/α],andtheunknownvariableΘ=Δ[xyzQ].
Now considering a noisy environment as in (3), (8) is replaced by(10)ω=GΘ-H,
where the range error vector is ω=Δ[ω1,ω2,…,ωN]T. By minimizing the range error ω, the least square (LS) solution gives an estimate of Θ as(11)Θ̂LSRSS=argmin⁡Θ(GΘ-H)T(GΘ-H)=(GTG)-1GTH.
To enhance its estimation accuracy, a weighted matrix W is used in (11) so that Θ̂LSRSS satisfies the fundamental constraint given in (7), and (11) is replaced by(12)Θ̂WLSRSS=argmin⁡Θ(GΘ-H)TWRSS(GΘ-H)=(GTWRSSG)-1GTWRSSH,
where WRSS=SJS,S=Δdiag[d1RSS2-α/α,d2RSS2-α/α,…,dNRSS2-α/α],J=Δdiag[σ1RSS2,σ2RSS2,…,σNRSS2], and σiRSS2, i=1,2,…,N, is the estimation variance using RSS.


## 2.2.2. TOA

The 2D LS TOA is derived in [19]. Appendix 5 extends it to the 3D WLS TOA case. In Appendix 5, WLS TOA is used instead of LS TOA in order to enhance its estimation accuracy in a LOS environment.


## 2.3. CRLB for the Proposed Scheme

The CRLB for 2D TOA positioning is derived in [19]. By using a similar derivation, the CRLB for 3D TOA positioning can be obtained as(13a)CRLBWLSTOA(x̂)={E[∑i=1N(x-xi)22σi2di2∑i=1N(x-xi)(y-yi)2σi2di2∑i=1N(x-xi)(z-zi)2σi2di2∑i=1N(x-xi)(y-yi)2σi2di2∑i=1N(y-yi)22σi2di2∑i=1N(y-yi)(z-zi)2σi2di2∑i=1N(x-xi)(z-zi)2σi2di2∑i=1N(y-yi)(z-zi)2σi2di2∑i=1N(z-zi)22σi2di2]}-1,and for each dimension, we have(13b)CRLBWLSTOA(x)=[CRLBWLSTOA(x̂)]11,(13c)CRLBWLSTOA(y)=[CRLBWLSTOA(x̂)]22,(13d)CRLBWLSTOA(z)=[CRLBWLSTOA(x̂)]33,
while the CRLB for 3D RSS positioning can be obtained as(14)CRLBWLSRSS(x̂)={E[∑i=1N(x-xi)22σi2di2(2-α)∑i=1N(x-xi)(y-yi)2σi2di2(2-α)∑i=1N(x-xi)(z-zi)2σi2di2(2-α)∑i=1N(x-xi)(y-yi)2σi2di2(2-α)∑i=1N(y-yi)22σi2di2(2-α)∑i=1N(y-yi)(z-zi)2σi2di2(2-α)∑i=1N(x-xi)(z-zi)2σi2di2(2-α)∑i=1N(y-yi)(z-zi)2σi2di2(2-α)∑i=1N(z-zi)22σi2di2(2-α)]}-1.Assume in our indoor experiment that there are M LOS BSs out of N BSs, for example, the ith BS, 1≤i≤M, corresponds to a LOS link between itself and MS, while the jth BS, M+1≤j≤N, corresponds to a NLOS link. In other words, our proposed method has M TOA measurements and N-M RSS measurements. Therefore, using a similar derivation as [19], the CRLB for the proposed HWLS RSS/TOA can be expressed as(15a)CRLBWLSRSS(x̂)={E[∂2ln⁡(p(r∣x))∂x2∂2ln⁡(p(r∣x))∂x∂y∂2ln⁡(p(r∣x))∂x∂z∂2ln⁡(p(r∣x))∂x∂y∂ln⁡(p(r∣x))∂y2∂2ln⁡(p(r∣x))∂y∂z∂2ln⁡(p(r∣x))∂x∂z∂2ln⁡(p(r∣x))∂y∂z∂ln⁡(p(r∣x))∂z2]}-1,where(15b)∂2ln⁡(p(r∣x))∂x2=∑i=1M(x-xi)22σi2di2+∑i=M+1N(x-xi)22σi2di2(2-α),(15c)∂2ln⁡(p(r∣x))∂y2=∑i=1M(y-yi)22σi2di2+∑i=M+1N(y-yi)22σi2di2(2-α),(15d)∂2ln⁡(p(r∣x))∂z2=∑i=1M(z-zi)22σi2di2+∑i=M+1N(z-zi)22σi2di2(2-α),(15e)∂2ln⁡(p(r∣x))∂x∂y=∑i=1M(x-xi)(y-yi)2σi2di2+∑i=M+1N(x-xi)(y-yi)2σi2di2(2-α),(15f)∂2ln⁡(p(r∣x))∂x∂z=∑i=1M(x-xi)(z-zi)2σi2di2+∑i=M+1N(x-xi)(z-zi)2σi2di2(2-α),(15g)∂2ln⁡(p(r∣x))∂y∂z=∑i=1M(y-yi)(z-zi)2σi2di2+∑i=M+1N(y-yi)(z-zi)2σi2di2(2-α).


## 3. Experimental Site and Setup

### 3.1. Experimental Site

To validate the proposed scheme in a realistic environment, an indoor experiment is conducted in the iRadio Lab, in the ICT building, at the University of Calgary. The 3D experimental site is shown in Figure1. Its size is approximately equal to 18.0 m × 7.0 m × 3.5 m. The Lab is divided into three sections: (1) Section 1 is a small room, that is, 3.0 m × 4.0 m, located at the entrance of the laboratory, (2) Section 2 is its main body, and (3) Section 3 is a small conference room partitioned from the lab’s main body.Figure 1
Experiment site in the iRadio lab, ICT Building at University of Calgary.The outer wall of the lab is glass-walled, which is a big source for signal reflecting. There are a couple of metallic heating pipes located around the room ceiling. The lab is equipped with many electronics instruments, computers, WiFi devices, scanners, printers, furniture, and metal cabinets. Polished furniture can generate signal reflection. For the LOS scenario, all BSs and MS are placed in the 2nd part of the lab, while for the NLOS scenario, the NLOS BSs are placed in the 1st Section and the 3rd Section. In Section3, the receiver is kept where a NLOS situation is guaranteed. In Figure 1, all BSs and MS are marked with red circles.


### 3.2. Experimental Setup

The whole system architecture of this experiment is depicted in Figure2. Each BS is specified by using the Agilent vector signal generator (VSG) N5182A capable of generating probing signal up to 4 GHz. The output parameter of the VSG is selected based on the technical specification of the analog front-end board MAX2830 Evaluation Kit from Maxim, and an agile digital signal processing (DSP) unit. The overall technical specification used in this study is summarized in Table 1.Table 1
Technical specification for analog front-end board from maxim (max 2830 evaluation kit) used in experiment.
SpecificationTechnical detailAnalog front end band2.4 GHz–2.5 GHzRx noise figure3.3 dBRx sensitivity−75 dBmRxI/Q gain0.1 dBRx phase imbalance0.35Based band gain control range33–62 dBFigure 2
Test bed for 3D indoor location.MATLAB/ADS is used to generate an OFDM signal, which is then upload to the VSG at 2.45 GHz carrier frequency, a 5 MHz bandwidth, and an 89.6 MHz sampling frequency. The VSG is connected to an RF switch with 5 outputs, each of which is further connected to a high gain omnidirection antenna using a coaxis cable. All BSs are clock-synchronized. The receiving antenna is directly connected to the Maxim board. The front-end is composed of a variable gain low-noise amplifier (LNA) followed by anI/Q zero-IF down converter, a programmable lowpass filter (LPF), and a variable gain amplifier (VGA) as shown in Figure 2. After downconversion, the baseband signal is sampled and quantified by an analog-to-digital converter. A MATLAB program is then used to capture both the I component and the Q component from the front end. The estimated time delay is obtained at the output of a matched filter, where the known time-multiplexed at the RF switch and the extra cable delay are considered for time delay compensation.


## 3.1. Experimental Site

To validate the proposed scheme in a realistic environment, an indoor experiment is conducted in the iRadio Lab, in the ICT building, at the University of Calgary. The 3D experimental site is shown in Figure1. Its size is approximately equal to 18.0 m × 7.0 m × 3.5 m. The Lab is divided into three sections: (1) Section 1 is a small room, that is, 3.0 m × 4.0 m, located at the entrance of the laboratory, (2) Section 2 is its main body, and (3) Section 3 is a small conference room partitioned from the lab’s main body.Figure 1
Experiment site in the iRadio lab, ICT Building at University of Calgary.The outer wall of the lab is glass-walled, which is a big source for signal reflecting. There are a couple of metallic heating pipes located around the room ceiling. The lab is equipped with many electronics instruments, computers, WiFi devices, scanners, printers, furniture, and metal cabinets. Polished furniture can generate signal reflection. For the LOS scenario, all BSs and MS are placed in the 2nd part of the lab, while for the NLOS scenario, the NLOS BSs are placed in the 1st Section and the 3rd Section. In Section3, the receiver is kept where a NLOS situation is guaranteed. In Figure 1, all BSs and MS are marked with red circles.


## 3.2. Experimental Setup

The whole system architecture of this experiment is depicted in Figure2. Each BS is specified by using the Agilent vector signal generator (VSG) N5182A capable of generating probing signal up to 4 GHz. The output parameter of the VSG is selected based on the technical specification of the analog front-end board MAX2830 Evaluation Kit from Maxim, and an agile digital signal processing (DSP) unit. The overall technical specification used in this study is summarized in Table 1.Table 1
Technical specification for analog front-end board from maxim (max 2830 evaluation kit) used in experiment.
SpecificationTechnical detailAnalog front end band2.4 GHz–2.5 GHzRx noise figure3.3 dBRx sensitivity−75 dBmRxI/Q gain0.1 dBRx phase imbalance0.35Based band gain control range33–62 dBFigure 2
Test bed for 3D indoor location.MATLAB/ADS is used to generate an OFDM signal, which is then upload to the VSG at 2.45 GHz carrier frequency, a 5 MHz bandwidth, and an 89.6 MHz sampling frequency. The VSG is connected to an RF switch with 5 outputs, each of which is further connected to a high gain omnidirection antenna using a coaxis cable. All BSs are clock-synchronized. The receiving antenna is directly connected to the Maxim board. The front-end is composed of a variable gain low-noise amplifier (LNA) followed by anI/Q zero-IF down converter, a programmable lowpass filter (LPF), and a variable gain amplifier (VGA) as shown in Figure 2. After downconversion, the baseband signal is sampled and quantified by an analog-to-digital converter. A MATLAB program is then used to capture both the I component and the Q component from the front end. The estimated time delay is obtained at the output of a matched filter, where the known time-multiplexed at the RF switch and the extra cable delay are considered for time delay compensation.


## 4. Empirical Result

### 4.1. Demonstration of Nakagami Distribution for LOS/NLOS Identification

As mentioned in [15], Nakagami distribution outperforms Rician distribution and Rayleigh distribution for LOS/NLOS identification because it is a more general form [15]. The following experiment verifies the availability of the Nakagami distribution for LOS identification, which is the first of our three objectives.The probability density function (PDF) of Nakagami distribution is given by(13)f(r)=2mmr2m-1Γ(m)Ωmexp⁡{-mr2Ω}r≥0,
where Γ(m) is the Gamma function, Ω, which is equal to E{r2}2, denotes the second moment, and m≥1/2, which is referred to as a shape factor, is defined as(14)m=E2{r2}E{[r2-E(r2)]2}.
The corresponding cumulative density function (CDF) can be expressed as(15)F(r)=P(mr2Ω,m),
where P(·) is the incomplete Gamma function. In the special case, when m=1, the Nakagami distribution behaves as a Rayleigh distribution. When m>1, it behaves as a Rician distribution. When m=0.5, it behaves as a one-sided Gaussian distribution.In our experiment, two BSs are assumed to be in a NLOS environment while three BSs are assumed to be in a LOS environment. 100 observations are captured from each BS, followed by anm-detector test to identify whether each observation belongs to a LOS environment or a NLOS one. One m is estimated from these 100 observations. For instance, Figure 3 shows several typical PDFs, their corresponding m factor values, and their corresponding channel types, that is, LOS or NLOS, where one can see that when m≤1, the channel is NLOS, and when m>1, the channel is LOS.Figure 3
PDF of LOS/NLOS channels and their corresponding Nakagami-m factors.To demonstrate the availability of the Nakagami distribution for LOS/NLOS identification, both a Kolmogorov-Smirnov (K-S) test and Chi-square test are performed at 5% significance level. Table2 illustrates their passing rate. One can see that the Nakagami m-factor has a 75% above passing rate for both tests. However, the passing rate of a LOS channel is shown to be higher than a NLOS one for both tests. This demonstrates that the LOS/NLOS channel can be properly identified using a Nakagami distribution for our indoor environment.Table 2
Passing rate of chi-square and K-S hypothesis tests.
MeasurementsChi-squareK-S TestLOS85.641276.06NLOS82.793375.71


### 4.2. Path Loss Model Estimation

Another objective in our experiment is to estimate the path loss exponent and shadowing effects of the experimental site. To generate this model, we measure the distance between each BS and MS as well as the received signal power. For a LOS test, we place all BSs in the lab’s Section2, select fifteen different positions for MS, and measure their received signal powers and their distances between each BS and MS. At each MS position, we stored a minimum of 20 signal power measurements and then move MS to another position. The same process is repeated for tests in a NLOS environment.Figure4 shows all received power measurements and its line of fit for pass loss exponent estimation in a LOS environment. Similarly, Figure 5 shows all received power measurements and its line of fit for pass loss exponent estimation in a NLOS environment. Table 3 illustrates the estimated path loss exponent and shadowing effect for both LOS and NLOS environments, where one can see that in a LOS environment, the path loss exponent and shadowing effect are found to be equal to 1.7 and 2.2, respectively; in a NLOS environment, the path loss exponent and shadowing effect are 3.2 and 2.7, respectively.Table 3
Estimated path loss exponent and shadowing effect.
Indoor EnvironmentPath Loss ExponentShadowing EffectLOS1.72.2NLOS3.22.7Figure 4
Path loss exponent estimation in a LOS environment.Figure 5
Path Loss Exponent Estimation in a LOS and NLOS combined environment.


### 4.3. Verification of Our Proposed HWLS RSS/TOA

The third objective of our experiment is to verify our proposed HWLS RSS/TOA for the indoor NLOS environments. For this purpose, (1) MS is placed in the lab’s 2nd Section; (2) two BSs are placed in the lab’s 1st Section and the 3rd Section, leading to a NLOS propagation channel and; three BSs are placed in the lab’s 2nd Section, leading to a LOS propagation channel. 100 observations are captured for this scenario. Using the derivation given in ((13a), (13b), (13c), and (13d)) through (15f), together with [19], all CRLBs for positioning estimation using 3D TOA alone, 3D RSS alone, 3D hybrid RSS/TOA, and our proposed 3D HWLS RSS/TOA are plotted in Figure 6 for comparison. One can see that the proposed 3D HWLS RSS/TOA outperforms all other three methods. Furthermore, the CDF of the root mean square error (RMSE) for positioning estimation using 3D TOA alone, 3D RSS alone, 3D hybrid RSS/TOA, and our proposed 3D HWLS RSS/TOA is shown in Figure 7, which is obtained from an indoor LOS/NLOS environments. Figure 7 indicates that the proposed HWLS RSS/TOA provides the best positioning accuracy compared to all other three methods. The overall performance of the considered methods for 3D positioning estimation is reported in Table 4, where a 2.02 m RMSE can be obtained by the proposed HWLS RSS/TOA, compared to 4.23 m by Hybrid RSS/TOA, 5.75 m by TOA alone, and 6.10 m by RSS alone.Table 4
Rmse for 3D positioning error using TOA alone, RSS alone, hybrid RSS/TOA and hwls RSS/TOA.
MethodCoordinatesVariance (m)Min Error (m)Mean Error (m)RMSE (m)TOAX1.501.0231.045.75Y0.701.0681.05Z1.801.121.38RSSX1.751.341.806.1Y1.201.901.72Z2.011.402.40Hybrid TOA/RSSX1.440.981.004.23Y0.500.970.99Z1.501.011.20HWLS TOA/RSSX0.0130.430.452.02Y0.0120.150.20Z0.280.1190.72Figure 6
CRLB for positioning estimation using TOA alone, RSS alone, hybrid RSS/TOA and the proposed HWLS RSS/TOA.Figure 7
CDF of positioning error using TOA alone, RSS alone, hybrid RSS/TOA, and the proposed HWLS RSS/TOA.


## 4.1. Demonstration of Nakagami Distribution for LOS/NLOS Identification

As mentioned in [15], Nakagami distribution outperforms Rician distribution and Rayleigh distribution for LOS/NLOS identification because it is a more general form [15]. The following experiment verifies the availability of the Nakagami distribution for LOS identification, which is the first of our three objectives.The probability density function (PDF) of Nakagami distribution is given by(13)f(r)=2mmr2m-1Γ(m)Ωmexp⁡{-mr2Ω}r≥0,
where Γ(m) is the Gamma function, Ω, which is equal to E{r2}2, denotes the second moment, and m≥1/2, which is referred to as a shape factor, is defined as(14)m=E2{r2}E{[r2-E(r2)]2}.
The corresponding cumulative density function (CDF) can be expressed as(15)F(r)=P(mr2Ω,m),
where P(·) is the incomplete Gamma function. In the special case, when m=1, the Nakagami distribution behaves as a Rayleigh distribution. When m>1, it behaves as a Rician distribution. When m=0.5, it behaves as a one-sided Gaussian distribution.In our experiment, two BSs are assumed to be in a NLOS environment while three BSs are assumed to be in a LOS environment. 100 observations are captured from each BS, followed by anm-detector test to identify whether each observation belongs to a LOS environment or a NLOS one. One m is estimated from these 100 observations. For instance, Figure 3 shows several typical PDFs, their corresponding m factor values, and their corresponding channel types, that is, LOS or NLOS, where one can see that when m≤1, the channel is NLOS, and when m>1, the channel is LOS.Figure 3
PDF of LOS/NLOS channels and their corresponding Nakagami-m factors.To demonstrate the availability of the Nakagami distribution for LOS/NLOS identification, both a Kolmogorov-Smirnov (K-S) test and Chi-square test are performed at 5% significance level. Table2 illustrates their passing rate. One can see that the Nakagami m-factor has a 75% above passing rate for both tests. However, the passing rate of a LOS channel is shown to be higher than a NLOS one for both tests. This demonstrates that the LOS/NLOS channel can be properly identified using a Nakagami distribution for our indoor environment.Table 2
Passing rate of chi-square and K-S hypothesis tests.
MeasurementsChi-squareK-S TestLOS85.641276.06NLOS82.793375.71


## 4.2. Path Loss Model Estimation

Another objective in our experiment is to estimate the path loss exponent and shadowing effects of the experimental site. To generate this model, we measure the distance between each BS and MS as well as the received signal power. For a LOS test, we place all BSs in the lab’s Section2, select fifteen different positions for MS, and measure their received signal powers and their distances between each BS and MS. At each MS position, we stored a minimum of 20 signal power measurements and then move MS to another position. The same process is repeated for tests in a NLOS environment.Figure4 shows all received power measurements and its line of fit for pass loss exponent estimation in a LOS environment. Similarly, Figure 5 shows all received power measurements and its line of fit for pass loss exponent estimation in a NLOS environment. Table 3 illustrates the estimated path loss exponent and shadowing effect for both LOS and NLOS environments, where one can see that in a LOS environment, the path loss exponent and shadowing effect are found to be equal to 1.7 and 2.2, respectively; in a NLOS environment, the path loss exponent and shadowing effect are 3.2 and 2.7, respectively.Table 3
Estimated path loss exponent and shadowing effect.
Indoor EnvironmentPath Loss ExponentShadowing EffectLOS1.72.2NLOS3.22.7Figure 4
Path loss exponent estimation in a LOS environment.Figure 5
Path Loss Exponent Estimation in a LOS and NLOS combined environment.


## 4.3. Verification of Our Proposed HWLS RSS/TOA

The third objective of our experiment is to verify our proposed HWLS RSS/TOA for the indoor NLOS environments. For this purpose, (1) MS is placed in the lab’s 2nd Section; (2) two BSs are placed in the lab’s 1st Section and the 3rd Section, leading to a NLOS propagation channel and; three BSs are placed in the lab’s 2nd Section, leading to a LOS propagation channel. 100 observations are captured for this scenario. Using the derivation given in ((13a), (13b), (13c), and (13d)) through (15f), together with [19], all CRLBs for positioning estimation using 3D TOA alone, 3D RSS alone, 3D hybrid RSS/TOA, and our proposed 3D HWLS RSS/TOA are plotted in Figure 6 for comparison. One can see that the proposed 3D HWLS RSS/TOA outperforms all other three methods. Furthermore, the CDF of the root mean square error (RMSE) for positioning estimation using 3D TOA alone, 3D RSS alone, 3D hybrid RSS/TOA, and our proposed 3D HWLS RSS/TOA is shown in Figure 7, which is obtained from an indoor LOS/NLOS environments. Figure 7 indicates that the proposed HWLS RSS/TOA provides the best positioning accuracy compared to all other three methods. The overall performance of the considered methods for 3D positioning estimation is reported in Table 4, where a 2.02 m RMSE can be obtained by the proposed HWLS RSS/TOA, compared to 4.23 m by Hybrid RSS/TOA, 5.75 m by TOA alone, and 6.10 m by RSS alone.Table 4
Rmse for 3D positioning error using TOA alone, RSS alone, hybrid RSS/TOA and hwls RSS/TOA.
MethodCoordinatesVariance (m)Min Error (m)Mean Error (m)RMSE (m)TOAX1.501.0231.045.75Y0.701.0681.05Z1.801.121.38RSSX1.751.341.806.1Y1.201.901.72Z2.011.402.40Hybrid TOA/RSSX1.440.981.004.23Y0.500.970.99Z1.501.011.20HWLS TOA/RSSX0.0130.430.452.02Y0.0120.150.20Z0.280.1190.72Figure 6
CRLB for positioning estimation using TOA alone, RSS alone, hybrid RSS/TOA and the proposed HWLS RSS/TOA.Figure 7
CDF of positioning error using TOA alone, RSS alone, hybrid RSS/TOA, and the proposed HWLS RSS/TOA.


## 5. Conclusion

In this paper, a novel HWLS RSS/TOA method has been proposed for target localization in an indoor NLOS environment. An experiment is conducted in the iRadio lab, in the ICT building at the University of Calgary, to demonstrate our proposed algorithm. Experimental results indicated that HWLS exhibits a much better performance than TOA alone, RSS alone, and hybrid RSS/TOA in terms of positioning accuracy in an indoor LOS/NLOS environment.



---

*Source: 503707-2012-03-01.xml*