COMPOSITE TENSOR BEAMFORMING METHOD FOR ELECTROMAGNETIC VECTOR COPRIME PLANAR ARRAY

20230048116 · 2023-02-16

Assignee

Inventors

Cpc classification

International classification

Abstract

The present invention belongs to the field of array signal processing and relates to a composite tensor beamforming method for an electromagnetic vector coprime planar array. The method includes: building an electromagnetic vector coprime planar array; performing tensor modeling of an electromagnetic vector coprime planar array receiving signal; designing a three-dimensional weight tensor corresponding to a coprime sparse uniform sub-planar array; forming a tensor beam power pattern of the coprime sparse uniform sub-planar array; and performing electromagnetic vector coprime planar array tensor beamforming based on coprime composite processing of the sparse uniform sub-planar array. Starting from the principles of receiving signal tensor spatial filtering of two sparse uniform sub-planar arrays that compose the electromagnetic vector coprime planar array, the present invention forms a coprime composite processing method based on a sparse uniform sub-planar array output signal.

Claims

1. A composite tensor beamforming method for an electromagnetic vector coprime planar array, comprising: step 1: building the electromagnetic vector coprime planar array; step 2: performing tensor modeling of an electromagnetic vector coprime planar array receiving signal; step 3: designing a three-dimensional weight tensor corresponding to a coprime sparse uniform sub-planar array; step 4: forming a tensor beam power pattern of the coprime sparse uniform sub-planar array; and step 5: performing electromagnetic vector coprime planar array tensor beamforming based on coprime composite processing of the sparse uniform sub-planar array.

2. The composite tensor beamforming method for the electromagnetic vector coprime planar array as claimed in claim 1, wherein the step 1 comprises: structuring a pair of sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2 on a plane coordinate system xoy of a receiving end, custom-character.sub.1 and custom-character.sub.2 respectively comprising custom-character×custom-character and custom-character×custom-character antenna elements, wherein custom-character, custom-character and custom-character, custom-character are respectively a pair of coprime integers; intervals of antenna elements of the sparse uniform sub-planar array custom-character.sub.1 in x axis and y axis directions being respectively custom-characterd and custom-characterd, wherein a unit interval is d=λ/2, and λ denotes a signal wavelength; similarly, intervals of the antenna elements of the sparse uniform sub-planar array custom-character.sub.2 in the x axis and the y axis directions being respectively custom-characterd and custom-characterd, wherein in custom-character.sub.1, positions of the (custom-character).sup.th antenna element in the x axis and the y axis directions are respectively x 1 ( m 1 ) = ( m 1 - 1 ) M 2 d and y 1 ( n 1 ) = ( n 1 - 1 ) N 2 d , wherein custom-character=1, 2, . . . , custom-character=1, 2, . . . , custom-character; in custom-character.sub.2; in custom-character.sub.2, positions of the (custom-character).sup.th antenna element in the x axis and the y axis directions are respectively x 2 ( m 2 ) = ( m 2 - 1 ) M 1 d and y 2 ( n 2 ) = ( n 2 - 1 ) N 1 d , wherein custom-character=1, 2, . . . , custom-character=1, 2, . . . , custom-character; and combining sub-arrays in a manner of superimposing antenna elements custom-character=custom-character=custom-character=custom-character=0) at a position of an origin point of the coordinate system for custom-character.sub.1 and custom-character.sub.2, thereby obtaining an electromagnetic vector coprime planar array actually comprising custom-character+custom-character−1 antenna elements, wherein each of the antenna elements uses three mutually orthogonal electric doublets and three mutually orthogonal magnetic dipoles to realize sensing of electromagnetic field, thereby possessing a six-path output.

