SPATIAL SPECTRUM ESTIMATION METHOD WITH ENHANCED DEGREE-OF-FREEDOM BASED ON BLOCK SAMPLING TENSOR CONSTRUCTION FOR COPRIME PLANAR ARRAY

Abstract

Disclosed is a spatial spectrum estimation method with enhanced degree-of-freedom based on block sampling tensor construction for coprime planar array, which mainly solves the multi-dimensional information loss in signals and degree-of-freedom limitation in the existing methods and which is implemented by the following steps: constructing a coprime planar array; modeling block sampling tensors of the coprime planar array; deducing coarray statistics based on the block sampling cross-correlation tensor; obtaining block sampling coarray signals of a virtual uniform array; constructing a three-dimensional block sampling coarray tensor and its fourth-order auto-correlation statistics; constructing signal and noise subspaces based on fourth-order auto-correlation tensor decomposition; estimating a tensor spatial spectrum with enhanced degrees-of-freedom. In the present disclosure, the block sampling tensors of the coprime planar array is constructed, where a coarray tensor is deduced, to realize tensor spatial spectrum estimation with enhanced degrees-of-freedom by extracting signal-to-signal subspace features from the four-order self-correlation tensor.

Claims

1. A spatial spectrum estimation method with enhanced degree-of-freedom based on block sampling tensor construction for coprime planar array, characterized by comprising steps of: (1) constructing, by a receiving end, an architecture using 4M.sub.xM.sub.y+N.sub.xN.sub.y−1 physical antenna array elements according to the structure of the coprime planar array; wherein, M.sub.x, N.sub.x and M.sub.y, N.sub.y are a pair of coprime integers, respectively, and M.sub.x<N.sub.x, M.sub.y<N.sub.y; and the coprime array can be decomposed into two sparse uniform sub-arrays .sub.1 and .sub.2; (2) assuming that there are K far-field narrowband incoherent sources from direction of {(θ.sub.1,φ.sub.1), (θ.sub.2,φ.sub.2), . . . , (θ.sub.K,φ.sub.K)}, taking L sample snapshots as one block sample, denoted as T.sub.r (r=1, 2, . . . , R), R denoting a number of the block samples; wherein, within a sampling range of each block, the received signals of a sparse sub-array .sub.1 of the coprime planar array can be represented by one three-dimensional tensor .sub.1.sup.(r)∈.sup.2M.sup.x.sup.×2M.sup.y.sup.×L (r=1, 2, . . . , R) as: $1_{\left(r\right)}^{}=\underset{k=1}{\overset{K}{.Math.}}{a}_{Mx}\left({\theta }_k,{\phi }_k\right)\circ a_{My}\left({\theta }_k,{\phi }_k\right)\circ s_k+1_{},$ wherein, s.sub.k=s.sub.k,1, s.sub.k,2, . . . , s.sub.k,L].sup.T is a multiple-snapshot sampling signal waveform corresponding to a k.sup.th incident source, [⋅].sup.T represents a transposition operation, .Math. represents a vector outer product, .sub.1 is a noise tensor independent of each signal source, a.sub.Mx(θ.sub.k,φ.sub.k) and a.sub.My(θ.sub.k,φ.sub.k) are steering vectors of .sub.1 in x-axis and y-axis directions, respectively, which correspond to the k.sup.th source with a direction-of-arrival (θ.sub.k,φ.sub.k) and are expressed as: $\begin{array}{cc}{a}_{Mx}\left({\theta }_k,{\phi }_k\right)={\left[1,e^{-j\pi u_1^{\left(2\right)}\sin \left({\phi }_k\right)\cos \left({\theta }_k\right)},\mathrm{.Math.},e^{-j\pi u_1^{\left(2M_x\right)}\sin \left({\phi }_k\right)\cos \left({\theta }_k\right)}\right]}^T,a_{My}\left({\theta }_k,{\phi }_k\right)={\left[1,e^{-j\pi v_1^{\left(2\right)}\sin \left({\phi }_k\right)\sin \left({\theta }_k\right)},\mathrm{.Math.},e^{-j\pi v_1^{\left(2M_y\right)}\sin \left({\phi }_k\right)\sin \left({\theta }_k\right)}\right]}^T,& \end{array}$ wherein, u.sub.1.sup.(i.sup.1.sup.) (i.sub.1=1, 2, . . . , 2M.sub.x) and v.sub.1.sup.(i.sup.2.sup.) (i.sub.2=1, 2, . . . , 2M.sub.y) represent actual positions of the i.sub.1.sup.th and i.sub.2.sup.th physical antenna array elements of the sparse sub-array .sub.1 in the x-axis and y-axis directions, respectively, and u.sub.1.sup.(1)=0, v.sub.1.sup.(1)=0, j=√{square root over (−1)}; within the sampling range of each block, the received signals of the sparse sub-array .sub.2 can be represented by another three-dimensional tensor .sub.2.sup.(r)∈.sup.N.sup.x.sup.×N.sup.y.sup.×L (r=1, 2, . . . , R) as: $2_{\left(r\right)}^{}=\underset{k=1}{\overset{K}{.Math.}}{a}_{Nx}\left({\theta }_k,{\phi }_k\right)\circ a_{Ny}\left({\theta }_k,{\phi }_k\right)\circ s_k+2_{},$ wherein, .sub.2 is a noise tensor independent of each signal source, a.sub.Nx(θ.sub.k,φ.sub.k) and a.sub.Ny(θ.sub.k,φ.sub.k) are steering vectors of the sparse sub-array .sub.2 in the x-axis and y-axis directions, respectively, which correspond to the k.sup.th source with a direction-of-arrival (θ.sub.k,φ.sub.k) and are expressed as: $\begin{array}{cc}{a}_{\mathrm{Nx}}\left({\theta }_k,{\phi }_k\right)={\left[1,e^{-j\pi u_2^{\left(2\right)}\sin \left({\phi }_k\right)\cos \left({\theta }_k\right)},\mathrm{.Math.},e^{-j\pi u_2^{\left(N_x\right)}\sin \left({\phi }_k\right)\cos \left({\theta }_k\right)}\right]}^T,a_{\mathrm{Ny}}\left({\theta }_k,{\phi }_k\right)={\left[1,e^{-j\pi v_2^{\left(2\right)}\sin \left({\phi }_k\right)\sin \left({\theta }_k\right)},\mathrm{.Math.},e^{-j\pi v_2^{\left(N_y\right)}\sin \left({\phi }_k\right)\sin \left({\theta }_k\right)}\right]}^T,& \end{array}$ wherein, u.sub.2.sup.(i.sup.3.sup.) (i.sub.3=1, 2, . . . , N.sub.x) and v.sub.2.sup.(i.sup.4.sup.) (i.sub.4=1, 2, . . . , N.sub.y) represent actual positions of the i.sub.3.sup.th and i.sub.4.sup.th physical antenna array elements of the sparse sub-array .sub.2 in the x-axis and y-axis directions, respectively, and u.sub.2.sup.(1)=0, v.sub.2.sup.(1)=0; for one block sample T.sub.r (r=1, 2, . . . , R), a second-order cross-correlation tensor .sub.r∈.sup.2M.sup.x.sup.×2M.sup.y.sup.×N.sup.x.sup.×N.sup.y (r=1, 2, . . . , R) of the received tensor signals .sub.1.sup.(r) and .sub.2.sup.(r) (r=1, 2, . . . , R) of the sub-arrays .sub.1 and .sub.2 within the block sampling range is calculated, which is expressed as: $r_{}=\frac{1}{L}\underset{l=1}{\overset{L}{.Math.}}{1}_{\left(r\right)}^{}\left(l\right)\circ 2_{{\left(r\right)}^{*}}^{}\left(l\right),$ wherein, .sub.1.sup.(r)(l) and .sub.2.sup.(r)(l) respectively represent the l.sup.th slice in a direction of a third dimension (i.e., snapshot dimension), and (⋅)* represents a conjugation operation; (3) obtaining an augmented non-uniform virtual array from the cross-correlation tensor .sub.r, wherein a position of each virtual element is expressed as:
={(−M.sub.xn.sub.xd+N.sub.xm.sub.xd,−M.sub.yn.sub.yd+N.sub.ym.sub.yd)|0≤n.sub.x<N.sub.x−1,0≤m.sub.x≤2M.sub.x−1,0≤n.sub.y≤N.sub.y−1,0≤m.sub.y≤2M.sub.y−1} wherein, a unit spacing d is taken as half of an incident narrowband signal wavelength λ, that is, d=λ/2; dimensional sets .sub.1={1, 3} and .sub.2={2, 4} are defined, then the equivalent signals U.sub.r∈.sup.2M.sup.x.sup.N.sup.x.sup.×2M.sup.y.sup.N.sup.y (r=1, 2, . . . , R) of the augmented virtual array can be obtained by modulo {.sub.1, .sub.2} PARAFAC-based unfolding on the ideal value .sub.r (noise-free scene) of the cross-correlation tensor .sub.r, which is ideally expressed as: $U_r\stackrel{\Delta }{=}{r_{{\left\{{\mathrm{��}}_1,{\mathrm{��}}_2\right\}}_{}}}_{}=\underset{k=1}{\overset{K}{.Math.}}{\sigma }_k^2{a}_x\left({\theta }_k,{\phi }_k\right)\circ a_y\left({\theta }_k,{\phi }_k\right)$ wherein, a.sub.x(θ.sub.k,φ.sub.k)=a.sub.Nx*(θ.sub.k,φ.sub.k)⊕a.sub.Mx(θ.sub.k,φ.sub.k) and a.sub.y(θ.sub.k,φ.sub.k)=a.sub.Ny*(θ.sub.k,φ.sub.k)⊕a.sub.My(θ.sub.k,φ.sub.k) are steering vectors of the augmented virtual array in the x-axis and y-axis directions, which correspond to the k.sup.th source with a direction-of-arrival (θ.sub.k,φ.sub.k); σ.sub.k.sup.2 represents the power of the k.sup.th incident source; wherein, ⊕ represents a Kronecker product; and the tensor subscripts represent PARAFAC-based tensor unfolding; (4) comprising a continuous uniform virtual array with x-axis distribution from (−N.sub.x+1)d to (M.sub.xN.sub.x+M.sub.x−1)d and y-axis distribution from (−N.sub.y+1)d to (M.sub.yN.sub.y+M.sub.y−1)d in, wherein there are a total of V.sub.x×V.sub.y virtual array elements in , where V.sub.x=M.sub.xN.sub.x+M.sub.x+N.sub.x−1, V.sub.y=M.sub.yN.sub.y+M.sub.y+N.sub.y−1, is expressed as:
={(x,y)|x=p.sub.xd,y=p.sub.yd,−N.sub.x+1≤p.sub.x≤M.sub.xN.sub.x+M.sub.x−1,−N.sub.y+1≤p.sub.y≤M.sub.yN.sub.y+M.sub.y−1} by selecting the elements in the coarray signals U.sub.r corresponding to the positions of the virtual elements of , the block sampling equivalent signals Ũ.sub.r∈.sup.V.sup.x.sup.×V.sup.y (r=1, 2, . . . , R) of the virtual uniform array is obtained and expressed as: ${\tilde{U}}_r=\underset{k=1}{\overset{K}{.Math.}}{\sigma }_k^2{b}_x\left({\theta }_k,{\phi }_k\right)\circ b_y\left({\theta }_k,{\phi }_k\right),$ where b.sub.x(θ.sub.k,φ.sub.k)=[e.sup.−jπ(−N.sup.x.sup.+1)sin(φ.sup.k.sup.)cos(θ.sup.k.sup.), e.sup.−jπ(−N.sup.x.sup.+2)sin(φ.sup.k.sup.)cos(θ.sup.k.sup.), . . . , e.sup.−jπ(M.sup.x.sup.N.sup.x.sup.+M.sup.x.sup.−1)sin(φ.sup.k.sup.)cos(θ.sup.k.sup.)] and b.sub.y(θ.sub.k,φ.sub.k)=[e.sup.−jπ(−N.sup.y.sup.+1)sin(φ.sup.k.sup.)sin(θ.sup.k.sup.), e.sup.−jπ(−N.sup.y.sup.+2)sin(φ.sup.k.sup.)cos(θ.sup.k.sup.), . . . , e.sup.−jπ(M.sup.y.sup.N.sup.y.sup.+M.sup.y.sup.−1)sin(φ.sup.k.sup.)sin(θ.sup.k.sup.)] are steering vectors of the virtual uniform array in x-axis and y-axis directions, which correspond to the k.sup.th source with the direction-of-arrival (θ.sub.k,φ.sub.k); (5) according to the foregoing steps, taking R block samples T.sub.r (r=1, 2, . . . , R) to correspondently obtain R coarray signals Ũ.sub.r (r=1, 2, . . . , R), and superimposing the R coarray signals Ũ.sub.r (r=1, 2, . . . , R) in the third dimension to obtain a coarray tensor ∈.sup.V.sup.x.sup.×V.sup.y.sup.×R in which the third dimension represents equivalent sampling snapshots; calculating a fourth-order auto-correlation tensor ∈.sup.V.sup.x.sup.×V.sup.y.sup.V.sup.x.sup.×V.sup.y of the block sampling coarray tensor and expressing it as: $=\frac{1}{R}\underset{r=1}{\overset{R}{.Math.}}\left(r\right)\circ {*}^{}\left(r\right),$ wherein, (r) represents the r.sup.th slice of in a direction of the third dimension (i.e., the equivalent sampling snapshot dimension represented by block sampling); (6) performing CANDECOMP/PARACFAC decomposition on the fourth-order auto-correlation coarray tensor to extract multi-dimension features, the results of which are expressed as follows:
=Σ.sub.k=1.sup.K{tilde over (b)}.sub.x(θ.sub.k,φ.sub.k).Math.b.sub.y(θ.sub.k,φ.sub.k).Math.{tilde over (b)}.sub.x*(θ.sub.k,φ.sub.k).Math.{tilde over (b)}.sub.y*(θ.sub.k,φ.sub.k), wherein, {tilde over (b)}.sub.x(θ.sub.k,φ.sub.k) (k=1, 2, . . . , K) and {tilde over (b)}.sub.y(θ.sub.k,φ.sub.k) (k=1, 2, . . . , K) are factor vectors obtained by CANDECOMP/PARACFAC decomposition, which represent x-axis direction spatial information and y-axis direction spatial information, respectively; at this time, a theoretical maximum of the number K of the signal sources, which are distinguishable by the auto-correlation CANDECOMP/PARACFAC decomposition, exceeds the actual number of physical array elements; further, a noise subspace .sub.s∈.sup.V.sup.x.sup.V.sup.y.sup.×K is constructed and expressed as:
=orth([{tilde over (b)}.sub.x(θ.sub.1,φ.sub.1).Math.{tilde over (b)}.sub.y(θ.sub.1,φ.sub.1),{tilde over (b)}.sub.x(θ.sub.2,φ.sub.2).Math.{tilde over (b)}.sub.y(θ.sub.2,φ.sub.2), . . . ,{tilde over (b)}.sub.x(θ.sub.K,φ.sub.K).Math.{tilde over (b)}.sub.y(θ.sub.K,φ.sub.K)]), wherein, orth(⋅) represents a matrix orthogonalization operation; further, .sub.n∈.sup.V.sup.x.sup.V.sup.y.sup.×(V.sup.x.sup.V.sup.y.sup.−K) represents a noise subspace, then .sub.s and .sub.n have a following relationship:
.sub.n.sub.n.sup.H=I−.sub.s.sub.s.sup.H, wherein, I represents a unit matrix; (⋅).sup.H represents a conjugate transposition operation; (7) constructing a tensor spatial spectrum function with enhanced degree-of-freedom according to the obtained signal subspace and the noise subspace, to obtain the spatial spectrum estimation corresponding to the two-dimensional direction-of-arrival.

