[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN115586495A - MR-FDA-MIMO radar interference suppression method - Google Patents

MR-FDA-MIMO radar interference suppression method Download PDF

Info

Publication number
CN115586495A
CN115586495A CN202211119206.7A CN202211119206A CN115586495A CN 115586495 A CN115586495 A CN 115586495A CN 202211119206 A CN202211119206 A CN 202211119206A CN 115586495 A CN115586495 A CN 115586495A
Authority
CN
China
Prior art keywords
interference
snapshot data
virtual
array
suppression
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.)
Pending
Application number
CN202211119206.7A
Other languages
Chinese (zh)
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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN202211119206.7A priority Critical patent/CN115586495A/en
Publication of CN115586495A publication Critical patent/CN115586495A/en
Pending legal-status Critical Current

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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures

Landscapes

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

Abstract

The invention provides an MR-FDA-MIMO radar interference suppression method, which comprises the steps of transmitting signals according to a signal transmission model by using a transmitting array element, and acquiring multiple times of snapshot data received by a receiving array element; compensating the snapshot data to obtain compensated snapshot data; suppressing side lobe suppression type interference on the compensated snapshot data in a receiving dimension by utilizing a generalized side lobe cancellation GSC (generalized minimum cost) method to obtain snapshot data after side lobe interference suppression; solving a covariance matrix for the equivalent emission dimension snapshot data after sidelobe interference suppression, expanding the virtual snapshot freedom degree by utilizing the idea of a difference common matrix, and then suppressing multi-main lobe interference by utilizing a non-adaptive method to obtain a true target. Compared with the prior art, the invention has the advantages that the deception interference with the number higher than that of the transmitting array elements can be inhibited under the condition that the transmitting array elements and the receiving array elements are limited, the number of the inhibited interference is effectively increased, and the capacity of the frequency diversity array for inhibiting the number of the interference is improved.

Description