3. The composite tensor beamforming method for the electromagnetic vector coprime planar array as claimed in claim 2, wherein the step 2 comprises: setting a far-field narrow-band desired signal that is incident to the electromagnetic vector coprime planar array from a (θ, φ) direction, wherein θ and φ respectively denote an azimuth angle and a pitch angle of the desired signal and θϵ[−π/2, π/2], φϵ[−π, π]; the six-path output of each of the elements in the electromagnetic vector coprime planar array simultaneously comprises Direction of Arrival (DOA) information U(θ, φ)∈custom-character.sup.6×2 and polarized state information g(γ, η)∈custom-character.sup.2, wherein γ∈[0, 2π] and η∈[−π, π] respectively denote a polarized auxiliary angle and a polarized phase difference, and a DOA matrix U(θ, φ) and a polarized state vector g(γ, η) are defined as: U ( θ , φ ) = [ - sin θ cos φcosθ cos θ cos φsinθ 0 - sin φ cos φcosθ sin θ cos φsinθ - cos θ - sin φ 0 ] , g ( γ , η ) = [ cos γ sin γ e j η ] , wherein j=√{square root over (−1)}, and correspondingly, output of each of the elements in the electromagnetic vector coprime planar array is denoted with a spatial electromagnetic response vector pϵcustom-character.sup.6 as follows:
p=U(θ, φ)g(γ, η); when G non-relevant interfering signals exist simultaneously in a space, the DOA matrix, the polarized state vector and the spatial electromagnetic response vector thereof are respectively denoted by Ū(θ.sub.g, φ.sub.g), g(γ.sub.g, η.sub.g) and p.sub.g=Ū(θ.sub.g, φ.sub.g)g(γ.sub.g,η.sub.g), wherein g=1, 2, . . . , G; reserving three-dimensional spatial information of a receiving signal of the sparse uniform sub-planar array custom-character.sub.i (i=1, 2) at time t, i.e. DOA information and spatial electromagnetic response information in the x axis direction and the y axis direction, which are denoted with one three-dimensional tensor, and superimposing a three-dimensional signal tensor snapped by the collected T sampling blocks on a time dimension as a fourth dimension, thereby constituting a receiving signal tensor χ i M i × N i × 6 × T corresponding to the sparse uniform sub-planar array custom-character.sub.i, the receiving signal tensor χ i M i × N i × 6 × T being denoted as follows: χ i = a x ( i ) a y ( i ) p s + .Math. g = 1 G a _ xg ( i ) a _ yg ( i ) p _ g s _ g + i , wherein a x ( i ) = [ 1 , e - j 2 π λ x i ( 2 ) μ , .Math. , e - j 2 π λ x i ( M i ) μ ] T M i and a y ( i ) = [ 1 , e - j 2 π λ x i ( 2 ) v , .Math. , e - j 2 π λ x i ( N i ) v ] T N i respectively denote a desired signal guiding vector of the electromagnetic vector coprime planar array in the x axis and the y axis directions, and μ=sin φ cos θ and v=sin φ sin θ, s=[s(1), s(2), . . . , s(T)].sup.T∈custom-character.sup.T is a signal waveform of the desired signal, ∘ denotes an outer product of vectors, (⋅).sup.T denotes an transposition operation, and i M P i × N i × 6 × T is an independent co-distributed additive white Gaussian noise tensor; and then a ¯ x g ( i ) M i and a ¯ y g ( i ) N i respectively denote guiding vectors of the electromagnetic vector coprime planar array in the x axis and the y axis directions, corresponding to the g.sup.th interfering signal, and s.sub.g∈custom-character.sup.T denotes a signal waveform of the g.sup.th interfering signal.

4. The composite tensor beamforming method for the electromagnetic vector coprime planar array as claimed in claim 3, wherein the step 3 comprises: for a receiving signal tensor χ i ( t ) M i × N i × 6 of two sparse uniform sub-planar arrays that compose the electromagnetic vector coprime planar array at the time t, setting a three-dimensional weight tensor 𝒲 i M i × N i × 6 matching multi-dimensional structure information thereof, performing spatial filtering on custom-character(t) through custom-character, and forming a beam directivity in the DOA corresponding to the desired signal, thereby obtaining an output signal custom-character(t), which is denoted as follows:
custom-character(t)=<custom-character(t), custom-character>, t=1, 2, . . . , T, wherein <⋅>denotes an inner product of tensors, (⋅)* denotes a conjugation operation; then minimizing an average output power of a tensor beamformer and performing optimization processing such that the DOA of the desired signal and a response corresponding to a polarized state thereof should not be distorted, thereby obtaining a tensor beamformer custom-character corresponding to two sparse uniform sub-planar arrays, the optimization processing expression being as follows: min 𝒲 i E [ .Math. "\[LeftBracketingBar]" y i ( t ) .Math. "\[RightBracketingBar]" 2 ] s . t . < 𝒲 1 * , i ( θ , φ , γ , η ) > = 1 , wherein i ( θ , φ , γ , η ) = a x ( i ) a y ( i ) p M i × N i × 6 denotes a three-dimensional space manifold tensor of the sparse uniform sub-planar array custom-character.sub.i corresponding to a DOA (θ, φ) and a polarized state (γ, η) of a desired signal, |⋅| denotes a modulo operation of complex number, and E[⋅] denotes an expectation-taking operation; through solving, three-dimensional weight tensors custom-character and custom-character respectively corresponding to sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2 are obtained and output signals custom-character(t) and custom-character(t) are generated; wherein each space dimension information of the three-dimensional weight tensors custom-character and custom-character(t) corresponds to each other, custom-character decomposed in a manner of CANDECOMP/PARAFAC is denoted as an outer product of a beamforming weight vector corresponding to DOA information w 1 ( i ) M i in the x axis, DOA information w 2 ( i ) N i in the y axis and spatial electromagnetic response information custom-charactercustom-character:
custom-character=custom-character then, an output signal custom-character(t) of the sparse uniform sub-planar array custom-character.sub.i, at the time t can be denoted as follows: y i ( t ) = χ i ( t ) × 1 w 1 ( i ) * × 2 w 2 ( i ) * × 3 w 3 ( i ) * , wherein ×.sub.r denotes an inner product of a tensor and a matrix along the r.sup.th dimension; a weight tensor custom-character corresponding to a receiving signal tensor custom-character(t) is weighted to be equivalently denoted as multi-dimensional weight of the above three beamforming weight vector custom-character, r=1, 2, 3, for custom-character(t), and a corresponding optimization problem can be denoted as follows: min w r ( i ) ( r = 1 , 2 , 3 ) w r ( i ) H R r ( i ) w r ( i ) s . t . w 1 ( i ) H a x ( i ) = 1 , w 2 ( i ) H a y ( i ) = 1 , w 3 ( i ) H p = 1 , wherein R r ( i ) = E [ x r ( i ) ( t ) ( x r ( i ) ( t ) ) H ] , x r ( i ) denotes an output signal of the sparse uniform sub-planar array custom-character.sub.i at the r.sup.th dimension, and a beamforming weight vector of remaining two dimensions other than the r.sup.th dimension is obtained after custom-character(t) is weighted, and is denoted as follows: x 1 ( i ) ( t ) = χ i ( t ) × 2 w 2 ( i ) * × 3 w 3 ( i ) * M i , x 2 ( i ) ( t ) = χ i ( t ) × 1 w 1 ( i ) * × 3 w 3 ( i ) * N i , x 3 ( i ) ( t ) = χ i ( t ) × 1 w 1 ( i ) * × 2 w 2 ( i ) * 6 , wherein (⋅).sup.H denotes conjugation and transposition operations, and Lagrangian multiplier method is used to solve in order six sub-optimization problems corresponding to sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2, and their respective three beamforming weight vectors custom-character (r=1, 2, 3) and custom-character (r=1, 2, 3), with closed-form solutions thereof as follows: w 1 ( 1 ) = R 1 ( 1 ) - 1 a x ( 1 ) a x ( 1 ) H R 1 ( 1 ) - 1 a x ( 1 ) , w 1 ( 2 ) = R 1 ( 2 ) - 1 a x ( 2 ) a x ( 2 ) H R 1 ( 2 ) - 1 a x ( 2 ) , w 2 ( 1 ) = R 2 ( 1 ) - 1 a y ( 1 ) a y ( 1 ) H R 2 ( 1 ) - 1 a y ( 1 ) , w 2 ( 2 ) = R 2 ( 2 ) - 1 a y ( 2 ) a y ( 2 ) H R 2 ( 2 ) - 1 a y ( 2 ) , w 3 ( 1 ) = R 3 ( 1 ) - 1 p p H R 3 ( 1 ) - 1 p , w 3 ( 2 ) = R 3 ( 2 ) - 1 p p H R 3 ( 2 ) - 1 p .

