Open AccessArticle

Enhancing the Accuracy and Robustness of a Compressive Sensing Based Device-Free Localization by Exploiting Channel Diversity

Dongping Yu

Yan Guo

^*,

Ning Li

and

Xiaoqin Yang

College of Communications Engineering, Army Engineering University of PLA, Nanjing 210007, China

Author to whom correspondence should be addressed.

Sensors 2019, 19(8), 1828; https://doi.org/10.3390/s19081828

Submission received: 8 March 2019 / Revised: 13 April 2019 / Accepted: 15 April 2019 / Published: 17 April 2019

(This article belongs to the Special Issue Multi-Sensor Systems for Positioning and Navigation)

Download

Browse Figures

Figure 1
An example of the CS-based DFL. "> Figure 2
Spatial impact area of a wireless link. "> Figure 3
Architecture of the ComDec method. "> Figure 4
Graphical model for the joint sparse recovery. "> Figure 5
Factor graph representation of the joint PDF (22). "> Figure 6
Simulation flowchart. "> Figure 7
The performance of ComDec when channel number varies from 1 to 25. "> Figure 8
The values of the dictionary atoms: (a) <math display="inline"><semantics> <mrow> <mi mathvariant="bold-italic">ϕ</mi> <msubsup> <mo>.</mo> <mrow> <mn>10</mn> </mrow> <mn>5</mn> </msubsup> </mrow> </semantics></math>; (b) <math display="inline"><semantics> <msubsup> <mi mathvariant="bold-italic">ϕ</mi> <mrow> <mn>42</mn> </mrow> <mn>12</mn> </msubsup> </semantics></math>, when <math display="inline"><semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics></math>, 5, and 10. "> Figure 9
Effect of environmental changes on the average localization error <math display="inline"><semantics> <mrow> <mo stretchy="false">(</mo> <mi>A</mi> <mi>v</mi> <mi>g</mi> <mo>.</mo> <mi>E</mi> <mi>r</mi> <mi>r</mi> <mi>o</mi> <mi>r</mi> <mo stretchy="false">)</mo> </mrow> </semantics></math> of multiple DFL approaches. "> Figure 10
Comparison of the original location vector and the reconstructed sparse vectors of different DFL approaches. (a) the original location vector; The reconstructed sparse vector corresponding to (b) E-HIPA; (c) LCS; (d) DR-DFL; (e) ComDec (<math display="inline"><semantics> <mrow> <mi>F</mi> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math>); (f) ComDec (<math display="inline"><semantics> <mrow> <mi>F</mi> <mo>=</mo> <mn>15</mn> </mrow> </semantics></math>). "> Figure 11
Comparison of the original target positions and the estimated target positions of different DFL approaches. "> Figure 12
Average localization error with different SNR. "> Figure 13
Average localization error with different K. ">

Versions Notes

Abstract

As an emerging and promising technique, device-free localization (DFL) estimates target positions by analyzing their shadowing effects. Most existing compressive sensing (CS)-based DFL methods use the changes of received signal strength (RSS) to approximate the shadowing effects. However, in changing environments, RSS readings are vulnerable to environmental dynamics. The deviation between runtime RSS variations and the data in a fixed dictionary can significantly deteriorate the performance of DFL. In this paper, we introduce ComDec, a novel CS-based DFL method using channel state information (CSI) to enhance localization accuracy and robustness. To exploit the channel diversity of CSI measurements, the DFL problem is formulated as a joint sparse recovery problem that recovers multiple sparse vectors with common support. To solve this problem, we develop a joint sparse recovery algorithm under the variational Bayesian inference framework. In this algorithm, dictionaries are parameterized based on the saddle surface model. To adapt to the environmental changes and different channel characteristics, dictionary parameters are modelled as tunable parameters. Simulation results verified the superior performance of ComDec as compared with other state-of-the-art CS-based DFL methods.

Keywords:

changing environment; channel diversity; variational Bayesian inference; joint sparse recovery; device-free localization

1. Introduction

As a main piece of the context information, location information is essential for the location-based services (LBS) (all the used abbreviations are explained at the end of the article) in many context-aware applications. Nowadays, localization techniques have been an active field in pervasive and ubiquitous computing. Generally, the target localization methods can be classified as: device-free [1,2,3] and device-based [4,5,6]. Device-based methods need objects to attach assistant wireless devices for signal transmitting or receiving. However, in many applications, such as emergency rescue, intruder detection, and smart homes, it is difficult to attach objects with any transceivers. In this case, the device-based method will be infeasible.

Device-free localization (DFL) methods do not need to attach objects with any assistant devices. This approach has become a crucial component in many context-aware applications. As the transceiver-free objects cannot be directly perceived, DFL methods estimate their positions by analyzing their influences on surrounding radio environments. Among current DFL methods, some rely on specialized hardware, such as radar-based methods [7], camera-based methods [8], and infrared-based methods [9]. These methods need an extensive deployment of dedicated devices in the monitoring area and may involve privacy issues. In recent years, DFL methods based on existing infrastructures (e.g., WSNs/WiFi) have attracted a lot of research interests. They do not require dedicated hardware and only make use of the measurement information available from the already deployed wireless devices.

As a common link measurement, the received signal strength (RSS) is nearly ubiquitously available from the standard radio transceivers [10,11]. The RSS measurement information has been widely used in DFL. When targets located at different positions, the RSS readings will be different. To exploit the location dependence of RSS, RSS readings on multiple wireless links are recorded before and after targets entering into the monitoring area. To establish the relationship between RSS variations and target locations, the fingerprinting technique has been introduced [12]. Typically, the fingerprinting-based DFL consists of an online and offline phase. In the offline training phase, by gathering the RSS variations caused by a target at every possible location, a radio map can be built. In the online locating phase, target locations are estimated by matching current RSS variations with the fingerprints in radio map. However, as a major drawback of the fingerprinting-based DFL methods, gathering fingerprints to build a radio map is labor-intensive and time-consuming.

As an alternative to fingerprinting-based methods, model-based methods have been widely used in DFL. They map RSS variations to target positions by theoretical or empirical shadowing models. Based on these models, the dictionary (a.k.a. radio map) can be built without site survey. Unfortunately, in changing environments, RSS is extremely sensitive to environmental dynamics, such as temperature, humidity, electromagnetic characteristics, and pedestrians around [13]. In this case, shadowing effects cannot be well approximated by RSS variations. The spatial and temporal environmental dynamics may result in mismatches between runtime RSS measurements and the data in a fixed dictionary. In [14], a dictionary refinement algorithm is introduced, which alleviates the dictionary mismatches by iteratively refining the model-based dictionary. However, the real-time optimization of dictionary parameters leads to a high computational complexity. To avoid this, in ComDec, the dictionary parameters are modelled as tunable parameters to adapt to the changes of environment, and we do not need to explicitly estimate their values.

Recent years, in target localization, the channel state information (CSI) [15] has been exploited, which is a measurement information from PHY layer. As a fine-grained value, CSI depicts the channel quality on multiple orthogonal subcarriers. In wireless environments, due to the frequency-independent attenuation and frequency-selective fading, each channel will exhibit a unique amplitude and phase. Hence, the dictionaries corresponding to different channels are different. Over the past few years, in many device-based localization methods, CSI measurements has been leveraged [16]. Recently, more efforts have been paid in device-free technique. In [17], the fine-grained subcarrier information is exploited based on the multitask Bayesian compressive sensing (MBCS). However, the DFL method proposed in [17] do not provide a solution to the dictionary mismatch problem caused by environmental dynamics.

Compared with RSS measurements, CSI measurements are more robust and suitable for being utilized in DFL. In this paper, a novel ComDec method is proposed. It can enhance the accuracy and robustness of the CS-based DFL in changing environments. In ComDec, to enrich the measurement information, CSI measurements are collected from multiple frequency bands. Moreover, to reduce wireless sensor nodes, the compressive sensing (CS) theory [18] is applied in ComDec by taking advantage of the spatial sparsity of target localization. The CS-based DFL problem in multi-channel scenario is formulated as a joint sparse recovery problem. Furthermore, to bypass the dictionary training and retraining works, the dictionaries are built based on the saddle surface model [19]. We treat dictionary parameters as tunable parameters to adapt the changes of environment and different channel characteristics. Afterwards, to recover the sparse vectors of multiple channels, we develop a joint sparse recovery algorithm under the variational Bayesian inference framework [20]. The main contributions of this work are as follows:

To enhance the localization accuracy and robustness of CS-based DFL, a novel ComDec method is proposed, which leverages the channel diversity of CSI measurements. In ComDec, the CS-based DFL problem is extended to multi-channel scenario. It is formulated as a joint sparse recovery problem that recovers multiple jointly sparse vectors over two known dictionaries.
To simultaneously recover the jointly sparse vectors, we develop a novel joint sparse recovery algorithm. The joint sparsity of the sparse vectors is induced by a novel two-layer hierarchical prior model. Then, the common support set of the sparse vectors is estimated by inferring the posteriors of the hidden variables that defined in the proposed prior model.
To mitigate the influence of environmental dynamics in changing environments, the dictionary parameters with respect to multiple channels are modelled as tunable parameters to adapt the environmental changes and different channel characteristics. In this way, the dictionary mismatch problem can be solved without the need of explicitly estimating the dictionary parameters.
To reduce the computational complexity, we introduce four methods in the proposed joint sparse recovery algorithm. Among them, the grid pruning method can improve the convergence speed of the proposed joint sparse recovery algorithm.

The remainder of this paper is organized as follows. Section 2 presents an overview of the related works on multi-target DFL. Section 3 gives the signal model and formulates the CS-based DFL problem as a joint sparse recovery problem. Section 4 proposes a novel joint sparse recovery algorithm to recover the jointly sparse vectors. Section 5 validates the proposed ComDec method with extensive numerical simulations. Finally, Section 6 concludes this paper.

