CN108957419A - Asynchronous interference suppressing method based on notch filter processing - Google Patents
Asynchronous interference suppressing method based on notch filter processing Download PDFInfo
- Publication number
- CN108957419A CN108957419A CN201810548240.3A CN201810548240A CN108957419A CN 108957419 A CN108957419 A CN 108957419A CN 201810548240 A CN201810548240 A CN 201810548240A CN 108957419 A CN108957419 A CN 108957419A
- Authority
- CN
- China
- Prior art keywords
- radar
- doppler
- airborne radar
- data matrix
- vector
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/38—Jamming means, e.g. producing false echoes
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
The invention discloses a kind of asynchronous interference suppressing methods based on notch filtering processing, its main thought are as follows: determine airborne radar, there are targets in the airborne radar detection range, the echo-signal that airborne radar emits signal into its detection range and receives after target reflects is denoted as radar raw radar data matrix;According to the radar raw radar data matrix, distance-Doppler numeric field data matrix is obtained;Then main-lobe clutter is determined;The Doppler frequency of main-lobe clutter is calculated, and obtains the column vector g that length is R;Wherein R is the positive integer greater than 1;It determines and refers to thresholdingSuccessively obtain modification treated result gproWith the power and vector G after normalization, and then notch filter weight coefficient vector F is obtained;Slide window processing is carried out to the radar raw radar data matrix using notch filter weight coefficient vector F, it is after obtaining slide window processing as a result, the result after the slide window processing is the asynchronous interference suppression result being filtered based on notch.
Description
Technical field
The invention belongs to Radar Technology field, in particular to a kind of asynchronous interference inhibition side based on notch filtering processing
Method is suitable for airborne radar echo asynchronous interference and inhibits.
Background technique
Airborne radar is considered as the strategic force for capableing of left and right situation of battlefield by the military of various countries with its unique characteristics of operation
Device.Can AF panel performance be the principal element that influences airborne radar and normally detect, therefore, airborne radar interference mitigation technology
Receive the attention of various countries researcher.
In radar signal environment, interfere being constantly present;Common interference can be mainly divided into Deceiving interference, obstruction
Formula interference, point frequency (nicking) interference and four kinds of asynchronous interference;Wherein asynchronous interference (Asynchronous Interference)
Mainly from the electromagnetic radiation of the equipment of industrial product, communication equipment, other radars etc., the characteristics of asynchronous interference is with burst pulse shape
Formula occurs at random, and amplitude is much larger than noise floor, therefore, also known as singular value (Singular Value);From origin
Say itself there should be certain periodicity, but the work of its frequency and radar that occur is out of step;Therefore, asynchronous dry
Disturb the time occurred in radar receiver be it is unfixed, show very big randomness.
On the other hand, the amplitude of asynchronous interference is very big, is much larger than signal and noise level, can even reach clutter sometimes
It is horizontal;Asynchronous interference shows as having one fixed width, more in distance domain in pulse Doppler (PD, Pulse Doppler) figure
The general horizontal stripe strangling domain and being full of;It, would generally since its width is similar to target echo signal in the signal detection of radar
It is detected as target;In automatic detection, it will usually Plot coherence technology is used, to reduce the quantity of original point mark;And due to
The intensity of asynchronous interference is too strong, can be due to Plot coherence algorithm, so that asynchronous interference will if there is target appearance in its vicinity
The target in its contiguous range is appeared in covering, causes target that indescribable loss can occur at certain moment or go out
Existing mistake (be embodied in detections of radar and go out singular value), to influence radar to the detection performance of target;In addition, radar is often adopted
The probability of detection that target is improved with accumulative means can make making an uproar for each channel after accumulating due to the presence of asynchronous interference signal
Sound substrate significantly improves, this will reduce the detection probability of target again.
Asynchronous interference is very high in the probability that certain frequency ranges (especially metric wave), certain occasions occur, and intensity is very big, and another
Some environment then influence smaller outside, or even do not occur;Therefore, it is necessary to adaptive method be designed in Radar Signal Processing to locate
Reason asynchronous interference carries out inhibition elimination to it when asynchronous interference occurs, and in no asynchronous interference, operates without inhibiting, with
Reduce signal processing loss.
Summary of the invention
In view of the above-mentioned problems of the prior art, it is an object of the invention to propose it is a kind of based on notch filtering processing
Asynchronous interference suppressing method, this kind based on the asynchronous interference suppressing method that notch is filtered can be adaptive inhibition it is asynchronous dry
It disturbs, when asynchronous interference occurs, is adaptively calculated out corresponding filtering weight coefficient and it is eliminated, in no asynchronous interference
When, it is operated without inhibiting, to reduce signal processing loss.
Realize the object of the invention main thought: asynchronous dry using radar raw data matrix before carrying out pulse compression
The characteristics of showing as a filament of single or several range gates after progress adding window Fourier transformation is tieed up in pulse is disturbed, is calculated adaptive
Notch filter weight coefficient is answered to inhibit it.
To reach above-mentioned technical purpose, the present invention is realised by adopting the following technical scheme.
A kind of asynchronous interference suppressing method based on notch filtering processing, comprising the following steps:
Step 1, airborne radar is determined, there are target in the airborne radar detection range, airborne radar detects model to it
The echo-signal enclosing interior transmitting signal and receiving after target reflects, is denoted as radar raw radar data matrix;According to described
Radar raw radar data matrix obtains distance-Doppler numeric field data matrix;Then main-lobe clutter is determined;
Step 2, the Doppler frequency of main-lobe clutter is calculated, and according to the distance-Doppler numeric field data matrix, is grown
Degree is the column vector g of R;Wherein R is the positive integer greater than 1;
Step 3, the column vector g for being R according to length is determined and is referred to thresholding
Step 4, according to the column vector for being R with reference to thresholding and the length, modification is obtained treated result gpro;
Step 5, according to modification treated result gproWith reference thresholdingPower and vector G after being normalized;
Step 6, according to the power and vector G after normalization, notch filter weight coefficient vector F is obtained;
Step 7, slide window processing is carried out to the radar raw radar data matrix using notch filter weight coefficient vector F,
It is after obtaining slide window processing as a result, result after the slide window processing is that asynchronous interference be filtered based on notch inhibits to tie
Fruit.
Compared with the prior art, the present invention has the following advantages:
First, the present invention can adaptively inhibit the asynchronous interference of different distance and frequency, or else there is asynchronous interference
Data when being handled adaptive-filtering weight coefficient be all 1, to ensure that treatment loss will not be caused to signal data.
Second, existing asynchronous interference suppressing method is offseted using two pulse delays, offsets result modulus and the frame that is delayed;
FAR processing and foundation and update apart from clutter map;Threshold value was detected, its distance unit position is recorded, and inserts one and moves into one's husband's household upon marriage
Limit table;Thresholding table was examined or check, determines singular value position;The method that signal at asynchronous interference is replaced with adjacent signal interpolation,
Calculating process is complicated, and time-consuming;Method of the invention directly carries out initial data to draw window processing, calculation amount by adaptively weighing
Few, process is simple, time-consuming shorter.
Detailed description of the invention
Invention is further described in detail with reference to the accompanying drawings and detailed description.
Fig. 1 is a kind of asynchronous interference suppressing method flow chart based on notch filtering processing of the invention;
Fig. 2 is the range Doppler spectrogram of radar raw radar data matrix;
Fig. 3 is that radar raw radar data matrix ties up the result signal obtain after adding window Fourier transformation along pulse
Figure;
Fig. 4 is the adaptive weight figure obtained after calculating;
Fig. 5 is the result signal for carrying out obtaining after pulse compression and range Doppler processing again after the present invention inhibits
Figure.
Specific embodiment
It referring to Fig.1, is a kind of asynchronous interference suppressing method flow chart based on notch filtering processing of the invention;Wherein institute
State the asynchronous interference suppressing method based on notch filtering processing, comprising the following steps:
Step 1, airborne radar is determined, there are target in the airborne radar detection range, airborne radar detects model to it
The echo-signal enclosing interior transmitting signal and receiving after target reflects, is denoted as radar raw radar data matrix BN×M'×R, thunder
Up to raw radar data matrix BN×M'×RFor array number N multiplied by umber of pulse M' multiplied by the data block of range gate number R, and radar is original
Echo data matrix BN×M'×RIt is N × M' × R three-dimensional matrice, N indicates the array element total number that airborne radar includes, M' table
Show that the pulse of the interior airborne radar transmitting of coherent processing inteval (CPI, Coherent Processing Interval) is always a
Number, R indicate airborne radar its detection range is divided after include range gate total number,PRF indicates pulse
Repetition rate, B indicate that airborne radar emits the bandwidth of signal into its detection range.
Directly to radar raw radar data matrix B without process of pulse-compressionN×M'×RIt is tieed up along pulse
It carries out plus the Fast Fourier Transform (FFT) (FFT) of Chebyshev window converts, obtain Fast Fourier Transform (FFT) (FFT) transformed N × M
The three-dimensional matrice B of × Rceil;N × M × R three-dimensional matrice B after the Fast Fourier Transform (FFT) (FFT)ceilIn each array element it is corresponding
The two-dimensional matrix that one size is R × M, that is, have:
Bceil=[Bceil(1) Bceil(2) … Bceil(i) … Bceil(N)]
Wherein, i=1,2 ..., N, Bceil(i) three-dimensional of Fast Fourier Transform (FFT) (FFT) transformed N × M × R is indicated
Matrix BceilIn the corresponding size of i-th of array element be R × M two-dimensional matrix, and Fast Fourier Transform (FFT) (FFT) transformed N
The three-dimensional matrice B of × M × RceilIn the corresponding size of i-th of array element be R × M two-dimensional matrix BceilIt (i) include R × M number
According to expression formula are as follows:
Wherein, b11(i) the three-dimensional matrice B of Fast Fourier Transform (FFT) (FFT) transformed N × M × R is indicatedceilIn i-th
The corresponding size of array element is the two-dimensional matrix B of R × Mceil(i) the 1st range gate in, the data at the 1st Doppler channel, b1M
(i) the three-dimensional matrice B of Fast Fourier Transform (FFT) (FFT) transformed N × M × R is indicatedceilIn the corresponding size of i-th of array element
For the two-dimensional matrix B of R × Mceil(i) the 1st range gate in, the data at m-th Doppler channel, bR1(i) quick Fu is indicated
In leaf transformation (FFT) transformed N × M × R three-dimensional matrice BceilIn the corresponding size of i-th of array element be R × M Two-Dimensional Moment
Battle array Bceil(i) the R range gate in, the data at the 1st Doppler channel, bRM(i) indicate that Fast Fourier Transform (FFT) (FFT) becomes
N × M × R three-dimensional matrice B after changingceilIn the corresponding size of i-th of array element be R × M two-dimensional matrix Bceil(i) R in
Data at range gate, m-th Doppler channel.
By N × M × R three-dimensional matrice B after the Fast Fourier Transform (FFT) (FFT)ceilIn it is each array element corresponding one big
The small two-dimensional matrix for R × M adds up, and then obtains transformed block of frequency domain data, is denoted as distance-Doppler numeric field data square
Battle array BFFT, calculation expression are as follows:
Distance-Doppler numeric field data matrix B at this timeFFTIt is the two-dimensional matrix of a R × M, and distance-Doppler numeric field data
Matrix BFFTIn include R × M distance-Doppler numeric field data, M expression distance-Doppler numeric field data matrix BFFTIncluding it is how general
Strangle channel total number, and with a coherent processing inteval (CPI, Coherent Processing Interval) interior airborne radar
The pulse total number M' value of transmitting is equal;R indicate airborne radar its detection range is divided after include range gate it is total
Number,PRF indicates pulse recurrence frequency, and B indicates that airborne radar emits the bandwidth of signal into its detection range.
Greatest irradiation beam definitions airborne radar emitted into its detection range in signal are main lobe wave beam, by main lobe
The direction of illumination of wave beam be defined as main beam direction, by main lobe beam to ground and through ground return generation echo-signal
It is defined as main-lobe clutter;Signal is transformed into frequency domain, main-lobe clutter meeting by time domain due to having carried out Fast Fourier Transform (FFT)
Assembled by compression, be in some Doppler frequency in the form of expression in distance-Doppler figure (reaction is in Doppler channel)
It is upper to form the vertical line with one fixed width, as in Fig. 2 at No. 30 Doppler channel.
Step 2, the Doppler frequency f of main-lobe clutter is calculatedd:
Wherein, v is the flying speed of airborne radar carrier aircraft, and λ is the wave that airborne radar emits signal into its detection range
It is long, φ0For the angle of main beam position and airborne radar carrier aircraft flying speed direction, closed according to airborne radar carrier aircraft flight geometry
System can obtain:
Wherein,For the azimuth of main beam position, θ0For the pitch angle of main beam position, sin indicates SIN function,
Cos indicates cosine function.
Obtain the Doppler frequency f of main-lobe clutterdLater, it needs in distance-Doppler numeric field data matrix BFFTMiddle removal is more
General Le frequency is in fdNeighbouring clutter data obtains the data of circle of good definition part to eliminate the influence of main-lobe clutter energy;Here
Along Doppler-pulse frequency domain data matrix BFFTDoppler direction find out the Doppler frequency f of main-lobe clutterd, and choose with master
The Doppler frequency f of valve clutterdCentered on, length beWidth is the region of R, is denoted as
Two-dimensional matrix, will be describedTwo-dimensional matrix inA distance-Doppler numeric field data is all rejected, and will
Distance-Doppler numeric field data matrix BFFTDescribed inThe remaining area Liang Ge after a distance-Doppler numeric field data whole rejecting
Domain is successively spliced, i.e., willTwo-dimensional matrix andTwo-dimensional matrix sequential concatenation, obtain
It arrivesTwo-dimensional matrix B, will be describedTwo-dimensional matrix B be denoted as clear area, and it is describedTwo
Tieing up each distance-Doppler numeric field data in matrix B is clear area data.
Again to describedTwo-dimensional matrix B in every a lineA Doppler-pulse frequency domain data is added, by phase
Result after adding be denoted as a range gate power and, and then obtain the power of R range gate with and by the function of R range gate
Rate and be denoted as length be R column vector g;Wherein, the range gate for including after R expression airborne radar divides its detection range
Total number.
Step 3, according to analysis before it is understood that asynchronous interference in radar signal amplitude be much larger than noise floor, because
This is embodied in the singular point one by one in clutter data, it is therefore desirable to remove singular point, specific practice is as follows:
To the power of R range gate in the column vector g that length is R with from small to large sorted, will sort from small to large
Obtained result afterwards is denoted as column vector g after the sequence that length is Rsort, and remove column vector g after the sequence that length is RsortIn
The excessive point of power, since length is column vector g after the sequence of RsortIt is the power of the R range gate after sorting from small to large
With, therefore only need the column vector g after length is the sequence of RsortThe data of middle selection suitable position.
Specific practice are as follows: by length be R sequence after column vector gsortIn the 1st range gate power and toIt is a
The power of range gate and Yi JiThe power of a range gate and to the R range gate power and all reject, will be remainingThe power of a range gate be successively added, and will add up rear result divided byAnd then it is flat to obtain statistics
Mean value, using the assembly average as with reference to thresholdingIts calculation expression are as follows:
Wherein,gsort(i') indicate length for column vector g after the sequence of RsortIn i-th ' a range gate
Power and, R indicate airborne radar its detection range is divided after include range gate total number,PRF table
Show pulse recurrence frequency, B indicates that airborne radar emits the bandwidth of signal into its detection range.
Remove length herein as column vector g after the sequence of RsortIn the 1st range gate power and toA range gate
Power and be to prevent too small sample from interfering overall data sample to guarantee accuracy.
Step 4, it takes with reference to thresholdingCertain multiple k after obtainAs carry out later detection judgement threshold value, this
In k be setting scale parameter, 1 < k < 10, k value is 4 in the present embodiment;It can be modified processing according to different situations,
Length is less than to be all in the column vector g of RRange gate power and replace all withThe column vector g for being R by length
In all be greater than or equal toRange gate power and remain unchanged;And then obtain modification treated result gpro, at this time
Modification treated result gproIt is the column vector that a length is R;Wherein, R indicates that airborne radar draws its detection range
The range gate total number for including after point.
Step 5, to modification treated result gproAbout the reference thresholding acquired in step 3It is normalized,
Power and vector G after being normalized, calculation expression are as follows:
Power and vector G after normalizing at this time are the column vectors that a length is R, and R indicates that airborne radar detects it
Range divided after include range gate total number.
Step 6, to after the normalization acquired in step 5 power and vector G asked and operated, will ask down operation after
To result be denoted as notch filter weight coefficient vector F, calculation expression are as follows:
F=1/G
At this time notch filter weight coefficient vector F be length be R column vector, R indicate airborne radar to its detection range into
The range gate total number for including after row division.
Step 7, using sliding window algorithm with the notch filter weight coefficient vector F that is acquired in step 6 to radar original echo number
According to matrix BN×M'×RAlong pulse tie up carry out slide window processing, it is after obtaining slide window processing as a result, the result after the slide window processing i.e.
For the asynchronous interference suppression result being filtered based on notch;Wherein, M' indicates coherent processing inteval (CPI, a Coherent
Processing Interval) transmitting of interior airborne radar pulse total number.
Advantages of the present invention can be further illustrated by following l-G simulation test.
(1) experiment parameter and experiment condition
The parameter that this experiment uses is as follows:
1) for airborne radar antenna using 2 rows × 16 column planar array, array element spacing is airborne radar transmitted waveform half-wave
Long, then available size is N × M × R radar return data after pitching filtering;Radar front tilt filter battle array is placed.
2) emit 101 coherent accumulation pulses in the same coherent processing inteval CPI, pulse recurrence frequency is
2.203kHz;Distance samples frequency is 2MHz;It is 176 ° that main beam, which is directed toward with carrier aircraft head angle, and yaw angle is 5 °;Carrier aircraft is high
Degree is 8.3 kilometers, and level flies at a constant speed, speed 149m/s;Earth radius is 6378 kilometers.
(2) experiment content and interpretation of result
A. this experiment first carries out at normal pulse compression and pulse Doppler radar raw radar data matrix
Reason, processing result are as shown in Figure 2;Wherein, abscissa indicates that Doppler's port number of signal, ordinate indicate the range gate of signal
Number, figure it is seen that having a large amount of apparent horizontal stripes at No. 50-150 and 670-770 range gate, while 300
There is a small amount of weaker horizontal stripe at number range gate, these horizontal stripes are asynchronous interference.
B. radar return data process according to the invention is handled;Fig. 3 is radar raw radar data matrix edge
Pulse dimension carries out the result schematic diagram obtained after adding window Fourier transformation, and comparison diagram 2 can see the energy of asynchronous interference at this time
50,70,90,270 and No. 680 range gates are concentrated on, the thin horizontal line that a plurality of energy is concentrated is formed;Fig. 4 is to obtain after calculating
Adaptive weight figure, it can be seen that except interference concentrated area foring adaptive trap recess to inhibit to interfere, and its
Remaining part fraction value is 1, will not change original radar return data.
C. Fig. 5 is the result signal for carrying out obtaining after pulse compression and range Doppler processing again after the present invention inhibits
Figure, comparison diagram 2 illustrate that asynchronous interference has obtained this it appears that the horizontal stripe of Fig. 2 corresponding part has disappeared in Fig. 5
It effectively inhibits;From the point of view of the result of Fig. 5, method of the invention can effectively inhibit asynchronous interference, and inhibitory effect is very good.
In conclusion emulation experiment demonstrates correctness of the invention, validity and reliability.
Obviously, various changes and modifications can be made to the invention without departing from essence of the invention by those skilled in the art
Mind and range;In this way, if these modifications and changes of the present invention belongs to the range of the claims in the present invention and its equivalent technologies
Within, then the present invention is also intended to include these modifications and variations.
Claims (9)
1. a kind of asynchronous interference suppressing method based on notch filtering processing, which comprises the following steps:
Step 1, airborne radar is determined, there are targets in the airborne radar detection range, and airborne radar is into its detection range
Transmitting signal and the echo-signal received after target reflects, are denoted as radar raw radar data matrix;According to the radar
Raw radar data matrix obtains distance-Doppler numeric field data matrix;Then main-lobe clutter is determined;
Step 2, the Doppler frequency of main-lobe clutter is calculated, and according to the distance-Doppler numeric field data matrix, obtaining length is
The column vector g of R;Wherein R is the positive integer greater than 1;
Step 3, the column vector g for being R according to length is determined and is referred to thresholding
Step 4, according to the column vector for being R with reference to thresholding and the length, modification is obtained treated result gpro;
Step 5, according to modification treated result gproWith reference thresholdingPower and vector G after being normalized;
Step 6, according to the power and vector G after normalization, notch filter weight coefficient vector F is obtained;
Step 7, slide window processing is carried out to the radar raw radar data matrix using notch filter weight coefficient vector F, obtained
It is after slide window processing as a result, result after the slide window processing is the asynchronous interference suppression result being filtered based on notch.
2. a kind of asynchronous interference suppressing method based on notch filtering processing as described in claim 1, which is characterized in that in step
In rapid 1, the radar raw radar data matrix is BN×M'×R, and radar raw radar data matrix BN×M'×RIt is N × M' × R
Three-dimensional matrice, N indicates the airborne radar array element total number that includes, and M' indicates airborne radar transmitting in a coherent processing inteval
Pulse total number, R indicate airborne radar its detection range is divided after include range gate total number,
PRF indicates pulse recurrence frequency, and B indicates that airborne radar emits the bandwidth of signal into its detection range;
The distance-Doppler numeric field data matrix, obtains process are as follows:
To radar raw radar data matrix BN×M'×RIt is carried out along pulse dimension plus the Fast Fourier Transform (FFT) of Chebyshev window converts,
Obtain the transformed N × M of Fast Fourier Transform (FFT) × R three-dimensional matrice Bceil, expression formula are as follows:
Bceil=[Bceil(1) Bceil(2) … Bceil(i) … Bceil(N)]
Wherein, i=1,2 ..., N, Bceil(i) N × M × R three-dimensional matrice B after Fast Fourier Transform (FFT) is indicatedceilIn i-th
The corresponding size of array element is the two-dimensional matrix of R × M;M indicates distance-Doppler numeric field data matrix BFFTIncluding Doppler channel
Total number, and it is equal with the pulse total number M' value that airborne radar in a coherent processing inteval emits;
By N × M after the Fast Fourier Transform (FFT) × R three-dimensional matrice BceilIn the corresponding size of each array element be R × M's
Two-dimensional matrix adds up, and then obtains transformed block of frequency domain data, is denoted as distance-Doppler numeric field data matrix BFFT, meter
Operator expression formula are as follows:
The distance-Doppler numeric field data matrix BFFTIt is the two-dimensional matrix of a R × M, and distance-Doppler numeric field data matrix
BFFTIn include R × M Doppler-pulse frequency domain data.
3. a kind of asynchronous interference suppressing method based on notch filtering processing as described in claim 1, which is characterized in that in step
In rapid 1, the main-lobe clutter, determination process are as follows:
Greatest irradiation beam definitions airborne radar emitted into its detection range in signal are main lobe wave beam, by main lobe wave beam
Direction of illumination be defined as main beam direction, by main lobe beam to ground and through ground return generate echo-signal definition
For main-lobe clutter.
4. a kind of asynchronous interference suppressing method based on notch filtering processing as claimed in claim 2 or claim 3, which is characterized in that
In step 2, the Doppler frequency of the main-lobe clutter is fd, calculation expression are as follows:
Wherein, v is the flying speed of airborne radar carrier aircraft, and λ is the wavelength that airborne radar emits signal into its detection range, φ0
For the angle of main beam position and airborne radar carrier aircraft flying speed direction;
The length is the column vector of R, obtains process are as follows:
Along Doppler-pulse frequency domain data matrix BFFTDoppler direction find out the Doppler frequency f of main-lobe clutterd, and choose
With the Doppler frequency f of main-lobe clutterdCentered on, length beWidth is the region of R, is denoted asTwo-dimensional matrix, will be describedTwo-dimensional matrix inA distance-Doppler numeric field data is all picked
It removes, and by distance-Doppler numeric field data matrix BFFTDescribed inIt is surplus after a distance-Doppler numeric field data whole rejecting
Remaininging two regions successively splices, i.e., willTwo-dimensional matrix andTwo-dimensional matrix it is suitable
Sequence splicing, obtainsTwo-dimensional matrix B;
Again to describedTwo-dimensional matrix B in every a lineA Doppler-pulse frequency domain data is added, after will add up
As a result be denoted as a range gate power and, and then obtain the power of R range gate with and by the power and note of R range gate
The column vector g for being R for length;Wherein, the range gate for including after R expression airborne radar divides its detection range is always a
Number.
5. a kind of asynchronous interference suppressing method based on notch filtering processing as claimed in claim 4, which is characterized in that in step
It is described to refer to thresholding in rapid 3Its determination process are as follows:
It is column vector g after the sequence of R by lengthsortIn the 1st range gate power and toThe power of a range gate and, with
And theThe power of a range gate and to the R range gate power and all reject, will be remainingA range gate
Power be successively added, and will add up rear result divided byAnd then assembly average is obtained, the statistics is flat
Mean value, which is used as, refers to thresholdingIts calculation expression are as follows:
Wherein,gsort(i') indicate length for column vector g after the sequence of RsortIn i-th ' a range gate function
Rate and, R indicate airborne radar its detection range is divided after include range gate total number,PRF indicates arteries and veins
Repetition rate is rushed, B indicates that airborne radar emits the bandwidth of signal into its detection range.
6. a kind of asynchronous interference suppressing method based on notch filtering processing as described in claim 1, which is characterized in that in step
In rapid 4, modification treated the result gpro, obtain process are as follows:
Length is less than to be all in the column vector g of RRange gate power and replace all withThe column for being R by length
It is all in vector g to be greater than or equal toRange gate power and remain unchanged;And then obtain modification treated result gpro,
Modification treated the result gproIt is the column vector that a length is R;Wherein, k indicates the scale parameter of setting, 1 < k < 10;
R indicate airborne radar its detection range is divided after include range gate total number,PRF indicates that pulse repeats
Frequency, B indicate that airborne radar emits the bandwidth of signal into its detection range.
7. a kind of asynchronous interference suppressing method based on notch filtering processing as described in claim 1, which is characterized in that in step
In rapid 5, power and vector G after the normalization obtain process are as follows:
To modification treated result gproAbout reference thresholdingIt is normalized, power and vector after being normalized
G, calculation expression are as follows:
Power and vector G after the normalization are the column vectors that a length is R, and R indicates airborne radar to its detection range
The range gate total number for including after being divided,PRF indicates pulse recurrence frequency, and B indicates that airborne radar is examined to it
Survey the bandwidth of transmitting signal in range.
8. a kind of asynchronous interference suppressing method based on notch filtering processing as described in claim 1, which is characterized in that in step
In rapid 6, the notch filter weight coefficient vector F, calculation expression are as follows:
F=1/G
The notch filter weight coefficient vector F is the column vector that length is R, and R indicates that airborne radar draws its detection range
The range gate total number for including after point,PRF indicates pulse recurrence frequency, and B indicates airborne radar to its detection range
The bandwidth of interior transmitting signal.
9. a kind of asynchronous interference suppressing method based on notch filtering processing as claimed in claim 6, which is characterized in that in step
It is after the slide window processing as a result, being specifically with the notch filter weight coefficient vector F using sliding window algorithm to radar in rapid 7
Raw radar data matrix BN×M'×RThe result for carrying out obtaining after slide window processing is tieed up along pulse;Wherein, M' indicates a coherent
Reason is spaced the pulse total number of interior airborne radar transmitting.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810548240.3A CN108957419B (en) | 2018-05-31 | 2018-05-31 | Asynchronous interference suppression method based on notch filtering processing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810548240.3A CN108957419B (en) | 2018-05-31 | 2018-05-31 | Asynchronous interference suppression method based on notch filtering processing |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108957419A true CN108957419A (en) | 2018-12-07 |
CN108957419B CN108957419B (en) | 2022-08-19 |
Family
ID=64492848
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810548240.3A Active CN108957419B (en) | 2018-05-31 | 2018-05-31 | Asynchronous interference suppression method based on notch filtering processing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108957419B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110618411A (en) * | 2019-10-23 | 2019-12-27 | 电子科技大学 | Airborne radar clutter real-time signal generation method |
CN111007466A (en) * | 2019-12-03 | 2020-04-14 | 西安电子科技大学 | DCAR interference suppression method and system based on introduced range gate degree of freedom |
CN112014806A (en) * | 2020-08-14 | 2020-12-01 | 西安电子科技大学 | Method for suppressing unintentional interference of airborne radar in complex interference scene |
RU2751532C1 (en) * | 2021-04-06 | 2021-07-14 | Акционерное общество "Всероссийский научно-исследовательский институт радиотехники" (АО "ВНИИРТ") | Inter-review device for mapping asynchronous impulse interference for pulse-doppler radar stations and method for its implementation |
CN114690143A (en) * | 2022-04-19 | 2022-07-01 | 珠海微度芯创科技有限责任公司 | Method and device for suppressing radar clutter, radar and medium |
CN116224249A (en) * | 2023-03-29 | 2023-06-06 | 上海雷骥电子科技有限公司 | Doppler frequency width acquisition method for main lobe clutter region of airborne radar |
CN112014806B (en) * | 2020-08-14 | 2024-04-23 | 西安电子科技大学 | Unintentional interference suppression method for airborne radar under complex interference scene |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102904604A (en) * | 2012-09-10 | 2013-01-30 | 北京大学 | Narrow-band interference suppression method and device |
CN103116155A (en) * | 2012-09-28 | 2013-05-22 | 北京理工大学 | Homotype radar same frequency interference suppression method used for ship formation condition |
CN103439691A (en) * | 2013-08-30 | 2013-12-11 | 西安电子科技大学 | Method for broadband networking radar to restrict narrow-band interference |
US20160072547A1 (en) * | 2014-03-31 | 2016-03-10 | King Fahd University Of Petroleum And Minerals | Evaluation of compressed sensing in uwb systems with nbi |
CN105929371A (en) * | 2016-04-22 | 2016-09-07 | 西安电子科技大学 | Airborne radar clutter suppression method based on covariance matrix estimation |
CN106546965A (en) * | 2016-10-31 | 2017-03-29 | 西安电子科技大学 | Based on radar amplitude and the space-time adaptive processing method of Doppler-frequency estimation |
-
2018
- 2018-05-31 CN CN201810548240.3A patent/CN108957419B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102904604A (en) * | 2012-09-10 | 2013-01-30 | 北京大学 | Narrow-band interference suppression method and device |
CN103116155A (en) * | 2012-09-28 | 2013-05-22 | 北京理工大学 | Homotype radar same frequency interference suppression method used for ship formation condition |
CN103439691A (en) * | 2013-08-30 | 2013-12-11 | 西安电子科技大学 | Method for broadband networking radar to restrict narrow-band interference |
US20160072547A1 (en) * | 2014-03-31 | 2016-03-10 | King Fahd University Of Petroleum And Minerals | Evaluation of compressed sensing in uwb systems with nbi |
CN105929371A (en) * | 2016-04-22 | 2016-09-07 | 西安电子科技大学 | Airborne radar clutter suppression method based on covariance matrix estimation |
CN106546965A (en) * | 2016-10-31 | 2017-03-29 | 西安电子科技大学 | Based on radar amplitude and the space-time adaptive processing method of Doppler-frequency estimation |
Non-Patent Citations (5)
Title |
---|
A.M. KUZMINSKIY等: ""Asynchronous interference cancellation with an antenna array"", 《THE 13TH IEEE INTERNATIONAL SYMPOSIUM ON PERSONAL, INDOOR AND MOBILE RADIO COMMUNICATIONS》 * |
KAZUMA ANDO等: ""Interference detection performance using asynchronous MU-MIMO and self-interference cancellation technique"", 《2016 INTERNATIONAL SYMPOSIUM ON ANTENNAS AND PROPAGATION (ISAP)》 * |
付启众等: ""一种基于信号处理的异步干扰消除方法"", 《雷达科学与技术》 * |
周昆正: ""基于实时优化的搜索雷达点迹提取方法"", 《现代导航》 * |
安红等: "《雷达电子战系统建模与仿真》", 31 December 2017 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110618411A (en) * | 2019-10-23 | 2019-12-27 | 电子科技大学 | Airborne radar clutter real-time signal generation method |
CN110618411B (en) * | 2019-10-23 | 2022-11-04 | 电子科技大学 | Airborne radar clutter real-time signal generation method |
CN111007466A (en) * | 2019-12-03 | 2020-04-14 | 西安电子科技大学 | DCAR interference suppression method and system based on introduced range gate degree of freedom |
CN111007466B (en) * | 2019-12-03 | 2023-05-26 | 西安电子科技大学 | DCAR interference suppression method and system based on introduced range gate degree of freedom |
CN112014806A (en) * | 2020-08-14 | 2020-12-01 | 西安电子科技大学 | Method for suppressing unintentional interference of airborne radar in complex interference scene |
CN112014806B (en) * | 2020-08-14 | 2024-04-23 | 西安电子科技大学 | Unintentional interference suppression method for airborne radar under complex interference scene |
RU2751532C1 (en) * | 2021-04-06 | 2021-07-14 | Акционерное общество "Всероссийский научно-исследовательский институт радиотехники" (АО "ВНИИРТ") | Inter-review device for mapping asynchronous impulse interference for pulse-doppler radar stations and method for its implementation |
CN114690143A (en) * | 2022-04-19 | 2022-07-01 | 珠海微度芯创科技有限责任公司 | Method and device for suppressing radar clutter, radar and medium |
CN114690143B (en) * | 2022-04-19 | 2023-01-06 | 珠海微度芯创科技有限责任公司 | Method and device for suppressing radar clutter, radar and medium |
CN116224249A (en) * | 2023-03-29 | 2023-06-06 | 上海雷骥电子科技有限公司 | Doppler frequency width acquisition method for main lobe clutter region of airborne radar |
CN116224249B (en) * | 2023-03-29 | 2023-11-21 | 上海雷骥电子科技有限公司 | Doppler frequency width acquisition method for main lobe clutter region of airborne radar |
Also Published As
Publication number | Publication date |
---|---|
CN108957419B (en) | 2022-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108957419A (en) | Asynchronous interference suppressing method based on notch filter processing | |
US9157992B2 (en) | Knowledge aided detector | |
CN104914415B (en) | Single-pulse radar coherent jamming method based on target range profile template matching | |
US6750804B2 (en) | System and method for detecting and estimating the direction of near-stationary targets in monostatic clutter using phase information | |
US8970426B1 (en) | Automatic matched Doppler filter selection | |
CN105785340B (en) | The angle measurement and recognition methods that monopulse radar interferes target in main lobe and bait | |
RU2280263C1 (en) | Method for selecting air decoys | |
CN103018727A (en) | Sample-training-based non-stationary clutter suppression method of vehicle-mounted radar | |
CN101881822B (en) | Method for inhibiting same frequency interference of shared-spectrum radars | |
CN108614245B (en) | Radio frequency stealth-based interference method while tracking | |
CN106383339B (en) | A kind of mirror target suppressing method of multi-site radar signal grade joint-detection | |
US9482744B1 (en) | Staggered pulse repetition frequency doppler processing | |
RU2704789C1 (en) | Method for adaptive signal processing in survey coherent-pulse radar stations | |
CN109283497B (en) | Bistatic FDA-MIMO radar range deception jamming identification method | |
CN104535972A (en) | Coherent transferring type jamming inhibition method for airborne radar | |
CN102288949A (en) | Optimal processor based maneuvering target detection method | |
Geng et al. | Ground moving target detection using beam-Doppler image feature recognition | |
WO2010039299A1 (en) | Counter target acquisition radar and acoustic adjunct for classification | |
RU2419107C1 (en) | Method of selecting moving targets in pulse-wise carrier frequency tuning | |
US9995819B1 (en) | System and method for spherical object identification using radar | |
CN109581366B (en) | Discrete sidelobe clutter identification method based on target steering vector mismatch | |
US9903944B2 (en) | Target detection system and method | |
CN108508413B (en) | Target detection method based on probability statistics under low signal-to-noise ratio condition | |
CN106054142B (en) | A kind of airborne MIMO radar main lobe smart munition suppressing method and system | |
Aboutanios et al. | Evaluation of the single and two data set STAP detection algorithms using measured data |
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 |