Disclosure of Invention
In order to solve the problems, the invention aims to provide a polarization information driven microwave correlation imaging method, which solves the technical problems that in the prior art, the research on polarization information is very blank and the research on the application of polarization information in the whole imaging process is lacking, can image a target with high resolution and high robustness by using the polarization information, and can perform instantaneous polarization measurement while imaging, has simple engineering realization and strong engineering application potential, can be popularized to actual scenes such as map reconnaissance, target identification and the like, and can rapidly and accurately acquire the imaging of a target area and the target polarization information.
In order to achieve the above object, the technical scheme of the present invention is as follows:
the invention provides a polarization information driven microwave correlation imaging method, which comprises the following steps:
s1, constructing a full polarization frequency hopping radar matrix;
s2, forming horizontal polarization and vertical polarization polarizing space-time random radiation fields;
s3, iterative imaging is performed by utilizing an orthogonal complement space of the random radiation field.
Further, the step S1 includes the following steps:
S11, arranging antennas of the frequency hopping radar at equal intervals along the wing direction of the aircraft to form a row of frequency hopping radar array in a multi-transmitting single-receiving mode;
S12, taking the middle full-polarization receiving antenna as a boundary, respectively distributing the horizontal polarization antenna and the vertical polarization antenna on two sides of the middle full-polarization receiving antenna;
s13, changing a target imaging area of the frequency hopping radar from a square area to a diamond area.
Further, the step S2 includes the following steps:
S21, setting a vector from the center of the radar array to a target scattering point position as Dividing M full polarization receiving and transmitting antennas into an H antenna and a V antenna, wherein the H antenna and the V antenna are respectively arranged at two sides of an antenna array;
S22, acquiring the incident field intensity of the whole antenna array to an imaging area;
S23, independently transmitting random frequency hopping signals through antennas with different polarization directions and randomly sampling the receiving antennas to form horizontal polarization space-time random radiation fields and vertical polarization space-time random radiation fields.
Further, the step S22 includes the following steps:
S221, the incident field intensity of the ith transmitting antenna to the target area is as follows:
wherein f Hi (t) and f Vi (t) are the random frequencies of the H antenna and the V antenna, respectively, γ i (t) is the random horizontal polarization angle, η i (t) is the random vertical polarization angle, Is an antenna position vector;
S222, the incident field intensity of the whole antenna array to an imaging area is as follows:
Further, the step S23 includes the following steps:
s231 formation of target region the random radiation field is of the formula:
s232, in the microwave correlation forward-looking imaging process, the random radiation field irradiates the target, then reflects to form an echo, and is received by a receiving antenna to obtain a target echo E Hr、EVr as follows:
S233, obtaining target scattering characteristics and target polarization scattering characteristics through inversion:
Further, the step S3 includes the following steps:
S31, if polarization information is not considered, TSVD processing is carried out on the space-time radiation field with insufficient randomness, then compression correlation operation is carried out on the echo and the space-time random radiation field, finally, a target scattering characteristic matrix is obtained, correlation imaging is achieved, and a formula of the target scattering characteristic matrix is simplified to E r=Es σ under the condition that noise is not considered;
s32, if the polarization scattering matrix of the target is considered, solving the instantaneous polarization scattering matrix of the target by using a microwave correlation imaging method.
Further, the instantaneous polarization scattering matrix echo in S32 is expressed as:
EHr=EHsSHH+EVsSHV
EVr=EHsSHV+EVsSVV
if the above formula is to be solved, the radiation field of the H antenna and the radiation field of the V antenna are orthogonal, namely:
EVsEHs=0
Further, the step S32 includes the following steps:
S321, obtaining a horizontal polarization radiation field orthogonal complementary matrix and a vertical polarization radiation field orthogonal complementary matrix through QR decomposition, and carrying out cross multiplication on the orthogonal complementary matrix and echo information;
s322, forming a new association type and iterative imaging by the radiation field after the matrix obtained in the step S321 is processed in the same way.
Further, the step S321 includes the following steps:
S3211, performing QR decomposition on a radiation field matrix of a target, and calculating to obtain orthogonal complement spaces E Hsb and E Vsb of the radiation field matrices E Hs and E Vs;
S3212, multiplying the target HV echo by the HV radiation field orthogonal complement space respectively to obtain:
EHsbEHr=EHsbEHsSHH+EHsbEVsSHV
EVsbEHr=EVsbEHsSHH+EVsbEVsSHV
EHsbEVr=EHsbEHsSVH+EHsbEVsSVV
EVsbEVr=EVsbEHsSVH+EVsbEVsSVV
the characteristics of the orthogonal complement matrix are as follows:
EHsbEHs=0
EVsbEVs=0
Namely:
By adopting the technical scheme, the invention has the following advantages:
The invention provides a polarization information driven microwave associated imaging method, which divides an antenna of a frequency hopping radar array into an H antenna and a V antenna, utilizes random frequency modulation and random sampling of each antenna to respectively generate a horizontal polarization radiation field and a vertical polarization radiation field, utilizes orthogonal complement space of the HV antenna radiation field to carry out cross multiplication processing on echo information, then forms new associated and iterative imaging with the radiation field after the same processing, can carry out high-resolution and high-robustness imaging on a target by utilizing the polarization information, carries out instantaneous polarization measurement while imaging, has better effect than the traditional instantaneous polarization measurement method under the condition of low time-frequency product, improves the anti-interference capability engineering implementation of the microwave associated imaging, has strong engineering application potential, and can be popularized to actual scenes such as map reconnaissance, target identification and the like to rapidly and accurately acquire the imaging of a target area and the target polarization information.
Detailed Description
In the following detailed description of the embodiments of the present invention, reference is made to the accompanying drawings, in which it is to be noted that, in this document, relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.
Fig. 1 shows a flowchart of a method for polarization information driven microwave correlation imaging according to the present invention, and the method according to the present invention is further explained below with reference to fig. 1. The specific steps and effects of the invention are as follows:
FIG. 2 shows a geometrical imaging schematic of a polarization information driven microwave-dependent imaging application scenario of the present invention;
s1, constructing a full polarization frequency hopping radar matrix;
s1 comprises the following specific steps:
S11. as shown in figure 2, a row of antennas (1 2....n.) to form a frequency hopping radar array, a multiple transmit-receive form forming a column;
s12, the middle full-polarization receiving antenna A is used as a boundary, and the horizontal polarization antenna H and the vertical polarization antenna V are respectively distributed on two sides, so that left and right blurring can be effectively restrained, and the device is simple in structure and high in practicability;
s13, changing the target imaging area B from a general square area to a diamond area, so that the space correlation of random radiation fields generated by radar arrays arranged facing the aircraft wing in each scattering unit can be effectively reduced, and the imaging effect is improved.
S2, forming horizontal polarization and vertical polarization polarizing space-time random radiation fields;
S2 comprises the following specific steps:
S21, setting a vector from the center of the radar array to a target scattering point position as Dividing M full-polarization receiving and transmitting antennas into an H antenna and a V antenna, wherein the H antenna and the V antenna are respectively arranged at two sides of an antenna array, and the random frequencies of the two antennas are different;
S22, acquiring the incident field intensity of the whole antenna array to an imaging area;
The step S22 includes the steps of:
S221. Random frequencies of the H antenna and the V antenna are f Hi (t) and f Vi (t) respectively, a random horizontal polarization angle is gamma i (t), a random vertical polarization angle is eta i (t), and an antenna position vector is The incident field strength of the ith transmit antenna to the target area is as follows:
S222, the incident field intensity of the whole antenna array to an imaging area is as follows:
S23, independently transmitting random frequency hopping signals through antennas with different polarization directions and randomly sampling the receiving antennas to form horizontal polarization space-time random radiation fields and vertical polarization space-time random radiation fields.
S23 comprises the following specific steps:
s231 formation of target region the random radiation field is of the formula:
s232, in the microwave correlation forward-looking imaging process, the random radiation field irradiates the target, then reflects to form an echo, and is received by a receiving antenna to obtain a target echo E Hr、EVr as follows:
S233, obtaining target scattering characteristics and target polarization scattering characteristics through inversion:
s3, iterative imaging is performed by utilizing an orthogonal complement space of the random radiation field.
S31, if polarization information is not considered, TSVD processing is carried out on the space-time radiation field with insufficient randomness, then compression correlation operation is carried out on the echo and the space-time random radiation field, and finally a target scattering characteristic matrix is obtained, so that correlation imaging is realized. Equation (11) can be reduced to without regard to noise:
Er=Esσ (12)
S32, if the polarization scattering matrix of the target is taken into consideration, the polarization scattering matrix of the target is obtained while the target is imaged, the above formula solution cannot be obtained, and the characteristics of different HV antenna radiation fields must be taken into consideration for calculation. Solving the polarization scattering matrix of the target is equivalent to instantaneous polarization measurement, namely, can be understood as solving the instantaneous polarization scattering matrix of the target by using a microwave correlation imaging method.
The echo at this time can be expressed as:
EHr=EHsSHH+EVsSHV (13)
EVr=EHsSHV+EVsSVV (14)
if the above equation is to be solved, the best effect is that the radiation field of the H antenna is orthogonal to the radiation field of the V antenna, even if:
EVsEHs=0 (15)
However, for microwave correlated imaging, the radiation field irradiated to the target is randomly uncontrollable, and even if the radiation field is directly generated by using a mask or other tools, the HV radiation field at the target cannot be ensured to be mutually orthogonal, so that only another auxiliary matrix can be selected to deform the HV radiation field during calculation.
S32 comprises the following specific steps:
S321, obtaining a horizontal polarization radiation field orthogonal complementary matrix and a vertical polarization radiation field orthogonal complementary matrix through QR decomposition, and carrying out cross multiplication on the orthogonal complementary matrix and echo information;
s321 includes the steps of:
S3211, performing QR decomposition on a radiation field matrix of a target, and calculating to obtain orthogonal complement spaces E Hsb and E Vsb of the radiation field matrices E Hs and E Vs;
S3212, multiplying the target HV echo by the HV radiation field orthogonal complement space respectively to obtain:
EHsbEHr=EHsbEHsSHH+EHsbEVsSHV (16)
EVsbEHr=EVsbEHsSHH+EVsbEVsSHV (17)
EHsbEVr=EHsbEHsSVH+EHsbEVsSVV (18)
EVsbEVr=EVsbEHsSVH+EVsbEVsSVV (19)
the characteristics of the orthogonal complement matrix are as follows:
Namely:
s322, forming a new association type and iterative imaging by the radiation field after the matrix obtained in the step S321 is processed in the same way.
Equation (21) can be solved by solving equation (12) to measure the instantaneous polarization scattering matrix of the target while imaging the target.
In order to verify the effect of the present invention, the method proposed by the present invention is initially verified in combination with simulation. Fig. 3 shows a simulation of polarization information driven microwave correlation imaging of the present invention, it should be understood that the simulation described in fig. 3 is merely illustrative of the present invention and is not intended to limit the present invention. Assuming that the bandwidth b=1 GHz of the frequency hopping signal, the pulse time width T p =0.2 us, the carrier frequency f 0 =5 GHz, the sampling rate f s =10 GHz, the number of antennas is 24, the plane distance is 500m, and the imaging area is set to be a plane of 20×20 meters, simulation results and analysis can be obtained as follows.
Fig. 3 (a) shows the state of the simulation target, the polarization distribution of the target scattering is shown, the space distribution of the target scattering is shown in the row and column of the simulation target, different target types are shown through lines and dots, the polarization scattering matrix of the target is randomly generated through simulation each time, fig. 3 (b) shows the space correlation of the random radiation field of one pulse, the space-time radiation field correlation of one sampling is shown in fig. 3 (b), the whole space correlation of the radiation field is spike-shaped, the extremely high autocorrelation of the space-time radiation field is represented, the radiation field condition of microwave correlation imaging is met, fig. 3 (c) shows the microwave correlation imaging result, namely the transverse cutting condition during one-time microwave correlation imaging, the imaging precision and resolution of the space-time distribution are visually shown in fig. 3 (c), the microwave correlation imaging result driven by the whole polarization information is shown in fig. 3 (d), and the microwave correlation imaging basically reaches the expected comparison with the target.
The target polarization scattering matrix measurements are shown in table 1 below,
TABLE 1 target polarization scattering matrix measurement results
It can be seen from table 1 that the measured target polarization scattering matrix is slightly erroneous, but has been substantially consistent, indicating that this method can be implemented to measure the instantaneous polarization matrix while imaging the target.
Fig. 4 is a graph comparing polarized microwave correlation imaging of the present invention with conventional instantaneous polarization measurement PSL and I, and comparing the polarized microwave correlation imaging method driven by polarization information of the present invention with conventional instantaneous polarization measurement method in order to further embody the advantages of the method in terms of measurement of polarization at any time. The instantaneous polarization measurement mainly uses the difference of H, V antenna two signals, so the difference degree of the two signals is an important basis for comparing the effect of the instantaneous polarization measurement, the difference is often measured by the autocorrelation and the cross correlation of the two signals, giuli defines two parameters in the research to describe the two signals, namely Isolation (I) and peak sidelobe level ratio (Peak Sidelobe Level, PSL), wherein the Isolation I is used for describing the signal orthogonality degree by the cross correlation of the signals, and the peak sidelobe level ratio PSL is used for describing the influence of the autocorrelation function sidelobes of the signals on the polarization scattering matrix measurement, so the instantaneous polarization measurement performance is evaluated by comparing PSL and I of the method with positive and negative linear frequency modulation signals used for general instantaneous polarization measurement. PSL and I are specifically defined as follows:
In the above description, i and j represent HV two polarized signal transmitting channels, The method is an autocorrelation function of signals, two paths of radiation fields of the HV antenna are calculated according to the above formula under the conditions that the time width T=0.2 mu s and the bandwidth is from 0.5GHz to 2.5GHz, the obtained results are shown in fig. 4, the change of PSL and I along with the time-frequency product of the signals is shown in fig. 4, a line represents PSL-polarized microwave correlation imaging, b line represents I-polarized microwave correlation imaging, c line represents I-traditional transient polarization measurement, and d line represents PSL-traditional transient polarization measurement. As can be seen from the comparison of PSL and I in two modes of polarized microwave correlation imaging and traditional instantaneous polarization measurement in FIG. 4, when the time-frequency product is smaller (Bτ < 340), peak Sidelobe Level (PSL) and isolation (I) of two paths of signals of an HV antenna of the polarized microwave correlation imaging method are higher than those of positive and negative chirp signals used in a common method by about 13dB and 3dB on average, so that the polarized microwave correlation imaging method has better instantaneous polarization measurement effect in a frequency product section, and can also be used for instantaneous polarization measurement by utilizing the polarized microwave correlation method, which is different from the common instantaneous polarization measurement method, the PSL has weak correlation with the time-frequency product of signals, which is determined by the randomness of radiation fields in the microwave correlation imaging, the microwave correlation imaging itself has high requirements on the randomness of the radiation fields, and is in a highly random state in the time domain and the frequency domain, so that the change of the time-frequency product has small influence degree on the randomness of the radiation fields, and the change amplitude of the time-frequency product is smaller than 0.4dB in the total time-frequency product section (100 Bτ < 500), thereby having the advantages of instantaneous polarization measurement can be obtained under the aspect of more instantaneous polarization measurement than the traditional method.
In order to further embody the anti-interference capability of the polarized information driven microwave correlation imaging method of the invention, the anti-interference capability is compared with that of the traditional microwave correlation imaging method under the condition of the same parameter (SNR=10dB), fig. 5 shows an anti-interference effect comparison graph of the polarized microwave correlation imaging method of the invention and the traditional microwave correlation imaging method, the simulation result is fig. 5, fig. 5 (a) is a simulation target, fig. 5 (b) is a traditional microwave correlation imaging result, fig. 5 (c) is a polarized microwave correlation imaging result, fig. 5 (d) is an error comparison of the traditional microwave correlation imaging method and the polarized microwave correlation imaging method of the invention, wherein an e-curve of fig. 5 (d) represents an NMSE-dependent SNR-dependent curve of the traditional microwave correlation imaging method, and an f-curve represents an E-dependent SNR-dependent curve of the polarized microwave correlation imaging method of the invention, and as can be seen from fig. 5, the anti-interference capability of the polarized microwave correlation imaging method is stronger, imaging quality is higher because the correlation imaging part calculates a polarization scattering matrix of the target while imaging, and the correlation effect is equivalent to the calculation of the correlation result is prevented from being effectively influenced by multiple times. Meanwhile, the simulation does not consider the influence of the polarized antenna on noise interference, and in the practical application process, the anti-interference capability of the polarized antenna is only stronger as the polarized antenna can only receive signals in the preset polarization direction.
Finally, it is pointed out that while the invention has been described with reference to a specific embodiment thereof, it will be understood by those skilled in the art that the above embodiments are provided for illustration only and not as a definition of the limits of the invention, and various equivalent changes or substitutions may be made without departing from the spirit of the invention, therefore, all changes and modifications to the above embodiments shall fall within the scope of the appended claims.