# Two-Dimensional DOA Estimation for Monostatic MIMO Radar with Electromagnetic Vector Received Sensors

Authors: Guimei Zheng; Jun Tang
Journal: International Journal of Antennas and Propagation
 (2016)

---

## Abstract

We study two-dimensional direction of arrival (2D-DOA) estimation problem of monostatic MIMO radar with the receiving array which consists of electromagnetic vector sensors (EMVSs). The proposed angle estimation algorithm can be applied to the arbitrary and unknown array configuration, which can be summarized as follows. Firstly, EMVSs in the receiver of a monostatic MIMO radar are used to measure all six electromagnetic-field components of an incident wavefield. The vector sensor array with the six unknown electromagnetic-field components is divided into six spatially identical subarrays. Secondly, ESPRIT is utilized to estimate the rotational invariant factors (RIFs). Parts of the RIFs are picked up to restore the source’s electromagnetic-field vector. Last, a vector cross product operation is performed between electric field and magnetic field to obtain the Pointing vector, which can offer the 2D-DOA estimation. Prior knowledge of array elements’ positions and angle searching procedure are not necessary for the proposed 2D-DOA estimation method. Simulation results prove the validity of the proposed method.


---

## Body

## 1. Introduction

MIMO radar has many merits in target detection and parameter estimation [1]. Therefore, it becomes a hot topic for researchers in the last decade [2–10]. The parameter estimation (such as direction of arrival (DOA) or polarization state angle (PSA)) is one of the main tasks in a radar system. Reference [2] separates the orthogonal signal received in receiving array by using the method of matched filtering to get a virtual array. It has more degrees of freedom of radar system than that of a conventional phased array radar with one version of signal waveform. So it can improve the DOA estimation performance. Reference [3] illuminates that the maximum number of localizable targets of MIMO radar is many times more than that of a counterpart phased array radar. Then, according to the characteristic of MIMO radar, many new parameter estimation algorithms have been developed. Firstly, classic Capon, APES, and CAPES algorithms are applied to monostatic MIMO radar to estimate DOA and amplitude of signals [4]. According to rotational invariance of transmitting and receiving array in bistatic MIMO radar [5], ESPRIT algorithm is used to estimate direction of departure (DOD) and DOA, but the algorithm needs additional pairing program between DODs and DOAs. So real-valued processing is introduced into bistatic MIMO radar to avoid this problem, and it also can improve estimation performance in the case of low snapshots [6]. In [7], by using the property of Kronecker product, reduced-dimensional MUSIC algorithm is presented, which employs one-dimensional angle searching to avoid the high computational cost compared with conventional two-dimensional angle searching algorithm. In addition, alternating projection based maximum likelihood estimation can provide more accurate estimations of target’s DOD and DOA [8], which outperforms the existing ESPRIT and MUSIC algorithms. Reference [9] presents a low-complexity algorithm with local one-dimensional searching for joint angle and mutual coupling estimation in bistatic MIMO radar. Reference [10] proposes the Jacobi-Anger expansion to transform the steering vector of uniform circular array into a steering vector with a Vandermonde structure, and then polynomial root finding algorithm is used to estimate DOA and DOD of the targets.Combined with the advantages of polarization diversity offered by vector antenna and waveform diversity offered by colocated MIMO radar, Jiang et al. [11] proposes polarimetric MIMO radar. It uses scalar antenna in transmitting array and cross-dipole vector antenna at receiving array. Then ESPRIT algorithm is used to estimate DOD, DOA, and PSA, but the method needs additional pairing program. ESPRIT-Root MUSIC algorithm [12] is extended to polarimetric MIMO radar in order to achieve automatic pairing. In addition, it has low complexity because of the absence of angle searching [13]. Reference [14] introduces parallel factor analysis to polarimetric MIMO radar, which can improve spatial resolution and achieve the goal of three-dimensional positioning. Note that the cross-dipole antennas are used in the above two papers. Therefore, the vector antennas only acquire two-component electric field of source. Reference [15] uses six-component electromagnetic vector sensor (EMVS) [16] to increase the information of a single vector antenna. It generally uses three orthogonal electric dipoles and magnetic loops to measure the three electric-field and three magnetic-field parameters of the incident signals, respectively. So there are more advantages in DOA estimation because of using more information. Unfortunately, the method in [15] is only applied to the array configuration with a single received electromagnetic vector sensor.In this paper, we solve the problem of two-dimensional (2D) DOA estimation of monostatic MIMO radar, whose received array consists of many EMVS antennas. We will show how to find the rotational invariance of spatial-polarimetric domain (SPD) in our signal model. Then a corresponding ESPRIT-like method is proposed to estimate 2D-DOA and PSA in monostatic MIMO radar. The proposed 2D-DOA and PSA estimation method does not need angle searching or any information about array configuration. Simulation results illuminate the advantages of the proposed method. Note that the proposed method can be seen as the extensional version of [17].Notation. Superscripts ·
∗, (
·
)
T, and (
·
)
H denote complex conjugation, transpose, and conjugate transpose, respectively. ⊙, ×, ⊗, and ⊕ denote Hadamard-Schur product, vector cross product, Kronecker product, and Khatri-Rao product (columnwise Kronecker product), respectively. R
e
{
·
} represents the real-value part of the entity inside. I
M
N is M
N
×
M
N identity matrix. T
r
{
·
} denotes the trace of the entity inside. ∠
z is the angle of complex z, which is in (
-
π
,
π
]. · denotes Frobenius norm. · denotes absolute value.


## 2. Signal Model

Consider a monostatic MIMO radar; the transmitting array hasM antennas. The receiving array is composed of N six-component EMVSs, as shown in Figure 1. Define the location of the mth sensor in the transmitting array as (
x
t
m
,
y
t
m
,
z
t
m
)
,
m
=
1
,
…
,
M and the location of the nth EMVS in the receiving array as (
x
r
n
,
y
r
n
,
z
r
n
)
,
n
=
1
,
…
,
N.Figure 1
Monostatic MIMO radar with arbitrarily spaced electromagnetic vector received sensors at unknown locations.TheM transmitting antennas transmit M orthogonal signals, which have the same bandwidth and center frequency; that is, s
(
τ
)
=
[
s
1
(
τ
)
,
…
,
s
M
(
τ
)
]
T
∈
C
M
×
1, where s
i
(
τ
) denotes ith transmitted signal; τ
=
1
,
…
,
Γ is the sampling index in fast time; therefore 1
/
Γ
∑
τ
=
1
Γ
s
(
τ
)
s
(
τ
)
H
=
I
M. Assume that there are K far-field point targets. The 2D-DOA of target are azimuth angle ϕ
∈
[
0,2
π
) and elevation angle θ
∈
[
0
,
π
], which are measured from positive x-axis and positive z-axis, respectively. PSA of target can be expressed as auxiliary polarization angle γ
∈
[
0
,
π
/
2
] and polarization phase difference η
∈
(
-
π
,
π
] [16], respectively. Then received signal can be represented as(1)
x
t
,
τ
=
A
r
θ
,
ϕ
⊕
A
p
o
l
θ
,
ϕ
,
γ
,
η
diag
⁡
b
t
A
t
T
θ
,
ϕ
s
τ
+
w
t
,
τ
∈
C
6
N
×
1
,where t denotes slow-time index, namely, snapshot index. A
r
(
θ
,
ϕ
)
=
[
a
r
(
θ
1
,
ϕ
1
)
,
…
,
a
r
(
θ
K
,
ϕ
K
)
]
∈
C
N
×
K denotes receiving spatially steering matrix. A
t
(
θ
,
ϕ
)
=
[
a
t
(
θ
1
,
ϕ
1
)
,
…
,
a
t
(
θ
K
,
ϕ
K
)
]
∈
C
M
×
K denotes transmitting spatially steering matrix. A
p
o
l
(
θ
,
ϕ
,
γ
,
η
)
=
[
a
p
o
l
(
θ
1
,
ϕ
1
,
γ
1
,
η
1
)
,
…
,
a
p
o
l
(
θ
K
,
ϕ
K
,
γ
K
,
η
K
)
] denotes SPD steering matrix of a single collocated EMVS. a
r
(
θ
k
,
ϕ
k
), a
t
(
θ
k
,
ϕ
k
), and a
p
o
l
(
θ
k
,
ϕ
k
,
γ
k
,
η
k
) are, respectively, given as follows: (2)
a
r
θ
k
,
ϕ
k
=
exp
⁡
-
j
2
π
λ
x
r
1
sin
⁡
θ
k
cos
⁡
ϕ
k
+
y
r
1
sin
⁡
θ
k
sin
⁡
ϕ
k
+
z
r
1
cos
⁡
θ
k
⁡
,
…
,
exp
⁡
-
j
2
π
λ
x
r
N
sin
⁡
θ
k
cos
⁡
ϕ
k
+
y
r
N
sin
⁡
θ
k
sin
⁡
ϕ
k
+
z
r
N
cos
⁡
θ
k
T
,
(3)
a
t
θ
k
,
ϕ
k
=
exp
⁡
-
j
2
π
λ
x
t
1
sin
⁡
θ
k
cos
⁡
ϕ
k
+
y
t
1
sin
⁡
θ
k
sin
⁡
ϕ
k
+
z
t
1
cos
⁡
θ
k
,
…
,
exp
⁡
-
j
2
π
λ
x
t
M
sin
⁡
θ
k
cos
⁡
ϕ
k
+
y
t
M
sin
⁡
θ
k
sin
⁡
ϕ
k
+
z
t
M
cos
⁡
θ
k
T
,
(4)
a
p
o
l
θ
k
,
ϕ
k
,
γ
k
,
η
k
=
e
h
=
a
p
o
l
,
1
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
2
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
3
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
4
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
5
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
6
θ
k
,
ϕ
k
,
γ
k
,
η
k
=
cos
⁡
ϕ
k
cos
⁡
θ
k
-
sin
⁡
ϕ
k
sin
⁡
ϕ
k
cos
⁡
θ
k
cos
⁡
ϕ
k
-
sin
⁡
θ
k
0
-
sin
⁡
ϕ
k
-
cos
⁡
ϕ
k
cos
⁡
θ
k
cos
⁡
ϕ
k
-
sin
⁡
ϕ
k
cos
⁡
θ
k
0
sin
⁡
θ
k
︸
Θ
θ
k
,
ϕ
k
sin
⁡
γ
k
e
j
η
k
cos
⁡
γ
k
︸
g
γ
k
,
η
k
,where e (the first three elements of a
p
o
l
(
θ
k
,
ϕ
k
,
γ
k
,
η
k
)) and h (the last three elements of a
p
o
l
(
θ
k
,
ϕ
k
,
γ
k
,
η
k
)) denote electric-field vector and magnetic-field vector, respectively. b
(
t
)
=
[
γ
1
(
t
)
,
…
,
γ
K
(
t
)
]
T
∈
C
K
×
1 denotes reflection coefficient vector, where γ
k
(
t
)
=
α
k
exp
⁡
(
j
2
π
f
k
t
). The amplitude of reflection coefficient α
k is determined by target’s radar cross section (RCS). The phase of reflection coefficient is determined by Doppler frequency of target. w
(
t
,
τ
)
∈
C
6
N
×
1 denotes additive white Gaussian noise with zero mean and covariance matrix σ
n
2
I
6
N.Multiplyx
(
t
,
τ
) by 1
/
Γ
s
m
(
τ
)
H
,
m
=
1
,
…
,
M; that is, perform matched filtering using the M transmitting signals; the output can be expressed as (5)
X
t
=
1
Γ
∑
τ
=
1
Γ
x
t
,
τ
s
τ
H
=
A
r
θ
,
ϕ
⊕
A
p
o
l
θ
,
ϕ
,
γ
,
η
diag
⁡
b
t
A
t
T
θ
,
ϕ
1
Γ
∑
τ
=
1
Γ
s
τ
s
τ
H
+
N
t
=
A
r
θ
,
ϕ
⊕
A
p
o
l
θ
,
ϕ
,
γ
,
η
diag
⁡
b
t
A
t
T
θ
,
ϕ
+
N
t
,where N
(
t
)
=
1
/
Γ
∑
τ
=
1
Γ
w
(
t
,
τ
)
s
(
τ
)
H. Vectorizing (5) yields the following 6
M
N-dimensional column vector: (6)
y
t
=
v
e
c
X
t
=
A
θ
,
ϕ
,
γ
,
η
b
t
+
n
t
∈
C
6
M
N
×
1
,where A
(
θ
,
ϕ
,
γ
,
η
)
=
A
t
(
θ
,
ϕ
)
⊕
A
r
(
θ
,
ϕ
)
⊕
A
p
o
l
(
θ
,
ϕ
,
γ
,
η
)
∈
C
6
M
N
×
K, whose columns are defined as virtual array steering vector a
(
θ
,
ϕ
,
γ
,
η
)
=
a
t
(
θ
,
ϕ
)
⊗
a
r
(
θ
,
ϕ
)
⊗
a
p
o
l
(
θ
,
ϕ
,
γ
,
η
). And the noise is transformed to n
(
t
)
=
v
e
c
[
N
(
t
)
]. The noise component corresponding to mth transmitting antenna and the n
t
h receiving antenna is n
n
,
m
=
w
n
s
m
H. Then the elements of corresponding noise covariance matrix are equal to (7)
E
n
n
2
,
m
2
H
·
n
n
1
,
m
1
=
E
w
n
2
s
m
2
H
H
·
w
n
1
s
m
1
H
=
E
s
m
2
w
n
2
H
·
w
n
1
s
m
1
H
=
0
else
σ
n
2
,
n
1
=
n
2
,
m
1
=
m
2
.It implies that n
(
t
) is still a complex Gaussian white noise with zero mean and covariance matrix σ
n
2
I
6
M
N.


## 3. Proposed Algorithm

For virtual array, after matched filtering, we extract the component with the same orient in EMVSs to consist a subarray. Then we can obtain six spatially identical subarrays. Here we take some EMVSs located iny-axis, for example, as shown in Figure 2.Figure 2
Six spatially identical subarrays offered by EMVS array.According to the SPD rotational invariance of the six subarrays, we can getC
6
2
=
15 rotational invariant factors (RIFs) by using ESPRIT algorithm. According to the knowledge of [17], we can calculate the incident electromagnetic wave’s electric field and magnetic field by only selecting five RIFs. Then we can perform vector cross product between electric field and magnetic field to obtain the Pointing vector of source. Last, 2D-DOA can be obtained by the triangle operation with the Pointing vector. The detailed derivation is given in the rest of paper.ForL snapshots, covariance matrix of data (6) is calculated by means of maximum likelihood estimation: R
^
=
1
/
L
∑
t
=
1
L
y
(
t
)
·
y
(
t
)
H
,
t
=
1
,
…
,
L. Then eigendecomposition is performed with the covariance matrix R
^ to obtain R
^
=
E
S
∑
S
E
S
H
+
σ
n
2
E
N
E
N
H, where E
S denotes 6
M
N
×
K signal subspace. Without noise, signal subspace and steering vectors {
a
(
θ
k
,
ϕ
k
,
γ
k
,
η
k
)
,
k
=
1
,
…
,
K
} have same space; that is, E
S
=
A
T. T denotes a nonsingular matrix. In the following, we will show how to estimate the 2D-DOA and PSA according to signal subspace E
S by means of ESPRIT-like algorithm.We firstly choose the array manifold ofi
t
h subarray by using selection matrix J
i. The following relationship equation with the array manifold will hold for kth target:(8)
J
i
a
θ
k
,
ϕ
k
,
γ
k
,
η
k
=
a
p
o
l
,
i
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
t
θ
k
,
ϕ
k
⊗
a
r
θ
k
,
ϕ
k
.Selection matrix can be expressed as J
i
=
[
I
M
N
⊗
e
i
]
T, where vector e
i
=
[
0
,
…
,
0
︸
i
-
1
,
1
,
0
,
…
,
0
︸
6
-
i
]
T
,
i
=
1
,
…
,
6. According to (8), we can obtain that the following rotational invariance relationship between i
t
h subarray and jth subarray holds: (9)
a
t
θ
k
,
ϕ
k
⊗
a
r
θ
k
,
ϕ
k
=
J
i
a
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
i
θ
k
,
ϕ
k
,
γ
k
,
η
k
=
J
j
a
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
j
θ
k
,
ϕ
k
,
γ
k
,
η
k
,where i
,
j
=
1
,
…
,
6. According to (9), we can obtain (10)
J
i
a
θ
k
,
ϕ
k
,
γ
k
,
η
k
=
J
j
a
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
i
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
j
θ
k
,
ϕ
k
,
γ
k
,
η
k
.Consider all of K targets; (10) can be transformed to matrix form:(11)
J
i
A
θ
,
ϕ
,
γ
,
η
=
J
j
A
θ
,
ϕ
,
γ
,
η
Φ
i
,
j
,where Φ
i
,
j
=
diag
⁡
[
a
p
o
l
,
i
θ
1
,
ϕ
1
,
γ
1
,
η
1
/
a
p
o
l
,
j
θ
1
,
ϕ
1
,
γ
1
,
η
1, …
,
a
p
o
l
,
i
θ
K
,
ϕ
K
,
γ
K
,
η
K
/
a
p
o
l
,
j
θ
K
,
ϕ
K
,
γ
K
,
η
K
]
T. For convenient presentation, we define RIF as χ
k
i
,
j
=
d
e
f
a
p
o
l
,
i
(
θ
k
,
ϕ
k
,
γ
k
,
η
k
)
/
a
p
o
l
,
j
(
θ
k
,
ϕ
k
,
γ
k
,
η
k
). Substituting E
S
=
A
T into (11), we can obtain rotational invariant equation about signal subspace E
S: (12)
J
i
E
S
=
J
j
E
S
Ψ
i
,
j
,where Ψ
i
,
j
=
(
T
i
,
j
)
-
1
Φ
i
,
j
T
i
,
j. In noise scenario, (12) can be solved by least squares or total least squares (TLS) method to obtain Ψ
i
,
j. Then, performing eigendecomposition with Ψ
i
,
j, RIF χ
k
i
,
j estimation can be obtained.Note that the RIFsχ
k
i
,
j
i
,
j
=
1
,
…
,
6 in different Φ
i
,
j have different orders. But there is no mismatch between χ
k
i
,
j and its corresponding eigenvector. Thus, different χ
k
i
,
j can be paired for the same source by pairing the orthogonal rows of T
i
,
j. Take, for example, χ
k
2,1 and χ
k
3,2. Let k denote the row index of the matrix element with the largest absolute value in the fth column of the matrix {
T
2,1
·
(
T
3,2
)
-
1
}. Then the kth row of T
2,1 must correspond to the fth row of T
3,2. Now, {
χ
k
2,1
,
k
=
1
,
…
,
K
} and χ
k
3,2
,
k
=
1
,
…
,
K have been correctly paired. All the RIFs χ
k
i
,
j
i
,
j
=
1
,
…
,
6 can be paired by performing the similar processing.There areC
6
2
=
15 choices for RIFs χ
k
i
,
j because of i
,
j
=
1
,
…
,
6. We set i
=
j
+
1; that is, 5 RIFs χ
k
j
+
1
,
j
,
j
=
1
,
…
,
5, are selected. In the following, we will show how to use the five RIFs to calculate the Pointing vector. According to the definition of χ
k
i
,
j, the following equations hold: (13)
a
p
o
l
,
2
θ
k
,
ϕ
k
,
γ
k
,
η
k
=
a
p
o
l
,
1
θ
k
,
ϕ
k
,
γ
k
,
η
k
·
χ
k
2,1
,
a
p
o
l
,
3
θ
k
,
ϕ
k
,
γ
k
,
η
k
=
a
p
o
l
,
1
θ
k
,
ϕ
k
,
γ
k
,
η
k
·
χ
k
2,1
·
χ
k
3,2
,
a
p
o
l
,
4
θ
k
,
ϕ
k
,
γ
k
,
η
k
=
a
p
o
l
,
1
θ
k
,
ϕ
k
,
γ
k
,
η
k
·
χ
k
2,1
·
χ
k
3,2
·
χ
k
4,3
,
a
p
o
l
,
5
θ
k
,
ϕ
k
,
γ
k
,
η
k
=
a
p
o
l
,
1
θ
k
,
ϕ
k
,
γ
k
,
η
k
·
χ
k
2,1
·
χ
k
3,2
·
χ
k
4,3
·
χ
k
5,4
,
a
p
o
l
,
6
θ
k
,
ϕ
k
,
γ
k
,
η
k
=
a
p
o
l
,
1
θ
k
,
ϕ
k
,
γ
k
,
η
k
·
χ
k
2,1
·
χ
k
3,2
·
χ
k
4,3
·
χ
k
5,4
·
χ
k
6,5
.According to Maxwell equation, we know that the result of cross-multiplication between the electric-field vector and the magnetic-field vector is equal to the Pointing vector as follows: (14)
a
p
o
l
,
1
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
2
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
3
θ
k
,
ϕ
k
,
γ
k
,
η
k
×
a
p
o
l
,
4
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
5
θ
k
,
ϕ
k
,
γ
k
,
η
k
a
p
o
l
,
6
θ
k
,
ϕ
k
,
γ
k
,
η
k
=
e
×
h
=
Γ
.Consider Γ
=
d
e
f
u
v
w is the Pointing vector, where u
=
sin
⁡
θ
cos
⁡
ϕ, v
=
sin
⁡
θ
sin
⁡
ϕ, and w
=
cos
⁡
θ are the three direction cosines aligning to x-axis, y-axis, and z-axis, respectively. Substituting (13) into (14), we have (15)
a
p
o
l
,
1
θ
k
,
ϕ
k
,
γ
k
,
η
k
2
1
χ
k
2,1
χ
k
2,1
·
χ
k
3,2
×
χ
k
2,1
·
χ
k
3,2
·
χ
k
4,3
χ
k
2,1
·
χ
k
3,2
·
χ
k
4,3
·
χ
k
5,4
χ
k
2,1
·
χ
k
3,2
·
χ
k
4,3
·
χ
k
5,4
·
χ
k
6,5
=
u
v
w
.We can compute the vector cross product operation of (15) according to the paired estimations χ
k
j
+
1
,
j
,
j
=
1
,
…
,
5. Then the estimations of direction cosines u
^
k, v
^
k, and w
^
k can be obtained by normalized processing with the result of vector cross product. Last, azimuth and elevation angle estimations can be obtained by the following triangle operations:(16)
ϕ
^
k
=
arctan
⁡
v
^
k
u
^
k
,
if
u
^
k
≥
0
,
ϕ
^
k
=
arctan
⁡
v
^
k
u
^
k
+
π
,
if
u
^
k
<
0
,
θ
^
k
=
arcsin
⁡
u
^
k
2
+
v
^
k
2
,
if
w
^
k
≥
0
,
θ
^
k
=
-
arcsin
⁡
u
^
k
2
+
v
^
k
2
+
π
,
if
w
^
k
<
0
,
k
=
1
,
…
,
K
.From the above processing, we can see that the information of antennas’ locations (
x
t
m
,
y
t
m
,
z
t
m
)
,
m
=
1
,
…
,
M, and (
x
r
n
,
y
r
n
,
z
r
n
)
,
n
=
1
,
…
,
N, is not utilized. It can be seen from (4) that(17)
g
γ
k
,
η
k
=
Θ
H
θ
k
,
ϕ
k
Θ
θ
k
,
ϕ
k
-
1
Θ
θ
k
,
ϕ
k
H
a
^
p
o
l
θ
k
,
ϕ
k
,
γ
k
,
η
k
.Substituting (
θ
^
k
,
ϕ
^
k
) into (17), we can obtain the estimations of g
(
γ
k
,
η
k
). Then we can calculate the corresponding PSA by the following formula:(18)
γ
^
k
=
a
r
c
t
a
n
g
^
γ
k
,
η
k
1
g
^
γ
k
,
η
k
2
η
^
k
=
∠
g
^
γ
k
,
η
k
1
-
∠
g
^
γ
k
,
η
k
2
.It can be seen from the above analysis that the proposed algorithm has the following advantages:(1)
Because we do not use the information of array configuration, it can be suitable for arbitrary array.
(2)
It requires no knowledge of the locations of antennas.
(3)
It has low complexity, since it requires no angle searching.
(4)
It has whole sphere validity region for 2D-DOA estimation.
(5)
It can be applied to EMVSs array in received side, but the method in [15] is only applied to a single EMVS in received side.Based on the above theoretical analysis, the implementation of the 2D-DOA estimation for monostatic MIMO radar with electromagnetic vector received sensors is summarized as follows.Step 1.
Perform matched filtering processing by using (5) and vectorization by using (6) according to the received data of (1).Step 2.
Compute the estimation of virtual array covariance matrixR
^ according to received data of (6). Compute the eigenvectors of R
^ to form the signal subspace E
S.Step 3.
Form the spatially rotational invariant equation of (12) according to the signal subspace E
S. Calculate the solution of (12) by using TLS, that is, estimation of {
Ψ
i
,
j
}, and then compute the eigenvalues of {
Ψ
i
,
j
} to obtain the RIFs estimations χ
k
j
+
1
,
j
,
j
=
1
,
…
,
5.Step 4.
Complete the pairing of the estimations{
χ
k
2,1
,
χ
k
3,2
,
χ
k
4,3
,
χ
k
5,4
,
χ
k
6,5
,
(
k
=
1
,
…
,
K
)
} for the same incident signal by using the relationship of eigenvalues and their corresponding eigenvectors.Step 5.
Perform the vector cross product operation of (15) according to the paired estimations {
χ
k
2,1
,
χ
k
3,2
,
χ
k
4,3
,
χ
k
5,4
,
χ
k
6,5
,
(
k
=
1
,
…
,
K
)
}. Then perform normalized processing to obtain estimations of direction cosines u
^
k, v
^
k, and w
^
k.Step 6.
Last, we can obtain azimuth and elevation angle estimations by (16). If necessary, we can calculate auxiliary polarization angle and polarization phase difference according to (18).


## 4. Complexity Analysis and Maximum Number of Localizable Targets

The computational complexity of the proposed algorithm is dominated by (1) formation of the covariance matrix of virtual MIMO array radar, (2) calculation of the estimation of array covariance matrix and its eigendecomposition to obtain signal subspace, (3) calculation of five RIFs, and (4) pairing scheme for the five RIFs. Process (1) requiresO
[
6
M
N
Γ
] flops; process (2) requires O
[
L
(
6
M
N
)
2
] flops; process (3) requires O
[
(
6
M
N
)
3
] flops; process (4) requires O
[
5
K
2
M
N
)
] flops; process (5) requires O
[
5
K
3
] flops. Therefore, the computational load of proposed algorithm is mainly O
[
6
M
N
Γ
+
L
(
6
M
N
)
2
+
(
6
M
N
)
3
+
5
K
2
M
N
)
+
5
K
3
] flops.The maximum number of localizable targets of the proposed algorithm is dominated by the choice of subarrays for ESPRIT algorithm. According to the choice of subarrays for ESPRIT in Figure2, we can obtain that the maximum number of targets that can be localized is M
N.


## 5. Computer Simulations

In this section, several numerical examples are presented to verify the effectiveness of the proposed algorithm. The root mean squared error (RMSE) of 2D-DOA and PSA estimations is defined asR
M
S
E
=
1
/
K
∑
k
=
1
K
E
(
ρ
^
k
-
ρ
k
)
2
,
ρ
=
θ
,
ϕ
,
γ
,
η. In following simulations, we consider a MIMO radar system with M
=
6 and N
=
6. We set the transmitting antennas located at x-axis with half-wavelength spacing and the receiving antennas located in y-axis with half-wavelength spacing. So the MIMO radar consists of L-shape array. Note that you can set MIMO radar with arbitrary configuration. There exists K
=
2 uncorrelated targets, which are located at the angles (
θ
1
,
ϕ
1
)
=
(
30
°
,
40
°
)
,
(
θ
2
,
ϕ
2
)
=
(
60
°
,
70
°
). The incident targets’ PSA are (
γ
1
,
η
1
)
=
(
45
°
,
90
°
)
,
(
γ
2
,
η
2
)
=
(
45
°
,
-
90
°
). The first target is left-circularly polarized, whereas the second target is right-circularly polarized. Note that the two targets have identical auxiliary polarization angles.Figure3 shows the histograms of estimation results with 500 Monte Carlo trials. In this simulation, assume that S
N
R
=
5 dB and the number of snapshots is 100. As can be seen in the figure, the 2D-DOA of the two targets are well localized and the PSA of the two targets are also well estimated.Figure 3
The estimation results of two targets with 500 trials.Figure4 shows the RMSE of 2D-DOA and PSA versus SNR, compared with CRB. In this simulation, we set L
=
200 and perform 1000 Monte Carlo trials. It is indicated in Figures 4(a)–4(d) that the performance of 2D-DOA estimation of the proposed algorithm becomes better with SNR increasing and there is a gap between proposed algorithm and CRB. We can also see that the performance of PSA estimations of the proposed algorithm is close to CRB. It is well known that DOA estimation performance is critically dependent on array aperture. Large array aperture can offer very accurate estimation. This result occurs because the proposed method does not make use of array aperture but only uses the information inside of the vector sensor. The polarization estimation performance is only dependent on the information inside of the vector sensor [17] and we fully make use of it.Figure 4
RMSE versus SNR.
(a)
Azimuth
(b)
Elevation
(c)
Auxiliary polarization angle
(d)
Polarization phase differenceFigure5 shows the RMSE of 2D-DOA and PSA versus the number of snapshots, compared with CRB. In this simulation, we set S
N
R
=
20dB and the other simulation conditions are the same as the above simulation ones. It is indicated in Figures 5(a)–5(d) that the performance of 2D-DOA estimation of the proposed algorithm becomes better in collaboration with L increasing and there is a gap between proposed algorithm and CRB. And the performance of PSA estimation of the proposed algorithm is close to CRB. The reason is the same as the above simulation.Figure 5
RMSE versus number of snapshots.
(a)
Azimuth
(b)
Elevation
(c)
Auxiliary polarization angle
(d)
Polarization phase difference


## 6. Conclusion

This paper has studied the 2D-DOA estimation problem for monostatic MIMO radar with arbitrary array. ESPRIT-like algorithm is used to calculate RIFs between the six subarrays within six-component EMVS. Then the electric-field and magnetic-field vectors are restored according to the five RIFs. Last, the 2D-DOA is obtained by the Pointing vector offered by the vector cross product between the electric-field and magnetic-field vectors. The proposed algorithm can be suitable for arbitrary array with electromagnetic vector received sensors at unknown locations. This advantage alleviates the nontrivial problem of spatial calibration. Moreover, it has low complexity because of no angle searching and the 2D-DOA estimation validity region is whole sphere.



---

*Source: 2952539-2016-08-17.xml*