2. The spatial spectrum estimation method with enhanced degree-of-freedom based on block sampling tensor construction for coprime planar array of claim 1, wherein a structure of the coprime planar array described in step (1) can be described as: a pair of spare uniform planar sub-arrays .sub.1 and .sub.2 are constructed on a planar coordinate system xoy, wherein .sub.1 contains 2M.sub.x×2M.sub.y antenna array elements, inter-element spacings in the x-axis direction and the y-axis direction are N.sub.xd and N.sub.yd, respectively, the position coordinates of which on xoy are {(N.sub.xdm.sub.x, N.sub.ydm.sub.y), m.sub.x=0, 1, . . . , 2M.sub.x−1, m.sub.y=0, 1, . . . , 2M.sub.y−1}; .sub.2 contains N.sub.x×N.sub.y antenna array elements, inter-element spacings in the x-axis direction and the y-axis direction are M.sub.xd and M.sub.yd, respectively, the position coordinates of which on xoy are {(M.sub.xdn.sub.x, M.sub.ydn.sub.y), n.sub.x=0, 1, . . . , N.sub.x−1, n.sub.y=0, 1, . . . , N.sub.y−1; wherein, M.sub.x, N.sub.x and M.sub.y, N.sub.y are a pair of coprime integers, respectively, and M.sub.x<N.sub.x, M.sub.y<N.sub.y; .sub.1 and .sub.2 are combined in sub-arrays by means of overlapping array elements at the coordinate (0,0), to obtain a coprime planar array that actually contains 4M.sub.xM.sub.y+N.sub.xN.sub.y−1 physical antenna array elements.

3. The spatial spectrum estimation method with enhanced degree-of-freedom based on block sampling tensor construction for coprime planar array of claim 1, wherein the cross-correlation tensor .sub.r described in step (3) can be ideally modeled as (noise-free scene):
=Σ.sub.k=1.sup.Kσ.sub.k.sup.2a.sub.Mx(θ.sub.k,φ.sub.k).Math.a.sub.My(θ.sub.k,φ.sub.k).Math.a.sub.Nx*(θ.sub.k,φ.sub.k).Math.a.sub.Ny*(θ.sub.k,φ.sub.k), wherein, in .sub.r, a.sub.mx(θ.sub.k,φ.sub.k).Math.a.sub.Nx*(θ.sub.k,φ.sub.k) is equivalent to an augmented coarray along the x-axis; a.sub.My(θ.sub.k,φ.sub.k).Math.a.sub.Ny*(θ.sub.k,φ.sub.k) is equivalent to an augmented coarray along the y-axis, such that a non-uniform virtual array can be obtained.

4. The spatial spectrum estimation method with enhanced degree-of-freedom based on block sampling tensor construction for coprime planar array of claim 1, wherein as described in step (5), the coarray signals Ũ.sub.r (r=1, 2, . . . , R) corresponding to R block samples T.sub.r (r=1, 2, . . . , R) is constructed, and Ũ.sub.r (r=1, 2, . . . , R) is superimposed along the third dimension to obtain a coarray tensor ∈.sup.V.sup.x.sup.×V.sup.y.sup.×R, the first two dimensions of the coarray tensor represent the spatial information of the virtual uniform array in x-axis and y-axis directions; the third dimension represents the equivalent sampling snapshot constructed by block sampling; the coarray tensor has the same structure as that of the actual received tensor signals .sub.1.sup.(r) and .sub.2.sup.(r) of the coprime planar array; for the coarray tensor , the fourth-order auto-correlation tensor can be directly calculated without need to introduce spatial smoothing process to compensate for a rank deficiency problem caused by a single snapshot of the coarray signals.

5. The spatial spectrum estimation method with enhanced degree-of-freedom based on block sampling tensor construction for coprime planar array of claim 1, wherein the CANDECOMP/PARACFAC decomposition for the fourth-order auto-correlation tensor described in step (6) follows a uniqueness condition as follows:
.sub.rank({tilde over (B)}.sub.x)+.sub.rank({tilde over (B)}.sub.y)+.sub.rank({tilde over (B)}.sub.x*)+.sub.rank({tilde over (B)}.sub.y*)≥2K+3, wherein, .sub.rank(⋅) represents Kruskal rank of the matrix, {tilde over (B)}.sub.x=[{tilde over (b)}.sub.x(θ.sub.1,φ.sub.1), {tilde over (b)}.sub.x(θ.sub.2,φ.sub.2), . . . , {tilde over (b)}.sub.x(θ.sub.K,φ.sub.K)] and {tilde over (B)}.sub.y=[{tilde over (b)}.sub.y(θ.sub.1,φ.sub.1), {tilde over (b)}.sub.y(θ.sub.2,φ.sub.2), . . . , {tilde over (b)}.sub.y(θ.sub.K,φ.sub.K)] represent factor sub-matrices, and .sub.rank({tilde over (B)}.sub.x)=min(V.sub.x, K), .sub.rank({tilde over (B)}.sub.y)=min(V.sub.y, K), .sub.rank({tilde over (B)}.sub.x*)=min(V.sub.x, K), .sub.rank({tilde over (B)}.sub.y*)=min(V.sub.y, K), min(⋅) represents minimum taking operation; therefore, the uniqueness condition for the CANDECOMP/PARACFAC decomposition is transformed into:
2min(V.sub.x,K)+2min(V.sub.y,K)≥2K+3, according to the above inequality, the number K of the distinguishable sources is greater than the number of the actual physical array elements, the maximum value of K is $.Math.\frac{2\left(V_x+V_y\right)-3}{2}.Math.,$ and └⋅┘ represents a rounding operation.

6. The spatial spectrum estimation method with enhanced degree-of-freedom based on block sampling tensor construction for coprime planar array of claim 1, wherein the signal and noise subspaces obtained by the fourth-order auto-correlation coarray tensor CANDECOMP/PARACFAC decomposition are utilized to construct the tensor spatial spectrum function in step (7); a two-dimensional direction-of-arrival ({tilde over (θ)}, {tilde over (φ)}), {tilde over (θ)}∈[−90°, 90°], {tilde over (φ)}∈[0°,180° ] for spectrum peak search are defined at first, and the steering information ({tilde over (θ)}, {tilde over (φ)})∈.sup.V.sup.x.sup.V.sup.y corresponding to the virtual uniform array is constructed, which is expressed as:
({tilde over (θ)},{tilde over (φ)})=b.sub.x({tilde over (θ)},{tilde over (φ)}).Math.b.sub.y({tilde over (θ)},{tilde over (φ)}), the tensor spatial spectrum function ({tilde over (θ)}, {tilde over (φ)}) based on the noise subspace is expressed as follows: $\left(\tilde{\theta },\tilde{\phi }\right)=,$ thus, the tensor spatial spectrum with enhanced degree-of-freedom corresponding to the two-dimensional search direction-of-arrival ({tilde over (θ)}, {tilde over (φ)}) is obtained.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0044] FIG. 1 is a block diagram of the overall flow of the present disclosure.

[0045] FIG. 2 is a schematic diagram of the structure of the coprime planar array in the present disclosure.

[0046] FIG. 3 is a schematic diagram of the structure of an augmented virtual array derived by the present disclosure.

DESCRIPTION OF EMBODIMENTS

[0047] Hereinafter, the technical solution of the present disclosure will be further described in detail with reference to the accompanying drawings.

[0048] In order to solve the problems of loss of signal multi-dimensional spatial structural information and limited degree-of-freedom performance in existing methods, the present disclosure provides a spatial spectrum estimation method with enhanced degree-of-freedom based on the coprime planar array block sampling tensor construction. Through the statistical analysis of the block sampling tensor of the coprime planar array, the coarray statistics based on the block sampling tensor statistics are derived, and the coarray tensor with equivalent sampling snapshots is constructed; the fourth-order auto-correlation coarray tensor is decomposed by CANDECOMP/PARACFAC to obtain the signal and noise subspaces without need to introduce a spatial smoothing process, thereby constructing the tensor spatial spectrum function with enhanced degree-of-freedom. Referring to FIG. 1, the implementation steps of the present disclosure are as follows:

[0049] Step 1: constructing a coprime planar array. At a receiving end, 4M.sub.xM.sub.y+N.sub.xN.sub.y−1 physical antenna array elements are used to construct a coprime planar array, as shown in FIG. 2: a pair of sparse uniform plane sub-arrays .sub.1 and .sub.2 are constructed on the plane coordinate system xoy, where .sub.1 contains 2M.sub.x×2M.sub.y antenna array elements, and the inter-elements spacing in the x-axis direction and the y-axis direction are N.sub.xd and N.sub.yd, respectively, and the position coordinates on xoy thereof are {(N.sub.xdm.sub.x, N.sub.ydm.sub.y), m.sub.x=0, 1, . . . , 2M.sub.x=1, m.sub.y=0, 1, . . . , 2M.sub.y−1}; .sub.2 contains N.sub.x×N.sub.y antenna array elements, and the inter-elements spacings in the x-axis direction and the y-axis direction are M.sub.xd and M.sub.yd, respectively, and the position coordinates on xoy thereof are {(M.sub.xdn.sub.x, M.sub.ydn.sub.y), n.sub.x=0, 1, . . . , N.sub.x−1, n.sub.y=0, 1, . . . , N.sub.y−1}; wherein, M.sub.x, N.sub.x and M.sub.y, N.sub.y are respectively a pair of coprime integers, and M.sub.x<N.sub.x, M.sub.y<N.sub.y; the unit spacing d is taken as half of the incident narrowband signal wavelength λ, that is, λ/2; .sub.1 and .sub.2 are combined in sub-arrays by means of overlapping array elements at the coordinate (0,0), to obtain a coprime planar array that actually contains 4M.sub.xM.sub.y+N.sub.xN.sub.y−1 of physical antenna array elements;

[0050] Step 2: modeling block sampling tensors of a coprime planar array; assuming that there are K far-field narrowband incoherent sources from the direction of {(θ.sub.p,φ.sub.1), (θ.sub.2,φ.sub.2), . . . , (θ.sub.K,φ.sub.K)} and taking L continue time sampling snapshots as a block sample, denoted as T.sub.r (r=1, 2, . . . , R), wherein R is the number of block samples; within the sampling range of each block, the sampling snapshot signals of the sparse sub-array .sub.1 of the coprime planar array are superimposed in the third dimension to obtain a three-dimensional block sampling tensor .sub.1.sup.(r)∈.sup.2M.sup.x.sup.×2M.sup.y.sup.×L (r=1, 2, . . . , R), which is expressed as:

[00009]$1_{\left(r\right)}^{}=\underset{k=1}{\overset{K}{.Math.}}{a}_{Mx}\left({\theta }_k,{\phi }_k\right)\circ a_{My}\left({\theta }_k,{\phi }_k\right)\circ s_k+1_{},$

[0051] wherein, s.sub.k=[s.sub.k,1, s.sub.k,2, . . . , s.sub.k,L].sup.T is the multi-snapshot signal waveform corresponding to the k.sup.th incident signal source, [⋅].sup.T represents a transposition operation, and .Math. represents the outer vector product, .sub.1 is the noise tensor independent of each source, a.sub.Mx(θ.sub.k,φ.sub.k) and a.sub.My(θ.sub.k,φ.sub.k) are respectively the steering vectors of .sub.1 in the x-axis and y-axis directions, corresponding to the source with direction-of-arrival (θ.sub.k,φ.sub.k), which are expressed as:

[00010]$\begin{array}{cc}{a}_{Mx}\left({\theta }_k,{\phi }_k\right)={\left[1,e^{-j\pi u_1^{\left(2\right)}\sin \left({\phi }_k\right)\cos \left({\theta }_k\right)},\mathrm{.Math.},e^{-j\pi u_1^{\left(2M_x\right)}\sin \left({\phi }_k\right)\cos \left({\theta }_k\right)}\right]}^T,a_{My}\left({\theta }_k,{\phi }_k\right)={\left[1,e^{-j\pi v_1^{\left(2\right)}\sin \left({\phi }_k\right)\sin \left({\theta }_k\right)},\mathrm{.Math.},e^{-j\pi v_1^{\left(2M_y\right)}\sin \left({\phi }_k\right)\sin \left({\theta }_k\right)}\right]}^T,& \end{array}$

[0052] wherein u.sub.1.sup.(i.sup.1.sup.) (i.sub.1=1, 2, . . . , 2M.sub.x) and v.sub.1.sup.(i.sup.2.sup.) (i.sub.2=1, 2, . . . , 2M.sub.y) represent the actual positions of the i.sub.1.sup.th and i.sub.2.sup.th physical antenna array elements of the sparse sub-array .sub.1 in the x-axis and y-axis directions, respectively, and u.sub.i.sup.(1)=0, v.sub.1.sup.(1)=0, j=√{square root over (−1)}.

[0053] Similarly, one block sampling signal of the sparse sub-array .sub.2 can be represented by another three-dimensional tensor .sub.2.sup.(r)∈.sup.N.sup.x.sup.×N.sup.y.sup.×L (r=1, 2, . . . , R) as:

[00011]$2_{\left(r\right)}^{}=\underset{k=1}{\overset{K}{.Math.}}{a}_{Nx}\left({\theta }_k,{\phi }_k\right)\circ a_{Ny}\left({\theta }_k,{\phi }_k\right)\circ s_k+2_{},$

[0054] wherein, .sub.2 is the noise tensor independent of each source, a.sub.Nx(θ.sub.k,φ.sub.k) and a.sub.Ny(θ.sub.k,φ.sub.k) are the steering vectors of the sparse sub-array .sub.2 in the x-axis and y-axis directions respectively, corresponding to the source with direction-of-arrival (θ.sub.k,φ.sub.k), which are expressed as:

[00012]$\begin{array}{cc}{a}_{Nx}\left({\theta }_k,{\phi }_k\right)={\left[1,e^{-j\pi u_2^{\left(2\right)}s\mathrm{in}\left({\phi }_k\right)\cos \left({\theta }_k\right)},\mathrm{.Math.},e^{-j\pi u_2^{\left(N_x\right)}\sin \left({\phi }_k\right)\cos \left({\theta }_k\right)}\right]}^T,a_{Ny}\left({\theta }_k,{\phi }_k\right)={\left[1,e^{-j\pi v_2^{\left(2\right)}\sin \left({\phi }_k\right)\sin \left({\theta }_k\right)},\mathrm{.Math.},e^{-j\pi v_2^{\left(N_y\right)}\sin \left({\phi }_k\right)\sin \left({\theta }_k\right)}\right]}^T,& \end{array}$

[0055] wherein u.sub.2.sup.(i.sup.3.sup.) (i.sub.3=1, 2, . . . , N.sub.x) and v.sub.2.sup.(i.sup.4.sup.) (i.sub.4=1, 2, . . . , N.sub.y) represent the actual positions of the i.sub.3.sup.th and i.sub.4.sup.th physical antenna array elements of the sparse sub-array .sub.2 in the x-axis and y-axis directions, respectively, and u.sub.2.sup.(1)=0, v.sub.2.sup.(1)=0.

[0056] For a block sample T.sub.r (r=1, 2, . . . , R), the cross-correlation statistics of the tensor signals .sub.1.sup.(r) and .sub.2.sup.(r) (r=1, 2, . . . , R) of the sub-arrays .sub.1 and .sub.2 within the block sampling range are calculated to obtain a second-order cross-correlation tensor .sub.r∈.sup.2M.sup.x.sup.×2M.sup.y.sup.×N.sup.x.sup.×N.sup.y (r=1, 2, . . . , R) with four-dimensional space information, which is expressed as:

[00013]$\begin{array}{cc}{r}_{}=\frac{1}{L}\underset{l=1}{\overset{L}{.Math.}}{1}_{\left(r\right)}^{}\left(l\right)\circ 2_{{\left(r\right)}^{*}}^{}\left(l\right),& \end{array}$

[0057] wherein, .sub.1.sup.(r)(l) and .sub.2.sup.(r)(l) respectively represent the l.sup.th slice of .sub.1.sup.(r) and .sub.2.sup.(r) in the direction of the third dimension (i.e., snapshot dimension), and (⋅)* represents a conjugation operation;

[0058] Step 3: deducing coarray signals based on the cross-correlation statistics of the block sampling tensor signals. The second-order cross-correlation tensor .sub.r of the bock sampling received tensor signal of the two sub-arrays in the coprime planar array can be ideally modeled (noise-free scene) as:

.sub.r=Σ.sub.k=1.sup.Kσ.sub.k.sup.2a.sub.Mx(θ.sub.k,φ.sub.k).Math.a.sub.My(θ.sub.k,φ.sub.k).Math.a.sub.Nx*(θ.sub.k,φ.sub.k).Math.a.sub.Ny*(θ.sub.k,φ.sub.k),

[0059] wherein, σ.sub.k.sup.2 represents the power of the k.sup.th incident signal source; at this time, in .sub.r, a.sub.Mx(θ.sub.k,φ.sub.k).Math.a.sub.Nx*(θ.sub.k,φ.sub.k) is equivalent to an augmented virtual domain along the x-axis, and a.sub.My(θ.sub.k,φ.sub.k).Math.a.sub.Ny*(θ.sub.k,φ.sub.k) is equivalent to an augmented virtual domain along the y-axis, thus an augmented non-uniform virtual array S can be obtained as shown in FIG. 3, where the position of each virtual array element is expressed as:

={(−M.sub.xn.sub.xd+N.sub.xm.sub.xd,−M.sub.yn.sub.yd+N.sub.ym.sub.yd)|0≤n.sub.x≤N.sub.x−1,0≤m.sub.x≤2M.sub.x−1,0≤n.sub.y≤N.sub.y−1,0≤m.sub.y≤2M.sub.y−1}.

[0060] In order to obtain the equivalent signals corresponding to the augmented virtual array, the first and third dimensions representing the spatial information of x-axis direction in the cross-correlation tensor .sub.r are merged into one dimension, and the second and fourth dimensions representing the spatial information of y-axis direction are merged into another dimension. The dimensional merging of tensors can be realized through the PARAFAC-based unfolding. Taking a four-dimensional tensor ∈.sup.I.sup.1.sup.×I.sup.2.sup.×I.sup.3.sup.×I.sup.4=Σ.sub.p=1.sup.Pb.sub.11.Math.b.sub.12.Math.b.sub.21.Math.b.sub.22 as an example, the dimension sets .sub.1={1, 2} and .sub.2={3, 4} are defined, then the modulo{.sub.1, .sub.2} PARAFAC-based unfolding of is as follows:

[00014]$B\in {ℂ}^{I_1{I}_2\times I_3{I}_4}\stackrel{\Delta }{=}{\left\{{\mathrm{��}}_1,{\mathrm{��}}_2\right\}}_{}=\underset{p=1}{\overset{P}{.Math.}}{b}_1\circ b_2,$

[0061] wherein, the tensor subscript represents the tensor PARAFAC-based unfolding, b.sub.1=b.sub.12.Math.b.sub.11 and b.sub.2=b.sub.2.Math.b.sub.21 represent the factor vectors of the two dimensions after unfolding; wherein, .Math. represents the Kronecker product. Therefore, the dimensional sets .sub.1={1, 3} and .sub.2={2, 4} are defined, and an equivalent received signal U.sub.r∈.sup.2M.sup.x.sup.N.sup.x.sup.×2M.sup.y.sup.N.sup.y (r=1, 2, . . . , R) of the augmented virtual arrays can be obtained by modulo {} PARAFAC-based unfolding on the cross-correlation tensor .sub.r, which is expressed as:

U.sub.r=Σ.sub.k=1.sup.Kσ.sub.k.sup.2a.sub.x(θ.sub.k,φ.sub.k).Math.a.sub.y(θ.sub.k,φ.sub.k),

[0062] wherein, a.sub.x(θ.sub.k,φ.sub.k)=a.sub.Nx*(θ.sub.k,φ.sub.k).Math.a.sub.Mx(θ.sub.k,φ.sub.k) and a.sub.y(θ.sub.k,φ.sub.k)=a.sub.Ny*(θ.sub.k,φ.sub.k).Math.a.sub.My(θ.sub.k,φ.sub.k) are steering vectors of the augmented virtual array along in x-axis and y-axis directions, which correspond to the k.sup.th signal source with direction-of-arrival (θ.sub.k,φ.sub.k);

[0063] Step 4: obtaining the block sampling coarray signals of a virtual uniform array. includes a virtual uniform array with x-axis distribution from (−N.sub.x+1)d to (M.sub.xN.sub.x+M.sub.x−1)d and y-axis distribution from (−N.sub.y+1)d to (M.sub.yN.sub.y+M.sub.y−1)d. There are a total of V.sub.x×V.sub.y virtual elements in , where V.sub.x=M.sub.xN.sub.x+M.sub.x+N.sub.x−1, V.sub.y=M.sub.yN.sub.y+M.sub.y+N.sub.y−1; the structure of the virtual uniform array is shown in the dotted box in FIG. 3, and is expressed as:

={(x,y)|x=p.sub.xd,y=p.sub.yd,−N.sub.x+1≤p.sub.x≤M.sub.xN.sub.x+M.sub.x−1,

−N.sub.y+1≤p.sub.y≤M.sub.yN.sub.y+M.sub.y−1}.

[0064] By selecting the elements in the equivalent signals U.sub.r corresponding to the positions of the virtual elements of of the augmented virtual array , the block sampling equivalent signals Ũ.sub.r∈.sup.V.sup.x.sup.×V.sup.y (r=1, 2, . . . , R) of the virtual uniform array can be obtained:

Ũ.sub.r=Σ.sub.k=1.sup.Kσ.sub.k.sup.2b.sub.x(θ.sub.k,φ.sub.k).Math.b.sub.y(θ.sub.k,φ.sub.k),

[0065] wherein,

[0066] b.sub.x(θ.sub.k,φ.sub.k)=[e.sup.−jπ(−N.sup.x.sup.+1)sin(φ.sup.k.sup.)cos(θ.sup.k.sup.), e.sup.−jπ(−N.sup.x.sup.+2)sin(φ.sup.k.sup.)cos(θ.sup.k.sup.), . . . , e.sup.−jπ(M.sup.x.sup.N.sup.x.sup.+M.sup.x.sup.−1)sin(φ.sup.k.sup.)cos(θ.sup.k.sup.)] and b.sub.y(θ.sub.k,φ.sub.k)=[e.sup.−jπ(−N.sup.y.sup.+1)sin(φ.sup.k.sup.)sin(θ.sup.k.sup.),

[0067] e.sup.−jπ(−N.sup.y.sup.+2)sin(φ.sup.k.sup.)cos(θ.sup.k.sup.), . . . , e.sup.−jπ(M.sup.y.sup.N.sup.y.sup.+M.sup.y.sup.−1)sin(φ.sup.k.sup.)sin(θ.sup.k.sup.)] are steering vectors of the virtual uniform array in the x-axis and y-axis directions, corresponding to the k.sup.th source with the direction-of-arrival (θ.sub.k,φ.sub.k);

[0068] Step 5: constructing a three-dimensional block sampling coarraytensor and its fourth-order auto-correlation statistics. According to the foregoing steps, R block samples T.sub.r (r=1, 2, . . . , R) are taken to correspondently obtain R coarray signals Ũ.sub.r (r=1, 2, . . . , R) and these R coarray signals Ũ.sub.r are superimposed in the third dimension to obtain a three-dimensional tensor ∈.sup.V.sup.x.sup.×V.sup.y.sup.×R. The first two dimensions of the coarray tensor represent the spatial information of the virtual uniform array in the x-axis and y-axis directions, and the third dimension represents the equivalent snapshots constructed by block sampling. It can be seen that the coarray tensor has the same structure as that of the coprime planar array which actually receives tensor signals .sub.1.sup.(r) and .sub.2.sup.(r). For the coarray tensor , the fourth-order auto-correlation tensor can be directly calculated, without the need to introduce a spatial smoothing process to compensate for the rank deficiency problem caused by the single-block shooting of the coarray signals. The fourth-order auto-correlation tensor ∈.sup.V.sup.x.sup.×V.sup.y.sup.×V.sup.x.sup.×V.sup.y of the block sampling coarray tensor is calculated and expressed as:

[00015]$=\frac{1}{R}\underset{r=1}{\overset{R}{.Math.}}\left(r\right)\circ {*}^{}\left(r\right),$

[0069] wherein, (r) represents the r.sup.th slice of in the direction of the third dimension (i.e., the equivalent snapshot dimension represented by block sampling);

[0070] Step 6: constructing the signal-to-noise subspace based on the fourth-order auto-correlation coarray tensor decomposition. In order to construct the tensor spatial spectrum, the fourth-order auto-correlation tensor is subjected to CANDECOMP/PARACFAC decomposition to extract multi-dimensional features, and the result is expressed as follows:

=Σ.sub.k=1.sup.K{tilde over (b)}.sub.x(θ.sub.k,φ.sub.k).Math.{tilde over (b)}.sub.y(θ.sub.k,φ.sub.k).Math.{tilde over (b)}.sub.x*(θ.sub.k,φ.sub.k).Math.{tilde over (b)}.sub.y(θ.sub.k,φ.sub.k),

[0071] wherein, {tilde over (b)}.sub.x(θ.sub.k,φ.sub.k) (k=1, 2, . . . , K) and {tilde over (b)}.sub.y(θ.sub.k,φ.sub.k) (k=1, 2, . . . , K) are the factor vectors obtained by CANDECOMP/PARACFAC decomposition, which respectively represent the spatial information in the x-axis direction and the y-axis direction; {tilde over (B)}.sub.x=[{tilde over (b)}.sub.x(θ.sub.1,φ.sub.1), {tilde over (b)}.sub.x(θ.sub.2,φ.sub.2), . . . , {tilde over (b)}.sub.x(θ.sub.K,φ.sub.K)] and {tilde over (B)}.sub.y=[{tilde over (b)}.sub.y(θ.sub.1,φ.sub.1), {tilde over (b)}.sub.y(θ.sub.2,φ.sub.2), . . . , {tilde over (b)}.sub.y(θ.sub.K,φ.sub.K)] represent the factor sub-matrices. At this time, CANDECOMP/PARACFAC decomposition follows the uniqueness condition as follows:

.sub.rank({tilde over (B)}.sub.x)+.sub.rank({tilde over (B)}.sub.y)+.sub.rank({tilde over (B)}.sub.x)+.sub.rank({tilde over (B)}.sub.y*)≥2K+3,

[0072] wherein, .sub.rank(⋅) represents the Kruskal rank of the matrix, and .sub.rank({tilde over (B)}.sub.x)=min(V.sub.x, K), .sub.rank({tilde over (B)}.sub.y)=min(V.sub.y,K), .sub.rank({tilde over (B)}.sub.x*)=min(V.sub.x, K), rank({tilde over (B)}.sub.y*)=min(V.sub.y,K), and min(⋅) represents the minimum operation. Therefore, the above unique decomposition conditions can be transformed into:

2min(V.sub.x,K)+2min(V.sub.y,K)≥2K+3.

[0073] It can be seen from the above inequality that the number of distinguishable incident sources K of the method proposed in the present disclosure is greater than the number of actual physical array elements, and the maximum value of K is

[00016]$.Math.\frac{2\left(V_x+V_y\right)-3}{2}.Math.,$

and └⋅┘ represents a rounding operation. Furthermore, the multi-dimensional features obtained by tensor decomposition are used to construct the signal subspace ∈.sup.V.sup.x.sup.V.sup.y.sup.×K:

=orth([{tilde over (b)}.sub.x(θ.sub.1,φ.sub.1).Math.{tilde over (b)}.sub.y(θ.sub.1,φ.sub.1),{tilde over (b)}.sub.x(θ.sub.2,φ.sub.2).Math.{tilde over (b)}.sub.y(θ.sub.2,φ.sub.2), . . . ,{tilde over (b)}.sub.x(θ.sub.K,φ.sub.K).Math.{tilde over (b)}.sub.y(θ.sub.K,φ.sub.K)]),

[0074] wherein, orth(⋅) represents the matrix orthogonalization operation; .sub.n∈.sup.V.sup.x.sup.V.sup.y.sup.×(V.sup.x.sup.V.sup.y.sup.−K) represents the noise subspace, then .sub.s and .sub.n have the following relationship:

.sub.n.sub.n.sup.H=I−.sub.s.sub.s.sup.H,

[0075] wherein, I represents the unit matrix; (⋅).sup.H represents the conjugate transposition operation;

[0076] Step 7: estimating a tensor spatial spectrum with enhanced degree-of-freedom. The two-dimensional directions of arrival ({tilde over (θ)},{tilde over (φ)}), {tilde over (θ)}∈[−90°, 90°], {tilde over (φ)}∈[0°,180° ] for spectrum peak search are defined and the steering information ({tilde over (θ)},{tilde over (φ)})∈.sup.V.sup.x.sup.V.sup.y corresponding to the virtual uniform array is constructed, which is expressed as:

({tilde over (θ)},φ)=b.sub.x({tilde over (θ)},{tilde over (φ)}).Math.b.sub.y({tilde over (θ)},{tilde over (φ)}).

[0077] The tensor spatial spectrum function ({tilde over (θ)}, {tilde over (φ)}) based on the noise subspace is expressed as follows:

[00017]$\left(\tilde{\theta },\tilde{\phi }\right)=$

[0078] thus, the tensor spatial spectrum with enhanced degree-of-freedom corresponding to the two-dimensional search direction-of-arrival ({tilde over (θ)}, {tilde over (φ)}) is obtained.

[0079] In summary, the present disclosure fully considers the multi-dimensional information structure of the coprime planar array signal, uses block sampling tensor signal modeling, constructs a virtual domain tensor signal with equivalent sampling time sequence information, and further uses tensor decomposition to extract the multi-dimensional feature of the fourth-order statistics of the block sampling coarray tensor to construct a signal-to-noise subspace based on the block sampling coarray tensor, and establish the correlation between the block sampling coarray tensor signal and the tensor spatial spectrum of the coprime planar array; at the same time, the present disclosure obtains a coarray tensor with a three-dimensional information structure through the block sample construction, thereby avoiding the need of the introduction of a spatial smoothing process in order to solve the rank deficiency problem resulting from single-block shooting of the coarray signals; therefore, the advantages of the degree-of-freedom brought by the virtual domain of the coprime planar array are sufficiently utilized and the multi-source tensor spatial spectrum estimation with enhanced degree of freedom is realized.

[0080] The above are only the preferred embodiments of the present disclosure. Although the present disclosure has been disclosed as above in preferred embodiments, it is not intended to limit the present disclosure. Anyone skilled in the art, without departing from the scope of the technical solution of the present disclosure, can use the methods and technical content disclosed above to make many possible changes and modifications to the technical solution of the present disclosure, or modify it into equivalent changes. Therefore, all simple variations, equivalent changes and modifications made to the above embodiments based on the technical essence of the present disclosure without departing from the content of the technical solution of the present disclosure still fall within the protection scope of the technical solution of the present disclosure.