5. The composite tensor beamforming method for the electromagnetic vector coprime planar array as claimed in claim 4, wherein the step 4 comprises: denoting the tensor beam power pattern custom-character({acute over (θ)}, {acute over (φ)}) of the sparse uniform sub-planar array tensor beamformer custom-character equivalently as follows through a CANDECOMP/PARAFAC decomposition form substituted into custom-character: 𝒫 i ( θ , φ ) = .Math. "\[LeftBracketingBar]" < 𝒲 1 * , 1 ( θ , φ , γ , η ) > .Math. "\[RightBracketingBar]" 2 = .Math. "\[LeftBracketingBar]" ( a x ( i ) ( θ , φ ) w 1 ( i ) H ) ( a y ( i ) ( θ , φ ) w 2 ( i ) H ) ( p ( θ , φ , γ , η ) w 3 ( i ) H ) .Math. "\[RightBracketingBar]" 2 , wherein {acute over (θ)}ϵ[−π/2, π/2] and {acute over (φ)}ϵ[−π, π]; when DOA is in a direction of a desired signal, i.e. {acute over (θ)}=θ and {acute over (φ)}=φ, a tensor beam power value of custom-character({acute over (θ)}, {acute over (φ)}) reaches a maximum, which is regarded as a main lobe; at a two-dimensional DOA plane, virtual peaks exist in both tensor beam power patterns custom-character({acute over (θ)}, {acute over (φ)}) and custom-character({acute over (θ)}, {acute over (φ)}) of sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2 and virtual peak positions (custom-character) and (custom-character) respectively corresponding thereto do not overlap each other, i.e. custom-charactercustom-charactercustom-character.

6. The composite tensor beamforming method for the electromagnetic vector coprime planar array as claimed in claim 5, wherein the step 5 comprises: performing coprime composite processing on output signals of two sparse uniform sub-planar arrays, the virtual peak positions of which do not overlap each other, thereby realizing virtual-peak restrained electromagnetic vector coprime planar array tensor beamforming, wherein the coprime composite processing comprises coprime composite processing based on multiplicative rules and coprime composite processing based on power minimization rules.

7. The composite tensor beamforming method for the electromagnetic vector coprime planar array as claimed in claim 6, wherein principles of the coprime composite processing based on multiplicative rules are as follows: when, in a two-dimensional DOA (custom-character) a tensor beam power pattern custom-character({acute over (θ)}, {acute over (φ)}) of custom-character.sub.1 corresponds to a virtual peak, and a tensor beam power pattern custom-character({acute over (θ)}, {acute over (φ)}) of custom-character.sub.2 does not correspond to a virtual peak, thus at a position of (custom-character), tensor beam power of custom-character({acute over (θ)}, {acute over (φ)}) and custom-character({acute over (θ)}, {acute over (φ)}) is multiplied and the virtual peak is retrained; similarly, when, in a two-dimensional DOA (custom-character), a tensor beam power pattern custom-character({acute over (θ)}, {acute over (φ)}) of custom-character.sub.2 corresponds to a virtual peak, and a tensor beam power pattern custom-character({acute over (θ)}, {acute over (φ)}) of custom-character.sub.1 does not correspond to a virtual peak, tensor beam power of custom-character({acute over (θ)}, {acute over (φ)}) and custom-character({acute over (θ)}, {acute over (φ)}) is multiplied and the virtual peak corresponding to the position can also be restrained; and an electromagnetic vector coprime planar array output signal y.sub.mul(t) based on multiplicative rules is obtained by multiplying output signals custom-character(t) and custom-character(t) of sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2 at the time t and is denoted as follows:
y.sub.mul(t)=custom-character(t)*custom-character(t), correspondingly, the tensor beam power pattern of the electromagnetic vector coprime planar array is an arithmetic square root of a product of tensor beam power patterns of two sparse uniform sub-planar arrays: 𝒫 mul ( θ , φ ) = 𝒫 1 ( θ , φ ) * 𝒫 2 ( θ , φ ) 2 = .Math. "\[LeftBracketingBar]" < 𝒲 1 * , 1 ( θ , φ , γ , η ) > * < 𝒲 2 * , 2 ( θ , φ , γ , η ) > .Math. "\[RightBracketingBar]" .

