CN109901146B - Hemispherical array and differential angle decoupling method based on spiral distribution - Google Patents

Hemispherical array and differential angle decoupling method based on spiral distribution Download PDF

Info

Publication number
CN109901146B
CN109901146B CN201910199805.6A CN201910199805A CN109901146B CN 109901146 B CN109901146 B CN 109901146B CN 201910199805 A CN201910199805 A CN 201910199805A CN 109901146 B CN109901146 B CN 109901146B
Authority
CN
China
Prior art keywords
azimuth
substep
angle
sum
array
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910199805.6A
Other languages
Chinese (zh)
Other versions
CN109901146A (en
Inventor
朱圣棋
周季峰
曾操
许京伟
刘婧
廖桂生
李军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201910199805.6A priority Critical patent/CN109901146B/en
Publication of CN109901146A publication Critical patent/CN109901146A/en
Application granted granted Critical
Publication of CN109901146B publication Critical patent/CN109901146B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a semispherical array and differential angle decoupling method based on spiral distribution, which comprises the following steps: firstly, distributing N array elements on a hemispherical surface with the radius of R according to a spiral distribution arrangement mode; then redefining the azimuth angle of the signal coming from a certain direction
Figure DDA0001996998090000011
And a pitch angle theta to obtain a redefined slow vector r; then, the azimuth deviation is obtained
Figure DDA0001996998090000012
And four beams with the pitching direction deviation of +/-Delta theta, respectively obtaining a sum beam, an azimuth direction difference beam and a pitching direction difference beam for the four beams, and then obtaining an azimuth direction difference sum ratio and a pitching direction difference sum ratio; then setting an azimuth observation window
Figure DDA0001996998090000013
And a pitch viewing window theta win So as to obtain the angle identifying curves of the uncoupled azimuth direction and the uncoupled elevation direction in the observation window. The method reduces the pre-stored data quantity, has smaller angle measurement error per se and high angle measurement speed, and ensures the requirements of speed and precision of angle measurement of the radar array.

Description

Hemispherical array and differential angle decoupling method based on spiral distribution
Technical Field
The invention belongs to the technical field of radar signal processing, and particularly relates to a hemispherical array and differential angle decoupling method, which is suitable for a hemispherical radar array in spiral distribution, and solves the problem that the data volume required to be prestored is huge due to the azimuth direction and elevation direction coupling existing in the process of measuring the sum and differential angles.
Background
Early radar array, main wired array and area array etc. because the array distribution is more regular, antenna element's directional unanimity, consequently it is easier to carry out theoretical analysis, and the shortcoming is that the receiving visual angle of array is limited, and hemispherical distributed's radar visual angle has advantages such as directive property is good, interference clutter suppression ability is strong, detection distance is far away and wide beam visual angle, but hemispherical array's antenna element directional nonconformity leads to hemispherical array theoretical analysis more difficult.
The sum and difference angle measurement technology has an early origin and obvious advantages, can measure the angle information of a target in a pulse, and has great application in engineering practice along with the research of a radar algorithm and the great improvement of the precision of the sum and difference angle measurement, particularly the sum and difference angle measurement technology of a linear array, which is developed to date and has a double-beam pointing method, a direct weighting method, a symmetrical negation method and the like. The main problem of area array sum and difference angle measurement exists in the information coupling of azimuth and elevation, the main solving method comprises two methods of coordinate system transformation, guide vector transformation and the like, when the method is applied to a hemispherical array, because the directions of all antenna units are inconsistent in the hemispherical array based on directional array elements, the theoretical analysis is difficult; meanwhile, azimuth information and elevation information are coupled, and a decoupling method of an area array is not completely applicable, so that a sum-difference angle measurement technology based on a hemispherical array has difficulty.
Disclosure of Invention
In order to solve the problems, the invention discloses a sum and difference angle measurement decoupling method based on a hemispherical array, which redefines the azimuth angle and the pitch angle of a coordinate system on the basis of the hemispherical array in spiral distribution, respectively synthesizes four beams, respectively obtains the difference and the ratio of the pitch direction and the pitch direction of an expected direction, approximately considers that the difference and the ratio of the azimuth direction are irrelevant to the pitch angle and the difference and the ratio of the pitch direction are irrelevant to the azimuth angle within a certain small angle range of the expected direction, and thus achieves the purpose of decoupling the azimuth direction and the pitch direction.
In order to achieve the purpose, the invention adopts the following technical scheme to realize the purpose:
a sum and difference angle decoupling method based on a hemispherical array comprises the following steps:
step 1: distributing N array elements on a hemispherical surface with the radius of R according to a spiral distribution arrangement mode;
step 2: redefining the azimuth angle of a signal arriving from a certain direction
Figure BDA0001996998070000021
And a pitch angle theta to obtain a redefined slow vector r;
and 3, step 3: assuming a target beam with a desired direction of
Figure BDA0001996998070000022
On the basis of step 1 and step 2, the respective synthetic beam directions are respectively
Figure BDA0001996998070000023
And
Figure BDA0001996998070000024
four beams, denoted P respectively 1 、P 2 、P 3 And P 4 (ii) a Wherein Δ θ and
Figure BDA0001996998070000025
angular deviations in azimuth and pitch directions, respectively;
and 4, step 4: separately calculating sum beam P sum Elevation difference beam P ele Azimuth difference beam P azi Differential pitch and ratio r ele Sum of azimuth differences r azi
And 5: set up the azimuth observation window
Figure BDA0001996998070000026
And a pitch viewing window theta win Respectively taking an azimuth viewInspection window
Figure BDA0001996998070000027
And a pitch viewing window theta win Inner pitch sum ratio r ele Sum of azimuth differences r azi And angle identifying curves of the azimuth direction and the pitching direction are obtained, so that the decoupling purpose is achieved.
Compared with the prior art, the invention has the following advantages:
(1) The method effectively solves the problem of coupling of the azimuth direction and the pitch direction, thereby effectively reducing the data volume of the angle identification curve required to be stored by the sum and difference angle measurement algorithm, reducing the engineering implementation complexity of the angle measurement algorithm and reducing the data volume;
(2) The method obtains the angle identification curve on the basis of approximation, reduces the data volume, does not increase errors obviously, and improves the comprehensive performance of the algorithm.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a schematic flow chart of a hemispherical array and differential angle decoupling method based on spiral distribution;
FIG. 2 is a graph of the array element coordinate of the mth array element in the hemispherical array;
FIG. 3 is a slow vector definition diagram in a rectangular coordinate system;
FIG. 4 is a spiral distribution diagram of a hemispherical array for simulation experiments;
FIG. 5 is a beam pattern for the pointing (0, 0) direction of the simulation experiment;
FIG. 6 shows a beam P of a simulation experiment 1 、P 2 、P 3 、P 4 A drawing;
FIG. 7 is a view in the X-Y dimension of FIG. 6;
FIG. 8 is an enlarged view of a portion of FIG. 7;
FIG. 9 is a sum beam diagram of a simulation experiment;
FIG. 10 is an azimuth difference beam pattern of a simulation experiment;
FIG. 11 is a view of elevation difference beams from a simulation experiment;
FIG. 12 is a plot of the azimuthal difference sum ratio within the observation window for a simulation experiment;
FIG. 13 is an X-Z view of the azimuthal difference sum ratio within the observation window of a simulation experiment;
FIG. 14 is a graph of the fitting results of the azimuthal curves in the observation window of the simulation experiment;
FIG. 15 is a plot of the sum of the pitching differences within the observation window of a simulation experiment;
FIG. 16 is a Y-Z dimensional view of the sum of the pitching differences within the observation window of the simulation experiment;
FIG. 17 is a graph of the results of fitting the azimuthal curves in the observation window for the simulation experiment;
FIG. 18 is a plot of SNR versus azimuthal estimation error for a simulation experiment;
FIG. 19 is a plot of SNR versus pitch estimation error for a simulation experiment;
FIG. 20 is a slow vector definition plot in the conventional array coordinate system for comparative experiments;
FIG. 21 is an original view of four beam patterns in a conventional coordinate system of a comparative experiment;
FIG. 22 is an X-Y view of FIG. 21;
FIG. 23 is an original view of the azimuthal difference and ratio within the observation window of the comparative experiment;
FIG. 24 is an X-Z view of FIG. 23;
FIG. 25 is an original view of the sum of the pitching differences within the observation window of a comparative experiment;
fig. 26 is a Y-Z view of fig. 25.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Fig. 1 is a schematic flowchart of a hemispherical array based on spiral distribution and a differential angle decoupling method according to an embodiment of the present invention.
As shown in fig. 1, the hemispherical array based on spiral distribution and the differential angle decoupling method provided in the embodiment of the present invention include the following steps:
step 1: distributing N array elements on a hemispherical surface with the radius of R according to a spiral distribution arrangement mode;
n array elements are distributed on the hemispherical spherical surface, and the whole hemisphere is collectively called a radar array, namely the radar array comprises the N array elements.
Sub-step 1a, n array elements, where the position of the m-th array element in the coordinate system as shown in fig. 2 is denoted as point p m Remember q m Is a point p m A projection point on a plane XOY, wherein
Figure BDA0001996998070000051
∠ZOp m =θ m Then, according to the definition of the helical distribution:
θ m =arccos(h m ),
Figure BDA0001996998070000052
wherein m is more than or equal to 1 and less than or equal to N;
Figure BDA0001996998070000053
wherein m is more than or equal to 2 and less than or equal to N-1,
Figure BDA0001996998070000054
substep 1b, as can be obtained from substep 1a, the XYZ coordinate of the m-th array element in the rectangular coordinate system as shown in the figure is:
Figure BDA0001996998070000055
Figure BDA0001996998070000056
z m =R sinθ m
and a substep 1c, according to the array element coordinates of the m-th array element obtained in the substep 1b, distributing the N array elements on a hemisphere with a radius of R according to a spiral.
Step 2: redefining the azimuth angle of a signal arriving in a certain direction
Figure BDA0001996998070000057
And a pitch angle θ, thereby obtaining a redefined slow vector r; the method specifically comprises the following substeps:
substep 2a, assume an incoming signal as in FIG. 3
Figure BDA0001996998070000058
Q is a projection point of a p point on a plane XOZ, black solid dots in the figure are array elements, and the angle qOZ = theta is defined,
Figure BDA0001996998070000059
and a substep 2b, obtaining a slow vector as follows according to the definition of the slow vector:
Figure BDA0001996998070000061
and step 3: assuming a target beam with a desired direction of
Figure BDA0001996998070000062
On the basis of step 1 and step 2, the synthesized beam is respectively pointed to
Figure BDA0001996998070000063
And
Figure BDA0001996998070000064
four beams, respectively denoted as P 1 、P 2 、P 3 And P 4 (ii) a Wherein, Δ θ and
Figure BDA0001996998070000065
angular deviations in azimuth and pitch directions, respectively;
the method specifically comprises the following substeps:
substep 3a of obtaining the pointing direction according to the definition of hemispherical array beam forming
Figure BDA0001996998070000066
The beam of (a) is:
Figure BDA0001996998070000067
wherein the content of the first and second substances,
Figure BDA0001996998070000068
and theta is as defined in figure 3,
Figure BDA0001996998070000069
angle of target azimuth, θ 0 Is the angle of a target pitch angle, N is the total number of array elements in the radar array, H represents the conjugate transpose,
Figure BDA00019969980700000610
denotes Heander Meng Deji, L i Is the spatial position of the ith array element in the radar array, and L i =[x i y i z i ] T I is more than or equal to 1 and less than or equal to N, r is a slow vector,
Figure BDA00019969980700000611
for spatial orientation of the ith array element
Figure BDA00019969980700000612
Weighted value of, and
Figure BDA00019969980700000613
j is an imaginary factor, λ is the wavelength of the incoming signal,
Figure BDA00019969980700000614
is the ith array element directional diagram function, and
Figure BDA00019969980700000615
p can be obtained according to the above definition 1 、P 2 、P 3 And P 4 Respectively as follows:
Figure BDA00019969980700000616
Figure BDA00019969980700000617
Figure BDA00019969980700000618
Figure BDA00019969980700000619
and 4, step 4: computing sum beam P sum Elevation difference beam P ele Azimuth difference beam P azi Differential pitch and ratio r ele Sum-azimuth ratio r azi
The method specifically comprises the following substeps:
substep 4a of obtaining a sum beam P based on the definition of the sum beam sum The calculation formula of (c) is:
P sum =P 1 +P 2 +P 3 +P 4
substep 4b, obtaining the elevation difference beam P according to the definition of the elevation difference beam ele The calculation formula of (2) is as follows:
P ele =P 1 +P 2 -P 3 -P 4
substep 4c, according to the definition of the azimuth difference beam, obtaining the azimuth difference beam P azi The calculation formula of (2) is as follows:
P azi =P 1 -P 2 +P 3 -P 4
substep 4d, obtaining the pitching difference sum ratio r according to the definition of the pitching difference sum ratio ele The calculation formula of (c) is:
Figure BDA0001996998070000071
substep 4e, obtaining the azimuth difference sum ratio r according to the definition of the azimuth difference sum ratio azi The calculation formula of (2) is as follows:
Figure BDA0001996998070000072
thereby obtaining a sum beam P sum Elevation difference beam P ele Azimuth difference beam P azi Differential pitch and ratio r ele Sum of azimuth differences r azi
And 5: set up the direction observation window
Figure BDA0001996998070000073
And a pitch viewing window theta win Respectively taking the azimuth viewing window
Figure BDA0001996998070000074
And a pitch viewing window theta win Inner pitch sum ratio r ele Sum-azimuth ratio r azi And angle identifying curves of the azimuth direction and the pitching direction are obtained, so that the decoupling purpose is achieved.
Specifically, the azimuth observation window is arranged at 1-2 degrees
Figure BDA0001996998070000075
And a pitch viewing window theta win . The sum of the differences in pitch direction r obtained in the step 4 ele Sum of azimuth differences r azi In fact one about θ and
Figure BDA0001996998070000081
by setting a small range of azimuthal viewing windows
Figure BDA0001996998070000082
And a pitch viewing window theta win Approximately, the sum of the pitching differences r is considered ele And azimuth angle
Figure BDA0001996998070000083
Independently, the sum of the azimuthal differences r azi Independent of the pitch angle theta, thereby achieving the purpose of decoupling. Then the angle curve of the azimuth direction is the azimuth direction difference and the ratio curve in theta o The elevation angle curve of the section is the elevation difference and the specific surface
Figure BDA0001996998070000084
The slice of (6).
The effect of the method provided by the embodiment of the invention is verified through a simulation experiment as follows:
1. simulation experiment environment and data:
the experimental environment is as follows: inter (R) Core (TM) i5-6500 CPU@3.20HGz,64 bits Windows operating system and MATLAB 2016b simulation software.
Experimental parameters: the total number N of array elements of the radar array is 200, and the azimuth angle
Figure BDA0001996998070000085
The value range of the pitch angle theta is-90 degrees, and the target direction is
Figure BDA0001996998070000086
Is (0 deg. ), a pitching declination angle delta theta =1 deg., an azimuth declination angle
Figure BDA0001996998070000087
Elevation direction observation window theta win =1 °, squarePosition observation window
Figure BDA0001996998070000088
2. Simulation experiment results:
the spiral distribution of the hemispherical array obtained from step 1 is shown in fig. 4;
the beam pointing to the target azimuth (0 ° ) obtained according to step 3 is shown in fig. 5;
beam P 1 、P 2 、P 3 And P 4 As shown in fig. 6, fig. 7 is an X-Y dimensional view of fig. 6, and fig. 8 is a partially enlarged view of fig. 7. As can be seen in fig. 8, the four beams are directed at (-1 ° ), (-1 °,1 °), (1 °, -1 °) and (1 ° ), respectively.
As can be seen in fig. 9, the sum beam is directed (0 ° ). As can be seen from fig. 10 and 11, the azimuth difference beam and the elevation difference beam have a peak and a trough respectively beside the (0 ° ) point, except that the peak and the trough beside the azimuth difference beam (0 ° ) are
Figure BDA0001996998070000089
In this dimension, the peaks and troughs next to the elevation difference beam (0 ° ) are in the dimension θ.
From fig. 12 and 13, it can be seen that the azimuthal difference sum ratio in the viewing window is only related to the azimuth angle
Figure BDA0001996998070000091
In this regard, regardless of the pitch angle θ, the present invention recognizes that the azimuthal difference and ratio approximate the azimuth angle within the observation window
Figure BDA0001996998070000092
Is a linear relation, the curve obtained by fitting the linear relation is shown in figure 14, as can be seen from figure 14, the fitting effect is good, and the expression of the fitting straight line is
Figure BDA0001996998070000093
From fig. 15 and 16, it can be seen that the sum of the difference in pitch in the observation windowRatio is related only to pitch angle θ, to azimuth angle
Figure BDA0001996998070000094
Regardless, in the observation window, the invention considers that the pitch difference sum ratio approximation and the pitch angle theta are in a linear relation, the curve obtained by fitting the linear relation is shown in figure 17, as can be seen from figure 17, the fitting effect is good, and the expression of the fitting straight line is r ele =-0.1175θ。
On the basis of the angle profile shown in fig. 14 and 17, the sum and difference angle measurement and error analysis are performed to obtain an error-SNR curve as shown in fig. 18 and 19.
As can be seen from fig. 18 and 19, on the basis of the angle identification curves obtained in fig. 14 and 17, the sum and difference angles are measured by using the array with hemispherical spiral distribution, the obtained error is about 0.07, the error is in a relatively small range, and the influence of the SNR on the error is not large, which is mainly caused by the accumulation of signals before the sum beam and the difference beam are ratioed.
As shown in fig. 18 and 19, the size of the azimuth observation window and the pitch observation window are set to be 2 ° respectively, because the azimuth observation window and the pitch observation window are already decoupled, and therefore the pre-stored straight line data is 180 straight lines; however, if there is no decoupling, the amount of data to be prestored is about thousands of planes, and each plane is assumed to be composed of several hundreds of straight lines, and the straight lines to be stored is about 10 5 And (3) strips.
3. And (3) comparing the experimental results:
the conventional array coordinate system is defined as shown in FIG. 20, where p is the space-electromagnetic wave vector, q is the projection of p on XOY plane, where
Figure BDA0001996998070000095
Angle ZOp = θ. Then the slow vector is:
Figure BDA0001996998070000096
assuming that the spherical radius R =10 λ and the angle of pointing of the center beam is (10 ° ), beams pointing in four directions of (11 °,9 °), (9 °,11 °), (9 ° ) and (11 °,11 °) are synthesized, respectively, as shown in fig. 21 and 22:
FIG. 21 is an original view of four beam patterns in a conventional coordinate system, and FIG. 22 is an X-Y view of FIG. 21; it can be seen that: only two beams in elevation are shown. This is because the beams are relatively wide in azimuth but have small differences in pointing angles, so that the two beams in azimuth are coincident, and therefore only the two beams in elevation can be seen in the figure. In summary, it can be concluded that: in a conventional coordinate system, four beam patterns distributed spirally are completely asymmetric in elevation.
Assumed azimuth observation window
Figure BDA0001996998070000101
Elevation direction observation window theta win The calculated azimuth difference sum ratio and the calculated pitch difference sum ratio in the observation window are shown in fig. 23 to 26. FIG. 23 is an original view of the azimuthal difference sum ratio within the viewing window, FIG. 24 is an X-Z view of FIG. 23, FIG. 25 is an original view of the pitch difference sum ratio within the viewing window, and FIG. 26 is a Y-Z view of FIG. 25. As can be seen from the figure, the azimuth difference sum ratio and the pitch difference sum ratio are not only related to the azimuth angle, but also related to the pitch angle, so that under the conventional coordinate system, the azimuth direction and the pitch direction of the spirally distributed hemispherical array are coupled, and the angle is measured, the pre-stored data will be much.
As can be seen from the results of the simulation experiment and the comparison experiment, the method effectively solves the problem of the coupling of the azimuth direction and the pitch direction existing in the process of the spiral distribution of the hemispherical array and the angle difference measurement in a partial range on the basis of the redefinition and approximation of coordinates, greatly reduces the data quantity required to be prestored in the angle difference measurement, and ensures that the error of the angle difference measurement is in a smaller range.
Those of ordinary skill in the art will understand that: all or part of the steps for realizing the method embodiments can be completed by hardware related to program instructions, the program can be stored in a computer readable storage medium, and the program executes the steps comprising the method embodiments when executed; and the aforementioned storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
The above description is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art can easily conceive of the changes or substitutions within the technical scope of the present invention, and all the changes or substitutions should be covered within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the appended claims.

Claims (6)

1. A semispherical array and differential angle decoupling method based on spiral distribution is characterized by comprising the following steps:
step 1: distributing N array elements on a hemispherical surface with the radius of R according to a spiral distribution arrangement mode;
step 2: redefining the azimuth angle of a signal arriving from a certain direction
Figure FDA0003869622910000011
And a pitch angle theta to obtain a redefined slow vector r;
and step 3: assuming a target beam with a desired direction of
Figure FDA0003869622910000012
On the basis of step 1 and step 2, the synthesized beam is respectively pointed to
Figure FDA0003869622910000013
And
Figure FDA0003869622910000014
four beams, respectively denoted as P 1 、P 2 、P 3 And P 4 (ii) a Wherein, Δ θ and
Figure FDA0003869622910000015
angular deviations in azimuth and pitch directions, respectively;
and 4, step 4: computing sum beam P sum Elevation difference beam P ele Azimuth difference beam P azi Differential pitch and ratio r ele Sum of azimuth differences r azi
And 5: set up the direction observation window
Figure FDA0003869622910000016
And a pitch viewing window theta win Respectively taking the azimuth direction observation window
Figure FDA0003869622910000017
And a pitch viewing window theta win Inner pitch sum ratio r ele Sum of azimuth differences r azi And angle identifying curves of the azimuth direction and the pitching direction are obtained, so that the decoupling purpose is achieved.
2. The hemispherical array and differential angle decoupling method based on spiral distribution as claimed in claim 1, wherein in step 1, the N array elements are distributed on a hemispherical surface with a radius R according to the arrangement of spiral distribution, which includes the following sub-steps:
sub-step 1a, N array elements, where the position of the m-th array element in the coordinate system is denoted as point p m Remember q m Is a point p m A projection point on a plane XOY, wherein
Figure FDA0003869622910000018
Then, according to the definition of the helical distribution, we get:
Figure FDA0003869622910000019
wherein m is more than or equal to 1 and less than or equal to N;
Figure FDA0003869622910000021
wherein m is more than or equal to 2 and less than or equal to N-1,
Figure FDA0003869622910000022
Substep 1b, according to substep 1a, the XYZ coordinate of the m-th array element in the rectangular coordinate system is:
Figure FDA0003869622910000023
Figure FDA0003869622910000024
z m =Rsinθ m
and a substep 1c, according to the array element coordinates of the m-th array element obtained in the substep 1b, distributing the N array elements on a hemisphere with a radius of R according to a spiral.
3. The hemispherical array and differential angle decoupling method based on spiral distribution as claimed in claim 1, wherein in step 2, the step of obtaining the redefined slow vector r comprises the following sub-steps:
substep 2a, assuming an incoming signal
Figure FDA0003869622910000025
q is the projection point of the point p on the plane XOZ, defines ≈ qOZ = theta,
Figure FDA0003869622910000026
and a substep 2b, obtaining a slow vector as follows according to the definition of the slow vector:
Figure FDA0003869622910000027
4. according to claim 1The hemispherical array based on spiral distribution and the differential angle decoupling method are characterized in that in step 3, four beams P 1 、P 2 、P 3 And P 4 The calculation algorithm of (2) comprises the following substeps:
substep 3a, deriving the pointing direction according to the definition of hemispherical array beamforming
Figure FDA0003869622910000028
The beam of (a) is:
Figure FDA0003869622910000029
wherein the content of the first and second substances,
Figure FDA00038696229100000210
angle of target azimuth, θ 0 Is the angle of a target pitch angle, N is the total number of array elements in the radar array, H represents conjugate transposition,
Figure FDA00038696229100000211
denotes Heander Meng Deji, L i Is the spatial position of the ith array element in the radar array, and L i =[x i y i z i ] T I is more than or equal to 1 and less than or equal to N, r is a slow vector,
Figure FDA0003869622910000031
for spatial orientation of the ith array element
Figure FDA0003869622910000032
Weighted value of, and
Figure FDA0003869622910000033
j is an imaginary factor, λ is the wavelength of the incoming signal,
Figure FDA0003869622910000034
is the ith array element directionA graph function, and
Figure FDA0003869622910000035
substep 3b, according to substep 3a, obtaining the orientation
Figure FDA0003869622910000036
Directional beam P 1 Comprises the following steps:
Figure FDA0003869622910000037
substep 3c, according to substep 3a, obtaining the orientation
Figure FDA0003869622910000038
Directional beam P 2 Comprises the following steps:
Figure FDA0003869622910000039
substep 3d, according to substep 3a, obtaining the orientation
Figure FDA00038696229100000310
Directional beam P 3 Comprises the following steps:
Figure FDA00038696229100000311
substep 3e, according to substep 3a, obtaining the orientation
Figure FDA00038696229100000312
Directional beam P 4 Comprises the following steps:
Figure FDA00038696229100000313
5. the method according to claim 1, wherein the calculating in step 4 comprises the following sub-steps:
substep 4a of obtaining a sum beam P based on the definition of the sum beam sum The calculation formula of (2) is as follows:
P sum =P 1 +P 2 +P 3 +P 4
wherein, P 1 Is directed to
Figure FDA00038696229100000314
A directional beam; p 2 Is directed to
Figure FDA00038696229100000315
A directional beam; p 3 Is directed to
Figure FDA00038696229100000316
A directional beam; p 4 Is directed to
Figure FDA00038696229100000317
A directional beam;
substep 4b, obtaining the elevation difference beam P according to the definition of the elevation difference beam ele The calculation formula of (2) is as follows:
P ele =P 1 +P 2 -P 3 -P 4
substep 4c, according to the definition of the azimuth difference beam, obtaining the azimuth difference beam P azi The calculation formula of (2) is as follows:
P azi =P 1 -P 2 +P 3 -P 4
substep 4d, obtaining the pitching difference sum ratio r according to the definition of the pitching difference sum ratio ele The calculation formula of (2) is as follows:
Figure FDA0003869622910000041
a substep 4e of obtaining a sum of azimuth differences r according to the definition of the sum of azimuth differences r azi The calculation formula of (2) is as follows:
Figure FDA0003869622910000042
6. the method for decoupling hemispherical array and differential angle based on spiral distribution according to claim 1, wherein in step 5, the obtaining of the angle identifying curves of azimuth and elevation specifically comprises the following sub-steps:
substep 5a, setting an azimuth observation window of 1-2 DEG
Figure FDA0003869622910000043
And a pitch viewing window theta win Viewing in the azimuth viewing window
Figure FDA0003869622910000044
And a pitch viewing window theta win Azimuthal difference within range and ratio r azi And the sum of the pitching differences r ele
Substep 5b, consider the sum of the pitching differences r ele And azimuth angle
Figure FDA0003869622910000045
Independently, the sum of the azimuthal differences r azi Independent of the pitch angle theta, obtaining azimuth angle curve and pitch angle curve respectively, namely the azimuth angle curve is azimuth difference and ratio curve surface at theta o The elevation angle curve of the section is the elevation difference and the specific surface
Figure FDA0003869622910000046
The slice of (6).
CN201910199805.6A 2019-03-15 2019-03-15 Hemispherical array and differential angle decoupling method based on spiral distribution Active CN109901146B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910199805.6A CN109901146B (en) 2019-03-15 2019-03-15 Hemispherical array and differential angle decoupling method based on spiral distribution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910199805.6A CN109901146B (en) 2019-03-15 2019-03-15 Hemispherical array and differential angle decoupling method based on spiral distribution

Publications (2)

Publication Number Publication Date
CN109901146A CN109901146A (en) 2019-06-18
CN109901146B true CN109901146B (en) 2023-03-14

Family

ID=66952666

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910199805.6A Active CN109901146B (en) 2019-03-15 2019-03-15 Hemispherical array and differential angle decoupling method based on spiral distribution

Country Status (1)

Country Link
CN (1) CN109901146B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB208087A (en) * 1923-05-28 1923-12-13 Francis Joseph Stawell Jones Improvements in and connected with internal combustion engine fuel tanks having automatically controlled feeds
US5130714A (en) * 1991-05-23 1992-07-14 Hughes Aircraft Company Stretch and chirp waveform format for reduced generating and receiving hardware complexity
CN202091337U (en) * 2011-04-21 2011-12-28 江苏南方轴承股份有限公司 Decoupling appliance core shaft capable of preventing force transmission interference and decoupling appliance of same
WO2012099704A2 (en) * 2011-01-17 2012-07-26 Novita Therapeutics, Llc Blockstent device and methods of use
CN103035994A (en) * 2011-08-22 2013-04-10 英飞凌科技股份有限公司 Microstrip coupler combining transmit-receive signal separation and differential to single ended conversion
CN103235292A (en) * 2013-05-08 2013-08-07 西安电子科技大学 Full-dimension and difference angle measurement method for zero setting conformal calibration of a planar phased array

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB208087A (en) * 1923-05-28 1923-12-13 Francis Joseph Stawell Jones Improvements in and connected with internal combustion engine fuel tanks having automatically controlled feeds
US5130714A (en) * 1991-05-23 1992-07-14 Hughes Aircraft Company Stretch and chirp waveform format for reduced generating and receiving hardware complexity
WO2012099704A2 (en) * 2011-01-17 2012-07-26 Novita Therapeutics, Llc Blockstent device and methods of use
CN202091337U (en) * 2011-04-21 2011-12-28 江苏南方轴承股份有限公司 Decoupling appliance core shaft capable of preventing force transmission interference and decoupling appliance of same
CN103035994A (en) * 2011-08-22 2013-04-10 英飞凌科技股份有限公司 Microstrip coupler combining transmit-receive signal separation and differential to single ended conversion
CN103235292A (en) * 2013-05-08 2013-08-07 西安电子科技大学 Full-dimension and difference angle measurement method for zero setting conformal calibration of a planar phased array

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Measuring Leaf Angle Distribution Using Terrestrial Laser Scanning in a European Beech Forest;Jing Liu et al.;《IGARSS 2018》;20180727;第8977-8980页 *
矢量推进解耦球面并联机构动力学研究;张荣敏 等;《农业机械学报》;20150616;第46卷(第6期);第319-325页 *

Also Published As

Publication number Publication date
CN109901146A (en) 2019-06-18

Similar Documents

Publication Publication Date Title
CN108549059B (en) Low-altitude target elevation angle estimation method under complex terrain condition
CN105445718B (en) A kind of DOA estimation method of the distributed overloading warship over-the-horizon radar based on array reconfiguration
CN108535698B (en) Meter-wave radar low elevation angle estimation method based on beam space
CN107991659B (en) Method for measuring height of low-elevation target of meter wave radar based on dictionary learning
CN105306123A (en) Robust beamforming method with resistance to array system errors
CN109188342A (en) Low complex degree arrival direction estimation method under conformal circle battle array
CN107229032A (en) A kind of method and apparatus for building four array element solid arrays
CN103885045A (en) Sub-array division based circulation combined adaptive beam forming method
CN103513225B (en) Sparse planar formation optimization method based on spatial gain
CN113314832A (en) Millimeter wave vehicle-mounted MIMO radar antenna array device and design method
CN111398914B (en) Radar tracking beam angle measurement method
CN112612010A (en) Meter-wave radar low elevation height measurement method based on lobe splitting pretreatment
CN109633638B (en) Large-angle SAR super-resolution imaging method based on electromagnetic scattering model
CN109901146B (en) Hemispherical array and differential angle decoupling method based on spiral distribution
CN114779236A (en) Improved meter-wave radar low-elevation height measurement method based on spatial smoothing MUSIC
Frid et al. Determining direction-of-arrival accuracy for installed antennas by postprocessing of far-field data
CN106990385B (en) Non-homogeneous face battle array design method based on minimum-redundancy linear arrays
CN111812607B (en) Meter wave MIMO radar low elevation angle estimation method based on beam space
CN108051792B (en) Distribution description method of scattering source coupled between target and rough surface
CN113671485B (en) ADMM-based two-dimensional DOA estimation method for meter wave area array radar
CN110389319B (en) MIMO radar DOA estimation method based on low-altitude multipath condition
CN111983599A (en) Target two-dimensional DOA estimation method based on azimuth-pitch dictionary
CN107144813B (en) Method and device for constructing four-array-element three-dimensional array
CN113820654B (en) S-band radar target low elevation DOA estimation method based on beam domain dimension reduction
CN114488027A (en) Wave beam zero setting and widening method of conformal array antenna

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant