阅读说明：本技术 一种mimo雷达近场三维成像方法 (MIMO radar near-field three-dimensional imaging method ) 是由 刘楠 杨明磊 王晓冉 陈晓玲 黎鹏斌 于 2020-11-02 设计创作，主要内容包括：本发明公开了一种MIMO雷达近场三维成像方法,包括：S1：获取目标回波信号；S2：对目标回波信号进行去斜、重构以及加窗处理,得到虚拟阵列加窗后的差拍信号；S3：按照三维形式将成像区域划分为若干网格,并计算雷达的不同阵元到不同方位向及距离向位置上的网格点的时延；S4：根据网格点的时延,得到虚拟阵列的相位补偿向量；S5：根据虚拟阵列的相位补偿向量对虚拟阵列加窗后的差拍信号进行补偿,并对网格点进行相干叠加和加窗叠加,得到该网格点处的回波强度值；S6：重复步骤S4-S5,遍历整个成像区域,以完成图像的重建。本发明提供的MIMO雷达近场三维成像方法避免了孔径渡越效应的影响,保证了近场三维成像结果的正确性。(The invention discloses a near-field three-dimensional imaging method for an MIMO radar, which comprises the following steps: s1: acquiring a target echo signal; s2: performing deskewing, reconstruction and windowing on the target echo signal to obtain a beat signal after the virtual array is windowed; s3: dividing an imaging area into a plurality of grids according to a three-dimensional form, and calculating time delay from different array elements of the radar to grid points in different azimuth directions and distance direction positions; s4: obtaining a phase compensation vector of the virtual array according to the time delay of the grid point; s5: compensating the windowed beat signal of the virtual array according to the phase compensation vector of the virtual array, and performing coherent superposition and windowed superposition on the grid points to obtain the echo intensity value of the grid points; s6: repeating steps S4-S5, traversing the entire imaging region to complete the reconstruction of the image. The MIMO radar near-field three-dimensional imaging method provided by the invention avoids the influence of the aperture transit effect and ensures the correctness of the near-field three-dimensional imaging result.)

1. A MIMO radar near-field three-dimensional imaging method is characterized by comprising the following steps:

s1: acquiring a target echo signal; the target echo signals are frequency-modulated continuous wave target echo signals respectively sent by a plurality of transmitting array elements according to a time-sharing mode;

s2: performing deskewing, reconstruction and windowing on the target echo signal to obtain a beat signal after virtual array windowing;

s3: dividing an imaging area into a plurality of grids according to a three-dimensional form, and calculating time delay from different array elements of the radar to grid points in different azimuth directions and distance direction positions;

s4: obtaining a phase compensation vector of the virtual array according to the time delay of the grid point;

s5: compensating the windowed beat signal of the virtual array according to the phase compensation vector of the virtual array, and performing coherent superposition and windowed superposition on the grid point to obtain an echo intensity value at the grid point;

s6: repeating steps S4-S5, traversing the entire imaging region to complete the reconstruction of the image.

2. The MIMO radar near-field three-dimensional imaging method of claim 1, wherein the expression of the target echo signal is as follows:

wherein s (-) represents the target echo signal, T_pP represents the P-th transmitting array element, P is more than or equal to 1 and less than or equal to P, P represents the number of the transmitting array elements, n_iI represents the ith scattering point, I is more than or equal to 1 and less than or equal to I, I represents the number of the scattering points, R_qRepresents the Q-th receiving array element, Q is more than or equal to 1 and less than or equal to Q, Q represents the number of the receiving array elements,representing fast time, A representing echo intensity, T representing echo time delay, f_cRepresenting the operating frequency of the radar, gamma representing the chirp rate, T_FMIndicating the frequency modulation period.

3. The MIMO radar near-field three-dimensional imaging method of claim 2, wherein the step S2 comprises:

s21: performing deskewing processing on the target echo signal and the reference signal to obtain a deskewed beat signal;

s22: reconstructing the deskewed beat signals according to the relation between the transmitting array elements to obtain beat signals corresponding to the virtual array;

s23: and windowing the beat signal corresponding to the virtual array to obtain the windowed beat signal of the virtual array.

4. The MIMO radar near-field three-dimensional imaging method according to claim 3, wherein the beat signals corresponding to the virtual array are represented as:

wherein the content of the first and second substances,representing a transmitting array element T_pReceive beat signals received by the array during transmission, an Representing the beat signal of the target echo signal after the target echo signal is subjected to the deskew.

5. The MIMO radar near-field three-dimensional imaging method according to claim 3, wherein the virtual array windowed beat signal is represented as:

wherein the content of the first and second substances,representing the beat signal to which the virtual array corresponds,representing the corresponding windowed weight vector for the virtual array.

6. The MIMO radar near-field three-dimensional imaging method of claim 1, wherein the time delay of the grid points is expressed as:

wherein the content of the first and second substances,representing grid points b_wC represents the speed of light,andrespectively expressed as:

representing grid points b_wIs determined by the three-dimensional coordinates of (a),three-dimensional coordinates representing the q-th receiving array element,representing the three-dimensional coordinates of the p-th transmit array element.

7. The MIMO radar near-field three-dimensional imaging method of claim 1, wherein the step S4 comprises:

s41: obtaining a phase compensation vector corresponding to each grid point according to the time delay of the grid point;

s42: and recombining the phase compensation vectors corresponding to each grid point according to the relation among the transmitting array elements to obtain the phase compensation vectors corresponding to the virtual array.

8. The MIMO radar near-field three-dimensional imaging method of claim 7, wherein the phase compensation vectors corresponding to the grid points are represented as:

wherein the content of the first and second substances,to representThe time transmitting array element is T_pPoly (A), poly (B)Focus is grid point b_wThe corresponding phase compensation vector of the array is received,indicating the transmitting array element as T_pThe focus point is a grid point b_wThe receiving array element is R_qA corresponding phase compensation vector, and:

wherein, P is more than or equal to 1 and less than or equal to P, P represents the number of transmitting array elements, Q is more than or equal to 1 and less than or equal to Q, Q represents the number of receiving array elements,representing fast time,. tau.representing echo time delay, f_cIndicating the operating frequency of the radar and gamma indicating the chirp rate.

9. The MIMO radar near-field three-dimensional imaging method of claim 8, wherein the phase compensation vector corresponding to the virtual array is represented as:

wherein the content of the first and second substances,to representGrid point at time b_wAnd (4) compensating the vector corresponding to the virtual array.

10. The MIMO radar near-field three-dimensional imaging method of claim 1, wherein the coherent superposition and the windowed superposition are performed on the grid points to obtain pixel values of the grid points, and the method comprises the following steps:

and carrying out coherent superposition on the grid points to obtain a coherent superposition result, wherein the coherent superposition result is as follows:

wherein the content of the first and second substances,to representTime grid point b_wAs a result of the coherent addition of (a),to representGrid point at time b_wThe phase compensation vector corresponding to the virtual array,to representA beat signal after windowing of the time virtual array;

and performing windowing superposition on the coherent superposition result in the distance direction to obtain the pixel value of the grid point, wherein the pixel value is as follows:

wherein, K represents the number of sampling times,indicating the windowing weight corresponding to the kth sampling instant,and representing the coherent superposition result corresponding to the kth sampling moment.

Technical Field

The invention belongs to the technical field of radars, and particularly relates to a near-field three-dimensional imaging method for an MIMO radar.

Background

In recent years, radar near-field detection technology is gradually becoming one of the hot problems in radar application research, and near-field radar can be applied to different occasions, for example; wall-through detection, gesture recognition, near-distance nondestructive detection and the like. At present, when the radar works under the near-field condition, detection is mainly used basically, and imaging application is less.

In the field of array signal processing, the computation complexity is often reduced by far-field assumption due to the computation complexity, but with the diversification of imaging scenes, the requirement on imaging accuracy under a near-field condition is continuously increased, and the processing according to a far-field signal model is not suitable. The existing near-field imaging algorithm based on beam forming mainly adopts a method of compressing echo signals in a distance direction first and then forming near-field beams to realize near-field imaging.

However, large aperture arrays are required to achieve high resolution imaging results. When the aperture of the array is large, the aperture transit effect can occur during imaging by adopting a method of compressing the array in the distance direction first and then forming near-field beams, so that the imaging result is staggered with the real position of a target, and the imaging result is distorted. In addition, the existing near-field three-dimensional imaging algorithm needs to respectively complete two steps of distance direction compression and near-field beam forming, and the algorithm has more steps and more complex process.

Disclosure of Invention

In order to solve the problems in the prior art, the invention provides a near-field three-dimensional imaging method for an MIMO radar. The technical problem to be solved by the invention is realized by the following technical scheme:

