CN114879138B - Millimeter wave radar two-dimensional angle calculation method and device and related equipment - Google Patents

Millimeter wave radar two-dimensional angle calculation method and device and related equipment Download PDF

Info

Publication number
CN114879138B
CN114879138B CN202210818680.2A CN202210818680A CN114879138B CN 114879138 B CN114879138 B CN 114879138B CN 202210818680 A CN202210818680 A CN 202210818680A CN 114879138 B CN114879138 B CN 114879138B
Authority
CN
China
Prior art keywords
target
angle
array elements
virtual array
virtual
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
CN202210818680.2A
Other languages
Chinese (zh)
Other versions
CN114879138A (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.)
Guangdong Dawan District Aerospace Information Research Institute
Original Assignee
Guangdong Dawan District Aerospace Information Research Institute
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 Guangdong Dawan District Aerospace Information Research Institute filed Critical Guangdong Dawan District Aerospace Information Research Institute
Priority to CN202210818680.2A priority Critical patent/CN114879138B/en
Publication of CN114879138A publication Critical patent/CN114879138A/en
Application granted granted Critical
Publication of CN114879138B publication Critical patent/CN114879138B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14Systems for determining direction or deviation from predetermined direction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/93Radar or analogous systems specially adapted for specific applications for anti-collision purposes
    • G01S13/931Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/415Identification of targets based on measurements of movement associated with the target
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The application discloses a millimeter wave radar two-dimensional angle calculation method, a millimeter wave radar two-dimensional angle calculation device and related equipment, wherein the method comprises the following steps: selecting a plurality of virtual array elements from a target virtual array to form three groups of target array elements; determining an angle set of each target in the corresponding direction based on an echo data matrix of each virtual array element in each group of target array elements to obtain a first angle set of each target in the horizontal direction, a second angle set of each target in the vertical direction and a third angle set of each target in the oblique direction; and finally, determining the angle of each target based on a preset calculation formula of the echo signal, the first angle set, the second angle set and the third angle set. According to the method and the device, only part of virtual array elements are adopted for calculation, calculation of unnecessary information is reduced to a great extent, the calculated amount cannot be increased due to the increase of the number of the array elements in the target virtual array, the calculation process is simple, and two-dimensional angle estimation of the target angle can be completed efficiently.

Description

Millimeter wave radar two-dimensional angle calculation method and device and related equipment
Technical Field
The application relates to the technical field of millimeter wave radars, in particular to a millimeter wave radar two-dimensional angle calculation method, a millimeter wave radar two-dimensional angle calculation device and related equipment.
Background
The automobile industry is developing towards a more intelligent direction, and from an Advanced Driver Assistance System (ADAS) to a higher-level fully automatic driving, the implementation of these functions is not necessary to provide various sensors for sensing surrounding environment information for a decision-making System. The millimeter wave radar has the advantages of long detection distance, capability of measuring speed, capability of working all day long and all weather, low cost and the like, and becomes one of the most core sensors of the automobile.
The traditional millimeter wave radar sensor only has the angle measuring function of horizontal dimension, the requirement on the sensor is higher due to the improvement of the automatic driving level, and then the 4D millimeter wave radar is promoted. The target object which can be perceived by the 4D radar comprises information of four dimensions of distance, speed, horizontal angle and vertical angle. Compared with the traditional vehicle-mounted millimeter wave radar, the 4D radar uses a two-dimensional area array, has the capability of vertically measuring the angle, and has higher angle measurement resolution of two dimensions.
The problem to be solved by two-dimensional angle measurement is angle matching of the same target in the horizontal and vertical directions, and the traditional 2D-DOA (2-division Direction of Arrival) method is to perform 2D-FFT (2-division Fast Fourier Transform) or beam forming on all row elements and array elements of a planar array and then obtain a target angle by methods such as finding an extremum and the like so as to avoid the problem of angle matching. This method requires a large amount of computation, not only needs to perform FFT or beamforming on all rows and columns, but also finds extreme values in the whole matrix, and the amount of computation increases linearly as the number of array elements increases.
Disclosure of Invention
In view of this, the present application provides a millimeter wave radar two-dimensional angle calculation method, device and related equipment, so as to efficiently implement calculation of a target angle.
In order to achieve the above object, a first aspect of the present application provides a millimeter wave radar two-dimensional angle calculation method, including:
selecting a plurality of virtual array elements from a target virtual array to form three groups of target array elements, wherein the target virtual array is uniformly distributed in the horizontal direction and the vertical direction respectively, the first group of target array elements comprises at least two virtual array elements in the same horizontal direction, the second group of target array elements comprises at least two virtual array elements in the same vertical direction, and the third group of target array elements comprises at least two virtual array elements which are not in the same horizontal direction and not in the same vertical direction;
determining a first angle set of each target in the horizontal direction based on an echo data matrix of each virtual array element in the first group of target array elements, determining a second angle set of each target in the vertical direction based on an echo data matrix of each virtual array element in the second group of target array elements, and determining a third angle set of each target in the oblique direction based on an echo data matrix of each virtual array element in the third group of target array elements, wherein the first angle set, the second angle set and the third angle set contain equal number of elements and are the number of targets;
determining the angle of each target based on the calculation formula of the echo signal, the first angle set, the second angle set and the third angle set.
Preferably, the process of determining the angle of each target based on the calculation formula of the echo signal, the first angle set, the second angle set and the third angle set includes:
matching a first angle in the first angle set and a second angle in the second angle set according to the following equation to obtain a matched angle combination:
Figure 953431DEST_PATH_IMAGE001
determining each angle combination as an angle of a target, wherein each angle combination comprises a first angle and a second angle;
wherein,
Figure 393640DEST_PATH_IMAGE002
the distance of adjacent target array elements in the third set of target array elements in the horizontal direction,
Figure 56702DEST_PATH_IMAGE003
is the third set of targetsThe distance of the adjacent target array elements in the vertical direction,
Figure 520045DEST_PATH_IMAGE004
for a first angle in the first set of angles,
Figure 926755DEST_PATH_IMAGE005
for a second angle in the second set of angles,
Figure 487050DEST_PATH_IMAGE006
is a third angle in the third set of angles.
Preferably, the process of matching a first angle in the first angle set and a second angle in the second angle set according to the following equation includes:
construction matrix
Figure 993117DEST_PATH_IMAGE007
Each element in the matrix is calculated by the following equation:
Figure 336898DEST_PATH_IMAGE008
constructing a list of vectors
Figure 281721DEST_PATH_IMAGE009
The vector list is calculated by the following equation
Figure 696521DEST_PATH_IMAGE009
Each vector of (a):
Figure 373490DEST_PATH_IMAGE010
for each
Figure 873742DEST_PATH_IMAGE011
Determining the same as
Figure 622255DEST_PATH_IMAGE011
Target with the closest value of
Figure 891562DEST_PATH_IMAGE012
And target the object
Figure 5012DEST_PATH_IMAGE012
In (1)
Figure 992559DEST_PATH_IMAGE013
And
Figure 276255DEST_PATH_IMAGE014
determining the angle combinations matched with each other;
wherein,
Figure 665648DEST_PATH_IMAGE015
as to the number of the objects in question,
Figure 12315DEST_PATH_IMAGE013
for a first angle in the first set of angles,
Figure 221580DEST_PATH_IMAGE014
for a second angle in the second set of angles,
Figure 249579DEST_PATH_IMAGE016
is a third angle in the third set of angles;
Figure 493478DEST_PATH_IMAGE002
is prepared by reacting with
Figure 948730DEST_PATH_IMAGE013
Corresponding target array element, and
Figure 910870DEST_PATH_IMAGE014
the distance of the corresponding target array element in the horizontal direction;
Figure 539298DEST_PATH_IMAGE003
is prepared by reacting with
Figure 637704DEST_PATH_IMAGE013
Corresponding target array element, and
Figure 263857DEST_PATH_IMAGE014
the corresponding target array element, the distance in the vertical direction.
Preferably, the process of determining the first set of angles of the target in the horizontal direction based on the echo data matrix of each virtual array element in the first group of target array elements includes:
performing two-dimensional fast Fourier transform (2D-FFT) on each echo data matrix of each virtual array element in the first group of target array elements to obtain a distance-Doppler spectrum matrix;
performing non-correlation accumulation processing on the distance-Doppler spectrum matrix to obtain a distance-Doppler spectrum matrix after accumulation;
performing constant false alarm CFAR detection on the accumulated distance-Doppler spectrum matrix to obtain a distance and a speed index of a target;
a first set of angles of the target is determined based on the range-doppler spectrum matrix and the range and velocity indices of the target.
Preferably, the target virtual array is a rectangular array, and the process of selecting a plurality of virtual array elements from the target virtual array to form three groups of target array elements includes:
forming a first group of target array elements by all virtual array elements in any row in the target virtual array;
forming a second group of target array elements by each virtual array element in any column in the target virtual array;
and forming a third group of target array elements by each virtual array element on the diagonal line or the sub diagonal line in the target virtual array.
Preferably, the target virtual array is a non-rectangular array, and the process of selecting a plurality of virtual array elements from the target virtual array to form three groups of target array elements includes:
determining a row with the maximum number of virtual array elements from a target virtual array, and forming a first group of target array elements by all the virtual array elements in the row;
determining a column with the maximum number of virtual array elements from a target virtual array, and forming a second group of target array elements by all the virtual array elements in the column;
determining a plurality of virtual array elements which are not in the same row and are not in the same column from the target virtual array, wherein the centers of the plurality of virtual array elements are connected into a line, the intervals of the adjacent array elements in the plurality of virtual array elements in the horizontal direction are consistent, and the intervals of the adjacent array elements in the plurality of virtual array elements in the vertical direction are consistent.
Preferably, the distance between adjacent target array elements in the target virtual array in the horizontal direction
Figure 450643DEST_PATH_IMAGE002
Satisfying the following equation:
Figure 882762DEST_PATH_IMAGE017
wherein,
Figure 38937DEST_PATH_IMAGE018
being the wavelength of the millimeter wave used,
Figure 898308DEST_PATH_IMAGE019
the angle measurement range is a preset maximum non-fuzzy angle measurement range.
The second aspect of the present application provides a millimeter wave radar two-dimensional angle calculation apparatus, including:
the target array element selecting unit is used for selecting a plurality of virtual array elements from a target virtual array to form three groups of target array elements, wherein the target virtual array elements are uniformly distributed in the horizontal direction and the vertical direction respectively, the first group of target array elements comprises at least two virtual array elements in the same horizontal direction, the second group of target array elements comprises at least two virtual array elements in the same vertical direction, and the third group of target array elements comprises at least two virtual array elements which are not in the same horizontal direction and not in the same vertical direction;
a single angle determining unit, configured to determine a first angle set of each target in a horizontal direction based on an echo data matrix of each virtual array element in the first group of target array elements, determine a second angle set of each target in a vertical direction based on an echo data matrix of each virtual array element in the second group of target array elements, and determine a third angle set of each target in an oblique direction based on an echo data matrix of each virtual array element in the third group of target array elements, where the numbers of elements included in the first angle set, the second angle set, and the third angle set are equal and are the numbers of targets;
an angle matching calculation unit, configured to determine an angle of each target based on a calculation formula of an echo signal, the first angle set, the second angle set, and the third angle set.
The third aspect of the present application provides a millimeter wave radar two-dimensional angle calculation device, including: a memory and a processor;
the memory is used for storing programs;
the processor is configured to execute the program, and implement each step of the millimeter wave radar two-dimensional angle calculation method.
A fourth aspect of the present application provides a storage medium having a computer program stored thereon, where the computer program, when executed by a processor, implements the steps of the millimeter wave radar two-dimensional angle calculation method as described above.
According to the technical scheme, firstly, a plurality of virtual array elements are selected from the target virtual array to form three groups of target array elements. The target virtual arrays are uniformly distributed in the horizontal direction and the vertical direction respectively, that is, the intervals of the virtual array elements of the target virtual arrays in the same horizontal direction are consistent, and the intervals of the virtual array elements of the target virtual arrays in the same vertical direction are consistent. Then, based on the echo data matrix of each virtual array element in the first group of target array elements, determining a first angle set of the target in the horizontal direction; determining a second angle set of the target in the vertical direction based on the echo data matrix of each virtual array element in the second group of target array elements; and determining a third angle set of the target in the oblique direction based on the echo data matrix of each virtual array element in the third group of target array elements. The number of elements contained in the first angle set, the second angle set and the third angle set is equal to the number of the target elements. Because the first group of target array elements comprises at least two virtual array elements in the same horizontal direction, the second group of target array elements comprises at least two virtual array elements in the same vertical direction, and the third group of target array elements comprises at least two virtual array elements which are not in the same horizontal direction and not in the same vertical direction, therefore, each first angle in the first angle set embodies the angle information of the target in the horizontal direction, each second angle in the second angle set embodies the angle information of the target in the vertical direction, and each third angle in the third angle set embodies the associated information of the target in the horizontal direction and the vertical direction. And finally, determining the angle of each target based on a preset calculation formula of the echo signal, the first angle set, the second angle set and the third angle set. According to the method and the device, only part of virtual array elements are adopted for calculation, calculation of unnecessary information is reduced to a great extent, the calculated amount cannot be increased due to the increase of the number of the array elements in the target virtual array, the calculation process is simple, and the two-dimensional angle estimation of the target angle can be completed efficiently.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, it is obvious that the drawings in the following description are only embodiments of the present application, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts.
FIG. 1 illustrates a schematic diagram of a transmit receive array and a virtual array as disclosed in embodiments of the present application;
fig. 2 is a schematic diagram of a two-dimensional angle calculation method for a millimeter wave radar disclosed in the embodiment of the present application;
FIG. 3 illustrates a selection method for each group of target array elements disclosed in the embodiments of the present application;
FIG. 4 illustrates another selection method for each group of target array elements disclosed in the embodiments of the present application;
FIG. 5 illustrates various groups of target array elements in a non-rectangular array as disclosed in embodiments of the present application;
fig. 6 is a schematic diagram of a two-dimensional angle calculation apparatus for a millimeter wave radar disclosed in an embodiment of the present application;
fig. 7 is a schematic diagram of a millimeter wave radar two-dimensional angle calculation device disclosed in the embodiment of the present application.
Detailed Description
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, 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 application.
For convenience of understanding, the following first describes an antenna array used in the millimeter wave radar two-dimensional angle calculation method provided in the present application. Referring to the left part of fig. 1, the antenna array (actual transceiving array) provided in the embodiment of the present application adopts a uniform area array design, which includes 5 transmitting antennas (Tx 1-Tx 5) uniformly arranged in the vertical direction and 5 receiving antennas (Rx 1-Rx 5) uniformly arranged in the horizontal direction, and is applied to a millimeter wave radar system in a Multiple-Input Multiple-Output (MIMO) mode. The multiple-sending and multiple-receiving method is characterized in that the antenna panel of the radar uses a plurality of transmitting antennas and a plurality of receiving antennas which are arranged according to a certain design to complete the transmission and the reception of electromagnetic wave signals, and the method can realize higher resolution in a mode of generating a virtual aperture so as to expand the aperture of the antenna.
Referring to the right part of fig. 1, 5 transmitting antennas and 5 receiving antennas in the antenna array can generate 25 virtual array elements, assuming that the distance between adjacent receiving antennas in the horizontal direction is
Figure 569461DEST_PATH_IMAGE020
With adjacent transmitting antennas at a vertical distance of
Figure 477374DEST_PATH_IMAGE021
Then, the maximum unambiguous angle measurement in the horizontal direction
Figure 550372DEST_PATH_IMAGE022
Can be calculated from the following equation:
Figure 518328DEST_PATH_IMAGE023
(1)
maximum unambiguous angle measurement in the vertical direction
Figure 676777DEST_PATH_IMAGE024
Can be calculated from the following equation:
Figure 450698DEST_PATH_IMAGE025
(2)
wherein,
Figure 643782DEST_PATH_IMAGE018
is the wavelength of the millimeter wave employed. Suppose that
Figure 782639DEST_PATH_IMAGE020
Value of 0.5
Figure 183313DEST_PATH_IMAGE018
Then the corresponding maximum unambiguous angle measurement in the horizontal direction
Figure 760925DEST_PATH_IMAGE022
In the range of
Figure 746198DEST_PATH_IMAGE026
The two-dimensional angle calculation method for the millimeter wave radar provided by the embodiment of the application is described below. Referring to fig. 2, a method for calculating a two-dimensional angle of a millimeter wave radar according to an embodiment of the present application may include the following steps:
step S101, selecting a plurality of virtual array elements from a target virtual array to form three groups of target array elements.
The target virtual arrays are evenly distributed in the horizontal direction and the vertical direction respectively, the first group of target array elements comprises at least two virtual array elements in the same horizontal direction, the second group of target array elements comprises at least two virtual array elements in the same vertical direction, and the third group of target array elements comprises at least two virtual array elements which are not in the same horizontal direction and are not in the same vertical direction. Illustratively, based on the target virtual array shown in fig. 1, a target array element may be selected as shown in fig. 3 or fig. 4. It is understood that the antenna array provided by the embodiments of the present application is merely exemplary, and in fact, the target virtual array generated by the antenna array only satisfies the following condition: the millimeter wave radar two-dimensional angle calculation method can be applied to the millimeter wave radar two-dimensional angle calculation method provided by the embodiment of the application.
Step S102, based on the echo data matrix of each virtual array element in each group of target array elements, determining an angle set of each target in the direction corresponding to the group of target array elements.
Specifically, a first angle set of each target in the horizontal direction is determined based on an echo data matrix of each virtual array element in a first group of target array elements; determining a second angle set of each target in the vertical direction based on the echo data matrix of each virtual array element in the second group of target array elements; and determining a third angle set of each target in the oblique direction based on the echo data matrix of each virtual array element in the third group of target array elements. It is understood that the first angle set, the second angle set and the third angle set contain equal number of elements, and all are the number of targets.
The vehicle-mounted millimeter wave radar generally uses a chirp continuous wave (chirp), for original echo data of the radar, how many virtual array elements are, how many echo data matrixes are corresponding to the original echo data, for example, the corresponding virtual array elements in fig. 1, 25 radar echo data matrixes can be obtained under radar single-frame observation, and the size of each radar echo data matrix is equal to that of the corresponding virtual array element
Figure 118274DEST_PATH_IMAGE027
Wherein
Figure 985736DEST_PATH_IMAGE028
for the number of sampling points of the ADC (analog to digital conversion) in each receiving antenna for a single chirp signal,
Figure 39142DEST_PATH_IMAGE029
the number of chirp transmitted by a single antenna under a single frame (generally, a transmitting antenna may transmit a number of chirps consecutively at a certain period, such as M chirps, which constitute a frame).
Because the first group of target array elements are virtual array elements in the target virtual array in the same horizontal direction, the echo data matrix of each virtual array element in the first group of target array elements contains angle information in the horizontal direction, namely the first angle is the angle of the target in the horizontal direction; similarly, since the second group of target array elements are virtual array elements in the target virtual array in the same vertical direction, the echo data matrix of each virtual array element in the second group of target array elements includes angle information in the vertical direction, that is, the second angle is an angle of the target in the vertical direction; because the third group of target array elements are virtual array elements which are not in the same horizontal direction and not in the same vertical direction in the target virtual array, the echo data matrix of each virtual array element in the third group of target array elements contains joint information in the horizontal direction and the vertical direction, that is, the third angle is the angle of the target in the direction of the straight line formed by connecting each virtual array element in the third group of target array elements. Based on this, by combining the echo data matrix of each virtual array element in each group of target array elements, the angle set of each target in the direction corresponding to the group of target array elements can be calculated by using an FFT or beam forming method.
Step S103, determining the angle of each target based on the calculation formula of the echo signal, the first angle set, the second angle set and the third angle set.
Because the third angle set combines the information of each target in the horizontal direction and the vertical direction, the first angle set and the second angle set respectively represent the angle information of each target in the horizontal direction and the vertical direction, so that the echo signals of the third group of target array elements are described by using the first angle and the second angle, the angle matching condition of the first angle and the second angle can be determined, and the angle of the target can be determined according to the well matched first angle and second angle.
The method comprises the steps of firstly selecting a plurality of virtual array elements from a target virtual array to form three groups of target array elements. The target virtual arrays are respectively and uniformly arranged in the horizontal direction and the vertical direction, namely, the intervals of all virtual array elements of the target virtual arrays in the same horizontal direction are consistent, and the intervals of all virtual array elements of the target virtual arrays in the same vertical direction are consistent. Then, based on the echo data matrix of each virtual array element in the first group of target array elements, determining a first angle set of the target in the horizontal direction; determining a second angle set of the target in the vertical direction based on the echo data matrix of each virtual array element in the second group of target array elements; and determining a third angle set of the target in the oblique direction based on the echo data matrix of each virtual array element in the third group of target array elements. The number of elements contained in the first angle set, the second angle set and the third angle set is equal to the number of the target elements. Because the first group of target array elements comprises at least two virtual array elements in the same horizontal direction, the second group of target array elements comprises at least two virtual array elements in the same vertical direction, and the third group of target array elements comprises at least two virtual array elements which are not in the same horizontal direction and not in the same vertical direction, therefore, each first angle in the first angle set embodies the angle information of the target in the horizontal direction, each second angle in the second angle set embodies the angle information of the target in the vertical direction, and each third angle in the third angle set embodies the associated information of the target in the horizontal direction and the vertical direction. And finally, determining the angle of each target based on a preset calculation formula of the echo signal, the first angle set, the second angle set and the third angle set. According to the method and the device, only part of virtual array elements are adopted for calculation, calculation of unnecessary information is reduced to a great extent, the calculated amount cannot be increased due to the increase of the number of the array elements in the target virtual array, the calculation process is simple, and the two-dimensional angle estimation of the target angle can be completed efficiently.
In some embodiments of the present application, the target virtual array is a rectangular array, and the step S101 of selecting a plurality of virtual array elements from the target virtual array to form three groups of target array elements may include:
s1, a first group of target array elements is formed by each virtual array element in any row of the target virtual array.
S2, a second group of target array elements is formed by each virtual array element in any column of the target virtual array.
And S3, forming a third group of target array elements by each virtual array element on the diagonal line or the sub diagonal line in the target virtual array.
In order to utilize the virtual array elements in the antenna array to the maximum, the virtual array elements can be taken as many as possible in the same horizontal direction, the same vertical direction and the diagonal line or the secondary diagonal line respectively, the more the selected virtual array elements are, the larger the aperture is, and the better the resolution is. That is, the scheme of fig. 3 may have a higher resolution than the selection of the target array elements of each group shown in fig. 4.
Specifically, for a uniformly arranged two-dimensional transmit-receive array like that shown in fig. 1, the angular resolution calculation formula:
Figure 941239DEST_PATH_IMAGE030
(3)
wherein,
Figure 421899DEST_PATH_IMAGE028
as to the number of virtual array elements on a row/column,Dis the distance between the two array elements,
Figure 42236DEST_PATH_IMAGE031
is the pore size, here equals
Figure 696072DEST_PATH_IMAGE032
Figure 718254DEST_PATH_IMAGE018
The wavelength corresponding to the carrier frequency is taken as the wavelength corresponding to 77GHz,
Figure 435062DEST_PATH_IMAGE033
at an angle to the normal of the radar (typically taking the value 0, i.e. 0)
Figure 214799DEST_PATH_IMAGE034
). Then, it can be seen from equation (3) that the larger the aperture, the higher the angular resolution, i.e. the closer the angle is to the target can be resolved.
In some embodiments of the present application, the above-mentioned target virtual array may not be a rectangular array, and as shown in fig. 5, the target virtual array at least includes two horizontally oriented virtual array elements to perform angle measurement in the horizontal direction, two vertically oriented virtual array elements to perform angle measurement in the vertical direction, and two non-horizontally oriented virtual array elements to obtain information related to two directions.
Based on this, in some embodiments of the present application, the target virtual array is a non-rectangular array, and the step S101 of selecting a plurality of virtual array elements from the target virtual array to form three groups of target array elements may include:
s1, a row with the largest number of virtual array elements is determined from the target virtual array, and each virtual array element in the row constitutes a first group of target array elements.
S2, a column with the largest number of virtual array elements is determined from the target virtual array, and each virtual array element in the column constitutes a second group of target array elements.
And S3, determining a plurality of virtual array elements which are not in the same row and are not in the same column from the target virtual array, wherein the centers of the plurality of virtual array elements are connected into a line, the intervals of the adjacent array elements in the plurality of virtual array elements in the horizontal direction are consistent, and the intervals of the adjacent array elements in the plurality of virtual array elements in the vertical direction are consistent.
In some embodiments of the present application, the step S102 of determining the first angle set of the target in the horizontal direction based on the echo data matrix of each virtual array element in the first group of target array elements may include:
and S1, performing 2D-FFT on each echo data matrix of each virtual array element in the first group of target array elements to obtain each range-Doppler spectrum matrix.
S2, performing Non-correlation accumulation (Non-coherent combining) processing on each range-doppler spectrum matrix to obtain a range-doppler spectrum matrix after accumulation.
S3, performing CFAR detection (Constant False Alarm Rate) on the accumulated range-doppler spectrum matrix to obtain a range and velocity index of the target.
S4, a first set of angles of the target is determined based on the respective range-Doppler spectrum matrices and the range and velocity indices of the target.
Illustratively, the first set of angles of the target may be determined by finding extrema or the like based on the range-doppler spectrum matrix and the range and velocity indices of the target.
Generally, since the echo data matrices are obtained from a single range and velocity gate, the range resolution and velocity resolution can be made small enough by setting parameters such as bandwidth, for example, only about 3 targets can be allowed under the same range and velocity gate, and then the angle corresponding to the first group of target array elements can be found based on the criterion and a set threshold (for example, -3dB of the maximum energy value).
Likewise, with reference to the above method, a second set of angles of the target in the vertical direction and a third set of angles in the diagonal direction may be determined.
It is understood that, in the process of determining the first angle set, the second angle set and the third angle set, the above steps S1 to S3 may be performed in combination, that is, in S1, 2D-FFT is performed on each echo data matrix of each virtual array element in the entire target virtual array, and after performing S2 and S3, the distance and velocity index of the target is obtained. In S4, the corresponding elements in each range-doppler spectrum matrix are selected to determine each angle set. Aiming at the first angle set, selecting elements corresponding to a first group of target array elements in each range-Doppler spectrum matrix, and determining the first angle set of the target by combining the range and speed indexes of the target; aiming at the second angle set, selecting elements corresponding to a second group of target array elements in each range-Doppler spectrum matrix, and determining the second angle set of the target by combining the distance and the speed index of the target; and aiming at the third angle set, selecting elements corresponding to a third group of target array elements in each distance-Doppler spectrum matrix, and determining the third angle set of the target by combining the distance and the speed index of the target.
In some embodiments of the present application, the step of determining the first angle set of the target based on the range-doppler spectrum matrices and the range and velocity indexes of the target at S4 may include:
s41, each element of the range and velocity index corresponding to the target is selected from each range-doppler spectrum matrix.
S42, new line data is constructed based on the elements.
Wherein the position of each element in the row coincides with the position of the element in the original range-doppler spectrum matrix.
S43, performing FFT or beam forming on the line data, finding the extremum to obtain the possible angles of the target, and forming the first angle set of the target by the angles.
In some embodiments of the present application, the step S103 of determining the angle of each target based on the calculation formula of the echo signal, the first angle set, the second angle set, and the third angle set may include:
s1, matching the first angle in the first angle set and the second angle in the second angle set according to equation (3), to obtain a matched angle combination:
Figure 672325DEST_PATH_IMAGE035
(4)
each angle combination comprises a first angle and a second angle, and each angle combination corresponds to one target.
Figure 549015DEST_PATH_IMAGE002
The distance of the adjacent target array elements in the third group of target array elements in the horizontal direction,
Figure 433794DEST_PATH_IMAGE003
the distance of the adjacent target array elements in the third group of target array elements in the vertical direction is shown. For example, for the third set of target array elements shown in fig. 3 or 4,
Figure 700827DEST_PATH_IMAGE036
Figure 962044DEST_PATH_IMAGE037
Figure 693240DEST_PATH_IMAGE004
is a first angle in the first set of angles,
Figure 686604DEST_PATH_IMAGE005
is a second angle in the second set of angles,
Figure 503250DEST_PATH_IMAGE006
is a third angle in the third set of angles.
S2, determining each angle combination as the angle of a target.
Assume that the first angle set, the second angle set and the third angle set each include 3 angles, which are
Figure 568158DEST_PATH_IMAGE038
Figure 91543DEST_PATH_IMAGE039
Figure 315195DEST_PATH_IMAGE040
After the matching process, three groups of angle combinations are obtained, namely
Figure 556821DEST_PATH_IMAGE041
Figure 159841DEST_PATH_IMAGE042
And
Figure 600049DEST_PATH_IMAGE043
then these three sets of angle combinations respectively correspond to a target.
Wherein the above equation is obtained based on the following phase relationship: for a particular target, let its horizontal direction angle with respect to the radar be
Figure 263112DEST_PATH_IMAGE004
At a vertical angle of
Figure 726454DEST_PATH_IMAGE005
The angle measured by the array where the third group of target array elements is
Figure 133165DEST_PATH_IMAGE006
Then, the echo signals received by each array element in the third group of target array elements have two equivalent expression modes:
Figure 427880DEST_PATH_IMAGE044
(5)
Figure 199527DEST_PATH_IMAGE045
(6)
corresponding to the virtual matrix shown in figure 3,
Figure 212482DEST_PATH_IMAGE036
Figure 94987DEST_PATH_IMAGE037
Figure 575035DEST_PATH_IMAGE018
is the wavelength of the carrier frequency and,
Figure 704534DEST_PATH_IMAGE046
corresponding to each array element, for the virtual array shown in fig. 3,
Figure 204785DEST_PATH_IMAGE028
and (5). Since equation (5) and equation (6) are equivalent, equation (4) can be obtained by combining these two equations.
In some embodiments of the present application, the process of matching the first angle in the first angle set and the second angle in the second angle set according to equation (4) at S1 may include:
s1, constructing a matrix
Figure 953298DEST_PATH_IMAGE007
Each element in the matrix is calculated by the following equation:
Figure 222606DEST_PATH_IMAGE008
s2, constructing a vector list
Figure 398372DEST_PATH_IMAGE009
The vector list is calculated by the following equation
Figure 382990DEST_PATH_IMAGE009
Each vector of (a):
Figure 669615DEST_PATH_IMAGE010
s3, for each
Figure 59008DEST_PATH_IMAGE011
Determining the same as
Figure 140096DEST_PATH_IMAGE011
Target with the closest value of
Figure 614940DEST_PATH_IMAGE012
And target the object
Figure 705256DEST_PATH_IMAGE012
In (1)
Figure 949155DEST_PATH_IMAGE013
And
Figure 466724DEST_PATH_IMAGE014
are determined as angle combinations matched with each other.
Wherein,
Figure 100968DEST_PATH_IMAGE015
is the number of the targets that are to be processed,
Figure 732325DEST_PATH_IMAGE013
is a first angle in the first set of angles,
Figure 96310DEST_PATH_IMAGE014
is a second angle in the second set of angles,
Figure 784781DEST_PATH_IMAGE016
is a third angle in the third set of angles;
Figure 906320DEST_PATH_IMAGE002
is prepared by reacting with
Figure 338439DEST_PATH_IMAGE013
Corresponding target array element, and
Figure 556930DEST_PATH_IMAGE014
the distance of the corresponding target array element in the horizontal direction;
Figure 416302DEST_PATH_IMAGE003
is prepared by reacting with
Figure 25138DEST_PATH_IMAGE013
Corresponding target array element, and
Figure 995368DEST_PATH_IMAGE014
the corresponding target array element, the distance in the vertical direction. Illustratively, corresponding to the virtual matrix shown in FIG. 3,
Figure 68366DEST_PATH_IMAGE036
Figure 36322DEST_PATH_IMAGE037
in practical applications, due to the influences of insufficient angle measurement accuracy, system errors and the like, the relational expressions cannot be strictly equal, and therefore certain errors are allowed.
The matching process in the above embodiment may be represented by a code as follows:
first, based on
Figure 191841DEST_PATH_IMAGE013
Figure 965762DEST_PATH_IMAGE014
And equation (3) constructs one
Figure 893267DEST_PATH_IMAGE047
Size matrix to be compared
Figure 32124DEST_PATH_IMAGE048
Based on
Figure 677869DEST_PATH_IMAGE016
And equation (3) constructs a length of
Figure 255481DEST_PATH_IMAGE015
To be compared vector
Figure 37492DEST_PATH_IMAGE009
Figure 347251DEST_PATH_IMAGE049
Then pair
Figure 480292DEST_PATH_IMAGE009
Traversing the values in the vector and comparing the values with the values obtained
Figure 533699DEST_PATH_IMAGE048
And comparing elements in the matrix, finding the element closest to the value, wherein the values of alpha and beta corresponding to the element are the angles in the horizontal and vertical directions corresponding to a certain target, and thus completing pairing.
Figure 435796DEST_PATH_IMAGE050
The List array obtained after the matching comprises the matching results of all targets: each row corresponds to a target, and the second column of each row corresponds to the angle corresponding to the horizontal direction of the target
Figure 981702DEST_PATH_IMAGE013
The third column corresponds to the angle corresponding to the target vertical direction
Figure 602039DEST_PATH_IMAGE014
The millimeter wave radar two-dimensional angle calculation device provided by the embodiment of the present application is described below, and the millimeter wave radar two-dimensional angle calculation device described below and the millimeter wave radar two-dimensional angle calculation method described above may be referred to in a mutually corresponding manner.
Referring to fig. 6, the millimeter wave radar two-dimensional angle calculating device provided in the embodiment of the present application may include:
a target array element selecting unit 21, configured to select a plurality of virtual array elements from a target virtual array to form three groups of target array elements, where the target virtual array is uniformly arranged in a horizontal direction and a vertical direction, a first group of target array elements includes at least two virtual array elements in the same horizontal direction, a second group of target array elements includes at least two virtual array elements in the same vertical direction, and a third group of target array elements includes at least two virtual array elements that are not in the same horizontal direction and are not in the same vertical direction;
a single angle determining unit 22, configured to determine a first angle set of each target in the horizontal direction based on an echo data matrix of each virtual array element in the first group of target array elements, determine a second angle set of each target in the vertical direction based on an echo data matrix of each virtual array element in the second group of target array elements, and determine a third angle set of each target in the oblique direction based on an echo data matrix of each virtual array element in the third group of target array elements, where the numbers of elements included in the first angle set, the second angle set, and the third angle set are equal and are the numbers of targets;
an angle matching calculation unit 23, configured to determine an angle of each target based on a calculation formula of an echo signal, the first angle set, the second angle set, and the third angle set.
In some embodiments of the present application, the process that the target virtual array is a rectangular array, and the selecting unit 21 selects a plurality of virtual array elements from the target virtual array to form three groups of target array elements may include:
forming a first group of target array elements by all virtual array elements in any row in the target virtual array;
forming a second group of target array elements by each virtual array element in any column in the target virtual array;
and forming a third group of target array elements by each virtual array element on the diagonal line or the sub diagonal line in the target virtual array.
In some embodiments of the present application, the target virtual array is a non-rectangular array, and the process of selecting a plurality of virtual array elements from the target virtual array by the target array element selecting unit 21 to form three groups of target array elements may include:
determining a row with the maximum number of virtual array elements from a target virtual array, and forming a first group of target array elements by all the virtual array elements in the row;
determining a column with the maximum number of virtual array elements from a target virtual array, and forming a second group of target array elements by all the virtual array elements in the column;
determining a plurality of virtual array elements which are not in the same row and are not in the same column from the target virtual array, wherein the centers of the plurality of virtual array elements are connected into a line, the intervals of the adjacent array elements in the plurality of virtual array elements in the horizontal direction are consistent, and the intervals of the adjacent array elements in the plurality of virtual array elements in the vertical direction are consistent.
In some embodiments of the present application, the distance of adjacent target array elements in the target virtual array in the horizontal direction
Figure 193558DEST_PATH_IMAGE002
Satisfying the following equation:
Figure 950161DEST_PATH_IMAGE017
wherein,
Figure 929618DEST_PATH_IMAGE018
being the wavelength of the millimeter wave used,
Figure 771673DEST_PATH_IMAGE019
the angle measurement range is a preset maximum non-fuzzy angle measurement range.
In some embodiments of the present application, the process of determining, by the single angle determining unit 22, the first angle set of the target in the horizontal direction based on the echo data matrix of each virtual array element in the first group of target array elements may include:
performing two-dimensional fast Fourier transform (2D-FFT) on each echo data matrix of each virtual array element in the first group of target array elements to obtain a distance-Doppler spectrum matrix;
performing non-correlation accumulation processing on the distance-Doppler spectrum matrix to obtain a distance-Doppler spectrum matrix after accumulation;
performing constant false alarm CFAR detection on the accumulated distance-Doppler spectrum matrix to obtain a distance and speed index of the target;
a first set of angles of the target is determined based on the range-doppler spectrum matrix and the range and velocity indices of the target.
In some embodiments of the present application, the process of determining the angle of each target by the angle matching calculation unit 23 based on the calculation formula of the echo signal, the first angle set, the second angle set, and the third angle set may include:
matching a first angle in the first angle set and a second angle in the second angle set according to the following equation to obtain a matched angle combination:
Figure 166882DEST_PATH_IMAGE001
determining each angle combination as an angle of a target, wherein each angle combination comprises a first angle and a second angle;
wherein,
Figure 43571DEST_PATH_IMAGE002
is the distance of the adjacent target array elements in the third group of target array elements in the horizontal direction,
Figure 928350DEST_PATH_IMAGE003
is the distance of the adjacent target array elements in the third group of target array elements in the vertical direction,
Figure 195384DEST_PATH_IMAGE004
for a first angle in the first set of angles,
Figure 456601DEST_PATH_IMAGE005
is a second angle of the second set of angles,
Figure 208304DEST_PATH_IMAGE006
is a third angle in the third set of angles.
In some embodiments of the present application, the process of matching the first angle in the first angle set and the second angle in the second angle set by the angle matching calculation unit 23 according to the following equation may include:
construction matrix
Figure 263985DEST_PATH_IMAGE007
Each element in the matrix is calculated by the following equation:
Figure 18314DEST_PATH_IMAGE008
constructing a vector list
Figure 817643DEST_PATH_IMAGE009
The vector list is calculated by the following equation
Figure 403345DEST_PATH_IMAGE009
Each vector of (a):
Figure 833189DEST_PATH_IMAGE010
for each
Figure 137132DEST_PATH_IMAGE011
Determining the same as
Figure 474572DEST_PATH_IMAGE011
Target with the closest value of
Figure 180360DEST_PATH_IMAGE012
And target the object
Figure 515526DEST_PATH_IMAGE012
In
Figure 306765DEST_PATH_IMAGE013
And
Figure 450826DEST_PATH_IMAGE014
determining the angle combinations matched with each other;
wherein,
Figure 948803DEST_PATH_IMAGE015
as to the number of the objects in question,
Figure 782767DEST_PATH_IMAGE013
for a first angle in the first set of angles,
Figure 795723DEST_PATH_IMAGE014
for a second angle in the second set of angles,
Figure 740545DEST_PATH_IMAGE016
is a third angle in the third set of angles;
Figure 93029DEST_PATH_IMAGE002
is prepared by reacting with
Figure 97894DEST_PATH_IMAGE013
Corresponding target array element, and
Figure 535828DEST_PATH_IMAGE014
the distance of the corresponding target array element in the horizontal direction;
Figure 18762DEST_PATH_IMAGE003
is prepared by reacting with
Figure 553649DEST_PATH_IMAGE013
Corresponding target array element, and
Figure 667098DEST_PATH_IMAGE014
the corresponding target array element, the distance in the vertical direction.
The millimeter wave radar two-dimensional angle calculating device provided by the embodiment of the application can be applied to millimeter wave radar two-dimensional angle calculating equipment such as a computer. Optionally, fig. 7 shows a block diagram of a hardware structure of the millimeter wave radar two-dimensional angle calculation device, and referring to fig. 7, the hardware structure of the millimeter wave radar two-dimensional angle calculation device may include: at least one processor 31, at least one communication interface 32, at least one memory 33 and at least one communication bus 34.
In the embodiment of the present application, the number of the processor 31, the communication interface 32, the memory 33 and the communication bus 34 is at least one, and the processor 31, the communication interface 32 and the memory 33 complete the communication with each other through the communication bus 34;
the processor 31 may be a central processing unit CPU, or an application Specific Integrated circuit asic, or one or more Integrated circuits configured to implement embodiments of the present application, etc.;
the memory 33 may include a high-speed RAM memory, and may further include a non-volatile memory (non-volatile memory) or the like, such as at least one disk memory;
wherein the memory 33 stores a program and the processor 31 may call the program stored in the memory 33 for:
selecting a plurality of virtual array elements from a target virtual array to form three groups of target array elements, wherein the target virtual array is uniformly distributed in the horizontal direction and the vertical direction respectively, the first group of target array elements comprises at least two virtual array elements in the same horizontal direction, the second group of target array elements comprises at least two virtual array elements in the same vertical direction, and the third group of target array elements comprises at least two virtual array elements which are not in the same horizontal direction and not in the same vertical direction;
determining a first angle set of each target in the horizontal direction based on an echo data matrix of each virtual array element in the first group of target array elements, determining a second angle set of each target in the vertical direction based on an echo data matrix of each virtual array element in the second group of target array elements, and determining a third angle set of each target in the oblique direction based on an echo data matrix of each virtual array element in the third group of target array elements, wherein the first angle set, the second angle set and the third angle set contain equal number of elements and are the number of targets;
determining the angle of each target based on the calculation formula of the echo signal, the first angle set, the second angle set and the third angle set.
Alternatively, the detailed function and the extended function of the program may be as described above.
An embodiment of the present application further provides a storage medium, where the storage medium may store a program adapted to be executed by a processor, where the program is configured to:
selecting a plurality of virtual array elements from a target virtual array to form three groups of target array elements, wherein the target virtual array is uniformly distributed in the horizontal direction and the vertical direction respectively, the first group of target array elements comprises at least two virtual array elements in the same horizontal direction, the second group of target array elements comprises at least two virtual array elements in the same vertical direction, and the third group of target array elements comprises at least two virtual array elements which are not in the same horizontal direction and not in the same vertical direction;
determining a first angle set of each target in the horizontal direction based on an echo data matrix of each virtual array element in the first group of target array elements, determining a second angle set of each target in the vertical direction based on an echo data matrix of each virtual array element in the second group of target array elements, and determining a third angle set of each target in the oblique direction based on an echo data matrix of each virtual array element in the third group of target array elements, wherein the first angle set, the second angle set and the third angle set contain the same number of elements, and the number of elements is the same as that of the targets;
determining the angle of each target based on the calculation formula of the echo signal, the first angle set, the second angle set and the third angle set.
Alternatively, the detailed function and the extended function of the program may refer to the above description.
In summary, the following steps:
the method comprises the steps of firstly selecting a plurality of virtual array elements from a target virtual array to form three groups of target array elements. The target virtual arrays are respectively and uniformly arranged in the horizontal direction and the vertical direction, namely, the intervals of all virtual array elements of the target virtual arrays in the same horizontal direction are consistent, and the intervals of all virtual array elements of the target virtual arrays in the same vertical direction are consistent. Then, based on the echo data matrix of each virtual array element in the first group of target array elements, determining a first angle set of the target in the horizontal direction; determining a second angle set of the target in the vertical direction based on the echo data matrix of each virtual array element in the second group of target array elements; and determining a third angle set of the target in the oblique direction based on the echo data matrix of each virtual array element in the third group of target array elements. The number of elements contained in the first angle set, the second angle set and the third angle set is equal to the number of the target elements. Because the first group of target array elements comprises at least two virtual array elements in the same horizontal direction, the second group of target array elements comprises at least two virtual array elements in the same vertical direction, and the third group of target array elements comprises at least two virtual array elements which are not in the same horizontal direction and not in the same vertical direction, therefore, each first angle in the first angle set embodies the angle information of the target in the horizontal direction, each second angle in the second angle set embodies the angle information of the target in the vertical direction, and each third angle in the third angle set embodies the associated information of the target in the horizontal direction and the vertical direction. And finally, determining the angle of each target based on a preset calculation formula of the echo signal, the first angle set, the second angle set and the third angle set. According to the method and the device, only part of virtual array elements are adopted for calculation, calculation of unnecessary information is reduced to a great extent, the calculated amount cannot be increased due to the increase of the number of the array elements in the target virtual array, the calculation process is simple, and two-dimensional angle estimation of the target angle can be completed efficiently.
Finally, it should also be noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other identical elements in a process, method, article, or apparatus that comprises the element.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, the embodiments may be combined as needed, and the same and similar parts may be referred to each other.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present application. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the application. Thus, the present application is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (10)

1. A millimeter wave radar two-dimensional angle calculation method is characterized by comprising the following steps:
selecting a plurality of virtual array elements from a target virtual array to form three groups of target array elements, wherein the target virtual array is uniformly distributed in the horizontal direction and the vertical direction respectively, the first group of target array elements comprises at least two virtual array elements in the same horizontal direction, the second group of target array elements comprises at least two virtual array elements in the same vertical direction, and the third group of target array elements comprises at least two virtual array elements which are not in the same horizontal direction and are not in the same vertical direction;
determining a first angle set of each target in the horizontal direction based on an echo data matrix of each virtual array element in the first group of target array elements, determining a second angle set of each target in the vertical direction based on an echo data matrix of each virtual array element in the second group of target array elements, and determining a third angle set of each target in the oblique direction based on an echo data matrix of each virtual array element in the third group of target array elements, wherein the first angle set, the second angle set and the third angle set contain equal number of elements and are the number of targets;
determining the angle of each target based on the calculation formula of the echo signal, the first angle set, the second angle set and the third angle set.
2. The method of claim 1, wherein the determining the angle of each target based on the calculation formula of the echo signal, the first angle set, the second angle set, and the third angle set comprises:
matching a first angle in the first angle set and a second angle in the second angle set according to the following equation to obtain a matched angle combination:
Figure 502055DEST_PATH_IMAGE001
determining each angle combination as an angle of a target, wherein each angle combination comprises a first angle and a second angle;
wherein,
Figure 763272DEST_PATH_IMAGE002
is the distance of the adjacent target array elements in the third group of target array elements in the horizontal direction,
Figure 432151DEST_PATH_IMAGE003
is the distance of the adjacent target array elements in the third group of target array elements in the vertical direction,
Figure 487831DEST_PATH_IMAGE004
for a first angle in the first set of angles,
Figure 304478DEST_PATH_IMAGE005
for a second angle in the second set of angles,
Figure 103806DEST_PATH_IMAGE006
is a third angle in the third set of angles.
3. The method of claim 2, wherein matching a first angle of the first set of angles to a second angle of the second set of angles is accomplished according to the following equation, comprising:
construction matrix
Figure 689508DEST_PATH_IMAGE007
Each element in the matrix is calculated by the following equation:
Figure DEST_PATH_IMAGE008
constructing a list of vectors
Figure 243986DEST_PATH_IMAGE009
The vector list is calculated by the following equation
Figure 547929DEST_PATH_IMAGE009
Each vector of (a):
Figure 823052DEST_PATH_IMAGE010
for each
Figure DEST_PATH_IMAGE011
Determining the same as
Figure 591157DEST_PATH_IMAGE011
Target with the closest value of
Figure 988640DEST_PATH_IMAGE012
And target the object
Figure 779879DEST_PATH_IMAGE012
In (1)
Figure 921010DEST_PATH_IMAGE013
And
Figure 543621DEST_PATH_IMAGE014
determining the angle combinations matched with each other;
wherein,
Figure 377585DEST_PATH_IMAGE015
as to the number of the objects in question,
Figure 390541DEST_PATH_IMAGE013
is a first angle of the first set of angles,
Figure 335363DEST_PATH_IMAGE014
for a second angle in the second set of angles,
Figure 750164DEST_PATH_IMAGE016
is a third angle in the third set of angles.
4. The method of claim 1, wherein determining the first set of angles of the target in the horizontal direction based on the echo data matrix of each virtual array element in the first set of target array elements comprises:
performing two-dimensional fast Fourier transform (2D-FFT) on each echo data matrix of each virtual array element in the first group of target array elements to obtain a distance-Doppler spectrum matrix;
performing non-correlation accumulation processing on the distance-Doppler spectrum matrix to obtain a distance-Doppler spectrum matrix after accumulation;
performing constant false alarm CFAR detection on the accumulated distance-Doppler spectrum matrix to obtain a distance and speed index of the target;
a first set of angles of the target is determined based on the range-doppler spectrum matrix and the range and velocity indices of the target.
5. The method of claim 1, wherein the target virtual array is a rectangular array, and the selecting a plurality of virtual array elements from the target virtual array to form three groups of target array elements comprises:
forming a first group of target array elements by all virtual array elements in any row in the target virtual array;
forming a second group of target array elements by each virtual array element in any column in the target virtual array;
and forming a third group of target array elements by each virtual array element on the diagonal line or the sub diagonal line in the target virtual array.
6. The method of claim 1, wherein the target virtual array is a non-rectangular array, and the selecting a number of virtual array elements from the target virtual array to form three groups of target array elements comprises:
determining a row with the maximum number of virtual array elements from a target virtual array, and forming a first group of target array elements by all the virtual array elements in the row;
determining a column with the maximum number of virtual array elements from a target virtual array, and forming a second group of target array elements by all the virtual array elements in the column;
determining a plurality of virtual array elements which are not in the same row and are not in the same column from the target virtual array, wherein the centers of the plurality of virtual array elements are connected into a line, the intervals of the adjacent array elements in the plurality of virtual array elements in the horizontal direction are consistent, and the intervals of the adjacent array elements in the plurality of virtual array elements in the vertical direction are consistent.
7. The method of claim 1, wherein the distance between adjacent target array elements in the target virtual array in the horizontal direction
Figure 755029DEST_PATH_IMAGE002
Satisfying the following equation:
Figure 255280DEST_PATH_IMAGE017
wherein,
Figure 675897DEST_PATH_IMAGE018
for the wavelength of the millimeter wave used,
Figure 273101DEST_PATH_IMAGE019
the angle measurement range is a preset maximum non-fuzzy angle measurement range.
8. A millimeter wave radar two-dimensional angle calculation device, comprising:
the target array element selecting unit is used for selecting a plurality of virtual array elements from a target virtual array to form three groups of target array elements, wherein the target virtual array is uniformly distributed in the horizontal direction and the vertical direction respectively, the first group of target array elements comprises at least two virtual array elements in the same horizontal direction, the second group of target array elements comprises at least two virtual array elements in the same vertical direction, and the third group of target array elements comprises at least two virtual array elements which are not in the same horizontal direction and are not in the same vertical direction;
a single angle determining unit, configured to determine a first angle set of each target in a horizontal direction based on an echo data matrix of each virtual array element in the first group of target array elements, determine a second angle set of each target in a vertical direction based on an echo data matrix of each virtual array element in the second group of target array elements, and determine a third angle set of each target in an oblique direction based on an echo data matrix of each virtual array element in the third group of target array elements, where the numbers of elements included in the first angle set, the second angle set, and the third angle set are equal and are the numbers of targets;
an angle matching calculation unit, configured to determine an angle of each target based on a calculation formula of an echo signal, the first angle set, the second angle set, and the third angle set.
9. A millimeter wave radar two-dimensional angle calculation device, comprising: a memory and a processor;
the memory is used for storing programs;
the processor is configured to execute the program to implement the steps of the millimeter wave radar two-dimensional angle calculation method according to any one of claims 1 to 7.
10. A storage medium having stored thereon a computer program, wherein the computer program, when executed by a processor, implements the steps of the millimeter wave radar two-dimensional angle calculation method according to any one of claims 1 to 7.
CN202210818680.2A 2022-07-13 2022-07-13 Millimeter wave radar two-dimensional angle calculation method and device and related equipment Active CN114879138B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210818680.2A CN114879138B (en) 2022-07-13 2022-07-13 Millimeter wave radar two-dimensional angle calculation method and device and related equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210818680.2A CN114879138B (en) 2022-07-13 2022-07-13 Millimeter wave radar two-dimensional angle calculation method and device and related equipment

Publications (2)

Publication Number Publication Date
CN114879138A CN114879138A (en) 2022-08-09
CN114879138B true CN114879138B (en) 2022-09-27

Family

ID=82683080

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210818680.2A Active CN114879138B (en) 2022-07-13 2022-07-13 Millimeter wave radar two-dimensional angle calculation method and device and related equipment

Country Status (1)

Country Link
CN (1) CN114879138B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521472A (en) * 2012-01-04 2012-06-27 电子科技大学 Method for constructing thinned MIMO (Multiple Input Multiple Output) planar array radar antenna
CN104777485A (en) * 2015-04-20 2015-07-15 西安交通大学 Three-dimensional wide-beam small-region rapid cavitating and imaging method of ultrasonic two-dimensional planar array
CN104931919A (en) * 2015-06-08 2015-09-23 西安电子科技大学 Direction finding method for two dimensional plane annular array
CN105785315A (en) * 2016-04-20 2016-07-20 西北工业大学 Angle measurement method of linear array, covariance matrix and main diagonal reconstruction
CN107918108A (en) * 2017-11-14 2018-04-17 重庆邮电大学 A kind of uniform circular array 2-d direction finding method for quick estimating
CN108957391A (en) * 2018-07-24 2018-12-07 北京理工大学 A kind of estimating two-dimensional direction-of-arrival method of the inverted-L antenna battle array based on nested array
CN109188386A (en) * 2018-08-30 2019-01-11 河海大学 Based on the MIMO radar high-resolution method for parameter estimation for improving two dimensional ESPRIT algorithm
CN113126021A (en) * 2021-04-19 2021-07-16 电子科技大学 Single-snapshot two-dimensional DOA estimation method based on three parallel linear arrays
CN113314832A (en) * 2021-06-15 2021-08-27 东南大学 Millimeter wave vehicle-mounted MIMO radar antenna array device and design method

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8354972B2 (en) * 2007-06-06 2013-01-15 Fractus, S.A. Dual-polarized radiating element, dual-band dual-polarized antenna assembly and dual-polarized antenna array
KR20150100909A (en) * 2013-01-28 2015-09-02 후지쯔 가부시끼가이샤 Feedback method for channel state information, transmission method for channel state information reference signal, user equipment and base station

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521472A (en) * 2012-01-04 2012-06-27 电子科技大学 Method for constructing thinned MIMO (Multiple Input Multiple Output) planar array radar antenna
CN104777485A (en) * 2015-04-20 2015-07-15 西安交通大学 Three-dimensional wide-beam small-region rapid cavitating and imaging method of ultrasonic two-dimensional planar array
CN104931919A (en) * 2015-06-08 2015-09-23 西安电子科技大学 Direction finding method for two dimensional plane annular array
CN105785315A (en) * 2016-04-20 2016-07-20 西北工业大学 Angle measurement method of linear array, covariance matrix and main diagonal reconstruction
CN107918108A (en) * 2017-11-14 2018-04-17 重庆邮电大学 A kind of uniform circular array 2-d direction finding method for quick estimating
CN108957391A (en) * 2018-07-24 2018-12-07 北京理工大学 A kind of estimating two-dimensional direction-of-arrival method of the inverted-L antenna battle array based on nested array
CN109188386A (en) * 2018-08-30 2019-01-11 河海大学 Based on the MIMO radar high-resolution method for parameter estimation for improving two dimensional ESPRIT algorithm
CN113126021A (en) * 2021-04-19 2021-07-16 电子科技大学 Single-snapshot two-dimensional DOA estimation method based on three parallel linear arrays
CN113314832A (en) * 2021-06-15 2021-08-27 东南大学 Millimeter wave vehicle-mounted MIMO radar antenna array device and design method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"The coefficient compensation of damaged element in phase-array radar based on region search method and total weight optimization method";Shujin Gao等;《2019 European Microwave Conference in Central Europe (EuMCE)》;20191017;46-49 *
"超宽带宽角覆盖阵列天线技术研究";付桂林;《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》;20210515;I136-126 *

Also Published As

Publication number Publication date
CN114879138A (en) 2022-08-09

Similar Documents

Publication Publication Date Title
JP7369852B2 (en) Radar systems and vehicles
WO2021129581A1 (en) Signal processing method and apparatus
CN114879139B (en) Joint angle measurement method and device for vehicle-mounted 4D millimeter wave radar and related equipment
CN109521426B (en) Method and device for obtaining angle of target based on automobile radar
CN115327473B (en) Equivalent multi-snapshot 4D millimeter wave radar angle measurement method and device and related equipment
CN113109781B (en) Direction-of-arrival estimation method, radar and mobile device
CN112098970B (en) Speed ambiguity resolving algorithm for traffic microwave detection and related equipment
CN112612020B (en) Novel millimeter wave radar signal processing method
CN110865364A (en) Target resolving method of radar and terminal equipment
KR20200047411A (en) Radar and antenna built in radar
CN112098968A (en) Target positioning method and device for radar
CN116569062A (en) Positioning method and device
CN115097447A (en) MIMO radar monitoring system and monitoring method based on MIMO radar monitoring system
US20210323560A1 (en) Vehicle speed calculation method, system, device, and storage medium
CN112710998B (en) Speed ambiguity-resolving algorithm for microwave detection and related equipment
CN114879138B (en) Millimeter wave radar two-dimensional angle calculation method and device and related equipment
CN117368922A (en) Vehicle-mounted radar angle measurement method, terminal equipment and storage medium
JP2020020696A (en) Arrival direction estimating device and arrival direction estimating method
CN116106847B (en) Millimeter wave radar two-dimensional combined super-resolution angle measurement method, device and storage medium
CN110515066B (en) Vehicle-mounted millimeter wave radar and target height measuring method thereof
CN116500620A (en) Data processing method and device of millimeter wave radar, storage medium and unmanned vehicle
CN114779183B (en) Self-adaptive three-dimensional angle Doppler compensation method based on FDA-MIMO radar
CN118068313B (en) Millimeter wave radar and angle-of-arrival measurement method and device
CN113406615B (en) Target tracking method and device of binary phase modulation array radar
US20240361444A1 (en) Method for radar angle estimation

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