Search

Distributed or decentralized estimation of covariance, and distributed principal component analysis have been introduced and studied in the signal processing community in recent years, and applications in array processing have been indicated in some detail. Inspired by these, this paper provides a detailed development of several distributed algorithms for array processing. New distributed algorithms are proposed for DOA estimation methods like root-MUSIC, total least squares-ESPRIT, and FOCUSS. Other contributions include distributed design of the Capon beamformer from data, and distributed implementation of the spatial smoothing method for coherent sources. A distributed implementation of a recently proposed beamspace method called the convolutional beamspace (CBS) is also proposed. The proposed algorithms are fully distributed – an average consensus (AC) is used to avoid the need for a fusion center. The algorithms are based on a recently reported finite-time version of AC which converges to the exact solution in a finite number of iterations. Numerical examples are given throughout the paper to show the effectiveness of the proposed algorithms. [ABSTRACT FROM AUTHOR]

The main aim of this paper is to extend the semiparametric inference methodology, recently investigated for Real Elliptically Symmetric (RES) distributions, to Complex Elliptically Symmetric (CES) distributions. The generalization to the complex field is of fundamental importance in all practical applications that exploit the complex representation of the acquired data. Moreover, the CES distributions has been widely recognized as a valuable and general model to statistically describe the non-Gaussian behaviour of datasets originated from a wide variety of physical measurement processes. The paper is divided in two parts. In the first part, a closed form expression of the constrained semiparametric Cramér-Rao Bound (CSCRB) for the joint estimation of complex mean vector and complex scatter matrix of a set of CES-distributed random vectors is obtained by exploiting the so-called Wirtinger or $\mathbb {C}\mathbb {R}$ -calculus. The second part deals with the derivation of the semiparametric version of the Slepian-Bangs formula in the context of the CES model. Specifically, the proposed semiparametric Slepian-Bangs (SSB) formula provides us with a useful and ready-to-use expression of the semiparametric Fisher Information Matrix (SFIM) for the estimation of a parameter vector parametrizing the complex mean and the complex scatter matrix of a CES-distributed vector in the presence of unknown, nuisance, density generator. Furthermore, we show how to exploit the derived SSB formula to obtain the semiparametric counterpart of the Stochastic CRB for Direction of Arrival (DOA) estimation under a random signal model assumption. Simulation results are also provided to clarify the theoretical findings and to demonstrate their usefulness in common array processing applications. [ABSTRACT FROM AUTHOR]

This paper utilizes sparse array motion to increase the numbers of achievable both degrees of freedom (DOFs) and consecutive lags in direction-of-arrival (DOA) estimation problems. We use commonly employed environment-independent sparse array configurations. The design of these arrays is not dependent on the sources in the field of view, but rather aims at achieving desirable difference co-arrays. They include structured coprime and nested arrays, minimum redundancy array (MRA), minimum hole array (MHA), and sparse uniform linear array (SULA). Array motion can fill the holes in the spatial autocorrelation lags associated with a fixed platform and, therefore, increases the number of sources detectable by the same number of array sensors. Quasi-stationarity of the environment is assumed where the source locations and waveforms are considered invariant over array motion of half wavelength. Closed-form expressions of the number of DOFs and consecutive spatial correlation lags for coprime and nested arrays as well as SULA, due to array translation motion, are derived. The number of DOFs and consecutive lags for the specific cases of MRA an 5 avaluated. We show the respective DOA estimation performance based on sparse reconstruction techniques. [ABSTRACT FROM AUTHOR]

For stationary signals in time, the weak law of large numbers (WLLN) states that ensemble and realization averages are within $\epsilon$ of each other with a probability of order ${\mathcal O}(1/N\epsilon ^2)$ when considering $N$ signal components. The graph WLLN introduced in this paper shows that the same is essentially true for signals supported on graphs. However, the notions of stationarity, ensemble mean, and realization mean are different. Recent papers have defined graph stationary signals as those that satisfy a form of invariance with respect to graph diffusion. The ensemble mean of a graph stationary signal is not a constant but a node-varying signal whose structure depends on the spectral properties of the graph. The realization average of a graph signal is defined here as an average of successive weighted averages of local signal values with signal values of neighboring nodes. The graph WLLN shows that these two node-varying signals are within $\epsilon$ of each other with probability of order ${\mathcal O}(1/N\epsilon ^2)$ in at least some nodes. In stationary time signals, the realization average is not only a consistent estimator of the ensemble mean but also optimal in terms of mean squared error (MSE). This is not true for graph signals. Optimal MSE graph filter designs are also presented. An example problem concerning the estimation of the mean of a Gaussian random field is presented. [ABSTRACT FROM AUTHOR]

Effective non-parametric density estimation is a key challenge in high-dimensional multivariate data analysis. In this paper, we propose a novel approach that builds upon tensor factorization tools. Any multivariate density can be represented by its characteristic function, via the Fourier transform. If the sought density is compactly supported, then its characteristic function can be approximated, within controllable error, by a finite tensor of leading Fourier coefficients, whose size depends on the smoothness of the underlying density. This tensor can be naturally estimated from observed and possibly incomplete realizations of the random vector of interest, via sample averaging. In order to circumvent the curse of dimensionality, we introduce a low-rank model of this characteristic tensor, which significantly improves the density estimate especially for high-dimensional data and/or in the sample-starved regime. By virtue of uniqueness of low-rank tensor decomposition, under certain conditions, our method enables learning the true data-generating distribution. We demonstrate the very promising performance of the proposed method using several toy, measured, and image datasets. [ABSTRACT FROM AUTHOR]

The state estimation problem of heterogeneous multi-agent systems with random transport protocol is investigated in this paper. Due to the dependency of the agent dynamics and the random sparse structure induced by the random transport protocol, the optimal state estimation design becomes complex and challenging. An optimal state estimator is successfully designed by the Hadamard product and gradient method. Based on the analysis of the matrix functions, a sufficient condition is established to guarantee that the average estimate error covariance is limited. Finally, a numerical example and a smart grid model are utilized to demonstrate the effectiveness of the deigned estimator. [ABSTRACT FROM AUTHOR]

This paper focuses on optimal time-of-arrival (TOA) sensor placement for multiple target localization simultaneously. In previous work, different solutions only using non-shared sensors to localize multiple targets have been developed. Those methods localize different targets one-by-one or use a large number of mobile sensors with many limitations, such as low effectiveness and high network complexity. In this paper, firstly, a novel optimization model for multi-target localization incorporating shared sensors is formulated. Secondly, the systematic theoretical results of the optimal sensor placement are derived and concluded using the A-optimality criterion, i.e., minimizing the trace of the inverse Fisher information matrix (FIM), based on rigorous geometrical derivations. The reachable optimal trace of Cramér-Rao lower bound (CRLB) is also derived. It can provide optimal conditions for many cases and even closed form solutions for some special cases. Thirdly, a novel numerical optimization algorithm to quickly find and calculate the (sub-)optimal placement and achievable lower bound is explored, when the model becomes complicated with more practical constraints. Then, a hybrid method for solving the most general situation, integrating both the analytical and numerical solutions, is proposed. Finally, the correctness and effectiveness of the proposed theoretical and mathematical methods are demonstrated by several simulation examples. [ABSTRACT FROM AUTHOR]

The problem of range estimation based on noisy measurements with phase ambiguity can be solved using signals across multiple frequencies by existing algorithms based on either the Chinese Remainder Theorem or lattice theory. The performances of these algorithms are constrained by the trade-off between SNR of measurements and required estimation accuracy. In this paper, we firstly propose an algorithm to find a set of wavelengths satisfying the conditions to apply the closed-form lattice algorithm from a given interval. Then we propose an enhanced lattice algorithm for range estimation based on existing lattice approaches. In the conventional lattice algorithm, the selection of an integer vector by the estimator often fails to be the ground truth because of noise. As a result, the true value can be a neighbouring lattice point of the estimated point where neighbours are defined in terms of a distance depending on the noise characteristics. In this paper, all neighbouring lattice points of the estimated point obtained by a conventional lattice algorithm are regarded as candidates for the new estimate. Then an efficient estimator is presented to choose one of those candidates to be the new estimate. The performance of the proposed algorithm, i.e. the successful probability of reconstruction, is derived mathematically. Additionally, a computation-reduced algorithm is presented to reduce the overhead of computation. The calculated probability and simulation results demonstrate the efficiency of the proposed algorithm. [ABSTRACT FROM AUTHOR]

This paper proposes a low complexity n-dimensional (nD) FastICA algorithm and architecture by introducing the concept of coordinate rotation where n\geq 2. The proposed algorithm can merge the two key steps of conventional FastICA—preprocessing and update and is therefore capable of reducing the hardware complexity of the conventional FastICA significantly as demonstrated in this paper. Hardware implementation can further be simplified due to the recursive nature of the proposed algorithm where the same 2D hardware module can be used as the fundamental core to implement nD architecture. Together with the algorithm formulation, its functionality is also validated and hardware complexity is analyzed and compared with the conventional nD FastICA. [ABSTRACT FROM PUBLISHER]

The conventional direction of arrival (DOA) estimation methods are performed with multiple receiving channels. In this paper, a changeling DOA estimation problem is addressed in a different scenario with only one full-functional receiving channel. A new unmanned aerial vehicle (UAV) swarm system using multiple lifted reconfigurable intelligent surface (RIS) is proposed for the DOA estimation. The UAV movement degrades the DOA estimation performance significantly, and the existing atomic norm minimization (ANM) methods cannot be used in the scenario with array perturbation. Specifically, considering the position perturbation of UAVs, a new atomic norm-based DOA estimation method is proposed, where an atomic norm is defined with the parameter of the position perturbation. Then, a customized semi-definite programming (SDP) method is derived to solve the atomic norm-based method, where different from the traditional SDP method, an additional transforming matrix is formulated. Moreover, a gradient descent method is applied to refine the estimated DOA and the position perturbation further. Simulation results show that the proposed method achieves much better DOA estimation performance in the RIS-aided UAV swarm system with only one receiving channel than various benchmark schemes. [ABSTRACT FROM AUTHOR]

Random matrices are widely used to estimate the extent of an elliptically contoured object. Usually, it is assumed that the measurements follow a normal distribution, with its standard deviation being proportional to the object’s extent. However, the random matrix approach can filter the center of gravity and the covariance matrix of measurements independently of the measurement model. This work considers the whole chain from data acquisition to the linear Kalman Filter with extension estimation as a reference plant. The input is the (unknown) ground truth (position and extent). The output is the filtered center of gravity and the filtered covariance matrix of the measurement distribution. A virtual measurement model emulates the behavior of the reference plant. The input of the virtual measurement model is adapted using the proposed algorithm until the output parameters of the virtual measurement model match the result of the reference plant. After the adaptation, the input to the virtual measurement model is considered an estimation for position and extent. The main contribution of this paper is the reference model concept and an adaptation algorithm to optimize the input of the virtual measurement model. [ABSTRACT FROM AUTHOR]

In this paper, we propose multi-input multi-output (MIMO) beamforming designs towards joint radar sensing and multi-user communications. We employ the Cramér-Rao bound (CRB) as a performance metric of target estimation, under both point and extended target scenarios. We then propose minimizing the CRB of radar sensing while guaranteeing a pre-defined level of signal-to-interference-plus-noise ratio (SINR) for each communication user. For the single-user scenario, we derive a closed form for the optimal solution for both cases of point and extended targets. For the multi-user scenario, we show that both problems can be relaxed into semidefinite programming by using the semidefinite relaxation approach, and prove that the global optimum can be generally obtained. Finally, we demonstrate numerically that the globally optimal solutions are reachable via the proposed methods, which provide significant gains in target estimation performance over state-of-the-art benchmarks. [ABSTRACT FROM AUTHOR]

This paper concerns the sparse Bayesian learning (SBL) problem for group sparse signals. Group sparsity means that the signal coefficients can be divided into groups and that the entries in one group are simultaneously zero or nonzero. In SBL, each group is controlled by a hyperparameter, which is estimated by solving the marginal likelihood maximization (MLM) problem. MLM is used to maximize the marginal likelihood of a given hyper-parameter by fixing all others. The main contribution of this paper is to extend the fast marginal likelihood maximization (FMLM) method to the group sparsity case. The key process of realizing FMLM in the group sparsity case is solved by finding the roots of a polynomial. Hence, the local maximum of the marginal likelihood can be found exactly, and the final estimation is obtained more efficiently. Furthermore, most large matrix inverses involved in MLM are replaced with the eigenvalue decompositions of much smaller matrices, which substantially reduces the computational complexity. Numerical results with both simulated and real data demonstrate the excellent performance and relatively high computational efficiency of the proposed method. [ABSTRACT FROM PUBLISHER]

This paper proposes a novel proximal-gradient algorithm for a decentralized optimization problem with a composite objective containing smooth and nonsmooth terms. Specifically, the smooth and nonsmooth terms are dealt with by gradient and proximal updates, respectively. The proposed algorithm is closely related to a previous algorithm, PG-EXTRA (W. Shi, Q. Ling, G. Wu, and W. Yin, “A proximal gradient algorithm for decentralized composite optimization,” IEEE Trans. Signal Process., vol. 63, no. 22, pp. 6013–6023, 2015), but has a few advantages. First of all, agents use uncoordinated step-sizes, and the stable upper bounds on step-sizes are independent of network topologies. The step-sizes depend on local objective functions, and they can be as large as those of the gradient descent. Second, for the special case without nonsmooth terms, linear convergence can be achieved under the strong convexity assumption. The dependence of the convergence rate on the objective functions and the network are separated, and the convergence rate of the new algorithm is as good as one of the two convergence rates that match the typical rates for the general gradient descent and the consensus averaging. We provide numerical experiments to demonstrate the efficacy of the introduced algorithm and validate our theoretical discoveries. [ABSTRACT FROM AUTHOR]

This paper studies linear distributed estimation of an unknown random parameter vector in a bandwidth-constrained multisensor network. To meet the bandwidth limitations, each sensor converts its observation into a low-dimensional datum via a suitable linear transformation. Then, the fusion center estimates the parameter vector by linearly combining all the received low-dimensional data, aiming at minimizing the estimation mean square error. The main purpose of this paper is to jointly determine the compression dimension of each sensor (referred to as dimension assignment) and design the corresponding compression matrix when the total compression dimensions is limited. Such a joint design problem can be formulated as a rank-constrained optimization problem and it is shown to be NP-hard for the first time. In addition, successive quadratic upper-bound minimization (SQUM), SQUM-block coordinate descent (SQUM-BCD) and nuclear norm regularization (NNR) methods are developed to solve it approximately. Furthermore, we show that any accumulation point of the sequence generated by the SQUM method satisfies the Karush-Kuhn-Tucker conditions of the rank-constrained optimization problem, and the Phase II algorithm of the SQUM-BCD and NNR methods (both are two-phase algorithms and have the same Phase II algorithm) guarantees convergence at least to a stationary point. Numerical experiments illustrate the advantages of the proposed methods compared with the existing method. [ABSTRACT FROM AUTHOR]

In the practical radar with multiple antennas, the antenna imperfections degrade the system performance. In this paper, the problem of estimating the direction of arrival (DOA) in a multiple-input and multiple-output (MIMO) radar system with unknown mutual coupling effect between antennas is investigated. To exploit the target sparsity in the spatial domain, the compressed sensing based methods have been proposed by discretizing the detection area and formulating the dictionary matrix, so an off-grid gap is caused by the discretization processes. In this paper, different from the present DOA estimation methods, both the off-grid gap due to the sparse sampling and the unknown mutual coupling effect between antennas are considered at the same time, and a novel sparse system model for DOA estimation is formulated. Then, a novel sparse Bayesian learning (SBL)-based method named sparse Bayesian learning with the mutual coupling (SBLMC) is proposed, where an expectation-maximum-based method is established to estimate all the unknown parameters including the noise variance, the mutual coupling vectors, the off-grid vector, and the variance vector of scattering coefficients. Additionally, the prior distributions for all the unknown parameters are theoretically derived. With regard to the DOA estimation performance, the proposed SBLMC method can outperform state-of-the-art methods in the MIMO radar with unknown mutual coupling effect, while keeping the acceptable computational complexity. [ABSTRACT FROM AUTHOR]

This paper has a twofold goal. The first aim is to provide a deeper understanding of the family of the real elliptically symmetric (RES) distributions by investigating their intrinsic semiparametric nature. The second aim is to derive a semiparametric lower bound for the estimation of the parametric component of the model. The RES distributions represent a semiparametric model, where the parametric part is given by the mean vector and by the scatter matrix, while the non-parametric, infinite-dimensional, part is represented by the density generator. Since, in practical applications, we are often interested only in the estimation of the parametric component, the density generator can be considered as nuisance. The first part of this paper is dedicated to conveniently place the RES distributions in the framework of the semiparametric group models. In the second part, building on the mathematical tools previously introduced, the constrained semiparametric Cramér–Rao bound (CSCRB) for the estimation of the mean vector and of the constrained scatter matrix of a RES distributed random vector is introduced. The CSCRB provides a lower bound on the mean squared error of any robust $M$ -estimator of mean vector and scatter matrix when no a priori information on the density generator is available. A closed-form expression for the CSCRB is derived. Finally, in simulations, we assess the statistical efficiency of the Tyler's and Huber's scatter matrix $M$ -estimators with respect to the CSCRB. [ABSTRACT FROM AUTHOR]

Frequency agile radar (FAR) is known to have excellent electronic counter-countermeasures performance and the potential to realize spectrum sharing in dense electromagnetic environments. Many compressed sensing (CS) based algorithms have been developed for joint range and Doppler estimation in FAR. This paper considers theoretical analysis of FAR via CS algorithms. In particular, we analyze the properties of the sensing matrix, which is a highly structured random matrix. We then derive bounds on the number of recoverable targets. Numerical simulations and field experiments validate the theoretical findings and demonstrate the effectiveness of CS approaches to FAR. [ABSTRACT FROM AUTHOR]

Linear estimation of signals is often based on covariance matrices estimated from training, which can perform poorly if the training data are limited and the estimated covariance matrices are ill-conditioned. Shrinking the covariance matrix toward a scaled identity matrix can improve the robustness against the model uncertainty provided the shrinkage factor is appropriately chosen. This paper introduces several cross-validation schemes for choosing the shrinkage factors in applications where the covariance matrices are replaced with sample covariance matrices or constructed from least squares estimates of the linear model parameters. For cases where the training and out-of-training data are identically distributed, we derive leave-one-out cross-validation (LOOCV) schemes that repeatedly split the training data with respect to time to determine the optimal shrinkage factors for either model fitting or signal estimation. For cases where they are distributed differently, we develop alternative LOOCV schemes that repeatedly split the out-of-training observations with respect to space. We derive computationally efficient implementations of those schemes and provide examples to demonstrate their performance. [ABSTRACT FROM AUTHOR]

In this paper, we address the problem of distributed adaptive estimation of node-specific signals for signal enhancement or noise reduction in wireless sensor networks with multi-sensor nodes. The estimation is performed by a multi-channel Wiener filter (MWF) in which a low-rank approximation based on a generalized eigenvalue decomposition (GEVD) is incorporated. In non-stationary or low-SNR conditions, this GEVD-based MWF has been demonstrated to be more robust than the original MWF. In a centralized realization where a fusion center has access to all the nodes’ sensor signal observations, the network-wide sensor signal correlation matrices and the low-rank approximation can be directly estimated and used to compute the network-wide GEVD-based MWF. However, in this paper, we aim to avoid centralizing the sensor signal observations, in which case the network-wide sensor signal correlation matrices cannot be estimated. To this end, we start from the so-called distributed adaptive node-specific signal estimation (DANSE) algorithm, and include GEVD-based low-rank approximations in the per-node local computations. Remarkably, the new algorithm is able to significantly compress the signal observations transmitted between the nodes, while still converging to the network-wide GEVD-based MWF as if each node would have access to all sensor signal observations, even though the low-rank approximations are applied locally at each node. We provide a theoretical convergence analysis, which shows that the algorithm converges to the network-wide GEVD-based MWF under conditions that are less strict than in the original DANSE algorithm. The convergence and performance of the algorithm are further investigated via numerical simulations. [ABSTRACT FROM AUTHOR]

The estimation of mutual information (MI) or conditional mutual information (CMI) from a set of samples is a long-standing problem. A recent line of work in this area has leveraged the approximation power of artificial neural networks and has shown improvements over conventional methods. One important challenge in this new approach is the need to obtain, given the original dataset, a different set where the samples are distributed according to a specific product density function. This is particularly challenging when estimating CMI. In this paper, we introduce a new technique, based on k nearest neighbors (k-NN), to perform the resampling and derive high-confidence concentration bounds for the sample average. Then the technique is employed to train a neural network classifier and the CMI is estimated accordingly. We propose three estimators using this technique and prove their consistency, make a comparison between them and similar approaches in the literature, and experimentally show improvements in estimating the CMI in terms of accuracy and variance of the estimators. [ABSTRACT FROM AUTHOR]

Inferring network topologies from observed graph-structured data (also known as graph signals) is a crucial task in many applications of network science. Existing papers on network topology inference typically assume that the observations at all nodes are available. However, there are many situations where only partial observations can be collected due to application-specific constraints. To handle the missing data problem, we propose a framework that relies on heterogeneous incomplete data from a collection of related networks to identify multiple network topologies simultaneously. This work advocates the Gaussian graphical model (GGM) and casts the topology inference problem in terms of estimating the precision matrix that has a form of graph Laplacian. Firstly, an unbiased estimator for the covariance matrix of incomplete data is established and then algorithms based on the alternating direction method of multipliers (ADMM) are developed to jointly estimate graph topologies, rather than estimate each graph topology separately. So that we can borrow information across the multiple related networks to eliminate the impact of missing data. Moreover, non-asymptotic statistical analysis is provided, which proves the consistency of the graph estimator and enables us to investigate the effect of several key factors on the graph estimation error bound. Furthermore, based on the consistent graph estimator, an adaptive algorithm that utilizes the reweighting scheme is proposed to improve the estimation accuracy of the graph-edge structure. Finally, we evaluate our method on both real and synthetic datasets, and the experimental results demonstrate the advantage of our method in comparison with benchmarking algorithms. [ABSTRACT FROM AUTHOR]

The paper presents a communication efficient distributed algorithm, $\mathcal {CIRFE}$ of the consensus+innovations type, to estimate a high-dimensional parameter in a multi-agent network, in which each agent is interested in reconstructing only a few components of the parameter. This problem arises, for example, when monitoring the high-dimensional distributed state of a large-scale infrastructure with a network of limited capability sensors and where each sensor is tasked with estimating some local components of the state. At each observation sampling epoch, each agent updates its local estimate of the parameter components in its interest set by simultaneously processing the latest locally sensed information (innovations) and the parameter estimates from agents (consensus) in its communication neighborhood given by a time-varying possibly sparse graph. Under minimal conditions, on the interagent communication network and the sensing models, almost sure convergence of the estimate sequence at each agent to the components of the true parameter in its interest set is established. Furthermore, the paper establishes the performance of $\mathcal {CIRFE}$ in terms of asymptotic covariance of the estimate sequences and specifically characterizes the dependencies of the component wise asymptotic covariance in terms of the number of agents tasked with estimating it. Finally, simulation experiments demonstrate the efficacy of $\mathcal {CIRFE}$. [ABSTRACT FROM AUTHOR]

Detecting periodicity in a short sequence is an important problem, with many applications across science and engineering. Several efficient algorithms have been proposed for this over the years. There is a wide choice available today in terms of the tradeoff between algorithmic complexity and estimation accuracy. In spite of such a rich history, one particular aspect of period estimation has received very little attention from a fundamental perspective. Namely, given a discrete time periodic signal and a list of candidate integer periods, what is the absolute minimum datalength required to estimate its integer period? Notice that the answer we seek must be a fundamental bound, i.e., independent of any particular period estimation technique. Common intuition suggests the minimum datalength as twice the largest expected period. However, this is true only under some special contexts. This paper derives the exact necessary and sufficient bounds to this problem. The above-mentioned question is also extended to the case of mixtures of periodic signals. First, a careful mathematical formulation discussing the unique identifiability of the component periods (hidden integer periods) is presented. Once again, a rigorous theoretical framework in this regard is missing in the existing literature but is a necessary platform to derive precise bounds on the minimum necessary datalength. The bounds given here are generic, that is, independent of the algorithms used. Specific algorithm-dependent bounds are also presented in the end for the case of dictionary-based integer period estimation reported in recent years. [ABSTRACT FROM AUTHOR]

Recent papers have formulated the problem of learning graphs from data as an inverse covariance estimation problem with graph Laplacian constraints. While such problems are convex, existing methods cannot guarantee that solutions will have specific graph topology properties (e.g., being a tree), which are desirable for some applications. The problem of learning a graph with topology properties is in general non-convex. In this paper, we propose an approach to solve these problems by decomposing them into two sub-problems for which efficient solutions are known. Specifically, a graph topology inference (GTI) step is employed to select a feasible graph topology. Then, a graph weight estimation (GWE) step is performed by solving a generalized graph Laplacian estimation problem, where edges are constrained by the topology found in the GTI step. Our main result is a bound on the error of the GWE step as a function of the error in the GTI step. This error bound indicates that the GTI step should be solved using an algorithm that approximates the data similarity matrix by another matrix whose entries have been thresholded to zero to have the desired type of graph topology. The GTI stage can leverage existing methods, which are typically based on minimizing the total weight of removed edges. Since the GWE stage is an inverse covariance estimation problem with linear constraints, it can be solved using existing convex optimization methods. We demonstrate that our approach can achieve good results for both synthetic and texture image data. [ABSTRACT FROM AUTHOR]