MR-FDA-MIMO radar interference suppression method
Technical Field
The invention belongs to the technical field of radar detection, and particularly relates to an MR-FDA-MIMO radar interference suppression method.
Background
With the development of information technology, the electromagnetic environment in a battlefield is more and more complex, the types and the number of electronic interferences are increased all the time, especially main lobe deceptive interferences, and the deceptive interferences radiate electromagnetic waves similar to target echoes to induce a radar system to treat false targets as true targets, so that the parameter estimation and tracking of the false targets on the true targets are influenced, and the battlefield viability of the radar is greatly threatened.
Indeed, conventional radars lack effective system freedom to distinguish true targets from false targets in the main lobe direction. Therefore, how to combat the deceptive jamming, especially the active deceptive jamming entering from the main lobe of the radar antenna, has become an urgent issue to be solved. Unlike a conventional phased array, a Frequency Diversity Array (FDA) radar uses small frequency increments in array elements to increase the degree of freedom in the distance dimension. In order to fully utilize the distance dimension information, FDA and Multiple Input Multiple Output (MIMO) technology are combined, and the method is favorable for always-mainlobe deceptive interference. In general, the number of main lobe deception jamming suppressions is fundamentally limited by the number of radar transmitting array elements.
In the prior art, the distance dimension degree of freedom of the FDA is utilized, a true target and deceptive interference are distinguished in a transmitting airspace frequency domain, and meanwhile, a robust distance angle two-dimensional beam forming method based on a direct data domain is provided. Likewise, the number of spoofed interference suppressions for multiple-input multiple-output (MIMO) radar is limited by the number of its transmitting elements.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides an MR-FDA-MIMO radar interference suppression method. The technical problem to be solved by the invention is realized by the following technical scheme:
the invention provides an MR-FDA-MIMO radar interference suppression method, which comprises the following steps:
step 1: establishing a signal transmitting model;
step 2: transmitting signals by using the transmitting array elements according to the signal transmitting model, and acquiring multiple times of snapshot data received by the receiving array elements;
the multi-time snapshot data comprises snapshot data of a true target, sidelobe pressing mode interference which is not greater than the number of the receiving array elements and main lobe deception interference which is greater than the number of the transmitting array elements;
and step 3: compensating the snapshot data to obtain compensated snapshot data;
and 4, step 4: suppressing side lobe suppression type interference on the compensated snapshot data in a receiving dimension by using a generalized side lobe cancellation GSC (generalized side lobe cancellation) method to obtain the snapshot data after side lobe interference suppression;
and 5: and (3) suppressing main lobe deceptive jamming on the snapshot data after the side lobe jamming is suppressed in an equivalent emission dimension by using a non-adaptive method, and obtaining virtual snapshot data of a real target.
The invention has the beneficial effects that:
the invention provides an MR-FDA-MIMO radar interference suppression method, which comprises the steps of transmitting signals according to a signal transmission model by using a transmitting array element, and acquiring multiple times of snapshot data received by a receiving array element; compensating the snapshot data to obtain compensated snapshot data; suppressing side lobe suppression type interference on the compensated snapshot data in a receiving dimension by using a generalized side lobe cancellation GSC (generalized side lobe cancellation) method to obtain the snapshot data after side lobe interference suppression;
and by using a non-adaptive method, the virtual degree of freedom of the transmission dimension is expanded, main lobe deceptive interference of the snapshot data after the sidelobe interference suppression is suppressed in the receiving dimension, and a true target is obtained. Compared with the prior art, the invention has the advantages that the deception interference with the number higher than that of the transmitting array elements can be inhibited under the condition that the transmitting array elements and the receiving array elements are limited, the number of the inhibited interference is effectively increased, and the capacity of the frequency diversity array for inhibiting the number of the interference is improved.
The present invention will be described in further detail with reference to the accompanying drawings and examples.
Drawings
Fig. 1 is a schematic flowchart of an MR-FDA-MIMO radar interference suppression method according to an embodiment of the present invention;
FIG. 2 is a diagram showing 4 array elements;
FIG. 3 is a MR-FDA-MIMO radar target and interference power spectrum;
FIG. 4 is a graph showing the results of side lobe interference suppression;
FIG. 5 is a graph of output signal-to-noise ratio (SINR) versus input signal-to-noise ratio (SNR) after 300 Monte Carlo tests;
FIG. 6 is a graph of multi-mainlobe deception jamming suppression results comparing a 7-array element Minimum Variance Distortionless Response (MVDR), a virtual array adaptive approach, and a virtual array non-adaptive approach;
FIG. 7 is a graph of signal-to-noise ratio for multiple spoofing interference;
FIG. 8 is a graph of the SNR loss performance of the 7-element Minimum Variance Distortionless Response (MVDR) method compared to the virtual array non-adaptive method for multiple decoy cases.
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.
As shown in fig. 1, an MR-FDA-MIMO radar interference suppression method provided by the present invention includes:
step 1: establishing a signal transmitting model;
wherein, the signal emission model is:
Figure RE-GDA0003975953580000031
where E is the transmit energy, M is the number of transmit array elements, T P Is the pulse width, f m =f 0 +ΔfD m Is the carrier frequency of the m-th antenna element, D m The position corresponding to the array element of the minimum redundant array is shown, and delta f < f 0 Phi is the envelope of the transmitted signal, the envelope of the transmitted signal satisfying
Figure RE-GDA0003975953580000032
M '=1,2, \ 8230, M ≠ M', τ isTime delay, m is the mth transmit array element, and m' is the mth transmit array element.
Step 2: transmitting signals by using the transmitting array elements according to the signal transmitting model, and acquiring multiple times of snapshot data received by the receiving array elements;
the multi-time snapshot data comprise snapshot data of a true target, sidelobe pressing mode interference which is not greater than the number of the receiving array elements and main lobe deception interference which is greater than the number of the transmitting array elements;
for a position in the far field (R) 00 ) The echo received by the nth receiving array element can be expressed as
Figure RE-GDA0003975953580000041
Under the assumption of narrow band,. Phi m (t-τ m,n )≈Φ m (t-τ 0 ) Beta is the scattering coefficient of the point target, tau m,n The time delay of the signal transmitted by the mth antenna and received by the nth array element can be written as
Figure RE-GDA0003975953580000042
Wherein is d R Spacing of receiving elements, τ 0 =2R 0 And c is the common delay time and c is the speed of light.
The signal is down-converted, matched filtered in the receiver. After the three signal processing flows, the received data can be expressed as the following parts:
Figure RE-GDA0003975953580000043
wherein
Figure RE-GDA0003975953580000044
λ 0 =c/f 0 Is the reference wavelength.
Thus, the snapshot data of true targets received in MR-FDA-MIMO is represented as
Figure RE-GDA0003975953580000045
Wherein,
Figure RE-GDA0003975953580000046
and
Figure RE-GDA0003975953580000047
respectively, a transmit steering vector and a receive steering vector, respectively, denoted as
Figure RE-GDA0003975953580000051
Figure RE-GDA0003975953580000052
Wherein, y MN The receiving array element receives the snapshot data of the real target after the Mth transmitting array element transmits signals,
Figure RE-GDA0003975953580000053
in order to guide the vector for the distance,
Figure RE-GDA0003975953580000054
an angle steering vector, \ indicates a Hadamard product, (R) 00 ) Indicating the position of a true object in the far field, c is the speed of light, λ 0 =c/f 0 Is a reference wavelength, d R The receive array element spacing, d the transmit element spacing, typically a half wavelength,
Figure RE-GDA0003975953580000055
is a complex set.
The transmitting frequency and the receiving frequency corresponding to the true target are respectively
Figure RE-GDA0003975953580000056
Figure RE-GDA0003975953580000057
Suppose Q are from theta q Q =1,2, \ 8230, and the interference signal in the Q direction has side-lobe suppression type interference of
Figure RE-GDA0003975953580000058
Wherein ξ q Is a uniform complex Gaussian random variable with a mean value of zero and a variance of
Figure RE-GDA0003975953580000059
E {. Is the mathematical expectation operation;
Figure RE-GDA00039759535800000510
and
Figure RE-GDA00039759535800000511
respectively representing the transmit and receive steering vectors of the interference of the respective noise classes, n, due to the statistical nature and the obtrusiveness of the noise of the interference T Satisfying a zero mean Gaussian noise distribution, i.e. n T ~CN(0,I M );
Suppose that in the far field of space (R) jj ) The interference machine intercepts radar signals and delays and forwards one or more pulse periods to form a false target, and the interference radar identifies and tracks the true target. The main lobe deception jamming contained in the received echo signal of the n array element is
Figure RE-GDA0003975953580000061
Wherein Q' is the number of decoys, beta q′ Is Q ' (Q ' =1, 2., Q ') Scattering coefficient of a false target, tau q′ Is the delay of the qth' th decoy, τ j Is a two-way time delay, the jammer is located in the far field at (R) jj ),
Figure RE-GDA0003975953580000062
The qth' th decoy can be expressed as
Figure RE-GDA0003975953580000063
Wherein
Figure RE-GDA0003975953580000064
The equivalent distance of the qth' th decoy is shown. The q' th decoy has a transmit and receive space-time frequency of
Figure RE-GDA0003975953580000065
Figure RE-GDA0003975953580000066
The snapshot data finally received by the receiving array element not only comprises the snapshot data of a true target, main lobe interference and side lobe interference, but also is expressed as
x=x s +x J +x I +x n (16)
Wherein x is n Is a mean of 0 and a variance of σ 2 White gaussian noise.
And 3, step 3: compensating the snapshot data to obtain compensated snapshot data;
the emission space frequency of the invention is distance-dependent, which needs further compensation, and the compensated snapshot data is
Figure RE-GDA0003975953580000067
Where p represents the number of false target delay pulses,
Figure RE-GDA0003975953580000068
for the purpose of the transmit-receive joint compensation vector,
Figure RE-GDA0003975953580000071
is a compensation vector of the emission domain, r a Is the primary value distance.
The transmission frequencies of the true target and the false target after compensation are respectively
Figure RE-GDA0003975953580000072
Figure RE-GDA0003975953580000073
Wherein R is u =c/(2f r ) Is the maximum unambiguous distance.
And 4, step 4: suppressing side lobe suppression type interference on the compensated snapshot data in a receiving dimension by utilizing a generalized side lobe cancellation GSC (generalized minimum cost) method to obtain snapshot data after side lobe interference suppression;
the side lobe interference is different from the target angle, and is suppressed by a generalized side lobe cancellation (GSC) method in a receiving dimension. The receiving dimension is the number dimension of the receiving array elements, and the snapshot data after the sidelobe interference suppression is
Figure RE-GDA0003975953580000074
Wherein w R =w 0 -w a ,w a =(C H RC) -1 (C H Rw 0 ),
Figure RE-GDA0003975953580000075
The covariance matrix is
Figure RE-GDA0003975953580000076
x H comp (l) Blocking matrix
Figure RE-GDA0003975953580000077
|θ-θ 0 | ≦ ε, θ is the range of the target that is likely to be in the main lobe, ε is the angle uncertainty parameter range.
And 5: and (3) suppressing main lobe deceptive jamming on the snapshot data after the side lobe jamming is suppressed in an equivalent emission dimension by using a non-adaptive method, and obtaining the snapshot data of the true target.
False target suppression is carried out, and the method can be divided into two steps: (1) Expanding equivalent transmitting freedom according to the idea of difference common array; (2) Adopting a non-adaptive beam forming method to restrain a false target, specifically:
step 5-1: obtaining the covariance matrix R of the snapshot data after the sidelobe interference suppression in the equivalent transmission dimension t
Wherein, covariance matrix R t Is composed of
Figure RE-GDA0003975953580000081
Wherein K is the number of information sources, L is the number of snapshots, x t And (5) carrying out snapshot on equivalent transmission dimension data after sidelobe interference suppression. l is the first snapshot, a k For the k-th source steering vector,
Figure RE-GDA0003975953580000082
is the noise power, I M Is an M-dimensional identity matrix.
The equivalent virtual signal of the minimum redundant array is constructed by the second order statistic calculated by equation (19). The minimum redundancy array is exemplified by four array elements, the covariance matrix R of which t Is composed of
Figure RE-GDA0003975953580000083
Thus the covariance matrix R t The ith row and jth column elements of (1) are:
Figure RE-GDA0003975953580000084
wherein i, j =1,2, \ 8230;, M, f comp Is the transmission frequency after the true and false target compensation.
Second order statistic r i,j The exponential function term in (1) can be regarded as a guide vector of an equivalent virtual array element, and the position of each element in the virtual array element is as follows:
Figure RE-GDA0003975953580000085
then
Figure RE-GDA0003975953580000086
Each individual element in (a) may be equivalently represented using a differential co-array of the lowermost redundant array,
namely:
Figure RE-GDA0003975953580000087
the differential common array V of the minimum redundant array at least comprises 2D M With M physical array elements, 2D can be obtained M An equivalent virtual array element. Thus, for covariance matrix R t Vectorization is carried out, and equivalent virtual signals corresponding to the virtual array V can be obtained.
Step 5-2: the covariance matrix R t Each row of the virtual array signals is superposed into a vector to obtain an equivalent virtual array signal z;
wherein the virtual array signal z is
Figure RE-GDA0003975953580000091
Wherein,
Figure RE-GDA0003975953580000092
is a virtual steering vector matrix corresponding to the virtual domain V,
Figure RE-GDA0003975953580000093
a virtual guide vector is generated by a virtual guide vector,
Figure RE-GDA0003975953580000094
contains K sources, e = vec (I) M )。
Step 5-3: utilizing the G matrix to remove redundancy and rearrange the virtual array signal z to obtain a rearranged virtual array signal
Figure RE-GDA0003975953580000095
The values indicate that: the virtual array signal z contains redundancy, the data is subjected to redundancy removal and rearrangement by using the G matrix, and the rearranged virtual array signal
Figure RE-GDA0003975953580000096
Is composed of
Figure RE-GDA0003975953580000097
Wherein,
Figure RE-GDA0003975953580000098
a matrix of virtual steering vectors is formed,
Figure RE-GDA0003975953580000099
the virtual guide vector is a vector of the virtual guide,
Figure RE-GDA00039759535800000910
according to the difference of the minimum redundant array element number, the structure of the selection matrix G is different according to the difference of the covariance matrix element position. Taking 4 array elements as an example, G can be expressed as
Figure RE-GDA00039759535800000911
Step 5-4: predicting a mainlobe deceptive interference location in a transmit frequency dimension to construct a covariance matrix R containing only virtual interference and noise v
Since the main lobe deceptive interference can be predicted in the transmit frequency dimension, a virtual interference can be constructed.
The covariance matrix thus constructed is expressed as
Covariance matrix R v Is composed of
Figure RE-GDA0003975953580000101
Wherein R is N And R I For the constructed noise covariance matrix and interference covariance matrix,
Figure RE-GDA0003975953580000102
in order to be able to measure the power of the noise,
Figure RE-GDA0003975953580000103
in order to be a virtual interference power,
Figure RE-GDA0003975953580000104
and a steering vector corresponding to the q virtual interference.
Step 5-5: the covariance matrix R is calculated using the Minimum Variance Distortionless Response (MVDR) criterion v A non-adaptive weight vector of (a);
according to the Minimum Variance Distortionless Response (MVDR) criterion, the non-adaptive weight vector is
Figure RE-GDA0003975953580000105
Wherein,
a v0 )=[exp(-j2πd sin(θ 0 )D M0 ),…,1,…,exp(j2πd sin(θ 0 )D M0 )]is a true targetA virtual steering vector;
and 5-6: using the non-adaptive weight vector to the rearranged virtual array signal
Figure RE-GDA0003975953580000107
And carrying out self-adaptive main lobe deception jamming suppression in a receiving dimension to obtain virtual snapshot data of the true target.
Wherein the true target is
Figure RE-GDA0003975953580000106
Where H denotes the transpose of the matrix.
Simulation experiments can further prove the beneficial effects of the invention.
Simulation experiment I:
first, table 1 provides the array element positions of the minimum redundant array, as shown in fig. 2, and is illustrated by taking 4 array elements as an example, for example, the array elements of the quaternary uniform array are located at {0,1,4,6} d, and the array elements of the quaternary MRA are located at {0,1,4,6} d, where d represents the unit pitch, typically half wavelength, and the array of the four-array element MRA is shown in fig. 1.
TABLE 1 MRA array element distribution
Figure RE-GDA0003975953580000111
The simulation sets that the number of the transmitting array elements and the number of the receiving array elements are 7, the target angle is 0 degrees, the SNR =20dB, the side lobe pressing type interference angles are 30 degrees, -30 degrees, the dry-to-noise ratio JNRn =40dB, the pulse repetition period PRF =10kHz, and the reference carrier frequency f 0 =16GHz, wavelength λ 0 =0.0187m, assuming there are 2 decoys, the dry-to-noise ratio INR =20db, INR =25db.
FIG. 3 is a MR-FDA-MIMO radar target and interference power spectrogram. As can be seen from the graph (b) in fig. 3, true and false targets can be distinguished in the transmit dimension, and false targets can be suppressed in the transmit dimension. As can be seen from the (c) diagram in fig. 3, the suppressed interference is distinguishable in the receiving dimension, and can be suppressed in the receiving dimension.
Fig. 4 is a diagram of a side lobe suppression result, which is distorted by a Minimum Variance Distortionless Response (MVDR) method because of a target signal component existing in snapshot sample data, and this greatly affects the performance of suppressing side lobe suppression interference. While the Generalized Sidelobe Canceling (GSC) method can effectively suppress the suppressed interference, with performance close to the ideal case with an accurate interference plus noise covariance matrix.
Fig. 5 shows a plot of output signal-to-noise ratio (SINR) versus input signal-to-noise ratio (SNR) after 300 monte carlo tests. As can be seen from the graph (a) in fig. 5, the Generalized Sidelobe Cancellation (GSC) method effectively suppresses sidelobe interference, and the performance approaches the upper bound in the case of a relatively small input signal-to-noise ratio. When the signal-to-noise ratio is greater than 20dB, the signal-to-noise ratio output of both generalized side lobe cancellation (GSC) and Minimum Variance Distortionless Response (MVDR) decreases due to the presence of the object in the sample. In particular, generalized Sidelobe Cancellation (GSC) performs better than Minimum Variance Distortionless Response (MVDR). As shown in fig. 5 (b), both the Generalized Sidelobe Cancellation (GSC) and Minimum Variance Distortionless Response (MVDR) methods are ineffective in the presence of mainlobe spoofing interference. Note that the Generalized Sidelobe Cancellation (GSC) curve differs from the ideal curve by approximately 27dB, which is exactly the sum of the dry-to-noise ratios of the two main lobe spoofing interferers.
A second simulation experiment: assuming that there are 14 decoys, the interference-to-noise ratio INR (1-7) =20db, INR (2-14) =25dB. The rest parameters are the same as those of the simulation experiment.
Fig. 6 compares the multi-mainlobe deception jamming suppression results of 7-array element Minimum Variance Distortionless Response (MVDR), virtual array adaptive method, and virtual array non-adaptive method. The 7-array Minimum Variance Distortionless Response (MVDR) method has only 7 array elements, so 14 decoys cannot be suppressed. The virtual array self-adaptive method is only
Figure RE-GDA0003975953580000121
One snapshot sample data has errors in covariance estimation, and a plurality of false targets cannot be suppressed. The virtual array non-adaptive method is not influenced by the number of snapshot samplesWith sufficient degrees of freedom, multiple spurious objects can be suppressed.
In order to comprehensively evaluate the performance of the MR-FDA-MIMO radar, the variation curve of the output signal-to-interference-and-noise ratio with the input signal-to-noise ratio under the condition of multiple deceptive interferences is simulated. As can be seen from fig. 7, when there are 14 decoy targets, the output signal-to-noise ratio performance of the 7-element Minimum Variance Distortionless Response (MVDR) method is severely degraded, below the ideal curve of 32dB, because the number of spoofing interferers is greater than the number of transmitting elements. For a virtual array non-adaptive method, the degree of freedom is expanded, the number of virtual array elements is greater than the number of deception interferences, a plurality of interferences are completely inhibited, and the performance is close to an ideal curve.
Fig. 8 shows the signal to interference plus noise ratio (snr) loss performance comparison of the 7-array element Minimum Variance Distortionless Response (MVDR) method and the virtual array non-adaptive method for the multi-decoy case. It can be seen that the 7-array Minimum Variance Distortionless Response (MVDR) method loses robustness completely when the number of interferers is greater than the number of array elements. The virtual array non-adaptive method is not affected, a plurality of interferences are completely inhibited, and a performance curve is close to an ideal curve.
Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of the present invention, "a plurality" means two or more unless specifically defined otherwise.
While the present application has been described in connection with various embodiments, other variations to the disclosed embodiments can be understood and effected by those skilled in the art in practicing the claimed application, from a review of the drawings, the disclosure, and the appended claims. In the claims, the word "comprising" does not exclude other elements or steps, and the word "a" or "an" does not exclude a plurality.
The foregoing is a further detailed description of the invention in connection with specific preferred embodiments and it is not intended to limit the invention to the specific embodiments described. 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.

