Scholarly article on topic 'A Hybrid TOA-Fingerprinting Based Localization of Mobile Nodes Using UWB Signaling for Non Line-Of-Sight Conditions'

A Hybrid TOA-Fingerprinting Based Localization of Mobile Nodes Using UWB Signaling for Non Line-Of-Sight Conditions Academic research paper on "Electrical engineering, electronic engineering, information engineering"

Share paper

Academic research paper on topic "A Hybrid TOA-Fingerprinting Based Localization of Mobile Nodes Using UWB Signaling for Non Line-Of-Sight Conditions"

Sensors 2012, 12, 11187-11204; doi:10.3390/s120811187



ISSN 1424-8220


A Hybrid TOA-Fingerprinting Based Localization of Mobile Nodes Using UWB Signaling for Non Line-Of-Sight Conditions

Md. Humayun Kabir * and Ryuji Kohno

Division of Physics, Electrical and Computer Engineering, Yokohama National University,

79-5 Tokiwadai, Hodogaya, Yokohama 240-8501, Japan; E-Mail:

* Author to whom correspondence should be addressed; E-Mail:;

Tel.: +81-45-339-4116; Fax: +81-45-338-1157.

Received: 2 July 2012; in revised form: 3 August 2012 /Accepted: 7 August 2012 /

Published: 10 August 2012

Abstract: Recently, Impulse Radio Ultra Wideband (IR-UWB) signaling has become popular for providing precise location accuracy for mobile and wireless sensor node localization in the indoor environment due to its large bandwidth and high time resolution while providing ultra-high transmission capacity. However, the Non-line-of-sight (NLOS) error mitigation has considerable importance in localization of wireless nodes. In order to mitigate NLOS errors in indoor localization this paper proposes and investigates a novel approach which creates a hybrid combination of channel impulse response (CIR)-based fingerprinting (FP) positioning and an iterative Time of Arrival (TOA) real time positioning method using Ultra Wideband (UWB) signaling. Besides, to reduce the calculation complexities in FP method, this paper also introduces a unique idea for the arrangement of reference nodes (or tags) to create a fingerprinting database. The simulation results confirm that the proposed hybrid method yields better positioning accuracies and is much more robust in NLOS error mitigation than TOA only and FP only and a conventional iterative positioning method.

Keywords: UWB signaling; node localization; target tracking; fingerprinting; TOA; NLOS

1. Introduction

Lately, the Impulse Radio (IR)-Ultra Wideband (UWB) has been applied for the wireless personal area network (PAN), body area network (BAN), and RFID standards IEEE802.15.4a, 802.15.6, and

802.15.4f, respectively, because IR-UWB can perform ultra-high resolution of ranging, positioning, and high speed data transmission. Among the variety of prospective applications of IR-UWB, one of the most promising is in wireless sensor networks (WSNs), which require both robust communications and high-precision ranging capabilities [1]. Moreover, IR-UWB has some inherent properties like low power consumption, low cost, little interference to other systems, resistance to severe multipath conflicts and jamming, and has very high time domain resolution allowing for precise location and tracking, making it well suited for WSN applications [1].

Since IR-UWB signals have very short duration pulses, they can provide very accurate ranging and positioning capability in short range indoor radio propagation environments. Besides, the high time resolution characteristic of the UWB signal makes the time of arrival (TOA) method a good choice for location estimation in UWB communications. However, UWB localization and communication have to face lots of challenges like signal acquisition, multipath effects, multi-user interference and NLOS propagation. Among them NLOS propagation introduces positive bias for ranging estimation. Hence, it is difficult to recognize the impact of NLOS conditions on localization systems and to develop techniques that mitigate their effects [2].

Although, with enhanced sampling rates UWB receivers can accomplish high TOA resolution, their operating ranges are typically limited in reduced interference [3]. On the other hand, a suitable alternative to this real-time geometric positioning method is RF fingerprinting (FP) or pattern matching method, where a pre-calculated location database is used to match the properties of the received signal for position estimation. Moreover, the FP method is well suited in positioning under dense multipath conditions.

Besides the FP method, to achieve higher accuracy in localization of mobile nodes, many conventional iterative methods like Gauss-Newton, Steepest Descent, Levenberg-Marquardt have been proposed in conjunction with TOA or TDOA methods for mitigating NLOS errors. Among them, the Steepest Descent algorithm shows the slowest convergence in the final iteration steps, but for maintaining higher accuracy and low complexity the Levenberg-Marquardt method is the most suitable one to estimate a node location among the iterative methods [4]. Moreover, these iterative methods require initial positioning guesses. If the guesses are not accurate, convergence problems may occur. Alternately, a non-iterative positioning method can provide good initial positions for an iterative positioning method [5]. This idea also leads to an energy efficient iterative localization scheme for wireless sensor nodes. In [6], a hybrid combination of pattern matching (PM) and TDOA positioning methods is proposed for a CDMA network, although the described method has low accuracy in an indoor environment where the GPS signal is not satisfactory.

In our proposal, we introduce a hybrid positioning technique which combines a Channel Impulse Response (CIR)-based FP method with polygonal arrangement (which will be discussed in a later section) of reference nodes (or tags) and an iterative-TOA positioning method using UWB signaling for wireless ad hoc networks such as PAN, BAN, and RFID as well as wireless sensor nodes. The proposed method assures significant improvement in positioning accuracy compared to TOA only, FP only and a conventional iterative-TOA positioning methods by reducing NLOS errors effectively in the indoor environment. Moreover, we employ a deterministic wave propagation model based on Geometrical Optics to predict indoor radio coverage and CIRs.

In this paper, we would like to concentrate on the successful implementations of our proposed hybrid method in an indoor radio propagation environment such as a typical office room. We shall discuss on how to utilize the positioning results of the CIR-based FP method with different polygonal arrangements (which will be discussed in Section 3) of reference nodes in an iterative-TOA positioning method and to get the mutual benefits of the TOA and FP methods in a hybrid method. Besides, we would like to limit our discussion on the positioning performances of FP and proposed hybrid localization considering various NLOS incidences experienced by access points (AP)s. We also compare the performances of proposed hybrid method with TOA only, FP only and a conventional iterative-TOA method in LOS cases as well as worst cases of NLOS. Besides, the detailed mathematical analysis of computational complexity of the proposed hybrid method is beyond the scope of this paper and hopefully would be provided in a future issue.

The remaining of the paper is organized as follows: Section 2 illustrates the preparations for our proposal for hybrid positioning algorithm. The detailed analysis for our proposal for hybrid algorithm is provided in Section 3. Section 4 summarizes the model's performance using computer simulations. Finally, a brief conclusion is made in Section 5.

2. Preparations for Our Proposal for Hybrid Positioning Algorithm

Before introducing our hybrid positioning algorithm the following three subsections are the necessary preparations for our proposal. Section 2.1 discusses and proposes a deterministic channel model constructed by a Ray Tracing technique for indoor radio wave propagation. Then, the CIR-based FP positioning method for an indoor propagation environment is explained in Section 2.2. Section 2.3 provides an overview on UWB TOA-based accurate ranging in a multipath environment.

2.1. Deterministic Channel Model for Indoor Radio Wave Propagation

We consider a typical well-furnished office room that provides an efficient means of simulating a short range indoor UWB channel. Due to reflection, transmission, refraction, and diffraction by obstacles signal propagation suffers from severe multipath effects in an indoor environment [7]. In order to predict the accurate indoor radio coverage and channel impulse response we consider a deterministic (site-specific) propagation modeling technique, called Ray Tracing (RT). The RT technique utilizes the knowledge of the locations and electromagnetic properties of indoor objects, is used to predict path loss, time invariant impulse response, and the RMS delay spread [8-10].

A general model of band-limited complex channel impulse response (CIR) is expressed as [11]:

h(t) = X AW-T )ej (1)

where, An is the path attenuation, Tn is the time delay of the nth path, S is the Dirac delta function, 6n is the phase of the nth path. For the implementation of this model, it is necessary to identify the amplitudes, time delay and phases of the N number of components of the response.

We employed the RT technique to identify the components of the above mentioned channel model. RT technique follows the ray launching approach (based on geometrical optics), which involves a number of rays launched uniformly in space around the transmitter antenna (Figure 1). Each ray is

traced until it reaches the receiver or its amplitude falls below a specified limit. In our model, the specular reflections and transmissions are considered while diffraction and scattering are neglected, so that every wave component traverses one of more free space propagation paths between the transmitter and the receiver.

Figure 1. Ray tracing for indoor propagation modeling.

The complex electric field incident at the receiver due to the nth impinging ray can be represented as:

En = E0

2 ( Qn Rn ^

2 mn r

^ q=1 r=1 J

- j2Mdn / 2c


E0: transmitted electric field

Xc: wave length correspond to center frequency

dn: propagation path length of nth ray

Rn: number of reflections of nth ray

Qn: number of penetrations of nth ray

rr: coefficient of rth reflection of nth ray

Pq: coefficient of qth penetration of nth ray

In Equation (2), the term_c— is the freespace pathloss, the exponential term e~J2ndn Uc represents

phase offsets due to pathlength dn, assuming omni-directional antennas where the azimuthal-plane antenna gains have been ignored for simplicity.

2.2. CIR Based FP Positioning Technique for Indoor Environment

The estimated CIRs described in the previous subsection are utilized for FP positioning and TOA based ranging (will be described in Section 2.3). The contents of Sections 2.2 and 2.3 are the preparations for the hybrid algorithm described in Section 3.

The CIR-based FP positioning, is also known as pattern matching technique, is especially suited for positioning under dense multipath conditions [12]. This technique consists of two phases:

(1) In the training phase, location signatures (fingerprints) based on the received signal's CIRs' are measured by several access points (APs) for various reference locations of nodes in a region to construct a fingerprinting database.

(2) In the positioning phase, the run time estimation of node's location is performed by correlating the CIRs' of received signal with the fingerprints stored in the database.

Generally, the FP technique requires exhaustive calculations for fingerprinting pattern matching as well as creating a database for storing the location fingerprints (signatures). Therefore, the computational complexity in the FP technique depends greatly on the size of the fingerprinting database employed in pattern matching. As the fingerprints of the reference nodes constitute a database, the larger the number of reference nodes the larger the size of the fingerprinting database needed to cover the whole area, but employing a smaller number of reference nodes reduces the computational complexity, although it will not provide a good accuracy in positioning. In order to combat this positioning constraint we introduce a new idea to optimize the number of reference nodes to cover the whole area as well as to generate a fingerprinting database of optimum size.

In our new idea, we arrange the reference nodes in the room to form polygon shapes (as shown in Figures 2 and 3) instead of deploying them in a scattered way. The number of reference nodes required to cover the whole area depends on the size and quantity of such polygons. Similarly, the size of fingerprinting database also depends on the size of the polygons. Therefore, the optimum polygon size will reduce the computational complexity in fingerprinting pattern matching to a reasonable level while maintaining a certain degree of positioning accuracy in the FP method, so the main idea of arranging reference nodes (or tags) to form polygons is to reduce unnecessary deployment of reference nodes which increase the size of the fingerprinting database. Moreover, by varying the size of such polygons we can vary the size of the fingerprinting database as well as the complexity of the FP method calculations.

In the training phase of the FP method, signatures (fingerprints) h(l,t) are estimated and recorded for various reference locations l £ R, where R is the region of interest. Location of a node (or, tag) in 2D horizontal plane is defined as l(x,y) in Cartesian coordinates.

Figure 2. Size of polygon.

50 cm-Tag1

Figure 3. Arrangements of reference nodes in the room forming polygonal shapes.

In the positioning phase, an estimate l0of the true instantaneous position l0, is obtained by using the corresponding instantaneous estimated CIR, h( Z0,t). This estimated CIR is corrupted by channel noise. For, m = h(l, t) and n = h( Z0,t), the channel spatial correlation RT" is defined as:

E{mn} - E{m}E{n}

Rm,n _

^ E{\m\2} - |E{m}\2)(E{|n\2} - |E{n}|2)]

where E{.} denotes the expectation operator. The position of the node is estimated by maximizing the correlation coefficient Ri"'" , i.e.,

L _ arg max

Due to imperfect channel estimation at the receiver the ideal correlation cannot be attained. Hence, a correlation threshold Rth is used such that, Rf > Rth. The values of R?" crosses Rth are defined as

RT" = [Ri], R12,.....,Rik], where k is the number of values of Rf crosses Rth Then the values of RT" are

indexed with their corresponding location coordinates in descending order. Finally, the node location is calculated by taking a weighted average of the coordinates of the first three locations in the index, assuming these three locations are very close to the node location.

2.3. UWB TO A Based Ranging in Muitipath Propagation Environment

The TOA estimation technique used with UWB transmission can be used for conducting accurate ranging in indoor multipath environments. For the indoor multipath channel, the impulse response is usually defined as [13]:

h(t) _ x ake^kö(t-rk)

where Lp is the number of multipath components and ak, (pk and Tk are the amplitude, phase and propagation delay of the kth path, respectively. The range measurement between a transmitter and a receiver is typically corrupted by multipath fading, thermal noise, blockage of direct path (i.e., NLOS path) and direct path excess delay [14].

For TOA estimation, we employ a threshold-based first path detection method which utilizes the strongest path of CIR. This method makes use of an iterative search back algorithm, referred in [14], to calculate the noise floor to be used in the detection of the first path.

Figure 4. Search-back detection method.

First peak Strongest peak

- Threshold p ;

Noise Floor J j >

It Iff if

Delay time -

3. Proposal for a Hybrid Localization Algorithm

3.1. A Brief Idea of the Proposal for Hybrid Localization

The deterministic channel model described in Section 2.1, the FP positioning method explained in Section 2.2, and followed by UWB TOA-based ranging technique described in the Section 2.3, are the preparation for the hybrid localization that will be described in this section. Our proposal for hybrid localization is a combination of the CIR-based FP positioning and iterative-TOA positioning methods.

It is obvious that an iterative positioning algorithm needs an initial guess which will influence the iterative process. A good approximation of the initial value also reduces the calculation complexity as well as increases the positioning accuracies. In our proposed method, we initialize the iterative-TOA positioning with position location estimated by the CIR-based precalculated FP method. As the CIR- based FP method has better positioning accuracy than the TOA method in dense multipath indoor conditions, the location point estimated in of the FP method acts as a better initial value for our hybrid method. Moreover, by varying the size of the fingerprinting polygons (as discussed in the previous section), the positioning errors as well the complexities in the proposed hybrid method could be minimized. The brief idea of the proposed localization algorithm is provided in the flow chart in Figure 5.

3.2. Hybrid Positioning Algorithm Combining CIR Based FP and Iterative-TOA Positioning Methods

In hybrid positioning, the node location estimated in the FP method is considered as the accurate initial values for an iterative-TOA method. The Levenberg-Marquardt based iterative algorithm (will be described in the later part of this section) is used to mitigate NLOS error in the iterative process of the hybrid positioning.

Let M access points (AP)s with position coordinates (x, y), where i = 1, 2,..M. Assume (x0, y0) is the sensor node location to be determined and it is away from access points (AP)s. In the NLOS case, the measured range R is always larger than the LOS range. This inequality can be expressed as:

LOS range = ^ ^ i = 1,2,-m

Figure 5. Flow chart for the proposed positioning algorithm.

Let us define the range between position (x, y) and the AP, in LOS condition as:

di (xy)=^(X-xT+cy—yf ,

i = 1,2... M

As d(x, y) is a nonlinear function, the measurement of d(x, y) is erroneous in accordance with the severity of nonlinearity. This ranging error will be minimum at the immediate vicinity of the minimum of the cost function [15]. The nonlinear model of ranging can be represented as an optimum interpolation between a Taylor series method and a gradient method.

If (xes, yes) is an estimated position of the true position of node, expanding the function d(x, y) in a Taylor's series we get:

d,(x y) « dt ^ yes)+

ddi (x y)

(X - Xes) +

dd,(X; y)

=es (y - y )

=y„v yes/

In Equation (8), only the terms of zero-order and first-order are kept. Let 00 = [x0, yo] be the position vector and D(00) = [d1(^0), d2(00),..., dM(00)]T be the distance vectors. According to Equation (8), the distance vectors D(00) can be written as:

D (00 ) ~ D0) + J0 ) (00 - 0es) (9)

where, the Jacobian matrix J(0es) can be expressed as:

x - x, y - y,

es 1 y es y 1

3(0es ) =

du (Oss ) dM (6es )

Thus, the most M linear independent ranges are:

D(0) = [Di(fl), D2(0), ...Dm^)]1 and the corresponding observed ranges are given by:

R = [R;, R2,.....Rm]

Then, according to [16,17] the cost function is defined as:

°(O)=2 - D (O)]

(ii) (12)

which should be minimized for the optimum case. The minimization of the cost function is shown as:

9 = arg min O (9)

Then, we assemble the individual components of [Ri -Di(6)] from Equation (13) into a residual vector Rn^Rm defined by:

e(ff) = \^s1(dls2(d),

Using the notation || , ||2 for Euclidean norm, we can rewrite Equation (13) as:

O(ff) = VA\e\\\ (16)

Consider the linear case, where every £(ff), i = 1,2,...,M is linear and Jacobian J of s is constant. Then the cost function O is given by the quadratic as [16,17]:

o(ff) = / \\ J.e ^ £\\2 (17)

V Off) = JT ( J.ff +£) (18)

V2 O(ff) = JT J (19)

Setting V O(ff) = 0 yields 6U = - ( JT J)-1 JT£,which is the solution to the set of normal equations. For the non-linear case, as mentioned in [16,17],we have:

V O(ff) = JT £ (20)

V2 O(ff) = JT J (21)

In order to estimate the node position by reducing NLOS error, it is required to minimize the cost function O(ff). We applied the Levenberg-Marquardt [15,17] iterative algorithm that minimizes O(ff) and can be expressed as:

ff n+1 = ff n + ( JT J + 1kdiag [JT J] )1 JT (ff n) £"

= ff n +( Hk + 1kdiag(H) )-1 gk (22)

The Levenberg-Marquardt algorithm starts with an initial value, ff(0) = [x(0),y(0)]T which is obtained from the FP method. In Equation (22), 1 is the damping factor; whose value is reconfigured during iteration in order to reduce the influence of gradient descent. The value of 1 affects both direction and step-size of iterations. Therefore, the initial value of 1 should be chosen according to the size of the elements in A0 = J(ff°)T J(ff0), such that [18]:

10 = t . maxi { aif0)} (23)

Here, t is defined by the user (by a rule of thumb, but a small value should be chosen). Throughout iteration the size of 1 is updated and is controlled by the cost function O(ff). In an iteration, a smaller value of O(ff) indicates a good approximation to ff n+1, and we can keep the weights at their new values and decrease 1 by a certain (say, 10) factor so that the next Levenberg-Marquardt step is closer to the Gauss-Newton step. In contrast, if O(ff) is large, then ff n+1 is a poor approximation, and we should reset the weights at their previous values and increase 1 by same factor for getting closer to the direction of Steepest Descent and reduce the steps.

4. Simulation Models and Performance Evaluations of Proposed Hybrid Method

This section provides the performances of our proposed hybrid method for node positioning by computer simulations. The performances of FP positioning method with various fingerprinting polygons (discussed in Section 2.2) are also evaluated in this section. Moreover, a comparison of the proposed hybrid method with TOA only, FP only and a conventional iterative method is shown in the later part of the section.

Tables 1 and 2 show the parameter values for evaluating model's performances and the properties of materials we employed in the RT technique, respectively.

Table 1. Parameter values for evaluating model's performance.


Transmitted signal Channel bandwidth Sampling frequencies AWGN parameter Transmitted Power

Gaussian monopulse

2 GHz 30, 50, 80 GHz.

g = 0.5 m -60 dBm

Fingerprinting polygon size Number of Access Point Area

20 cm, 30 cm, 50 cm 4

8 m x 6 m

Table 2. Properties of the materials employed in RT technique.


g[S/m] gr Thickness[cm]

Concrete wall Table, chair, cabinet (wood) _Window (glass)_

0.01 9 7.5

10-5 13 3.0

10-12 7.6 3.0

To get the accurate results in the ray launching type RT technique, many rays have to be launched and only a fraction of these reach the receiver. The accuracy of this technique also depends on the radius of the reception sphere. If it is too small, rays will pass by. If it is too large, paths might be duplicated. We set the radius of the reception sphere at 1.5 cm. Moreover, we consider the reflection of rays up to the third order, since the received signal strength beyond third order ray reflection becomes very low and is neglected. CIRs collected by the RT technique are convoluted with a UWB signal with bandwidth of 2 GHz. To consider a sufficient amount of multipath effects in the CIR pattern for NLOS incidences we allow a wider observation window, i.e., 65 ns, since the probability that the direct path is further apart from the strongest path is higher in the NLOS than in LOS. Therefore, it could accommodate maximum of 20 paths in the CIR patterns of NLOS cases.

As we consider an indoor office room of size 8 m x 6 m for our simulations, we think four APs are enough to cover radio mapping of this limited area and adequate for CIR-based FP positioning. It is obvious that increasing the number of AP would increase the possibility of LOS incidences and thus improve the accuracy in TOA-based positioning. However, increasing the number of APs will also increase the size of the fingerprinting database and thus increase the computational complexity of positioning in the FP method.

To evaluate the performance of the FP method for positioning of the same test nodes which are used in an iterative TOA positioning, we place the reference nodes (or tags) in such a manner to form polygon shapes (as mentioned in Section 2.2). We employ three different sizes of polygon, 20 cm, 30 cm and 50 cm and evaluate the positioning accuracy by varying signal's sampling rate (i.e., 30 GHz, 50 GHz and 80 GHz) for each case. Here, Figures 6-8 show cumulative distribution function (CDF) of distance errors (distance between estimated and real positions of test nodes) in the FP method for LOS cases as well as worst cases of NLOS. The NLOS rates shown in the figures are the ratios of the number of access points (AP)s experiencing NLOS among all APs. For instance, a NLOS rate of 0.5 means the half of the total number of APs are experiencing NLOS.

Figure 6 depicts that, in LOS cases, 50% probability is achieved for the accuracy of less than 10 cm in cases with a polygon size of 20 cm, whereas, 35% and 15% probability is achieved for the cases with polygon size of 30 cm and 50 cm, respectively, for the accuracy of less than 10 cm. Hence, cases with polygon size of 20 cm have better positioning accuracies than other polygon cases in LOS. This is because in cases with smaller polygons the reference nodes (or tags) are in close proximity to the test nodes positions as compared to the cases with larger polygons. Therefore, in smaller polygons, there are fewer chances of errors in matching the CIR pattern of the test node with that of its close neighbouring reference nodes under LOS conditions.

Figure 6. CDF of the distance errors for FP method in different cases with LOS.

0.9 8.0-8 I 0.7

¿0.6 c

|o-3 □

¡0.2 0

FP method in LOS case

I M___

■¿-polygon 20 cm & fs:30 GHz.

Jfi* -0-polygon:20 cm & fs:50 GHz '"fr"polygon:20 cm & fs:80 GHz.

0.3 0.4 Distance error (m)

Figure 7. CDF of the distance errors for FP method in different cases with NLOS rate = 0.5.

Distance error (m) Distance error (m)

(a) (b)

In NLOS situations, unlike LOS conditions the position estimation mostly depends on the multipath components of the received signal. Therefore, the estimated position would be very far away from the real position of the test node in NLOS conditions. As the estimated position is shifted from its real position, it is more likely to match the CIR pattern of estimated location with that of a far location point rather than matching CIR pattern of the point near to its real position. Figures (7) and (8) show the CDF of distance error in the FP method for worst cases of NLOS. With a NLOS rate of 0.5, the case with polygon size of 50 cm (with sampling rate 50 GHz.) has a 35% probability to achieve accuracy of less than 1.5 m, whereas the cases with polygon sizes of 30 cm and 20 cm have a probability of 20% to achieve accuracy of less than 1.5 m. Therefore, the case with polygon size of 50 cm (with sampling rate 50 GHz.) performs slightly better than other cases with different polygons because when the NLOS rate is 0.5, two of the APs cannot see the test node and CIR patterns measured by these NLOS APs are also considered for node positioning. These two NLOS CIR patterns of the test node could have a better match with the CIR patterns of reference node points of a polygon

size of 50 cm (which are relatively far away from the real position of the test node) rather than matching with the CIR patterns of reference node points of polygon size 30 cm and 20 cm. Similarly, when the NLOS rate is 0.75, three out of four APs cannot see the test node and CIR patterns measured by these three NLOS APs are also considered for node positioning. Therefore, in these worst NLOS conditions, for the reasons stated above, the case with polygon size of 50 cm performs slightly better than other cases.

Besides, in the above results, the cases with higher received signal sampling rate sometimes provides slightly better positioning accuracies. Because, the higher sampling rate of the received signal reduces the time quantization error in the time delay estimation. Therefore, positioning accuracy increases with higher time resolutions.

It is obvious that, the smaller the polygon size the larger the fingerprinting database as well as the higher calculation complexity will be in the FP method. Alternatively, the cases with the larger polygon size provide lesser calculation complexity in FP positioning.

Figure 8. CDF of the distance errors for FP method (8a-c), Hybrid method (8d) in different cases with NLOS rate = 0.75.

2.4 2.5 Distance error (m)

23 24 25 Distance error (m)

FP method in NLOS rate of 0.75

-£-polygon:30 cm & fs:30 GHz -0-polygon:30 cm & fs:50 GHz polygon :30 cm&fs:80GHz

Hybrid method in NLOS rate of 0.

-A polygon;30 cm & fs;30 GHz, polygon:30 crn & fs:50 GHz. * polygon:30 cm & fs:80 GHz.

2.3 2.4 2.5 Distance error (m)

2.4 2.5 Distance error (m)

Figure 8(c,d) depicts the comparison of positioning accuracies in between FP and the proposed hybrid method for the cases with polygon size 30 cm when the NLOS rate is 0.75. Figures depict that the proposed hybrid method has better positioning accuracies by mitigating NLOS errors than the FP method. Figure 9 shows the performances of the hybrid positioning in root mean squared error (RMSE) for various cases with polygon sizes in LOS as well as different NLOS conditions. Moreover, the smaller polygon size in FP positioning provides better location accuracies in the hybrid method for lower rates of NLOS, whereas, for higher rates of NLOS, the positioning accuracies for the cases with various polygon sizes are almost identical.

Figure 9. Positioning accuracies in hybrid method for various cases with polygon sizes in LOS and different NLOS conditions.

In the proposed hybrid method, there is possibility to get mutual benefits from the TOA only and FP only methods. If any case in FP method does not have better positioning accuracy than the TOA only method, so the hybrid method can choose its initial value from the positioning result of TOA only method for iterations.

The positioning accuracies for different cases using the proposed hybrid method can be increased by increasing the iteration number. However, as depicted in Figure 10, the accuracy does not increase beyond 4 or 5 iterations for various NLOS cases. Hence, setting a maximum iteration number beyond 4 or 5 would not be necessary. Moreover, in Figure 10, there is a wide gap between the RMSE performances of NLOS rates of 0.25 and 0.5. This is because we have a better positioning accuracy in the FP method using a polygon size of 20 cm when the NLOS rate is 0.25. Therefore, in hybrid positioning, RMSE for the case of a NLOS rate of 0.25 decreases and it is closer to the positioning accuracy of the case with LOS where the polygon size is 20 cm. Besides, we don't have better positioning accuracy for the case using a polygon size of 20 cm than the case of polygon size of 50 cm in the FP method while NLOS rate is 0.5. Therefore, in hybrid positioning, RMSE performances in

each iteration for the case of NLOS rate of 0.5 are not good enough. Moreover, it is obvious that if we increase the number of test node samples in the FP method the positioning accuracy for the case of NLOS rate of 0.5 would also increase a bit.

Figure 10. Performance of the hybrid method in each iteration for different NLOS conditions, considering fingerprinting polygon size: 20 cm and signal sampling rate: 50 GHz.

Figure 11. Comparison on positioning accuracies for TOA only, FP only, conventional iterative TOA and hybrid methods in NLOS conditions.

Comparison of methods

gl □É

0 A TOA only t? FP with potygon:30 cm. ■0"Hybrid considering FP -O-Conventionallterative-TOA M

a ■ y/

fr?/ . FT > / ' x3C / /H y

.('■ / /

A / / / ✓ / ✓

.*>' ry ......

04 0.6

NLOS rale

Figure 11 shows that the proposed hybrid method has better positioning accuracies by mitigating NLOS error than TOA only, FP only and a conventional iterative method under different NLOS conditions. In a conventional iterative TOA method, we also applied the Levenberg-Marquardt iteration algorithm; however, initial values are random guesses unlike considering the positioning results of the FP method. Similarly, the proposed method also renders better positioning accuracies than the Gauss-Newton based iterative-TOA positioning method, reported in [19].

5. Conclusion

We have proposed a hybrid approach for localization of mobile or wireless sensor nodes using UWB signaling which combines the CIR-based FP positioning and iterative-TOA real time positioning methods. The proposed hybrid method follows a Levenb erg-Mar quardt based iterative algorithm, which effectively reduces NLOS errors. Our simulation results show that the proposed hybrid method yields better performances and is more robust in NLOS error mitigation than TOA only, FP only and a conventional iterative method. Moreover, the smaller size of polygons, resembling reference nodes (or tags) arrangement in the FP method, the bigger the size of fingerprinting database as well as the higher calculation complexity in positioning. Alternately, the larger size of polygon provides a smaller fingerprinting database compromising positioning accuracies. The optimum size of polygon could be predicted while evaluating the trade-off between accuracies in the proposed hybrid method and the complexities in positioning calculations, which is beyond the scope of this paper. A detailed analysis of the computational complexity of the proposed method will be provided in our next paper on this topic.


1. Zhang, J.; Orlik, P.V.; Sahinoglu, Z.; Molisch, A.F.; Kinney, P. UWB systems for wireless sensor networks. Proc. IEEE 2009, 97, 313-331.

2. Marano, S.; Gifford, W.M.; Wymeersch, H.; Win, M.Z. NLOS Identification and mitigation for localization based on UWB experimental data. IEEE J. Sel. Areas Commun. 2010, 28, 1026-1035.

3. Jin, Y.; Soh, W.-S.; Wong, W.C. Indoor localization with channel impulse response based fingerprint and nonparametric regression. IEEE Trans. Wirel. Commun. 2010, 9, 1120-1127.

4. Mensing, C.; Plass, S. Positioning Algorithms for Cellular Networks Using TDOA. In Proceddings of the IEEE International Conference on Acoustics, Speech and Signal Processing, Toulouse, France, 14-19 May 2006; pp. 1-4.

5. Xiong, J-Y.; Wang, W.; Zhu, Z.-L. A New TDOA Location Technique Based on Taylor Series Expansion in Cellular Networks. In Proceedings of the Fourth International Conference on Parallel and Distributed Computing, Applications and Technologies, Chengdu, China, 27-29 August 2003; pp. 378-381.

6. Akbarov, D.; Choi, H.; Park, Y.; Han, G.; Moon, J. Hybrid Pattern Matching/TDOA Positioning Method for CDMA Networks. In Proceeding of the 4th Workshop on Positioning, Navigation and Communication, Hannover, Germany, 22-22 March 2007; pp. 199-203.

7. Hashemi, H. The indoor radio propagation channel. Proc. IEEE 1993, 81, 943-968.

8. McKown, J.W.; Hamilton, R.L. Ray tracing as a design tool for radio networks. IEEENetw. Mag. 1991, 5, 27-30.

9. Schaubach, K.R.; Davis, N.J.; Rappaport, T.S. A Ray Tracing Method for Predicting Path Loss and Delay Spread in Microcellular Environments. In Proceedings of the 42nd IEEE Vehicular Technology Conference, Denver, CO, USA, 10-13 May 1992; pp. 932-935.

10. Diskin, J.; Brennan, C. Accelerated Ray-Tracing for Indoor Ultra-Wideband Propagation Modeling. In Proceedings of the 65th IEEE Vehicular Technology Conference, Dublin, Ireland, 22-25 April 2007; pp. 418-422.

11. Turin, G.L.; Clapp, F.D.; Johnston, T.L.; Fine, S.B.; Lavry, D. A statistical model of urban multipath propagation. IEEE Trans. Veh. Technol. 1972, 21, 1-9.

12. Malik, W.Q.; Allen, B. Wireless Sensor Positioning with Ultra Wideband Fingerprinting. In Proceeding of the 1st European Conference on Antennas and Propagation, Nice, France, 6-10 November 2006; pp. 1-5.

13. Alsindi, N.A.; Alavin, B.; Pahlavan, K. Measurement and modeling of Ultrawideband TOA-based ranging in indoor multipath environments. IEEE Trans. Veh. Technol. 2009, 58, 1046-1058.

14. Dashti, M.; Ghoraishi, M.; Takada, J. High-Precision Time-of-Arrival Estimation for UWB Localizers in Indoor Multipath Channels. In Novel Applications of the UWB Technologies; Lembrikov, B.I., Ed.; InTech: Rijeka, Croatia, 2011; pp. 397-422.

15. Marquardt, D.W. An algorithm for least-squares estimation of nonlinear parameters. J. Soc. Ind. Appl. Math. 1963, 11, 431-441.

16. Nocedal, J.; Wright, S.J. Numerical Optimization; Springer-Verlag Inc.: New York, NY, USA, 1999.

17. Ranganathan, A. The Levenberg-Marquardt Algorithm, Honda Research Institute; 8 June 2004. Available online: (accessed on 15 June 2012).

18. Madsen, K.; Nielsen, H.B.; Tingleff, O. Methods for Non-Linear Least Squares Problems, 2nd ed.; Lecture notes, Informatics and Mathematical Modelling, Technical University of Denmark: Lyngby, Denmark, April 2004; p. 25.

19. Kabir, M.H.; Kohno, R. A Hybrid Positioning Approach by UWB Radio Communication Systems for Non Line-Of-Sight Conditions. In Proceedings of the IEEE Global Telecommunications Conference (GLOBECOM 2011), Houston, TX, USA, 5-9 December 2011; pp. 1-5.

© 2012 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (

Copyright of Sensors (14248220) is the property of MDPI Publishing and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.