A coprime sensor array (CSA) interleaves two undersampled uniform linear arrays with coprime undersampling factors and has recently found broad applications in signal detection and estimation. CSAs commonly use the product processor by multiplying the scanned responses of two colinear subarrays to estimate the spatial power spectral density (PSD) of the received signal. This paper proposes a new CSA processor, the CSAmin processor, which chooses the minimum between the two CSA subarray periodograms at each bearing to estimate the spatial PSD. The proposed CSAmin processor resolves the CSA subarray spatial aliasing equally well as the product processor. For an extended aperture CSA, the CSAmin reduces the peak sidelobe height and total sidelobe area over the product processor for the same CSA geometry. Moreover, unlike the PSD estimate from the product processor, the PSD estimate from the CSAmin is guaranteed to be positive semidefinite. This paper derives the probability density function, the complementary cumulative distribution function (CCDF, or tail distribution), and the first two moments of the CSAmin PSD estimator in closed form for Gaussian sources in white Gaussian noise. Numerical simulations verify the derived CSAmin statistics and demonstrate that the CSAmin improves the performance over the product processor in detecting narrowband Gaussian sources in the presence of loud interferers and noise. The CSAmin spatial PSD estimate achieves lower variance than the product processor estimate, and keeps the PSD estimate unbiased for white Gaussian processes and asymptotically unbiased for nonwhite Gaussian processes. [ABSTRACT FROM AUTHOR]

This paper studies distributed estimation of unstable dynamic random fields observed by a sparsely connected network of sensors. The field dynamics are globally detectable, but not necessarily locally detectable. We propose a consensus+innovations distributed estimator, termed Distributed Information Kalman Filter. We prove under what conditions this estimator is asymptotically unbiased with bounded mean-squared error, smaller than for other alternative distributed estimators. Monte Carlo simulations confirm our theoretical error asymptotic results. [ABSTRACT FROM PUBLISHER]

In this paper, we consider the problem of distributed inference in tree based networks. In the framework considered in this paper, distributed nodes make a 1-bit local decision regarding a phenomenon before sending it to the fusion center (FC) via intermediate nodes. We propose the use of coding theory based techniques to solve this distributed inference problem in such structures. Data is progressively compressed as it moves towards the FC. The FC makes the global inference after receiving data from intermediate nodes. Data fusion at nodes as well as at the FC is implemented via error correcting codes. In this context, we analyze the performance for a given code matrix and also design the optimal code matrices at every level of the tree. We address the problems of distributed classification and distributed estimation separately and develop schemes to perform these tasks in tree networks. The proposed schemes are of practical significance due to their simple structure. We study the asymptotic inference performance of our schemes for two different classes of tree networks: fixed height tree networks, and fixed degree tree networks. We derive the sufficient conditions under which the proposed schemes are asymptotically optimal. [ABSTRACT FROM PUBLISHER]

Carrier frequency offset (CFO) introduces spectrum misalignment of the transmit and receive filters, causing energy loss and distortion of received signal. A recent study presents an accurate signal model that incorporates these characteristics. However, CFO estimators in the existing literature were developed based on the signal models without these characteristics. This paper generalizes the accurate signal model with CFO to include timing offset and sampling offset and then studies what effects the accurate signal model bring to the CFO estimation and how to address them. We investigate different existing CFO estimators and find some of them lose their optimality under the accurate signal model. We also develop two new integer normalized CFO (ICFO) estimators under the accurate signal model using only one OFDM symbol; one of them is a pilot-aided estimator using both pilot and data and the other one is a blind estimator using data only. Compared with the existing approaches, both of the new estimators make better use of channel and energy loss information, thus yielding more accurate estimation. Furthermore, the proposed methods overcome the existing approaches’ limitation to phase shift keying modulation. Simulation results corroborate the advantages of the new estimators. Analytical performance results are also provided for the proposed estimators. [ABSTRACT FROM PUBLISHER]

An analytic solution to the multi-target Bayes recursion known as the \delta -Generalized Labeled Multi-Bernoulli (\delta -GLMB) filter has been recently proposed by Vo and Vo in [“Labeled Random Finite Sets and Multi-Object Conjugate Priors,” IEEE Trans. Signal Process., vol. 61, no. 13, pp. 3460–3475, 2014]. As a sequel to that paper, the present paper details efficient implementations of the \delta -GLMB multi-target tracking filter. Each iteration of this filter involves an update operation and a prediction operation, both of which result in weighted sums of multi-target exponentials with intractably large number of terms. To truncate these sums, the ranked assignment and K-th shortest path algorithms are used in the update and prediction, respectively, to determine the most significant terms without exhaustively computing all of the terms. In addition, using tools derived from the same framework, such as probability hypothesis density filtering, we present inexpensive (relative to the \delta -GLMB filter) look-ahead strategies to reduce the number of computations. Characterization of the L1-error in the multi-target density arising from the truncation is presented. [ABSTRACT FROM PUBLISHER]

This paper deals with the nonlinear minimum mean-squared error estimation problem by using a quaternion widely linear model. On the basis of the information supplied by a Gaussian signal and its square, a quaternion observation process is defined and then, by applying a widely linear processing, an optimal estimator for the continuous-time setting is provided. The special structure of the estimator proves its superiority over the complex-valued widely linear solution. The continuous-discrete version of the problem is also studied where the solution takes the form of a suboptimal estimator useful in practical applications. In addition, the particular case of signal plus noise is considered in which the suboptimal solution and its associated error can be implemented through an iterative algorithm. Two numerical simulation examples are presented showing the advantages of the proposed approach. [ABSTRACT FROM PUBLISHER]

We correct errors in step 6 of Algorithm 1 given in the paper ‘'Sparse-Group Lasso for Graph Learning from Multi-Attribute Data.’' Other parts of the paper including numerical results are unchanged as the algorithm was implemented correctly. [ABSTRACT FROM AUTHOR]

Covariance matrix estimation is a crucial task in adaptive signal processing applied to several surveillance systems, including radar and sonar. In this paper we propose a dynamic learning strategy to track both the covariance matrix of data and its structure (class). We assume that, given the class, the posterior distribution of the covariance is described through a mixture of inverse Wishart distributions, while the class evolves according to a Markov chain. Hence, we devise a novel and general filtering strategy, called multi-class inverse Wishart mixture filter, able to capitalize on previous observations so as to accurately track and estimate the covariance. Some case studies are provided to highlight the effectiveness of the proposed technique, which is shown to outperform alternative methods in terms of both covariance estimation accuracy and probability of correct model selection. Specifically, the proposed filter is compared with class-clairvoyant covariance estimators, e.g., the maximum likelihood and the knowledge-based recursive least square filter, and with the model order selection method based on the Bayesian information criterion. [ABSTRACT FROM AUTHOR]

This paper studies distributed diffusion adaptation over clustered multi-task networks in the presence of impulsive interferences and Byzantine attacks. We develop a robust resilient diffusion least mean Geman-McClure-estimation (RDLMG) algorithm based on the cost function used by the Geman-McClure estimator, which can reduce the sensitivity to large outliers and make the algorithm robust under impulsive interferences. Moreover, the mean sub-sequence reduced method, in which each node discards the extreme value information of cost contributions received from its neighbors, can make the network resilient against Byzantine attacks. In this regard, the proposed RDLMG algorithm ensures that all normal nodes converge to their ideal states with cooperation among nodes. A statistical analysis of the RDLMG algorithm is also carried out in terms of mean and mean-square performances. Numerical results evaluate the proposed RDLMG algorithm in applications to multi-target localization and multi-task spectrum sensing. [ABSTRACT FROM AUTHOR]

The recent criteria of a preferable linear sparse array for a robust direction of arrival (DOA) estimation are: the closed-form expression for its sensor locations; the large central uniform linear array (ULA) segment in its resulting co-array; and the fewer sensor pairs with small separations in its arrangement. These features, in turn, will lead to a straightforward array design process, more sources being resolved and strong helpful performance for the array within a real-world mutual coupling environment. This paper introduces a new array called the dilated nested array (DNA), which takes into account all these considerations. The sensors of the DNA are arranged in a way that yields an effective difference co-array with a central virtual ULA part larger than that of the co-prime array and equal to that of the nested array, and in an equivalent importance, a physical arrangement with smaller numbers of sensor pairs with small separations than those of the super nested array, so that the array can notably reduce the mutual coupling effects in comparison. A displaced arrangement obtained from the DNA geometry is also introduced to further increase the co-array total number of DOFs and at the same time to decrease the numbers of sensor pairs with small separations. Extensive numerical simulations are included to demonstrate the effectiveness of the proposed DNAs over other sparse arrays. [ABSTRACT FROM AUTHOR]

In this paper, we address the problem of joint direction-of-arrival (DoA) and range estimation using frequency diverse coprime array (FDCA). By incorporating the coprime array structure and coprime frequency offsets, a two-dimensional space-frequency virtual difference coarray corresponding to uniform array and uniform frequency offset is considered to increase the number of degrees-of-freedom (DoFs). However, the reconstruction of the doubly-Toeplitz covariance matrix is computationally prohibitive. To solve this problem, we propose an interpolation algorithm based on decoupled atomic norm minimization (DANM), which converts the coarray signal to a simple matrix form. On this basis, a relaxation-based optimization problem is formulated to achieve joint DoA-range estimation with enhanced DoFs. The reconstructed coarray signal enables application of existing subspace-based spectral estimation methods. The proposed DANM problem is further reformulated as an equivalent rank-minimization problem which is solved by cyclic rank minimization. This approach avoids the approximation errors introduced in nuclear norm-based approach, thereby achieving superior root-mean-square error which is closer to the Cramér-Rao bound. The effectiveness of the proposed method is confirmed by theoretical analyses and numerical simulations. [ABSTRACT FROM AUTHOR]

In this paper, we propose a federated online training framework for deep neural network (DNN)-based channel estimation (CE) in frequency-division duplexing (FDD) multi-user massive multiple-input multiple-output (MU-MIMO) systems. The proposed DNN consists of two stages, where Stage-I explores the partial common sparsity structure among the users, while Stage-II estimates the channel state information (CSI) with reduced pilot overhead leveraging the common support information provided by Stage-I. To realize online training, we first propose three axioms for a legitimate online loss function, based on which we develop the federated online training algorithm with convergence analysis. The proposed two-tier DNN is trained online at the base station (BS) based on real-time pilot measurements distributively fed back from the users without the need of true channel labels, and the estimates for the CSI at the transmitter (CSIT) can be simultaneously generated in real-time. Meanwhile, the weights of Stage-II can be broadcasted to the users for real-time estimation of the CSI at the receiver (CSIR) at each user. Simulation shows that the proposed solution achieves higher CE performance than traditional compressive sensing (CS)-based algorithms while enjoying much faster computation. The solution is also robust to the channel model mismatches induced by the change of propagation environment, and is able to track the time-varying channel model. [ABSTRACT FROM AUTHOR]

Covariance matrix tapers have a long history in signal processing and related fields. Examples of applications include autoregressive models (promoting a banded structure) or beamforming (widening the spectral null width associated with an interferer). In this paper, the focus is on high-dimensional setting where the dimension $p$ is high, while the data aspect ratio $n/p$ is low. We propose an estimator called Tabasco (TApered or BAnded Shrinkage COvariance matrix) that shrinks the tapered sample covariance matrix towards a scaled identity matrix. We derive optimal and estimated (data adaptive) regularization parameters that are designed to minimize the mean squared error (MSE) between the proposed shrinkage estimator and the true covariance matrix. These parameters are derived under the general assumption that the data is sampled from an unspecified elliptically symmetric distribution with finite 4th order moments (both real- and complex-valued cases are addressed). Simulation studies show that the proposed Tabasco outperforms all competing tapering covariance matrix estimators in diverse setups. An application to space-time adaptive processing (STAP) also illustrates the benefit of the proposed estimator in a practical signal processing setup. [ABSTRACT FROM AUTHOR]

In this paper, we propose and implement blind multiuser carrier frequency offsets (CFOs) estimation for multi-input multi-output (MIMO) single carrier frequency division multiple access (SC-FDMA) uplink system using orthogonal propagator method (OPM) over a noisy frequency-selective channel. The performance of typical OPM is limited at low signal-to-noise ratio (SNR) due to its noise-free signal model. The additive white Gaussian noise mainly impacts the diagonal elements of the data covariance matrix (DCM) of the received signal. First, the proposed fractional CFOs estimator minimizes the noise effect by reconstructing the denoised diagonal elements of the DCM for the received signal through cubic Hermite interpolation. Therefore, we use modified OPM to estimate multiple CFOs for MIMO SC-FDMA systems. Second, we have used an iterative method based on the first-order Taylor series approximation of the CFO vector to avoid peak ambiguity and reduce the number of iterations required in the grid search. The mean square error of the proposed CFO estimator considerably improves compared with the existing methods at a low SNR regime. The bit error rate performance of the proposed method for MIMO SC-FDMA is also compared with the existing methods. Finally, the proposed estimator is validated by implementing it on radio frequency hardware testbed over an indoor propagation environment. [ABSTRACT FROM AUTHOR]

This paper presents a low-complexity linear minimum mean-squared-error-based iterative soft interference cancellation (LMMSE-ISIC) scheme for multiple-input multiple-output (MIMO) systems. To avoid direct matrix inversion in the conventional LMMSE-ISIC scheme, expressions are derived to obtain the filtering vector and estimate of the conventional scheme; the direct inverse operation is then replaced by the vector-based operations in the proposed scheme. Thus, for the worst-case scenario in the proposed scheme, the complexity of estimating each transmit symbol is approximately proportional to the square of the number of receive antennas; in the conventional scheme, it is approximately proportional to the cube of the number of receive antennas. In addition, because there are no approximations for deriving the filtering vector and estimate, the proposed LMMSE-ISIC scheme achieves near-optimum performance identical to that of the conventional scheme, which is close to the matched filter bound of the channel. Further, the extension of the proposed scheme for iterative detection and decoding is developed for coded systems. The simulated results confirm that the conventional and proposed schemes outperform the approximated matrix inversion based schemes and achieve the identical error performance in both uncoded and coded systems, while the complexity order of the proposed scheme is similar to or even lower than those of the approximation schemes. Therefore, the proposed scheme can be considered an effective near-optimum LMMSE-based iterative detection approach for MIMO systems, especially for massive MIMO systems with a high load factor. [ABSTRACT FROM AUTHOR]

This paper considers stochastic convex optimization problems with two sets of constraints: (a) deterministic constraints on the domain of the optimization variable, which are difficult to project onto; and (b) deterministic or stochastic constraints that admit efficient projection. Problems of this form arise frequently in the context of semidefinite programming as well as when various NP-hard problems are solved approximately via semidefinite relaxation. Since projection onto the first set of constraints is difficult, it becomes necessary to explore projection-free algorithms, such as the stochastic Frank-Wolfe (FW) algorithm. On the other hand, the second set of constraints cannot be handled in the same way, and must be incorporated as an indicator function within the objective function, thereby complicating the application of FW methods. Similar problems have been studied before; however, they suffer from slow convergence rates. This work, equipped with momentum based gradient tracking technique, guarantees fast convergence rates on par with the best-known rates for problems without the second set of constraints. Zeroth-order variants of the proposed algorithms are also developed and again improve upon the state-of-the-art rate results. We further propose the novel trimmed FW variants that enjoy the same convergence rates as their classical counterparts, but are empirically shown to require significantly fewer calls to the linear minimization oracle speeding up the overall algorithm. The efficacy of the proposed algorithms is tested on relevant applications of sparse matrix estimation, clustering via semidefinite relaxation, and uniform sparsest cut problem. [ABSTRACT FROM AUTHOR]

We consider the problem of recovering random graph signals from nonlinear measurements. For this setting, closed-form Bayesian estimators are usually intractable and even numerical evaluation may be difficult to compute for large networks. In this paper, we propose a graph signal processing (GSP) framework for random graph signal recovery that utilizes information on the structure behind the data. First, we develop the GSP-linear minimum mean-squared-error (GSP-LMMSE) estimator, which minimizes the mean-squared-error (MSE) among estimators that are represented as an output of a graph filter. The GSP-LMMSE estimator is based on diagonal covariance matrices in the graph frequency domain, and thus, has reduced complexity compared with the LMMSE estimator. This property is especially important when using the sample-mean estimators that are based on a training dataset. We then state conditions under which the low-complexity GSP-LMMSE estimator coincides with the optimal LMMSE estimator. Next, we develop an approximate parametrization of the GSP-LMMSE estimator by graph filters. We present three implementations of the parametric GSP-LMMSE estimator for typical graph filters. These parametric graph filters are more robust to outliers and to network topology changes. In our simulations, we evaluate the performance of the proposed GSP-LMMSE estimators for the problem of state estimation in power systems, which can be interpreted as a graph signal recovery task. We show that the proposed sample-GSP estimators outperform the sample-LMMSE estimator for a limited training dataset and that the parametric GSP-LMMSE estimators are more robust to topology changes in the form of adding/removing vertices/edges. [ABSTRACT FROM AUTHOR]

Conventionalliterature on array signal processing (ASP) is based on the “capture first, process later” philosophy and to this end, signal processing algorithms are typically decoupled from the hardware. This poses fundamental limitations because if the sensors result in information loss, the algorithms may no longer be able to achieve their guaranteed performance. In this paper, our goal is to overcome the barrier of information loss via sensor saturation and clipping. This is a significant problem in application areas including physiological monitoring and extra-terrestrial exploration where the amplitudes may be unknown or larger than the dynamic range of the sensor. To overcome this fundamental bottleneck, we propose “computational arrays” which are based on a co-design approach so that a collaboration between the sensor array hardware and algorithms can be harnessed. Our work is inspired by the recently introduced unlimited sensing framework. In this context, our computational arrays encode the high-dynamic-range information by folding the signal amplitudes, thus introducing a new form of information loss in terms of the modulo measurements. On the decoding front, we develop mathematically guaranteed recovery algorithms for spatio-temporal array signal processing tasks that include DoA estimation, beamforming and signal reconstruction. Numerical examples corroborate the applicability of our approach and pave a path for the development of novel computational arrays for ASP. [ABSTRACT FROM AUTHOR]

Learning generative probabilistic models is a core problem in machine learning, which presents significant challenges due to the curse of dimensionality. This paper proposes a joint dimensionality reduction and non-parametric density estimation framework, using a novel estimator that can explicitly capture the underlying distribution of appropriate reduced-dimension representations of the input data. The idea is to jointly design a nonlinear dimensionality reducing auto-encoder to model the training data in terms of a parsimonious set of latent random variables, and learn a canonical low-rank tensor model of the joint distribution of the latent variables in the Fourier domain. The proposed latent density model is non-parametric and “universal,” as opposed to the predefined prior that is assumed in variational auto-encoders. Joint optimization of the auto-encoder and the latent density estimator is pursued via a formulation which learns both by minimizing a combination of the negative log-likelihood in the latent domain and the auto-encoder reconstruction loss. We demonstrate that the proposed model achieves very promising results on toy, tabular, and image datasets on regression tasks, sampling, and anomaly detection. [ABSTRACT FROM AUTHOR]

In this paper, two types of linear estimators are considered for three related estimation set-ups involving set-theoretic uncertainty pertaining to $ \mathcal {H}_{2}$ and $ \mathcal {H}_{\infty }$ balls of frequency-responses. The problems at stake correspond to robust $ \mathcal {H}_{2}$ and $ \mathcal {H}_{\infty }$ estimation in the face of non-parametric “channel-model’ uncertainty and to a nominal $ \mathcal {H}_{\infty }$ estimation problem. The estimators considered here are defined by minimizing the worst-case estimation error magnitude over the “uncertainty set” or by minimizing an average cost under the constraint that the worst-case error of any admissible estimator does not exceed a prescribed value. The main point is to explore the derivation of estimators which may be viewed as less conservative alternatives to minimax estimators, or in other words, that allow for trade-offs between worst-case performance and better performance over “large” subsets of the uncertainty set. The “average costs” over $ \mathcal {H}_{2}-$ signal balls are obtained as limits of averages over sets of finite impulse responses, as their length grows unbounded. The estimator design problems for the two types of estimators and the three set-ups addressed here are recast as semi-definite programming problems (SDPs, for short). These SDPs are solved in the case of simple examples to illustrate the potential of the “average cost/worst-case constraint” estimators to mitigate the inherent conservatism of the minimax estimators. [ABSTRACT FROM AUTHOR]

Direction of arrival (DOA) estimation in array processing using uniform/sparse linear arrays is concerned in this paper. While sparse methods via approximate parameter discretization have been popular in the past decade, the discretization may cause problems, e.g., modeling error and increased computations due to dense sampling. In this paper, an exact discretization-free method, named as sparse and parametric approach (SPA), is proposed for uniform and sparse linear arrays. SPA carries out parameter estimation in the continuous range based on well-established covariance fitting criteria and convex optimization. It guarantees to produce a sparse parameter estimate without discretization required by existing sparse methods. Theoretical analysis shows that the SPA parameter estimator is a large-snapshot realization of the maximum likelihood estimator and is statistically consistent (in the number of snapshots) under uncorrelated sources. Other merits of SPA include improved resolution, applicability to arbitrary number of snapshots, robustness to correlation of the sources and no requirement of user-parameters. Numerical simulations are carried out to verify our analysis and demonstrate advantages of SPA compared to existing methods. [ABSTRACT FROM PUBLISHER]

The problem of distributed inference with M-ary quantized data at the sensors is investigated in the presence of Byzantine attacks. We assume that the Byzantine nodes attack the inference network by modifying the symbol corresponding to the quantized data to one of the other symbols in the quantization alphabet-set and transmitting falsified symbol to the fusion center (FC). In this paper, we find the optimal Byzantine attack that blinds any distributed inference network. As the quantization alphabet size increases, a tremendous improvement in the security performance of the distributed inference network is observed. In addition to the perfect channel case, in Appendix A , we also analyze the optimal Byzantine attack when the channel between the nodes and the FC is noisy and is modelled as a discrete M-ary channel. We also investigate the optimal attack within the restricted space of highly-symmetric attack strategies, that maximally degrades the performance of the inference network in the presence of resource-constrained Byzantine attacks. A reputation-based scheme for identifying malicious nodes is also presented as the network's strategy to mitigate the impact of Byzantine threats on the inference performance of the distributed sensor network. We also provide asymptotic analysis to find the optimal reputation-based scheme as a function of the fraction of compromised nodes in the network. [ABSTRACT FROM PUBLISHER]

This paper addresses the problem of direction of arrival (DOA) estimation and blind source separation (BSS) for nonstationary signals in the underdetermined case. These two problems are strongly related to the mixing matrix estimation problem. To deal with the nonstationary characteristics of signals, this study uses high-resolution quadratic time-frequency distributions (TFDs) to reduce cross-terms while keeping a good resolution for the construction of spatial TFDs. The main contributions of this paper are two-fold. First, the formulation of a statistical test for the noise thresholding step improves robustness and avoids the use of empirical parameters; this test performs multisource selection of the time-frequency points where the signal of interest is present. Second, an algorithm based on image processing methods performs an auto-source selection for mixing matrix estimation. The results on simulated signals demonstrate an improvement of 10 dB in terms of normalized mean square error for BSS and 7% in terms of relative error for DOA estimation over standard methods. [ABSTRACT FROM PUBLISHER]

The product of the scanned responses of two colinear interleaved arrays estimates the spatial power spectral density of the received signal. This estimate has mathematical structure similar to the classical periodogram and hence is called the generalized periodogram. The generalized periodogram includes sparse arrays such as coprime arrays and two level nested arrays. This paper derives the analytical expected value of the generalized periodogram output for any signal. The generalized periodogram expected value is the convolution of the true spectral density with a windowing function. This result reveals that the expected product spectrum is not guaranteed to be positive definite. The paper also derives the covariance function of the generalized periodogram for a white Gaussian process. Finally, the paper illustrates the implicit equivalence between the product spectrum and estimation of the correlation function from the products of appropriately spaced pairs of sensors. [ABSTRACT FROM PUBLISHER]

The rapid proliferation of wireless devices offers new means to infer current conditions surrounding mobile users, their locations, and their actions. This is especially true for industrial, scientific, and medical (ISM) radio bands where communication protocols are often open and information about neighboring devices abounds. This paper examines efficient occupancy estimation based on Wi-Fi metadata, with an emphasis on algorithms attuned to directional antenna technology. Certain Wi-Fi interfaces can be switched into monitoring mode, an operating state where all local packets are observed and recorded. Using a network of such monitoring sensors, it is then possible to estimate the number of active devices within a specific area. The envisioned estimators take as input received signal strength indicators and media access control addresses. By using directional monitoring antennas, one can gain additional and/or more discriminating information about current conditions, thereby yielding enhanced occupancy estimates. This paper introduces novel estimation algorithms and characterizes the performance gains associated with RF-aware sensing devices. Experimental results based on a prototype implementation of this distributed monitoring system provide further supporting evidence for the proposed techniques. [ABSTRACT FROM PUBLISHER]