Claims (10)

1. An MR-FDA-MIMO radar interference suppression method, comprising:
step 1: establishing a signal transmitting model;
step 2: transmitting signals by using the transmitting array elements according to the signal transmitting model, and acquiring multiple times of snapshot data received by the receiving array elements;
the multi-time snapshot data comprise snapshot data of a true target, sidelobe pressing mode interference which is not greater than the number of the receiving array elements and main lobe deception interference which is greater than the number of the transmitting array elements;
and step 3: compensating the snapshot data to obtain compensated snapshot data;
and 4, step 4: suppressing side lobe suppression type interference on the compensated snapshot data in a receiving dimension by utilizing a generalized side lobe cancellation GSC (generalized minimum cost) method to obtain snapshot data after side lobe interference suppression;
and 5: and (3) suppressing main lobe deceptive jamming on the snapshot data after the side lobe jamming suppression in an equivalent transmission dimension by using a non-adaptive method, and obtaining virtual snapshot data of the true target.
2. The method of claim 1, wherein the signal transmission model in step 1 is:
Figure RE-FDA0003975953570000011
where E is the transmit energy, M is the number of transmit array elements, T P Is the pulse width, f m =f 0 +ΔfD m Is the carrier frequency of the m-th antenna element, D m The position corresponding to the array element of the minimum redundant array is shown, and delta f < f 0 Phi is the envelope of the transmitted signal, the envelope of the transmitted signal satisfying
Figure RE-FDA0003975953570000012
m′=1,2, Mm ≠ m ', τ being the time delay, m being the mth transmitting array element, m' being the mth transmitting array element.
3. The MR-FDA-MIMO radar interference mitigation method according to claim 2, wherein the snapshot data of the true target is represented as
Figure RE-FDA0003975953570000013
Wherein,
Figure RE-FDA0003975953570000014
and
Figure RE-FDA0003975953570000015
respectively, transmit and receive steering vectors, denoted as
Figure RE-FDA0003975953570000021
Figure RE-FDA0003975953570000022
Wherein, y MN The receiving array element receives the snapshot data of the true target after the Mth transmitting array element transmits the signal,
Figure RE-FDA0003975953570000023
in order to guide the vector for the distance,
Figure RE-FDA0003975953570000024
is an angle steering vector, < i > represents a Hadamard product, < i > R 00 ) Indicating the position of a true object in the far field, c is the speed of light, λ 0 =c/f 0 Is a reference wavelength, d R The receive array element spacing, d the transmit element spacing, typically a half wavelength,
Figure RE-FDA0003975953570000025
is a complex set;
the transmitting frequency and the receiving frequency corresponding to the true target are respectively
Figure RE-FDA0003975953570000026
Figure RE-FDA0003975953570000027
The side lobe suppression type interference is
Figure RE-FDA0003975953570000028
Wherein ξ q Is a uniform complex Gaussian random variable with a mean value of zero and a variance of
Figure RE-FDA0003975953570000029
E {. Is the mathematical expectation operation;
Figure RE-FDA00039759535700000210
and
Figure RE-FDA00039759535700000211
respectively representing the transmit and receive steering vectors of the interference of the respective noise classes, n, due to the statistical nature and the obtrusiveness of the noise of the interference T Satisfying zero mean value Gaussian noiseAcoustic distribution, i.e. n T ~CN(0,I M );
The main lobe deceptive jamming is
Figure RE-FDA00039759535700000212
Wherein Q' is the number of decoys, beta q′ Is the scattering coefficient, τ, of the Q ' (Q ' =1, 2.., Q ') th false target q 'is the delay of the q' th decoy, τ j Is a two-way time delay, the jammer is located in the far field at (R) jj ),
The snapshot data in the step 2 is
x=x s +x J +x I +x n
Wherein x is n Is a mean of 0 and a variance of σ 2 White gaussian noise.
4. The MR-FDA-MIMO radar interference mitigation method according to claim 3, wherein the snapshot data compensated in step 3 is
Figure RE-FDA0003975953570000031
Wherein p represents the number of false target delay pulses,
Figure RE-FDA0003975953570000032
for the purpose of a transmit-receive joint compensation vector,
Figure RE-FDA0003975953570000033
is a compensation vector of the emission domain, r a The distance is a distance of a main value,
the transmission frequencies of the true target and the false target after compensation are respectively
Figure RE-FDA0003975953570000034
Figure RE-FDA0003975953570000035
Wherein R is u =c/(2f r ) Is the maximum unambiguous distance.
5. The MR-FDA-MIMO radar interference suppression method according to claim 4, wherein the snapshot data after sidelobe interference suppression in the step 4 is
Figure RE-FDA0003975953570000036
Wherein w R =w 0 -w a ,w a =(C H RC) -1 (C H Rw 0 ),
Figure RE-FDA0003975953570000041
The covariance matrix is
Figure RE-FDA0003975953570000042
Blocking matrix
Figure RE-FDA0003975953570000043
Where | theta-theta 0 | ≦ ε, θ is the range of the target that is likely to be in the main lobe, ε is the angle uncertainty parameter range.
6. The MR-FDA-MIMO radar interference suppression method according to claim 5, wherein the step 5 comprises:
step 5-1: obtaining the covariance matrix R of the snapshot data after the sidelobe interference suppression in the equivalent transmission dimension t
Step 5-2: the covariance matrix R t Are superimposed into a vectorObtaining an equivalent virtual array signal z;
step 5-3: utilizing the G matrix to remove redundancy and rearrange the virtual array signal z to obtain a rearranged virtual array signal
Figure RE-FDA0003975953570000044
Step 5-4: predicting a mainlobe deceptive interference location in a transmit frequency dimension to construct a covariance matrix R containing only virtual interference and noise v
Step 5-5: the covariance matrix R is calculated using the Minimum Variance Distortionless Response (MVDR) criterion v A non-adaptive weight vector of (a);
and 5-6: using the non-adaptive weight vector to the rearranged virtual array signal
Figure RE-FDA0003975953570000045
And carrying out self-adaptive main lobe deception jamming suppression in the emission dimension to obtain a true target.
7. The MR-FDA-MIMO radar interference mitigation method according to claim 6, wherein the covariance matrix R in step 5-1 t Is composed of
Figure RE-FDA0003975953570000046
Wherein K is the number of information sources, L is the number of snapshots, x t For suppressing the number of the equivalent emission dimension snapshots after the side lobe interference system, l is the ith snapshot, a k The steering vector of the k-th source,
Figure RE-FDA0003975953570000047
is the power of the k-th source,
Figure RE-FDA0003975953570000051
is the noise power, I M Is an M-dimensional identity matrix.
8. The method of claim 7, wherein the method further comprises the step of,
the virtual array signal z in the step 5-2 is
Figure RE-FDA0003975953570000052
Wherein,
Figure RE-FDA0003975953570000053
is a virtual steering vector matrix corresponding to the virtual field V,
Figure RE-FDA0003975953570000054
a virtual guide vector is generated by a virtual guide vector,
Figure RE-FDA0003975953570000055
contains K sources, e = vec (I) M ),
Figure RE-FDA0003975953570000056
Is a real number set;
the rearranged virtual array signal in the step 5-3
Figure RE-FDA0003975953570000057
Is composed of
Figure RE-FDA0003975953570000058
Wherein,
Figure RE-FDA0003975953570000059
a matrix of virtual steering vectors is formed,
Figure RE-FDA00039759535700000510
virtual directorThe amount of the compound (A) is,
Figure RE-FDA00039759535700000511
9. the method as claimed in claim 8, wherein the covariance matrix R in step 5-4 is v Is composed of
Figure RE-FDA00039759535700000512
Wherein R is N And R I For the constructed noise covariance matrix and interference covariance matrix,
Figure RE-FDA00039759535700000513
in order to be able to measure the power of the noise,
Figure RE-FDA00039759535700000514
in order to be a virtual interference power,
Figure RE-FDA00039759535700000515
and a steering vector corresponding to the qth virtual interference.
10. The method of claim 9, wherein the non-adaptive weight vector in step 5-5 is
Figure RE-FDA00039759535700000516
Wherein,
a v0 )=[exp(-j2πd sin(θ 0 )D M0 ),…,1,…,exp(j2πd sin(θ 0 )D M0 )]a virtual steering vector that is a true target;
the true target in the step 5-6 is
Figure RE-FDA0003975953570000061
Where H denotes the transpose of the matrix.
CN202211119206.7A 2022-09-13 2022-09-13 MR-FDA-MIMO radar interference suppression method Pending CN115586495A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211119206.7A CN115586495A (en) 2022-09-13 2022-09-13 MR-FDA-MIMO radar interference suppression method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211119206.7A CN115586495A (en) 2022-09-13 2022-09-13 MR-FDA-MIMO radar interference suppression method

