IEEE COMMUNICATIONS LETTERS, VOL. 21, NO. 8, AUGUST 2017 1783
DOA Estimation for Noncircular Sources with
Multiple Noncoherent Subarrays
Fei Wen, Member, IEEE, Wei Xie, Xin Chen, Member, IEEE, and Peilin Liu, Member, IEEE
Abstract—This letter addresses the direction-of-arrival
estimation issue for noncircular sources with multiple noncoher-
ent subarrays. First, we present the Cramer–Rao lower bound
for this problem. Then, for multiple noncoherent subarrays,
we propose the asymptotically minimum variance second-order
estimators for the two cases of circular and noncircular sources,
respectively. Furthermore, we propose a computationally efficient
MUSIC-like algorithm for such arrays which can exploit the
noncircularity of the sources. Finally, numerical examples have
been provided to demonstrate the performance of the new
methods.
Index Terms— Noncircular signals, direction of arrival (DOA)
estimation, MUSIC algorithm, time varying arrays, distributed
sensor array networks.
I. INTRODUCTION
D
IRECTION-OF-ARRIVAL (DOA) estimation of multiple
sourcesusinganarrayisanessentialtaskinmanyappli-
cations, such as radar, sonar, communications, geophysics,
tracking and localization [1]. In this work, we consider
the DOA estimation problem using multiple noncoherent
subarrays. Such arrays arise in the many important applica-
tions, such as time-varying arrays [5]–[11], subarray sam-
pling [12], [13] and partially coherent arrays. It is also the
case of distributed sensor array networks, in which each
array is perfectly coherent locally but the correlation between
different arrays is irrelevant [2]–[4].
All of the above works [2]–[13] focused on either determin-
istic or stochastic circular signals, while this work considers
the case of noncircular signals. Noncircular signals are usually
encountered in the context of communications, such as AM,
MASK, BPSK or UQPSK signals. It has been shown that,
the noncircularity of complex signals can be exploited to
improve the performance of DOA estimation [14]–[19], [27],
beamforming [20]–[22] and time delay estimation [23], [24].
Existing DOA estimation methods for noncircular sig-
nals (e.g., [14]–[19] and the references therein) commonly
assume a single coherent array, and none of them applies
to multiple noncoherent subarrays. The goal of this work
is to develop DOA estimation methods for noncircular sig-
nals applicable to generalized multiple noncoherent subarrays.
Manuscript received March 10, 2017; revised April 12, 2017; accepted
May 7, 2017. Date of publication May 9, 2017; date of current version
August 10, 2017. This work was supported in part by the National Natural
Science Foundation of China (NSFC) under grants 61401501 and 61472442.
The associate editor coordinating the review of this paper and approving it
for publication was J. Prieto. (Corresponding author: Fei Wen.)
F. Wen, X. Chen, and P. Liu are with the Department of Electronic
Engineering, Shanghai Jiao Tong University, Shanghai 200240, China (e-mail:
wenfei@sjtu.edu.cn; xin.chen@sjtu.edu.cn; liupeilin@sjtu.edu.cn).
W. Xie is with the Department of Electronic Engineering, University
of Electronic Science and Technology of China, Chengdu 610054, China
(e-mail: raysheik@foxmail.com).
Digital Object Identifier 10.1109/LCOMM.2017.2703649
The main contributions of this work are as follows. First,
we present the CRB of DOA estimation for noncircular
complex Gaussian (NCG) sources with multiple noncoher-
ent subarrays. Then, the GLS estimator [10] is improved
to be an asymptotically minimum variance (AMV) second-
order (SO) estimator for circular sources, and the AMV
SO estimator for noncircular sources [15] is extended for
multiple noncoherent subarrays. For multiple sources, these
methods require multi-dimensional nonlinear optimization,
which prompts us to design a computationally simple MUSIC-
like method subsequently. Specifically, we propose a weighted
MUSIC algorithm for noncircular sources, which fuses the
spatial spectrum of the subarrays in a weighted summation
manner.
Notations: (·)
T
and (·)
H
denote the transpose and Hermitian
transpose, respectively. δ
l,k
is the Kronneck function. || · ||
2
F
,
(·) and (·) stand for the Frobenius norm, real and imaginary
part operators, respectively. vec(·) is the “vectorization” oper-
ator stacking the columns of the matrix one below another, and
v(·) denotes the operator obtained from vec(·) by eliminating
all supradiagonal elements of the matrix. ⊗ denotes the
Kronecker product. I
M
stands for an M × M identity matrix.
II. P
ROBLEM STATEMENT,MLESTIMATOR AND CRB
A. Signal Model
Consider an array consisting of L noncoherent subar-
rays receiving d narrowband far-field signals impinging with
unknown DOAs θ
1
, ··· ,θ
d
.Thel-th subarray has M
l
sensors.
The output of the l-th subarray at time t
l
at some local
reference point can be expressed as
y
l
(t
l
) = A
l
()x(t
l
) + n
l
(t
l
), t
l
= 1, ··· , N
l
(1)
where A
l
() =[a
l
(θ
1
) ··· a
l
(θ
d
)] is the full column rank
steering matrix with a
l
(θ) represents the complex manifold of
the l-th subarray and =[θ
1
··· θ
d
]
T
denotes the vector
of the unknown DOAs; x(t
l
) =[x
1
(t
l
) ··· x
d
(t
l
)]
T
and
n
l
(t
l
) are the complex signal amplitudes and additive white
sensor noise, respectively; x(t
l
) and n
l
(t
l
) are multivariate
independent zero-mean wide-sense stationary processes; n
l
(t
l
)
is assumed spatially white circular complex Gaussian (CG)
with E{n
l
(t
l
)n
H
l
(t
l
)}=σ
2
I
M
l
, whereas x(t
l
) is NCG with
P = E{x(t
l
)x
H
(t
l
)} and Q = E{x(t
l
)x
T
(t
l
)}.
From (1), the two covariance matrices corresponding to the
l-th subarray are
R
l
() = A
l
PA
H
l
+ σ
2
I
M
l
and
l
() = A
l
QA
T
l
(2)
where =[
T
T
σ
2
]
T
denotes the unknown para-
meter vector with being a real vector made from the
free real parameters P(i, i), {(P(i, j)), (P(i, j))}
j>i
and
{(Q(i, j)), (Q(i, j))}
j≥i
.
1558-2558 © 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission.
See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.