2. Related Work

Beginning with the initial papers of Youssef et al. [1] and Zhang et al. [2], numerous research works on DFL have been carried out [17,21]. Typically, there are four types of DFL methods: (1) geometry-based methods, (2) fingerprinting-based methods, (3) radio tomographic imaging (RTI)-based methods, and (4) CS-based methods. The geometry-based methods estimate target locations based on the geometry information of wireless links. They have a restriction on target spacing and require to know the deployment information of wireless nodes. Fingerprinting-based methods estimate target locations by matching runtime measurements with fingerprints in a radio map. They can achieve enhanced accuracy, but a site survey is needed for radio map building and retraining. RTI-based methods do not require offline training effort. They treat target locations as the attenuation images of distorted links and can achieve an improved performance. However, they need a dense deployment of wireless links. It may lead to high hardware cost and great energy consumption. CS-based methods formulate the DFL problem as a sparse signal reconstruction problem. Compared with the aforementioned DFL methods, the CS-based DFL methods can ensure a high accuracy with fewer measurements. LCS [22] is a CS-based DFL method, in which the model-based dictionary satisfies the restricted isometry property (RIP) with high probability. E-HIPA [23] is a representative CS-based DFL method. It adopts an adaptive orthogonal matching pursuit (OMP) algorithm to estimate the target number and location vector. DR-DFL [14] introduced a real-time dictionary refinement algorithm for CS-based DFL. It can mitigate the influence of environmental dynamics in changing environments. The dictionary refinement is realized by optimizing the environment-related dictionary parameters. However, explicitly estimating the dictionary parameters may lead to a high computational complexity.

To leverage CSI for DFL, Pilot [24] regards the correlations of CSI as fingerprints and uses the maximum a posteriori probability (MAP) estimator to estimate target location. MonoPHY [25] is a fingerprinting-based DFL method that leverages CSI to locate a target with only one stream. Gao et al. [26] transformed CSI measurements into a radio image and adopted the machine learning method to estimate the position of a person. The aforementioned DFL method leverage CSI by fingerprinting technique. They seek appropriate location-dependent CSI features for DFL and trying to build a robust and precise relationship between CSI measurements and target position. However, to make full use of the channel diversity of CSI measurements is still a challenging issue in CS-based DFL.

3. Preliminaries and Problem Formulation

3.1. Overview of Multi-Target Device-Free Localization

The proposed ComDec method is a CS-based DFL method that aims to estimate the number and locations of multiple transceiver-free objects in changing environments. Figure 1 shows an example of the CS-based DFL. As can be seen, multiple transceiver-free objects are randomly distributed in an

l \times l

two-dimensional (2D) monitoring area

A

. Owing to the inherent spatial sparsity of the targets in

A

, the position information of multiple targets can be considered as a sparse signal. To achieve this, we discretize

A

into N equal-sized grids. When a target is located in grid n, we regard the grid center as its position. In this case, the target number and position information can be encoded in a sparse location vector, i.e.,

θ = {[θ_{1}, θ_{2}, . . ., θ_{n}, . . ., θ_{N}]}^{T},

(1)

where

θ \in R^{N \times 1}

is a K-sparse vector.

θ_{n} \in {0, 1}

is its n-th component, which indicates whether a target is located in grid n. If there is a target,

θ_{n} = 1

; otherwise,

θ_{n} = 0

. Thus, the location information of K targets can be denoted as

L = \{(x_{n}, y_{n}) |θ_{n} = 1, n \in \{1, . . ., N\}\}

, where

(x_{n}, y_{n})

is the coordinate of grid n. Moreover, the target number can be denoted as

K = {∥ θ ∥}_{0}

. When the target number K is far less than the grid number N (

K ≪ N

), the sparsity of

θ

can be ensured. As a matter of fact, the grid number N is usually a big number, while there is only a small number of targets in the monitoring area.

As seen from Figure 1,

A

is surrounded by several wireless nodes. Between these nodes, M bidirectional wireless links are established to cover

A

. They travel through

A

and sense the target-induced shadowing effects in the electromagnetic field. When there exists a target, some of the wireless links will be shadowed, and the signal power or other features of radio signals will be affected due to the scattering, reflection, refraction, and absorption of the signals. Fortunately, the changes of the signal features are closely related to target positions. Therefore, by analyzing the target-induced shadowing effects, the location information of multiple targets can be inferred.

3.2. CSI Collection and Feature Extraction

RSS is a coarse-grained measurement information from MAC layer. It represents the overall signal power across all subchannels. Due to the multipath fading, RSS is unreliable and varies with time even in a static environment. This may result in limited localization accuracy, especially in changing environments. To enhance the localization performance of CS-based DFL, ComDec adopts CSI to characterize target-induced perturbations. CSI is a fine-grained measurement information from PHY layer. It can provide the phase and strength information of the signals on different subchannels. For link m, the CSI on each subchannel is a complex value which is defined as

H_{m}^{f} = | H_{m}^{f} | e^{j \sin \{∠ H_{m}\}}, f = 1, 2, . . ., F,

(2)

where

H_{m}^{f}

denotes the CSI measurement corresponding to the f-th channel of link m.

| H_{m}^{f} |

and

∠ H_{m}

denote the amplitude and phase response, respectively. F is the number of channels. The amplitude response

| H_{m}^{f} |

is the change of amplitude of link m on channel f. By converting it from linear space to logarithmic space, the corresponding power fading can be written as

{\tilde{H}}_{m}^{f} = 20 log (| H_{m}^{f} | / 10^{3})

(dBm) [27]. We collect the CSI measurements from F channels, and a set of power fading information with channel diversity can be obtained. They can provide redundancy information to alleviate the location ambiguity that is incurred by environmental dynamics.

3.3. Problem Formulation

To characterize shadowing effects, we measure the the change of power fading on each link. For each channel, the change of power fading on link m is

\begin{matrix} Δ {\tilde{H}}_{m}^{f} = {\tilde{H}}_{m}^{f} - {\tilde{H}}_{m}^{f} (0) \approx S_{m}^{f} + ϵ_{m}^{f}, \end{matrix}

(3)

where

{\tilde{H}}_{m}^{f}

denotes the current power fading,

{\tilde{H}}_{m}^{f} (0)

is the reference power fading recorded when

A

is vacant,

S_{m}^{f}

represents the target-induced shadowing loss, and

ϵ_{m}^{f}

denotes the measurement noise.

The changes of power fading are measured on M links at runtime. For each link, we can obtain F measurements from F channels. Thus, the measurement information can be represented as

Y = [y^{1}, y^{2}, . . ., y^{f}, . . ., y^{F}] = {[\begin{matrix} y_{1}^{1} & y_{1}^{2} & \dots & y_{1}^{F} \\ y_{2}^{1} & y_{2}^{2} & \dots & y_{2}^{F} \\ ⋮ & ⋮ & ⋱ & ⋮ \\ y_{M}^{1} & y_{M}^{2} & \dots & y_{M}^{F} \end{matrix}]}_{M \times F,}

(4)

where

Y \in R^{M \times F}

is the measurement matrix.

y^{f} \in R^{M \times 1}

is the f-th column of

Y

. Its m-th component

y_{m}^{f} = Δ {\tilde{H}}_{m}^{f}

If any two targets are located sparsely [23] in the monitoring area,

S_{m}^{f}

can be represented as the summation of attenuations that occur in each cell [28]. Therefore,

Δ {\tilde{H}}_{m}^{f}

can be expressed as

Δ {\tilde{H}}_{m}^{f} = \sum_{n = 1}^{N} θ_{n} \cdot {\tilde{h}}_{m, n}^{f} + ϵ_{m}^{f},

(5)

where

{\tilde{h}}_{m, n}^{f}

is the shadowing loss that caused by a target located in grid n. For M links, the measurement vector

y^{f}

can be written as

y^{f} = Φ^{f} θ + ϵ^{f},

(6)

where

Φ^{f} \in R^{M \times N}

is a dictionary. Its (

m, n

)-th element is

{\tilde{h}}_{m, n}^{f}

ϵ^{f} \in R^{M \times 1}

represents the noise vector of channel f, and its m-th component is

ϵ_{m}^{f}

. In ComDec, the DFL problem can be viewed as a problem of reconstructing the location vector

θ

from the measurements

Y

. Theoretically, this problem can be solved by existing joint sparse recovery methods. Nonetheless, as a common shortcoming, they require the true dictionaries

{Φ^{f}}_{f = 1}^{F}

to be known in advance. However, in changing environments, it is impossible for us to accurately estimate these dictionaries.

Due to the difference in channel characteristics, the dictionaries with respect to different channels are different. Furthermore, in changing environments, the dictionary for an individual channel may differ when observed at different times. In this case, the fingerprinting method will be infeasible because it may cumulate the effects of dictionary mismatches on multiple channels and deteriorate the localization performance significantly. As an alternative, we adopt the saddle surface model to characterize the shadowing effect and establish multiple model-based dictionaries for multiple channels. To adapt to the changes of environment and different channel characteristics, the environmental parameters in these dictionaries are considered as adjustable. For simplicity, we denote

ϕ_{m, n}^{f} = {\tilde{h}}_{m, n}^{f}

, and thus

Φ^{f}

can be expressed as

Φ^{f} = {[ϕ_{1}^{f}, ϕ_{2}^{f}, . . ., ϕ_{m}^{f}, . . ., ϕ_{M}^{f}]}^{T} = {[\begin{matrix} ϕ_{1, 1}^{f} & ϕ_{1, 2}^{f} & \dots & ϕ_{1, N}^{f} \\ ϕ_{2, 1}^{f} & ϕ_{2, 2}^{f} & \dots & ϕ_{2, N}^{f} \\ ⋮ & ⋮ & ⋱ & ⋮ \\ ϕ_{M, 1}^{f} & ϕ_{M, 2}^{f} & \dots & ϕ_{M, N}^{f} \end{matrix}]}_{M \times N .}