8. The composite tensor beamforming method for the electromagnetic vector coprime planar array as claimed in claim 6, wherein principles of the coprime composite processing based on power minimization rules are as follows: in a two-dimensional DOA (custom-character), a virtual peak response value custom-character(custom-character), of custom-character({acute over (θ)}, {acute over (φ)}) is greater than a response value custom-character(custom-character) corresponding to a non-virtual peak position of custom-character({acute over (θ)}, {acute over (φ)}) and the virtual peak is restrained by selecting a minimum value thereof; similarly, on (custom-character), a virtual peak response value custom-character(custom-character) of custom-character({acute over (θ)}, {acute over (φ)}) is greater than a non-virtual peak position response value custom-character(custom-character) of custom-character({acute over (θ)}, {acute over (φ)}) and the virtual peak is also restrained by selecting a minimum value thereof; and an output signal y.sub.min(t) of the electromagnetic vector coprime planar array based on power minimization rules is obtained by conducting minimization processing on power of output signals custom-character(t) and custom-character(t) of sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2 at the time t:
y.sub.min(t)=min(|custom-character(t)|.sup.2, |custom-character(t)|.sup.2), wherein min(⋅) denotes a minimum value taking operation; and correspondingly, the tensor beam power pattern of the electromagnetic vector coprime planar array is constituted by selecting a minimum value through comparison of tensor beam power of two sparse uniform sub-planar arrays in each two-dimensional DOA:
custom-character({acute over (θ)}, {acute over (φ)})=min(|<custom-character({acute over (θ)}, {acute over (φ)}, γ, η)>|.sup.2).

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0040] FIG. 1 is a block diagram for a total flow of the present invention.

[0041] FIG. 2 is a structural diagram of an electromagnetic vector coprime planar array in the present invention.

[0042] FIG. 3 is a block diagram of coprime composite processing flow based on multiplicative rules in the present invention.

[0043] FIG. 4 is a block diagram of coprime composite processing flow based on power minimization rules in the present invention.

[0044] FIG. 5a is a schematic diagram of a tensor beam power pattern effect based on multiplicative rules in the present invention.

[0045] FIG. 5b is a schematic diagram of a tensor beam power pattern effect based on power minimization rules in the present invention.

[0046] FIG. 6a is a performance comparison diagram of an output signal to noise ratio (SINR) varying with SNR in the present invention.

[0047] FIG. 6b is a performance comparison diagram of an output signal to noise ratio (SINR) varying with the number of snapshots of samples T.

DESCRIPTION OF THE EMBODIMENTS

[0048] To make the objectives, technical solutions, and technical effects of the present invention more comprehensible, the following describes the present invention in details with reference to accompanying drawings and embodiments.

[0049] As shown in FIG. 1, the present invention realizes electromagnetic vector coprime planar array tensor beamforming having capability of restraining a virtual peak and output performance improvement by spatial filtering of a coprime sparse uniform sub-planar array receiving signal tensor and matching coprime composite processing on a sub-planar array output signal of a feature that virtual peaks corresponding to coprime sparse uniform sub-planar arrays do not overlap each other, with specific realization steps comprising:

[0050] step 1: building an electromagnetic vector coprime planar array;

[0051] using custom-character+custom-character−1 electromagnetic vector antenna elements by a receiving end to construct an electromagnetic vector coprime planar array, wherein each of the antenna elements uses three mutually orthogonal electric doublets and three mutually orthogonal magnetic dipoles to realize sensing of electromagnetic field, thereby possessing a six-path output;

[0052] As shown in FIG. 2, a pair of sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2 are structured on a plane coordinate system xoy, custom-character.sub.1 and custom-character.sub.2 respectively comprising custom-character×custom-character and custom-character×custom-character antenna elements, wherein custom-character and custom-character are respectively a pair of coprime integers; intervals of antenna elements of the sparse uniform sub-planar array custom-character.sub.1 in x axis and y axis directions are respectively custom-characterd and custom-characterd, wherein a unit interval is d=λ/2, and λ denotes a signal wavelength; similarly, intervals of the antenna elements of the sparse uniform sub-planar array custom-character.sub.2 in x axis and y axis directions are respectively custom-characterd and custom-characterd, wherein in custom-character.sub.1, positions of the (custom-character).sup.th antenna element in x axis and y axis directions are respectively

[00018] x 1 ( m 1 ) = ( m 1 - 1 ) M 2 d and y 1 ( n 1 ) = ( n 1 - 1 ) N 2 d ,

wherein custom-character=1, 2, . . . , custom-character, custom-character=1, 2, . . . , custom-character; in custom-character.sub.2, positions of the (custom-character).sup.th antenna element in x axis and y axis directions are respectively

[00019] x 2 ( m 2 ) = ( m 2 - 1 ) M 1 d and y 2 ( n 2 ) = ( n 2 - 1 ) N 1 d ,

wherein custom-character=1, 2, . . . , custom-character, custom-character=1, 2, . . . , custom-character; and sub-arrays are combined in a manner of superimposing) elements (custom-character=custom-character=custom-character=custom-character=0) at a position of an origin point of the coordinate system for custom-character.sub.1 and custom-character.sub.2, thereby obtaining an electromagnetic vector coprime planar array actually comprising custom-character+custom-character−1 antenna elements;

[0053] step 2: performing tensor modeling of an electromagnetic vector coprime planar array receiving signal;

[0054] setting a far-field narrow-band desired signal that is incident to the electromagnetic vector coprime planar array from a (θ, φ) direction, wherein θ and φ respectively denote an azimuth angle and a pitch angle of the desired signal and θϵ[−π/2, π/2], φϵ[−π, π]; the six-path output of each of the elements in the electromagnetic vector coprime planar array simultaneously comprises Direction of Arrival (DOA) information U(θ, φ)∈custom-character.sup.6×2 and polarized state information g(γ, η)∈custom-character.sup.2, wherein γ∈[0, 2π] and η∈[−π, π] respectively denote a polarized auxiliary angle and a polarized phase difference, and a DOA matrix U(θ, φ) and a polarized state vector g(γ, η) are specifically defined as:

[00020] U ( θ , φ ) = [ - s in θ cos φ cos θ cos θ cos φ sin θ 0 - s in φ cos φ cos θ sin θ cos φ sin θ - cos θ - s in φ 0 ] , g ( γ , η ) = [ cos γ sin γ e j η ] ,

[0055] wherein j=√{square root over (−1)}, and correspondingly, output of each of the elements in the electromagnetic vector coprime planar array is denoted with a spatial electromagnetic response vector pϵcustom-character.sup.6 as follows:


p=U(θ, φ)g(γ, η).

[0056] when G non-relevant interfering signals exist simultaneously in a space, the DOA matrix, the polarized state vector and the spatial electromagnetic response vector thereof are respectively denoted by Ū(θ.sub.g, 100.sub.g), g(γ.sub.g, η.sub.g) and p.sub.g=Ū(θ.sub.g, φ.sub.g)g(γ.sub.g, η.sub.g), wherein g=1, 2, . . . , G;

[0057] reserving three-dimensional spatial information of a receiving signal of the sparse uniform sub-planar array custom-character.sub.i (i=1, 2) at time t, i.e. DOA information and spatial electromagnetic response information in x axis direction and y axis direction, which are denoted with one three-dimensional tensor, and superimposing a three-dimensional signal tensor snapped by the collected T sampling blocks on a time dimension as a fourth dimension, thereby constituting a receiving signal tensor

[00021] 𝓍 i M i × N i × 6 × T

corresponding to the sparse uniform sub-planar array custom-character.sub.i, the receiving signal tensor

[00022] 𝓍 i M i × N i × 6 × T

being denoted as follows:

[00023] 𝓍 i = a x ( i ) a y ( i ) 𝓅 𝓈 + .Math. g = 1 G a ¯ x g ( i ) a ¯ y g ( i ) 𝓅 _ g 𝓈 _ g + i , ( 2 ) wherein a x ( i ) = [ 1 , e - j 2 π λ x i ( 2 ) μ , .Math. , e - j 2 π λ x i ( M i ) μ ] T M i and a y ( i ) = [ 1 , e - j 2 π λ y i ( 2 ) υ , .Math. , e - j 2 π λ y i ( N i ) υ ] T N i

respectively denote a desired signal guiding vector of the electromagnetic vector coprime planar array in x axis and y axis directions, and μ=sin φ cos θ and v=sin φ sin θ, s=[s(1), s(2), . . . , s(T)].sup.T∈custom-character.sup.T is a signal waveform of the desired signal, ∘ denotes an outer product of vectors, (⋅).sup.T denotes an transposition operation, and

[00024] i M i × N i × 6 × T

is an independent co-distributed additive white Gaussian noise tensor; and then

[00025] a ¯ x g ( i ) M i and a ¯ y g ( i ) N i

respectively denote guiding vectors of the electromagnetic vector coprime planar array in x axis and y axis directions, corresponding to the g.sup.th interfering signal, and s.sub.g∈custom-character.sup.T denotes a signal waveform of the g.sup.th interfering signal.

[0058] step 3: designing a three-dimensional weight tensor corresponding to a coprime sparse uniform sub-planar array;

[0059] for a receiving signal tensor

[00026] 𝓍 i ( t ) M i × N i × 6

of two sparse uniform sub-planar arrays that compose the electromagnetic vector coprime planar array at time t, setting a three-dimensional weight tensor

[00027] 𝓌 i M i × N i × 6

matching multi-dimensional structure information thereof, performing spatial filtering on custom-character(t) through custom-character, and forming a beam directivity in the DOA corresponding to the desired signal, thereby obtaining an output signal custom-character(t), which is denoted as follows:


custom-character(t)=<custom-character(t), custom-character>, t=1, 2, . . . , T,

[0060] wherein <⋅> denotes an inner product of tensors, (⋅)* denotes a conjugation operation; then in order to obtain a tensor beamformer custom-character corresponding to two sparse uniform sub-planar arrays, minimizing an average output power of a tensor beamformer and performing optimization processing to ensure that the DOA of the desired signal and a response corresponding to a polarized state thereof should not be distorted, with an expression being as follows:

[00028] min 𝓌 i E [ .Math. "\[LeftBracketingBar]" y i ( t ) .Math. "\[RightBracketingBar]" 2 ] s . t . < 𝓌 i * , i ( θ , φ , γ , η ) >= 1 , wherein i ( θ , φ , γ , η ) = a x ( i ) a y ( i ) p M i × N i × 6

denotes a three-dimensional space manifold tensor of the sparse uniform sub-planar array custom-character.sub.i corresponding to a DOA (θ, φ) and a polarized state (γ, η) of a desired signal, |⋅| denotes a modulo operation of complex number, and E[⋅] denotes an expectation-taking operation; through solving, three-dimensional weight tensors custom-character and custom-character respectively corresponding to sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2 are obtained and output signals custom-character(t) and custom-character(t) are generated;

[0061] wherein each space dimension information of the three-dimensional weight tensors custom-character and custom-character(t) corresponds to each other, custom-character decomposed in a manner of CANDECOMP/PARAFAC is denoted as an outer product of a beamforming weight vector corresponding to DOA information in

[00029] 𝓌 1 ( i ) M i

in x axis, DOA information

[00030] 𝓌 2 ( i ) M i

in y axis and spatial electromagnetic response information custom-character:


custom-character=custom-character.

[0062] then, an output signal custom-character(t) of the sparse uniform sub-planar array custom-character.sub.i, at time t can be denoted as follows:


custom-character(t)=custom-character(t)×.sub.1custom-character×.sub.2custom-character×.sub.3custom-character,

[0063] wherein ×.sub.r denotes an inner product of a tensor and a matrix along the r.sup.th dimension; thus, a weight tensor custom-character corresponding to a receiving signal tensor custom-character(t) is weighted to be equivalently denoted as multi-dimensional weight of the above three beamforming weight vector custom-character, r=1, 2, 3, for custom-character(t), and a corresponding optimization problem can be denoted as follows:

[00031] min 𝓌 r ( i ) ( r = 1 , 2 , 3 ) 𝓌 r ( i ) H R r ( i ) 𝓌 r ( i ) s . t . 𝓌 1 ( i ) H a x ( i ) = 1 , 𝓌 2 ( i ) H a y ( i ) = 1 , w 3 ( i ) H p = 1 , wherein R r ( i ) = E [ 𝓍 r ( i ) ( t ) ( x r ( i ) ( t ) ) H ] , 𝓍 r ( i )

denotes an output signal of the sparse uniform sub-planar array custom-character.sub.i at the r.sup.th dimension, and a beamforming weight vector of remaining two dimensions other than the r.sup.th dimension is obtained after custom-character(t) is weighted, and is denoted as follows:

[00032] x 1 ( i ) ( t ) = 𝓍 i ( t ) × 2 𝓌 2 ( i ) * × 3 𝓌 3 ( i ) * M i , x 2 ( i ) ( t ) = 𝓍 i ( t ) × 1 𝓌 1 ( i ) * × 3 𝓌 3 ( i ) * N i , x 3 ( i ) ( t ) = 𝓍 i ( t ) × 1 𝓌 1 ( i ) * × 2 𝓌 2 ( i ) * 6 ,

[0064] wherein (⋅).sup.H denotes conjugation and transposition operations, and Lagrangian multiplier method is used to solve in order six sub-optimization problems corresponding to sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2, and their respective three beamforming weight vectors custom-character (r=1, 2, 3) and custom-character (r=1, 2, 3), with closed-form solutions thereof as follows:

[00033] 𝓌 1 ( 1 ) = R 1 ( 1 ) - 1 a x ( 1 ) a x ( 1 ) H R 1 ( 1 ) - 1 a x ( 1 ) , 𝓌 1 ( 2 ) = R 1 ( 2 ) - 1 a x ( 2 ) a x ( 2 ) H R 1 ( 2 ) - 1 a x ( 2 ) , 𝓌 2 ( 1 ) = R 2 ( 1 ) - 1 a y ( 1 ) a y ( 1 ) H R 2 ( 1 ) - 1 a y ( 1 ) , 𝓌 2 ( 2 ) = R 2 ( 2 ) - 1 a y ( 2 ) a y ( 2 ) H R 2 ( 2 ) - 1 a y ( 2 ) , 𝓌 3 ( 1 ) = R 3 ( 1 ) - 1 𝓅 𝓅 H R 3 ( 1 ) - 1 𝓅 , 𝓌 3 ( 2 ) = R 3 ( 2 ) - 1 𝓅 𝓅 H R 3 ( 2 ) - 1 𝓅

[0065] step 4: forming a tensor beam power pattern of the coprime sparse uniform sub-planar array;

[0066] denoting the tensor beam power pattern custom-character({acute over (θ)}, {acute over (φ)}) of the sparse uniform sub-planar array tensor beamformer custom-character equivalently as follows through a CANDECOMP/PARAFAC decomposition form substituted into custom-character:

[00034] 𝒫 i ( θ , φ ) = .Math. "\[LeftBracketingBar]" < 𝓌 i * , i ( θ , φ , γ , η ) > .Math. "\[RightBracketingBar]" 2 = .Math. "\[LeftBracketingBar]" ( a x ( i ) ( θ , φ ) 𝓌 1 ( i ) H ) ( a y ( i ) ( θ , φ ) 𝓌 2 ( i ) H ) ( 𝓅 ( θ , φ , γ , η ) 𝓌 3 ( i ) H ) .Math. "\[RightBracketingBar]" 2 ,

[0067] wherein {acute over (θ)}ϵ[−π/2, π/2] and {acute over (φ)}ϵ[−π, π]; when DOA is in a direction of a desired signal, i.e. {acute over (θ)}=θ and {acute over (φ)}=φ, a tensor beam power value of custom-character({acute over (θ)}, {acute over (φ)}) reaches a maximum, which is regarded as a main lobe. However, an interval of elements in the sparse uniform sub-planar array is greater than a half wavelength, thereby failing to satisfy a Nyquist sampling rate, resulting in that when ({acute over (θ)}, {acute over (φ)})=(custom-character), α=1, 2, . . . , A a virtual peak exists in custom-character({acute over (θ)}, {acute over (φ)}); and when ({acute over (θ)}, {acute over (φ)})=(custom-character), b=1, 2, . . . , B, a virtual peak exists in custom-character({acute over (θ)}, {acute over (φ)}). Since the arrangement of the elements of the sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2 along x axis direction and y axis direction satisfies a coprime feature, thus in a two-dimensional DOA plane, the virtual peak positions (custom-character) and (custom-character) respectively corresponding to the sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2 do not overlap each other, i.e. custom-charactercustom-charactercustom-character.

[0068] step 5: performing electromagnetic vector coprime planar array tensor beamforming based on coprime composite processing of the sparse uniform sub-planar array.

[0069] performing composite processing on output signals of coprime sparse uniform sub-planar arrays by using a feature that virtual peak positions of the two sparse uniform sub-planar arrays do not overlap each other, thereby realizing virtual-peak restrained electromagnetic vector coprime planar array tensor beamforming,

[0070] wherein the coprime composite processing of the sparse uniform sub-planar array output signal comprises coprime composite processing based on multiplicative rules and coprime composite processing based on power minimization rules.

[0071] Principles of the coprime composite processing based on multiplicative rules are as follows: when, in a two-dimensional DOA (custom-character), a tensor beam power pattern custom-character({acute over (θ)}, {acute over (φ)}) of custom-character.sub.1 corresponds to a virtual peak, and a tensor beam power pattern custom-character({acute over (θ)}, {acute over (φ)}) of custom-character.sub.2 does not correspond to a virtual peak, thus at a position of (custom-character), tensor beam power of custom-character({acute over (θ)}, {acute over (φ)}) and custom-character({acute over (θ)}, {acute over (φ)}) is multiplied and the virtual peak is retrained; similarly, when, in a two-dimensional DOA (custom-character), a tensor beam power pattern custom-character({acute over (θ)}, {acute over (φ)}) of custom-character.sub.2 corresponds to a virtual peak, and a tensor beam power pattern custom-character({acute over (θ)}, {acute over (φ)}) of custom-character.sub.1 does not correspond to a virtual peak, tensor beam power of custom-character({acute over (θ)}, {acute over (φ)}) and custom-character({acute over (θ)}, {acute over (φ)}) is multiplied and the virtual peak corresponding to the position can also be restrained. As shown in FIG. 3, an electromagnetic vector coprime planar array output signal based on multiplicative rules is obtained by multiplying output signals and of sparse uniform sub-planar arrays and at time and is denoted as follows: As shown in FIG. 3, an electromagnetic vector coprime planar array output signal y.sub.mul(t) based on multiplicative rules is obtained by multiplying output signals custom-character(t) and custom-character(t) of sparse uniform sub-planar array custom-character.sub.1 and custom-character.sub.2 at time t and is denoted as follows:


y.sub.mul(t)=custom-character(t)*custom-character(t).

[0072] correspondingly, the tensor beam power pattern thereof is an arithmetic square root of a product of tensor beam power patterns of two sparse uniform sub-planar arrays:

[00035] 𝒫 mul ( θ , φ ) = 𝒫 1 ( θ , φ ) * 𝒫 2 ( θ , φ ) 2 = .Math. "\[LeftBracketingBar]" < 𝒲 1 * , 1 ( θ , φ , γ , η ) > * < 𝒲 2 * , 2 ( θ , φ , γ , η ) > .Math. "\[RightBracketingBar]" .

[0073] Principles of the coprime composite processing based on power minimization rules are as follows: in a two-dimensional DOA (custom-character), a virtual peak response value custom-character(custom-character) of custom-character({acute over (θ)}, {acute over (φ)}) is greater than a response value custom-character(custom-character) corresponding to a non-virtual peak position of custom-character({acute over (θ)}, {acute over (φ)}) and the virtual peak is restrained by selecting a minimum value thereof; similarly, on (custom-character), a virtual peak response value custom-character(custom-character) of custom-character({acute over (θ)}, {acute over (φ)}) is greater than a non-virtual peak position response value custom-character(custom-character) of custom-character({acute over (θ)}, {acute over (φ)}) and the virtual peak is also restrained by selecting a minimum value thereof. As shown in FIG. 4, an output signal under the rules is obtained by conducting minimization processing on power of output signals custom-character(t) and custom-character(t) of sparse uniform sub-planar arrays custom-character.sub.1 and custom-character.sub.2 at time t.


y.sub.min(t)=min(|custom-character(t)|.sup.2, |custom-character(t)|.sup.2),

[0074] wherein min (⋅) denotes a minimum value taking operation; and correspondingly, the tensor beam power pattern thereof is constituted by selecting a minimum value through comparison of tensor beam power of two sparse uniform sub-planar arrays in each two-dimensional DOA:


custom-character.sub.min({acute over (θ)}, {acute over (φ)})=min(|<custom-character({acute over (θ)}, {acute over (φ)}, γ, η)>|.sup.2).

[0075] The effects of the present invention are further described below in combination with embodiments:

[0076] Embodiment 1: an electromagnetic vector coprime planar array is used to receive an incident signal and parameters thereof are selected as custom-character=custom-character=5 and custom-character=custom-character=4, that is, the electromagnetic vector coprime planar array of the architecture comprises custom-character+custom-character−1=40 antenna elements in total. It is assumed that a desired signal is located at (θ, φ)=(30°, 45°) and carries a polarized auxiliary angle γ=15° and a phase difference subangle η=−20°; an interfering signal is located at (θ.sub.1, φ.sub.1)=(10°, 15°); and γ.sub.1=−25°, η.sub.1=−35°.

[0077] When a signal-to-noise ratio (SNR) of a desired signal is 0 dB and the number of snapshots of the samples is T=300 , the tensor beam power patterns custom-character.sub.mul({acute over (θ)}, {acute over (φ)}) and custom-character.sub.min({acute over (θ)}, {acute over (φ)}) of the electromagnetic vector coprime planar array based on multiplicative rules and power minimization rules are drawn as shown in FIG. 5a and FIG. 5b. The tensor beam power pattern of the electromagnetic vector coprime planar array corresponds to one main lobe at a DOA position of a desired signal and a virtual peak does not exist at other positions. Thus, it can be seen that the composite tensor beamforming method of the electromagnetic vector coprime planar array as proposed restrains the virtual peak effectively.

[0078] Embodiment 2: further, the composite tensor beamforming method of the electromagnetic vector coprime planar array as proposed is compared with a signal-to-interference-plus-noise ratio (SINR) performance of the tensor signal processing method based on the electromagnetic vector uniform planar array. In order to ensure fairness of simulative comparison, 40 elements are arranged for the electromagnetic vector uniform planar array according to a structure with 5 rows and 8 columns. When the number of snapshots of the samples is T=300, a performance comparison curve of an output SINR varying with the SNR is drawn as shown in FIG. 6a. Under the condition of SNR=0 dB, a performance comparison curve of an output SINR varying with the number of snapshots of the samples T is drawn as shown in FIG. 6a. It can be seen from a comparison result of FIG. 6a and FIG. 6b that no matter in a SNR scene for different desired signals or in a different scene for the number of snapshots of samples T, the output SINR performance of the composite tensor beamforming method of the electromagnetic vector coprime planar array based on multiplicative rules and power minimization rules as proposed is superior to the tensor signal processing method based on the electromagnetic vector uniform planar array. Benefited from a big-hole-diameter advantage brought about by sparse arrangement of elements in the electromagnetic vector coprime planar array and effective function of the method in restraining the virtual peak, the electromagnetic vector coprime planar array has a higher output SINR as compared with a uniform planar array. Meanwhile, since power minimization rules restrict response of a virtual peak to a maximum extent on a tensor beam power pattern and the electromagnetic coprime planar array tensor beamforming corresponding thereto is superior to the electromagnetic vector coprime planar array tensor beamforming based on multiplicative rules.

[0079] To sum up, the present invention matches structural space information covered in the multi-dimensional receiving signal of the electromagnetic vector coprime planar array, thereby forming principles of spatial filtering of a coprime sparse uniform sub-planar array receiving signal tensor. In addition, the present invention matches the coprime arrangement feature of the two sparse uniform sub-planar arrays to perform coprime composite processing on the output signal of the sparse uniform sub-planar array by using a feature that virtual peaks do not overlap each other in the tensor beam power patterns of the two sparse uniform sub-planar arrays, so as to realize electromagnetic vector coprime planar array tensor beamforming having a capability of restraining the virtual peak and improvement of output performance.

[0080] The above are merely preferred embodiments of the present invention. Although the preferred embodiments of the present invention are disclosed above, they are not used to restrict the present invention. Any person skilled in the art who is familiar with the field can make many likely changes and modifications to the technical solution of the present invention with the method and technical contents disclosed above or modify them as equivalent embodiments of equivalent change without departing from the scope of the technical solution of the present invention. Therefore, any content without departing from the technical solution of the present invention, any simple change, equivalent change and modification made to the above embodiments according to the technical substance of the present invention, belong to the protection scope of the technical solution of the present invention.