a MIMO radar near-field three-dimensional imaging method comprises the following steps:

s1: acquiring a target echo signal; the target echo signals are frequency-modulated continuous wave target echo signals respectively sent by a plurality of transmitting array elements according to a time-sharing mode;

s2: performing deskewing, reconstruction and windowing on the target echo signal to obtain a beat signal after virtual array windowing;

s3: dividing an imaging area into a plurality of grids according to a three-dimensional form, and calculating time delay from different array elements of the radar to grid points in different azimuth directions and distance direction positions;

s4: obtaining a phase compensation vector of the virtual array according to the time delay of the grid point;

s5: compensating the windowed beat signal of the virtual array according to the phase compensation vector of the virtual array, and performing coherent superposition and windowed superposition on the grid point to obtain an echo intensity value at the grid point;

s6: repeating steps S4-S5, traversing the entire imaging region to complete the reconstruction of the image.

In an embodiment of the present invention, the expression of the target echo signal is:

wherein s (-) represents the target echo signal, T_pP represents the P-th transmitting array element, P is more than or equal to 1 and less than or equal to P, P represents the number of the transmitting array elements, n_iI represents the ith scattering point, I is more than or equal to 1 and less than or equal to I, I represents the number of the scattering points, R_qRepresents the Q-th receiving array element, Q is more than or equal to 1 and less than or equal to Q, Q represents the number of the receiving array elements,representing fast time, A representing echo intensity, T representing echo time delay, f_cRepresenting the operating frequency of the radar, gamma representing the chirp rate, T_FMIndicating the frequency modulation period.

In one embodiment of the present invention, step S2 includes:

s21: performing deskewing processing on the target echo signal and the reference signal to obtain a deskewed beat signal;

s22: reconstructing the deskewed beat signals according to the relation between the transmitting array elements to obtain beat signals corresponding to the virtual array;

s23: and windowing the beat signal corresponding to the virtual array to obtain the windowed beat signal of the virtual array.

In one embodiment of the present invention, the beat signal corresponding to the virtual array is represented as:

wherein the content of the first and second substances,representing a transmitting array element T_pDuring transmission, the beat signals received by the array are received, andrepresenting the beat signal of the target echo signal after the target echo signal is subjected to the deskew.

In one embodiment of the present invention, the virtual array windowed beat signal is represented as:

wherein the content of the first and second substances,representing the beat signal to which the virtual array corresponds,representing the corresponding windowed weight vector for the virtual array.

In one embodiment of the present invention, the time delay of the mesh point is expressed as:

wherein the content of the first and second substances,representing grid points b_wC represents the speed of light,andrespectively expressed as:

representing grid points b_wIs determined by the three-dimensional coordinates of (a),three-dimensional coordinates representing the q-th receiving array element,representing the three-dimensional coordinates of the p-th transmit array element.

In one embodiment of the present invention, step S4 includes:

s41: obtaining a phase compensation vector corresponding to each grid point according to the time delay of the grid point;

s42: and recombining the phase compensation vectors corresponding to each grid point according to the relation among the transmitting array elements to obtain the phase compensation vectors corresponding to the virtual array.

In one embodiment of the present invention, the phase compensation vector corresponding to the grid point is represented as:

wherein the content of the first and second substances,to representThe time transmitting array element is T_pThe focus point is a grid point b_wThe corresponding phase compensation vector of the array is received,indicating the transmitting array element as T_pThe focus point is a grid point b_wThe receiving array element is R_qA corresponding phase compensation vector, and:

wherein, P is more than or equal to 1 and less than or equal to P, P represents the number of transmitting array elements, Q is more than or equal to 1 and less than or equal to Q, Q represents the number of receiving array elements,representing fast time,. tau.representing echo time delay, f_cIndicating the operating frequency of the radar and gamma indicating the chirp rate.

In one embodiment of the present invention, the phase compensation vector corresponding to the virtual array is represented as:

wherein the content of the first and second substances,to representGrid point at time b_wAnd (4) compensating the vector corresponding to the virtual array.

In an embodiment of the present invention, performing coherent superposition and windowed superposition on the grid points to obtain pixel values of the grid points, includes:

and carrying out coherent superposition on the grid points to obtain a coherent superposition result, wherein the coherent superposition result is as follows:

wherein the content of the first and second substances,to representTime grid point b_wAs a result of the coherent addition of (a),to representGrid point at time b_wThe phase compensation vector corresponding to the virtual array,to representA beat signal after windowing of the time virtual array;

and performing windowing superposition on the coherent superposition result in the distance direction to obtain the pixel value of the grid point, wherein the pixel value is as follows:

wherein, K represents the number of sampling times,indicating the windowing weight corresponding to the kth sampling instant,and representing the coherent superposition result corresponding to the kth sampling moment.

The invention has the beneficial effects that:

1. according to the MIMO radar near-field three-dimensional imaging method, the azimuth direction is processed firstly, and then the distance direction is processed, so that the influence of the aperture transit effect is avoided, and the correctness of a near-field three-dimensional imaging result is ensured;

2. the MIMO radar near-field three-dimensional imaging method provided by the invention combines the processing of the azimuth direction and the processing of the distance direction into the phase compensation, simplifies the algorithm steps and improves the operation efficiency of the algorithm.

The present invention will be described in further detail with reference to the accompanying drawings and examples.

Drawings

Fig. 1 is a schematic flow chart of a MIMO radar near-field three-dimensional imaging method according to an embodiment of the present invention;

fig. 2 is a diagram of a MIMO radar front architecture used in simulations 1 and 2;

FIG. 3 is an enlarged view of a portion of FIG. 2;

FIG. 4 is a spatial position distribution diagram of 5 scattering points set in simulation 1;

FIG. 5 is a slice of the imaging results obtained in simulation 1 using the method of the present invention at a fixed distance of 5m from the wavefront reference array element;

FIG. 6 is a slice of the result of imaging using the method of the present invention in simulation 1 with a fixed X coordinate of 0 m;

FIG. 7 is a slice of the imaging results obtained in simulation 1 using the method of the present invention at a fixed Y coordinate of 0 m;

fig. 8 is a diagram of the positional form of the MIMO radar front in space used in simulation 2;

FIG. 9 is a schematic diagram of the geometrical relationship between the array elements and the targets set in simulation 2;

FIG. 10 is a schematic diagram of the geometrical relationship between the array elements and the target in simulation 2, assuming that the shielded part has no scattering points;

FIG. 11 is a cloud of imaging points from simulation 2 using the method of the present invention in combination with a center of gravity ranging method;

FIG. 12 is a diagram of the array format used to verify the effect of aperture transit on different imaging algorithms;

fig. 13 is a slice diagram of the imaging result obtained when the array shown in fig. 12 is adopted in simulation 3, and the array element is referenced to 5m according to the fixed distance front plane obtained by the existing near-field imaging algorithm based on beam forming;

FIG. 14 is a slice of the imaging results obtained by the method of the present invention at a fixed wavefront reference array element of 5m after the array of FIG. 12 was used in simulation 3;

FIG. 15 is a diagram of a larger array format of apertures used to further verify the effect of aperture transit on different imaging algorithms;

fig. 16 is a slice diagram of the imaging result obtained when the array shown in fig. 15 is adopted in simulation 3, and the array element is referenced to 5m according to the fixed distance front plane obtained by the existing near-field imaging algorithm based on beam forming;

fig. 17 is a slice of the imaging results obtained by the method according to the invention at a fixed distance front with reference to the array element of 5m after the array of fig. 15 has been used in simulation 3.

Detailed Description

The present invention will be described in further detail with reference to specific examples, but the embodiments of the present invention are not limited thereto.

Example one

Referring to fig. 1, fig. 1 is a schematic flow chart of a MIMO radar near-field three-dimensional imaging method according to an embodiment of the present invention, including:

s1: acquiring a target echo signal; the target echo signals are frequency modulated continuous waves which are respectively sent by a plurality of transmitting array elements according to a time-sharing mode.

The near-field imaging method provided by the embodiment is applied to the MIMO radar, the signal transmission modes of a plurality of transmission array elements are firstly adjusted so as to facilitate the subsequent near-field imaging,

specifically, assuming that the radar has P transmitting array elements and Q receiving array elements, the P transmitting array elements are set to respectively send out frequency modulated continuous waves in a time-sharing manner, and signals transmitted by the transmitting array elements can be represented as:

wherein, T_FMWhich represents the period of the frequency modulation,representing fast time, f_cWhich represents the operating frequency of the radar, i.e. the carrier frequency, gamma represents the chirp rate.

The received target echo signal may be expressed as:

wherein s (-) represents the target echo signal, T_pP represents the P-th transmitting array element, P is more than or equal to 1 and less than or equal to P, P represents the number of the transmitting array elements, n_iI represents the ith scattering point, I is more than or equal to 1 and less than or equal to I, I represents the number of the scattering points, R_qRepresents the Q-th receiving array element, Q is more than or equal to 1 and less than or equal to Q, Q represents the number of the receiving array elements,representing fast time, A representing echo intensity, τ representing echo time delay, and

wherein the content of the first and second substances,representing the p-th transmitting element to the p-thThe distance between the i scattering points is,represents the distance from the q-th receiving array element to the i-th scattering point, and c represents the speed of light.

S2: the target echo signal is subjected to deskew, reconstruction and windowing processing to obtain a beat signal after virtual array windowing, and the method comprises the following steps:

s21: the target echo signal and the reference signal are subjected to deskew processing to obtain a deskewed beat signal, and the expression is as follows:

wherein the content of the first and second substances,representing a transmitting array element T_pThe emitted frequency-modulated signal passes through a scattering point n_iBack to the receiving array element R_qSubject the target echo signal to a deskewed beat signal,representing the conjugate of the transmitted signal.

S22: and reconstructing the deskewed beat signals according to the relation between the transmitting array elements to obtain the beat signals corresponding to the virtual array.

Specifically, assuming that the array has P transmit elements and Q receive elements, for a certain time instantArray element T_pThe beat signal received by the receiving array at the time of transmission can be expressed as:

the beat signal corresponding to the virtual array can be expressed as:

s23: and windowing the beat signal corresponding to the virtual array to obtain the beat signal after the virtual array is windowed.

Specifically, assume time of dayThe windowing weight vector corresponding to the virtual array is:

the virtual array windowed beat signal can be represented as:

s3: and dividing the imaging area into a plurality of grids according to a three-dimensional form, and calculating time delay from different array elements of the radar to grid points in different azimuth directions and distance directions.

Specifically, the imaging region is divided into a plurality of grids according to a three-dimensional form, and a plurality of grid points, that is, pixel points, are obtained.

Suppose a certain grid point b_wHas three-dimensional coordinates ofThe p-th transmitting array element transmits signals through b_wAnd reflecting, when the q-th receiving array element receives the signal, the time delay is as follows:

wherein the content of the first and second substances,representing grid points b_wC represents the speed of light,andrespectively expressed as:

representing grid points b_wIs determined by the three-dimensional coordinates of (a),three-dimensional coordinates representing the q-th receiving array element,representing the three-dimensional coordinates of the p-th transmit array element.

S4: obtaining a phase compensation vector of the virtual array according to the time delay of the grid point, wherein the phase compensation vector comprises the following steps:

s41: and obtaining a phase compensation vector corresponding to each grid point according to the time delay of each grid point.

Specifically, assuming that the array has P transmit elements and Q receive elements, for a certain time instantWhen array element T₁Transmitting and focusing point as grid point b_wThen, the phase compensation vector corresponding to the receiving array is:

correspondingly, when array element T_pTransmitting and focusing point as grid point b_wWhen the temperature of the water is higher than the set temperature,

wherein the content of the first and second substances,indicating the transmitting array element as T_pThe focus point is a grid point b_wThe receiving array element is R_qA corresponding phase compensation vector, and:

s42: and recombining the phase compensation vectors corresponding to each grid point according to the relation among the transmitting array elements to obtain the phase compensation vectors corresponding to the virtual array.

Specifically, the phase compensation vector corresponding to the entire virtual array can be expressed as:

s5: and compensating the windowed beat signal of the virtual array according to the phase compensation vector of the virtual array, and performing coherent superposition and windowed superposition on the grid points to obtain the echo intensity value at the grid points.

Firstly, correspondingly compensating the beat signals of the virtual array according to the phase compensation vectors obtained in the last step.

For a certain moment of timeWhen the focus point is a netGrid point b_wThen, the phase compensation vector corresponding to the virtual array isThe beat signal after the virtual array windowing isCorrespondingly compensating the beat signals of the virtual array, and performing coherent superposition on each pixel point:

wherein the content of the first and second substances,to representTime grid point b_wAnd (4) coherent superposition results.

The coherent addition results are then added in the distance direction.

