Noise resilient solution and its analysis for multistatic localization using differential arrival times

We can only extract the differential time delay (DTD) measurements between the direct and reflected paths in multistatic localization when there is no synchronization in time between the transmitter and receiver and among the receivers. This paper first addresses the problem of multistatic localization of a fixed object when the transmitter position is not available by using the DTD measurements. We propose a two-step optimization method for jointly estimating the object and transmitter positions. In the first step, we formulate a non-convex constrained weighted least squares (CWLS) problem by transforming the DTD measurement model and introducing nuisance variables. Such a non-convex CWLS problem is then relaxed to a tractable convex semidefinite programming (SDP) problem by applying semidefinite relaxation. In the second step, the error coming from relaxation and approximation in the SDP solution is reduced iteratively through solving a generalized trust region subproblem (GTRS) in each iteration. If the receivers are synchronized such that the time difference of arrival (TDOA) measurements can be acquired in addition to DTD, we formulate a different CWLS problem by utilizing both DTD and TDOA measurements, which is solved by convex relaxation as well. The relaxed SDP problem can achieve the optimal solution of the CWLS problem, and further refinement is no longer needed. We conduct the mean square error (MSE) analysis to validate that both proposed methods are able to achieve the Cramer-Rao Lower Bound (CRLB) performance under small Gaussian noise, which is also validated by simulations.

Semidefinite relaxation method for unified near-Field and far-Field localization by AOA

In this paper, we address the motion parameter estimation problem using time difference of arrival (TDOA) measurements, where a mobile source starts from an initial position with constant velocity. The problem is formulated as three non-convex constrained weighted least squares (CWLS) problems. Owing to their non-convex nature, local convergence may occur when solving them by an iterative algorithm, implying good estimates are not guaranteed. We propose to solve these CWLS problems by applying semidefinite relaxation (SDR) to generate three different semidefinite programming (SDP) problems, namely, the motion unconstrained semidefinite programming (MU-SDP), motion constrained semidefinite programming (MC-SDP), and motion parameter direct semidefinite programming (MPD-SDP). The MU-SDP, MC-SDP, and MPD-SDP methods are then extended to the scenario of unknown propagation speed (PS), which is common in underwater and underground localizations. Specifically, for this scenario, the two-step MU-SDP, two-step MC-SDP, and two-step MPD-SDP methods are designed to keep the relaxed SDP problems tight by introducing the penalty function. Simulation results show that MU-SDP performs worse than MC-SDP and MPD-SDP owing to the ignorance of the motion equation, and MC-SDP and MPD-SDP are able to reach the Cramér-Rao lower bound (CRLB) accuracy.

Reactive power planning problem is the key to secure and economic operation of power systems. Optimal management of existing reactive sources leads to loss minimization and economic performance of the system. Because of the nonlinear inter-relation between the physical parameters of the electric grid, this problem is a highly nonlinear and nonconvex constrained optimization problem. The application of semidefinite programming (SDP) to power system problems has recently gained considerable research attention. A recent SDP formulation uses a convex relaxation to the nonconvex optimal power flow problem under some technical conditions. This paper proposes a novel equivalent convex optimization formulation for the optimal reactive power dispatch (ORPD) problem and presents a new framework for finding the global optimum. Numerical results for the IEEE 30-bus and 118-bus test systems show that the proposed scheme obtains the optimum operation point and outperforms various state-of-the-art methods significantly.

Visible light positioning (VLP) has drawn great attention in the field of indoor positioning as light communication has been popularized in low-cost Internet-of-Things (IOT) devices. In this paper, we investigate the VLP problem using the received signal strength (RSS) and by only considering the line-of-slight (LOS) propagation. The RSS-based VLP problem is highly nonlinear, and its solutions may be trapped in local optima without a good initial guess. To circumvent this difficulty, we propose closed-form solutions of the VLP problem considering a known or unknown user orientation. By applying the weighted least squares (WLS) method, the closed-form solutions are divided into two stages. In the stage-one WLS solution, the nonlinear VLP problem is transformed into a pseudo-linear form by introducing some auxiliary variables, which are considered to be independent of each other. The estimates of the stage-one WLS solution are further refined in the stage-two WLS solution by exploiting the constrained relationships among these defined variables. The simulation results show that the stage-two WLS solution provides good estimates for the user position and orientation. The proposed stage-two WLS solution outperforms the existing methods especially at a high signal-to-noise ratio (SNR).

Efficient semidefinite solutions for TDOA-based source localization under unknown PS

Asymptotically efficient GNSS trilateration

This paper addresses the multistatic localization of a moving object in position and velocity using time delay (TD) and Doppler frequency shift (DFS) measurements, where the position of the transmitter is unknown and has not yet been synchronized with the receivers. Based on the TD and DFS measurements from the direct and indirect paths, we formulate a non-convex weighted least squares (WLS) minimization problem, and then apply semidefinite relaxation (SDR) to relax the WLS problem into a convex semidefinite program. Simulation results show that the proposed SDR method is able to achieve the Cramer-Rao lower bound accuracy under mild Gaussian noise condition and outperforms the existing method.

A vast majority of machine learning algorithms train their models and perform inference by solving optimization problems. In order to capture the learning and prediction problems accurately, structural constraints such as sparsity or low rank are frequently imposed or else the objective itself is designed to be a non-convex function. This is especially true of algorithms that operate in high-dimensional spaces or that train non-linear models such as tensor models and deep networks. The freedom to express the learning problem as a non-convex optimization problem gives immense modeling power to the algorithm designer, but often such problems are NP-hard to solve. A popular workaround to this has been to relax non-convex problems to convex ones and use traditional methods to solve the (convex) relaxed optimization problems. However this approach may be lossy and nevertheless presents significant challenges for large scale optimization. On the other hand, direct approaches to non-convex optimization have met with resounding success in several domains and remain the methods of choice for the practitioner, as they frequently outperform relaxation-based techniques - popular heuristics include projected gradient descent and alternating minimization. However, these are often poorly understood in terms of their convergence and other properties. This monograph presents a selection of recent advances that bridge a long-standing gap in our understanding of these heuristics. We hope that an insight into the inner workings of these methods will allow the reader to appreciate the unique marriage of task structure and generative models that allow these heuristic techniques to (provably) succeed. The monograph will lead the reader through several widely used nonconvex optimization techniques, as well as applications thereof. The goal of this monograph is to both, introduce the rich literature in this area, as well as equip the reader with the tools and techniques needed to analyze these simple procedures for non-convex problems.

In cognitive relay networks, the interferences from secondary users (SUs) and relays to primary users are constrained to be lower than a threshold. The interference constraints are difficult to satisfy when the channel state information (CSI) is imperfect. In this paper, we propose a robust beamforming scheme for the multi-antenna non-regenerative cognitive relay network where the multi-antenna relay with imperfect CSIs helps the communication of single-antenna SU. Our objective is to design a robust beamforming scheme which maximizes the system capacity subject to transmit power constraint and interference constraints. The bounded channel uncertainties are modeled using the worst-case model. The robust beamforming problem, neglecting the correlation of channel uncertainties, is reformulated as a convex semidefinite programming (SDP) by rank-one relaxation. This convex SDP is related with the worst-case relay transmit power minimization problem, which is further reformulated as a convex SDP, whose rank-one solution is proved to exist. Thus, we propose the suboptimal solution to the robust beamforming problem which is found effectively by solving two convex SDPs. Simulation results are provided to demonstrate the effectiveness of the proposed scheme.

Source localization in a wireless sensor network (WSN) is sensitive to the sensors' positions. In practice, due to mobility, the receivers' positions may be known inaccurately, leading to non-negligible degradation in source localization estimation performance. The goal of this paper is to develop a semidefinite programming (SDP) method using time-difference-of arrival (TDOA) and frequency-difference-of-arrival (FDOA) by taking the sensor position uncertainties into account. Specifically, we transform the commonly used maximum likelihood estimator (MLE) problem into a convex optimization problem to obtain an initial estimation. To reduce the coupling between position and velocity estimator, we also propose an iterative method to obtain the velocity and position, by using weighted least squares (WLS) method and SDP method, respectively. Simulations show that the method can approach the Cramér-Rao lower bound (CRLB) under both mild and high noise levels.

We revisit the problem of protein structure determination from geometrical restraints from NMR, using convex optimization. It is well-known that the NP-hard distance geometry problem of determining atomic positions from pairwise distance restraints can be relaxed into a convex semidefinite program (SDP). However, often the NOE distance restraints are too imprecise and sparse for accurate structure determination. Residual dipolar coupling (RDC) measurements provide additional geometric information on the angles between atom-pair directions and axes of the principal-axis-frame. The optimization problem involving RDC is highly non-convex and requires a good initialization even within the simulated annealing framework. In this paper, we model the protein backbone as an articulated structure composed of rigid units. Determining the rotation of each rigid unit gives the full protein structure. We propose solving the non-convex optimization problems using the sum-of-squares (SOS) hierarchy, a hierarchy of convex relaxations with increasing complexity and approximation power. Unlike classical global optimization approaches, SOS optimization returns a certificate of optimality if the global optimum is found. Based on the SOS method, we proposed two algorithms-RDC-SOS and RDC-NOE-SOS, that have polynomial time complexity in the number of amino-acid residues and run efficiently on a standard desktop. In many instances, the proposed methods exactly recover the solution to the original non-convex optimization problem. To the best of our knowledge this is the first time SOS relaxation is introduced to solve non-convex optimization problems in structural biology. We further introduce a statistical tool, the Cramér-Rao bound (CRB), to provide an information theoretic bound on the highest resolution one can hope to achieve when determining protein structure from noisy measurements using any unbiased estimator. Our simulation results show that when the RDC measurements are corrupted by Gaussian noise of realistic variance, both SOS based algorithms attain the CRB. We successfully apply our method in a divide-and-conquer fashion to determine the structure of ubiquitin from experimental NOE and RDC measurements obtained in two alignment media, achieving more accurate and faster reconstructions compared to the current state of the art.

Efficient closed-form solution for target localization using TDOA measurements in the presence of sensor position errors

It is known that the design of optimal transmit beamforming vectors for cognitive radio multicast transmission can be formulated as indefinite quadratic optimization programs. Given the challenges of such nonconvex problems, the conventional approach in literature is to recast them as convex semidefinite programs (SDPs) together with rank-one constraints. Then, these nonconvex and discontinuous constraints are dropped allowing for the realization of a pool of relaxed candidate solutions, from which various randomization techniques are utilized with the hope to recover the optimal solutions. However, it has been shown that such approach fails to deliver satisfactory outcomes in many practical settings, wherein the determined solutions are found to be unacceptably far from the actual optimality. On the contrary, we in this contribution tackle the aforementioned optimal beamforming problems differently by representing them as SDPs with additional reverse convex (but continuous) constraints. Nonsmooth optimization algorithms are then proposed to locate the optimal solutions of such design problems in an efficient manner. Our thorough numerical examples verify that the proposed algorithms offer almost global optimality whilst requiring relatively low computational load.

In this paper, we address the issue of optimal transmitter design to achieve physical layer security for a multiple-input single-output (MISO) cognitive radio network (CRN), in which a secondary user transmitter (SU-Tx) sends confidential information to a SU receiver (SU-Rx) on the same frequency band with a primary user (PU) in the presence of an eavesdropper receiver (ED-Rx). It is assumed that all the channel state information (CSI) of the secondary, primary and eavesdropper channels is not perfectly known at the SU-Tx. The optimal transmitter design, under the restriction of Gaussian signaling without preprocessing of information, involves a nonconvex semiinfinite optimization problem which maximizes the rate of the secondary link while avoiding harmful interference to the PU and keeping the eavesdropper totally ignorant of the messages sent regardless of the uncertainties in the CSI. We propose two approaches to solve this challenging optimization problem. The first one relates the original problem to a sequence of semiinfinite capacity-achieving transmitter design problems in an auxiliary CRN without any eavesdropper, which can then be solved through transformations and using convex semidefinite programs (SDPs). The second approach explores the hidden convexity of the problem and hence transforms it into a single SDP, which significantly reduces the computational complexity. Furthermore, a few heuristic beamforming solutions for the ease of implementation are also introduced. Finally, simulation results are presented to evaluate the performance of the proposed optimal and suboptimal solutions.