
We propose superresolution MIMO channel estimators for millimeterwave
(mmWave) systems that employ hybrid analog and digital beamforming and
generalized spatial modulation, respectively. Exploiting the inherent sparsity
of mmWave channels, the channel estimation problem is formulated as an atomic
norm minimization that enhances sparsity in the continuous angles of departure
and arrival. Both pilotassisted and dataaided channel estimators are
developed, with the former one formulated as a convex problem and the latter as
a nonconvex problem. To solve these formulated channel estimation problems, we
develop a computationally efficient conjugate gradient descent method based on
nonconvex factorization which restricts the search space to lowrank matrices.
Simulation results are presented to illustrate the superior channel estimation
performance of the proposed algorithms for both types of mmWave systems
compared to the existing compressedsensingbased estimators with finely
quantized angle grids.

Most existing approaches to coexisting communication/radar systems assume
that the radar and communication systems are coordinated, i.e., they share
information, such as relative position, transmitted waveforms and channel
state. In this paper, we consider an uncoordinated scenario where a
communication receiver is to operate in the presence of a number of radars, of
which only a subset may be active, which poses the problem of estimating the
active waveforms and the relevant parameters thereof, so as to cancel them
prior to demodulation. Two algorithms are proposed for such a joint waveform
estimation/data demodulation problem, both exploiting sparsity of a proper
representation of the interference and of the vector containing the errors of
the data block, so as to implement an iterative joint interference removal/data
demodulation process. The former algorithm is based on classical ongrid
compressed sensing (CS), while the latter forces an atomic norm (AN)
constraint: in both cases the radar parameters and the communication
demodulation errors can be estimated by solving a convex problem. We also
propose a way to improve the efficiency of the ANbased algorithm. The
performance of these algorithms are demonstrated through extensive simulations,
taking into account a variety of conditions concerning both the interferers and
the respective channel states.

We study the problem of estimating $\beta \in \mathbb{R}^p$ from its noisy
linear observations $y= X\beta+ w$, where $w \sim N(0, \sigma_w^2 I_{n\times
n})$, under the following highdimensional asymptotic regime: given a fixed
number $\delta$, $p \rightarrow \infty$, while $n/p \rightarrow \delta$. We
consider the popular class of $\ell_q$regularized least squares (LQLS)
estimators, a.k.a. bridge, given by the optimization problem: \begin{equation*}
\hat{\beta} (\lambda, q ) \in \arg\min_\beta \frac{1}{2} \yX\beta\_2^2+
\lambda \\beta\_q^q, \end{equation*} and characterize the almost sure limit
of $\frac{1}{p} \\hat{\beta} (\lambda, q ) \beta\_2^2$. The expression we
derive for this limit does not have explicit forms and hence are not useful in
comparing different algorithms, or providing information in evaluating the
effect of $\delta$ or sparsity level of $\beta$. To simplify the expressions,
researchers have considered the ideal "nonoise" regime and have characterized
the values of $\delta$ for which the almost sure limit is zero. This is known
as the phase transition analysis.
In this paper, we first perform the phase transition analysis of LQLS. Our
results reveal some of the limitations and misleading features of the phase
transition analysis. To overcome these limitations, we propose the study of
these algorithms under the low noise regime. Our new analysis framework not
only sheds light on the results of the phase transition analysis, but also
makes an accurate comparison of different regularizers possible.

The focus of this paper is on coexistence between a communication system and
a pulsed radar sharing the same bandwidth. Based on the fact that the
interference generated by the radar onto the communication receiver is
intermittent and depends on the density of scattering objects (such as, e.g.,
targets), we first show that the communication system is equivalent to a set of
independent parallel channels, whereby precoding on each channel can be
introduced as a new degree of freedom. We introduce a new figure of merit,
named the {\em compound rate}, which is a convex combination of rates with and
without interference, to be optimized under constraints concerning the
signaltointerferenceplusnoise ratio (including {\em signaldependent}
interference due to clutter) experienced by the radar and obviously the powers
emitted by the two systems: the degrees of freedom are the radar waveform and
the aforementioned encoding matrix for the communication symbols. We provide
closedform solutions for the optimum transmit policies for both systems under
two basic models for the scattering produced by the radar onto the
communication receiver, and account for possible correlation of the
signalindependent fraction of the interference impinging on the radar. We also
discuss the region of the achievable communication rates with and without
interference. A thorough performance assessment shows the potentials and the
limitations of the proposed coexisting architecture.

The millimeterwave (mmWave) fulldimensional (FD) MIMO system employs planar
arrays at both the base station and user equipment and can simultaneously
support both azimuth and elevation beamforming. In this paper, we propose
atomicnormbased methods for mmWave FDMIMO channel estimation under both
uniform planar arrays (UPA) and nonuniform planar arrays (NUPA). Unlike
existing algorithms such as compressive sensing (CS) or subspace methods, the
atomicnormbased algorithms do not require to discretize the angle spaces of
the angle of arrival (AoA) and angle of departure (AoD) into grids, thus
provide much better accuracy in estimation. In the UPA case, to reduce the
computational complexity, the original largescale 4D atomic norm minimization
problem is approximately reformulated as a semidefinite program (SDP)
containing two decoupled twolevel Toeplitz matrices. The SDP is then solved
via the alternating direction method of multipliers (ADMM) where each iteration
involves only closedform computations. In the NUPA case, the atomicnormbased
formulation for channel estimation becomes nonconvex and a
gradientdecentbased algorithm is proposed to solve the problem. Simulation
results show that the proposed algorithms achieve better performance than the
CSbased and subspacebased algorithms.

Feature aided tracking can often yield improved tracking performance over the
standard multiple target tracking (MTT) algorithms with only kinematic
measurements. However, in many applications, the feature signal of the targets
consists of sparse Fourierdomain signals. It changes quickly and nonlinearly
in the time domain, and the feature measurements are corrupted by missed
detections and misassociations. These two factors make it hard to extract the
feature information to be used in MTT. In this paper, we develop a
featureaided nearest neighbour joint probabilistic data association filter
(NNJPDAF) for joint MTT and feature extraction in dense target environments.
To estimate the rapidly varying feature signal from incomplete and corrupted
measurements, we use the atomic norm constraint to formulate the sparsity of
feature signal and use the $\ell_1$norm to formulate the sparsity of the
corruption induced by misassociations. Based on the sparse representation, the
feature signal are estimated by solving a semidefinite program (SDP) which is
convex. We also provide an iterative method for solving this SDP via the
alternating direction method of multipliers (ADMM) where each iteration
involves closedform computation. With the estimated feature signal,
refiltering is performed to estimate the kinematic states of the targets,
where the association makes use of both kinematic and feature information.
Simulation results are presented to illustrate the performance of the proposed
algorithm in a radar application.

We consider the problem of superresolving the line spectrum of a
multisinusoidal signal from a finite number of samples, some of which may be
completely corrupted. Measurements of this form can be modeled as an additive
mixture of a sinusoidal and a sparse component. We propose to demix the two
components and superresolve the spectrum of the multisinusoidal signal by
solving a convex program. Our main theoretical result is that up to
logarithmic factors this approach is guaranteed to be successful with high
probability for a number of spectral lines that is linear in the number of
measurements, even if a constant fraction of the data are outliers. The result
holds under the assumption that the phases of the sinusoidal and sparse
components are random and the line spectrum satisfies a minimumseparation
condition. We show that the method can be implemented via semidefinite
programming and explain how to adapt it in the presence of dense perturbations,
as well as exploring its connection to atomicnorm denoising. In addition, we
propose a fast greedy demixing method which provides good empirical results
when coupled with a local nonconvexoptimization step.

In this paper, we consider the problem of joint delayDoppler estimation of
moving targets in a passive radar that makes use of orthogonal
frequencydivision multiplexing (OFDM) communication signals. A compressed
sensing algorithm is proposed to achieve supperresolution and better accuracy,
using both the atomic norm and the $\ell_1$norm. The atomic norm is used to
manifest the signal sparsity in the continuous domain. Unlike previous works
which assume the demodulation to be error free, we explicitly introduce the
demodulation error signal whose sparsity is imposed by the $\ell_1$norm. On
this basis, the delays and Doppler frequencies are estimated by solving a
semidefinite program (SDP) which is convex. We also develop an iterative method
for solving this SDP via the alternating direction method of multipliers (ADMM)
where each iteration involves closedform computation. Simulation results are
presented to illustrate the high performance of the proposed algorithm.

In many application areas we are faced with the following question: Can we
recover a sparse vector $x_o \in \mathbb{R}^N$ from its undersampled set of
noisy observations $y \in \mathbb{R}^n$, $y=A x_o+w$. The last decade has
witnessed a surge of algorithms and theoretical results addressing this
question. One of the most popular algorithms is the $\ell_p$regularized least
squares (LPLS) given by the following formulation: \[ \hat{x}(\gamma,p )\in
\arg\min_x \frac{1}{2}\y  Ax\_2^2+\gamma\x\_p^p, \] where $p \in [0,1]$.
Despite the nonconvexity of these problems for $p<1$, they are still appealing
because of the following folklores in compressed sensing: (i) $\hat{x}(\gamma,p
)$ is closer to $x_o$ than $\hat{x}(\gamma,1)$. (ii) If we employ iterative
methods that aim to converge to a local minima of LPLS, then under good
initialization these algorithms converge to a solution that is closer to $x_o$
than $\hat{x}(\gamma,1)$. In spite of the existence of plenty of empirical
results that support these folklore theorems, the theoretical progress to
establish them has been very limited.
This paper aims to study the above folklore theorems and establish their
scope of validity. Starting with approximate message passing algorithm as a
heuristic method for solving LPLS, we study the impact of initialization on the
performance of AMP. Then, we employ the replica analysis to show the connection
between the solution of AMP and $\hat{x}(\gamma, p)$ in the asymptotic
settings. This enables us to compare the accuracy of $\hat{x}(\gamma,p)$ for $p
\in [0,1]$. In particular, we will characterize the phase transition and noise
sensitivity of LPLS for every $0\leq p\leq 1$ accurately. Our results in the
noiseless setting confirm that LPLS exhibits the same phase transition for
every $0\leq p <1$ and this phase transition is much higher than that of LASSO.