Assuming that the total sampling is performed for K times, the difference frequency signals of K time moments are obtained, and the coherent superposition result obtained at the K sampling time moment isAt this time, the corresponding windowing weight value isThe results from the K sub-samples are then superimposed as:

the superposition result is the pixel value of the grid point (pixel point).

S6: repeating steps S4-S5, traversing the entire imaging region to complete the reconstruction of the image.

Specifically, the imaging region is divided into meshes in a three-dimensional form, each mesh point b_wAll as focus points, and are repeated according to the steps S4 and S5 respectively to obtain the value out (b) of the corresponding pixel point_w) Thereby completing the reconstruction of the image.

According to the MIMO radar near-field three-dimensional imaging method provided by the embodiment, the azimuth direction processing is performed firstly, and then the distance direction processing is performed, so that the influence of the aperture transit effect is avoided, and the correctness of the near-field three-dimensional imaging result is ensured. Meanwhile, in the subsequent processing process, the processing in the azimuth direction and the processing in the distance direction are combined to phase compensation, so that the steps of the algorithm are simplified, and the operation efficiency of the algorithm is improved.

The radar near-field three-dimensional imaging method provided by the invention can obtain an accurate three-dimensional image, avoids the influence of an aperture transit effect, does not generate the offset of an imaging position, and can realize higher imaging resolution under the condition of the same array element number aiming at the MIMO array.

Example two

The beneficial effects of the present invention are further verified and explained by the simulation experiment.

Simulation 1: the method of the present invention is used to simulate the imaging of multiple point targets.

1.1 simulation conditions:

the signal adopts a sawtooth frequency modulation continuous wave form, the carrier frequency of a transmitting signal is set to be 300GHz, the signal bandwidth is set to be 10GHz, the frequency modulation period is set to be 100us, and the sampling rate is set to be 16 MHz. The distance between receiving array elements of the MIMO array is half wavelength, namely 0.5mm, the distance between transmitting array elements is 10.5mm, the virtual array formed in this way is non-porous, specific array form diagrams are shown in fig. 2 and fig. 3, fig. 2 is a structure diagram of the MIMO radar front provided by the embodiment of the invention, and fig. 3 is a partial enlarged view of fig. 2. Wherein, 5 scattering point targets are respectively located: [0m,0m,3m ], [ -0.12m,0m,3m ], [0.12m,0m,3m ], [0m, -0.12m,3m ], [0m,0.12m,3m ], wherein the 1 st number represents the X-axis coordinate of the scattering point projected on the wavefront, the 2 nd number represents the Y-axis coordinate of the scattering point projected on the wavefront, the 3 rd number represents the distance between the scattering point and the reference array element, and the spatial position distribution is shown in fig. 4.

1.2 simulation content and result analysis:

under the simulation condition of the 1.1, the method of the invention is used for carrying out three-dimensional imaging simulation on the 5 scattering point targets, and the obtained imaging result is a three-dimensional data matrix. In order to more intuitively show the imaging result, the imaging result is displayed in a slicing mode. Referring to fig. 5 to 7, fig. 5 is a slice diagram of an imaging result corresponding to a fixed distance array plane reference array element of 5m according to an embodiment of the present invention, fig. 6 is a slice diagram of an imaging result corresponding to a fixed X coordinate of 0m according to an embodiment of the present invention, and fig. 7 is a slice diagram of an imaging result corresponding to a fixed Y coordinate of 0m according to an embodiment of the present invention. As can be seen from fig. 5, 6, and 7, the imaging result can correspond to the real position of the target, and the correctness of the three-dimensional imaging algorithm provided by the present invention is verified.

Simulation 2: the method of the invention is used for imaging simulation of a cylinder.

2.1 simulation conditions:

the signal adopts a sawtooth frequency modulation continuous wave form, the carrier frequency of a transmitting signal is set to be 300GHz, the signal bandwidth is set to be 10GHz, the frequency modulation period is set to be 100us, and the sampling rate is set to be 16 MHz. The distance between receiving array elements of the MIMO array is half wavelength, namely 0.5mm, the distance between transmitting array elements is 10.5mm, and the virtual array formed by the method is non-porous. A specific array format is shown in fig. 2, and its position in space is shown in fig. 8, where the symbol "x" represents the transmitting array element and the symbol "o" represents the receiving array element.

The object imaged by the simulation experiment is a cylindrical model and a background surface, the background surface is parallel to the array surface, is a rectangular surface with the length of 4m by 4m, the height is-5 m, the value range of an x axis is (-2m, 2m), and the value range of a y axis is (-2m, 2 m). The center coordinates of the upper bottom surface of the cylinder are (0, 0-4.5 m), the radius is 1m, the height of the cylinder is 0.5m, and the geometrical relationship between the array elements and the target is shown in FIG. 9.

Assuming that the portion of the background surface shielded by the cylindrical surface has no scattering points, the side surface of the cylindrical surface has no scattering points, and the area projected from the transmitting array element to the upper bottom surface to the background surface also has no scattering points, the relationship between the array element and the scattering points is shown in fig. 10.

2.2 simulation content and result analysis:

under the simulation condition of the 2.1, the method is used for carrying out three-dimensional imaging simulation on the imaging object, then the center of gravity of the output result of each wave position is calculated, a three-dimensional measurement point is given according to the distance corresponding to the center of gravity and the angle corresponding to the wave position, all three-dimensional measurement points are drawn, and the imaging point cloud picture is obtained.

The resulting imaged dot cloud results are shown in fig. 11. As can be seen from the figure, the imaging result can correspond to the actual situation, and the correctness of the three-dimensional imaging algorithm provided by the invention is verified.

Simulation 3: the method is compared with the existing near-field imaging algorithm based on beam forming in a simulation mode.

3.1, simulation conditions:

because the aperture transit effect is mainly related to the aperture size of the array, the array is in a form of one-shot multiple-shot in order to simplify the verification step and reduce the calculation amount.

The signal adopts a sawtooth frequency modulation continuous wave form, the carrier frequency of a transmitting signal is set to be 300GHz, the signal bandwidth is set to be 10GHz, the frequency modulation period is set to be 100us, and the sampling rate is set to be 16 MHz. The point target is located at a position having an azimuth angle of 10 degrees, a pitch angle of 30 degrees, and a distance of 5 m.

3.2 simulation content and result analysis:

under the simulation condition of 3.1, the imaging effect of the method of the present invention is compared with that of the existing near-field imaging algorithm based on beam forming, the array forms are respectively shown in fig. 12 and fig. 15, when the array shown in fig. 12 is adopted, the imaging result slice corresponding to the fixed distance front reference array element 5m obtained according to the existing near-field imaging algorithm based on beam forming is shown in fig. 13, and the imaging result slice corresponding to the fixed distance front reference array element 5m obtained according to the method of the present invention is shown in fig. 14. When the array shown in fig. 15 is used, the slice image of the imaging result obtained by the conventional near-field imaging algorithm based on beam forming when the distance front is fixed to the reference array element by 5m is shown in fig. 16, and the slice image of the imaging result obtained by the method of the present invention when the distance front is fixed to the reference array element by 5m is shown in fig. 17. Wherein, the imaging results are the imaging results of the azimuth-elevation dimension.

In all the figures representing the imaging results, the symbol "+" represents the true position of the target. As can be seen from fig. 13, because of the influence of the aperture transit effect, the imaging result obtained by using the existing near-field imaging algorithm based on beam forming has a deviation from the actual position of the target, because the aperture ratio of the array shown in fig. 12 is small, and the influence of the aperture transit effect is not so large, so the deviation is not obvious. In fig. 14, it can be seen that the imaging result obtained by using the method of the present invention can be perfectly aligned with the actual position of the target without deviation and without being affected by the aperture transit effect.

To be able to show the aperture transit effect more clearly, an array with larger apertures is used, as shown in fig. 15. In fig. 16, it can be seen that, due to the effect of the aperture transit effect, the imaging result obtained by using the existing near-field imaging algorithm based on beamforming has a large deviation from the actual position of the target. In fig. 17, it can be seen that the imaging result obtained by using the method of the present invention can be perfectly aligned with the actual position of the target without being affected by the aperture transit effect. Meanwhile, by comparing fig. 16 and fig. 13, it can be demonstrated that the influence of the aperture transit effect is larger and larger as the aperture of the array is increased.

In conclusion, the method can avoid the influence of the aperture transit effect and ensure the correctness of the near-field three-dimensional imaging result. Meanwhile, the processing of the azimuth direction and the processing of the distance direction are combined on the phase compensation, so that the steps of the algorithm are simplified, and the operation efficiency of the algorithm is improved.

The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.