CN111665501B - MIMO radar two-dimensional imaging method based on improved CBP - Google Patents
MIMO radar two-dimensional imaging method based on improved CBP Download PDFInfo
- Publication number
- CN111665501B CN111665501B CN202010603678.4A CN202010603678A CN111665501B CN 111665501 B CN111665501 B CN 111665501B CN 202010603678 A CN202010603678 A CN 202010603678A CN 111665501 B CN111665501 B CN 111665501B
- Authority
- CN
- China
- Prior art keywords
- imaging
- mimo radar
- algorithm
- area
- cbp
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses an MIMO radar two-dimensional imaging method based on improved CBPs, and relates to the technical field of millimeter wave MIMO radar imaging. The method comprises the following steps: establishing a near-field millimeter wave MIMO radar echo model; designing a group of matched filter functions to compensate phase errors introduced by an equivalent phase center principle; dividing an imaging area, detecting a target, and mapping the target to an imaging unit; extracting pixel units according to the probability to form an imaging pixel unit set; and sequentially imaging the sampling set pixel units by using a CBP algorithm, and reconstructing the rest pixel units according to a non-precise augmented Lagrange IALM algorithm to finish the two-dimensional imaging of the MIMO radar. The method greatly shortens the imaging time of the CBP algorithm, and can be widely applied to MIMO radar near-field imaging.
Description
Technical Field
The invention relates to the technical field of millimeter wave MIMO radar imaging, in particular to an MIMO radar two-dimensional imaging method based on improved CBP.
Background
The millimeter wave radar has the advantages of weather resistance and all-weather working capability, shows a series of advantages in the application of the millimeter wave in the fields of engineering and military, and plays a very important role in electromagnetic wave resources. The millimeter wave technology is mature day by day and is widely applied to the fields of communication radar, remote sensing, accurate guidance and automobile radar. The millimeter wave radar has stronger anti-interference capability due to the selectivity of the narrower bandwidth of the millimeter wave radar to atmospheric propagation and the larger Doppler frequency. With the continuous progress of unmanned technology, the performance requirements on radar are gradually improved, the military field comprises military reconnaissance, target detection/tracking and the like, the functions of distance measurement, speed measurement, angle measurement and the like of targets can be realized, the application value of high-resolution two-dimensional and three-dimensional imaging is particularly huge, large ground and sea surface scene information under a working environment can be provided, and the detection and identification precision of the military targets is greatly improved.
Disclosure of Invention
The invention aims to solve the technical problem of how to provide a two-dimensional imaging method of an MIMO radar, which can improve the imaging speed.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows: a MIMO radar two-dimensional imaging method based on improved CBP is characterized by comprising the following steps:
establishing a near-field millimeter wave MIMO radar echo model;
designing a group of matched filter functions to compensate phase errors introduced by an equivalent phase center principle;
dividing an imaging area, detecting a target, and mapping the target to an imaging unit;
extracting pixel units according to the probability to form an imaging pixel unit set;
and sequentially imaging the sampling set pixel units by using a CBP algorithm, and reconstructing the rest pixel units according to a non-precise augmented Lagrange IALM algorithm to finish the two-dimensional imaging of the MIMO radar.
The technical scheme is that the method for establishing the near-field millimeter wave MIMO radar echo model comprises the following steps:
let At(xt,yt) And Ar(xr,yr) Respectively, coordinates of a pair of transmitting and receiving antennas, for an object M (x) in the object arean,yn) At respective distances R from the transmitting antenna and the receiving antennat,n、Rr,nWherein:
the emission signals are:
s(t)=exp[-j2πft] (3)
wherein f ═ f0+(n-1)Δf,f0Is the initial frequency;
the FMCW radar mixing signal is:
Sn(xt,yt;xr,yr;Kω)=σ(xn,yn)exp[-jKω(Rt,n+Rr,n)] (4)
wherein KwSince 2 pi f/c is the spatial frequency, the signals received by the transmitting and receiving antenna in the target area are:
the further technical scheme is that a group of matched filter functions is designed, and a method for compensating phase errors introduced by an equivalent phase center principle comprises the following steps:
is equivalent to a self-transmitting and self-receiving antenna, and the position of the antenna is Ae(xe,ye) And then:
the compensation function is expressed as:
the compensated signals are:
the further technical solution is that the method for dividing the imaging area and detecting the target and mapping to the imaging unit is as follows:
establishing MIMO radar space coordinates, taking an MIMO radar array as an x axis, taking the array direction as a y axis vertically, and dividing an imaging area, wherein the MIMO radar is provided with Mt transmitting antennas and Mr receiving antennas, and the imaging area is divided into (x) receiving antennasmin,ymin)、(xmin,ymax)、(xmax,ymin)、(xmax,ymax) Four points are enclosed.
The further technical scheme is that the CBP algorithm is used for imaging the sampling set pixel units in sequence, and the method for reconstructing the rest pixel units according to the non-precise augmented Lagrangian IALM algorithm is as follows:
the distance up integration is done for equation (8) and summed along the azimuth direction:
if the designed MIMO array has Q sets of transceiving combinations in the azimuth direction, the imaging result of the nth pixel point is:
the matrix filling non-precise augmentation Lagrangian IALM algorithm solves the following problems:
the element in constraint E is 0, and the augmented lagrange function can be expressed as:
the iterative formula is as follows:
when the X is updated,
when the E is updated, the data of the E is updated,
the Y is updated, and the Y is updated,
Yk+1=Yk+μk(M-Xk+1-Ek+1) (14)
if the elements of the whole area are imaged, the operation time is too long, pixel points in a certain proportion are randomly extracted from all pixel points corresponding to the imaging area, the extracted pixel points are imaged by using a convolution inverse projection algorithm to obtain imaging results of the pixel points, and then the values of all the pixel points are restored through a matrix filling algorithm, namely the imaging results of all the pixel points are finally obtained.
Adopt the produced beneficial effect of above-mentioned technical scheme to lie in: according to the method, the MIMO radar two-dimensional imaging CBP algorithm and the IALM algorithm are combined, two-dimensional imaging is carried out on the randomly extracted pixel points, and then the IALM algorithm is used for restoring the whole imaging area, so that the imaging time of the CBP algorithm is greatly shortened, and the method can be widely applied to near-field imaging of the MIMO radar.
Drawings
The present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
FIG. 1 is a flow chart of a method according to an embodiment of the invention;
FIG. 2 is a schematic diagram of the method of the present invention for introducing phase error using the principle of equivalent phase center;
FIG. 3 is a flow chart of the combination of CBP algorithm and IALM algorithm in the method according to the embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention are clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, but the present invention may be practiced in other ways than those specifically described and will be readily apparent to those of ordinary skill in the art without departing from the spirit of the present invention, and therefore the present invention is not limited to the specific embodiments disclosed below.
As shown in fig. 1, the embodiment of the invention discloses a MIMO radar two-dimensional imaging method based on improved CBP, and the algorithm includes the following steps:
step 1): establishing a near-field millimeter wave MIMO radar echo model;
step 2): designing a group of matched filter functions for counteracting errors introduced by an equivalent phase center principle, and compensating the phase errors;
step 3): dividing an imaging area, detecting a target, and mapping the target to an imaging unit;
step 4): extracting pixel units according to the probability to form an imaging pixel unit set;
step 5): and sequentially imaging the sampling set pixel units by using a CBP algorithm, and reconstructing the rest pixel units according to an Inaccurate Augmented Lagrange (IALM) algorithm to finish the MIMO radar two-dimensional imaging, as shown in figure 3.
The above steps are described below with reference to specific methods:
the method for establishing the near-field millimeter wave MIMO radar echo model in the step 1) comprises the following steps:
let At(xt,yt) And Ar(xr,yr) Respectively, coordinates of a pair of transmitting and receiving antennas, for an object M (x) in the object arean,yn) At respective distances R from the transmitting antenna and the receiving antennat,n、Rr,nWherein:
the emission signals are:
s(t)=exp[-j2πft] (3)
wherein f ═ f0+(n-1)Δf,f0Is the initial frequency.
The FMCW radar mixing signal is:
Sn(xt,yt;xr,yr;Kω)=σ(xn,yn)exp[-jKω(Rt,n+Rr,n)] (4)
wherein Kw2 pi f/c is the spatial frequency. Therefore, the target area signal received by the transceiving antenna is:
as shown in fig. 2, a schematic diagram of introducing a phase error by using an equivalent phase center principle, and the method for introducing a matched filter function by using the step 2) of eliminating the equivalent phase center principle is as follows:
is equivalent to a self-transmitting and self-receiving antenna, and the position of the antenna is Ae(xe,ye) And then:
the compensation function can be expressed as:
the compensated signals are:
the step 3) of dividing the target area, detecting and mapping the target area to the imaging unit comprises the following steps:
and establishing a space coordinate of the MIMO radar, taking the MIMO radar array as an x axis, and taking the array direction as a y axis vertically. Dividing an imaging area, wherein the MIMO radar has Mt transmitting antennas and Mr receiving antennas, and the imaging area is formed by (x)min,ymin)、(xmin,ymax)、(xmax,ymin)、(xmax,ymax) Four points are enclosed.
And 4, extracting certain point target pixels according to the probability, forming a pixel set so as to improve the speed of an imaging algorithm, and selecting a specific point target in a target area for imaging.
As shown in fig. 3, in step 5, imaging is performed by using a CBP algorithm, and an imaging matrix is reconstructed by using an Imprecise Augmented Lagrangian (IALM) algorithm, and a method for obtaining a target image finally is as follows:
and 5.1) integrating the distance upwards on the formula (8) and summing along the azimuth direction. Assuming that the designed MIMO array has Q sets of transceiving combinations in the azimuth direction, the imaging result of the nth pixel point is:
step 5.2) matrix filling non-precise augmentation Lagrange (IALM) algorithm, which solves the following problems:
the element in constraint E is 0, and the augmented lagrange function can be expressed as:
the iterative formula is as follows:
when the X is updated,
when the E is updated, the data of the E is updated,
the Y is updated, and the Y is updated,
Yk+1=Yk+μk(M-Xk+1-Ek+1) (14)
if the elements of the whole area are imaged, the operation time is too long, pixel points in a certain proportion are randomly extracted from all pixel points corresponding to the imaging area, the extracted pixel points are imaged by using a convolution inverse projection algorithm to obtain imaging results of the pixel points, and then the values of all the pixel points are restored through a matrix filling algorithm, namely the imaging results of all the pixel points are finally obtained.
Claims (1)
1. A MIMO radar two-dimensional imaging method based on improved CBP is characterized by comprising the following steps:
establishing a near-field millimeter wave MIMO radar echo model;
designing a group of matched filter functions to compensate phase errors introduced by an equivalent phase center principle;
dividing an imaging area, detecting a target, and mapping the target to an imaging unit;
extracting pixel units according to the probability to form an imaging pixel unit set;
sequentially imaging the sampling set pixel units by using a CBP algorithm, and reconstructing the rest pixel units according to a non-precise augmented Lagrange IALM algorithm to complete the two-dimensional imaging of the MIMO radar;
the method for establishing the near-field millimeter wave MIMO radar echo model comprises the following steps:
let At(xt,yt) And Ar(xr,yr) Respectively, coordinates of a pair of transmitting and receiving antennas, for an object M (x) in the object arean,yn) At respective distances R from the transmitting antenna and the receiving antennat,n、Rr,nWherein:
the emission signals are:
s(t)=exp[-j2πft] (3)
wherein the transmission frequency f ═ f0+(n-1)Δf,f0Is the initial frequency;
the FMCW radar mixing signal is:
Sn(xt,yt;xr,yr;Kω)=σ(xn,yn)exp[-jKω(Rt,n+Rr,n)] (4)
wherein KwSince 2 pi f/c is the spatial frequency, the signals received by the transmitting and receiving antenna in the target area are:
the method for designing a group of matched filter functions and compensating the phase error introduced by the equivalent phase center principle comprises the following steps:
is equivalent to a self-transmitting and self-receiving antenna, and the position of the antenna is Ae(xe,ye) And then:
the compensation function is expressed as:
the compensated signals are:
the method for dividing the imaging area and detecting the target and mapping to the imaging unit is as follows:
establishing MIMO radar space coordinates, taking an MIMO radar array as an x axis, taking the array direction as a y axis vertically, and dividing an imaging area, wherein the MIMO radar is provided with Mt transmitting antennas and Mr receiving antennas, and the imaging area is divided into (x) receiving antennasmin,ymin)、(xmin,ymax)、(xmax,ymin)、(xmax,ymax) Area surrounded by four points, xminIs the minimum value of the abscissa, xmaxIs the maximum value of the abscissa, yminIs the ordinate minimum, ymaxIs the maximum value of the ordinate;
the method for imaging the sampling set pixel units in sequence by using the CBP algorithm and reconstructing the rest pixel units according to the non-precise augmented Lagrange IALM algorithm comprises the following steps:
the distance up integration is done for equation (8) and summed along the azimuth direction:
if the designed MIMO array has Q sets of transceiving combinations in the azimuth direction, the imaging result of the nth pixel point is:
the matrix filling non-precise augmentation Lagrangian IALM algorithm solves the following problems:
wherein X is a band filling matrix, and E is a constraint matrix;
the element in constraint E is 0, and the augmented lagrange function can be expressed as:
wherein Y is a Lagrange multiplier, M is an original data matrix in coordinates, and mu is a penalty factor more than 0;
the iterative formula is as follows:
when the X is updated,
when the E is updated, the data of the E is updated,
the Y is updated, and the Y is updated,
Yk+1=Yk+μk(M-Xk+1-Ek+1) (14)
if the elements of the whole area are imaged, the operation time is too long, pixel points in a certain proportion are randomly extracted from all pixel points corresponding to the imaging area, the extracted pixel points are imaged by using a convolution inverse projection algorithm to obtain imaging results of the pixel points, and then the values of all the pixel points are restored through a matrix filling algorithm, namely the imaging results of all the pixel points are finally obtained.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010603678.4A CN111665501B (en) | 2020-06-29 | 2020-06-29 | MIMO radar two-dimensional imaging method based on improved CBP |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010603678.4A CN111665501B (en) | 2020-06-29 | 2020-06-29 | MIMO radar two-dimensional imaging method based on improved CBP |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111665501A CN111665501A (en) | 2020-09-15 |
CN111665501B true CN111665501B (en) | 2021-08-24 |
Family
ID=72390242
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010603678.4A Active CN111665501B (en) | 2020-06-29 | 2020-06-29 | MIMO radar two-dimensional imaging method based on improved CBP |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111665501B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113917461B (en) * | 2021-10-21 | 2022-10-28 | 苏州威陌电子信息科技有限公司 | MIMO radar imaging method and system |
CN115574711B (en) * | 2022-11-24 | 2023-02-28 | 天津理工大学 | Measuring method and system for environmental engineering |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103969642A (en) * | 2014-05-04 | 2014-08-06 | 中国电子科技集团公司第四十一研究所 | Phase compensating method used for multi-probe array imaging |
CN105652271A (en) * | 2015-12-29 | 2016-06-08 | 电子科技大学 | Super-resolution processing method for augmented Lagrangian real-beam radar angle |
CN105676219A (en) * | 2016-01-11 | 2016-06-15 | 桂林电子科技大学 | Quadrature-phase modulation-based MIMO radar three-dimensional imaging method |
CN108919229A (en) * | 2018-10-09 | 2018-11-30 | 中国科学院沈阳自动化研究所 | A kind of matrix reconstruction imaging method based on convolution inverse projection |
CN209821372U (en) * | 2018-12-29 | 2019-12-20 | 清华大学 | Electromagnetic imaging device for active microwave millimeter wave security inspection equipment and security inspection equipment |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11002828B2 (en) * | 2018-01-12 | 2021-05-11 | Tiejun Shan | Method of using a multi-input and multi-output antenna (MIMO) array for high-resolution radar imaging and wireless communication for advanced driver assistance systems (ADAS) and autonomous driving |
-
2020
- 2020-06-29 CN CN202010603678.4A patent/CN111665501B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103969642A (en) * | 2014-05-04 | 2014-08-06 | 中国电子科技集团公司第四十一研究所 | Phase compensating method used for multi-probe array imaging |
CN105652271A (en) * | 2015-12-29 | 2016-06-08 | 电子科技大学 | Super-resolution processing method for augmented Lagrangian real-beam radar angle |
CN105676219A (en) * | 2016-01-11 | 2016-06-15 | 桂林电子科技大学 | Quadrature-phase modulation-based MIMO radar three-dimensional imaging method |
CN108919229A (en) * | 2018-10-09 | 2018-11-30 | 中国科学院沈阳自动化研究所 | A kind of matrix reconstruction imaging method based on convolution inverse projection |
CN209821372U (en) * | 2018-12-29 | 2019-12-20 | 清华大学 | Electromagnetic imaging device for active microwave millimeter wave security inspection equipment and security inspection equipment |
Non-Patent Citations (2)
Title |
---|
ISAR快速CBP成像算法研究;舒明敏等;《微波学报》;20141231;第30卷(第6期);第31-35页 * |
块关联匹配与低秩矩阵超分辨融合的图像修复;马爽等;《计算机辅助设计与图形学学报》;20150228;第27卷(第2期);第271-278页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111665501A (en) | 2020-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101369018B (en) | Satellite machine combined double-base synthetic aperture radar frequency domain imaging method | |
CN111443339A (en) | Bistatic SAR space-variant correction imaging method, device, equipment and storage medium | |
CN108427115B (en) | Method for quickly estimating moving target parameters by synthetic aperture radar | |
CN104833972B (en) | A kind of bistatic CW with frequency modulation synthetic aperture radar frequency becomes mark imaging method | |
CN104749570B (en) | It is a kind of to move constant airborne biradical synthetic aperture radar target localization method | |
CN110632594B (en) | Long-wavelength spaceborne SAR imaging method | |
CN111665501B (en) | MIMO radar two-dimensional imaging method based on improved CBP | |
CN105223572B (en) | A kind of positive forward sight Bistatic SAR image processing method based on PFA algorithms | |
KR102151362B1 (en) | Image decoding apparatus based on airborn using polar coordinates transformation and method of decoding image using the same | |
CN109597076B (en) | Data processing method and device for ground-based synthetic aperture radar | |
CN111190151A (en) | System parameter design and emission power optimization method of multi-mode small satellite SAR in scanning mode | |
CN108132466B (en) | Airborne array antenna downward-looking three-dimensional imaging method and system | |
An et al. | Geosynchronous spaceborne–airborne bistatic SAR imaging based on fast low-rank and sparse matrices recovery | |
CN113702974A (en) | Method for quickly optimizing airborne/missile-borne synthetic aperture radar image | |
CN104820221B (en) | The target 3-D positioning method of many base synthetic aperture radar | |
CN109143236B (en) | Bistatic bunching SAR large-scene imaging method suitable for complex flight trajectory | |
An et al. | Simultaneous moving and stationary target imaging for geosynchronous spaceborne-airborne bistatic SAR based on sparse separation | |
Liu et al. | Ambiguities Suppression for Azimuth Multichannel SAR Based on ${L_ {2, q}} $ Regularization With Application to Gaofen-3 Ultra-Fine Stripmap Mode | |
CN108732555B (en) | Automatic driving array microwave imaging motion compensation method | |
CN112255596B (en) | Space-borne SAR deception jamming signal generation method based on spatial frequency interpolation | |
CN113670301A (en) | Airborne SAR motion compensation method based on inertial navigation system parameters | |
Pingping et al. | MIMO-SAR imaging technology for helicopter-borne based on ARC antenna array | |
CN116136595A (en) | Collaborative detection double-base forward-looking SAR imaging processing method based on two-stage scale fine adjustment | |
CN116559905A (en) | Undistorted three-dimensional image reconstruction method for moving target of bistatic SAR sea surface ship | |
CN116148851A (en) | Multi-mode-based multi-transmission multi-reception synthetic aperture radar wide swath imaging method |
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 |