CN117784026B - Space-time-frequency domain combined active anti-composite interference method and device - Google Patents
Space-time-frequency domain combined active anti-composite interference method and device Download PDFInfo
- Publication number
- CN117784026B CN117784026B CN202410208694.1A CN202410208694A CN117784026B CN 117784026 B CN117784026 B CN 117784026B CN 202410208694 A CN202410208694 A CN 202410208694A CN 117784026 B CN117784026 B CN 117784026B
- Authority
- CN
- China
- Prior art keywords
- interference
- data
- channel data
- representing
- array
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 239000002131 composite material Substances 0.000 title claims abstract description 54
- 230000001629 suppression Effects 0.000 claims abstract description 54
- 230000035485 pulse pressure Effects 0.000 claims abstract description 38
- 230000003044 adaptive effect Effects 0.000 claims abstract description 27
- 238000005259 measurement Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 18
- 230000008859 change Effects 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 6
- 230000002401 inhibitory effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 20
- 238000005070 sampling Methods 0.000 description 8
- 238000001228 spectrum Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000002452 interceptive effect Effects 0.000 description 4
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention relates to the technical field of signal processing, and provides a space-time-frequency domain combined active anti-composite interference method and device. The method comprises the following steps: finding a distance gate unit where the slice interference is located from the channel data, constructing an adaptive weight for multi-channel array data in the distance gate unit where the slice interference is located, and obtaining adaptive channel data by using the adaptive weight; pulse pressure is carried out on the self-adaptive channel data, and a starting position R t of a target echo signal is found out from the pulse pressure; selecting first signal data from the starting position R t from the sum and difference channel data, and performing interference suppression on the first signal data in a fractional order domain; and then performing sum and difference angle measurement to obtain a target angle. The invention utilizes the difference of the composite interference and the target echo in three different dimensions of a space domain, a time domain and a fractional order domain to inhibit the composite interference, thereby effectively inhibiting the composite interference of main lobe slice interference in a pulse or between pulses and obtaining an accurate target angle.
Description
Technical Field
The invention relates to the technical field of signal processing, in particular to a space-time-frequency domain combined active anti-composite interference method and device.
Background
Some researches on slice interference suppression have been performed in the prior art, but the following problems still exist in the prior art: firstly, only certain slice interference suppression is considered, and a scene that an actual interference environment is generally composite interference is not considered; secondly, most models of the existing slice interference suppression method are built on the basis of a single-channel model, and spatial correlation among main stream phased array radar antenna units is not considered. Therefore, the effect of suppressing the slice interference is poor in practical use.
In view of this, overcoming the drawbacks of the prior art is a problem to be solved in the art.
Disclosure of Invention
The invention aims to provide a space-time-frequency domain combined active anti-composite interference method.
The invention adopts the following technical scheme:
In a first aspect, the present invention provides a space-time-frequency domain combined active anti-composite interference method, wherein a radar transmits inter-pulse combined agile pulses, the method includes:
from radar received multi-channel array data Generating sum channel data and difference channel data;
Finding a distance gate unit where the slice interference is located from the sum channel data, constructing an adaptive weight for multi-channel array data in the distance gate unit where the slice interference is located, and multiplying the adaptive weight with the multi-channel array data to obtain adaptive channel data;
pulse pressure is carried out on the self-adaptive channel data to obtain self-adaptive channel data after pulse pressure, and a position corresponding to a power peak value is found from the self-adaptive channel data after pulse pressure, namely, a starting position R t of a target echo signal;
Selecting first signal data with the length of a preset length from a starting position R t from the sum channel data, and performing interference suppression on the first signal data in a fractional order domain to obtain pulse pressure and channel data after interference suppression;
selecting second signal data with the length of a preset length from a starting position R t from the difference channel data, and performing interference suppression on the second signal data in a fractional order domain to obtain pulse pressure difference channel data after interference suppression;
And carrying out sum and difference angle measurement according to the interference-suppressed pulse pressure, a third position corresponding to the power peak value in the channel data and a fourth position corresponding to the power peak value in the interference-suppressed pulse pressure difference channel data to obtain a target angle.
Preferably, the multi-channel array data received according to radarGenerating sum channel data and difference channel data, specifically including:
Computing the multi-channel array data Angle of spatial source/>; Wherein/>Representing conjugate transpose operation,/>Representing airspace angle;
According to the angle of the space information source Judging whether the side lobe interference exists or not, and if so, judging that the side lobe interference exists, and according to the angle/>, of the side lobe interferenceConstruction of an orthographic projection matrix/>; Wherein/>The unit matrix representing M multiplied by M has the elements on opposite corner lines of all 1 and the elements at other positions of all 0;
Using the orthogonal projection matrix For the multichannel array data/>Performing sidelobe interference suppression to obtain data/>, after the sidelobe interference suppression；
Using the sidelobe interference suppressed dataSum channel data and difference channel data are generated.
Preferably, said using said sidelobe interference suppressed dataGenerating sum channel data and difference channel data, specifically including:
using the sidelobe interference suppressed data Formed sum channel data is/>；
Using the sidelobe interference suppressed dataThe resulting bad channel data is/>; Wherein/>Is a vector of length M x 1, the first M/2 elements being 1 and the last element being-1.
Preferably, the distance gate unit where the slice interference is found from the sum channel data specifically includes:
finding a first position and a second position with abrupt change of power in the sum channel data, and judging whether the average power between the first position and the second position is higher than a preset power or not;
And if the average power between the first position and the second position is higher than the preset power, taking the area between the first position and the second position as a distance gate unit where the slice interference exists.
Preferably, the constructing the adaptive weight for the multichannel array data in the range gate unit where the slice interference exists specifically includes:
Selecting multichannel array data in a range gate unit where slice interference is located ; Wherein,Sta is the starting position of the distance gate unit where the slice interference is located, and end is the ending position of the distance gate unit where the slice interference is located;
According to Construction of adaptive weights/>。
Preferably, the performing interference suppression on the first signal data in the fractional order domain to obtain pulse pressure and channel data after interference suppression specifically includes:
And calculating the average power of the first signal data, setting zero at a position higher than the average power to obtain processed first signal data, converting the processed first signal data into a fractional order domain, performing interference suppression on the processed first signal data in the fractional order domain, and performing FrFT inverse conversion to obtain pulse pressure and channel data after interference suppression.
Preferably, the array of the radar is a uniform linear array composed of M array elements, and the signals incident to the array comprise k t target echo signals, k jm main lobe slice interferences and k js side lobe noise suppression interferences. The guiding vector of the linear array is; Wherein/>Represents the direction of arrival of the source, d represents the array element spacing,Representing wavelength,/>Representing a transpose operation.
Preferably, the multi-channel array data; Wherein,，，；/>、/>And/>DOA,/>, respectively representing the kth t target, the kth jm main lobe slice disturbance and the kth js side lobe noise suppression disturbance、/>AndRespectively representing a target echo signal, a main lobe slice interference signal and a side lobe noise suppression interference signal received by the array,/>Representing a complex matrix of x rows and y columns, L representing the number of snapshots and N representing noise subject to gaussian distribution.
Preferably, the data is received in a multi-channel array according to radarBefore the sum channel data and the difference channel data are generated, band-pass filtering is further carried out on radar receiving signals so as to filter interference in a non-working bandwidth of the radar in the current period, and multichannel array data/>, is obtained。
In a second aspect, the present invention further provides a space-time-frequency domain combined active anti-composite interference device, for implementing the space-time-frequency domain combined active anti-composite interference method in the first aspect, where the device includes:
At least one processor; and a memory communicatively coupled to the at least one processor; the memory stores instructions executable by the at least one processor, where the instructions are executable by the processor to perform the space-time-frequency domain joint active anti-composite interference method of the first aspect.
In a third aspect, the present invention further provides a non-volatile computer storage medium, where computer executable instructions are stored, where the computer executable instructions are executed by one or more processors to perform the space-time-frequency domain joint active anti-composite interference method according to the first aspect.
According to the method, the composite interference is restrained by utilizing the difference of the composite interference and the target echo in three different dimensions of a space domain, a time domain and a fractional order domain, so that the composite interference including main lobe slice interference in a pulse or among pulses can be restrained effectively, and an accurate target angle can be obtained.
Drawings
In order to more clearly illustrate the technical solution of the embodiments of the present invention, the drawings that are required to be used in the embodiments of the present invention will be briefly described below. It is evident that the drawings described below are only some embodiments of the present invention and that other drawings may be obtained from these drawings without inventive effort for a person of ordinary skill in the art.
Fig. 1 is a schematic flow chart of a space-time-frequency domain combined active anti-composite interference method provided by an embodiment of the invention;
fig. 2 is a schematic diagram of spatial spectrum in a space-time-frequency domain combined active anti-composite interference method according to an embodiment of the present invention;
Fig. 3 is a schematic diagram of sum channel data after side lobe interference is suppressed in a space-time-frequency domain combined active anti-composite interference method provided by an embodiment of the present invention;
fig. 4 is a time-frequency diagram of a target echo signal in a practical application scenario according to an embodiment of the present invention;
fig. 5 is a time-frequency diagram of cross-period intermittent sampling direct interference forwarding in a first space-time-frequency domain combined active anti-composite interference method provided by an embodiment of the present invention;
fig. 6 is a time-frequency diagram of cross-period intermittent sampling direct interference forwarding in a second space-time-frequency domain combined active anti-composite interference method according to an embodiment of the present invention;
fig. 7 is a time-frequency diagram of cross-period intermittent sampling direct interference forwarding in a third space-time-frequency domain combined active anti-composite interference method according to an embodiment of the present invention;
Fig. 8 is a time-frequency diagram of cross-period intermittent sampling direct interference forwarding in a fourth space-time-frequency domain combined active anti-composite interference method according to an embodiment of the present invention;
Fig. 9 is a flow chart of another method for space-time-frequency domain joint active anti-composite interference according to an embodiment of the present invention;
fig. 10 is a schematic diagram of adaptive channel data after pulse pressure in a space-time-frequency domain combined active anti-composite interference method according to an embodiment of the present invention;
fig. 11 is a schematic diagram of sum channel data after all interference is suppressed in a space-time-frequency domain combined active anti-composite interference method according to an embodiment of the present invention;
Fig. 12 is a schematic diagram of difference channel data after all interference is suppressed in a space-time-frequency domain combined active anti-composite interference method according to an embodiment of the present invention;
Fig. 13 is a schematic diagram of pulse pressure and channel data after all interference is suppressed in a space-time-frequency domain combined active anti-composite interference method according to an embodiment of the present invention;
fig. 14 is a schematic diagram of channel data of pulse pressure difference after all interference is suppressed in a space-time-frequency domain combined active anti-composite interference method according to an embodiment of the present invention;
FIG. 15 is a schematic diagram of a target angle estimation result of a Monte Carlo simulation method according to an embodiment of the present invention;
fig. 16 is a schematic diagram of an architecture of a space-time-frequency domain combined active anti-composite interference device according to an embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
The terms "first," "second," and the like herein 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 defining "a first", "a second", etc. may explicitly or implicitly include one or more such feature. In the description of the present application, unless otherwise indicated, the meaning of "a plurality" is two or more.
In addition, the technical features of the embodiments of the present invention described below may be combined with each other as long as they do not collide with each other.
Example 1:
The embodiment 1 of the invention provides a space-time-frequency domain combined active anti-composite interference method, which is used for transmitting pulse of inter-pulse combined agility in a pulse by a radar, wherein the inter-pulse agility is realized by transmitting sub-pulses of a plurality of linear frequency modulation (Linear Frequency Modulation, abbreviated as LFM) signals with different slopes, the inter-pulse agility is realized by transmitting different carrier frequencies, and the frequency spectrum of the transmitted signals in adjacent heavy weeks is ensured not to overlap, as shown in figure 1, and the method comprises the following steps:
in step 201, multi-channel array data received from radar Sum channel data and difference channel data are generated.
In step 202, a distance gate unit where the slice interference is located is found from the sum channel data, an adaptive weight is constructed on the multi-channel array data in the distance gate unit where the slice interference is located, and the adaptive weight is multiplied by the multi-channel array data to obtain adaptive channel data.
In step 203, pulse pressure is performed on the adaptive channel data to obtain adaptive channel data after pulse pressure, and a position corresponding to a power peak is found from the adaptive channel data after pulse pressure, that is, the starting position R t of the target echo signal.
In step 204, from the sum channel data, selecting first signal data with a length of a preset length from a starting position R t, and performing interference suppression on the first signal data in a fractional order domain to obtain pulse pressure and channel data after interference suppression; the preset length is determined by a person skilled in the art according to the length of the transmitting signal, and in an actual application scene, the preset length is the length of the transmitting signal.
In step 205, second signal data with a length being a preset length from a start position R t is selected from the difference channel data, and interference suppression is performed on the second signal data in a fractional order domain, so as to obtain pulse pressure difference channel data after interference suppression.
In step 206, the sum and difference angle is measured according to the interference suppressed pulse pressure and the third position corresponding to the power peak in the channel data and the interference suppressed pulse pressure and the fourth position corresponding to the power peak in the channel data, so as to obtain the target angle. And selecting the numerical value of the pulse pressure and channel data at the third position after interference suppression and the numerical value of the pulse pressure difference channel data at the fourth position after interference suppression, and performing sum and difference angle measurement.
According to the method, the device and the system, the composite interference is restrained by utilizing the difference of the composite interference and the target echo in three different dimensions of a space domain, a time domain and a fractional order domain, so that the composite interference including main lobe slice interference in a pulse or among pulses can be restrained effectively, and an accurate target angle can be obtained.
In an actual application scene, if the array of the radar is a uniform linear array composed of M array elements, the signals incident to the array comprise k t target echo signals, k jm main lobe slice interference and k js side lobe noise suppression interference. The guiding vector of the linear array is; Wherein/>Representing the direction of arrival of the source, d representing the array element spacing,/>Representing wavelength,/>Representing a transpose operation.
The multi-channel array data; Wherein,，，；/>、/>And/>Directions of arrival (Direction Of Arrival, abbreviated as DOA) representing the kth t target, the kth jm main lobe slice interference, and the kth js side lobe noise suppression interference, respectively,/>、/>And/>Respectively representing a target echo signal, a main lobe slice interference signal and a side lobe noise suppression interference signal received by the array,/>Representing a complex matrix of x rows and y columns, L representing the number of snapshots and N representing noise subject to gaussian distribution.
The multi-channel array data received according to radarGenerating sum channel data and difference channel data, specifically including:
Computing the multi-channel array data Angle of spatial source/>; Wherein/>Representing conjugate transpose operation,/>Representing airspace angle; the number of the airspace angles can be selected by a person skilled in the art according to the needs, and the more the number is, the higher the angle estimation accuracy is.
According to the angle of the space information sourceJudging whether side lobe interference exists or not, if so, constructing an orthogonal projection matrix/>, according to the angle of the side lobe interference; Wherein/>The unit matrix representing M multiplied by M has the elements on opposite corner lines of all 1 and the elements at other positions of all 0; the process of calculating the spatial source may be understood as a process of converting the multichannel array data X into the spatial domain, and the determining whether the side lobe interference exists is determined according to the number of the sources in the spatial domain, as in the spatial spectrum shown in fig. 2, it can be seen that there are currently 3 angles of-4 degrees, 3.5 degrees, 15 degrees, and the like, and after the angle of the spatial source is obtained by calculation, filtering may be performed according to the angle of the spatial source to obtain the interference angle.
Using the orthogonal projection matrixFor the multichannel array data/>Performing sidelobe interference suppression to obtain data/>, after the sidelobe interference suppression。
Using the sidelobe interference suppressed dataSum channel data and difference channel data are generated.
In an actual application scenario, the data after the sidelobe interference suppression is usedGenerating sum channel data and difference channel data, specifically including:
using the sidelobe interference suppressed data Formed sum channel data is/>; Data/>, using the sidelobe interference suppressedThe resulting bad channel data is/>; Wherein/>Is the main beam direction,/>Is a vector with length of M multiplied by1, the first M/2 elements are 1, and the latter elements are-1; symbol/>Representing dot product.
In an optional embodiment, the distance gate unit where the slice interference is found from the sum channel data specifically includes:
finding a first position and a second position with abrupt change of power in the sum channel data, and judging whether the average power between the first position and the second position is higher than a preset power or not; the preset power is obtained by an empirical analysis by a person skilled in the art. The abrupt change is determined by the power difference value of the previous data and the next data, and when the power difference value of the previous data and the next data is larger than a preset difference value, the abrupt change is considered to occur; wherein the preset difference is obtained by empirical analysis by a person skilled in the art.
And if the average power between the first position and the second position is higher than the preset power, taking the area between the first position and the second position as a distance gate unit where the slice interference exists.
As shown in fig. 3, the first position is 900, the second position is 2399, and the formed range gate unit c= [900, 901, …,2399].
The constructing the self-adaptive weight for the multichannel array data in the range gate unit where the cut-off interference is located specifically comprises the following steps: selecting multichannel array data in a range gate unit where slice interference is located; Wherein,Sta is the starting position of the distance gate unit where the slice interference is located, and end is the ending position of the distance gate unit where the slice interference is located; according to/>Construction of adaptive weights/>。
The step of performing interference suppression on the first signal data in the fractional order domain to obtain pulse pressure and channel data after interference suppression specifically includes:
Calculating the average power of the first signal data, setting zero at a position higher than the average power to obtain processed first signal data, converting the processed first signal data into a fractional order domain, performing interference suppression on the processed first signal data in the fractional order domain, and performing fractional Fourier transform (Fractional Fourier Transform, abbreviated as FrFT) inverse transform to obtain pulse pressure and channel data after interference suppression. The process of processing the first signal data to obtain the processed first signal data may also be understood as time-domain concealment.
In a preferred embodiment, the data is received in accordance with radar in a multi-channel arrayBefore the sum channel data and the difference channel data are generated, band-pass filtering is further carried out on radar receiving signals so as to filter interference in a non-working bandwidth of the radar in the current period, and multichannel array data/>, is obtained。
Example 2:
the invention is based on the method described in embodiment 1, and combines specific application scenes, and the implementation process in the characteristic scene of the invention is described by means of technical expression in the relevant scene.
The embodiment provides a space-time-frequency domain combined active anti-composite interference method, which comprises the following steps:
(1) The radar transmits pulses of inter-pulse combined agility, the inter-pulse agility is realized by transmitting sub-pulses of a plurality of linear frequency modulation signals with different slopes, the inter-pulse agility is realized by transmitting different carrier frequencies, and the frequency spectrums of the transmitted signals in adjacent heavy cycles are ensured not to be overlapped.
(2) And carrying out band-pass filtering to filter interference in a non-working bandwidth in the current period of the radar.
(3) And estimating the angle of the space information source by adopting a conventional beam forming algorithm for multi-channel array data received by the radar.
(4) Judging whether side lobe interference exists, if so, performing the step (5), and if not, performing the step (6).
(5) Based on the angle estimation value of the sidelobe interference, an orthogonal projection matrix is constructed, and the sidelobe interference is restrained by using an orthogonal projection algorithm.
(6) And forming sum and difference channel data according to the main beam direction.
(7) And roughly estimating a range gate unit where the main lobe slice interference is located according to the sum channel data.
(8) And obtaining self-adaptive channel data by adopting a least mean square undistorted response (Minimum Variance Distortionless Response, MVDR) algorithm according to the interference sample distance gate unit estimated in the last step.
(9) And estimating the starting position R t of the target echo signal according to the self-adaptive channel data.
(10) From the sum and difference channel data, signal data having the same length as the transmission signal is extracted from the R t th distance gate unit, and time domain concealment processing is performed on the data.
(11) And converting the pulse data after the time domain concealment into a fractional order domain, performing interference suppression in the fractional order domain, and performing FrFT inverse conversion to obtain sum and difference channel data of the time domain.
(12) And carrying out pulse pressure processing on the sum and difference channel data, and detecting a distance gate unit where the target is located according to the sum channel data after pulse pressure.
(13) Taking the numbers in the sum and difference channel data at the peak value to realize the sum and difference angle of the target.
Assuming that the array is a uniform linear array consisting of M array elements, the signals incident to the array comprise k t target echo signals, k jm main lobe slice disturbances and k js side lobe noise suppression disturbances. The guiding vector of the linear array is as follows:
Wherein, Representing the direction of arrival (Direction Of Arrival, DOA) of the source, d representing the array element spacing,/>Representing wavelength,/>Representing a transpose operation. The array received data can be expressed as:
Wherein, ；；；/>、/>And/>DOA,/>, respectively representing the kth t target, the kth jm main lobe slice disturbance and the kth js side lobe noise suppression disturbance、/>AndRespectively representing a target echo signal, a main lobe slice interference signal and a side lobe noise suppression interference signal received by the array,/>Representing a complex matrix of x rows and y columns, L representing the number of snapshots and N representing noise subject to gaussian distribution. The time-frequency diagram of the target real echo signal is shown in fig. 4, and it can be seen that the radar transmission signal is composed of a plurality of LFM sub-pulse signals with different frequency modulation slopes.
Assuming that the signal of the previous period transmitted by the radar is forwarded by the interfering party, the start range gate of the interference can lead or lag the start range gate of the target real echo signal, the time-frequency diagram of the slice interference directly forwarded by the 1 st main lobe intermittent sampling is shown in fig. 5, and the time-frequency diagram of the slice interference repeatedly forwarded by the 2 nd main lobe intermittent sampling is shown in fig. 6. It can be seen that the first difference is that the two main lobe slice interference is different from the carrier frequency of the target real echo signal, because the radar transmit signal carrier frequency is agile between pulses; the second difference is that the frequency modulation slope change rule of the two main lobe slice interferences and the target real echo signal is different.
Assuming that the signal of the current period of radar transmission is forwarded by the interfering party, the start range gate of interference must be the start range gate of the real echo signal of the lagging target, because the slice forwarding of the interfering party needs a certain time, at this time, the time-frequency diagram of the slice interference directly forwarded by the 1 st main lobe intermittent sampling is shown in fig. 7, and the time-frequency diagram of the slice interference repeatedly forwarded by the 2 nd main lobe intermittent sampling is shown in fig. 8. It can be seen that the first difference is that the two main lobe slice interference and the target real echo signal initially appear at different times, because the interfering party needs a certain time to perceive the radar and forward the interference signal; the second difference is that the frequency modulation slope change rule of the two main lobe slice interferences and the target real echo signal is different.
Fig. 9 is a schematic flow chart of a space-time-frequency domain combined active anti-composite interference method according to the present invention, the method includes:
in step 301, the radar transmits pulses with inter-pulse joint agility, the inter-pulse agility is realized by transmitting a plurality of LFM sub-pulse signals with different slopes, the inter-pulse agility is realized by transmitting different carrier frequencies, and it is ensured that the frequency spectrums of the transmission signals in adjacent heavy cycles are not overlapped.
In step 302, bandpass filtering is performed to filter out interference in the non-operating bandwidth of the radar in the current period.
In step 303, estimating an angle of a spatial source for multi-channel array data received by the radar; the specific implementation method is that; Wherein/>Representing conjugate transpose operation,/>The number of the airspace angles can be selected according to the requirement, and the angle estimation precision is higher as the number is larger. Assuming that two main lobe slices interfere in the air and 1 side lobe suppresses interfere, a spatial spectrum as shown in fig. 2 can be obtained according to the angle calculation formula of the spatial information source, and it can be seen that the information source exists at 3 angles of-4 degrees, 3.5 degrees, 15 degrees and the like.
In step 304, it is determined whether or not there is a side lobe disturbance, if there is a side lobe disturbance, step 305 is performed, and if there is no side lobe disturbance, step 306 is performed directly.
From the source angle estimated in the last step, it is easy to know that there is a side lobe interference at 15 degrees, and the angle of the side lobe interference is recorded as。
In step 305, an orthogonal projection matrix is constructed based on the angle estimation value of the sidelobe interference, and the sidelobe interference is suppressed by using an orthogonal projection algorithm; the method for constructing the orthogonal projection matrix comprises the following steps: ; wherein/> The unit matrix representing M multiplied by M has the elements on opposite corner lines of all 1 and the elements at other positions of all 0; and then utilizing an orthogonal projection matrix to inhibit side lobe interference to obtain/>; Wherein Y represents the data after suppression of side lobe interference.
In step 306, sum and difference channel data are formed according to the main beam direction; assuming the main beam is directed toThe method for forming and channel data comprises the following steps: /(I); The method for forming the difference channel data comprises the following steps: /(I); Wherein/>Is a vector of length M x 1, the first M/2 elements being 1 and the last element being-1.
In step 307, the range gate unit where the main lobe slice interference is located is roughly estimated according to the sum channel data; and channel data results as shown in fig. 3, it can be seen that the range gate of slice disturbance occurs in the range from 900 to 2399 range gates, which is denoted as c= [900, 901, …,2399].
In step 308, adaptive channel data is obtained by using MVDR algorithm according to the interference sample distance gate unit estimated in the previous step.
The method for taking the data of the interference sample distance gate unit comprises the following steps ofFor example, when c= [900, 901, …,2399], the 900 th to 2399 th column of data X of the array reception data is selected; then get the construction adaptive weight。
The adaptive channel data is obtained by multiplying the adaptive weight w with the array received data X, and then pulse pressure is performed on the adaptive channel data, and the result is shown in fig. 10.
In step 309, the starting position of the target echo signal is estimated from the adaptive channel data; from fig. 10, it can be seen that the abscissa 1000 corresponding to the maximum peak is the target echo start range gate position, and this abscissa is labeled R t.
In step 310, from the sum and difference channel data, signal data having the same length as the transmission signal is extracted from the R t th distance gate unit, and time domain concealment is performed on the data; and taking out signal data with the same length as the transmitted signal from the R t th distance gate unit in the sum and difference channel data, averaging the signal data, and setting zero for data larger than the average value.
In step 311, the pulse data after the time-domain concealment is transformed into a fractional domain, interference suppression is performed in the fractional domain, and then FrFT inverse transformation is performed, so as to obtain sum and difference channel data of the time domain, as shown in fig. 11 and 12.
In step 312, pulse pressure processing is performed on the sum and difference channel data, and as shown in fig. 13 and 14, the distance gate unit where the target is located is detected according to the pulse pressure-processed sum channel data, and the maximum peak value can be seen at the 1000 th distance gate unit.
In step 313, the sum and difference angles of the target are obtained by taking the values in the sum and difference channel data at the peak, and fig. 15 shows the target angle estimated by 10 monte carlo simulations, so that the target estimated angle is very close to the real angle by 0.2 degrees, which indicates that the method of the embodiment has strong anti-interference capability and high target angle estimation accuracy.
The embodiment utilizes the difference of the composite interference and the target echo in three different dimensions of a space domain, a time domain and a fractional order domain to inhibit the composite interference, thereby being capable of effectively inhibiting the composite interference of main lobe slice interference in a pulse or between pulses and obtaining an accurate target angle.
Example 3:
Fig. 16 is a schematic diagram of an architecture of a space-time-frequency domain combined active anti-composite interference device according to an embodiment of the present invention. The space-time-frequency domain joint active anti-composite interference device of the present embodiment includes one or more processors 21 and a memory 22. In fig. 16, a processor 21 is taken as an example.
The processor 21 and the memory 22 may be connected by a bus or otherwise, which is illustrated in fig. 16 as a bus connection.
The memory 22 is used as a non-volatile computer readable storage medium for storing non-volatile software programs and non-volatile computer executable programs, such as the space-time-frequency-domain joint active anti-composite interference method of embodiment 1. The processor 21 executes the space-time-frequency-domain joint active anti-composite interference method by running non-volatile software programs and instructions stored in the memory 22.
The memory 22 may include high-speed random access memory, and may also include non-volatile memory, such as at least one magnetic disk storage device, flash memory device, or other non-volatile solid-state storage device. In some embodiments, memory 22 may optionally include memory located remotely from processor 21, which may be connected to processor 21 via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The program instructions/modules are stored in the memory 22, which when executed by the one or more processors 21, perform the space-time-frequency domain joint active immunity method of embodiment 1 described above.
It should be noted that, because the content of information interaction and execution process between modules and units in the above-mentioned device and system is based on the same concept as the processing method embodiment of the present invention, specific content may be referred to the description in the method embodiment of the present invention, and will not be repeated here.
Those of ordinary skill in the art will appreciate that all or a portion of the steps in the various methods of the embodiments may be implemented by a program that instructs associated hardware, the program may be stored on a computer readable storage medium, the storage medium may include: read Only Memory (ROM), random access Memory (RAM, random Access Memory), magnetic or optical disk, and the like.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather is intended to cover all modifications, equivalents, and alternatives falling within the spirit and principles of the invention.
Claims (10)
1. A space-time-frequency domain combined active anti-composite interference method is characterized in that a radar transmits inter-pulse combined agile pulses, and the method comprises the following steps:
from radar received multi-channel array data Generating sum channel data and difference channel data;
Finding a distance gate unit where the slice interference is located from the sum channel data, constructing an adaptive weight for multi-channel array data in the distance gate unit where the slice interference is located, and multiplying the adaptive weight with the multi-channel array data to obtain adaptive channel data;
pulse pressure is carried out on the self-adaptive channel data to obtain self-adaptive channel data after pulse pressure, and a position corresponding to a power peak value is found from the self-adaptive channel data after pulse pressure, namely, a starting position R t of a target echo signal;
Selecting first signal data with the length of a preset length from a starting position R t from the sum channel data, and performing interference suppression on the first signal data in a fractional order domain to obtain pulse pressure and channel data after interference suppression;
selecting second signal data with the length of a preset length from a starting position R t from the difference channel data, and performing interference suppression on the second signal data in a fractional order domain to obtain pulse pressure difference channel data after interference suppression;
And carrying out sum and difference angle measurement according to the interference-suppressed pulse pressure, a third position corresponding to the power peak value in the channel data and a fourth position corresponding to the power peak value in the interference-suppressed pulse pressure difference channel data to obtain a target angle.
2. The method for space-time-frequency domain joint active anti-composite interference according to claim 1, wherein the multi-channel array data received by the radarGenerating sum channel data and difference channel data, specifically including:
Computing the multi-channel array data Angle of spatial source/>; Wherein,Representing conjugate transpose operation,/>Representing airspace angle,/>Representing a airspace search guide vector;
According to the angle of the space information source Judging whether the side lobe interference exists or not, and if so, judging that the side lobe interference exists, and according to the angle/>, of the side lobe interferenceConstruction of an orthographic projection matrix/>; Wherein/>Representing spatial domain interference guide vector,/>The unit matrix representing M multiplied by M has the elements on opposite corner lines of all 1 and the elements at other positions of all 0;
Using the orthogonal projection matrix For the multichannel array data/>Side lobe interference suppression is carried out to obtain data/>, after side lobe interference suppression；
Data after using the side lobe interference suppressionSum channel data and difference channel data are generated.
3. The method for space-time-frequency domain joint active anti-composite interference according to claim 2, wherein the data after side lobe interference suppression is usedGenerating sum channel data and difference channel data, specifically including:
Data after using the side lobe interference suppression Formed sum channel data is/>；
Data after using the side lobe interference suppressionThe resulting bad channel data is/>; Wherein,Is the main beam direction,/>Is a vector of length M x 1, the first M/2 elements being 1 and the last M/2 elements being-1.
4. The method for space-time-frequency domain joint active anti-composite interference according to claim 1, wherein the finding the distance gate unit where the slice interference is located from the sum channel data specifically comprises:
finding a first position and a second position with abrupt change of power in the sum channel data, and judging whether the average power between the first position and the second position is higher than a preset power or not;
And if the average power between the first position and the second position is higher than the preset power, taking the area between the first position and the second position as a distance gate unit where the slice interference exists.
5. The method for space-time-frequency domain combined active anti-composite interference according to claim 1, wherein the method is characterized by constructing an adaptive weight for multichannel array data in a range gate unit where the cut-off interference is located, and specifically comprising:
Selecting multichannel array data in a range gate unit where slice interference is located ; Wherein,Sta is the starting position of the distance gate unit where the slice interference is located, and end is the ending position of the distance gate unit where the slice interference is located;
According to Construction of adaptive weights/>Wherein/>Representing the main beam pointing steering vector,/>Representing the number of shots.
6. The method of claim 1, wherein the performing interference suppression on the first signal data in the fractional order domain to obtain pulse pressure and channel data after interference suppression specifically comprises:
And calculating the average power of the first signal data, setting zero at a position higher than the average power to obtain processed first signal data, converting the processed first signal data into a fractional order domain, performing interference suppression on the processed first signal data in the fractional order domain, and performing FrFT inverse conversion to obtain pulse pressure and channel data after interference suppression.
7. The space-time-frequency domain combined active anti-composite interference method according to claim 1, wherein the radar array is a uniform linear array composed of M array elements, and the signals incident on the array comprise k t target echo signals, k jm main lobe slice interferences and k js side lobe suppression interferences; the guiding vector of the linear array is; Wherein/>Representing the direction of arrival of the source, d representing the array element spacing,/>Representing wavelength,/>Representing a transpose operation.
8. The method for space-time-frequency domain joint active anti-composite interference according to claim 1, wherein the multichannel array data; Wherein/>，，；/>、/>And/>DOA,/>, respectively representing the kth t target, the kth jm main lobe slicing disturbance and the kth js side lobe suppression disturbance、/>AndRespectively representing a target echo signal, a main lobe slice interference signal and a side lobe suppression interference signal received by the array,/>Representing a complex matrix of x rows and y columns, L representing the number of snapshots and N representing noise subject to gaussian distribution.
9. The method for space-time-frequency domain joint active anti-composite interference according to claim 1, wherein the data is received in a multi-channel array according to radarBefore the sum channel data and the difference channel data are generated, band-pass filtering is further carried out on radar receiving signals so as to filter interference in a non-working bandwidth of the radar in the current period, and multichannel array data/>, is obtained。
10. The utility model provides a space-time frequency domain unites initiative anti-composite interference device which characterized in that includes:
At least one processor; and a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor for performing the space-time-frequency domain joint active anti-composite interference method of any one of claims 1-9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410208694.1A CN117784026B (en) | 2024-02-26 | 2024-02-26 | Space-time-frequency domain combined active anti-composite interference method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410208694.1A CN117784026B (en) | 2024-02-26 | 2024-02-26 | Space-time-frequency domain combined active anti-composite interference method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117784026A CN117784026A (en) | 2024-03-29 |
| CN117784026B true CN117784026B (en) | 2024-05-03 |
Family
ID=90379974
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410208694.1A Active CN117784026B (en) | 2024-02-26 | 2024-02-26 | Space-time-frequency domain combined active anti-composite interference method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117784026B (en) |
Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5371506A (en) * | 1993-07-19 | 1994-12-06 | General Electric Co. | Simultaneous multibeam approach for cancelling multiple mainlobe jammers while preserving monopulse angle estimation accuracy on mainlobe targets |
| JP2002243824A (en) * | 2001-02-15 | 2002-08-28 | Mitsubishi Electric Corp | Angle measuring device and angle measuring method |
| JP2009162613A (en) * | 2008-01-07 | 2009-07-23 | Mitsubishi Electric Corp | Radar system |
| CN101907703A (en) * | 2010-07-20 | 2010-12-08 | 中国人民解放军空军雷达学院 | Sky-wave beyond visual range radar adaptive space-time joint interference-resistant method |
| KR20160012284A (en) * | 2014-07-23 | 2016-02-03 | 국방과학연구소 | Method and Apparatus for suppressing jammer signals and estimating Angle Of Arrival of original signal using orthogonal of transmitting signal waveform |
| KR20170021149A (en) * | 2015-08-17 | 2017-02-27 | 국방과학연구소 | Apparatus and Method for adaptive side lobe cancelation applicable to interference environment |
| CN109946664A (en) * | 2019-03-06 | 2019-06-28 | 西安电子科技大学 | A kind of array radar target seeker Monopulse estimation method under major lobe suppression |
| CN112881975A (en) * | 2021-01-08 | 2021-06-01 | 电子科技大学 | Single pulse sum and difference beam angle measurement method based on subarray characteristic matrix joint diagonalization |
| CN112949100A (en) * | 2020-11-06 | 2021-06-11 | 中国人民解放军空军工程大学 | Main lobe interference resisting method for airborne radar |
| CN113238192A (en) * | 2021-03-22 | 2021-08-10 | 西安电子科技大学 | Four-channel main lobe interference suppression angle measurement method based on feature projection preprocessing |
| CN113835068A (en) * | 2021-09-22 | 2021-12-24 | 南京信息工程大学 | Blind source separation real-time main lobe interference resisting method based on independent component analysis |
| CN114935738A (en) * | 2022-05-06 | 2022-08-23 | 中国人民解放军空军预警学院 | Radar sum and difference angle measurement method based on sample identification under time domain discrete main lobe interference |
| CN115932750A (en) * | 2022-12-23 | 2023-04-07 | 电子科技大学 | Radar space-time-frequency-energy multi-domain joint intelligent active anti-interference method |
| CN116068502A (en) * | 2023-04-06 | 2023-05-05 | 中国人民解放军空军预警学院 | Multi-domain combined anti-composite interference method, device and system |
| CN116068498A (en) * | 2023-02-02 | 2023-05-05 | 中国人民解放军63892部队 | Multi-interference-source suppression method based on main and auxiliary lobe airspace polarization characteristics |
| CN116125420A (en) * | 2023-01-12 | 2023-05-16 | 电子科技大学 | Polarization-space-time domain combined main lobe suppression interference resisting method based on tensor decomposition |
| CN116774162A (en) * | 2023-05-15 | 2023-09-19 | 西安电子科技大学 | Main and side lobe interference resistant self-adaptive monopulse angle measurement method for uniform linear array |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10140498C1 (en) * | 2001-08-17 | 2003-05-15 | Eads Deutschland Gmbh | Method for the suppression of jammer signals |
| US9473182B2 (en) * | 2014-11-07 | 2016-10-18 | Nokia Solutions And Networks Oy | Oportunistic interference cancellation of radar signals |
-
2024
- 2024-02-26 CN CN202410208694.1A patent/CN117784026B/en active Active
Patent Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5371506A (en) * | 1993-07-19 | 1994-12-06 | General Electric Co. | Simultaneous multibeam approach for cancelling multiple mainlobe jammers while preserving monopulse angle estimation accuracy on mainlobe targets |
| JP2002243824A (en) * | 2001-02-15 | 2002-08-28 | Mitsubishi Electric Corp | Angle measuring device and angle measuring method |
| JP2009162613A (en) * | 2008-01-07 | 2009-07-23 | Mitsubishi Electric Corp | Radar system |
| CN101907703A (en) * | 2010-07-20 | 2010-12-08 | 中国人民解放军空军雷达学院 | Sky-wave beyond visual range radar adaptive space-time joint interference-resistant method |
| KR20160012284A (en) * | 2014-07-23 | 2016-02-03 | 국방과학연구소 | Method and Apparatus for suppressing jammer signals and estimating Angle Of Arrival of original signal using orthogonal of transmitting signal waveform |
| KR20170021149A (en) * | 2015-08-17 | 2017-02-27 | 국방과학연구소 | Apparatus and Method for adaptive side lobe cancelation applicable to interference environment |
| CN109946664A (en) * | 2019-03-06 | 2019-06-28 | 西安电子科技大学 | A kind of array radar target seeker Monopulse estimation method under major lobe suppression |
| CN112949100A (en) * | 2020-11-06 | 2021-06-11 | 中国人民解放军空军工程大学 | Main lobe interference resisting method for airborne radar |
| CN112881975A (en) * | 2021-01-08 | 2021-06-01 | 电子科技大学 | Single pulse sum and difference beam angle measurement method based on subarray characteristic matrix joint diagonalization |
| CN113238192A (en) * | 2021-03-22 | 2021-08-10 | 西安电子科技大学 | Four-channel main lobe interference suppression angle measurement method based on feature projection preprocessing |
| CN113835068A (en) * | 2021-09-22 | 2021-12-24 | 南京信息工程大学 | Blind source separation real-time main lobe interference resisting method based on independent component analysis |
| CN114935738A (en) * | 2022-05-06 | 2022-08-23 | 中国人民解放军空军预警学院 | Radar sum and difference angle measurement method based on sample identification under time domain discrete main lobe interference |
| CN115932750A (en) * | 2022-12-23 | 2023-04-07 | 电子科技大学 | Radar space-time-frequency-energy multi-domain joint intelligent active anti-interference method |
| CN116125420A (en) * | 2023-01-12 | 2023-05-16 | 电子科技大学 | Polarization-space-time domain combined main lobe suppression interference resisting method based on tensor decomposition |
| CN116068498A (en) * | 2023-02-02 | 2023-05-05 | 中国人民解放军63892部队 | Multi-interference-source suppression method based on main and auxiliary lobe airspace polarization characteristics |
| CN116068502A (en) * | 2023-04-06 | 2023-05-05 | 中国人民解放军空军预警学院 | Multi-domain combined anti-composite interference method, device and system |
| CN116774162A (en) * | 2023-05-15 | 2023-09-19 | 西安电子科技大学 | Main and side lobe interference resistant self-adaptive monopulse angle measurement method for uniform linear array |
Non-Patent Citations (2)
| Title |
|---|
| 一种基于干扰重构的自适应干扰对消算法;周伟光;谭怀英;李归;;现代雷达;20180615(06);全文 * |
| 基于非均匀子阵的双和/三差通道同时抑制主副瓣干扰;孙晨伟;陶海红;郭晓双;宋嘉奇;;电子学报;20190415(04);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117784026A (en) | 2024-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | A method for finding best channels in beam-space post-Doppler reduced-dimension STAP | |
| Sira et al. | Waveform-agile sensing for tracking | |
| Patton et al. | Efficient design of radar waveforms for optimised detection in coloured noise | |
| US20230152424A1 (en) | Method for eliminating one-bit signal harmonic false target and related component | |
| Hamid et al. | Performance comparison of time-domain and frequency-domain beamforming techniques for sensor array processing | |
| CN105510887B (en) | To the active suppressing formula disturbance restraining method of airborne radar under a kind of clutter background | |
| CN107153178B (en) | Target detection method for external radiation source radar reference signal with multipath interference | |
| CN104459667B (en) | Sparse array DOA estimation method based on CLEAN | |
| CN104502904B (en) | Torpedo homing beam sharpening method | |
| CN103728595A (en) | Subspace-projection-based networked radar inhibition pressing type main lobe interference method | |
| CN104678372A (en) | Joint estimation method for super-resolution distance value and angle value by using orthogonal frequency division multiplexing radar | |
| Li et al. | Dimension-reduced space-time adaptive clutter suppression algorithm based on lower-rank approximation to weight matrix in airborne radar | |
| CN110531311A (en) | A kind of LTE external illuminators-based radar DOA estimation method based on matrix recombination | |
| Li et al. | Parameter estimation based on fractional power spectrum density in bistatic MIMO radar system under impulsive noise environment | |
| CN109655819A (en) | A kind of clutter recognition three-D imaging method based on real aperture Doppler beam sharpening | |
| Yang et al. | Enhanced knowledge-aided space–time adaptive processing exploiting inaccurate prior knowledge of the array manifold | |
| CN108490425B (en) | Angle measuring method of bistatic MIMO radar | |
| CN115508799A (en) | Distributed passive radar target detection method based on moment space | |
| Smith et al. | A new approach to moving target detection using unit circle roots constrained adaptive matched filter | |
| CN117784026B (en) | Space-time-frequency domain combined active anti-composite interference method and device | |
| CN102932034B (en) | Fast broadband coherent source direction estimation method | |
| CN116068502B (en) | Multi-domain combined anti-composite interference method, device and system | |
| CN105319538B (en) | Space-time adaptive clutter suppression method based on accessory channel | |
| CN109490859B (en) | Robust detector for orientation and Doppler disturbances in a partially homogeneous environment | |
| CN112014807A (en) | Self-adaptive clutter suppression method for frequency agile 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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |