NavigationEvents

IEEE Transactions on Signal ProcessingResource Optimization of NonAdditive Utility Functions in Localized SCFDMA SystemsIn this paper, we study the problem of resource allocation in SCFDMA systems. A sumutility maximization is considered where the utility of each user may be neither additive nor super additive. Unlike OFDMA, in addition to the restriction of allocating a subchannel to at most one user, the multiple subchannels allocated to a user in SCFDMA must be consecutive. This makes the resource allocation problem prohibitively difficult and challenging. We provide a fundamental complexity analysis of the optimization problem for general nonadditive utility functions and show that not only the problem is NPhard but also approximating it within a factor better than $ {{2011}over {2012}}$ is not possible unless ${rm P}={rm NP}$. An efficient cutting plane algorithm is presented and five suboptimal heuristics are also presented that achieve near optimal solution in different scenarios. Computational results of the cutting plane and heuristic algorithms are reported and a comparison between these heuristics is provided as well.
Categories: Scientifical publications
Rigid Body Localization Using Sensor NetworksA framework for joint position and orientation estimation of a rigid body using range measurements is proposed. We consider a setup in which a few sensors are mounted on a rigid body. The absolute position of the rigid body is not known. However, we know how the sensors are mounted on the rigid body, i.e., the sensor topology is known. The rigid body is localized using noisy range measurements between the sensors and a few anchors (nodes with known absolute positions), and without using any inertial measurements. We propose a leastsquares (LS), and a number of constrained LS estimators, where the constrained estimators solve an optimization problem on the Stiefel manifold. As a benchmark, we derive a unitarily constrained Cramér–Rao bound. Finally, the known topology of the sensors can be perturbed during fabrication or if the body is not entirely rigid. To take these perturbations into account, constrained totalleastsquares estimators are also proposed.
Categories: Scientifical publications
Multiuser MISO Beamforming for Simultaneous Wireless Information and Power TransferThis paper studies a multiuser multipleinput singleoutput (MISO) broadcast simultaneous wireless information and power transfer (SWIPT) system, where a multiantenna access point (AP) sends wireless information and energy simultaneously via spatial multiplexing to multiple singleantenna receivers each of which implements information decoding (ID) or energy harvesting (EH). We aim to maximize the weighted sumpower transferred to all EH receivers subject to given minimum signaltointerferenceandnoise ratio (SINR) constraints at different ID receivers. In particular, we consider two types of ID receivers (referred to as Type I and Type II, respectively) without or with the capability of cancelling the interference from (a priori known) energy signals. For each type of ID receivers, we formulate the joint information and energy transmit beamforming design as a nonconvex quadratically constrained quadratic program (QCQP). First, we obtain the globally optimal solutions for our formulated QCQPs by applying an optimization technique socalled semidefinite relaxation (SDR). It is shown via SDR that under the condition of independently distributed user channels, no dedicated energy beam is used for the case of Type I ID receivers to achieve the optimal solution; while for the case of Type II ID receivers, employing no more than one energy beam is optimal. Next, we establish a new form of the celebrated uplink–downlink duality for our studied downlink beamforming problems, and thereby develop alternative algorithms to obtain the same optimal solutions as by SDR. Finally, numerical results are provided to evaluate the performance of proposed optimal beamforming designs for MISO SWIPT systems.
Categories: Scientifical publications
Theory, Design and Application of Arbitrary Order Arbitrary Delay FilterbanksClassical design of filterbanks results in equal order inverses, thus resulting in filterbanks with equal analysis and synthesis filter length. While it is known that filterbanks may have unequal orders, there does not exist a simple, systematic method to design a $p$thorder analysis and $q$thorder synthesis filterbank where $p ne q$. Also, the design criteria of the $p$ thorder analysis having $q$ thorder synthesis filters ($p ne q$) with a flexibility to control the system delay has never been addressed concomitantly. In this work, we propose the use of unit nondiagonal matrix polynomials to design such filterbanks for given orders $p$, $q$ with arbitrary delay. We demonstrate how these orders depend on the sequence of unit matrices. Not only the achievable ranges of $p$ , $q$ are found, but construction methods for all such cases are also suggested. The proposed design has several desirable properties such as structural perfect reconstruction in its factorized form, completeness, and ability to control the resulting system delay of the filterbank. Several design examples and an application example for audio illustrate the flexibility and usefulness of the proposed design.
Categories: Scientifical publications
<formula formulatype="inline"><tex Notation="TeX">$R$</tex> </formula>Dimensional ESPRITType Algorithms for Strictly SecondOrder NonCircular Sources and Their Performance AnalysisHighresolution parameter estimation algorithms designed to exploit the prior knowledge about incident signals from strictly secondorder (SO) noncircular (NC) sources allow for a lower estimation error and can resolve twice as many sources. In this paper, we derive the $R$ ${mathchar"702D}{rm D}$ NC Standard ESPRIT and the $R$ ${mathchar"702D}{rm D}$ NC Unitary ESPRIT algorithms that provide a significantly better performance compared to their original versions for arbitrary source signals. They are applicable to shiftinvariant $R$ ${mathchar"702D}{rm D}$ antenna arrays and do not require a centrosymmetric array structure. Moreover, we present a firstorder asymptotic performance analysis of the proposed algorithms, which is based on the error in the signal subspace estimate arising from the noise perturbation. The derived expressions for the resulting parameter estimation error are explicit in the noise realizations and asymptotic in the effective signaltonoise ratio (SNR), i.e., the results become exact for either high SNRs or a large sample size. We also provide mean squared error (MSE) expressions, where only the assumptions of a zero mean and finite SO moments of the noise are required, but no assumptions about its statistics are necessary. As a main result, we analytically prove that the asymptotic performance of both $R$ ${mathchar"702D}{rm D}$ NC ESPRITtype algorithms is identical in the higheffective SNR regime. Finally, 
case study shows that no improvement from strictly noncircular sources can be achieved in the special case of a single source.
Categories: Scientifical publications
<formula formulatype="inline"><tex Notation="TeX">${cal H}_2$</tex></formula> Sampled—Data Filtering of Linear SystemsThis paper addresses the ${cal H}_2$ state estimation problem in the context of sampleddata systems, which is closely related to ${cal H}_2$ filtering design under limited bandwidth constraints. Our first approach considers a fixed datarate setting, for which a discretetime equivalent system is devised and, based on this concept, sufficient conditions for the ${cal H}_2$ filtering design are yielded. These results are then generalized to obtain a robust filter that is able to cope with uneven datarates. Academical examples illustrate the theoretical results and point out the main characteristics of the design techniques devised in this paper.
Categories: Scientifical publications
Convex Optimization Approaches for Blind Sensor Calibration Using SparsityWe investigate a compressive sensing framework in which the sensors introduce a distortion to the measurements in the form of unknown gains. We focus on blind calibration, using measures performed on multiple unknown (but sparse) signals and formulate the joint recovery of the gains and the sparse signals as a convex optimization problem. We divide this problem in 3 subproblems with different conditions on the gains, specifially i) gains with different amplitude and the same phase, ii) gains with the same amplitude and different phase and iii) gains with different amplitude and phase. In order to solve the first case, we propose an extension to the basis pursuit optimization which can estimate the unknown gains along with the unknown sparse signals. For the second case, we formulate a quadratic approach that eliminates the unknown phase shifts and retrieves the unknown sparse signals. An alternative form of this approach is also formulated to reduce complexity and memory requirements and provide scalability with respect to the number of input signals. Finally for the third case, we propose a formulation that combines the earlier two approaches to solve the problem. The performance of the proposed algorithms is investigated extensively through numerical simulations, which demonstrates that simultaneous signal recovery and calibration is possible with convex methods when sufficiently many (unknown, but sparse) calibrating signals are provided.
Categories: Scientifical publications
Superoscillations With Optimum Energy ConcentrationOscillations of a bandlimited signal at a rate faster than the bandlimit are called “superoscillations” and have applications e.g. in superresolution and superdirectivity. The synthesis of superoscillating signals is a numerically difficult problem. Minimum energy superoscillatory signals seem attractive for applications because (i) the minimumenergy solution is unique (ii) it has the smallest energy cost (iii) it may yield a signal of the smallest possible amplitude. On the negative side, superoscillating functions of minimumenergy depend heavily on cancellation and give rise to expressions that have very large coefficients. Furthermore, these coefficients have to be found by solving equations that are very illconditioned. Surprisingly, we show that by dropping the minimum energy requirement practicality can be gained rather than lost. We give a method of constructing superoscillating signals that leads to coefficients and condition numbers that are smaller by several orders of magnitude than the minimumenergy solution, yet yields energies close to the minimum. In contrast with the minimumenergy method, which builds superoscillations by linearly combining functions with an illconditioned Gram matrix, our method combines orthonormal functions, the Gram matrix of which is obviously the identity. Another feature of the method is that it yields the superoscillatory signal that maximises the energy concentration in a given set, which may or may not include the superoscillatory segment.
Categories: Scientifical publications
Alternating Projections and DouglasRachford for Sparse Affine FeasibilityThe problem of finding a vector with the fewest nonzero elements that satisfies an underdetermined system of linear equations is an NPcomplete problem that is typically solved numerically via convex heuristics or nicelybehaved nonconvex relaxations. In this work we consider elementary methods based on projections for solving a sparse feasibility problem without employing convex heuristics. It has been shown recently that, locally, the fundamental method of alternating projections must converge linearly to a solution to the sparse feasibility problem with an affine constraint. In this paper we apply different analytical tools that allow us to show global linear convergence of alternating projections under familiar constraint qualifications. These analytical tools can also be applied to other algorithms. This is demonstrated with the prominent DouglasRachford algorithm where we establish local linear convergence of this method applied to the sparse affine feasibility problem.
Categories: Scientifical publications
LowSamplingRate UltraWideband Channel Estimation Using EquivalentTime SamplingIn this paper, a lowsamplingrate scheme for ultrawideband channel estimation is proposed. The scheme exploits multiple observations generated by transmitting multiple pulses. In the proposed scheme, $P$ pulses are transmitted to produce channel impulse response estimates at a desired sampling rate, while the ADC samples at a rate that is $P$ times slower. To avoid loss of fidelity, the number of sampling periods (based on the desired rate) in the interpulse interval is restricted to be coprime with $P$ . This condition is affected when clock drift is present and the transmitted pulse locations change. To handle this case, and to achieve an overall good channel estimation performance, without using prior information, we derive an improved estimator based on the bounded data uncertainty (BDU) model. It is shown that this estimator is related to the Bayesian linear minimum mean squared error (LMMSE) estimator. Channel estimation performance of the proposed subsampling scheme combined with the new estimator is assessed in simulation. The results show that high reduction in sampling rate can be achieved. The proposed estimator outperforms the least squares estimator in almost all cases, while in the high SNR regime it also outperforms the LMMSE estimator. In addition to channel estimation, a synchronization method is also proposed that utilizes the same pulse sequence used for channel estimation.
Categories: Scientifical publications
Elliptic Localization: Performance Study and Optimum Receiver PlacementRange based localization often relies on TOA that creates a circular contour or TDOA that yields hyperbolic surface to determine the position of an emitting object. Another approach that has gained attention recently is elliptic localization. Elliptic localization uses an active transmitter to send out a signal and several receivers to acquire the reflected or relayed signal from an object to obtain its location. This paper performs a fundamental investigation on the performance of elliptic localization, for both the synchronous and asynchronous cases between the transmitter and the receivers. Their performance is also characterized with respect to the hyperbolic localization. The study is based on the Cramér–Rao lower bound (CRLB) for the object location under Gaussian measurement noise as well as Gaussian errors in the transmitter and receiver positions. When the receiver positions are controllable, we derive the optimum receiver placement for elliptic positioning to improve the localization accuracy. We conclude the study with simulation results to demonstrate the elliptic localization performance using the CRLBs and the Maximum Likelihood estimators.
Categories: Scientifical publications
An EmpiricalBayes Approach to Recovering Linearly Constrained NonNegative Sparse SignalsWe propose two novel approaches for the recovery of an (approximately) sparse signal from noisy linear measurements in the case that the signal is a priori known to be nonnegative and obey given linear equality constraints, such as a simplex signal. This problem arises in, e.g., hyperspectral imaging, portfolio optimization, density estimation, and certain cases of compressive imaging. Our first approach solves a linearly constrained nonnegative version of LASSO using the maxsum version of the generalized approximate message passing (GAMP) algorithm, where we consider both quadratic and absolute loss, and where we propose a novel approach to tuning the LASSO regularization parameter via the expectation maximization (EM) algorithm. Our second approach is based on the sum–product version of the GAMP algorithm, where we propose the use of a Bernoulli nonnegative Gaussianmixture signal prior and a Laplacian likelihood and propose an EMbased approach to learning the underlying statistical parameters. In both approaches, the linear equality constraints are enforced by augmenting GAMP's generalizedlinear observation model with noiseless pseudomeasurements. Extensive numerical experiments demonstrate the stateoftheart performance of our proposed approaches.
Categories: Scientifical publications
Construction of Doppler Resilient Complete Complementary Code in MIMO RadarOrthogonal waveform design has become an important research topic in recent years due to the performance advantages of MIMO radar incorporating orthogonal transmitted signals. Complete complementary code (CCC) has ideal ambiguity along the zero Doppler axis, which is desirable in MIMO radar. However, the main obstacle for the extensive usage of CCC is its sensitivity to nonzero Doppler shift. In this paper, we introduce a systematic method of constructing Doppler resilient CCC in MIMO radar, where the range sidelobes can vanish at modest Doppler shifts. In the construction of Doppler resilient waveforms, the generalized Prouhet–Thue–Morse (GPTM) sequence plays a key role. Numerical results are presented to demonstrate the effectiveness of the proposed method.
Categories: Scientifical publications
KnowledgeAided Parametric Adaptive Matched Filter With Automatic Combining for Covariance EstimationIn this paper, a knowledgeaided parametric adaptive matched filter (KAPAMF) is proposed that utilizing both observations (including the test and training signals) and a priori knowledge of the spatial covariance matrix. Unlike existing KAPAMF methods, the proposed KAPAMF is able to automatically adjust the combining weight of a priori covariance matrix, thus gaining enhanced robustness against uncertainty in the prior knowledge. Meanwhile, the proposed KAPAMF is significantly more efficient than its KA nonparametric counterparts when the amount of training signals is limited. One distinct feature of the proposed KAPAMF is the inclusion of both the test and training signals for automatic determination of the combining weights for the prior spatial covariance matrix and observations. Numerical results are presented to demonstrate the effectiveness of the proposed KAPAMF, especially in the limited training scenarios.
Categories: Scientifical publications
A Variational Bayes Framework for Sparse Adaptive EstimationRecently, a number of mostly $ell_1$ norm regularized leastsquarestype deterministic algorithms have been proposed to address the problem of sparse adaptive signal estimation and system identification. From a Bayesian perspective, this task is equivalent to maximum a posteriori probability estimation under a sparsity promoting heavytailed prior for the parameters of interest. Following a different approach, this paper develops a unifying framework of sparse variational Bayes algorithms that employ heavytailed priors in conjugate hierarchical form to facilitate posterior inference. The resulting fully automated variational schemes are first presented in a batch iterative form. Then, it is shown that by properly exploiting the structure of the batch estimation task, new sparse adaptive variational Bayes algorithms can be derived, which have the ability to impose and track sparsity during realtime processing in a timevarying environment. The most important feature of the proposed algorithms is that they completely eliminate the need for computationally costly parameter finetuning, a necessary ingredient of sparse adaptive deterministic algorithms. Extensive simulation results are provided to demonstrate the effectiveness of the new sparse adaptive variational Bayes algorithms against stateoftheart deterministic techniques for adaptive channel estimation. The results show that the proposed algorithms are numerically robust and exhibit in general superior estimation performance compared to their deterministic counterparts.
Categories: Scientifical publications
BlockWise QRDecomposition for the Layered and Hybrid Alamouti STBC MIMO Systems: Algorithms and Hardware ArchitecturesUnlike the channel matrix in the spatial division multiplexing (SDM) multipleinput multipleoutput (MIMO) communication system, the equivalent channel matrix in the layered Alamouti spacetime block coding (STBC) MIMO system comprised 2by2 Alamouti subblocks. One novel property, found by Sayed about the QRdecomposition (QRD) of this equivalent channel matrix is that the produced ${rm Q}$ and ${rm R}$matrices are also matrices with Alamouti subblocks. Taking advantage of this property, we propose a new blockwise complex Givens rotation (BCGR) based algorithm and a triangular systolic array (TSA) to compute the QRD of the equivalent channel matrix in an Alamouti block by block manner. Implementation results reveal that our new TSA can compute QRDs of 4by4 equivalent channel matrices faster than any architecture that has been developed for the SDM MIMO system. This property of fast QRD makes our TSA very attractive for the layered Alamouti STBC MIMO system combined with the orthogonal frequency division multiplexing. Our new BCGR based approach can also be applied to the hybrid Alamouti STBC MIMO system, which is also a system with equivalent channel matrix consisting of Alamouti subblocks.
Categories: Scientifical publications
IRCI Free Range Reconstruction for SAR Imaging With Arbitrary Length OFDM PulseOur previously proposed OFDM with sufficient cyclic prefix (CP) synthetic aperture radar (SAR) imaging algorithm is interrangecell interference (IRCI) free and achieves ideally zero range sidelobes for range reconstruction. In this OFDM SAR imaging algorithm, the minimum required CP length is almost equal to the number of range cells in a swath, while the number of subcarriers of an OFDM signal needs to be more than the CP length. This makes the length of a transmitted OFDM sequence at least almost twice of the number of range cells in a swath and for a wide swath imaging, the transmitted OFDM pulse length becomes long, which may cause problems in some radar applications. In this paper, we propose a CPbased OFDM SAR imaging with arbitrary pulse length, which has IRCI free range reconstruction and its pulse length is independent of a swath width. We then present a novel design method for our proposed arbitrary length OFDM pulses. Simulation results are presented to illustrate the performances of the OFDM pulse design and the arbitrary pulse length CPbased OFDM SAR imaging.
Categories: Scientifical publications
A Barycentric Coordinate Based Distributed Localization Algorithm for Sensor NetworksThis paper studies the problem of determining the sensor locations in a large sensor network using only relative distance (range) measurements. Based on a generalized barycentric coordinate representation, our work generalizes the DILOC algorithm to the localization problem under arbitrary deployments of sensor nodes and anchor nodes. First, a criterion and algorithm are developed to determine a generalized barycentric coordinate of a node with respect to its neighboring nodes, which do not require the node to be inside the convex hull of its neighbors. Next, for the localization problem based on the generalized barycentric coordinate representation, a necessary and sufficient condition for the localizability of a sensor network with a generic configuration is obtained. Finally, a new linear iterative algorithm is proposed to ensure distributed implementation as well as global convergence to the true coordinates.
Categories: Scientifical publications
MultiStage Robust Chinese Remainder TheoremIt is wellknown that the traditional Chinese remainder theorem (CRT) is not robust in the sense that a small error in a remainder may cause a large reconstruction error. A robust CRT was recently proposed for a special case when the greatest common divisor (gcd) of all the moduli is more than 1 and the remaining integers factorized by the gcd are coprime. It basically says that the reconstruction error is upper bounded by the remainder error level $tau$ if $tau$ is smaller than a quarter of the gcd of all the moduli. In this paper, we consider the robust reconstruction problem for a general set of moduli. We first present a necessary and sufficient condition on the remainder errors with a general set of moduli and also a corresponding robust reconstruction method. This can be thought of as a singlestage robust CRT. We then propose a twostage robust CRT by grouping the moduli into several groups as follows. First, the singlestage robust CRT is applied to each group. Then, with these robust reconstructions from all the groups, the singlestage robust CRT is applied again across the groups. This is easily generalized to multistage robust CRT. With this twostage robust CRT, the robust reconstruction holds even when the remainder error level $tau$ is above the quarter of the gcd of all the moduli, and an algorithm on how to group a set of moduli for a better reconstruction robustness is proposed in some special cases.
Categories: Scientifical publications
Hierarchical Interference Mitigation for Massive MIMO Cellular NetworksWe propose a hierarchical interference mitigation scheme for massive MIMO cellular networks. The MIMO precoder at each base station (BS) is partitioned into an inner precoder and an outer precoder. The inner precoder controls the intracell interference and is adaptive to local channel state information (CSI) at each BS (CSIT). The outer precoder controls the intercell interference and is adaptive to channel statistics. Such hierarchical precoding structure reduces the number of pilot symbols required for CSI estimation in massive MIMO downlink and is robust to the backhaul latency. We study joint optimization of the outer precoders, the user selection, and the power allocation to maximize a general concave utility which has no closedform expression. We first apply random matrix theory to obtain an approximated problem with closedform objective. Then using the hidden convexity of the problem, we propose an iterative algorithm to find the optimal solution for the approximated problem. We also obtain a low complexity algorithm with provable convergence. Simulations show that the proposed design has significant gain over various stateoftheart baselines.
Categories: Scientifical publications
