CN115421157A - Method and device for constructing radar array based on undirected adjacency graph - Google Patents

Method and device for constructing radar array based on undirected adjacency graph Download PDF

Info

Publication number
CN115421157A
CN115421157A CN202211298872.1A CN202211298872A CN115421157A CN 115421157 A CN115421157 A CN 115421157A CN 202211298872 A CN202211298872 A CN 202211298872A CN 115421157 A CN115421157 A CN 115421157A
Authority
CN
China
Prior art keywords
array
radar
sub
virtual
elements
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.)
Granted
Application number
CN202211298872.1A
Other languages
Chinese (zh)
Other versions
CN115421157B (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.)
Shenzhen University
Original Assignee
Shenzhen 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 Shenzhen University filed Critical Shenzhen University
Priority to CN202211298872.1A priority Critical patent/CN115421157B/en
Publication of CN115421157A publication Critical patent/CN115421157A/en
Application granted granted Critical
Publication of CN115421157B publication Critical patent/CN115421157B/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
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88Lidar systems specially adapted for specific applications
    • G01S17/89Lidar systems specially adapted for specific applications for mapping or imaging
    • 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/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • 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)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Mathematical Analysis (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Remote Sensing (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Electromagnetism (AREA)
  • Computing Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention relates to the technical field of radars, in particular to a method and a device for constructing a radar array based on a non-directional adjacency graph. The method comprises the steps of firstly simulating a receiving signal of the radar, then constructing a virtual array by using the simulated receiving signal of the radar, then processing the virtual array to obtain an undirected adjacent matrix, and judging whether undirected graphs of the undirected adjacent matrix are communicated or not, wherein the undirected graphs are communicated to show that the position of a hollow array element in the virtual array can be used for placing the radar. If the undirected graph is connected, the position of the hollow array element in the virtual array is suitable for placing the radar. And obtaining a subarray which is formed by the minimum array elements and can enable the undirected graph to be non-connected through traversal, and replacing the hollow array elements in the array to be supplemented with the non-hollow array elements to obtain the final target radar array. The invention can balance the problem between the aperture of the radar array and the accuracy of the reference data provided by the radar array.

Description

Method and device for constructing radar array based on undirected adjacency graph
Technical Field
The invention relates to the technical field of radars, in particular to a method and a device for constructing a radar array based on a non-directional adjacency graph.
Background
The realization of obstacle avoidance of unmanned driving, autonomous robots and unmanned aerial vehicles requires three-dimensional modeling of external complex environments. The unmanned system mainly relies on a laser radar point cloud image to construct an environment three-dimensional model. The laser radar has accurate distance measurement and high spatial resolution, and is an ideal environment sensor for automatic driving. However, as the wavelength range of the laser radar is in the nm level, when the laser radar encounters dust in the air, macromolecular aerosols such as water vapor and raindrops can generate refraction or scattering phenomena, and the laser radar cannot be used in severe weather such as rainstorm, snowstorm, haze and sand storm. In order to reduce the influence of the above natural phenomena on the lidar so that the lidar can provide necessary reference data for the unmanned technology, the prior art often adopts increasing the number of radars to increase the accuracy of the reference data provided by the array of radars. For example, an array of radars can provide reference data of the angle of arrival for the unmanned technology, and the resolution of the angle of arrival is larger when the aperture (the aperture is the size of the space occupied by the radars in the array) of the radar array is larger, so that the higher the resolution of the angle of arrival is, the safer the unmanned technology can perform. Conversely, the smaller the aperture of the radar array, the greater the resolution of the angle of arrival. And the radar array with large aperture increases the difficulty of using the radar. Therefore, it is necessary to balance the resolution of the radar array aperture and the angle of arrival, i.e. to improve the resolution of the angle of arrival (reference data) under the premise of limited radar array aperture, or to make a compromise between the radar array aperture and the quality of the reference data provided by the radar array difficult.
In summary, the prior art is difficult to balance the problem of radar array aperture and the quality of the reference data provided by the radar array.
Thus, there is a need for improvements and enhancements in the art.
Disclosure of Invention
In order to solve the technical problems, the invention provides a method and a device for constructing a radar array based on an undirected adjacency graph, and solves the problem that in the prior art, the aperture of the radar array and the quality of reference data provided by the radar array are difficult to balance.
In order to achieve the purpose, the invention adopts the following technical scheme:
in a first aspect, the present invention provides a method for constructing a radar array based on an undirected adjacency graph, wherein the method comprises:
generating a Hankel matrix by taking array elements in a virtual array as elements in the matrix and taking virtual receiving signals corresponding to the array elements as element values in the matrix, wherein the array elements are used for representing a radar, and the virtual receiving signals are used for simulating receiving signals of the radar;
setting a non-zero element value in the Hankel matrix as one and a zero element value as zero to obtain an undirected adjacency matrix;
when an undirected graph corresponding to the undirected adjacency matrix has connectivity, selecting the array elements from the virtual array to form a sub-array, and taking the sub-array formed by the minimum array elements corresponding to the undirected graph having the non-connectivity as an array to be supplemented;
and according to a convex optimization algorithm, filling the cavity array elements covered by the array to be filled with non-cavity array elements to obtain a target radar array, wherein the positions of the non-cavity array elements in the array to be filled before filling are used for placing the radar, and the positions of the cavity array elements in the array to be filled before filling are not placed with the radar.
In one implementation, the virtual array and the method for generating the virtual received signal include:
respectively endowing each array element in two preset arrays with simulation receiving signals, wherein the array elements are used for representing a radar, and the simulation receiving signals are used for simulating receiving signals of the radar;
and subtracting the array elements in the two preset arrays in pairs to obtain virtual arrays and virtual receiving signals corresponding to the virtual arrays, wherein the subtraction is to subtract the position of the array element of one of the preset arrays in one of the preset arrays by the position of the array element of the other preset array in the other preset array, and subtract the simulated receiving signals of the array elements of the one of the preset arrays by the simulated receiving signals of the array elements of the other preset array.
In one implementation manner, the total number of the array elements covered by one of the two preset arrays and the total number of the array elements covered by the other preset array are mutually prime numbers, where the array element pitch of one of the preset arrays is equal to the total number of the array elements of the other preset array multiplied by a wavelength, the array element pitch of the other preset array is equal to the total number of the array elements of one of the preset arrays multiplied by a wavelength, and the wavelength is a wavelength of a signal to be received by the radar.
In one implementation manner, the method for generating the simulated received signal includes:
setting a pitching angle when the radar receives a transmission signal of a signal source, wherein the signal source is used for transmitting the signal to the radar;
setting noise information of the environment where the radar is located;
and generating the simulation receiving signal according to the pitching angle, the transmitting signal, the noise information and the position of the array element in the preset array.
In one implementation, when an undirected graph corresponding to the undirected adjacency matrix has connectivity, selecting the array elements from the virtual array to form a sub-array, and using the sub-array formed by the minimum array elements corresponding to the undirected graph having the non-connectivity as an array to be complemented includes:
forming sub-arrays by using the array elements at the head end and the tail end of the virtual array, and if the undirected graphs corresponding to the sub-arrays have connectivity, adding the array elements except the array elements at the head end and the tail end to the sub-arrays one by one until the undirected graphs corresponding to the sub-arrays after the array elements are added have non-connectivity to obtain target sub-arrays;
and obtaining an array to be compensated according to the target subarray.
In an implementation manner, the forming a sub-array with the array elements at the head and the tail of the virtual array, and if the undirected graph corresponding to the sub-array has connectivity, adding the array elements except the array elements at the head and the tail to the sub-array one by one until the undirected graph corresponding to the sub-array after the addition of the array elements has non-connectivity, to obtain a target sub-array includes:
forming a first sub-array in the sub-arrays by using the array elements at the head end and the tail end of the virtual array;
forming the array elements except the array elements at the head end and the tail end in the virtual array into a second sub array in the sub arrays, wherein the second sub array and the first sub array are relatively prime arrays;
if the undirected graphs corresponding to the virtual sub-arrays formed by the first sub-array and the second sub-array have connectivity, the array elements of the virtual arrays are respectively added to the second sub-array and the first sub-array one by one until the undirected graphs corresponding to the virtual sub-arrays formed by the second sub-array and the first sub-array after the array elements are added have non-connectivity, and the second sub-array and the first sub-array after the array elements are added are taken as target sub-arrays.
In one implementation, the method for completing a hole array element covered by an array to be completed as a non-hole array element according to a convex optimization algorithm to obtain a target radar array, where a position of the non-hole array element in the array to be completed before completion is used for placing the radar, and a position of the hole array element in the array to be completed before completion is not placed with the radar includes:
adding unknown parameters to the non-cavity array elements in the array to be compensated to obtain virtual non-cavity array elements, wherein the unknown parameters are used for representing the cavity array elements in the array to be compensated;
generating the Hankel matrix corresponding to the virtual non-cavity array element according to the position of the virtual non-cavity array element in the array to be compensated and the simulation received signal corresponding to the virtual non-cavity array element;
generating an estimation matrix corresponding to the virtual non-hollow array element according to the Hankel matrix corresponding to the virtual non-hollow array element;
establishing a constraint condition according to the F norm of the difference between the estimation matrix and the Hankel matrix corresponding to the virtual non-cavity array element being less than a set value;
calculating a parameter value corresponding to the unknown parameter when the kernel norm of the estimation matrix takes the minimum value under the constraint condition;
obtaining the cavity array elements to be completed according to the parameter values;
and placing the radar at the position of the cavity array element to be completed in the array to be completed to obtain a target radar array.
In a second aspect, an embodiment of the present invention further provides an apparatus for constructing a radar array based on an undirected adjacency graph, where the apparatus includes the following components:
the device comprises a Hankel matrix generation module, a receiving module and a processing module, wherein the Hankel matrix generation module is used for generating a Hankel matrix by taking array elements in a virtual array as elements in the matrix and taking virtual receiving signals corresponding to the array elements as element values in the matrix, the array elements are used for representing a radar, and the virtual receiving signals are used for simulating receiving signals of the radar;
the undirected adjacency matrix construction module is used for setting a non-zero element value in the Hankel matrix as one and a zero element value as zero to obtain an undirected adjacency matrix;
a to-be-supplemented array construction module, configured to select the array elements from the virtual array to form a sub-array when an undirected graph corresponding to the undirected adjacency matrix has connectivity, and use the sub-array formed by the minimum array elements corresponding to the undirected graph having the non-connectivity as a to-be-supplemented array;
and the array element completion module is used for completing the cavity array elements covered by the array to be completed into non-cavity array elements according to a convex optimization algorithm to obtain a target radar array, the positions of the non-cavity array elements in the array to be completed before completion are used for placing the radars, and the positions of the cavity array elements in the array to be completed before completion are not placed with the radars.
In a third aspect, an embodiment of the present invention further provides a terminal device, where the terminal device includes a memory, a processor, and a program stored in the memory and executable on the processor and used for constructing a radar array based on a directed adjacency graph, and when the processor executes the program for constructing a radar array based on a directed adjacency graph, the steps of the method for constructing a radar array based on a directed adjacency graph are implemented.
In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium, where a program for constructing a radar array based on a non-directional adjacency graph is stored on the computer-readable storage medium, and when the program for constructing a radar array based on a non-directional adjacency graph is executed by a processor, the steps of the method for constructing a radar array based on a non-directional adjacency graph described above are implemented.
Has the advantages that: the method comprises the steps of firstly simulating a receiving signal of a radar, then constructing a virtual array by using the simulated receiving signal of the radar, then processing the virtual array to obtain an undirected adjacent matrix, and judging whether undirected graphs of the undirected adjacent matrix are communicated, wherein the undirected graphs are communicated to show that the position of a cavity array element in the virtual array can be used for placing the radar, namely the cavity array element is expanded into a non-cavity array element. If the undirected graph is connected, the position of the hollow array element in the virtual array is suitable for placing the radar. Traversing the subarrays of the virtual arrays communicated by the undirected graph, finding the subarrays (to-be-complemented arrays) which can make the undirected graph be non-communicated and consist of the minimum array elements, and then replacing the hollow array elements in the to-be-complemented arrays with the non-hollow array elements to obtain the final target radar array. The invention adopts a one-by-one traversal method to obtain the minimum radar number (the radar number corresponds to the radar aperture) required by the array when the undirected graph is not connected, and the radar array corresponding to the undirected graph can provide accurate reference data when the undirected graph is not connected. Therefore, the invention can balance the problem between the aperture of the radar array and the accuracy of the reference data provided by the radar array, namely the invention can provide the reference data with draft precision under the condition of smaller aperture of the radar array.
Drawings
FIG. 1 is an overall flow chart of the present invention;
FIG. 2 is a schematic diagram of a virtual array in an embodiment of the invention;
FIG. 3 is a co-prime array in an embodiment of the invention;
FIG. 4 is a schematic non-communicating view in an embodiment of the present invention;
FIG. 5 is a schematic illustration of communication in an embodiment of the present invention;
FIG. 6 is a diagram illustrating the physical array element positions of MCCA in an embodiment of the present invention;
FIG. 7 is a schematic diagram of a virtual array element location according to an embodiment of the present invention;
FIG. 8 is a graph comparing FFT spatial spectra in an embodiment of the present invention;
FIG. 9 is a graph illustrating deviation from a noise-free full aperture received signal in accordance with an embodiment of the present invention;
fig. 10 is a schematic block diagram of an internal structure of a terminal device according to an embodiment of the present invention.
Detailed Description
The technical scheme of the invention is clearly and completely described below by combining the embodiment and the attached drawings of the specification. 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.
Research shows that the realization of obstacle avoidance of unmanned driving, autonomous robots and unmanned aerial vehicles requires three-dimensional modeling of external complex environments. The unmanned system mainly relies on a laser radar point cloud image to construct an environment three-dimensional model. The laser radar has accurate distance measurement and high spatial resolution, and is an ideal environment sensor for automatic driving. However, as the wavelength range of the laser radar is in the nm level, when the laser radar encounters dust in the air, macromolecular aerosols such as water vapor and raindrops can generate refraction or scattering phenomena, and the laser radar cannot be used in severe weather such as rainstorm, snowstorm, haze and sand storm. In order to reduce the influence of the above natural phenomena on the lidar so that the lidar can provide necessary reference data for the unmanned technology, the prior art often adopts increasing the number of radars to increase the accuracy of the reference data provided by the array of radars. For example, an array of radars can provide reference data of the angle of arrival for the unmanned technology, and the resolution of the angle of arrival is larger when the aperture (the aperture is the size of the space occupied by the radars in the array) of the radar array is larger, so that the higher the resolution of the angle of arrival is, the safer the unmanned technology can perform. Conversely, the smaller the aperture of the radar array, the greater the resolution of the angle of arrival. And the radar array with large aperture increases the difficulty of using the radar. Therefore, it is necessary to equalize the radar array aperture and the resolution of the angle of arrival, i.e., to improve the resolution of the angle of arrival (reference data) with a limited radar array aperture.
In order to solve the technical problems, the invention provides a method and a device for constructing a radar array based on a non-directional adjacency graph, and solves the problem that in the prior art, the aperture of the radar array and the accuracy of reference data provided by the radar array are difficult to balance. During specific implementation, a virtual array is constructed by using a radar and a virtual receiving signal thereof, then the virtual array is converted into a Hankel matrix, the Hankel matrix is converted into an undirected adjacent matrix, whether undirected graphs of the undirected adjacent matrix have connectivity or not is judged, if the undirected graphs of the undirected adjacent matrix have the connectivity, array elements are selected from the virtual array to form a sub-array, the sub-array is the minimum array element required by enabling the undirected graphs corresponding to the sub-array to be non-connected, and finally, the screened hole array elements serving as an array to be complemented are complemented into non-hole array elements. The embodiment can take account of the problem between the radar array aperture and the accuracy of the reference data which can be provided by the radar array.
For example, in the virtual array shown in fig. 2, the position of each solid circle in fig. 2 represents that radar should be placed at the position, and each dotted circle represents that radar is not placed at the position temporarily. Generating a Hankel matrix according to each radar position in the virtual array and a radar signal, then converting the Hankel matrix into a non-directional adjacent matrix, judging whether the non-directional adjacent matrix has connectivity, when the non-directional adjacent matrix has connectivity, selecting radar to form a sub-array from the virtual array of FIG. 2, sequentially converting each sub-array into a corresponding Hankel matrix, a non-directional adjacent matrix and a non-directional graph, if the sub-array A formed by 20 array elements of P-12 (P represents the position), P-11, P-10, P-7 to P7, P11 and P12 in FIG. 2 is a non-connected graph, and if the number of the array elements contained in the sub-array corresponding to other non-connected graphs is more than 20, the sub-array A is a completeable array (to-be-completed array), finally completing the radar at P-11 in the sub-array A (adding a radar at the position) through a convex optimization algorithm, and finally forming 19 radar positions by the radar array of P-12, P-11, P-10, P-7 to P-7 and P12 (without P11) and the radar signal) to form a target.
Exemplary method
The method for constructing the radar array based on the undirected adjacency graph can be applied to terminal equipment, and the terminal equipment can be terminal products with a computing function, such as a computer and the like. In this embodiment, as shown in fig. 1, the method for constructing a radar array based on a non-directional adjacency graph specifically includes the following steps:
s100, generating a Hankel matrix by taking array elements in a virtual array as elements in the matrix and taking virtual receiving signals corresponding to the array elements as element values in the matrix, wherein the array elements are used for representing a radar, and the virtual receiving signals are used for simulating receiving signals of the radar.
The virtual array in this embodiment is composed of non-cavity array elements and cavity array elements, the cavity array elements are without radars, and the non-cavity array elements are with radars. And simulating a real radar array by using the virtual array, and constructing the real radar array by establishing the virtual array. The radar array constructed by the embodiment can provide high-quality reference data (DOA resolution).
In one embodiment, the method for constructing the virtual array and the simulation received signal includes the following steps S101 and S102:
and S101, respectively endowing each array element in the two preset arrays with simulation receiving signals, wherein the array elements are used for representing the radar, and the simulation receiving signals are used for simulating the receiving signals of the radar.
As shown in FIG. 3, the number of array elements of the first predetermined array is
Figure 680986DEST_PATH_IMAGE001
The number of array elements of the second predetermined array is
Figure 438727DEST_PATH_IMAGE002
Figure 945931DEST_PATH_IMAGE001
And
Figure 689897DEST_PATH_IMAGE002
are prime numbers of each other, and the spacing between the array elements in the first predetermined array is
Figure 333368DEST_PATH_IMAGE003
The spacing between the array elements in the second predetermined array is
Figure 12742DEST_PATH_IMAGE004
Figure 741663DEST_PATH_IMAGE005
The position of the array element is the position where the radar is supposed to be located, so that the radar actually needs to receive the half wavelength of the signal. The purpose of giving each array element an artificial receiving signal is to make the array element simulate a radar receiving signal.
In one embodiment, the array element position is expressed as follows:
Figure 289319DEST_PATH_IMAGE006
in this embodiment, generating the simulated received signal includes the steps of:
s1011, setting the pitching angle when the radar receives the emission signal of the signal source
Figure 52876DEST_PATH_IMAGE007
And the signal source is used for transmitting signals to the radar.
And S1012, setting noise information of the environment where the radar is located.
And S1013, generating the simulation receiving signal according to the pitching angle, the transmitting signal, the noise information and the position of the array element in the preset array.
If a radar can be from different angles during the period t
Figure 886840DEST_PATH_IMAGE007
(Angle relative to horizontal, i.e. pitch angle)
Figure 368637DEST_PATH_IMAGE008
Receiving
Figure 454404DEST_PATH_IMAGE009
A (a)
Figure 72468DEST_PATH_IMAGE009
Is relative to
Figure 326600DEST_PATH_IMAGE010
The number of the smaller number is,
Figure 30114DEST_PATH_IMAGE011
then each preset array
Figure 919573DEST_PATH_IMAGE012
Each array element in
Figure 392142DEST_PATH_IMAGE013
Simulated received signal of
Figure 833488DEST_PATH_IMAGE014
Figure 758719DEST_PATH_IMAGE015
In the formula (I), the compound is shown in the specification,
Figure 451868DEST_PATH_IMAGE016
is a pitch angle of the radar of
Figure 44523DEST_PATH_IMAGE007
When a received transmission signal (a reception signal from the perspective of radar) is sent from a source,
Figure 407503DEST_PATH_IMAGE017
is as follows
Figure 820030DEST_PATH_IMAGE012
In a predetermined array
Figure 51291DEST_PATH_IMAGE013
The position of the array elements (relative to one of the ends of the predetermined array),
Figure 498453DEST_PATH_IMAGE018
in order to emit the wavelength of the signal,
Figure 281601DEST_PATH_IMAGE019
which is noise information of the environment in which the radar is located (white gaussian noise in this embodiment).
And S102, subtracting the array elements in the two preset arrays by two to obtain virtual arrays and virtual receiving signals corresponding to the virtual arrays, wherein the subtraction is to subtract the position of the array element of one preset array in one preset array by the position of the array element of the other preset array in the other preset array, and subtract the simulated receiving signal of the array element of the one preset array by the simulated receiving signal of the array element of the other preset array.
The first positions of the two predetermined arrays in this embodiment are both zero, and the spacing between adjacent positions in each array is half a wavelength
Figure 447003DEST_PATH_IMAGE005
(in this embodiment is
Figure 216376DEST_PATH_IMAGE005
Is 1), such as the first position of the first predetermined array in FIG. 3
Figure 518044DEST_PATH_IMAGE020
The second position is
Figure 721361DEST_PATH_IMAGE021
The third position is
Figure 108480DEST_PATH_IMAGE022
Up to the fifteenth position
Figure 681544DEST_PATH_IMAGE023
Figure 103298DEST_PATH_IMAGE002
Is a mixture of a water-soluble polymer and a water-soluble polymer, and is 5,
Figure 228249DEST_PATH_IMAGE002
multiplication by
Figure 102664DEST_PATH_IMAGE005
Equal to 5, i.e. the spacing between array elements in the first predetermined array is 5, there are 3 array elements in the first predetermined array, 3x5 is equal to 15), likewise obtaining each position of the second predetermined array, the maximum position of the second predetermined array being
Figure 213840DEST_PATH_IMAGE024
A first position of a first predetermined array
Figure 490100DEST_PATH_IMAGE020
Subtracting the first position of the second predetermined array
Figure 536685DEST_PATH_IMAGE020
Equal to 0 (position 0 in fig. 2), due to the first presetArray and second predetermined array
Figure 898396DEST_PATH_IMAGE025
There is radar, so the radar should be set at the position of the virtual array 0 in fig. 2. The first predetermined array being at position 8 in FIG. 2
Figure 813262DEST_PATH_IMAGE026
By subtracting a second predetermined array
Figure 537505DEST_PATH_IMAGE025
Due to the first predetermined array
Figure 941941DEST_PATH_IMAGE026
And a second predetermined array
Figure 463052DEST_PATH_IMAGE025
There is no radar, so no radar is set at the position of the virtual array 8 in fig. 2 for the time being. The first predetermined array being at position 9 in FIG. 2
Figure 243927DEST_PATH_IMAGE027
By subtracting a second predetermined array
Figure 806364DEST_PATH_IMAGE025
Due to the first predetermined array
Figure 381702DEST_PATH_IMAGE027
There is radar, so the radar should be set at the position of the virtual array 9 in fig. 2 (the solid line in fig. 2 indicates that the radar is set, and the dotted line indicates that the radar is not set at all). Subtracting the other positions of the first preset array and the second preset array in turn can obtain whether all the positions of the virtual array in fig. 2 are provided with radars (only a part of the virtual array in the embodiment is shown in fig. 2).
After the virtual array is obtained, the virtual received signal of the cavity array element (the dotted line is the cavity array element, corresponding to the position where the radar is not set for the moment) in the virtual array is zero, and the virtual received signal of the non-cavity array element (the solid line is the non-cavity array element, corresponding to the position where the radar is set) is the result of subtracting the simulated received signals of the two corresponding array elements in the two preset arrays.
In one embodiment, generating a Hankel matrix (Hankel matrix) from a virtual array is prior art, provided that the virtual array comprises sequentially arranged array elements
Figure 390109DEST_PATH_IMAGE028
Then the hankel matrix is:
Figure 709095DEST_PATH_IMAGE029
the generation of the virtual array and the hankerr matrix in step S100 is based on the following principle:
setting a sub-array element number as
Figure 876771DEST_PATH_IMAGE001
And
Figure 623010DEST_PATH_IMAGE030
the array elements are respectively at intervals of
Figure 384293DEST_PATH_IMAGE031
And
Figure 241390DEST_PATH_IMAGE032
of a mutually prime linear array of
Figure 279885DEST_PATH_IMAGE005
Is half wavelength. All physical array element positions of the co-prime array
Figure 197025DEST_PATH_IMAGE033
Figure 180024DEST_PATH_IMAGE034
(1)
Physical array element arrangement of co-prime arrayThe schematic diagram is shown in fig. 1, two sub-arrays share a first array element and share a common array element
Figure 840813DEST_PATH_IMAGE035
A physical array element. The difference covariance array is a virtual array which performs difference extension expression on received array element signals, and the specific array element position relationship can be expressed as follows:
Figure 983081DEST_PATH_IMAGE036
(2)
from the properties of the co-prime numbers and equation (2), it can be inferred that the aperture of the difference co-array (virtual array) of the co-prime array is
Figure 71123DEST_PATH_IMAGE037
. E.g. the number of elements of a sub-array is respectively
Figure 338156DEST_PATH_IMAGE038
Figure 740319DEST_PATH_IMAGE039
A schematic diagram of the co-prime array and its difference co-array (virtual array) is shown in fig. 2. As can be seen from fig. 2, if only the continuous array element part in the middle of the virtual array is used, part of the aperture will be lost and the maximum aperture cannot be reached. And then, the holes in the virtual array are approximately completed through a matrix completion method.
Next assume that the co-prime array is from a different angle
Figure 251940DEST_PATH_IMAGE040
Receive to
Figure 510883DEST_PATH_IMAGE041
Figure 468475DEST_PATH_IMAGE041
Is a relatively small number L) of mutually incoherent narrow-band far-field target signals, and a difference co-array is generated, representing the data at the holes in the virtual array as 0, which signal isCan be expressed as:
Figure 736645DEST_PATH_IMAGE042
(3)
wherein
Figure 587927DEST_PATH_IMAGE043
Represented as a collection of received signals, and,
Figure 17771DEST_PATH_IMAGE044
to
Figure 462659DEST_PATH_IMAGE045
The individual steering vectors can be expressed as:
Figure 268941DEST_PATH_IMAGE046
(4)
will have a size of
Figure 725461DEST_PATH_IMAGE047
Received signal of
Figure 591786DEST_PATH_IMAGE048
Folded into a Hankel matrix along an inverse diagonal
Figure 258391DEST_PATH_IMAGE049
Figure 868364DEST_PATH_IMAGE050
(5)
Wherein
Figure 428658DEST_PATH_IMAGE051
Figure 669146DEST_PATH_IMAGE052
Indicating an upward integer.
S200, setting the non-zero element value in the Hankel matrix H as one and the zero element value as zero to obtain an undirected neighborConnection matrix
Figure 885364DEST_PATH_IMAGE053
(adjacency matrix of undirected graph):
Figure 345033DEST_PATH_IMAGE054
(6)
s300, when the undirected graph corresponding to the undirected adjacency matrix has connectivity, selecting the array elements from the virtual array to form a sub-array, and taking the sub-array formed by the minimum array elements corresponding to the undirected graph having the non-connectivity as an array to be supplemented.
In this embodiment, it is the prior art to generate an undirected graph from an undirected adjacency matrix.
For example, the following non-directional adjacency matrix
Figure 963096DEST_PATH_IMAGE055
And undirected adjacency matrix
Figure 108907DEST_PATH_IMAGE056
The connectivity of an undirected graph is illustrated by way of example:
Figure 812421DEST_PATH_IMAGE057
matrix array
Figure 826513DEST_PATH_IMAGE055
In (1)
Figure 299083DEST_PATH_IMAGE058
(first row, third column) mapped to the first and third points in the undirected graph (FIG. 4), and
Figure 615795DEST_PATH_IMAGE058
is a non-zero value, so mapping into an undirected graph is: there is a connection line between the first point and the third point, since there is a connection point in fig. 4 that is not connected to any one point, there is a connection line between the first point and the third pointThe non-directional adjacency matrix
Figure 150812DEST_PATH_IMAGE055
Undirected graph of (a) is unconnected.
By way of example, FIG. 4 is
Figure 578382DEST_PATH_IMAGE055
Undirected graph of (2), generating
Figure 171038DEST_PATH_IMAGE055
Is based on the following principle:
an undirected graph only considers whether the elements in the undirected adjacency matrix that are off-diagonal are zero, e.g.
Figure 517706DEST_PATH_IMAGE055
In (1)
Figure 195812DEST_PATH_IMAGE058
(first row and third column) are mapped to the first and third points in the undirected graph, since
Figure 427073DEST_PATH_IMAGE058
Is non-zero and therefore there is a line between the first point and the third point. While
Figure 874235DEST_PATH_IMAGE055
The first and second points in the undirected graph mapped by the elements of the first row and the second column in
Figure 906650DEST_PATH_IMAGE055
The elements in the first row and the second column are zero, so that no connecting line exists between the first point and the second point in the undirected graph, and the same can be seen
Figure 806473DEST_PATH_IMAGE055
There is no connection line between the second point and the third point in the undirected graph. Therefore, the temperature of the molten metal is controlled,
Figure 575846DEST_PATH_IMAGE055
undirected graphs are unconnected.
FIG. 5 is
Figure 143094DEST_PATH_IMAGE059
Generating an undirected graph of
Figure 97143DEST_PATH_IMAGE059
Is based on the following principle:
matrix array
Figure 687525DEST_PATH_IMAGE059
In (1)
Figure 322905DEST_PATH_IMAGE060
(first row, second column) are mapped to the first and second points in the undirected graph (FIG. 5), and
Figure 557709DEST_PATH_IMAGE060
is a non-zero value, so mapping into an undirected graph is: there is a connection between the first point and the second point, and similarly, a non-zero value
Figure 620343DEST_PATH_IMAGE061
(second row, third column) maps the second and third points in the undirected graph (fig. 5), so the mapping into the undirected graph is: there is a connection line between the second point and the third point, and in summary, there is no directional adjacency matrix
Figure 229178DEST_PATH_IMAGE059
Undirected graph of (1) is connected.
In one embodiment, when the undirected graph corresponding to the undirected adjacency matrix does not have connectivity, the virtual array at this time is the array to be complemented (the complemented array).
In another embodiment, when the undirected graph corresponding to the undirected adjacency matrix has connectivity, an array element needs to be selected from the virtual array to form each sub-array, and then a sub-array which can be used as an array to be complemented is selected from each sub-array. The embodiment includes the following steps S301, S302, S303, and S304:
s301, the array elements at the head end and the tail end of the virtual array form a first sub-array in the sub-arrays.
S302, the array elements except the array elements at the head end and the tail end in the virtual array form a second sub array in the sub arrays, and the second sub array and the first sub array are relatively prime arrays.
S303, if the undirected graph corresponding to the virtual sub-array formed by the first sub-array and the second sub-array has connectivity, adding the array elements of the virtual array to the second sub-array and the first sub-array one by one until the undirected graph corresponding to the virtual sub-array formed by the second sub-array and the first sub-array after the addition of the array elements has non-connectivity, and using the second sub-array and the first sub-array after the addition of the array elements as a target sub-array MCCA (minimum complementary co-prime array).
And S304, obtaining an array to be compensated according to the target subarray.
The detailed process of S301 to S304 is illustrated:
if the undirected graph corresponding to the virtual array in FIG. 2 is not connected, p is added 12 And p -12 Is (p) 12 And p -12 At the head and tail ends) to form a first sub-array, and then randomly selecting p 1 、p 9 、p 10 And the array elements form a second sub-array, whether undirected graphs corresponding to a new virtual array formed by the first sub-array and the second sub-array are communicated or not is calculated, traversal is completed once at the moment, and the number of the array elements contained in the first sub-array and the second sub-array are mutually equal in traversal each time. If the undirected graph obtained by the traversal is communicated, the next traversal is carried out, an array element is added to the first subarray and the second subarray traversed last time in the next traversal, whether the undirected graph corresponding to the virtual array generated by the first subarray and the second subarray after the array element is added has connectivity is judged, and the undirected graph is generated until the undirected graph is formed by the first subarray and the second subarray after the array element is added one by oneThe undirected graph corresponding to the pseudo array does not have connectivity, and the first sub-array and the second sub-array after the array elements are added are finally needed target sub-arrays MCCA. And generating a final array to be compensated by the two target sub-arrays. (two target sub-arrays are differenced to obtain an array to be complemented, the form of the array to be complemented is the same as that of the virtual array in the figure 2, and the array comprises a cavity array element and a non-cavity array element).
The generation of the two target sub-arrays MCCA in steps S301 to S304 is based on the following principle:
in order not to destroy the maximum aperture of the difference covariance array (virtual array), the first and last elements of the physical array are fixed. Then, the number of antennas (array elements) is from small to large (because the antennas at two fixed ends are included, the number of the MCCA antennas is increased from 3 upwards), the other antennas except the antennas at two ends are subjected to position combination traversal to generate sub-arrays of a co-prime array, the connectivity of an undirected graph corresponding to the co-prime array is judged, and the minimum number of antennas is obtained and conforms to the non-connectivity MCCA.
And S400, according to a convex optimization algorithm, completing the hollow array elements covered by the array to be completed into non-hollow array elements to obtain a target radar array, wherein the positions of the non-hollow array elements in the array to be completed before completing are used for placing the radars, and the positions of the hollow array elements in the array to be completed before completing are not used for placing the radars.
Step S300 obtains a to-be-supplemented array formed by two target sub-arrays MCCA, and when the to-be-supplemented array includes a cavity array element, it needs to be considered whether to supplement each cavity array element included in the array, that is, whether to set a radar at a position where the cavity array element is located, so that the supplemented radar array can provide high-quality reference data (such as a wave arrival angle) for the unmanned technology. Step S400 includes the following steps:
s401, adding unknown parameters to the non-cavity array elements in the array to be compensated to obtain virtual non-cavity array elements, wherein the unknown parameters are used for representing the cavity array elements in the array to be compensated.
Because it is not known which cavity array element and how many cavities need to be completedThe hole array elements can enable the radar array corresponding to the array after completion to provide high-quality wave arrival angles, and therefore the hole array elements needing completion are used as unknown parameters to be solved. The method comprises the steps of firstly assuming a hollow array element as a non-hollow array element, and forming a set serving as an observation element by the assumed non-hollow array element and an original real non-hollow array element of an array to be compensated
Figure 605933DEST_PATH_IMAGE062
The virtual non-hole array element.
S402, generating the Hankel matrix H corresponding to the virtual non-hollow array element according to the position of the virtual non-hollow array element in the array to be compensated and the simulation received signal corresponding to the virtual non-hollow array element.
In this embodiment, the hankel matrix H corresponding to the virtual non-cavity array element is calculated by using formula (5).
And S403, generating an estimation matrix X corresponding to the virtual non-hollow array element according to the Hankel matrix corresponding to the virtual non-hollow array element.
S404, establishing a constraint condition by using the F norm of the difference between the estimation matrix and the Hankel matrix corresponding to the virtual non-cavity array element to be less than a set value:
Figure 944511DEST_PATH_IMAGE063
Figure 178046DEST_PATH_IMAGE064
represents the F norm (Frobenius-norm).
Set value in the present embodiment
Figure 477440DEST_PATH_IMAGE065
Noise information characterizing the environment in which the radar is located, i.e. the greater the noise,
Figure 454623DEST_PATH_IMAGE065
the larger and the opposite
Figure 896975DEST_PATH_IMAGE065
The smaller.
In one embodiment, the error parameter
Figure 301412DEST_PATH_IMAGE065
Related to noise power and array sparsity.
S405, calculating the parameter value corresponding to the unknown parameter when the kernel norm of the estimation matrix X is the minimum value under the constraint condition.
And S406, obtaining the cavity array element to be completed according to the parameter value.
S407, placing the radar at the position of the cavity array element to be completed in the array to be completed to obtain a target radar array.
Namely, it is
Figure 88102DEST_PATH_IMAGE066
When the virtual non-hollow array elements are established, after the specific virtual non-hollow array elements are obtained, the known non-hollow array elements which are the hollow array elements corresponding to the unknown parameters are removed, and the cavity array elements which need to be complemented into the non-hollow array elements are solved.
In one embodiment, the constraints are derived using the following principles
Figure 868976DEST_PATH_IMAGE067
And an objective function
Figure 916567DEST_PATH_IMAGE066
Hankel matrix formed by received signals
Figure 491904DEST_PATH_IMAGE068
Is a rank of
Figure 765891DEST_PATH_IMAGE069
According to the theory of completion by matrix:
Figure 84877DEST_PATH_IMAGE070
wherein
Figure 3285DEST_PATH_IMAGE071
Which represents the operation of determining the rank of the matrix,
Figure 749524DEST_PATH_IMAGE072
is composed of
Figure 510807DEST_PATH_IMAGE068
The total set of (a) and (b),
Figure 367905DEST_PATH_IMAGE073
is an approximation of the generation
Figure 390087DEST_PATH_IMAGE072
Of the estimation matrix, and
Figure 572807DEST_PATH_IMAGE074
for observing sets of elements, i.e.
Figure 555806DEST_PATH_IMAGE075
Is composed of
Figure 793759DEST_PATH_IMAGE073
The observation projection matrix of (2):
Figure 873710DEST_PATH_IMAGE076
due to the fact that
Figure 899435DEST_PATH_IMAGE077
The minimization matrix rank of (c) belongs to the non-convex problem of NP-Hard, so the optimization problem is relaxed:
Figure 432047DEST_PATH_IMAGE078
wherein
Figure 958844DEST_PATH_IMAGE079
Represents the nuclear norm (nuclear-norm). The matrix is destroyed to some extent due to the presence of noise
Figure 893302DEST_PATH_IMAGE080
Low rank characteristics, affecting the outcome of the optimization. In order to improve the optimization effect and remove noise, the distance with error parameters is characterized by the constraint conditions:
Figure 89928DEST_PATH_IMAGE081
the number of sub-arrays is
Figure 188465DEST_PATH_IMAGE082
Figure 191056DEST_PATH_IMAGE083
The target radar array constructed by the invention can provide high-quality arrival angles, namely the arrival angles have higher resolution.
The above co-prime array and the physical and virtual array element positions of 17 MCCAs are shown in fig. 6 and 7, wherein the 0 th line is a co-prime array. The number of virtual array elements of the co-prime array in fig. 7 is 45, and the number of virtual array elements of the MCCA in the other 17 rows is 23 except for row 2, which is 19. Virtual array
Figure 917704DEST_PATH_IMAGE084
The sparsity of (a) is defined as:
Figure 613127DEST_PATH_IMAGE085
wherein is
Figure 448228DEST_PATH_IMAGE086
The number of the arrays is represented by,
Figure 926614DEST_PATH_IMAGE087
representing a completion operation. The higher the sparsity, the more the number of array elements of the array is, and the more sparse the array is. The sparsity of the coprime array is about 0.26, while the sparsity of the second row of MCCA is about 0.69 and the sparsity of the remaining MCCAs is about 0.63.
Next, a simulation experiment link is carried out, and according to the power of noise and the sparsity of the array, the error parameters of the co-prime array low-rank matrix completion optimization algorithm are compensated
Figure 835664DEST_PATH_IMAGE088
Set to 0.8 and MCCA to 0.2. Taking the second row of MCCA with the highest sparsity and all virtual arrays (Coprime Array, CA) of the initial Coprime Array, continuous virtual apertures (CA-continuous) of the Coprime Array and Full-Aperture arrays (Full-Aperture) with noise as experimental objects, and respectively complementing the MCCA and the CA into MCCA-completed and CA-completed Full-Aperture arrays by using a low-rank matrix method.
Experiment one: fast Fourier Transform (FFT) spatial spectrum
Setting a single snapshot of dual target signals with array incidence angles of
Figure 747994DEST_PATH_IMAGE089
And
Figure 742495DEST_PATH_IMAGE090
the Signal-to-Noise Ratio (SNR) was set to 20dB. As shown in FIG. 8, MCCA-completed, CA-completed and Full-Aperture can form peaks accurately, and the other methods have different degrees of angle ambiguity.
Experiment two: deviation from noiseless full aperture received signal
MCCA-completed and CA-completed, full-Aperture and noiseless Full Aperture arrays
The variation trend of the deviation value of the received signal with SNR is shown in fig. 9, where the deviation is defined as:
Figure 24572DEST_PATH_IMAGE091
wherein
Figure 788128DEST_PATH_IMAGE092
Represents the first
Figure 153251DEST_PATH_IMAGE093
The received signal of the sub-monte carlo experiment,
Figure 307151DEST_PATH_IMAGE094
then it is a noise-free full-aperture received signal,
Figure 533864DEST_PATH_IMAGE095
representing a 1-norm operation. The experimental result is shown in fig. 9, the MCCA and CA after completion can eliminate partial noise, and the noise reduction effect of MCCA is better than that of CA.
In summary, the method includes firstly simulating a receiving signal of the radar, then constructing a virtual array by the simulated receiving signal of the radar, then processing the virtual array to obtain an undirected adjacent matrix, and judging whether undirected graphs of the undirected adjacent matrix are communicated or not, wherein the undirected graphs are communicated to show that the position of a cavity array element in the virtual array can be used for placing the radar, namely the cavity array element is expanded into a non-cavity array element. If the undirected graph is connected, the position of the hollow array element in the virtual array is suitable for placing the radar. Traversing the subarrays of the virtual arrays communicated by the undirected graph, finding the subarrays (to-be-complemented arrays) which can make the undirected graph be non-communicated and consist of the minimum array elements, and then replacing the hollow array elements in the to-be-complemented arrays with the non-hollow array elements to obtain the final target radar array. The invention adopts a one-by-one traversal method to obtain the minimum radar number (the radar number corresponds to the radar aperture) required by the array when the undirected graph is not connected, and the radar array corresponding to the undirected graph can provide accurate reference data when the undirected graph is not connected. Therefore, the invention can balance the problem between the aperture of the radar array and the accuracy of the reference data provided by the radar array, namely the invention can provide the reference data with draft precision under the condition of smaller aperture of the radar array.
In addition, the invention takes a vehicle-mounted millimeter wave radar as a research platform, and designs a minimum complementable co-prime linear array according to the connectivity of a covariance matrix of a received signal. Firstly, setting an initial co-prime linear array and generating a covariance matrix of a received signal; then, traversing and judging the subarray of the array through a discriminator, and screening a completeable subarray with the least array elements as MCCA; and finally, complementing the sparse array obtained in the last step into a uniform array by a convex optimization method, thereby realizing the aperture expansion of the array without the angle ambiguity.
Exemplary devices
The embodiment also provides a device for constructing a radar array based on the undirected adjacency graph, which comprises the following components:
the device comprises a Hankel matrix generation module, a receiving module and a processing module, wherein the Hankel matrix generation module is used for generating a Hankel matrix by taking array elements in a virtual array as elements in the matrix and taking virtual receiving signals corresponding to the array elements as element values in the matrix, the array elements are used for representing a radar, and the virtual receiving signals are used for simulating receiving signals of the radar;
the undirected adjacency matrix construction module is used for setting a non-zero element value in the Hankel matrix as one and a zero element value as zero to obtain an undirected adjacency matrix;
a to-be-supplemented array construction module, configured to select the array elements from the virtual array to form a sub-array when an undirected graph corresponding to the undirected adjacency matrix has connectivity, and use the sub-array formed by the minimum array elements corresponding to the undirected graph having the non-connectivity as a to-be-supplemented array;
and the array element completion module is used for completing the cavity array elements covered by the array to be completed into non-cavity array elements according to a convex optimization algorithm to obtain a target radar array, the positions of the non-cavity array elements in the array to be completed before completion are used for placing the radars, and the positions of the cavity array elements in the array to be completed before completion are not placed with the radars.
Based on the above embodiments, the present invention further provides a terminal device, and a schematic block diagram thereof may be as shown in fig. 10. The terminal equipment comprises a processor, a memory, a network interface, a display screen and a temperature sensor which are connected through a system bus. Wherein the processor of the terminal device is configured to provide computing and control capabilities. The memory of the terminal equipment comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The network interface of the terminal device is used for connecting and communicating with an external terminal through a network. The computer program is executed by a processor to implement a method of constructing a radar array based on an undirected adjacency graph. The display screen of the terminal equipment can be a liquid crystal display screen or an electronic ink display screen, and the temperature sensor of the terminal equipment is arranged in the terminal equipment in advance and used for detecting the operating temperature of the internal equipment.
It will be understood by those skilled in the art that the block diagram of fig. 10 is only a block diagram of a part of the structure related to the solution of the present invention, and does not constitute a limitation to the terminal device to which the solution of the present invention is applied, and a specific terminal device may include more or less components than those shown in the figure, or may combine some components, or have different arrangements of components.
In one embodiment, a terminal device is provided, where the terminal device includes a memory, a processor, and a program stored in the memory and executable on the processor for constructing a radar array based on a directed adjacency graph, and when the processor executes the program for constructing a radar array based on a directed adjacency graph, the following operation instructions are implemented:
taking an array element in a virtual array as an element in a matrix, and taking a virtual received signal corresponding to the array element as an element value in the matrix to generate a Hankel matrix, wherein the array element is used for representing a radar, and the virtual received signal is used for simulating a received signal of the radar;
setting a non-zero element value in the Hankel matrix as one and a zero element value as zero to obtain an undirected adjacency matrix;
when an undirected graph corresponding to the undirected adjacency matrix has connectivity, selecting the array elements from the virtual array to form a sub-array, and taking the sub-array formed by the minimum array elements corresponding to the undirected graph having the non-connectivity as an array to be supplemented;
and according to a convex optimization algorithm, complementing the cavity array elements covered by the array to be complemented into non-cavity array elements to obtain a target radar array, wherein the positions of the non-cavity array elements in the array to be complemented are used for placing the radar before complementing, and the positions of the cavity array elements in the array to be complemented are not placed into the radar before complementing.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, databases or other media used in the embodiments provided herein may include non-volatile and/or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically Programmable ROM (EPROM), electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous Link DRAM (SLDRAM), rambus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (10)

1. A method for constructing a radar array based on an undirected adjacency graph is characterized by comprising the following steps:
generating a Hankel matrix by taking array elements in a virtual array as elements in the matrix and taking virtual receiving signals corresponding to the array elements as element values in the matrix, wherein the array elements are used for representing a radar, and the virtual receiving signals are used for simulating receiving signals of the radar;
setting a non-zero element value in the Hankel matrix as one and a zero element value as zero to obtain an undirected adjacency matrix;
when an undirected graph corresponding to the undirected adjacency matrix has connectivity, selecting the array elements from the virtual array to form a sub-array, and taking the sub-array formed by the minimum array elements corresponding to the undirected graph having the non-connectivity as an array to be supplemented;
and according to a convex optimization algorithm, filling the cavity array elements covered by the array to be filled with non-cavity array elements to obtain a target radar array, wherein the positions of the non-cavity array elements in the array to be filled before filling are used for placing the radar, and the positions of the cavity array elements in the array to be filled before filling are not placed with the radar.
2. The method for constructing a radar array based on a non-directional adjacency graph as claimed in claim 1, wherein the virtual array and the method for generating the virtual received signal comprise:
respectively endowing each array element in two preset arrays with simulation received signals, wherein the array elements are used for representing a radar, and the simulation received signals are used for simulating received signals of the radar;
and subtracting the array elements in the two preset arrays in pairs to obtain virtual arrays and virtual receiving signals corresponding to the virtual arrays, wherein the subtraction is to subtract the position of the array element of one of the preset arrays in one of the preset arrays by the position of the array element of the other preset array in the other preset array, and subtract the simulated receiving signals of the array elements of the one of the preset arrays by the simulated receiving signals of the array elements of the other preset array.
3. The method as claimed in claim 2, wherein the total number of the array elements covered by one of the two predetermined arrays and the total number of the array elements covered by the other of the two predetermined arrays are mutually prime numbers, the array element pitch of one of the predetermined arrays is equal to the total number of the array elements of the other of the predetermined arrays multiplied by a wavelength, and the array element pitch of the other of the predetermined arrays is equal to the total number of the array elements of one of the predetermined arrays multiplied by a wavelength, and the wavelength is a wavelength of a signal to be received by the radar.
4. The method for constructing a radar array based on an undirected adjacency graph according to claim 2, wherein the method for generating the simulated received signals comprises:
setting a pitching angle when the radar receives a transmission signal of a signal source, wherein the signal source is used for transmitting the signal to the radar;
setting noise information of the environment where the radar is located;
and generating the simulation receiving signal according to the pitching angle, the transmitting signal, the noise information and the position of the array element in the preset array.
5. The method according to claim 1, wherein when an undirected graph corresponding to the undirected adjacency matrix has connectivity, selecting the array elements from the virtual array to form a sub-array, and using the sub-array formed by the fewest array elements corresponding to the undirected graph having non-connectivity as an array to be complemented, comprises:
forming sub-arrays by using the array elements at the head end and the tail end of the virtual array, and if the undirected graphs corresponding to the sub-arrays have connectivity, adding the array elements except the array elements at the head end and the tail end to the sub-arrays one by one until the undirected graphs corresponding to the sub-arrays after the array elements are added have non-connectivity to obtain target sub-arrays;
and obtaining an array to be compensated according to the target subarray.
6. The method for constructing a radar array based on the undirected adjacency graph as claimed in claim 5, wherein the forming of the array elements at the head and the tail ends of the virtual array into sub-arrays, and if the undirected graphs corresponding to the sub-arrays have connectivity, adding the array elements except the array elements at the head and the tail ends to the sub-arrays one by one until the undirected graphs corresponding to the sub-arrays after the addition of the array elements have non-connectivity, to obtain a target sub-array, comprises:
forming a first sub-array in the sub-arrays by using the array elements at the head end and the tail end of the virtual array;
forming the array elements except the array elements at the head end and the tail end in the virtual array into a second sub array in the sub arrays, wherein the second sub array and the first sub array are relatively prime arrays;
if the undirected graphs corresponding to the virtual sub-arrays formed by the first sub-array and the second sub-array have connectivity, the array elements of the virtual arrays are respectively added to the second sub-array and the first sub-array one by one until the undirected graphs corresponding to the virtual sub-arrays formed by the second sub-array and the first sub-array after the array elements are added have non-connectivity, and the second sub-array and the first sub-array after the array elements are added are taken as target sub-arrays.
7. The method according to claim 4, wherein the method for constructing a radar array based on the undirected adjacency graph includes the following steps of completing the hole array elements covered by the array to be completed into non-hole array elements according to a convex optimization algorithm to obtain a target radar array, wherein the positions of the non-hole array elements in the array to be completed before completing are used for placing the radars, and the positions of the hole array elements in the array to be completed before completing are not used for placing the radars, including:
adding unknown parameters to the non-cavity array elements in the array to be compensated to obtain virtual non-cavity array elements, wherein the unknown parameters are used for representing the cavity array elements in the array to be compensated;
generating the Hankel matrix corresponding to the virtual non-cavity array element according to the position of the virtual non-cavity array element in the array to be compensated and the simulation received signal corresponding to the virtual non-cavity array element;
generating an estimation matrix corresponding to the virtual non-hollow array element according to the Hankel matrix corresponding to the virtual non-hollow array element;
establishing a constraint condition according to the F norm of the difference between the estimation matrix and the Hankel matrix corresponding to the virtual non-cavity array element, which is smaller than a set value;
calculating a parameter value corresponding to the unknown parameter when the kernel norm of the estimation matrix takes the minimum value under the constraint condition;
obtaining the cavity array elements to be completed according to the parameter values;
and placing the radar at the position of the cavity array element to be completed in the array to be completed to obtain a target radar array.
8. An apparatus for constructing a radar array based on an undirected adjacency graph, the apparatus comprising:
the device comprises a Hankel matrix generation module, a radar detection module and a Hankel matrix generation module, wherein the Hankel matrix generation module is used for generating a Hankel matrix by taking array elements in a virtual array as elements in the matrix and taking virtual receiving signals corresponding to the array elements as element values in the matrix, the array elements are used for representing the radar, and the virtual receiving signals are used for simulating receiving signals of the radar;
the undirected adjacency matrix construction module is used for setting a non-zero element value in the Hankel matrix as one and a zero element value as zero to obtain an undirected adjacency matrix;
a to-be-supplemented array construction module, configured to select the array elements from the virtual array to form a sub-array when an undirected graph corresponding to the undirected adjacency matrix has connectivity, and use the sub-array formed by the minimum array elements corresponding to the undirected graph having the non-connectivity as a to-be-supplemented array;
and the array element completion module is used for completing the cavity array elements covered by the array to be completed into non-cavity array elements according to a convex optimization algorithm to obtain a target radar array, the positions of the non-cavity array elements in the array to be completed are used for placing the radar before completion, and the positions of the cavity array elements in the array to be completed are not placed in the radar before completion.
9. A terminal device, characterized in that the terminal device comprises a memory, a processor and a program stored in the memory and executable on the processor for constructing a radar array based on a directed adjacency graph, and the processor implements the steps of the method for constructing a radar array based on a directed adjacency graph according to any one of claims 1 to 7 when executing the program for constructing a radar array based on a directed adjacency graph.
10. A computer-readable storage medium, wherein the computer-readable storage medium has stored thereon a program for constructing a radar array based on a undirected adjacency graph, which program, when executed by a processor, carries out the steps of the method for constructing a radar array based on a undirected adjacency graph as claimed in any one of claims 1 to 7.
CN202211298872.1A 2022-10-24 2022-10-24 Method and device for constructing radar array based on undirected adjacency graph Active CN115421157B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211298872.1A CN115421157B (en) 2022-10-24 2022-10-24 Method and device for constructing radar array based on undirected adjacency graph

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211298872.1A CN115421157B (en) 2022-10-24 2022-10-24 Method and device for constructing radar array based on undirected adjacency graph

Publications (2)

Publication Number Publication Date
CN115421157A true CN115421157A (en) 2022-12-02
CN115421157B CN115421157B (en) 2023-03-24

Family

ID=84207616

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211298872.1A Active CN115421157B (en) 2022-10-24 2022-10-24 Method and device for constructing radar array based on undirected adjacency graph

Country Status (1)

Country Link
CN (1) CN115421157B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116698016A (en) * 2023-08-04 2023-09-05 西安交通大学 System and method for constructing miniature radar array in complex narrow space
CN117473687A (en) * 2023-12-28 2024-01-30 深圳大学 Construction method and equipment of redundancy-removing mutual radar array based on single-bit quantization

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108828551A (en) * 2018-08-28 2018-11-16 中国人民解放军空军工程大学 A kind of compressed sensing based flexible MIMO radar compound target DOA estimation method
CN108983168A (en) * 2018-04-27 2018-12-11 常熟理工学院 Compressed sensing MIMO radar object detection method based on aperture completion
CN109471082A (en) * 2018-11-07 2019-03-15 南京信息工程大学 Array element defect MIMO radar angle estimating method based on signal subspace reconstruct
CN109471086A (en) * 2018-10-18 2019-03-15 浙江大学 Relatively prime MIMO radar Wave arrival direction estimating method based on more sampling snap sum aggregate array signal discrete Fourier transforms
CN109471087A (en) * 2018-10-18 2019-03-15 浙江大学 Wave arrival direction estimating method based on relatively prime MIMO radar difference set sum aggregate signal Fast Fourier Transform (FFT)
CN109782243A (en) * 2018-12-26 2019-05-21 南京信息工程大学 Incomplete data MIMO radar angle estimating method based on block Hankel matrix filling
CN113219400A (en) * 2020-08-05 2021-08-06 哈尔滨工业大学(威海) Underdetermined direction-of-arrival estimation method based on void-free co-prime array structure
US20210311182A1 (en) * 2020-04-01 2021-10-07 The Board Of Trustees Of The University Of Alabama Sparse linear array approach in automotive radars using matrix completion
CN114740471A (en) * 2022-03-24 2022-07-12 中国人民解放军国防科技大学 Array radar foresight imaging method and device based on echo signal completion

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108983168A (en) * 2018-04-27 2018-12-11 常熟理工学院 Compressed sensing MIMO radar object detection method based on aperture completion
CN108828551A (en) * 2018-08-28 2018-11-16 中国人民解放军空军工程大学 A kind of compressed sensing based flexible MIMO radar compound target DOA estimation method
CN109471086A (en) * 2018-10-18 2019-03-15 浙江大学 Relatively prime MIMO radar Wave arrival direction estimating method based on more sampling snap sum aggregate array signal discrete Fourier transforms
CN109471087A (en) * 2018-10-18 2019-03-15 浙江大学 Wave arrival direction estimating method based on relatively prime MIMO radar difference set sum aggregate signal Fast Fourier Transform (FFT)
CN109471082A (en) * 2018-11-07 2019-03-15 南京信息工程大学 Array element defect MIMO radar angle estimating method based on signal subspace reconstruct
CN109782243A (en) * 2018-12-26 2019-05-21 南京信息工程大学 Incomplete data MIMO radar angle estimating method based on block Hankel matrix filling
US20210311182A1 (en) * 2020-04-01 2021-10-07 The Board Of Trustees Of The University Of Alabama Sparse linear array approach in automotive radars using matrix completion
CN113219400A (en) * 2020-08-05 2021-08-06 哈尔滨工业大学(威海) Underdetermined direction-of-arrival estimation method based on void-free co-prime array structure
CN114740471A (en) * 2022-03-24 2022-07-12 中国人民解放军国防科技大学 Array radar foresight imaging method and device based on echo signal completion

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JING WANG等: "A_New_Generalized_Khatri-Rao_Products_Approach_for_Correlated_Source_Number_Estimation", 《IEEE》 *
张伟科: "高分辨阵列信号DOA估计关键技术研究", 《中国博士学位论文全文数据库信息科技辑》 *
韦娟等: "基于互质阵虚拟阵列空间平滑的相干信号DOA估计方法", 《***工程与电子技术》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116698016A (en) * 2023-08-04 2023-09-05 西安交通大学 System and method for constructing miniature radar array in complex narrow space
CN116698016B (en) * 2023-08-04 2023-10-20 西安交通大学 System and method for constructing miniature radar array in complex narrow space
CN117473687A (en) * 2023-12-28 2024-01-30 深圳大学 Construction method and equipment of redundancy-removing mutual radar array based on single-bit quantization
CN117473687B (en) * 2023-12-28 2024-05-10 深圳大学 Construction method and equipment of redundancy-removing mutual radar array based on single-bit quantization

Also Published As

Publication number Publication date
CN115421157B (en) 2023-03-24

Similar Documents

Publication Publication Date Title
CN115421157B (en) Method and device for constructing radar array based on undirected adjacency graph
CN109100718B (en) Sparse aperture ISAR self-focusing and transverse calibration method based on Bayesian learning
JP6660857B2 (en) Method of generating image of region of interest using radar system
CN107340434B (en) Satellite-borne synthetic aperture radar antenna direction picture capturing method based on internal calibration data
CN108562897A (en) A kind of sparse imaging method of structure and device of MIMO through-wall radars
CN111190151B (en) System parameter design and emission power optimization method of multi-mode small satellite SAR in scanning mode
CN110609282A (en) Three-dimensional target imaging method and device
CN115291169A (en) Sound source imaging method, system, equipment and storage medium
CN113630197B (en) Antenna phase adjustment method and device, storage medium and electronic equipment
CN115184916A (en) Sea surface wind speed joint inversion method, device, medium and computing equipment
CN111965645B (en) Multi-baseline vegetation height inversion method and device considering geometric constraints
CN113835063A (en) Unmanned aerial vehicle array amplitude and phase error and signal DOA joint estimation method
CN114994674B (en) Intelligent microwave staring correlated imaging method and equipment and storage medium
CN106772368B (en) The super-resolution three-D imaging method of more random frequency radar arrays
CN108445458A (en) A kind of synthetic aperture radar orbit error removing method and device
CN111915570A (en) Atmospheric delay estimation method based on back propagation neural network
CN111259511A (en) Method and device for processing sea surface flow field
CN106918809B (en) Method for simulating original echo time domain of fast interferometric synthetic aperture sonar
CN114740470A (en) Microwave wavefront modulation foresight imaging method and device based on attribute scattering model
CN111539106B (en) Electronic reconnaissance equipment output data simulation method based on parametric modeling data
CN111123250B (en) Pulse Doppler radar based on pattern search algorithm and beam forming method
CN117473687B (en) Construction method and equipment of redundancy-removing mutual radar array based on single-bit quantization
CN112859074A (en) Multi-band multi-view ISAR fusion imaging method
KR102589829B1 (en) Validation apparatus and method for radar training device
CN113805147B (en) Vehicle-mounted radar level measurement angle self-correction method, device, medium and equipment

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