US 12,392,855 B2 · Two-dimensional direction-of-arrival estimation method for coprime surface array based on virtual domain tensor filling

Abstract

Disclosed in the present invention is a two-dimensional direction-of-arrival estimation method for a coprime surface array based on virtual domain tensor filling, which mainly solves the problems of the loss of multi-dimensional signal structural information and the inability to fully utilize virtual domain statistics in the existing method. The steps thereof are as follows: constructing a coprime surface array; modeling a tensor of a received signal of the coprime surface array; constructing an augmented non-continuous virtual surface array based on cross-correlation tensor transformation of the coprime surface array; deriving a virtual domain tensor based on mirror extension of the non-continuous virtual surface array; dispersing contiguous missing elements by reconstructing the virtual domain tensor; filling the virtual domain tensor based on the minimization of a tensor kernel norm; and decomposing a filled virtual domain tensor to obtain a direction-of-arrival estimation result.

A two-dimensional direction-of-arrival estimation method for a coprime surface array based on virtual domain tensor filling, wherein the method comprises the following steps:

(1) configuring a receiving end with a coprime surface array by using 4MxMy+NxNy−1 physical antenna array elements, wherein Mx, Nx and My, Ny are a pair of coprime integers respectively;
decomposing the coprime surface array into two sparse uniform sub-surface arrays and , wherein contains 2Mx×2My antenna array elements, array element spacings in an x axial direction and a y axial direction are respectively Nxd and Nyd, includes Nx×Ny antenna array elements, array element spacings in the x axial direction and the y axial direction are respectively Mxd and Myd, and an unit interval d is taken as half of wavelength λ of an incident narrowband signal;
(2) if there are K far-field narrowband uncorrelated signal sources from {(θ1, φ1), (θ2, φ2), . . . , θK, φK)} directions, θk and φk are respectively an azimuth angle and an elevation angle of a kth incident signal source, k=1, 2, . . . , K, utilizing a three-dimensional tensor
X
ℙ
1
∈
ℂ
2
⁢
M
x
×
2
⁢
M
y
×
T
to express T sampling snapshot signals of a sparse uniform sub-surface array as follows:
X
ℙ
1
=
∑
k
=
1
K
a
x
(
ℙ
1
)
(
θ
k
,
φ
k
)
∘
a
y
(
ℙ
1
)
(
θ
k
,
φ
k
)
∘
s
k
+
N
ℙ
1
,
wherein sk=[sk,1, sk,2, . . . , sk,T]T is a multi-snapshot sampling signal waveform corresponding to the kth incident signal source, [⋅]T represents a transpose operation, ∘ represents an outer product of a vector, is a noise tensor independent of each signal source,
a
x
(
ℙ
1
)
(
θ
k
,
φ
k
)
⁢
and
⁢
a
y
(
ℙ
1
)
(
θ
k
,
φ
k
)
are respectively steering vectors of in the x axial direction and the y axial direction, correspond to a signal source with an incoming wave direction of (θk, φk), and are expressed as follows:
a
x
(
ℙ
1
)
(
θ
k
,
φ
k
)
=
[
1
,
e
-
j
⁢
π
⁢
x
ℙ
1
(
2
)
⁢
μ
k
,
…
,
e
-
j
⁢
π
⁢
x
ℙ
1
(
2
⁢
M
x
)
⁢
μ
k
]
T
,
a
y
(
ℙ
1
)
(
θ
k
,
φ
k
)
=
[
1
,
e
-
j
⁢
π
⁢
y
ℙ
1
(
2
)
⁢
v
k
,
…
,
e
-
j
⁢
π
⁢
y
ℙ
1
(
2
⁢
M
y
)
⁢
v
k
]
T
,
and wherein
{
x
ℙ
1
(
1
)
,
x
ℙ
1
(
2
)
,
…
,
x
ℙ
1
(
2
⁢
M
x
)
}
{
y
ℙ
1
(
1
)
,
y
ℙ
1
(
2
)
,
…
,
y
ℙ
1
(
2
⁢
M
y
)
}
represent respectively actual positions of physical antenna elements of the sparse uniform sub-surface array in the x axial direction and the y axial direction, and
x
ℙ
1
(
1
)
=
0
,
y
ℙ
1
(
1
)
=
0
,
μ
k
=
sin
⁡
(
φ
k
)
⁢
cos
⁡
(
θ
k
)
,
v
k
=
sin
⁡
(
φ
k
)
⁢
sin
⁡
(
θ
k
)
,
j
=
-
1
;
expressing the T sampled snapshot signals of the sparse uniform sub-surface array by another three-dimensional tensor
X
ℙ
2
∈
ℂ
N
x
×
N
y
×
T
as follows:
X
ℙ
2
=
∑
k
=
1
K
a
x
(
ℙ
2
)
(
θ
k
,
φ
k
)
∘
a
y
(
ℙ
2
)
(
θ
k
,
φ
k
)
∘
s
k
+
N
ℙ
2
,
wherein is a noise tensor independent of each signal source,
a
x
(
ℙ
2
)
(
θ
k
,
φ
k
)
⁢
and
⁢
⁢
a
y
(
ℙ
2
)
(
θ
k
,
φ
k
)
are respectively steering vectors of in the x axial direction and the y axial direction, correspond to a signal source with an incoming wave direction of (θk, φk), and are expressed as follows:
a
x
(
ℙ
2
)
(
θ
k
,
φ
k
)
=
[
1
,
e
-
j
⁢
π
⁢
x
ℙ
2
(
2
)
⁢
μ
k
,
…
,
e
-
j
⁢
π
⁢
x
ℙ
2
(
N
x
)
⁢
μ
k
]
T
,
a
y
(
ℙ
2
)
(
θ
k
,
φ
k
)
=
[
1
,
e
-
j
⁢
π
⁢
y
ℙ
2
(
2
)
⁢
v
k
,
…
,
e
-
j
⁢
π
⁢
y
ℙ
2
(
N
y
)
⁢
v
k
]
T
,
and wherein
{
x
ℙ
2
(
1
)
,
x
ℙ
2
(
2
)
,
…
,
x
ℙ
2
(
N
x
)
}
{
y
ℙ
2
(
1
)
,
y
ℙ
2
(
2
)
,
…
,
y
ℙ
2
(
N
y
)
}
represent respectively actual positions of physical antenna elements of the sparse uniform sub-surface array in the x axial direction and the y axial direction, and
x
ℙ
2
(
1
)
=
0
,
y
ℙ
2
(
1
)
=
0
;
obtaining a second-order cross-correlation tensor
R
ℙ
1
⁢
ℙ
2
∈
ℂ
2
⁢
M
x
×
2
⁢
M
y
×
N
x
×
N
y
by solving cross-correlation statistic of the three-dimensional tensors and :
R
ℙ
1
⁢
ℙ
2
=
E
[
〈
X
ℙ
1
,
X
ℙ
2
*
〉
3
]
=
∑
k
=
1
K
σ
k
2
⁢
a
x
(
ℙ
1
)
(
θ
k
,
φ
k
)
∘
a
y
(
ℙ
1
)
(
θ
k
,
φ
k
)
∘
a
x
(
ℙ
2
)
*
(
θ
k
,
φ
k
)
∘
a
y
(
ℙ
2
)
*
(
θ
k
,
φ
k
)
+
N
ℙ
1
⁢
ℙ
2
,
wherein σk2=E[sksk*] represents power of a kth incident signal source,
N
ℙ
1
⁢
ℙ
2
=
E
[
〈
N
ℙ
1
,
N
ℙ
2
*
〉
3
]
represents a cross-correlation noise tensor, <⋅,⋅>r represents a tensor contraction operation of two tensors along a rth dimension, E[⋅] represents a mathematical expectation operation, and (⋅)* represents a conjugation operation;
the cross-correlation noise tensor only has an element with a value σn2 in the (1, 1, 1, 1)th position, wherein σn2 represents a noise power, and elements in other positions have the same value 0;
(3) defining dimension sets J1={1,3}, J2={2,4}, and obtaining a virtual domain signal
U
W
∈
ℂ
2
⁢
M
x
⁢
N
x
×
2
⁢
M
y
⁢
N
y
by performing a tensor transformation of dimension merging on the cross-correlation tensor :
U
W
=
△
R
ℙ
1
⁢
ℙ
2
{
J
1
,
J
2
}
=
∑
k
=
1
K
σ
k
2
⁢
b
x
(
k
)
∘
b
y
(
k
)
,
wherein by respectively forming a difference set array on an exponential term,
b
x
(
k
)
=
a
x
(
ℙ
2
)
*
(
θ
k
,
φ
k
)
⊗
a
x
(
ℙ
1
)
(
θ
k
,
φ
k
)
⁢
and
⁢
b
y
(
k
)
=
a
y
(
ℙ
2
)
*
(
θ
k
,
φ
k
)
⊗
a
y
(
ℙ
1
)
(
θ
k
,
φ
k
)
configure a two-dimensional augmented virtual surface array along the x axial direction and the y axial direction, ⊗ represents a Kronecker product;
therefore, UW corresponds to a non-continuous virtual surface array W of size JWx×JWy, JWx=3MxNx−Mx−Nx+1, JWy=3MyNy−My−Ny+1, and the non-continuous virtual surface array W contains holes in an entire row and an entire column;
(4) configuring a virtual surface array W that mirrors the non-continuous virtual surface array W about a coordinate axis, and superimposing the W and W on a third dimension into a three-dimensional non-continuous virtual cubic array
ℚ
⁢
of
⁢
size
⁢
J
ℚ
x
×
J
ℚ
y
×
J
ℚ
z
,
wherein
⁢
J
ℚ
x
=
J
W
x
,
J
ℚ
y
=
J
W
y
⁢
and
⁢
J
ℚ
z
=
2
;
correspondingly, rearranging elements in a conjugate transposed signal UW* of the virtual domain signal UW to correspond to positions of virtual array elements in W, so as to obtain a virtual domain signal UW corresponding to the virtual surface array W;
superimposing UW and UW in the third dimension to obtain a virtual domain tensor corresponding to the non-continuous virtual cubic array , which is represented as:
U
ℚ
=
∑
k
=
1
K
σ
k
2
⁢
b
~
x
∘
b
~
y
(
k
)
∘
c
k
,
wherein {tilde over (b)}x(k) and {tilde over (b)}y(k) are respectively steering vectors of the non-continuous virtual cubic array on the x axial direction and the y axial direction, and correspond to the signal source with the incoming wave direction (θk, φk);
due to existence of the holes in , {tilde over (b)}x(k) and {tilde over (b)}y(k) respectively correspond to elements in hole positions in in the x axial direction and the y axial direction which are set to be zero,
c
k
=
[
1
,
e
-
j
⁢
π
⁡
(
-
(
M
x
⁢
M
y
+
M
x
+
M
y
)
⁢
μ
k
-
(
N
x
⁢
N
y
+
N
x
+
N
y
)
⁢
v
k
)
]
T
represents a mirror transformation factor vector corresponding to W and W;
since the non-continuous virtual surface array W contains the holes in the entire row and the entire column, the non-continuous virtual cubic array obtained by superimposing W with a mirror image part thereof W contains contiguous holes, corresponds to virtual field tensor of the non-continuous virtual cubic array , and thus contains contiguous missing elements;
(5) designing a translation window of size Px×Py×2 to select a sub-tensor of the virtual domain tensor , wherein contains elements of which indices are (1:
Px−1), (1:
Py−1) and (1:2) respectively in three dimensions of ;
then, translating the translation window by one element in turn along the x axial direction and the y axial direction, dividing into Lx×Ly sub-tensors, expressed as , sx=1, 2, . . . , Lx, sy=1, 2, . . . , Ly, wherein a value range of a size of the translation window is as follows:
2
≤
P
x
≤
J
ℚ
x
-
1
,
2
≤
P
y
≤
J
ℚ
y
-
1
,
and Lx, Ly, Px, Py satisfy the following relationship:
P
x
+
L
x
-
1
=
J
ℚ
x
,
P
y
+
L
y
-
1
=
J
ℚ
y
,
superimposing the sub-tensors with the same index subscript sy in a fourth dimension to obtain Ly four-dimensional tensors with Px×Py×2×Lx dimensions;
further, superimposing the Ly four-dimensional tensors in a fifth dimension to obtain a five-dimensional virtual domain tensor
T
ℚ
∈
ℂ
P
x
×
P
y
×
2
×
L
x
×
L
y
,
wherein the five-dimensional virtual domain tensor contains spatial angle information in the x axial direction and the y axial direction, spatial mirror transformation information, and spatial translation information in the x axial direction and the y axial direction;
defining dimension sets K1={1, 2}, K2={3}, K3={4, 5}, and merging through the dimensions of to obtain a three-dimensional reconfigured virtual domain tensor
K
ℚ
∈
ℂ
P
x
⁢
P
y
×
L
x
⁢
L
y
×
2
;
K
ℚ
=
Δ
T
ℚ
{
K
1
,
K
2
,
K
3
}
wherein the three dimensions of respectively represent the spatial angle information, the spatial translation information and the spatial mirror transformation information, thus, the contiguous missing elements in the original virtual domain tensor are randomly distributed to the three spatial dimensions contained in ;
(6) designing a virtual domain tensor filling optimization problem based on tensor kernel norm minimization:
min
K
_
ℚ
_

K
_
ℚ
_

*
⁢
s
.
t
.
P
Ω
_
(
K
_
ℚ
_
)
=
P
Ω
_
(
K
ℚ
)
,
wherein optimization variable
K
¯
ℚ
¯
∈
ℂ
P
x
⁢
P
y
×
L
x
⁢
L
y
×
2
is a filled virtual domain tensor corresponding to the virtual uniform cubic array , ∥⋅|* represents a tensor kernel norm, Ω represents a position index set of non-missing elements in , PΩ(⋅) represents a mapping of a tensor on Ω, wherein the virtual domain tensor filling optimization problem is solved to obtain ;
(7) expressing the filled virtual domain tensor as follows:
K
_
ℚ
_
=
∑
k
=
1
K
σ
k
2
⁢
p
k
∘
q
k
∘
c
k
,
wherein pk=dy(k)⊗dx(k), qk=gy(k)⊗gx(k) are the space factors of ,
d
x
(
k
)
=
[
e
-
j
⁢
π
⁡
(
-
M
x
⁢
N
x
+
M
x
)
⁢
μ
k
,
e
-
j
⁢
π
⁡
(
-
M
x
⁢
N
x
+
M
x
+
1
)
⁢
μ
k
,
…
,
e
-
j
⁢
π
⁡
(
2
⁢
M
x
⁢
N
x
-
N
x
)
⁢
μ
k
]
T
,
d
y
(
k
)
=
[
e
-
j
⁢
π
⁡
(
-
M
y
⁢
N
y
+
M
y
)
⁢
v
k
,
e
-
j
⁢
π
⁡
(
-
M
y
⁢
N
y
+
M
y
+
1
)
⁢
v
k
,
…
,
e
-
j
⁢
π
⁡
(
2
⁢
M
y
⁢
N
y
-
N
y
)
⁢
v
k
]
T
,
respectively represent steering vectors of the virtual uniform cubic array along the x axial direction and the y axial directions,
g
x
(
k
)
=
[
1
,
e
-
j
⁢
πμ
k
,
…
,
e
-
j
⁢
π
⁡
(
L
x
-
1
)
⁢
μ
k
]
T
,
g
y
(
k
)
=
[
1
,
e
-
j
⁢
π
⁢
v
k
,
…
,
e
-
j
⁢
π
⁡
(
L
y
-
1
)
⁢
v
k
]
T
,
are respectively spatial translation factor vectors corresponding to the x axial direction and the y axial direction in a process of intercepting the sub-tensor by the translation window;
performing canonical polyadic decomposition on the filled virtual domain tensor to obtain an estimated values of three factor vectors pk, qk and ck, which are expressed as {circumflex over (p)}k, {circumflex over (q)}k and ĉk;
and extracting angle parameters contained in exponential terms of {circumflex over (p)}k and {circumflex over (q)}k to obtain a two-dimensional direction-of-arrival estimation result ({circumflex over (θ)}k, φk), wherein
the receiving end with the coprime surface array steers, based on the two-dimensional direction-of-arrival estimation result ({circumflex over (θ)}k, {circumflex over (φ)}k), vectors corresponding to the signal source with the incoming wave direction (θk, φk).

1,414 more elements

Two-dimensional direction-of-arrival estimation method for coprime surface array based on virtual domain tensor filling

Life of this patent

Two-dimensional direction-of-arrival estimation method for coprime surface array based on virtual domain tensor filling

Life of this patent

Abstract

Claims4