Publications (1)

Publication Number Publication Date
CN115586495A true CN115586495A (en) 2023-01-10

Family

ID=84773430

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211119206.7A Pending CN115586495A (en) 2022-09-13 2022-09-13 MR-FDA-MIMO radar interference suppression method

Country Status (1)

Country Link
CN (1) CN115586495A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116359857A (en) * 2023-06-02 2023-06-30 中国人民解放军空军预警学院 Space-time-frequency self-adaptive main lobe deception jamming prevention method and device for airborne early warning radar

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116359857A (en) * 2023-06-02 2023-06-30 中国人民解放军空军预警学院 Space-time-frequency self-adaptive main lobe deception jamming prevention method and device for airborne early warning radar
CN116359857B (en) * 2023-06-02 2023-09-01 中国人民解放军空军预警学院 Space-time-frequency self-adaptive main lobe deception jamming prevention method and device for airborne early warning radar

Similar Documents

Publication Publication Date Title
Lan et al. Suppression approach to main-beam deceptive jamming in FDA-MIMO radar using nonhomogeneous sample detection
CN109507649B (en) Method for resisting main lobe deception jamming of wave diversity array radar
CN104297734B (en) Deceiving interference suppressing method based on the MIMO radar of frequency diversity array
Melvin A stap overview
Wen et al. Enhanced three-dimensional joint domain localized STAP for airborne FDA-MIMO radar under dense false-target jamming scenario
Xu et al. Joint range and angle estimation using MIMO radar with frequency diverse array
US6252540B1 (en) Apparatus and method for two stage hybrid space-time adaptive processing in radar and communication systems
Yu et al. Polarimetric multiple-radar architectures with distributed antennas for discriminating between radar targets and deception jamming
CN115586495A (en) MR-FDA-MIMO radar interference suppression method
Yu et al. Mainbeam deceptive jammer suppression with joint element-pulse phase coding
Wu et al. Joint Beamforming Design and Power Control Game for a MIMO Radar System in the presence of Multiple Jammers
Wu et al. FDA-MRMO radar interference suppression
Adve et al. Adaptive space/frequency processing for distributed aperture radars
Qiu et al. Designing Constant Modulus Approximate Binary Phase Waveforms for Multi-Target Detection in MIMO Radar using LSTM Networks
CN106054142B (en) A kind of airborne MIMO radar main lobe smart munition suppressing method and system
CN110146854B (en) Robust anti-interference method for FDA-MIMO radar
CN113030877A (en) Precise intermittent sampling interference method for canceling narrow-band self-adaptive side lobe
Vijaykumar Mahamuni Space-time adaptive processing techniques (STAP) for mitigation of jammer interference and clutter suppression in airborne radar systems: A MATLAB implementation-based study
Mecca et al. Slow-time MIMO spacetime adaptive processing
CN113030874B (en) Close-range main lobe interference suppression method based on frequency control array radar
Mahamuni Space-Time Adaptive Processing (STAP) Techniques for Mitigation of Jammer Interference and Clutter Suppression in Airborne Radar Systems: A MATLAB Implementation-Based Study
Wu et al. Interference Suppression with MR-FDA-MIMO Radar Using Virtual Samples-based Beamforming
Zhang et al. Mainlobe Deceptive Jammer Suppression based on Secondary Frequency Compensation in FDA-MIMO Radar
Liu et al. Fast Repeater Deception Jamming Suppression with DFI-FDA-MIMO Radar
Wu et al. A Robust Interference Suppression Method Based on FDA-MIMO Radar

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