(7)

Figure 2 depicts the spatial impact area of a wireless link. As can be seen, the spatial impact area is an ellipse area, and we set up an U-V coordinate system in the area. Based on the saddle surface model,

ϕ_{m, n}^{f}

can be expressed as

\begin{matrix} ϕ_{m, n}^{f} = γ_{m}^{f} \cdot (\frac{1 - ρ_{m}^{f}}{λ_{1}^{2}} U_{m, n}^{2} + ρ_{m}^{f} \cdot (1 - \frac{V_{m, n}^{2}}{λ_{2}^{2}})), \\ s . t . \frac{U_{m, n}^{2}}{λ_{1}^{2}} + \frac{V_{m, n}^{2}}{λ_{2}^{2}} ⩽ 1, \end{matrix}

(8)

where

(U_{m, n}, V_{m, n})

denotes the coordinate of grid n.

λ_{1}

and

λ_{2}

denote the semi-major and semi-minor axes of the spatial impact area, respectively. In the saddle surface model,

γ_{m}^{f}

represents the maximum shadowing loss, and

ρ_{m}^{f} \in (0, 1]

denotes the shadowing ratio, which represents the normalized shadowing loss at the midpoint of LOS path.

Due to

A

being an isotropic free space and all links sharing a common link length, we have

γ_{m}^{f} = γ^{f}

and

ρ_{m}^{f} = ρ^{f}

. This means that the values of

γ^{f}

and

ρ^{f}

are only determined by the environmental characteristics. In this case, the shadowing loss

ϕ_{m, n}^{f}

can be rewritten as

ϕ_{m, n}^{f} = \frac{U_{m, n}^{2}}{λ_{1}^{2}} γ^{f} + γ^{f} ρ^{f} (1 - \frac{U_{m, n}^{2}}{λ_{1}^{2}} - \frac{V_{m, n}^{2}}{λ_{2}^{2}}) .

(9)

Let

ω^{f} = γ^{f} \cdot θ

and

υ^{f} = γ^{f} ρ^{f} \cdot θ

be the unknown sparse vectors. Based on (6) and (9),

y^{f}

can be rewritten as

y^{f} = Ψ ω^{f} + Ψ^{'} υ^{f} + ϵ^{f},

(10)

where

Ψ \in R^{M \times N}

and

Ψ^{'} \in R^{M \times N}

are known dictionaries. Their (

m, n

)-th elements are computed as

ψ_{m, n} = U_{m, n}^{2} / λ_{1}^{2}

and

ψ_{m, n}^{'} = (1 - U_{m, n}^{2} / λ_{1}^{2} - V_{m, n}^{2} / λ_{2}^{2})

, respectively. We can determine their values before the localization stage according to the deployment of wireless links. It is noteworthy that the location vector

θ

and sparse vectors

S = {ω^{f}, υ^{f}}_{f = 1}^{F}

share a common sparse support set

T \subseteq {1, 2, . . ., N}

. In other words, their nonzero entries are concentrated at some common locations.

As pointed out earlier, the values of

{γ^{f}, ρ^{f}}_{f = 1}^{F}

are closely related to the environmental characteristics, which are different for different channels and times. Let

γ = [γ^{1}, γ^{2}, . . ., γ^{F}]

and

ρ = [ρ^{1}, ρ^{2}, . . ., ρ^{F}]

. We define two matrices

Θ, Θ^{'} \in R^{N \times F}

, which are constructed as follows:

Θ = θ γ = [ω^{1}, ω^{2}, . . ., ω^{F}],

(11)

Θ^{'} = θ (γ ⊙ ρ) = [υ^{1}, υ^{2}, . . ., υ^{F}] .

(12)

Both

Θ

and

Θ^{'}

are K jointly sparse matrices. This means there are at most K rows in them that have nonzero elements. The column vectors of

Θ

and

Θ^{'}

share the common support

T

, which is the index set of the grids where a target exists. The cardinality of

T

K : = |supp Θ| = |supp Θ^{'}|

, which also denotes the target number. Based on (10), (11) and (12), the measurement matrix

Y

can be expressed as

Y = Ψ Θ + Ψ^{'} Θ^{'} + E,

(13)

where

E \in R^{M \times F}

is the matrix of measurement noises, whose f-th column is

ϵ^{f}

Θ

and

Θ^{'}

are the unknown matrices. By (13), the CS-based DFL problem is recast as a problem that needs to recover

Θ

and

Θ^{'}

simultaneously over two known dictionaries. It should be pointed out that the problem is different from the conventional joint sparse recovery problem, which only needs to reconstruct a single sparse matrix over a single sparsifying dictionary. However, in this problem, there are two different sparse matrices to be estimated. Therefore, existing joint sparse recovery algorithms cannot be directly applied in ComDec. In this context, the key issue of ComDec is to design a joint sparse recovery algorithm to estimate

Θ

and

Θ^{'}

simultaneously. It should be noted that, the channel diversity of CSI measurements can be exploited to enhance the accuracy and robustness of the joint sparse recovery.

Figure 3 shows an architectural overview of ComDec. The proposed ComDec consists of four main modules: dictionary construction, measurement information collection, joint sparse recovery, and location estimation. The process of target counting and localization is illustrated as follows: First, in the dictionary construction module, we establish two constant dictionaries

Ψ

and

Ψ^{'}

according to the saddle surface model. Then, at runtime, CSI measurements are collected from M links. Each stream contains the CSI readings of F channels. Afterwards, in the joint sparse recovery module, the posteriors of all hidden variables are inferred by a joint sparse recovery algorithm. Finally, with the knowledge of the posteriors, in the location estimation module, the common support set

T

can be estimated, and thus the estimated Cartesian coordinates of multiple targets can be obtained.

4. Target Localization via Variational Bayesian Inference

In this section, we develop a novel joint sparse recovery algorithm. First, to induce the joint sparsity of the jointly sparse vectors

S

, a two-layer hierarchical prior model is introduced. Then, by using the variational Bayesian inference technique, we infer the posteriors of the hidden variables that defined in the proposed hierarchical prior model. Finally, based on the posteriors of

S

, target counting and localization are implemented.

4.1. Hierarchical Prior Model

The joint sparsity of

S

is induced by a non-separable sparsity inducing prior model [29]. The graphical model for the joint sparse recovery is shown in Figure 4, which describes the dependencies between random variables. In the first layer, we regard

{ω^{f}, υ^{f}}_{f = 1}^{F}

as stochastic variables and define

ω = {ω^{f}}_{f = 1}^{F}

and

υ = {υ^{f}}_{f = 1}^{F}

. Moreover, the Gaussian-inverse-Gamma prior is imposed on each sparse vector to encourage its sparsity.

ω^{f}

is treated as a Gaussian random variable, whose prior distribution is can be given as

p (ω^{f} |α) = \prod_{n = 1}^{N} N (ω_{n}^{f} |0, α_{n}^{- 1}) = {(2 π)}^{- N / 2} {|Λ|}^{1 / 2} exp (\frac{1}{2} {(ω^{f})}^{T} Λ ω^{f}),

(14)

where

N (\cdot |0, α_{n}^{- 1})

denotes the Gaussian distribution with zero mean and variance of

α_{n}^{- 1}

α_{n}

is the inverse variance (precision) of

{ω_{n}^{f}}_{f = 1}^{F}

. We define

α = {[α_{1}, α_{2}, . . ., α_{N}]}^{T}

and

Λ = diag (α)

. As

υ^{f} = ρ^{f} ω^{f}

, where

ρ^{f}

is regarded as an unknown deterministic parameter, a multivariate Gaussian prior is also imposed on

υ^{f}

. The variance of

υ_{n}^{f}

can be given as

{(ρ^{f})}^{2} α_{n}^{- 1}

. We set

ρ^{f} = 1

to accommodate the worst case of the variance. In this situation, the prior of

υ^{f}

can be expressed as

p (υ^{f} |α) = \prod_{n = 1}^{N} N (υ_{n}^{f} |0, α_{n}^{- 1}) .

(15)

The variances of

υ_{n}^{f}

and

ω_{n}^{f}

are identical. Obviously, a large value of

α_{n}

will simultaneously drive

υ_{n}^{f}

and

ω_{n}^{f}

to zero, which is in correspondence with the joint sparsity between

υ^{f}

and

ω^{f}

. Intuitively, if most of the components of

α

have large values, the jointly sparse solutions will be obtained. Hence, we model the hyperparameters

{α_{n}}_{n = 1}^{N}

as hidden variables to allow the flexibility to learn and adapt to the true situation.

To allow conjugate-exponential analysis [20], an independent identically distributed (i.i.d.) Gamma distribution is imposed on

α

, which can be expressed as

p (α; c, d) = \prod_{n = 1}^{N} Gamma (α_{n} |c, d_{n}) = \prod_{n = 1}^{N} \frac{1}{Γ (c)} d_{n}^{c} α_{n}^{c - 1} exp (- d_{n} α_{n}),

(16)

where

Gamma (\cdot |c, d_{n})

denotes the Gamma distribution with parameters c and

d_{n}

Γ (c) = \int_{0}^{\infty} x^{c - 1} e^{- x} d x

is the Gamma function,

d = {[d_{1}, d_{2}, . . ., d_{N}]}^{T}

. We set the hyperparameters c and

d_{n}

to very small values (e.g.,

10^{- 6}

) to provide a non-informative hyperprior over

α_{n}

. The Gamma distribution is generally chosen as the prior for the inverse variance of a Gaussian distribution, because it is the conjugate prior of the Gaussian distribution. In this case, the associated Bayesian inference can be performed in closed form [4].

When a hierarchical Gaussian prior model imposed on

ω^{f}

, the true prior distribution of

ω^{f}

can be computed by marginalizing the parameter

α

, i.e.,

\begin{matrix} p (ω^{f}; c, d) & = \int p (ω^{f} | α) p (α; c, d) d α \\ = \int \prod_{n = 1}^{N} N (ω_{n}^{f} | 0, α_{n}^{- 1}) Gamma (α_{n} | c, d) d α_{n} \\ = \prod_{n = 1}^{N} St (ω_{n}^{f} | λ, v) . \end{matrix}

(17)

In this case, the true distribution of

ω^{f}

is a Student-t pdf,

St (x | κ, λ, v) = \frac{Γ ((v + 1) / 2)}{Γ (v / 2)} {(\frac{λ}{π v})}^{1 / 2} {(1 + \frac{λ {(x - κ)}^{2}}{v})}^{- (v + 1) / 2},

(18)

with mean

κ = 0

, parameter

λ = c / d

and degrees of freedom

v = 2 c

. According to the property of the Student-t distribution, when v is small, this distribution will exhibit very heavy tails. Thus, it favours sparse solutions, which include only few of the basis functions and prunes the remaining basis functions by setting the corresponding weights to very small values. In this case, a sparse vector

ω^{f}

can be induced with the hierarchical Gaussian prior model. As the case of

ω^{f}

, we can also induce the sparsity of

υ^{f}

by the proposed hierarchical Gaussian prior model.

As illustrated earlier,

ϵ^{f}

follows an i.i.d. Gaussian distribution. The prior distribution of

ϵ^{f}

is defined as

p (ϵ^{f} | β^{f}) = \prod_{n = 1}^{N} N (ϵ_{m}^{f} |0, {(β^{f})}^{- 1}) = {(\frac{β^{f}}{2 π})}^{\frac{M}{2}} exp (- \frac{β^{f}}{2} {∥ ϵ^{f} ∥}_{2}^{2}),

(19)

where

β^{f}

is the inverse variance of

ϵ^{f}

. We treat

β^{f}

as a hidden variable and define

β = {β^{f}}_{f = 1}^{F}

. As

ϵ^{f}

follows a Gaussian distribution, a Gamma prior is also imposed on each

β^{f}

, i.e.,

p (β^{f}; a^{f}, b^{f}) = Gamma (β^{f} | a^{f}, b^{f}) = \frac{{(b^{f})}^{a^{f}} {(β^{f})}^{a^{(f)} - 1}}{Γ (a^{f})} exp (- b^{f} β^{f}),

(20)

where

a^{f}

and

b^{f}

are deterministic parameters of the Gamma distribution. We denote

a = {a^{f}}_{f = 1}^{F}

and

b = {b^{f}}_{f = 1}^{F}

. To assume uninformative priors for

β

, the hyperparameters

a

and

b

are also set to very small values (e.g.,

10^{- 6}

In the proposed graphical model, the observed variables are

y = {y^{f}}_{f = 1}^{F}

, and the hidden variables are

z ≜ \{α, ω, υ, β\}

. To estimate the jointly sparse vectors, we need to infer the posterior distributions of

z

based on the predefined prior evidence and the measurement data. In addition, the deterministic parameters of the prior model are

Ω ≜ \{c, d, a, b\}

, which are fixed at small values to allow uninformative hyperpriors for

α

and

β

4.2. Variational Bayesian Inference

In the joint sparse recovery module, the key task is to infer the posteriors of

z

. Afterwards, based on these posteriors, the target number and locations can be estimated. For this objective, the variational Bayesian inference technique is adopted, which is applied due to it can deal with complicated Bayesian models [20]. Based on (10) and (19), the likelihood function of channel f can be written as

p (y^{f} | ω^{f}, υ^{f}, β^{f}) = {(\frac{β^{f}}{2 π})}^{\frac{M}{2}} exp (- \frac{β^{f}}{2} {∥ y^{f} - Ψ ω^{f} - Ψ^{'} υ^{f} ∥}_{2}^{2}) .

(21)

According to the chain rule of probability, the joint probability density function (PDF) of

y

and

z

can be written as

\begin{matrix} p (y, z) & = p (y |ω, υ, β) p (β) p (ω |α) p (υ |α) p (α) \\ = \prod_{f = 1}^{F} p (y^{f} | ω^{f}, υ^{f}, β^{f}) \prod_{f = 1}^{F} p (β^{f}) \prod_{f = 1}^{F} \prod_{n = 1}^{N} p (ω_{n}^{f} |α_{n}) \prod_{f = 1}^{F} \prod_{n = 1}^{N} p (υ_{n}^{f} |α_{n}) \prod_{n = 1}^{N} p (α_{n}) . \end{matrix}

(22)

Figure 5 illustrates the factor graph representation of the joint PDF. In this figure, the circle nodes represent the random variables that the complicated global function relied on, while the square nodes represent the local functions.

For an arbitrary density function

q (z)

, the evidence

p (y) = \int p (y, z) d z

can be decomposed as

ln p (y) = F (q; Ω) + K L (q ∥p),

(23)

where

F (q; Ω) = \int q (z) ln (\frac{p (y, z; Ω)}{q (z)}) d z,

(24)

K L (q ∥p) = - \int q (z) ln (\frac{p (z |y; Ω)}{q (z)}) d z .

(25)

F (q; Ω)

is a lower bound of

ln p (y)

, and

K L (q ∥p)

represents the Kullback-Leibler divergence (KLD) between the approximated posterior

q (z)

and the exact posterior

p (z |y; Ω)

. The proposed joint sparse recovery algorithm maximizes

ln p (y)

iteratively. At each iteration, we set

K L (q ∥p) = 0

to minimize the KLD and update

q (z)

accordingly. This will lead to the lower bound

F (q, Ω)

increasing to

ln p (y)

. Meanwhile, the updating of

q (z)

may enlarge

ln p (y)

and lead to a new non-negative KLD. We will minimize the new KLD and update

q (z)

in the next iteration. By doing so, the log-likelihood

ln p (y)

will be maximized, and the approximated posterior

q (z)

can be optimized iteratively.

However,

p (z |y; Ω)

cannot be computed analytically. Thereby, directly updating

q (z)

is intractable. To bypass this difficulty, we resort to the variational approximation method. It assumes that the posteriors of

z

are independent, i.e.,

q (z) = q (α) q (ω) q (υ) q (β) = q (α) \prod_{f = 1}^{F} q (ω^{f}) \prod_{f = 1}^{F} q (υ^{f}) \prod_{f = 1}^{F} q (β^{f}) .

(26)

By applying the above assumption, in each iteration, the log-posterior of

z_{i} \in z

can be approximated as the expectation of the joint PDF with respect to other hidden variables

(z - z_{i})

. More specifically, the log-posteriors of

z

are approximated as

\begin{matrix} ln q (ω^{f}) & = {〈ln p (y, z)〉}_{q (α) \prod_{i \neq f} q (ω^{i}) q (υ) q (β)} + ξ, \end{matrix}

(27)

\begin{matrix} ln q (υ^{f}) & = {〈ln p (y, z)〉}_{q (α) q (ω) \prod_{i \neq f} q (υ^{i}) q (β)} + ξ, \end{matrix}

(28)

\begin{matrix} ln q (β^{f}) & = {〈ln p (y, z)〉}_{q (α) q (ω) q (υ) \prod_{i \neq f} q (β^{i})} + ξ, \end{matrix}

(29)

\begin{matrix} ln q (α) & = {〈ln p (y, z)〉}_{q (ω) q (υ) q (β)} + ξ, \end{matrix}

(30)

where

ξ

denotes a normalizing constant. It makes the corresponding

q (\cdot)

a true PDF. The update rule for each hidden variable is derived below.

In (27), the terms independent of

ω^{f}

can be treated as a constant value. Keeping only the terms that are related to

ω^{f}

ln q (ω^{f})

can be given as

ln q (ω^{f}) \propto {〈ln p (y^{f} | ω^{f}, υ^{f}, β^{f}) + ln p (ω^{f} | α)〉}_{q (α) q (υ^{f}) q (β^{f}) .}

(31)

Note that

q (ω^{f})

and

q (υ^{f})

follow the multivariate Gaussian distribution. We assume they have the following forms after the updating

\begin{matrix} q (ω^{f}) & = N (ω^{f} | μ_{ω}^{f}, Σ_{ω}^{f}), \end{matrix}

(32)

\begin{matrix} q (υ^{f}) & = N (υ^{f} | μ_{υ}^{f}, Σ_{υ}^{f}), \end{matrix}

(33)

where

μ_{ω}^{f}

and

μ_{υ}^{f}

are the mean vectors,

Σ_{ω}^{f}

and

Σ_{υ}^{f}

are the covariance matrices. The update of the posterior distribution is equivalent to seeking appropriate values for the parameters in the approximated posterior distribution. Our goal is to learn the values of the mean vectors and covariance matrices based on the prior distributions and likelihood function. Substituting (14) and (21) into (31), after some rearrangement, yields

\begin{matrix} ln q (ω^{f}) \propto - \frac{1}{2} {(ω^{f})}^{T} {(Σ_{ω}^{f})}^{- 1} ω^{f} + {(ω^{f})}^{T} {(Σ_{ω}^{f})}^{- 1} μ_{ω}^{f}, \end{matrix}

(34)

where

\begin{matrix} Σ_{ω}^{f} & = {(〈 β^{f} 〉 Ψ^{T} Ψ + 〈Λ〉)}^{- 1}, \end{matrix}

(35)

\begin{matrix} μ_{ω}^{f} & = 〈 β^{f} 〉 Σ_{ω}^{f} Ψ^{T} (y^{f} - Ψ^{'} μ_{υ}^{f}) . \end{matrix}

(36)

In the same manner,

ln q (υ^{f})

can be given as

ln q (υ^{f}) \propto {〈ln p (y^{f} | ω^{f}, υ^{f}, β^{f}) + ln p (υ^{f} |α)〉}_{q (α) q (ω^{f}) q (β^{f}) .}

(37)

Substituting (15) and (21) into (37), after some rearrangement, yields

\begin{matrix} ln q (υ^{f}) \propto & - \frac{1}{2} {(υ^{f})}^{T} {(Σ_{υ}^{f})}^{- 1} υ^{f} + {(υ^{f})}^{T} {(Σ_{υ}^{f})}^{- 1} μ_{υ}^{f}, \end{matrix}

(38)

where

\begin{matrix} Σ_{υ}^{f} & = {(〈 β^{f} 〉 {Ψ^{'}}^{T} Ψ^{'} + 〈Λ〉)}^{- 1}, \end{matrix}

(39)

\begin{matrix} μ_{υ}^{f} & = 〈 β^{f} 〉 Σ_{υ}^{f} {Ψ^{'}}^{T} (y^{f} - Ψ μ_{ω}^{f}) . \end{matrix}

(40)

Keeping only the terms of (29) that are related to

β^{f}

ln q (β^{f})

can be given as

ln q (β^{f}) \propto {〈ln p (y^{f} | ω^{f}, υ^{f}, β^{f}) + ln p (β^{f})〉}_{q (ω^{f}) q (υ^{f}) .}

(41)

We assume that the posterior of

β^{f}

follows a Gamma distribution, i.e.,

q (β^{f}) = Gamma (β^{f} | {\tilde{a}}^{f}, {\tilde{b}}^{f}),

(42)

where

{\tilde{a}}^{f}

and

{\tilde{b}}^{f}

denote the deterministic parameters of the updated posterior distribution. To infer them, we substitute (20) and (21) into (41). After some rearrangement, the posterior can be given as

ln q (β^{f}) \propto ({\tilde{a}}^{f} - 1) ln β^{f} - {\tilde{b}}^{f} β^{f},

(43)

where

\begin{matrix} {\tilde{a}}^{f} = & a^{f} + M / 2, \end{matrix}

(44)

\begin{matrix} {\tilde{b}}^{f} = & b^{f} + \frac{1}{2} {∥y^{f} - Ψ μ_{ω}^{f} - Ψ^{'} μ_{υ}^{f}∥}_{2}^{2} + \frac{1}{2} (tr ({Ψ Σ}_{ω}^{f} Ψ^{T}) + tr ({Ψ^{'} Σ}_{υ}^{f} {Ψ^{'}}^{T})) . \end{matrix}

(45)

Similarly, we only keep the terms that are related to

α

in (30), and thereby

ln q (α)

can be given as

ln q (α) \propto {〈\sum_{f = 1}^{F} ln p (ω^{f} | α) + \sum_{f = 1}^{F} ln p (υ^{f} | α) + ln p (α)〉}_{\prod_{f = 1}^{F} q (ω^{f}) \prod_{f = 1}^{F} q (υ^{f}) .}

(46)

The approximated posterior of

α

is assumed to be a multivariate Gamma distribution, i.e.,

q (α) = \prod_{n = 1}^{N} Gamma (α_{n} | \tilde{c}, {\tilde{d}}_{n}),

(47)

where

\tilde{c}

and

{{\tilde{d}}_{n}}_{n = 1}^{N}

are deterministic parameters. Substituting (14), (15) and (16) into (46), after some rearrangement, yields

ln q (α) \propto (\tilde{c} - 1) \sum_{n = 1}^{N} ln α_{n} - \sum_{n = 1}^{N} {\tilde{d}}_{n} α_{n},

(48)

where

\begin{matrix} \tilde{c} & = c + F, \end{matrix}

(49)

\begin{matrix} {\tilde{d}}_{n} & = d_{n} + \frac{1}{2} \sum_{f = 1}^{F} ({[μ_{ω}^{f}]}_{n}^{2} + {[Σ_{ω}^{f}]}_{n, n} + {[μ_{υ}^{f}]}_{n}^{2} + {[Σ_{υ}^{f}]}_{n, n}) . \end{matrix}

(50)

The notation

{[\cdot]}_{n}

denotes the n-th component of the input vector, and

{[\cdot]}_{n, n}

denotes the

(n, n)

-th entry of the input matrix. Based on the results of posterior inference, the required expectations can be calculated as

\begin{matrix} 〈 β^{f} 〉 = {\tilde{a}}^{f} / {\tilde{b}}^{f}, f = 1, 2, . . ., F . \end{matrix}

(51)

\begin{matrix} 〈 Λ 〉 = diag (〈α_{1}〉, 〈α_{2}〉, . . ., 〈α_{N}〉), \end{matrix}

(52)

where

〈α_{n}〉 = \tilde{c} / {\tilde{d}}_{n}, n = 1, 2, . . ., N .

(53)

4.3. Joint Sparse Reconstruction

According to the above update rules, we can successively update the posteriors of hidden variables

z

. The jointly sparse vectors

S

can be reconstructed according to these posteriors. The algorithm of reconstructing

S

is summarized as follows:

For $f \in \{1, 2, . . ., F\}$ , update $q (ω^{f})$ by using (35) and (36); update $q (υ^{f})$ by using (39) and (40). $〈 β^{f} 〉$ and $〈Λ〉$ are obtained based on the current posteriors of $β^{f}$ and $α$ .
For $f \in \{1, 2, . . ., F\}$ , update $q (β^{f})$ according to (44), (45) and the current posteriors of $ω^{f}$ and $υ^{f}$ .
Update $q (α)$ according to (49), (50) and the current posteriors of $ω$ and $υ$ .
If a convergence criterion has been met, terminate and choose the posterior means of $ω$ and $υ$ as the estimation of $S$ . Otherwise, go to step $1)$ .

Based on the above joint sparse recovery algorithm, we can obtain a reconstructed sparse vector set

\hat{S} = {μ_{ω}^{f}, μ_{υ}^{f}}_{f = 1}^{F}

. The computational cost of the algorithm is dominated by the matrix inversion operations in (35) and (39), whose computational complexities are

O (N^{3} F)

. Moreover, the matrix-vector multiplications in (36) and (40), as well as the matrix multiplication operations in (45) can also incur heavy computational burden. Their computational complexities are

O (N^{2} M F)

. Thus, when applying ComDec in large-scale areas (where N is large), the joint sparse recovery algorithm will be computationally expensive. In this paper, we adopt the following four methods to alleviate the computational burden of the proposed algorithm.

In step 1, the covariance matrices

Σ_{ω}^{f}

and

Σ_{υ}^{f}

are computed by (35) and (39), which contain matrix inversion operations. Using the matrix inversion lemma [30], the covariance matrices can be evaluated as

\begin{matrix} Σ_{ω}^{f} & = {〈Λ〉}^{- 1} - {〈Λ〉}^{- 1} Ψ^{T} {(Ξ_{ω}^{f})}^{- 1} Ψ {〈Λ〉}^{- 1}, \end{matrix}

(54)

\begin{matrix} Σ_{υ}^{f} & = {〈Λ〉}^{- 1} - {〈Λ〉}^{- 1} {Ψ^{'}}^{T} {(Ξ_{υ}^{f})}^{- 1} Ψ^{'} {〈Λ〉}^{- 1}, \end{matrix}

(55)

where

\begin{matrix} Ξ_{ω}^{f} & = Ψ {〈Λ〉}^{- 1} Ψ^{T} + {〈 β^{f} 〉}^{- 1} I_{M}, \end{matrix}

(56)

\begin{matrix} Ξ_{υ}^{f} & = Ψ^{'} {〈Λ〉}^{- 1} {Ψ^{'}}^{T} + {〈 β^{f} 〉}^{- 1} I_{M} . \end{matrix}

(57)

With these matrix transformations, we only need to compute the inversions of

Ξ_{ω}^{f}

and

Ξ_{υ}^{f}

, whose computational complexities are

O (M^{3})

(M ≪ N)

. Note that,

〈Λ〉

is an

N \times N

diagonal matrix. We can easily obtain its inversion.

As mentioned before, the expectation evaluations of

μ_{ω}^{f}

and

μ_{υ}^{f}

can also lead to a high computational cost. To reduce their computational complexities, we reformulate (36) and (40) to cast them as a problem of solving the following linear systems of equations:

\begin{matrix} (Ψ^{T} Ψ + \frac{〈Λ〉}{〈 β^{f} 〉}) μ_{ω}^{f} & = Ψ^{T} (y^{f} - Ψ^{'} μ_{υ}^{f}), \end{matrix}

(58)

\begin{matrix} ({Ψ^{'}}^{T} Ψ^{'} + \frac{〈Λ〉}{〈 β^{f} 〉}) μ_{υ}^{f} & = {Ψ^{'}}^{T} (y^{f} - Ψ μ_{ω}^{f}) . \end{matrix}

(59)

The equations can be solved by the conjugate gradient (CG) algorithm. Theoretically, we can reach the exact solution after N iterations. However, this may incur a considerable computational burden. Fortunately, in practice, a few iterations is sufficient to obtain a good solution. If W

(W ≪ N)

iterations are required, the computational complexity of the CG algorithm will be

O (W N log N)

In step 2, the computational cost is mainly attributed to the matrix multiplications in the last term of (45). To mitigate this, (45) can be rewritten as

{\tilde{b}}^{f} = b^{f} + \frac{1}{2} {∥y^{f} - Ψ μ_{ω}^{f} - Ψ^{'} μ_{υ}^{f}∥}_{2}^{2} + \frac{1}{〈 β^{f} 〉} \sum_{n = 1}^{N} (1 - \frac{〈 α_{n} 〉}{2} ({[Σ_{ω}^{f}]}_{n, n} + {[Σ_{υ}^{f}]}_{n, n})) .

(60)

Here, matrix multiplications are replaced by simple element-wise operations. As a consequence, the computational complexity of step 2 can be reduced to

O (N M F)

. By using the CG method and the above reformulations, the computational cost per iteration of the proposed joint sparse recovery algorithm can be reduced to

O (N M^{2} F)

To further reduce the computational load and speed up the convergence, we conduct real-time grid pruning according to the posterior of

α

. In each iteration, when

α_{n}

is sufficiently large,

{ω_{n}^{f}, υ_{n}^{f}}_{f = 1}^{F}

will be driven to zero. This imply that the contribution of grid n to the signal power fading is negligible. In this case, we can remove grid n from the grid set

Π

, which is defined as the set of grids where a target possibly exists. With the reduction of grids in

Π

, the computational load of the next iteration will decrease.

We denote the initial grid set by

Π^{(0)} ≜ \{1, 2, . . ., N\}

. Then, the grid set in each iteration can be updated as

Π^{(τ)} = Π^{(τ - 1)} - {n |〈α_{n}〉 > α_{th}},

(61)

where

τ

denotes the iteration number and

α_{th}

is the threshold of

α_{n}

. For the reconstruction of

S

, the selection of

α_{th}

provides a trade-off between the localization accuracy and the convergence speed. After grid pruning, the dimension of

θ

can be reduced to

N = | Π^{(τ)} |

. In addition,

α

S

Ψ

, and

Ψ^{'}

also should be pruned accordingly.

By virtue of the local convergence property of variational Bayesian inference [20], the proposed joint sparse recovery algorithm is guaranteed to be convergent. The stop criterion is set as the residual error

R e s

smaller than a pre-determined threshold

r_{th}

, where the

R e s

is defined as

R e s = \sum_{f = 1}^{F} {∥y^{f} - Ψ μ_{ω}^{f} - Ψ^{'} μ_{υ}^{f}∥}_{2}^{2} .

(62)

Furthermore, to prevent the computational load from being excessively high, we set a maximum iteration number

τ_{max}

. In this situation, the iterative algorithm will stop when

R e s

becomes smaller than

r_{th}

τ

reaches

τ_{max}

4.4. Target Counting and Localization

Based on

\hat{S}

, the target number and locations can be estimated. However, these vectors are not strictly sparse. They may contain many negligible but nonzero coefficients. Thus, we estimate their common support by the following formulae:

\hat{T} = \{n |20 lg (\frac{{[μ_{ω}^{\hat{f}}]}_{n}}{max_{i \in Π} {[μ_{ω}^{\hat{f}}]}_{i}}) ⩾ μ_{th}, n \in Π\},

(63)

where

\hat{f} = arg min_{f} {∥y^{f} - Ψ μ_{ω}^{f} - Ψ^{'} μ_{υ}^{f}∥}_{2}^{2} .

(64)

The common support of the jointly sparse vectors is estimated based on

μ_{ω}^{\hat{f}}

, in which

\hat{f}

is the channel that has the minimum residual error. In (63),

μ_{th}

is the sparsity threshold. We use it to filter out the small coefficients in

μ_{ω}^{\hat{f}}

. After that, we can estimate the target number and locations as

\hat{K} = | \hat{T} |

and

\hat{L} = {(x_{n}, y_{n}) | n \in \hat{T}}

, respectively.

5. Numerical Evaluation

5.1. Simulation Setup

The numerical simulations are carried out in MATLAB R2015b 64 bit version running on a PC with i7-8550U and 16 GB memory. In our simulations, the channel state information is assumed to be collected from the wireless devices that implement an orthogonal frequency division multiplexing (OFDM) system. The wireless devices are assumed to mode at 2.4 GHz with 20 MHz bandwidth. Table 1 summarizes the default values of simulation parameters. The numerical simulations are conducted in a complex radio transmission environment. In this transmission environment, as the impact of environmental dynamics, the environmental parameters are changing with time. Moreover, the measurement of each channel is corrupted by addictive white Gaussian noise. The signal-to-noise ratio (SNR) of each channel is defined as

SNR = 10 lg \frac{∥ Φ^{f} θ ∥_{2}^{2}}{M {(σ^{f})}^{2}}, f = 1, 2, . . ., F,

(65)

where

σ^{f}

is the standard deviation of the measurement noise vector.

In our simulations, all results are averaged over

T = 10^{3}

Monte Carlo trials. For each trial, the localization error

E_{t}

is computed as

E_{t} = \frac{\sum_{k = 1}^{min {K, \hat{K}}} {∥ L_{k} - {\hat{L}}_{k} ∥}_{2}}{min {K, \hat{K}}},

(66)

where

L_{k}

and

{\hat{L}}_{k}

represent the true and estimated Cartesian coordinates of target k. We use the “average localization error” (

A v g . E r r o r

) and “root-mean-square error” (

R m s . E r r o r

) as the metrics to measure the localization accuracy. The

A v g . E r r o r

is defined as

A v g . E r r o r = \frac{\sum_{t = 1}^{T} E_{t}}{T} .

(67)

The

R m s . E r r o r

is defined as

R m s . E r r o r = \sqrt{\frac{\sum_{t = 1}^{T} E_{t}^{2}}{T}} .

(68)

To evaluate the target counting performance, we introduce another metric

P r o b . C o C

. It is the probability of correct counting (i.e.,

\hat{K} = K

). We compare the performances of ComDec with other three CS-based DFL approaches: (i) LCS with the GMP algorithm (LCS-GMP) as the CS recovery algorithm [22], (ii) E-HIPA with the adaptive OMP algorithm (E-HIPA-OMP) as the CS recovery algorithm [23], and (iii) DR-DFL with the VEM algorithm for dictionary refinement and sparse recovery (DR-DFL-VEM) [14]. Table 2 reports the computational complexities and localization accuracies of multiple DFL methods. As can be seen from this table, the ComDec method can achieve the lowest average localization error (

A v g . E r r o r

). It can enhance the accuracy and robustness of CS-based DFL by exploiting channel diversity. The computational complexity of the ComDec method is higher than the LCS and E-HIPA methods. The inferiority is mainly attributed to the estimation of the posterior distributions of sparse vectors. It noteworthy that the E-HIPA and LCS methods reconstruct the unknown sparse vector by using the OMP and GMP algorithms, respectively. They are greedy algorithms for sparse recovery, and can only provide a point estimation of the unknown sparse signal. In contrast, the proposed sparse recovery algorithm and VEM algorithm reconstruct sparse signal from a Bayesian perspective. They can provide a posterior belief (distribution function) for the values of the sparse signal, and therefore can achieve an improved accuracy. In real deployment, the sparse recovery algorithm is usually conducted on a server, where the computational resources are abundant. The server collects CSI measurements from all wireless nodes and estimates the target location accordingly. The simulation flowchart is shown in Figure 6. In the upper part of the flowchart, the measurement data is generated. In the lower part of the flowchart, the target locations can be estimated.

5.2. Impact of the Number of Channels

In this subsection, the impact of channel number on the performance of ComDec is evaluated. Theoretically, the accuracy of joint sparse recovery is closely related to the channel number. If we increase F, more useful information will be provided. Consequently, the target counting and localization performance will be improved.

In Figure 7, F varies from 1 to 25, and the

A v g . E r r o r

of ComDec decreases as F increases. It is noteworthy that, the

A v g . E r r o r

decreases dramatically when

F < 10

, while decreases slowly when

F > 10

. As in the case of

A v g . E r r o r

, the

P r o b . C o C

also increases as F increases, and its growth rate drops gradually. This is because the increase of F leads to a more complex hierarchical prior model, which contains (

3 F + 1

) hidden variables. When reconstructing

S

, the posterior distributions of these hidden variables should be inferred. Thus, a large F may deteriorate the accuracy of the posterior inference. From this perspective, increasing the channel number has a negative effect on sparse recovery. In fact, when

F > 20

, the negative effect of increasing F will offset the positive effect contributed by channel diversity. Hence, we set

F = 15

in our simulation to achieve a trade-off between model simplicity and the performance of target counting and localization.

5.3. Effectiveness of ComDec in Changing Environments

To mitigate the dictionary mismatches introduced by environmental dynamics, the proposed ComDec method combines the environmental parameters and the location vector to form a set of jointly sparse vectors. We treat these vectors as random variables and learn their values by variational Bayesian inference. This makes the dictionary parameters can adapt to the environmental changes and different channel characteristics. In this subsection, we test the effectiveness and robustness of ComDec in changing environments.

Firstly, we set up a simulated changing environment. To simulate the environmental dynamics in changing environments, we add Gaussian noises to the environmental parameters (

γ^{f}

and

ρ^{f}

). Based on these noisy environmental parameters and the saddle surface model, the dictionaries

{Φ^{f}}_{f = 1}^{F}

can be built. After that, we can obtain the measurement vectors

{y^{f}}_{f = 1}^{F}

according to (6). In our simulation, the values of

γ^{f}

and

ρ^{f}

are calculated by

\begin{matrix} γ_{t}^{f} & = γ_{0}^{f} + \sum_{i = 1}^{t} ϖ_{γ}^{i}, f = 1, 2, . . ., F, \end{matrix}

(69)

\begin{matrix} ρ_{t}^{f} & = ρ_{0}^{f} + \sum_{i = 1}^{t} ϖ_{ρ}^{i}, f = 1, 2, . . ., F, \end{matrix}

(70)

where

ϖ_{γ}^{i}

and

ϖ_{ρ}^{i}

are additive white Gaussian noises whose variances are

0.6

and

0.06

, respectively. The setting of the noise variances is according to the real data in changing environments. t represents the times of environmental changes and also denotes the number of additive noises added to the original environmental parameters (

γ_{0}^{f}

and

ρ_{0}^{f}

). In our simulations, we set

γ_{0}^{f} = 17.0831

and

ρ_{0}^{f} = 0.35122

. In addition, uncorrelated Gaussian noises are added to the dictionaries

{Φ^{f}}_{f = 1}^{F}

to produce an SNR of 40 dB. Figure 8 depicts the variations of two dictionary atoms in the simulated changing environment. As expected, the accumulation of additive noises can lead to a deviation of dictionary atoms, which is a simulation of the effect of environmental dynamics in real-world environments.

Secondly, in the simulated changing environment, the localization performances of multiple CS-based multi-target DFL methods are evaluated. Figure 9 compares the localization performances when t increases from 1 to 10. The proposed ComDec method that using the variational Bayesian inference technique for joint sparse recovery (ComDec-VBI) can achieve the minimum

A v g . E r r o r

. It should be noted that, its localization performance is robust to the environmental changes. In contrast, with the increase of t, the localization accuracies of other three DFL methods deteriorate gradually. This reveals that these methods are sensitive to environmental dynamics. The simulation results demonstrate the impact of environmental dynamics on

A v g . E r r o r

and verify the accuracy and robustness of the ComDec method.

Finally, we investigate the sparse reconstruction performances of the CS recovery algorithms that adopted in the CS-based DFL methods. We run E-HIPA, LCS, DR-DFL, and the ComDec method with

F = 5

and

F = 15

to reconstruct the target location vector. Figure 10 shows an example of the recovered sparse vectors of different DFL approaches. As can be seen, the recovered sparse vectors corresponding to E-HIPA and DR-DFL have many negligible but nonzero coefficients. Although the reconstructed sparse vector in Figure 10c does not have small coefficients, the indices of the nonzero components are different from those of the original location vector. Fortunately, the recovered sparse vector corresponding to our proposed ComDec method has a few small coefficients, and the indices of their significant coefficients are the same as those in the original location vector. Based on them, the target number and target locations can be estimated correctly.

Figure 11 shows the target positions estimated by multiple CS-based DFL methods. As can be seen, the ComDec method can accurately estimate target positions, whereas other DFL methods have 1 or 2 incorrect estimated positions. This is because the proposed ComDec method can eliminate the position ambiguities introduced by environmental dynamics by taking advantage of the channel diversity.

5.4. Localization Error vs. SNR

In this subsection, the localization performances of different DFL methods under different SNRs are evaluated. In cellular communication (e.g., LTE), the typical SNR value is below 0 dB. However, in short range communication protocols such as 802.11 b/g, the typical SNR value is from 10–15 dB (low signal) to 40 dB (excellent signal). In fact, the scenario that is chosen for our simulations is the short-range communication scenario. Thus, the SNR in this simulation is varied from 0 dB to 40 dB. Figure 12 shows the

A v g . E r r o r

of the ComDec method

(F = 10

, 15, and

20)

and other DFL methods. As can be seen, the

A v g . E r r o r

of all DFL methods decrease with the increase of SNR. Generally, the localization error is mainly caused by the measurement noise and the environmental dynamics in changing environments. As the ComDec and DR-DFL have countermeasures to mitigate the dictionary mismatches introduced by environmental dynamics, their

A v g . E r r o r

decrease dramatically when SNR increases from 0 dB to 25 dB.

Furthermore, as can be seen from Figure 12, by increasing the channel number, the accuracy improvement of ComDec in low SNR case (SNR < 20 dB) is much larger than the improvement in high SNR case (SNR > 20 dB). This demonstrates that the channel diversity can help to reduce the location uncertainty that introduced by measurement noises especially in the low SNR case. Moreover, the DR-DFL performs better than LCS in the high SNR case, while performs worse in low SNR cases. This is because, in the low SNR case, the measurement noise can significantly affect the dictionary refinement process. It may lead to a wrong estimation of the dictionary parameters, and the dictionary mismatches will deteriorate the localization performance seriously. Fortunately, in the high SNR case, the DR-DFL can correctly estimate the dictionary parameters and performs better than LCS.

5.5. Localization Error vs. Number of Targets

In this subsection, we examine the performances of multiple CS-based DFL methods when K varies from 2 to 20 at a step of 2. In this simulation, we set the environmental change time

t = 1

to simulate the environmental dynamics in changing environments. With the increase of K, the joint sparsity level of the jointly sparse vectors

S

changes accordingly. This will make the localization accuracies of all CS-based DFL decrease. Figure 13 plots the

A v g . E r r o r

as a function of K. As we can see, the

A v g . E r r o r

of all DFL methods increase as K increases, and the ComDec is more accurate than other DFL methods. Furthermore, we can also observe that the DR-DFL and ComDec perform better than the E-HIPA and LCS. This verifies that the Bayesian recovery algorithm can provide a more accurate sparse reconstruction than the greedy algorithm for CS recovery.

6. Conclusions and Future Work

In this paper, a novel multi-channel CS-based DFL method, ComDec is proposed. It can solve the dictionary mismatch problem caused by environmental dynamics in changing environments. The key novelty of ComDec is the making use of the channel diversity of CSI measurements in CS-based DFL. Moreover, the dictionary parameters of multiple channels are assumed to be adjustable. They can adapt to the changes of environment and different channel characteristics. In this manner, we can avoid the site survey process that is typically adopted in fingerprinting-based DFL and improve the robustness of DFL against environmental dynamics. We formulate the CS-based DFL problem as a joint sparse recovery problem, and develop a novel joint sparse recovery algorithm under the variational Batesian inference framework. Furthermore, four methods are presented to reduce its computational complexity. Finally, numerical simulation results validate the effectiveness of the ComDec method.

As future work, we will further investigate the dictionary mismatch problem in CS-based multi-target DFL caused by more realistic and complex environmental dynamics. Additionally, we intend to design and implement an efficient and accurate DFL framework that can utilize the phase information of CSI measurements.

Author Contributions

All authors contributed extensively to the work in this paper. D.Y., Y.G., and N.L. conceived the proposed scheme, derived inference procedure and conducted numerical simulation. D.Y. analyzed the data and wrote the manuscript. Y.G., N.L. and X.Y. provided valuable comments and contributed to the revision of the manuscript.

Funding

This work was supported in part by the Natural Science Foundation of Jiangsu Province under grant BK20171401; the National Natural Science Foundation of China under grant 61871400 and 61571463.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:

DFL	Device-free localization
VBI	Variational Bayesian inference
LBS	Location-based service
OFDM	Orthogonal Frequency Division Multiplexing
LTE	Long Term Evolution
SNR	Signal-to-noise ratio
CS	Compressive sensing
RIP	Restricted isometry property
RTI	Radio tomographic imaging
GMP	Greedy matching pursuit
LCS	Device-free localization based on compressive sensing
PDF	Probability density function
CG	Conjugate gradient
OMP	Orthogonal matching pursuit
E-HIPA	The energy-efficient framework for high-precision multi-target-adaptive device-free localization
VEM	Variational expectation-maximization
DR-DFL	Dictionary refinement based device-free localization
ComDec	Compressive sensing-based multi-target device-free localization
M-OMP	Multiple measurement vector orthogonal matching pursuit
MMV	Multiple measurement vector
COTS	Commercial off-the-shelf
RSS	Received signal strength
WSN	Wireless sensor networks
CSI	Channel state information
LASSO	Least absolute shrinkage and selection operator
LOS	Line-of-sight
M-SBL	Multiple sparse Bayesian learning
KLD	Kullback-Leibler divergence
M-BMP	Multiple measurement vector basic matching pursuit

References

Youssef, M.; Mah, M.; Agrawala, A. Challenges: Device-Free Passive Localization for Wireless Environments. In Proceedings of the 13th Annual ACM International Conference on Mobile Computing and Networking (MobiCom’07), Montreal, QC, Canada, 9–14 September 2007; pp. 222–229. [Google Scholar] [CrossRef]
Zhang, D.; Ma, J.; Chen, Q.; Ni, L. An RF-Based System for Tracking Transceiver-Free Objects. In Proceedings of the Fifth Annual IEEE International Conference on Pervasive Computing and Communications (PerCom’07), White Plains, NY, USA, 19–23 March 2007; pp. 135–144. [Google Scholar] [CrossRef]
Lei, Q.; Zhang, H.; Sun, H.; Tang, L. Fingerprint-Based Device-Free Localization in Changing Environments using Enhanced Channel Selection and Logistic Regression. IEEE Access 2018, 6, 2569–2577. [Google Scholar] [CrossRef]
Sun, B.; Guo, Y.; Li, N.; Fang, D. Multiple Target Counting and Localization Using Variational Bayesian EM Algorithm in Wireless Sensor Networks. IEEE Trans. Commun. 2017, 65, 2985–2998. [Google Scholar] [CrossRef]
Qian, P.; Guo, Y.; Li, N. Multitarget Localization with Inaccurate Sensor Locations via Variational EM Algorithm. IEEE Sens. Lett. 2019, 3, 1–4. [Google Scholar] [CrossRef]
Guo, Y.; Sun, B.; Li, N.; Fang, D. Variational Bayesian Inference-based Counting and Localization for Off-Grid Targets with Faulty Prior Information in Wireless Sensor Networks. IEEE Trans. Commun. 2018, 66, 1273–1283. [Google Scholar] [CrossRef]
Lin, A.; Ling, H. Doppler and Direction-Of-Arrival (DDOA) Radar for Multiple-Mover Sensing. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 1496–1509. [Google Scholar] [CrossRef]
Krumm, J.; Harris, S.; Meyers, B.; Brumitt, B.; Hale, M.; Shafer, S. Multi-Camera Multi-Person Tracking for EasyLiving. In Proceedings of the IEEE 3rd International Workshop on Visual Surveillance, Dublin, Ireland, 1 July 2000; pp. 3–10. [Google Scholar] [CrossRef]
Want, R.; Hopper, A.; Falco, V.; Gibbons, J. The Active Badge Location System. ACM Trans. Inf. Syst. 1992, 10, 91–102. [Google Scholar] [CrossRef]
Botta, M.; Simek, M. Adaptive Distance Estimation Based on RSSI in 802.15.4 Network. Radioengineering 2013, 22, 1162–1168. [Google Scholar]
Neburka, J.; Tlamse, Z.; Benes, V.; Polak, L.; Kaller, O.; Bolecek, L.; Sebesta, J.; Kratochvil, T. Study of the Performance of RSSI Based Bluetooth Smart Indoor Positioning. In Proceedings of the IEEE 26th International Conference Radioelektronika (RADIOELEKTRONIKA), Kosice, Slovakia, 19–20 April 2016; pp. 121–125. [Google Scholar] [CrossRef]
Zhang, D.; Liu, Y.; Guo, X.; Ni, L.M. RASS: A Real-Time, Accurate, and Scalable System for Tracking Transceiver-Free Objects. IEEE Trans. Parallel Distrib. Syst. 2013, 24, 996–1008. [Google Scholar] [CrossRef]
Mager, B.; Lundrigan, P.; Patwari, N. Fingerprint-Based Device-Free Localization Performance in Changing Environments. IEEE J. Sel. Areas Commun. 2015, 33, 2429–2438. [Google Scholar] [CrossRef]
Yu, D.; Guo, Y.; Li, N.; Fang, D. Dictionary Refinement for Compressive Sensing Based Device-Free Localization via the Variational EM Algorithm. IEEE Access 2016, 4, 9743–9757. [Google Scholar] [CrossRef]
Yang, Z.; Zhou, Z.; Liu, Y. From RSSI to CSI: Indoor Localization via Channel Response. ACM Comput. Surv. 2013, 46, 1–32. [Google Scholar] [CrossRef]
Wang, X.; Gao, L.; Mao, S.; Pandey, S. CSI-Based Fingerprinting for Indoor Localization: A Deep Learning Approach. IEEE Trans. Veh. Technol. 2017, 66, 763–776. [Google Scholar] [CrossRef]
Guo, Y.; Yu, D.; Li, N. Exploiting Fine-Grained Subcarrier Information for Device-Free Localization in Wireless Sensor Networks. Sensors 2018, 18, 3110. [Google Scholar] [CrossRef]
Candes, E.; Wakin, M. An Introduction to Compressive Sampling. IEEE Signal Process. Mag. 2008, 25, 21–30. [Google Scholar] [CrossRef]
Wang, J.; Gao, Q.; Pan, M.; Zhang, X.; Yu, Y.; Wang, H. Towards Accurate Device-Free Wireless Localization with a Saddle Surface Model. IEEE Trans. Veh. Technol. 2016, 65, 6665–6677. [Google Scholar] [CrossRef]
Tzikas, D.G.; Likas, A.C.; Galatsanos, N.P. The Variational Approximation for Bayesian Inference. IEEE Signal Process. Mag. 2008, 25, 131–146. [Google Scholar] [CrossRef]
Wang, J.; Zhang, X.; Gao, Q.; Ma, X.; Feng, X.; Wang, H. Device-Free Simultaneous Wireless Localization and Activity Recognition with Wavelet Feature. IEEE Trans. Veh. Technol. 2017, 66, 1659–1669. [Google Scholar] [CrossRef]
Wang, J.; Fang, D.; Chen, X.; Yang, Z.; Xing, T.; Cai, L. LCS: Compressive Sensing Based Device-Free Localization for Multiple Targets in Sensor Networks. In Proceedings of the IEEE INFOCOM 2013, Turin, Italy, 14–19 April 2013; pp. 145–149. [Google Scholar] [CrossRef]
Wang, J.; Fang, D.; Yang, Z.; Jiang, H.; Chen, X.; Xing, T.; Cai, L. E-HIPA: An Energy-Efficient Framework for High-Precision Multi-Target-Adaptive Device-Free Localization. IEEE Trans. Mob. Comput. 2017, 16, 716–729. [Google Scholar] [CrossRef]
Xiao, J.; Wu, K.; Yi, Y.; Wang, L.; Ni, L. Pilot: Passive Device-Free Indoor Localization Using Channel State Information. In Proceedings of the IEEE 33rd International Conference on Distributed Computing Systems, Philadelphia, PA, USA, 8–11 July 2013; pp. 236–245. [Google Scholar] [CrossRef]
Abdel-Nasser, H.; Samir, R.; Sabek, I.; Youssef, M. Monophy: Mono-Stream-Based Device-Free WLAN Localization via Physical Layer Information. In Proceedings of the IEEE Wireless Communication Network Conference (WCNC), Shanghai, China, 7–10 July 2013; pp. 4546–4551. [Google Scholar] [CrossRef]
Gao, Q.; Wang, J.; Ma, X.; Feng, X.; Wang, H. CSI-Based Device-Free Wireless Localization and Activity Recognition Using Radio Image Features. IEEE Trans. Veh. Technol. 2017, 66, 10346–10356. [Google Scholar] [CrossRef]
Wang, J.; Xiong, J.; Jiang, H.; Jamieson, K.; Chen, X.; Fang, D.; Wang, C. Low Human-Effort, Device-Free Localization with Fine-Grained Subcarrier Information. IEEE Trans. Mob. Comput. 2018, 17, 2550–2563. [Google Scholar] [CrossRef]
Wilson, J.; Patwari, N. Radio Tomographic Imaging with Wireless Networks. IEEE Trans. Mob. Comput. 2010, 9, 621–632. [Google Scholar] [CrossRef]
Wipf, D.; Rao, B.; Nagarajan, S. Latent Variable Bayesian Models for Promoting Sparsity. IEEE Trans. Inf. Theory 2011, 57, 6236–6255. [Google Scholar] [CrossRef]
Wipf, D.; Rao, B. Sparse Bayesian Learning for Basis Selection. IEEE Trans. Signal Process. 2004, 52, 2153–2164. [Google Scholar] [CrossRef]

Figure 1. An example of the CS-based DFL.

Figure 2. Spatial impact area of a wireless link.

Figure 3. Architecture of the ComDec method.

Figure 4. Graphical model for the joint sparse recovery.

Figure 5. Factor graph representation of the joint PDF (22).

Figure 6. Simulation flowchart.

Figure 7. The performance of ComDec when channel number varies from 1 to 25.

Figure 8. The values of the dictionary atoms: (a)

ϕ ._{10}^{5}

; (b)

ϕ_{42}^{12}

, when

t = 1

, 5, and 10.

Figure 8. The values of the dictionary atoms: (a)

ϕ ._{10}^{5}

; (b)

ϕ_{42}^{12}

, when

t = 1

, 5, and 10.

Figure 9. Effect of environmental changes on the average localization error

(A v g . E r r o r)

of multiple DFL approaches.

Figure 9. Effect of environmental changes on the average localization error

(A v g . E r r o r)

of multiple DFL approaches.

Figure 10. Comparison of the original location vector and the reconstructed sparse vectors of different DFL approaches. (a) the original location vector; The reconstructed sparse vector corresponding to (b) E-HIPA; (c) LCS; (d) DR-DFL; (e) ComDec (

F = 5

); (f) ComDec (

F = 15

F = 5

); (f) ComDec (

F = 15

Figure 11. Comparison of the original target positions and the estimated target positions of different DFL approaches.

Figure 12. Average localization error with different SNR.

Figure 13. Average localization error with different K.

Table 1. Default values of simulation parameters.

Parameters	Explains	Values
F	number of channels	15
K	number of targets	5
M	number of links	56
N	number of grids	784
l	link length	14 m
W	iteration number of CG algorithm	17
$α_{th}$	pruning threshold	10
$τ_{max}$	maximum number of iteration	600
$r_{th}$	residual error threshold	$10^{- 4}$
$μ_{th}$	sparsity threshold	$- 20$ dB

Table 2. Computational complexity and localization accuracy.

DFL Method	Sparse Recovery Algorithm	Computational Complexity	$Avg . Error$	$Rms . Error$
ComDec	VBI	$O (N M^{2} F)$	$0.296$ m	$0.612$ m
DR-DFL	VEM	$O (M N^{2})$	$0.981$ m	$1.420$ m
E-HIPA	OMP	$O (K M N + \log (2 N))$	$3.351$ m	$3.509$ m
LCS	GMP	$O (M N)$	$2.188$ m	$2.462$ m

© 2019 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 (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Yu, D.; Guo, Y.; Li, N.; Yang, X. Enhancing the Accuracy and Robustness of a Compressive Sensing Based Device-Free Localization by Exploiting Channel Diversity. Sensors 2019, 19, 1828. https://doi.org/10.3390/s19081828

AMA Style

Yu D, Guo Y, Li N, Yang X. Enhancing the Accuracy and Robustness of a Compressive Sensing Based Device-Free Localization by Exploiting Channel Diversity. Sensors. 2019; 19(8):1828. https://doi.org/10.3390/s19081828

Chicago/Turabian Style

Yu, Dongping, Yan Guo, Ning Li, and Xiaoqin Yang. 2019. "Enhancing the Accuracy and Robustness of a Compressive Sensing Based Device-Free Localization by Exploiting Channel Diversity" Sensors 19, no. 8: 1828. https://doi.org/10.3390/s19081828

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Enhancing the Accuracy and Robustness of a Compressive Sensing Based Device-Free Localization by Exploiting Channel Diversity

Abstract

1. Introduction

2. Related Work

3. Preliminaries and Problem Formulation

3.1. Overview of Multi-Target Device-Free Localization

3.2. CSI Collection and Feature Extraction

3.3. Problem Formulation

4. Target Localization via Variational Bayesian Inference

4.1. Hierarchical Prior Model

4.2. Variational Bayesian Inference

4.3. Joint Sparse Reconstruction

4.4. Target Counting and Localization

5. Numerical Evaluation

5.1. Simulation Setup

5.2. Impact of the Number of Channels

5.3. Effectiveness of ComDec in Changing Environments

5.4. Localization Error vs. SNR

5.5. Localization Error vs. Number of Targets

6. Conclusions and Future Work

Author Contributions

Funding

Conflicts of Interest

Abbreviations

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI