Journal Search Engine
Search Advanced Search Adode Reader(link)
Download PDF Export Citaion korean bibliography PMC previewer
ISSN : 1598-7248 (Print)
ISSN : 2234-6473 (Online)
Industrial Engineering & Management Systems Vol.16 No.1 pp.92-102

A Study on Particle Filter based on KLDResampling for Wireless Patient Tracking

Ly-Tu Nga*, Le-Tien Thuong, Mai Linh
School of Computer Science and Engineering, International University VNUHCM, Vietnam
Department of Electrical & Electronics Engineering, University of Technology-VNUHCM, Vietnam
Corresponding author :
April 29, 2016 October 10, 2016 December 6, 2016


In this paper, we consider a typical health care system via the help of Wireless Sensor Network (WSN) for wireless patient tracking. The wireless patient tracking module of this system performs localization out of samples of Received Signal Strength (RSS) variations and tracking through a Particle Filter (PF) for WSN assisted by multiple transmitpower information. We propose a modified PF, Kullback-Leibler Distance (KLD)-resampling PF, to ameliorate the effect of RSS variations by generating a sample set near the high-likelihood region for improving the wireless patient tracking. The key idea of this method is to approximate a discrete distribution with an upper bound error on the KLD for reducing both location error and the number of particles used. To determine this bound error, an optimal algorithm is proposed based on the maximum gap error between the proposal and Sampling Important Resampling (SIR) algorithms. By setting up these values, a number of simulations using the health care system’s data sets which contains the real RSSI measurements to evaluate the location error in term of various power levels and density nodes for all methods. Finally, we point out the effect of different power levels vs. different density nodes for the wireless patient tracking.



    The current challenge in health care system is to estimate physical parameters associated to the patient (Mann et al., 2004) because of increasingly aging population combined with nursing staff shortages and decreasingly the capacity of hospital (Ko et al., 2010). Wireless patient tracking/monitoring provides flexible and powerful patient surveillance through wearable devices anytime and anywhere by Ren et al. (2010).

    The core of health care system is an automatic supervision of patients within the nursing institute based on the lightweight system via Wireless Sensor Network (WSN). Patient’s supervision includes two major services: tracking the location of patients and monitoring the patient’s status. In the first service, patients need urgent staff intervention, the exact knowledge of their location plays an important factor to reach them in a short period of time. On the other hand, the last service, the status of patients in critical situation must be continuously available to the medical staff, when patients can roam around the hosting institution’s premises. Many kinds of information on the patient’s status are continually collected and automatically detected such as moving characteristics, heartbeat, etc.

    Several researchers such as Redondi et al. (2010), Lim et al. (2010) and Evennou et al. (2006) introduced the LAURA (LocAlization and Ubiquitous monitoRing of pAtients) system, with Personal Monitoring System (PMS) and Personal Localization and Tracking System (PLTS) to monitor/track the wireless patient localization based on Received Signal Strength (RSS). In case of moving targets, Lim et al. (2010) introduced gradient descent to determine the wireless patient position when the distances from the patient to the anchor nodes computed via Signalto- Distance Mapping (SDM).

    Most research in WSN had emphasized the use of PFs for solving many problems. Nguyen et al. (2013) derived a new method for evaluation of the number of particles in WSN. Li et al. (2013 and 2015) introduce KLD-resampling for solving problems in WSN under the benchmark model of maneuvering in a two-dimensional plane. Ly-Tu et al. (2014) introduced the ability of KLDresampling PF to evaluate the number of particles used for non-linear problems.

    Recently, by incorporating these PF methods and gradient descent to avoid wall crossing under known environment’s map for improving wireless patient tracking via the help of WSN. Redondi et al. (2010) and Evennou et al. (2006) introduced PF based on Sampling Importance Resampling (SIR) to enhance the quality of tracking. While Ly-Tu et al. (2015) introduced the KLDresampling PF for wireless patient tracking by the given upper bound error with fixed probability so as to create a sample set near the high-likelihood region. By setting up a specific upper bound error for the resampling parameter, we found that this technique further reduces the location error of patient but maintains the proper KLD-resampling (Li et al., 2015) for health care systems because the variation of RSS measurement values is the diminished. The upper bound error of the resampling parameter (Ly-Tu et al., 2015) is only a special case of Remark in Section 3.

    The main contributions of this paper are summarized as follows. First, we extend the underlying system (Ly-Tu et al., 2015) by proposing an algorithm to obtain the upper bound error of resampling parameter based on the maximum gap error between the proposed algorithm and SIR method for various power levels and density nodes. Second, given the strong contribution of the selected upper bound values of resampling parameter and weight function method, we use them in the proposal for assessing the impact of power levels of self RSSI vs. different density nodes to track wireless patient location.

    The paper is organized as follows. Introduction to a problem solving framework such as the patient localization based on gradient descent method, the improved localization via PF, and error criteria are given in Section 2. All related schemes, namely, SIR and our method are presented in Section 3. All experimental results based on MATLAB for wireless patient tracking are shown in Section 4. Finally, we conclude the paper in Section 5.


    As has been said, we focus on the wireless patient tracking in health care system. Lim et al. (2010) discussed about the problem context and its location data requirements are given in the first section. The wireless patient localization by gradient descent method is shown in the next section. Final, the last section describes PFs serving to enhance the wireless patient tracking by Evennou et al. (2006) and Rendodi et al. (2010).


    Evennou et al. (2006) introduced a typical health care system, as shown in Figure 1, to assemble the required data from the patients and distribute the same data via a central controller. This system contains the personal monitoring system (PMS), the personal localization and tracking system and a lightweight and flexible network architecture (NA). The PMS is composed of smallwearable devices geared with sensing and wireless transmission capabilities to monitor patient’s status and transmit the comprehensible information to the central controller via the help of the NA.

    The accuracy of localization calculates on the site of the anchor nodes and on their density, for instance, an average localization error of just a few meters is evaluated via real experiments.

    Regardless of the fundamental network infrastructure, localization is performed adopting analogous processing methods. Under fixed and known the placement of the anchor nodes, RSS measurements are gathered between pairs of nodes which related to the inter-node distances. Indeed, a RSS indicator (RSSI), which addressed by non-parametric or parametric methods, is supplied by off-the-shelf devices without the demand of ded- icated hardware.

    The first method, namely non-parametric method or fingerprinting which is introduced by Lim et al. (2010) and Evennou et al. (2006), is described to avoid the explicit modeling of the RSSI-distance relationship. The collection of fingerprint databases, which is implemented during a training phase, has many fingerprint vectors. Each vector consists of the RSSI values measured between a client node at a specific site and the anchor node. A client node forms a vector of RSSI measurements between itself and each anchor node and it compares this vector with the pre-computed fingerprints in the database. Final, its position is estimated based on the placement of the most similar fingerprints. This method suffers from a non-negligible setup cost which is related to the construction of the fingerprint database. Furthermore, it is not robust in face of changes of the environment due to invalid the affection of RF propagation and pre-computed fingerprints.

    The second method, which is called parametric method, explicitly assumes the propagation model’s knowledge that involves the RSSI metrics S (expressed in dBm units) and the inter-node distance d. Patwari et al. (2005) introduced an effortless dependency as the following model

    S = S 0 10 α log 10 d d 0 + v ,

    where S0 is the RSSI metrics measured between two nodes, d0 meters apart, and α indicates the factor in path loss for indoor environments. The effects of shadowfading in complex multipath environments is represented by the noise term v with a Gaussian random variable N ( 0 , σ v 2 ) , whereas the value of standard deviation σv depends on the characteristics of the specific environment. Barsocchi et al. (2009) introduced this method to solve the localization problem by computing the position that maximizes the likelihood with respect to the model in (1). The extension of this model is also introduced by Barsocchi et al. (2009) to take into account the signal attenuation because of walls. The advantage of this approach is to restrict the system-setup cost by estimating based on RSSI measurements between anchor nodes. Meanwhile, Lim et al. (2010) also introduced a similar approach via applied the propagation model in (1) under implicit assumptions, but the distance between a client and anchor node is modeled as a linear combination of the RSSI measurements between the client and all the nodes. The construction of these RSSI measurements based on this way is called the RSSI-to-distance linear mapping. The disadvantage of this approach is to setup a zero-configuration.

    Because the noise influents RSSI measurements, temporal averaging is used over a time window. The length of the window enables to define a trade-off between the RSSI measurement’s accuracy and the system delay. As a rule of thumb, every 200ms each node is collected and 3-5 RSSI measurements are updated.

    In case of the moving patient, the position of the wireless patient is tracked over time so as to improve the quality of localization. To this end, Evennou et al. (2006) applied a simply PF method which enables to merge the RSSI measurements together with the typical movement patterns of humans in indoor environments. This way can avoid wall crossing because of incorporating this a-priori information and PF method.

    2.2.Indoor Localization Via Gradient Descent

    Let us denote xi ∈ R2 , i = 1, …, M be the position of an anchor node, while M indicates M anchor nodes. Each anchor node obtains a vector s i =  [ s i1 ,  s i2 ,  ,  s iM ] T of RSSI measurements exchanging data with the other nodes, where sij is the RSSI relative to the signal emitted by the j-th anchor node. Then, let us define S=  [ s 1 ,  s 2 ,  ,  s M ] R MxM be the overall RSSI measurements among all anchor node pairs. In general, S is not symmetric in radio links; diagonal entries of S contain the self-RSSI values.

    Let us define d i  =   [ d i1 ,  d i2 ,  ,  d iM ] T with dij is the Euclidean distance between anchors i and j, in which case of the distance vector, and the corresponding matrix D =   [ d 1 ,  d 2 ,  ,  d M ] is symmetric and has zero diagonal entries. Note that these are the solely parameters that require to be manually obtained during the system calibration phase, as they depend on the specific hardware. A linear relationship between the RSSI measurements and the logarithm of the inter-node distance by Lim et al. (2010) is followed(2)

    log ( D ) = TS

    where T∈ℝMxM is the signal-to-distance mapping and each row vector is represented as a linear combination of the columns of S, weighted by the elements of the i-th row t i T . The matrix T is calculated by least square method as(3)

    T = log ( D ) S T ( SS T ) 1 .

    Once the signal-to-distance mapping is determined, the client node’s localization proceeds collecting RSSI measurements between itself and its adjacency anchor nodes in the vector s ^ . Then, the corresponding distance vector d ^ can be computed as

    d ^ = exp ( T s ^ )

    The location based on the obtained distance vector d ^ is estimated owing to apply a gradient descent algorithm by Lim et al. (2010). It is exploited in LAURA system to initialize tracking.

    2.3.Improving Localization Via Particle Filtering

    The paper focuses on improving the wireless patient tracking in LAURA PLT (Personal Localization and Tracking) module by Redondi et al. (2010 and 2013) to smooth down the functions of the RSSI variations due to PF.

    In case of static conditions and fixed client nodes, the localization algorithm described above achieves slightly perfect when robust RSSI measurements can be obtained by means of the time averaging. On the other hand, when the node is attached to the patient’s body, in which case of moving client nodes, severe short-time RSSI fluctuations are observed. As such, the accuracy of the estimated location decreases, and artificial motion discontinuities are detected.

    The localization accuracy is improved with a help of incorporating a-priori knowledge about the moving node and the geometry of the environment. In this system, a PF to track the position of the client node over time is used. At each time instant t, the state of the node is represented by the vector z ( t ) = [ x ( t ) T ,  v ( t ) T ] T , where x(t)∈R2 indicates the position of the node and v(t)∈ℝ2 its velocity.

    The sequence of previous states and all available observations are given. It is defined z p ( t ) = [ x p ( t ) T , v p ( t ) T ] T PF estimates the posteriori probability density function of the state z at time t with k = 1,…, N and wp(t) is a weighted set of particles. The state of the target node is then point-wise estimated as:

    { x ^ ( t ) = ( 1 α x ) ( x ^ ( t 1 ) + v ^ ( t 1 ) Δ T ) + α x w p ( t ) x ^ ( t ) , v ^ ( t ) = ( 1 α v ) v ^ ( t 1 ) + α v ( x ^ ( t ) x ^ ( t 1 ) ) ,

    where ΔT is the sampling period, and αx , αv are the adaption rates which are tuned experimentally based on the node dynamics.

    For the problem at hand, PF alternates the execution of two steps: prediction and update. Firstly, these assumptions of system will be discussed. Then, prediction step will be shown in the next part. Finally, how to update process will be mentioned in the last.


    The location of a target node x ^ ( 0 ) is the first thing at startup which can be computed by using the estimated vector d ^ through a gradient descent method which minimizes the following cost function:

    x ^ ( 0 ) = arg min x 1 2 i = 1 M θ i ( x-x i 2 d ^ i ) 2 ,

    where d ^ i is the estimated distance between the target node and the i-th anchor and θi is a weight factor computed as(7)

    θ i = d ^ i 2 i=1 M d ^ i 2

    By differentiating equation (6) with respect to x it is possible to obtain an iterative equation used to update the estimated x ^ :(8)

    x ^ ( k+1 ) = x ^ ( k ) + α i = 1 M θ i ( 1 d ^ i x ^ ( k ) x i 2 ) ( x ^ ( k ) x i ) ,

    with α is set to 0.1; the initial estimate x ^ ( 0 ) is equal to the location of the nearest anchor node, and then zp (0)=0[xp,(0,0)T  Tp .

    Secondly, since the initialization of the particle filter might potentially locate the target outside the building, the geometry leverage’s knowledge of the environment to constrain the initialization point to lie inside the building.


    At each step, the new state vector for the p-th particle is computed. The following dynamic model constructed upon kinematic equations is employed.(9)

    { x p ( t ) = x p ( t - 1 )  + v p ( t - 1 ) Δ T + ξ x v p ( t ) = v p ( t -​​  1 ) + ξ v ,

    with ξx and ξv are samples drawn from zero-mean Gaussian random variables with variance, respectively, σ ξ x 2 and σ ξ v 2 , which provide the system with a diversity of hypotheses. If the geometry of the environment is known, it can be incorporated in the dynamic model, e.g. by assigning wp = 0 to those particles that crossed a wall.


    After the prediction step, the weight of the p-th particle is updated based on the RSSI measurements collected at time t. Firstly, the vector d ^ as descried in (4) is obtained. Then, we assign a weight computed as

    w p ( t ) = w p ( t - 1 ) exp ( 1 2 i = 1 M w i ( x p ( t ) x i 2 d ^ i ) 2 ) ,

    where d ^ i is the estimated distance, and x p ( t ) x i 2 is the current distance of the particle from the i-th anchor.

    In order to avoid the critical situation, the weight distribution of the particles in (10) is tested. When all particles are trapped within a room while the client node had moved outside.(11)

    p = 1 N w p ( t ) > τ 1

    where τ1 is a threshold value that we set equal to 10-5 in our experiments. If the test fails, particles will not correct tracking the client node. Thus, the PF is re-initialized as described below. Otherwise, the particle weights are normalized such that p=1 N w p ( t ) = 1 and the client node’s position is determined based on (5).

    Working with PFs as customary, a resampling step will perform if the weight distribution is severely skewed, such that has just a few particles with non-negligible weight. In this paper, two kinds of resampling steps, namely SIR and KLD-resampling, will be discussed in the next step.

    2.4.Error Criteria

    Let us define Er to be the error for measuring the performance of the residuals which are considered the difference between the observed and predicted values.(12)

    E r = i = 1 K ( x ^ i x i )

    where x ^ i and xi are respectively the predicted location and the observed location at the i-th observation, and K is the number of observations.


    In this section, SIR PF is reviewed and our proposal (KLD-resampling PF) is presented.

    3.1.SIR Algorithm

    The concept of the auxiliary PF has been introduced by a number of researchers such as Redondi et al. (2010 and 2013) and Arulampalam et al. (2002). A PF is also known as the bootstrap filter, Monte Carlo technique, condensation algorithm, interacting particle approximations and survival of the fittest. The key idea is to represent the required posterior density function by a set of random samples (particles) with associated primary weights, and to compute the estimates based on the samples and primary weights. As the number of samples becomes very large, the Monte Carlo characterization is the closest equivalent for representation of the posterior probability function, and the solution approaches the optimal Bayesian. Schön et al. (2010) introduced the sequential importance sampling (SIS) algorithm, shown in Algorithm 1 for the PF, includes a resampling step at each instant as described in detail in the reference. The SIS algorithm uses the important density, which is a proposed density to represent another one that cannot be exactly computed, that is the sought posterior density in the present case. Hence, samples are drawn from the important density instead of the actual density. The degeneracy phenomenon is known as a common problem with the SIS PF, which all but one particle has negligible primary weight after a few states. Thus, a large computational effort is devoted to updating particles whose contribution to the approximation of the posterior density function is almost zero. This problem can be overcome by increasing the number of particles or more efficiently by approximately selecting the important density. In addition, the use of resampling technique is recommended to avoid the degeneracy of the particles by Arulampalam et al. (2002), as Algorithm 1. The pseudo-code of SIR PF is shown in Algorithm 2.

    3.2.Proposal Algorithm

    In this subsection, we discuss two problems. First, we describe our KLD-resampling for applied PF. Second, we present a way to determine the optimal bound values of resampling parameter for our method, as shown in Algorithm 4. Because our proposal which relates to KLDsampling, thus we briefly overview this method.

    Fox et al. (2003) introduced KLD-sampling to determine the number of samples such that, with probability 1−δ , the error between the true posterior and the samplebased approximation is less than ε . KLD is used to show how to determine the number of samples so that distance between the sample-based maximum likelihood estimate and the true posterior does not exceed a prespecified threshold.

    The KLD between the proposal (q) and (p) distributions is defined in discrete form as(13)

    d KL ( p q ) = x p ( x ) log ( p ( x ) q ( x ) ) = x W ( x ) q ( x ) log W ( x )

    with W(x) = p(x) / q(x). The required number Nr of samples is determined as follows

    N r = 1 2 ε χ l 1 , 1 δ 2

    where l is the number of bins with support, the quantizes of Chi-square distribution is computed as(15)

    P ( χ l 1 2 χ l 1 , 1 δ 2 )  =  1  -  δ

    Based on Wilson-Hilferty transformation by Johnson and Kotz (1994) to compute the approximation of χ l 1 , 1 δ 2 in (14) is expressed as

    N r  =  l 1 2 ε ( 1 2 9 ( l 1 ) + 2 9 ( l 1 ) z 1 δ ) 3

    where z1−δ is the upper quartile of the standard normal distribution.

    The undesirable results above are caused by the number of particles needed to approximate a discrete distribution with an upper bound ε on the KLD. In other words, the KLD-sampling method leads to statistical bounds of the approximation quality of samples that are actually drawn from the proposal rather than the true posterior distributions. Therefore, the mismatch between the true and the proposal distribution is ignored,


    To avoid the mismatch from KLD-sampling, the total number of required particles in (16) is applied in the resampling process better than sampling one. Li et al. (2013 and 2015) introduced this method which divides the particles of the posterior distribution into bins and count the number l of bins which at least one particle is resampled to determine the total number of particles to resample. Therefore, the required number Nr in (16) is rewritten as(17)

    N r, re  =min ( N max , c e i l ( l 1 2 ε ( 1 2 9 ( l 1 ) + 2 9 ( l 1 ) z 1 δ ) 3 ) )

    The Pseudo-code of proposal is shown in Algorithm 3.

    3.2.2.Finding the upper bound error

    We propose an algorithm to find the upper bound error as shown in Algorithm 4. Let us denote εi to be the upper bound error of i th, E r ε i Pr o and ErSIR which are respectively the error of proposal at εi and that of SIR; and εi* (line 11) is the set of the upper bound value that fulfils the condition of Remark (line 10). Let us define Δ E r = [ Δ E r s 1 , , Δ E R s Q ] with Δ E r ε i = E R SIR E r ε i Pr o is the gap of error value between SIR and proposal at εi . Algorithm 4 provides the optimal bound value of proposal.

    Remark: if ΔEr > 0 (lines 10-13), then a value εopt (line 16) exists to maximize the function ΔEr.


    We conducted a series of simulations via the data sets, which are used from publicly available by ANT Lab, to determine the optimal bound value of resampling parameter in (16) based on Algorithm 4. Next, the performance of tracking error, various power levels of self- RSSI by Hongliang et al. (2009) and the gap error for all methods are evaluated under four density nodes by Redondi et al. (2010 and 2013) and Ly-Tu et al. (2015). The network architecture described by Lim et al. (2010) which consists of MicaZ and IRIS motes operated by the TinyOS.

    All simulation results are run on PC Core i5-2400 @3.10GHz, 4.00GB RAM and MATLAB 2012a ( Furthermore, all parameters assignment for this system is given in Table 1.

    The localization system performance integrated PF is evaluated to describe in Section 2. A client node carried by a person is tracked that followed a known path in an indoor area of around 250m2. The node was placed at each test point for 30 seconds, so that 30 RSSI packets were sent to the PAN coordinator. To synchronize the RSSI message reception with the actual position of the patient node, the path was sampled each 90cm. For every RSSI packet, the node’s location is estimated using the algorithm which described in Section 2.

    The test was repeated for different density nodes and different power levels. To discard the bias related to a specific anchor node, the choosing estimation in random subset of anchor nodes was repeated for each density node.

    The results in Table 2 show the optimal bound values with various power levels of self-RSSI in term of four density nodes for our proposal. For example, at power level of self-RSSI (-15dBm), the bound values for four various density nodes are respectively about 0.3, 0.95, 0.75 and 0.9 (see column 7 in Table 2). First, we set up these values for our proposed method in Figure 3 to compare the cumulative density function of the tracking error for Gradient descent, SIR with 200 particles, SIR with 400 particles and proposal with 200 particles is shown in Figure 2. Final, all these values in Table 2 are applied to evaluate the performance of gap error vs. four various density nodes for all these methods and examine thoroughly the impact of density nodes on the whole power levels for our proposal in Figure 3 and Figure 4, respectively.

    The comparison of the cumulative density function of the tracking error for Gradient descent, SIR with 200 particles, SIR with 400 particles and proposal with 200 particles in case of power level (-15dBm) for four density nodes (5nodes, 10nodes, 15nodes and 20nodes) is illustrated in Figure 2. It is observed that the tracking error of each density node, at errors in 80% of the cases (Pr < 80%), for Gradient descent method achieves below 5m, around 4m and below 2.5m are respectively in Figure 2a, Figure 2b, and Figure 2cd. As enabling the PFs, the localization accuracy of our proposed is better than that of SIRs and Gradient descent methods. For example, at Pr = 0.8 in Figure 2c, it verifies that the curve of the localization accuracy for our proposal with 200 particles enhances respectively about 0.11 and around 0.37 that of SIR with 200 particles and SIR with 400 particles. It also verifies that our proposal not only reduces the number of used particles because of the selected upper bound error with fixed probability of resampling parameter, but also improves the average patient localization error compared with SIRs PF in clinical environments. Specially, as increasing the number of nodes our curve can improve the performance of the tracking error compared with others and satisfy the accuracy error from 3m to 5m by Farid et al. (2013).

    We further study the relationship between the gap error and various power levels for all methods under different density nodes is shown in Figure 3. At 5 nodes as illustrated in Figure 3a, overall particles used of proposal is almost equal to the number of particles used of SIR, it confirms the statistical gaps error of the proposal and SIR are respectively about 0.8, 0.38, 0.22, 0.2, 0.04, 0.26, and 0.36 when changing power levels of self-RSSI from -25dBm to 0dBm, respectively. In addition, as increasing the number of particles for SIR to twice, our method performs slightly about 0.03, 0.11, 0.04, 0.1, 0.06, 0.05, and 0.19 when compared this method. Similar to Figure 3a, when changing the density nodes from 0.02 nodes/m2 to 0.08 nodes/m2, the proposal also performs with traditional others.

    Before reaching a conclusion, it is important to examine thoroughly the impact various power levels on the tracking error for our proposal from 0.02nodes/m2 to 0.08nodes/m2 as illustrated in Figure 4. First, Figure 4d shows the performance of tracking error, the all curves of our proposal converge and reach below 2m at Pr = 80%. On the other hand, with 20 nodes or 0.08nodes/m2, the proposal can apply to any power levels or can use with a single transmit-power information. Moreover, the proposal in case of 15 nodes fulfils hard accuracy criteria only two power levels which are 0dBm and -3dBm. Figure 4a and Figure 4b illustrate to find power levels for improving tracking error from 3m to 5m in others density nodes.


    In this study, a KLD-resampling under the optimal bound error of resampling parameter is introduced to smooth down the fluctuations of RSS variations. Ultimately, the wireless patient tracking is enhanced using WSN with a PF system structure that is designed by these parameters. Furthermore, in using these bound error values, this technique reduces the number of particles used and produces accurate patient location in real environment. We conduct a large set of simulations to study the tracking error, the relationship between the gap error and various power levels, and the effect of different density nodes under the whole power levels. For future work, we will assess PFs based on an efficient KLD-resampling algorithm on board. Future works, we will analyse the problem under probability to show the great contribution in the slope of the mean number of particles used.


    This research is funded by Vietnam National University HoChiMinh City (VNU-HCM) under grant number C2015-28-03. The authors would like to thank the International University of the VNU-HCM for its support and the anonymous reviewers for their valuable comments.



    Functional architecture of system (Evennou et al., 2006).


    Pseudo-code of SIS


    Pseudo-code of SIR PF


    Pseudo-code of KLD-resampling (Proposal)


    Finding the upper bound error


    Cumulative density function of the tracking error in case of power level (-15dBm) (a) 5nodes (b) 10nodes (c) 15nodes (d) 20nodes.


    Gap error vs. various power levels of self-RSSI for all methods (a) 5nodes (b) 10nodes (c) 15nodes (d) 20nodes.


    Cumulative density function of the tracking error vs. various power levels for our proposal (a) 5nodes (b) 10nodes (c) 15nodes (d) 20nodes.


    Parameters assignment of LAURA system

    Performance of finding the optimal bound algorithm


    1. ANT L (2015) LAURA-Localization and Ubiquitous Monitoring of Patients for Health Care Support, available from: Accessed: Aug. 30
    2. Arulampalam M S , Maskell S , Gordon N , Clapp T (2002) A tutorial on particle filters for online nonlinear/non-gaussian bayesian tracking , IEEE Transactions on Signal Processing, Vol.50 (2) ; pp.174-188
    3. Barsocchi P , Lenzi S , Chessa S , Giunta G (2009) Virtual calibration for rssi-based indoor localization, with ieee 802154 , IEEE International Conference on Communications, ; pp.1-5
    4. Evennou F , Marx F (2006) Advanced integration of wifi and inertial navigation systems for indoor mobile positioning , EURASIP Journal on Applied Signal Processing, ; pp. 164-164
    5. Farid Z , Nordin R , Ismail M (2013) Recent advances in wireless indoor localization techniques and system , Journal of Computer Networks and Communications, ; pp.1-12
    6. Fox D (2003) Adapting the sample size in particle filters through KLD-sampling , The International Journal of Robotics Research, Vol.22 (12) ; pp.985-1003
    7. Hongliang R , Max Q-H (2009) Power adaptive localization algorithm for wireless sensor networks using particle filter , IEEE Trans. Vehicular Technology, Vol.58 (5) ; pp.2498-2508
    8. Johnson N , Kotz S , Balakrishnan N (1994) Continuous Unvariate Distribution, John Wiley & Sons, Publisers,
    9. Ko J , Gao T , Rothman R , Terzis A (2010) Wireless sensing systems in clinical environments: Improving the efficiency of the patient monitoring process , Engineering in Medicine and Biology Magazine, Vol.29 (2) ; pp.103-109
    10. Li T (2015) Resampling methods for particle filtering: classification, implementation, and strategies , IEEE Signal Processing Magazine, Vol.32 (2) ; pp.70-86
    11. Li T , Sun S , Sattar T P (2013) Adapting sample size in particle filters through KLD-resampling , Electronics Letters, Vol.49 (12) ; pp.1-2
    12. Lim H , Kung L-C , Hou J C , Luo H (2010) Zero-configuration indoor localization over ieee 80211 wireless infrastructure , Wireless Networks, Vol.16 (2) ; pp.405-420
    13. Ly-Tu N , Le-Tien T , Tran-Hoang-Thi Oanh , Huynh Kha Tu , Mai L (2014) Performance of sampling/resampling-based particle filters applied to non-linear problems , REV Journal on Electronics and Communications, Vol.4 (2-4) ; pp.75-84
    14. Ly-Tu N , Le-Tien Thuong , Mai L (2015) Kullback-leibler distance resampling based particle filter applied to health care system, in Proceedings on the 16th ASIA Pacific Industrial Engineering and Management Systems Conference , APIEMS, ; pp.301-307
    15. Mann W (2004) The aging population and its needs , Pervasive Computing, Vol.3 (2) ; pp.12-14
    16. Nguyen T D , Trinh H T , Nguyen N T , Nguyen T T , Koichiro W (2013) A new evaluation of particle filter algorithm and apply it to the wirelesssensor network , International Conference on Computing Management and Telecommunicaions, ; pp.169-174
    17. Patwari N , Ash J N , Kyperountas S , Hero A O , Moses R L , Correal N S (2005) Locating the nodes: Cooperative localization in wireless sensor networks , IEEE Signal Processing Magazine, Vol.22 (4) ; pp.54-69
    18. Redondi A , Chirico M , Borsani L , Cesana M , Tagliasacchhi M (2013) An integrated system based on wireless sensor networks monitoring, location and tracking , Ad Hoc Networks Elsevier, Vol.11 (1) ; pp.39-53
    19. Redondi A , Tagliasacchi M , Cesana M , Borsani L , Tarrio P , Salice F (2010) LAURA-Loca Alization and ubiquitous monitoring of pAtients for health care support, Symposium on Personal , Indoor and Mobile Radio Communications Workshops, ; pp.218-222
    20. Ren Y , Pazzi R , Boukerche A (2010) Monitoring patients via a secure and mobile health care system , Wireless Communications, Vol.17 (1) ; pp.59-65
    21. Schön T B (2010) Solving nonlinear state estimation problems using particle filters-an engineering perspective, Technical Report, Automatic Control at Linkopings Universitet