CN103605114A - Non-broadside array airborne radar short range clutter suppression method based on multiple frequencies - Google Patents
Non-broadside array airborne radar short range clutter suppression method based on multiple frequencies Download PDFInfo
- Publication number
- CN103605114A CN103605114A CN201310646050.2A CN201310646050A CN103605114A CN 103605114 A CN103605114 A CN 103605114A CN 201310646050 A CN201310646050 A CN 201310646050A CN 103605114 A CN103605114 A CN 103605114A
- Authority
- CN
- China
- Prior art keywords
- range
- data
- unit
- clutter
- short range
- 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
Images
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/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
- G01S7/2928—Random or non-synchronous interference pulse cancellers
-
- 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/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a non-broadside array airborne radar short range clutter suppression method based on multiple frequencies. The non-broadside array airborne radar short range clutter suppression method is implement in the following steps: (1) a distance unit range and a Doppler value range of strong short range clutter are acquired; (2) echo data of all distance units received by a radar are subjected to pitching sum and difference beam forming to obtain sum and difference beam data; (3) the sum and difference beam data within the distance unit range of the strong short range clutter are subjected to smoothing to obtain a smoothing result of the distance units within the distance unit range; (4) the sum and difference beam data beyond the distance unit range are subjected to pitching beam forming to obtain a smoothing result of the distance units beyond the distance unit range; (5) the smoothing result in the step (3) and the smoothing result in the step (4) are integrated to obtain a smoothing result of all the distance units, and then short range clutter suppression is completed. By means of the non-broadside array airborne radar short range clutter suppression method based on the multiple frequencies, the distance stability of a radar echo is improved, the clutter suppression performance of a radar system is improved, the amount of calculation is lowered, and the non-broadside array airborne radar short range clutter suppression method can be applied to target detection.
Description
Technical field
The invention belongs to signal processing technology field, further relate to radar clutter suppression method, the short range clutter distributing in the space-time two-dimensional plane that can be used for suppressing to cause when airborne radar antenna non-working side is settled.
Background technology
In airborne radar signal is processed, because the irradiation of radar signal main lobe can produce very strong land clutter, especially, when airborne radar antenna non-working side is settled, the Range-dependent that clutter presents in space-time two-dimensional distribution can cause stronger short range clutter.Specifically, between each range unit of short range, luffing angle changes greatly, cause the Doppler frequency of clutter can produce along with the variation of range unit comparatively strong variation, so just produced short range clutter, the destruction that the appearance of short range clutter is serious the stationarity of clutter, make to obtain abundant independent same distribution sample, cause traditional space-time adaptive to process the clutter rejection degradation of (STAP) technology.
The method of eliminating at present clutter Range-dependent mainly contains the full spectral domain algorithm of compensation class algorithm, change of scale class algorithm, interpolation class.
Compensation class algorithm is that each training sample is compensated respectively, and it mainly comprises: the Doppler effect correction that calendar year 2001 Kreyenkamp O and Klemm R propose; Calendar year 2001 Pearson F and Borsari G propose high-order Doppler effect correction; The self-adaptation angle Doppler effect correction that Melvin W and Davis M in 2007 proposes on periodical Aerospace and Electronic Systems; The angle Doppler effect correction that Colone F in 2011 proposes, this class algorithm all can compensate clutter Range-dependent at clutter without range ambiguity in the situation that effectively, but when clutter exists range ambiguity, this class algorithm is no longer applicable;
2011 by Ries P, the method that Lapierre F D and Verly J G adopts spectrum to estimate at the upper change of scale class algorithm proposing of periodical Aerospace and Electronic Systems (AES), the clutter covariance matrix of reconstruct unit to be detected, then according to reconstruct clutter covariance matrix out, the clutter in training sample is transformed in the clutter space of sample to be detected, this class algorithm can effectively compensate the Range-dependent of clutter when clutter exists range ambiguity, but this class algorithm operation quantity is larger, is not easy to practical application;
The full spectral domain algorithm of the interpolation class that Liu Jinhui etc. propose in electronics and information journal adopts least square technology the clutter of training sample to be transformed in the clutter space of sample to be detected, also can effectively compensate the Range-dependent of clutter, but steering vector in the time of need to estimating clutter empty, but due to the reason of error estimate empty time steering vector and actual when empty steering vector have the deviation that can not estimate, so this class algorithm clutter rejection when there is array error can decline greatly.
Summary of the invention
The object of the invention is to the deficiency for above-mentioned prior art, a kind of non-positive side battle array airborne radar short range clutter suppression method based on multiple frequency has been proposed, to avoid the Range-dependent of training sample, improve the distance stationarity of radar return, improve the clutter rejection of radar system, reduce operand.
Technical scheme of the present invention is achieved in that
One. know-why
For airborne radar, normally by launching continuously the waveform of several groups of different pulse repetition, solve the problem that target range is fuzzy, in certain time interval, because the short range clutter of same distance unit in these several groups of different pulse repetition echoes has higher correlativity, therefore can estimate the covariance matrix of clutter as training sample by the data of same distance unit in these several groups of different pulse repetition echoes.For the less problem of training sample, can, by adopting pitching and the poor adaptive processing technique of two-freedom, reduce the needed training sample number of accurate estimation clutter covariance matrix, by being filtered the inhibition of non-positive side battle array airborne radar clutter in pairs.
Two. technical scheme
The technical step of realizing the object of the invention comprises as follows:
(1) obtain range unit scope and Doppler's span of the strong clutter of short range:
1a) according to pulse repetition rate maximum in several pulse repetition raties of radar emission and the flying height of bandwidth and radar, calculate the not numerical value of fuzzy distance unit, ultimate range unit numerical value L as the strong clutter of short range, and calculating altitude line range unit H, the range unit scope that obtains the strong clutter of short range is H~L;
1b) calculate in the range unit scope H~L of the strong clutter of short range Doppler's span of the strong clutter of short range in each range unit main cone face;
(2) echo data of all range units of I group different pulse repetition radar being received carries out pitching and difference beam forms, obtain all pulse repetition raties all range units to be processed with difference beam data X
1, X
2... X
i... X
i, i=1 wherein, 2...I, I is the number of the data of the different pulse repetition that receives of radar;
(3) to carrying out filtering with difference beam data in the range unit scope H~L of the strong clutter of short range:
3a) from all range units to be processed with difference beam data X
1, X
2... X
i... X
iin, the data x of l range unit within the scope of taking-up H~L range unit
1, l, x
2, l... x
i,l... x
i,l;
3b) from step 1b) take out the doppler values in the middle of being positioned in Doppler's span of the strong clutter of short range of l range unit calculating, with this doppler values formation time domain steering vector v
1, l, v
2, l... v
i,l... v
i,l, then use this time domain steering vector respectively with the data x of l range unit
1, l, x
2, l... x
i,l... x
i,lmultiply each other, by the result multiplying each other, carry out the covariance matrix of short range clutter, then with the covariance matrix of this short range clutter, carry out the calculating of adaptive weight;
3c) use the data x of weights to l range unit
1, l, x
2, l... x
i,l... x
i,lcarry out filtering, obtain the filtering result of the data of l range unit: y
1, l, y
2, l... y
i,l... y
i,l;
3d) repeat 3b)-3c), until obtain all range units with difference beam data X
1, X
2... X
i... X
ibe positioned at the data filtering result of H~L range unit scope;
(4) to outside the range unit scope H~L of the strong clutter of short range to carry out pitching wave beam synthetic with difference beam data:
4a) from i all range units of pulse repetition rate to be processed with difference beam data X
ithe data of the extraneous range unit of middle taking-up H~L range unit, carry out pitching wave beam synthetic, obtain i pulse repetition rate with difference beam data X
iin be positioned at the filtering result of the extraneous range unit of H~L range unit, i=1 wherein, 2...I, I is the number of the data of the different pulse repetition that receives of radar;
4b) repeat 4a), until obtain all pulse repetition raties with difference beam data X
1, X
2... X
i... X
ibe positioned at the filtering result of the extraneous data of H~L range unit;
(5) integral data, completes short range clutter and suppresses:
5a) get the filtering result of i all range units of pulse repetition rate in step 3) and step 4), according to range unit sequence number, carry out Data Integration, obtain the integration kalman filter result Y of i pulse repetition rate
i;
5b) repeating step 5a), obtain the integration kalman filter result Y of all pulse repetition raties
1, Y
2... Y
i... Y
i, complete the inhibition to airborne radar short range clutter.
The present invention compared with prior art, has the following advantages:
1. the present invention compares along the method for choosing sample apart from dimension with common, because the echo data of choosing a plurality of pulse repetition rate same distance unit is as sample, can avoid the Range-dependent problem of training sample, can estimate comparatively exactly clutter covariance matrix, improve the distance stationarity of radar return, improve the clutter rejection of radar system.
2. the present invention is owing to adopting pitching and the poor adaptive processing technique of two-freedom to suppress the strong clutter of short range, can reduce the requirement to training sample number, reduce operand simultaneously, the difficulty of estimate covariance is further reduced, better suppressed short range clutter.
Accompanying drawing explanation
Fig. 1 is realization flow figure of the present invention;
Fig. 2 is not for being used the range Doppler image of partial distance unit of the radar return of the inventive method;
Fig. 3 is for being used the inventive method to process the range Doppler image of the partial distance unit of radar return later;
Fig. 4 is the partial distance cell power figure of target place Doppler's passage from Fig. 2 and Fig. 3.
Embodiment
With reference to Fig. 1, implementation step of the present invention is as follows:
Step 1, obtains range unit scope and Doppler's span of the strong clutter of short range.
1a) according to pulse repetition rate f maximum in several pulse repetition raties of radar emission
maxwith the flying height h of bandwidth B and radar, calculate the not numerical value of fuzzy distance unit
as the ultimate range unit numerical value L of the strong clutter of short range, and according to altitude line range unit H, the range unit scope that obtains the strong clutter of short range is H~L;
The strong clutter of short range that 1b) calculates l range unit in the strong clutter range unit of short range scope H~L is at pitching main cone face interior orientation angle θ
lspan be:
Wherein, θ
0the position angle of radar main beam pointing, φ
0for the angle of pitch of radar main beam pointing, θ
3dBfor main beam 3dB width,
be l the angle of pitch that range unit is corresponding;
1c) by the strong clutter of short range of l range unit at pitching main cone face interior orientation angle θ
lspan, the span of calculating l the strong clutter Doppler of short range in range unit main cone face is:
Wherein, v is the flying speed of radar, and λ is radar wavelength, φ
lthe angle of pitch of l range unit, θ
αfor the angle of antenna array axially and between carrier aircraft heading;
1d) repeat 1b)-1c), obtain the span of the strong clutter Doppler of short range in all range unit main cone faces in the range unit scope H~L of the strong clutter of short range.
Step 2, the echo data that radar is received carries out pitching and difference beam forms.
The echo data that takes out all range units of I group different pulse repetition from the data that radar receives, carries out pitching and difference beam and forms, obtain all pulse repetition raties all range units to be processed with difference beam data X
1, X
2... X
i... X
i, i=1 wherein, 2...I, I is the number of the data of the different pulse repetition that receives of radar;
Step 3, to carrying out filtering with difference beam data in the range unit scope H~L of the strong clutter of short range.
3a) from all range units to be processed with difference beam data X
1, X
2... X
i... X
iin, the data x of l range unit within the scope of taking-up H~L range unit
1, l, x
2, l... x
i,l... x
i,l, i=1,2...I;
3b) from step 1b) take out the doppler values in the middle of being positioned in Doppler's span of the strong clutter of short range of l range unit calculating, with this doppler values formation time domain steering vector v
1, l, v
2, l... v
i,l... v
i,l, each time domain steering vector is:
Wherein, f
ithe size of i pulse repetition rate, f
lthe doppler values that is positioned at centre taking out in Doppler's span of the strong clutter of short range of l range unit determining in step 1, ()
trepresent transposition, K
ithe echo data of value and i pulse repetition rate in coherent pulse number identical, I is the number of the data of the different pulse repetition that receives of radar;
3c) from the data x of l range unit
1, l, x
2, l... x
i,l... x
i,lthe data x of middle taking-up m array unit
1, l, m, x
2, l, m... x
i, l, m... x
i, l, m, m=1 wherein, 2...M, M is the columns of radar array element;
3d) use the data x of the m array unit taking out
1, l, m, x
2, l, m... x
i, l, m... x
i, l, mrespectively with time domain steering vector v
1, l, v
2, l... v
i,l... v
i,leach carry out correspondence and multiply each other, obtain the multiplied result c of m array unit
1, l, m, c
2, l, m... c
i, l, m... c
i, l, m, c wherein
i, l, m=x
i, l, m* v
i,l, i=1,2 ... I;
3e) utilize the multiplied result c of m array unit
1, l, m, c
2, l, m... c
i, l, m... c
i, l, m, calculate the covariance matrix R of the short range clutter of m array unit
m,l:
3f) take out the covariance matrix R of the short range clutter of m array unit
m,l, utilize following formula to calculate l range unit data x
1, l, x
2, l... x
i,l... x
i,lin the weight w of m array unit
m,l:
3g) take out the weight w of m array unit
m,l, respectively with the data x of m array unit
1, l, m, x
2, l, m... x
i, l, m... x
i, l, mmultiply each other, obtain m array metadata x
1, l, m, x
2, l, m... x
i, l, m... x
i, l, mfiltering result y
1, l, m, y
2, l, m... y
i, l, m... y
i, l, m;
3h) repetition 3c)-3g) obtain the data x of l range unit
1, l, x
2, l... x
i,l... x
i,lthe filtering result of all M row: y
1, l, y
2, l... y
i,l... y
i,l;
3i) repeat 3a)-3h), until obtain all range units with difference beam data X
1, X
2... X
i... X
ibe positioned at the data filtering result of H~L range unit scope.
Step 4, carries out pitching wave beam to the extraneous and difference beam data of H~L range unit synthetic.
4a) from i all range units of pulse repetition rate to be processed with difference beam data X
ithe data of the extraneous range unit of middle taking-up H~L range unit, carry out pitching wave beam synthetic, obtain i pulse repetition rate with difference beam data X
iin be positioned at the filtering result of the extraneous range unit of H~L range unit, i=1 wherein, 2...I, I is the number of the data of the different pulse repetition that receives of radar;
4b) repeating step 4a), until obtain all pulse repetition raties with difference beam data X
1, X
2... X
i... X
ibe positioned at the extraneous data filtering result of H~L range unit.
Step 5, integral data, completes short range clutter and suppresses.
5a) get the filtering result of i all range units of pulse repetition rate in step 3) and step 4), according to range unit sequence number, carry out Data Integration, obtain the integration kalman filter result Y of i pulse repetition rate
i;
5b) repeating step 5a), obtain the integration kalman filter result Y of all pulse repetition raties
1, Y
2... Y
i... Y
i, complete the inhibition to airborne radar short range clutter.
Effect of the present invention can be processed and be further illustrated by following emulated data:
1. simulation parameter
Carrier aircraft flying speed 100m/s, flying height 4150m, 32 of coherent pulse numbers, transmitted bandwidth 2.5MHz, operation wavelength 0.1m, 15 of orientation element number of array, 8 of pitching element number of array, array element interval, orientation and pitching array element interval are 0.05m, main beam azimuth 60 degree, the main beam angle of pitch 0 degree, five kinds of pulse repetition raties are used respectively 7337Hz, 6846Hz, 6533Hz, 5979Hz, 5200Hz.
According to systematic parameter emulation the land clutter data of five kinds of pulse repetition raties of forward sight battle array radar, the miscellaneous noise ratio of every array element, every pulse is 50dB.In main beam, add a target, apart from radar 393.48Km, the flying speed of radar is 18m/s relatively, and target is positioned at the 79th range unit, and its Doppler frequency is 360Hz.
2. emulation content and interpretation of result
2.1) echo data of above-mentioned 5 kinds of pulse repetition raties is directly carried out to wave beam and synthesize, then the data that are 7337Hz from pulse repetition rate, take out front 150 range unit data, draw its range Doppler figure, as Fig. 2.
2.2) use the inventive method to process the echo data of 5 kinds of pulse repetition raties, from process data later, take out front 150 range unit data that pulse repetition rate is the data of 7337Hz, draw its range Doppler figure, as Fig. 3.
2.3) from Fig. 2 and Fig. 3, take out respectively the range unit power of 150 points that the Doppler frequency at target place is 360Hz, draw out 150 whole power diagrams that range unit is corresponding, as Fig. 4.
Comparison diagram 2 and Fig. 3, can significantly see and in Fig. 3, use the short range clutter of the inventive method to be inhibited, and is subject to the impact of the strong clutter of short range in Fig. 2, cannot detect target; After short range clutter is inhibited in Fig. 3, can detect target.
As can see from Figure 4, after the inventive method filtering, short range clutter is effectively suppressed, and target reveals completely from clutter.
Claims (5)
1. the non-positive side battle array airborne radar short range clutter suppression method based on multiple frequency, comprises the steps:
(1) obtain range unit scope and Doppler's span of the strong clutter of short range:
1a) according to pulse repetition rate maximum in several pulse repetition raties of radar emission and the flying height of bandwidth and radar, calculate the not numerical value of fuzzy distance unit, ultimate range unit numerical value L as the strong clutter of short range, and calculating altitude line range unit H, the range unit scope that obtains the strong clutter of short range is H~L;
1b) calculate in the range unit scope H~L of the strong clutter of short range Doppler's span of the strong clutter of short range in each range unit main cone face;
(2) echo data of all range units of I group different pulse repetition radar being received carries out pitching and difference beam forms, obtain all pulse repetition raties all range units to be processed with difference beam data X
1, X
2... X
i... X
i, i=1 wherein, 2...I, I is the number of the data of the different pulse repetition that receives of radar;
(3) to carrying out filtering with difference beam data in the range unit scope H~L of the strong clutter of short range:
3a) from all range units to be processed with difference beam data X
1, X
2... X
i... X
iin, the data x of l range unit within the scope of taking-up H~L range unit
1, l, x
2, l... x
i,l... x
i,l;
3b) from step 1b) take out the doppler values in the middle of being positioned in Doppler's span of the strong clutter of short range of l range unit calculating, with this doppler values formation time domain steering vector v
1, l, v
2, l... v
i,l... v
i,l, then use this time domain steering vector respectively with the data x of l range unit
1, l, x
2, l... x
i,l... x
i,lmultiply each other, by the result multiplying each other, carry out the covariance matrix of short range clutter, then with the covariance matrix of this short range clutter, carry out the calculating of adaptive weight;
3c) use the data x of weights to l range unit
1, l, x
2, l... x
i,l... x
i,lcarry out filtering, obtain the filtering result of the data of l range unit: y
1, l, y
2, l... y
i,l... y
i,l;
3d) repeat 3b)-3c), until obtain all range units with difference beam data X
1, X
2... X
i... X
ibe positioned at the data filtering result of H~L range unit scope;
(4) to outside the range unit scope H~L of the strong clutter of short range to carry out pitching wave beam synthetic with difference beam data:
4a) from i all range units of pulse repetition rate to be processed with difference beam data X
ithe data of the extraneous range unit of middle taking-up H~L range unit, carry out pitching wave beam synthetic, obtain i pulse repetition rate with difference beam data X
iin be positioned at the filtering result of the extraneous range unit of H~L range unit, i=1 wherein, 2...I, I is the number of the data of the different pulse repetition that receives of radar;
4b) repeat 4a), until obtain all pulse repetition raties with difference beam data X
1, X
2... X
i... X
ibe positioned at the filtering result of the extraneous data of H~L range unit;
(5) integral data, completes short range clutter and suppresses:
5a) get the filtering result of i all range units of pulse repetition rate in step 3) and step 4), according to range unit sequence number, carry out Data Integration, obtain the integration kalman filter result Y of i pulse repetition rate
i;
5b) repeating step 5a), obtain the integration kalman filter result Y of all pulse repetition raties
1, Y
2... Y
i... Y
i, complete the inhibition to airborne radar short range clutter.
2. the non-positive side battle array airborne radar short range clutter suppression method based on multiple frequency according to claim 1, step 1b wherein) described calculating is in the range unit scope H~L of the strong clutter of short range, Doppler's span of the strong clutter of short range in each range unit main cone face, carry out as follows:
The strong clutter of short range that 1b1) calculates l range unit in the strong clutter range unit of short range scope H~L is at pitching main cone face interior orientation angle θ
lspan be:
Wherein, θ
0the position angle of radar main beam pointing, φ
0for the angle of pitch of radar main beam pointing, θ
3dBfor main beam 3dB width,
be l the angle of pitch that range unit is corresponding;
1b2) by the strong clutter of short range of l range unit at pitching main cone face interior orientation angle θ
lspan, the span of calculating l the strong clutter Doppler of short range in range unit main cone face is:
Wherein, v is the flying speed of radar, and λ is radar wavelength, φ
lthe angle of pitch of l range unit, θ
αfor the angle of antenna array axially and between carrier aircraft heading.
3. the non-positive side battle array airborne radar short range clutter suppression method based on multiple frequency according to claim 1, wherein step 3b) described in the middle doppler values of use form time domain steering vector v
1, l, v
2, l... v
i,l... v
i,l, by following formula, undertaken:
Wherein, f
ithe size of i pulse repetition rate, f
lthe doppler values that is positioned at centre taking out in Doppler's span of the strong clutter of short range of l range unit determining in step 1, ()
trepresent transposition, K
ithe echo data of value and i pulse repetition rate in coherent pulse number identical, I is the number of the data of the different pulse repetition that receives of radar.
4. the non-positive side battle array airborne radar short range clutter suppression method based on multiple frequency according to claim 1, wherein said step 3b) with time domain steering vector v
1, l, v
2, l... v
i,l... v
i,lrespectively with the data x of l range unit
1, l, x
2, l... x
i,l... x
i,lmultiply each other, by the result multiplying each other, carry out the covariance matrix of short range clutter, carry out as follows:
(3b1) from the data x of l range unit within the scope of H~L range unit
1, l, x
2, l... x
i,l... x
i,lthe data x of the m array unit of middle taking-up
1, l, m, x
2, l, m... x
i, l, m... x
i, l, m, m=1 wherein, 2...M, M is the columns of radar array element;
(3b2) use the data x of the m array unit taking out
1, l, m, x
2, l, m... x
i, l, m... x
i, l, mrespectively with time domain steering vector v
1, l, v
2, l... v
i,l... v
i,leach carry out correspondence and multiply each other, obtain the multiplied result c of m array unit
1, l, m, c
2, l, m... c
i, l, m... c
i, l, m, c wherein
i, l, m=x
i, l, m* v
i,l, i=1,2 ... I;
(3b3) utilize the multiplied result c of m array unit
1, l, m, c
2, l, m... c
i, l, m... c
i, l, m, calculate the covariance matrix R of the short range clutter of m array unit
m,l:
Wherein, ()
hrepresent conjugate transpose, i=1 wherein, 2 ... I, I is the number of the data of the different pulse repetition that receives of radar;
(3b4) repetition (3b1)-(3b3), obtain the data x of l range unit within the scope of H~L range unit
1, l, x
2, l... x
i,l... x
i,lin the covariance matrix R of short range clutter of all M array unit
1, l, R
2, l... R
m,l... R
m,l.
5. the non-positive side battle array airborne radar short range clutter suppression method based on multiple frequency according to claim 1, wherein said step 3b), with the covariance matrix of short range clutter, carry out the calculating of adaptive weight, carry out as follows:
(3ba) from the data x of l range unit within the scope of H~L range unit
1, l, x
2, l... x
i,l... x
i,lthe data x of the m array unit of middle taking-up
1, l, m, x
2, l, m... x
i, l, m... x
i, l, m, m=1 wherein, 2...M, M is the columns of radar array element;
(3bb) take out the covariance matrix R of the short range clutter of m array unit
m,l, utilize following formula to calculate l range unit data x
1, l, x
2, l... x
i,l... x
i,lin the weight w of m array unit
m,l:
(3bc) repetition (3ba)-(3bb), obtain the data x of l range unit
1, l, x
2, l... x
i,l... x
i,lin the weight w of all M array unit
1, l, w
2, l... w
m,l... w
m,l.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310646050.2A CN103605114B (en) | 2013-12-03 | 2013-12-03 | Non-broadside array airborne radar short range clutter suppression method based on multiple frequencies |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310646050.2A CN103605114B (en) | 2013-12-03 | 2013-12-03 | Non-broadside array airborne radar short range clutter suppression method based on multiple frequencies |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103605114A true CN103605114A (en) | 2014-02-26 |
CN103605114B CN103605114B (en) | 2015-07-15 |
Family
ID=50123354
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310646050.2A Active CN103605114B (en) | 2013-12-03 | 2013-12-03 | Non-broadside array airborne radar short range clutter suppression method based on multiple frequencies |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103605114B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103926572A (en) * | 2014-03-28 | 2014-07-16 | 西安电子科技大学 | Clutter rejection method of self-adaption subspace for non-sidelooking airborne array radar |
CN104635214A (en) * | 2015-02-12 | 2015-05-20 | 西安电子科技大学 | Airborne forward-looking frequency diversity array radar distance fuzzy clutter suppression method |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106940442B (en) * | 2017-03-15 | 2019-09-10 | 西北大学 | A kind of high speed missile-borne radar waveform design method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101561497A (en) * | 2009-05-22 | 2009-10-21 | 西安电子科技大学 | Airborne radar clutter suppression method |
WO2012067557A1 (en) * | 2010-11-19 | 2012-05-24 | Saab Ab | A method and radar system for repetition jammer and clutter suppression |
CN102721947A (en) * | 2012-06-14 | 2012-10-10 | 河海大学 | Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar |
CN102788970A (en) * | 2011-12-22 | 2012-11-21 | 河海大学 | High-efficiency self-adapting suppression method for short-range clutters of airborne digital array radar |
-
2013
- 2013-12-03 CN CN201310646050.2A patent/CN103605114B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101561497A (en) * | 2009-05-22 | 2009-10-21 | 西安电子科技大学 | Airborne radar clutter suppression method |
WO2012067557A1 (en) * | 2010-11-19 | 2012-05-24 | Saab Ab | A method and radar system for repetition jammer and clutter suppression |
CN102788970A (en) * | 2011-12-22 | 2012-11-21 | 河海大学 | High-efficiency self-adapting suppression method for short-range clutters of airborne digital array radar |
CN102721947A (en) * | 2012-06-14 | 2012-10-10 | 河海大学 | Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar |
Non-Patent Citations (1)
Title |
---|
李明等: "机载前视阵雷达距离模糊杂波抑制新方法", 《***仿真学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103926572A (en) * | 2014-03-28 | 2014-07-16 | 西安电子科技大学 | Clutter rejection method of self-adaption subspace for non-sidelooking airborne array radar |
CN103926572B (en) * | 2014-03-28 | 2016-06-29 | 西安电子科技大学 | A kind of clutter suppression method of battle array radar self adaptation subspace, airborne anon-normal side |
CN104635214A (en) * | 2015-02-12 | 2015-05-20 | 西安电子科技大学 | Airborne forward-looking frequency diversity array radar distance fuzzy clutter suppression method |
CN104635214B (en) * | 2015-02-12 | 2017-01-04 | 西安电子科技大学 | Air-borne Forward-looking frequency diversity array radar range ambiguity clutter suppression method |
Also Published As
Publication number | Publication date |
---|---|
CN103605114B (en) | 2015-07-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105738879B (en) | Radar clutter space-time adaptive pre-filtering method based on sparse recovery | |
Chen et al. | Radon-fractional ambiguity function-based detection method of low-observable maneuvering target | |
CN103399310B (en) | Method for detecting radar weak moving target based on PD (Phase Differentiation) RLVD (Radon-Lv Distribution) | |
CN106093870B (en) | The SAR-GMTI clutter suppression methods of hypersonic aircraft descending branch | |
CN102721947B (en) | Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar | |
CN104502898B (en) | The maneuvering target method for parameter estimation that modified R FT and amendment MDCFT are combined | |
CN104977571B (en) | Range ambiguity clutter suppression method based on pitching frequency diversity STAP | |
CN103954941B (en) | The multipulse cognitive clutter suppression method of a kind of airborne phased array radar two dimension | |
CN104833972B (en) | A kind of bistatic CW with frequency modulation synthetic aperture radar frequency becomes mark imaging method | |
CN104035095A (en) | Low-altitude wind shear wind speed estimating method based on space-time optimal processor | |
CN103969629A (en) | Airborne radar clutter self-adaption restraining method based on main-lobe clutter registering | |
CN104345300B (en) | The airborne non-working side battle array radar STAP method of clutter space-time spectrum linear compensation | |
CN105403864B (en) | Based on the two-dimentional boat-carrying high-frequency ground wave radar ocean clutter cancellation method for improving oblique projection | |
CN109581362B (en) | Signal processing method of synthetic aperture radar altimeter in variable pulse cluster mode | |
CN104849708B (en) | High speed machine moving target parameter estimation method based on the conversion of frequency domain polynomial-phase | |
CN103176168A (en) | Short-range cluster cancellation method for airborne non-side-looking array radar | |
CN112612006B (en) | Deep learning-based non-uniform clutter suppression method for airborne radar | |
CN103399309A (en) | Space-time two-dimensional clutter spectrum estimation method based on iterative weighted minimum variance | |
CN103412309B (en) | Move constant bistatic forward sight synthetic-aperture radar NLCS formation method | |
CN108931766A (en) | A kind of non-homogeneous STAP jamming target filtering method based on sparse reconstruct | |
CN103901416A (en) | Multichannel clutter suppression method based on robust principal component analysis | |
CN104076337A (en) | Airborne radar clutter suppression method based on array element amplitude and phase error correction | |
CN104950295A (en) | High-speed maneuvering target detecting method based on correlation functions and scale changes | |
CN104375128B (en) | Fast high maneuvering target accumulating and detecting method based on cross-correlation functions | |
CN103605114B (en) | Non-broadside array airborne radar short range clutter suppression method based on multiple frequencies |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |