CN102721947B - Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar - Google Patents

Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar Download PDF

Info

Publication number
CN102721947B
CN102721947B CN 201210195856 CN201210195856A CN102721947B CN 102721947 B CN102721947 B CN 102721947B CN 201210195856 CN201210195856 CN 201210195856 CN 201210195856 A CN201210195856 A CN 201210195856A CN 102721947 B CN102721947 B CN 102721947B
Authority
CN
China
Prior art keywords
pitching
pulse
clutter
array
dimensionality reduction
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.)
Expired - Fee Related
Application number
CN 201210195856
Other languages
Chinese (zh)
Other versions
CN102721947A (en
Inventor
沈明威
蒋德富
静大海
冯芸
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN 201210195856 priority Critical patent/CN102721947B/en
Publication of CN102721947A publication Critical patent/CN102721947A/en
Application granted granted Critical
Publication of CN102721947B publication Critical patent/CN102721947B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses an efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar, and belongs to the field of suppression of clutters of airborne radar. In the method, a cascading mode is adopted for the airborne non-front lateral planar array radar working in a medium and high pulse repetition frequency (PRF) system, and the method includes designing a robust adaptive digital beam former (ADBF) in a pitching dimension of airspace to suppress short-range clusters at first; and designing a two-dimensional dimension-reduction space-time adaptive processor (STAP) in a direction-pulse field to further suppress long-range clutters. By the aid of the efficient three-dimensional space-time adaptive clutter suppression method, the clutters at various fuzzy distances of the airborne non-front lateral planar array radar can be effectively filtered, the convergence speed is fast, computational complexity is low, and the method is suitable forengineering application.

Description

The efficient three-dimensional space-time adaptive clutter suppression method of airborne Digital Array Radar
Technical field
The present invention relates to the efficient three-dimensional space-time adaptive clutter suppression method of a kind of airborne Digital Array Radar, belong to airborne radar clutter and suppress the field.
Background technology
Look work under the airborne radar, the land clutter frequency spectrum is wide, intensity is big, and it is the core means that improve detection performance that clutter is effectively suppressed.Handling (STAP) technology based on the two-dimentional space-time adaptive of orientation, spatial domain and time domain impulse sampling is the effective ways that airborne radar clutter suppresses, and just adopts this technology as the E2-D airborne early warn ing radar of the U.S..STAP is as a kind of statistical adaptive Processing Algorithm, only satisfies under independent same distribution (IID) condition performance at clutter and just can give full play to list of references [1-6].
Airborne radar generally adopts pulse Doppler (PD) system.When radar pulse repetition frequency (PRF) is selected medium-high frequency, will there be range ambiguity, in a range unit, can receive the noise signal of a plurality of rang rings this moment simultaneously, wherein the clutter of fuzzy distance is nearest for the first time, the intensity maximum, list of references [4,14] is defined as the short range clutter, and its characteristic is conducted in-depth analysis.When airborne radar adopts the non-working side battle array, short range clutter and far field clutter spectrum distribution trajector deviation are very big, and spectrum was exceedingly fast with range conversion when the short range clutter was empty, caused the serious non-stationary of noise performance of different distance unit, and the clutter rejection that two-dimentional STAP handles sharply descends.
Compensate heterogeneity between the different distance door at the existing a lot of algorithms of distance non-stationary clutter, as lists of references [7-13] such as self-adaptation angle-Doppler effect correction, registration compensation.But under the range ambiguity condition, the far field clutter has constituted different clutter spectrum with the short range clutter, and backoff algorithm can't be taken into account the clutter of the fuzzy distance of each time simultaneously.Therefore, this type of algorithm can't improve the clutter rejection of follow-up two-dimentional STAP.
Summary of the invention
The present invention be directed under the high PRF system airborne non-positive side battle array Digital Array Radar and adopt traditional orientation-pulse domain two dimension STAP clutter rejection defect of insufficient, based on pitching-orientation-pulse three-dimensional territory the efficient three-dimensional space-time adaptive clutter suppression method of a kind of airborne Digital Array Radar has been proposed, at first suppress the short range clutter at spatial domain pitching dimension design robust ADBF processor, design two-dimentional dimensionality reduction STAP processor in orientation-pulse domain then and suppress long-range clutter.
The present invention adopts following technical scheme for solving its technical matters:
The efficient three-dimensional space-time adaptive clutter suppression method of a kind of airborne Digital Array Radar comprises the steps:
1) pitching robust ADBF adaptive beam forms: to each row pitching linear array, dimensionality reduction forms pitching and wave beam and trim wave beam respectively, namely
Figure 51572DEST_PATH_IMAGE001
（1）
M is antenna pitching dimension element number of array in the formula, Be listed as for every
Figure 240294DEST_PATH_IMAGE003
Individual array element is in the output signal of each pulse,
Figure 647004DEST_PATH_IMAGE004
Be pitching and the wave beam that forms behind M array element dimensionality reduction of every row,
Figure 879403DEST_PATH_IMAGE005
Be the trim wave beam that forms behind M array element dimensionality reduction of every row; Pitching and wave beam
Figure 713366DEST_PATH_IMAGE004
With the trim wave beam
Figure 740970DEST_PATH_IMAGE005
Correspond to pitching every alignment battle array output signal and dimensionality reduction matrix
Figure 623476DEST_PATH_IMAGE006
Product, that is:
Figure 38276DEST_PATH_IMAGE007
（2）
In the formula
Figure 43141DEST_PATH_IMAGE008
Be the output signal of each array element of the every row of pitching in each pulse;
Figure 543393DEST_PATH_IMAGE009
Be the conjugate transpose operational symbol;
Figure 229589DEST_PATH_IMAGE010
Be transpose operator; After the dimension-reduction treatment of pitching spatial domain, adopt first zero filling pulse as training sample in order to the estimation self-adaptive weights, each array element of digital array received in signal and only comprised short range clutter information this moment; When definition the First zero filling reception of impulse signal of the every row pitching of individual range unit antenna array linear array is
Figure 674663DEST_PATH_IMAGE012
, every array unit then
Figure 662211DEST_PATH_IMAGE011
Individual range unit first zero filling pulse pitching and wave beam
Figure 886518DEST_PATH_IMAGE013
Output signal and trim wave beam Output signal is:
Figure 359930DEST_PATH_IMAGE015
（3）
According to the secondary lobe SLC principle that disappears mutually, to
Figure 897090DEST_PATH_IMAGE011
Individual range unit offsets short range clutter in pitching and the wave beam, adaptive weight with the trim wave beam
Figure 925089DEST_PATH_IMAGE016
Satisfy:
Figure 293622DEST_PATH_IMAGE017
（4）
Wherein
Figure 748875DEST_PATH_IMAGE018
Be the mathematical expectation operational symbol,
Figure 445435DEST_PATH_IMAGE019
Be the conjugate transpose operational symbol, then
Figure 73863DEST_PATH_IMAGE020
（5）
Autocorrelation matrix in the formula
Figure 440777DEST_PATH_IMAGE021
Adopt after every row dimension-reduction treatment
Figure 253882DEST_PATH_IMAGE011
Individual range unit and vicinity
Figure 437738DEST_PATH_IMAGE022
The trim beamformer output signal of first zero filling pulse of individual range unit is estimated to obtain, and generally gets
Figure 869857DEST_PATH_IMAGE023
, namely
Figure 26031DEST_PATH_IMAGE024
,
Figure 885403DEST_PATH_IMAGE025
Be
Figure 556556DEST_PATH_IMAGE003
Individual range unit trim beamformer output signal; Cross-correlation matrix
Figure 464469DEST_PATH_IMAGE026
Adopt after every row dimension-reduction treatment
Figure 540397DEST_PATH_IMAGE011
Individual range unit and vicinity
Figure 570670DEST_PATH_IMAGE022
First zero filling pulse pitching of individual range unit and beamformer output signal and trim beamformer output signal are estimated to obtain, and generally get
Figure 666802DEST_PATH_IMAGE023
, namely
Figure 440723DEST_PATH_IMAGE027
,
Figure 368227DEST_PATH_IMAGE028
Be
Figure 569402DEST_PATH_IMAGE003
Individual range unit pitching and beamformer output signal;
By formula (2), (3), (4), (5), to
Figure 152830DEST_PATH_IMAGE011
Each row pitching linear array of individual range unit realizes the adaptive weight that the short range clutter suppresses For:
Figure 512453DEST_PATH_IMAGE030
（6）
Therefore, to
Figure 822211DEST_PATH_IMAGE011
Each alignment battle array of individual range unit adopts adaptive weight
Figure 955252DEST_PATH_IMAGE029
Pulse by Pulse carries out pitching filtering, just can effectively suppress the short range clutter in each pulse signal, avoids the far field echo signal to incur loss simultaneously; Can adopt above-mentioned flow process to calculate corresponding adaptive weight respectively to different range units, and then the first Pulse by Pulse of every array be carried out pitching filtering just can realize the short range clutter in each pulse signal of all range units is suppressed.
2) orientation-pulse domain two dimension dimensionality reduction STAP space-time adaptive is handled: after pitching robust ADBF suppressed the short range clutter, antenna array was equivalent to the orientation linear array; Linear array has N array element when the definition orientation, and each array element receives K pulse signal; When definition the
Figure 62187DEST_PATH_IMAGE011
Individual range unit receives signal
Figure 901967DEST_PATH_IMAGE031
, guiding vector when detecting target empty
Figure 444944DEST_PATH_IMAGE032
（7）
In the formula
Figure 65281DEST_PATH_IMAGE033
,
Figure 719116DEST_PATH_IMAGE034
,
Figure 413403DEST_PATH_IMAGE035
For the array element orientation to spacing,
Figure 392860DEST_PATH_IMAGE036
Be radar wavelength,
Figure 234914DEST_PATH_IMAGE037
Be the target Doppler frequency,
Figure 692440DEST_PATH_IMAGE038
Be the target azimuth incident angle,
Figure 506813DEST_PATH_IMAGE039
Be radar pulse recurrence interval PRI, For Kronecker amasss;
Further suppress to remain long-range clutter based on local Combined Treatment algorithm design orientation-pulse domain two dimension dimensionality reduction STAP processor, the spatial domain, orientation is pointed to both sides at the detection object beam and is respectively got 2 adjacent beams, time domain is respectively got 1 adjacent Doppler unit in detection both sides, Doppler unit, and then the dimensionality reduction matrix is
Figure 723872DEST_PATH_IMAGE041
（8）
In the formula
Figure 985089DEST_PATH_IMAGE042
Be spatial domain adjacent beams interval,
Figure 653968DEST_PATH_IMAGE043
Be the adjacent Doppler's unit interval of time domain,
Figure 709648DEST_PATH_IMAGE044
, PRF is radar pulse repetition frequency; Through matrix
Figure 526295DEST_PATH_IMAGE045
After the dimensionality reduction conversion, receive signal and target guiding vector and be respectively:
Figure 849009DEST_PATH_IMAGE047
（9）
Behind the dimensionality reduction, receive signal
Figure 341170DEST_PATH_IMAGE048
With the target guiding vector
Figure 645112DEST_PATH_IMAGE049
Degree of freedom be 15, therefore to the
Figure 920236DEST_PATH_IMAGE011
Individual range unit adopts the assorted covariance matrix of making an uproar after contiguous 50 range unit samples are estimated dimensionality reduction
Figure 626023DEST_PATH_IMAGE050
, that is:
Figure 961190DEST_PATH_IMAGE051
（10）
Wherein
Figure 755358DEST_PATH_IMAGE052
Be
Figure 896489DEST_PATH_IMAGE003
The output signal of individual range unit after the dimensionality reduction matrixing,
Figure 394467DEST_PATH_IMAGE053
For detecting the protected location number of range unit both sides, desirable
Figure 228431DEST_PATH_IMAGE054
The self-adaptation power of orientation-pulse domain two dimension dimensionality reduction STAP is:
Figure 241386DEST_PATH_IMAGE055
（11）
In the formula
Figure 186208DEST_PATH_IMAGE056
Be the normalization complex constant;
By formula (8), (9), (11), to
Figure 538692DEST_PATH_IMAGE011
Individual range unit orientation-pulse domain realizes the adaptive weight that long-range clutter suppresses
Figure 543557DEST_PATH_IMAGE057
For:
Figure 43809DEST_PATH_IMAGE058
（12）
Adopt above-mentioned flow process to calculate corresponding adaptive weight respectively to different range units, just can effectively suppress to influence the long-range clutter of target detection.
To sum up, the efficient three-dimensional space-time adaptive clutter suppression method of the present invention's proposition is made of pitching ADBF cascade orientation-pulse domain two dimension dimensionality reduction STAP; That is: to
Figure 464426DEST_PATH_IMAGE011
Individual range unit at first pitching dimension calculates by the row self-adaptation
Figure 64559DEST_PATH_IMAGE029
Suppress the short range clutter, then in the orientation-calculating of pulse domain self-adaptation
Figure 240325DEST_PATH_IMAGE057
Suppress long-range clutter.
Beneficial effect of the present invention is as follows:
The present invention is directed to the airborne non-positive side battle array Digital Array Radar that is operated under the middle high pulse repetition frequency and proposed a kind of efficient three-dimensional space-time adaptive clutter suppression method, comprise pitching robust ADBF, orientation-two steps of pulse domain two dimension dimensionality reduction STAP; With the two-dimentional dimensionality reduction STAP algorithm ratio of list of references [6] based on the local Combined Treatment, operand of the present invention is suitable substantially, but the pitching robust ADBF step that increases is the strong distance of filtering short range clutter heterogeneous effectively in advance, and then improves the rejection of follow-up orientation-residual spur of pulse domain two dimension STAP; Emulation experiment shows that operand of the present invention has only increased by 2.5% than the two-dimentional dimensionality reduction STAP algorithm based on the local Combined Treatment, but has improved about 15.5dB at the strong short range clutter district of the 400th range unit correspondence clutter rejection; Therefore the present invention can effectively improve the clutter inhibition of China's airborne non-positive side battle array Digital Array Radar, and fast convergence rate, and operand is little, is easy to engineering construction.
Description of drawings
Fig. 1 is the non-positive side array antenna geometric relationship figure of airborne radar.
Fig. 2 is forward sight battle array angle-Doppler's spectrogram.
Fig. 3 is forward sight battle array distance-Doppler spectrogram.
Fig. 4 (a) is the 400th range unit section two dimension dimensionality reduction STAP improvement factor synoptic diagram; Fig. 4 (b) is the 600th range unit section two dimension dimensionality reduction STAP improvement factor synoptic diagram.
Fig. 5 is pitching spatial domain dimensionality reduction exemplary plot.
Fig. 6 is that three-dimensional space-time adaptive clutter suppresses signal processing flow figure.
Fig. 7 is pitching and beam pattern figure.
Fig. 8 (a) is that the 400th the efficient three-dimensional space-time adaptive clutter of range unit section suppresses the improvement factor synoptic diagram; Fig. 8 (b) is that the 600th the efficient three-dimensional space-time adaptive clutter of range unit section suppresses the improvement factor synoptic diagram.
Embodiment
Below in conjunction with accompanying drawing the invention is described in further details.
The non-positive side of airborne radar battle array geometric relationship as shown in Figure 1.Suppose that antenna is the rectangle battle array, the rows of array elements column pitch is half-wavelength.If the normal direction of antenna plane is the x axle, the horizontal of antenna plane is the y axle, and the direction that makes progress perpendicular to ground is the z axle.The angle that the heading of carrier aircraft departs from y axle left side is The angle of radar main lobe beam position and array antenna transverse axis is
Figure 514498DEST_PATH_IMAGE060
The ground clutter of different distance is different to the angle of pitch of array antenna.The angle of pitch is defined as the angle of radar ray and antenna array normal, as shown in Figure 1
Figure 638312DEST_PATH_IMAGE061
,
Figure 922662DEST_PATH_IMAGE062
The pass of clutter Doppler frequency and distance is,
Figure 397506DEST_PATH_IMAGE063
（1）
In the formula
Figure 550139DEST_PATH_IMAGE064
Be carrier aircraft speed,
Figure 808687DEST_PATH_IMAGE065
Be the carrier aircraft height,
Figure 60676DEST_PATH_IMAGE066
Be radar wavelength,
Figure 22816DEST_PATH_IMAGE067
Be the clutter oblique distance,
Figure 588927DEST_PATH_IMAGE068
Be the clutter angle horizontal with respect to antenna plane,
Figure 15229DEST_PATH_IMAGE059
Crab angle for the relative y axle of carrier aircraft heading left side.For non-positive side battle array, the clutter two-dimensional spectrum of the fuzzy distance of each time is elliptic curve, and the clutter two-dimensional spectrum of different distance does not overlap.Angle-doppler spectral shown in Figure 2 is the clutter two dimension spectral line of the different fuzzy distances of same range unit under the forward sight battle array, and wherein the short range clutter of fuzzy distance and the clutter spectrum of the fuzzy distance of other each times of far field significantly depart from for the first time.Fig. 3 has provided the distance-doppler spectral of radar reception clutter, and wherein leftmost curve namely corresponds to the short range clutter of different distance unit.As seen, clutter slowly changes with distance at a distance, and distance variation more far away is more slow.
Conventional two-dimensional STAP is synthetic in pitching dimension with antenna array, forms the orientation linear array, then to the orientation-impulse sampling Design of Signal clutter self-adapting suppression filter of each range unit.Consider speed of convergence and the operand of adaptive algorithm, list of references [6] has proposed local Combined Treatment algorithm (JDL), and the guiding vector is respectively got in beam position both sides, spatial domain during namely according to the detection target empty
Figure 438120DEST_PATH_IMAGE069
Individual adjacent beams detects both sides, Doppler unit in time domain and respectively gets
Figure 887556DEST_PATH_IMAGE070
Individual adjacent Doppler unit, and then
Figure 257357DEST_PATH_IMAGE071
Local design dimensionality reduction STAP processor.Under simulation parameter, get
Figure 210270DEST_PATH_IMAGE072
,
Figure 72571DEST_PATH_IMAGE073
, utilize adjacent 50 range units to estimate two-dimentional dimensionality reduction STAP weights.Adopt improvement factor (IF) to weigh the clutter rejection, it is defined as
Figure 743724DEST_PATH_IMAGE074
（2）
In the formula
Figure 917216DEST_PATH_IMAGE075
Be the output letter miscellaneous noise ratio after the two-dimentional dimensionality reduction STAP processing,
Figure 990215DEST_PATH_IMAGE076
Receive the input letter miscellaneous noise ratio of signal for the array element individual pulse.Fig. 4 has provided the IF curve in the 400th range unit section and the 600th range unit section respectively.As seen, short range clutter power is big, seriously non-homogeneous down because airborne non-positive side battle array distance by radar blurs the 400th range unit section, the remarkable broadening of IF recess, and target detection performance descends at a slow speed; Because distance increase short range clutter power descends, and non-homogeneous intensity weakens, IF recess broadening slightly improves the 600th distance segment unit section.
The present invention is the improvement to conventional two-dimensional STAP algorithm, by increasing digital array spatial domain pitching dimension degree of freedom, in pitching-orientation-pulse domain studies three-dimensional space-time adaptive clutter suppression method, adopt cascade system, at first Adaptive Suppression is strong apart from short range clutter heterogeneous, and then designs two-dimentional dimensionality reduction STAP processor in orientation-pulse domain and suppress long-range clutter.
This paper scheme comprises following 2 steps: (1) pitching robust ADBF; (2) orientation-pulse domain two dimension dimensionality reduction STAP.Provide concrete analysis below.
(1) pitching robust ADBF
Can calculate the weight coefficient of each range unit pitching airspace filter device under the airborne numeral battle array radar ideal conditions according to radar parameter respectively.But because carrier aircraft height, the array element width of cloth equate error effect, the pitching wave filter recess of calculating can with short range clutter mismatch.List of references [14] has further been studied the complete space domain self-adapted algorithm of pitching on this basis, and its degree of freedom in system is pitching bay number.But under the range ambiguity, except short range clutter information, also comprised the far field clutter that receives in the training sample of choosing.Clutter its pitching dimension spatial domain guiding vector in far field is consistent with the pitching main beam direction, and this moment, the far field clutter was echo signal.Therefore, estimate that by the training sample that comprises the far field noise signal pitching ADBF weights that suppress the short range clutter can cause signal cancellation, the adaptive weight that calculates when serious can make the pitching main beam produce distortion.
The radar simulation parameter is as shown in table 1, and Fig. 4 has provided spatial domain and the wave beam corresponding figures of pitching routine and wave beam and different distance pitching adaptive weight.As seen, adaptive weight has all formed deep notch at short range clutter place to be suppressed it, but the pitching main beam has also produced distortion, the target signal to noise ratio loss that this must cause follow-up two-dimentional STAP to handle.
Table 1 radar system parameter
Parameter name Parameter values
Pulse repetition rate 5000Hz
The sampling bandwidth 5MHz
Line number 8
Columns 64
Carrier aircraft speed 150m/s
Pulse number in the CPI 64
The carrier aircraft flying height 8000m
Array element distance and wavelength ratio
Figure 20487DEST_PATH_IMAGE077
Several received pulses do not return as yet owing to far field signal (target and long-range clutter) before the Airborne PD Radar, generally as the zero filling pulse.If with the zero filling pulse as training sample estimation self-adaptive weights, the signal cancellation that can effectively avoid the far field clutter to cause then.PRF is in the radar system parameter During 5000Hz, the maximum distance of the 1st zero filling pulse correspondence is 30 kilometers, also has 20 ° deviation with the pitching main beam.Therefore, utilize first zero filling pulse to estimate short range clutter self-adapting weights, just can when suppressing the short range clutter, improve pitching main beam conformality.
Adopt first zero filling pulse to estimate short range clutter information, compare with list of references [14], number of training has just significantly reduced.The short range clutter requires to adopt too much range unit to come the estimate covariance matrix with the fast-changing non-stationary property of distance.Therefore, we further study pitching spatial domain dimensionality reduction scheme, to improve speed of convergence, reduce number of training.
Adopt first zero filling pulse estimate covariance matrix, each range unit has only short range clutter information.For contiguous some range units, its short range clutter spatial domain angle of pitch is approximate identical, and with the pitching main beam bigger differential seat angle is arranged, and therefore only needs two degree of freedom just can suppress it well.To phased array antenna, we handle respectively by row.To each row pitching linear array, dimensionality reduction forms pitching and wave beam and trim wave beam respectively as shown in Figure 5, namely
Figure 890540DEST_PATH_IMAGE001
（3）
Wherein M is antenna pitching dimension element number of array;
Figure 818045DEST_PATH_IMAGE002
Be listed as for every
Figure 19219DEST_PATH_IMAGE079
Individual array element is in the output signal of each pulse. Be pitching and the wave beam that forms behind M array element dimensionality reduction of every row,
Figure 914680DEST_PATH_IMAGE005
Be the trim wave beam that forms behind M array element dimensionality reduction of every row.In fact
Figure 965200DEST_PATH_IMAGE004
,
Figure 274959DEST_PATH_IMAGE005
Correspond to pitching every alignment battle array output signal and dimensionality reduction matrix
Figure 408000DEST_PATH_IMAGE006
Product, that is:
Figure 523723DEST_PATH_IMAGE007
（4）
Wherein
Figure 363503DEST_PATH_IMAGE008
Be the output of each array element of the every row of pitching in each pulse,
Figure 906480DEST_PATH_IMAGE009
Be the conjugate transpose operational symbol, Be transpose operator.After the dimension-reduction treatment of pitching spatial domain, adopt first zero filling pulse as training sample in order to the estimation self-adaptive weights, each array element of digital array received in signal and only comprised short range clutter information this moment.Suppose
Figure 118336DEST_PATH_IMAGE011
First zero filling reception of impulse signal of the every row pitching of individual range unit antenna array linear array is
Figure 874939DEST_PATH_IMAGE012
Signal, then every array unit
Figure 526500DEST_PATH_IMAGE011
Individual range unit first zero filling pulse pitching and wave beam
Figure 634134DEST_PATH_IMAGE013
Output signal and trim wave beam
Figure 91660DEST_PATH_IMAGE014
Output signal is:
Figure 640453DEST_PATH_IMAGE015
（5）
According to secondary lobe (SLC) principle that disappears mutually, to
Figure 793741DEST_PATH_IMAGE011
Individual range unit offsets short range clutter in pitching and the wave beam, adaptive weight with the trim wave beam
Figure 123091DEST_PATH_IMAGE016
Satisfy:
Figure 321991DEST_PATH_IMAGE017
（6）
Wherein
Figure 787608DEST_PATH_IMAGE018
Be the mathematical expectation operational symbol, Be the conjugate transpose operational symbol, then
Figure 863197DEST_PATH_IMAGE020
（7）
Autocorrelation matrix wherein Adopt after every row dimension-reduction treatment Individual range unit and vicinity
Figure 678072DEST_PATH_IMAGE022
The trim beamformer output signal of first zero filling pulse of individual range unit is estimated to obtain, and generally gets , namely
Figure 257138DEST_PATH_IMAGE024
,
Figure 962926DEST_PATH_IMAGE025
Be
Figure 351620DEST_PATH_IMAGE003
Individual range unit trim beamformer output signal.Cross-correlation matrix
Figure 142859DEST_PATH_IMAGE026
Adopt after every row dimension-reduction treatment
Figure 221673DEST_PATH_IMAGE011
Individual range unit and vicinity
Figure 781968DEST_PATH_IMAGE022
First zero filling pulse pitching of individual range unit and beamformer output signal and trim beamformer output signal are estimated to obtain, and generally get , namely
Figure 566570DEST_PATH_IMAGE027
,
Figure 573709DEST_PATH_IMAGE028
Be Individual range unit pitching and beamformer output signal.
By formula (4), (5), (6), (7), to
Figure 931058DEST_PATH_IMAGE011
Each row pitching linear array of individual range unit realizes the adaptive weight that the short range clutter suppresses
Figure 431310DEST_PATH_IMAGE029
For:
Figure 851927DEST_PATH_IMAGE030
（8）
Therefore, to Each row pitching linear array of individual range unit adopts adaptive weight
Figure 362247DEST_PATH_IMAGE029
Pulse by Pulse carries out filtering, just can effectively suppress the short range clutter in each pulse signal, avoids the far field echo signal to incur loss simultaneously.Can adopt above-mentioned flow process to calculate corresponding adaptive weight respectively to different range units, and then the first Pulse by Pulse of every array be carried out pitching filtering just can realize the short range clutter in each pulse signal of all range units is suppressed.
2) orientation-pulse domain two dimension dimensionality reduction STAP
After the short range clutter pitching filtering, the antenna array equivalence has formed an orientation linear array.Since elimination the short range clutter component, the clutter rejection of follow-up azimuth-range territory two dimension STAP will significantly improve.
The assumed position linear array has N array element, and each array element receives K pulse signal.Suppose Individual range unit receives signal
Figure 574102DEST_PATH_IMAGE031
, the guiding vector is when detecting target empty
Figure 963495DEST_PATH_IMAGE080
（9）
In the formula
Figure 310163DEST_PATH_IMAGE081
,
Figure 847324DEST_PATH_IMAGE082
,
Figure 674990DEST_PATH_IMAGE035
For the array element orientation to spacing,
Figure 856573DEST_PATH_IMAGE036
Be radar wavelength,
Figure 374142DEST_PATH_IMAGE037
Be the target Doppler frequency, Be the target azimuth incident angle,
Figure 902392DEST_PATH_IMAGE039
Be radar pulse recurrence interval (PRI),
Figure 798DEST_PATH_IMAGE040
For Kronecker amasss.
Further suppress to remain long-range clutter based on JDL algorithm design orientation-pulse domain two dimension dimensionality reduction STAP processor, the spatial domain, orientation is pointed to both sides at the detection object beam and is respectively got 2 adjacent beams, time domain is respectively got 1 adjacent Doppler unit in detection both sides, Doppler unit, and then the dimensionality reduction matrix is
Figure 689268DEST_PATH_IMAGE041
（10）
In the formula
Figure 76387DEST_PATH_IMAGE042
Be spatial domain adjacent beams interval,
Figure 242926DEST_PATH_IMAGE043
Be the adjacent Doppler's unit interval of time domain,
Figure 461418DEST_PATH_IMAGE083
, PRF is radar pulse repetition frequency; Through matrix
Figure 258473DEST_PATH_IMAGE045
After the dimensionality reduction conversion, receive signal and target guiding vector and be respectively
Figure 168365DEST_PATH_IMAGE047
（11）
Behind the dimensionality reduction, receive signal
Figure 179046DEST_PATH_IMAGE048
With the target guiding vector
Figure 209319DEST_PATH_IMAGE049
Degree of freedom be 15.Therefore to the
Figure 367768DEST_PATH_IMAGE011
Individual range unit adopts the assorted covariance matrix of making an uproar after contiguous 50 range unit samples are estimated dimensionality reduction
Figure 813792DEST_PATH_IMAGE050
, that is:
（12）
Wherein
Figure 208051DEST_PATH_IMAGE052
Be
Figure 791479DEST_PATH_IMAGE003
The output signal of individual range unit after the dimensionality reduction matrixing,
Figure 103511DEST_PATH_IMAGE053
For the protected location number of detecting unit both sides, desirable
Figure 151102DEST_PATH_IMAGE054
The self-adaptation power of orientation-pulse domain two dimension dimensionality reduction STAP processor is:
Figure 460860DEST_PATH_IMAGE055
（13）
In the formula
Figure 593901DEST_PATH_IMAGE056
Be the normalization complex constant.
By formula (10), (11), (13), to
Figure 724273DEST_PATH_IMAGE011
Individual range unit orientation-pulse domain realizes the adaptive weight that long-range clutter suppresses For:
Figure 107030DEST_PATH_IMAGE058
（14）
To sum up, the efficient three-dimensional space-time adaptive clutter suppression method of the present invention's proposition is made of pitching ADBF cascade orientation-pulse domain two dimension dimensionality reduction STAP; That is: to
Figure 461788DEST_PATH_IMAGE011
Individual range unit at first pitching dimension calculates by the row self-adaptation
Figure 318886DEST_PATH_IMAGE029
Suppress the short range clutter, then in the orientation-calculating of pulse domain self-adaptation
Figure 75489DEST_PATH_IMAGE057
Suppress long-range clutter, its signal processing flow as shown in Figure 6.
Simulation parameter is with table 1, and Fig. 7 calculates the spatial domain response diagram of adaptive weight respectively at 12.3Km and 18.3Km for this paper scheme step (1) pitching robust ADBF.Adaptive weight has all formed deep notch at the angle of pitch of short range clutter correspondence, and the pitching main beam is consistent with conventional formation main beam performance.
Come the performance of the efficient three-dimensional space-time adaptive clutter suppression method of quantitative analysis below with improvement factor (IF).Fig. 8 has provided the IF curve of two-dimentional dimensionality reduction STAP and this paper method respectively.Short range clutter district, the more two-dimentional dimensionality reduction STAP of the IF of this paper method is significantly increased, and wherein the 400th range unit section, the short range noise intensity is big, changes soon, and performance improvement is especially obvious, and this paper algorithm has improved about 15.5dB; The 600th range unit section, the short range noise intensity dies down, compare with two-dimentional dimensionality reduction STAP, and, this paper algorithm has improved about 6.5dB.
The present invention compares with the two-dimentional dimensionality reduction STAP of routine, has increased cascade pitching robust ADBF flow process, and its self-adaptive processing dimension only is 2.In simulation calculation, the self-adaptive processing dimension of two dimension dimensionality reduction STAP is 15, compare two-dimentional dimensionality reduction STAP, the operand of efficient three-dimensional space-time adaptive clutter suppression method has only increased about 0.25%, and when two-dimentional dimensionality reduction STAP self-adaptive processing dimension increased, it is littler that operand of the present invention increases ratio.Therefore, at in airborne non-positive side battle array Digital Array Radar under the high PRF system, the efficient three-dimensional space-time adaptive clutter suppression method that the present invention proposes is suitable with conventional two-dimensional dimensionality reduction STAP operand, and can significantly improve the clutter rejection of the fuzzy distance of each time.
List of references
[1] Brennan L E, Mallet J D, Reed I S. Theory of adaptive radar [J]. IEEE Trans. On AES., 1973, 9(2): 237-251。
[2] Wang Yongliang, Peng Yingning. the space-time adaptive signal is handled [M]. publishing house of Tsing-Hua University. 2000.
[3] Reed I S, Mallett J D, and Brennan L E. Rapid convergence rate in adaptive arrays[J]. IEEE Trans. AES. 1974, 10(6): 853-863。
[4] Meng Xiangdong. space-time two-dimensional Adaptive Signal Processing and moving-target detect [D]. Xi'an: Xian Electronics Science and Technology University, 2009.
[5] Shen Mingwei. handle moving target detection technique research [D] during airborne radar space. Nanjing: Nanjing Aero-Space University, 2008.
[6] Wang H, Cai L J. On Adaptive Spatial-Temporal Processing for Airborne Surveillance Radar System [J]. IEEE Trans. on AES. 1994, 30(3):660:670。
[7] Kreyenkamp O, Klemm R. Doppler compensation in forward-looking STAP radar [J]. IEE Proceedings-Radar Sonar and Naviqation. 2001,148(5): 253-258。
[8] Braham Himed, Zhang Yuhong, Abdelhak Hajjari. STAP with Angle-Doppler Compensation for Bistatic Airborne Radars [C]. Proceedings of the IEEE Radar Conference. 2002:311–317。
[9] Zatman M. Circular array STAP [J]. IEEE Trans. on AES, 2000, 36(2): 510~517。
[10] Braham H., Michel J.H., Zhang Y.H.. Bistatic STAP performance analysis in radar applications. Atlanta, GA, USA: Proceedings of the IEEE National Radar Conference, 2001: 198~203。
[11] Lapierre, F D, Verly J G, Van Droogenbroeck. New solutions to the problem of range dependence in bistatic STAP radars[C]. Piscataway NJ, USA: Proceedings of the IEEE International Radar Conference, 2003: 452-459。
[12] Xie W C, Wang Y L. New solution to Range-Dependence Problem in STAP Radar with HPRF[C]. Boston, MA: IEEE 2007 Radar Conference, 2007: 562-567。
[13] Colone F. Spectral Slope-based approach for mitigating bistatic Space-Time Adaptive Processing clutter dispersion [J]. IET Radar Sonar & Navigation. 2011, 5(5):593-603。
[14] Meng Xiangdong, Wang Tong etc. the pitching that airborne phased array radar short range clutter suppresses is to space domain self-adapted algorithm [J]. electronics and information journal, 2010,32 (4): 948-952.

Claims (1)

1. the efficient three-dimensional space-time adaptive clutter suppression method of airborne Digital Array Radar is characterized in that comprising the steps:
1) pitching robust ADBF adaptive beam forms: to each row pitching linear array, dimensionality reduction forms pitching and wave beam and trim wave beam respectively, namely
Figure 201210195856X100001DEST_PATH_IMAGE001
（1）
M is antenna pitching dimension element number of array in the formula, Be listed as for every
Figure 155000DEST_PATH_IMAGE004
Individual array element is in the output signal of each pulse,
Figure 201210195856X100001DEST_PATH_IMAGE005
Be pitching and the wave beam that forms behind M array element dimensionality reduction of every row,
Figure 552483DEST_PATH_IMAGE006
Be the trim wave beam that forms behind M array element dimensionality reduction of every row; Pitching and wave beam
Figure 343722DEST_PATH_IMAGE005
With the trim wave beam
Figure 422536DEST_PATH_IMAGE006
Correspond to pitching every alignment battle array output signal and dimensionality reduction matrix
Figure DEST_PATH_IMAGE007
Product, that is:
Figure 982830DEST_PATH_IMAGE008
（2）
In the formula Be the output signal of each array element of the every row of pitching in each pulse;
Figure 879111DEST_PATH_IMAGE010
Be the conjugate transpose operational symbol;
Figure DEST_PATH_IMAGE011
Be transpose operator; After the dimension-reduction treatment of pitching spatial domain, adopt first zero filling pulse as training sample in order to the estimation self-adaptive weights, each array element of digital array received in signal and only comprised short range clutter information this moment; When definition the
Figure DEST_PATH_IMAGE013
First zero filling reception of impulse signal of the every row pitching of individual range unit antenna array linear array is
Figure 892067DEST_PATH_IMAGE014
, every array unit then
Figure 839818DEST_PATH_IMAGE013
Individual range unit first zero filling pulse pitching and wave beam
Figure DEST_PATH_IMAGE015
Output signal and trim wave beam
Figure 316936DEST_PATH_IMAGE016
Output signal is:
Figure DEST_PATH_IMAGE017
（3）
According to the secondary lobe SLC principle that disappears mutually, to
Figure 321801DEST_PATH_IMAGE013
Individual range unit offsets short range clutter in pitching and the wave beam, adaptive weight with the trim wave beam Satisfy:
Figure DEST_PATH_IMAGE019
（4）
Wherein Be the mathematical expectation operational symbol,
Figure DEST_PATH_IMAGE021
Be the conjugate transpose operational symbol, then
Figure 839873DEST_PATH_IMAGE022
（5）
Autocorrelation matrix in the formula
Figure DEST_PATH_IMAGE023
Adopt after every row dimension-reduction treatment Individual range unit and vicinity
Figure 6117DEST_PATH_IMAGE024
The trim beamformer output signal of first zero filling pulse of individual range unit is estimated to obtain, and gets , namely
Figure 292742DEST_PATH_IMAGE026
, Be
Figure 416556DEST_PATH_IMAGE004
Individual range unit trim beamformer output signal; Cross-correlation matrix Adopt after every row dimension-reduction treatment
Figure 175750DEST_PATH_IMAGE013
Individual range unit and vicinity
Figure 266066DEST_PATH_IMAGE024
First zero filling pulse pitching of individual range unit and beamformer output signal and trim beamformer output signal are estimated to obtain, and get
Figure 447648DEST_PATH_IMAGE025
, namely , Be
Figure 661778DEST_PATH_IMAGE004
Individual range unit pitching and beamformer output signal;
By formula (2), (3), (4), (5), to
Figure 290205DEST_PATH_IMAGE013
Each row pitching linear array of individual range unit realizes the adaptive weight that the short range clutter suppresses
Figure DEST_PATH_IMAGE031
For:
（6）
2) orientation-pulse domain two dimension dimensionality reduction STAP space-time adaptive is handled: after pitching robust ADBF suppressed the short range clutter, antenna array was equivalent to the orientation linear array; Linear array has N array element when the definition orientation, and each array element receives K pulse signal; When definition the
Figure DEST_PATH_IMAGE033
Individual range unit receives signal , guiding vector when detecting target empty
（7）
In the formula
Figure 591764DEST_PATH_IMAGE036
,
Figure DEST_PATH_IMAGE037
,
Figure 23883DEST_PATH_IMAGE038
For the array element orientation to spacing,
Figure DEST_PATH_IMAGE039
Be radar wavelength,
Figure 976795DEST_PATH_IMAGE040
Be the target Doppler frequency, Be the target azimuth incident angle,
Figure 836167DEST_PATH_IMAGE042
Be radar pulse recurrence interval PRI,
Figure DEST_PATH_IMAGE043
For Kronecker amasss;
Further suppress to remain long-range clutter based on local Combined Treatment algorithm design orientation-pulse domain two dimension dimensionality reduction STAP processor, the spatial domain, orientation is pointed to both sides at the detection object beam and is respectively got 2 adjacent beams, time domain is respectively got 1 adjacent Doppler unit in detection both sides, Doppler unit, and then the dimensionality reduction matrix is
Figure 507320DEST_PATH_IMAGE044
（8）
In the formula
Figure DEST_PATH_IMAGE045
Be spatial domain adjacent beams interval, Be the adjacent Doppler's unit interval of time domain,
Figure DEST_PATH_IMAGE047
, PRF is radar pulse repetition frequency; Through matrix
Figure 869655DEST_PATH_IMAGE048
After the dimensionality reduction conversion, receive signal and target guiding vector and be respectively:
Figure DEST_PATH_IMAGE049
Figure 899928DEST_PATH_IMAGE050
（9）
Behind the dimensionality reduction, receive signal
Figure DEST_PATH_IMAGE051
With the target guiding vector
Figure 58377DEST_PATH_IMAGE052
Degree of freedom be 15, therefore to the Individual range unit adopts the assorted covariance matrix of making an uproar after contiguous 50 range unit samples are estimated dimensionality reduction
Figure DEST_PATH_IMAGE053
, that is:
Figure 759802DEST_PATH_IMAGE054
（10）
Wherein
Figure DEST_PATH_IMAGE055
Be
Figure 960976DEST_PATH_IMAGE004
The output signal of individual range unit after the dimensionality reduction matrixing,
Figure 671968DEST_PATH_IMAGE056
For detecting the protected location number of range unit both sides, get
Figure DEST_PATH_IMAGE057
The adaptive weight of orientation-pulse domain two dimension dimensionality reduction STAP is:
Figure 984001DEST_PATH_IMAGE058
（11）
In the formula
Figure DEST_PATH_IMAGE059
Be the normalization complex constant;
By formula (8), (9), (11), to
Figure 93908DEST_PATH_IMAGE013
Individual range unit orientation-pulse domain realizes the adaptive weight that long-range clutter suppresses
Figure 465984DEST_PATH_IMAGE060
For:
Figure DEST_PATH_IMAGE061
（12）
Adopt above-mentioned flow process to calculate corresponding adaptive weight respectively to different range units, just can effectively suppress to influence the long-range clutter of target detection.
CN 201210195856 2012-06-14 2012-06-14 Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar Expired - Fee Related CN102721947B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210195856 CN102721947B (en) 2012-06-14 2012-06-14 Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210195856 CN102721947B (en) 2012-06-14 2012-06-14 Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar

Publications (2)

Publication Number Publication Date
CN102721947A CN102721947A (en) 2012-10-10
CN102721947B true CN102721947B (en) 2013-09-25

Family

ID=46947763

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210195856 Expired - Fee Related CN102721947B (en) 2012-06-14 2012-06-14 Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar

Country Status (1)

Country Link
CN (1) CN102721947B (en)

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103136165B (en) * 2013-01-29 2016-01-20 西安电子工程研究所 A kind of method of the Adaptive Sidelobe Canceling weights based on FPGA
CN103176168B (en) * 2013-02-05 2015-09-09 西安电子科技大学 A kind of airborne non-working side battle array radar short range clutter cancellation method
CN103414482B (en) * 2013-06-17 2015-08-19 电子科技大学 The closed loop sidelobe cancellation method and system that modulus combines
CN103383448B (en) * 2013-06-25 2015-06-17 西安电子科技大学 Clutter suppression method suitable for high pulse repetition frequency (HPRF) waveform airborne radar
CN103364764B (en) * 2013-06-25 2015-06-17 西安电子科技大学 Airborne radar non-stationary clutter suppression method
CN103383449B (en) * 2013-07-14 2015-04-08 西安电子科技大学 ESPRIT algorithm based short-range clutter suppression method for airborne radar
CN103605114B (en) * 2013-12-03 2015-07-15 西安电子科技大学 Non-broadside array airborne radar short range clutter suppression method based on multiple frequencies
CN103728607A (en) * 2014-01-16 2014-04-16 西安电子科技大学 Space time code three-dimensional self-adaptation clutter cancelling method for onboard multiple input multiple output (MIMO) radar
CN104076343B (en) * 2014-06-25 2017-02-15 西安电子科技大学 Satellite-borne three-channel SAR-GMTI self-adaptive clutter suppression method
CN104535972A (en) * 2015-01-30 2015-04-22 西安电子科技大学 Coherent transferring type jamming inhibition method for airborne radar
CN104635214B (en) * 2015-02-12 2017-01-04 西安电子科技大学 Air-borne Forward-looking frequency diversity array radar range ambiguity clutter suppression method
CN106291476B (en) * 2016-07-29 2019-03-29 西安电子科技大学 The Radar Clutter acquisition methods of airborne three-dimensional isomery battle array
CN107678002B (en) * 2017-08-23 2020-11-17 西安电子工程研究所 Method for efficiently measuring ground clutter improvement factor of MTD radar
CN107797111A (en) * 2017-09-28 2018-03-13 中国人民解放军国防科技大学 Robust multi-channel SAR signal reconstruction method under non-uniform scattering coefficient scene
CN109375174B (en) * 2018-12-17 2023-03-21 西安八阵图电子科技有限公司 Two-dimensional electrical scanning three-coordinate radar method based on inverted T-shaped linear array
CN111398907B (en) * 2019-01-02 2023-04-07 中国人民解放军空军预警学院 Phased array radar space domain two-dimensional digital difference beam forming method
CN110133603A (en) * 2019-06-27 2019-08-16 哈尔滨工业大学 High-frequency ground wave radar ocean clutter cancellation method based on rooting Euclidean geometry center of gravity
CN111781603B (en) * 2020-06-09 2023-12-19 南京航空航天大学 Ground clutter suppression method for airborne weather radar
CN114966650A (en) * 2021-02-22 2022-08-30 华为技术有限公司 Signal processing method, signal processing device, storage medium and vehicle
CN113406615B (en) * 2021-05-25 2022-10-04 森思泰克河北科技有限公司 Target tracking method and device of binary phase modulation array radar
CN113156392B (en) * 2021-05-26 2022-04-08 哈尔滨工业大学 Clutter suppression method based on pitching domain self-adaptive processing
CN113740824B (en) * 2021-09-08 2023-09-15 中国科学院空天信息创新研究院 Three-dimensional combined beam dimension-reducing space-time self-adaptive processing method based on beam priori

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101414002A (en) * 2008-12-01 2009-04-22 西安电子科技大学 Method for counteracting airborne radar non-self-adapting clutter
CN101561497A (en) * 2009-05-22 2009-10-21 西安电子科技大学 Airborne radar clutter suppression method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101414002A (en) * 2008-12-01 2009-04-22 西安电子科技大学 Method for counteracting airborne radar non-self-adapting clutter
CN101561497A (en) * 2009-05-22 2009-10-21 西安电子科技大学 Airborne radar clutter suppression method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
和洁 等.机载MIMO雷达降维空时自适应处理算法.《电子科技大学学报》.2012,第41卷(第1期),第31-35页.
孟祥东 等.机载相控阵雷达近程杂波抑制的俯仰向空域自适应算法.《电子与信息学报》.2010,第32卷(第4期),第948-952页.
机载MIMO雷达降维空时自适应处理算法;和洁 等;《电子科技大学学报》;20120131;第41卷(第1期);第31-35页 *
机载相控阵雷达近程杂波抑制的俯仰向空域自适应算法;孟祥东 等;《电子与信息学报》;20100430;第32卷(第4期);第948-952页 *

Also Published As

Publication number Publication date
CN102721947A (en) 2012-10-10

Similar Documents

Publication Publication Date Title
CN102721947B (en) Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar
CN102156279B (en) Method for detecting moving target on ground by utilizing bistatic radar based on MIMO (Multiple Input Multiple Output)
CN101858976B (en) Moving target detection method based on multiple sub-apertures of single-channel SAR
CN103399303B (en) Airborne radar resists intensive deceiving jamming method and system
CN102288948B (en) High-speed platform high-speed air moving target detection method based on STAP (Spacetime Adaptive Processing)
CN103018727A (en) Sample-training-based non-stationary clutter suppression method of vehicle-mounted radar
CN106342236B (en) Based on the airborne radar clutter suppression method of difference beam
CN105510887B (en) To the active suppressing formula disturbance restraining method of airborne radar under a kind of clutter background
CN111913157B (en) Sea clutter suppression method based on radar signal space-time decorrelation model
CN106569212A (en) Multichannel SAR-GMTI range ambiguity clutter suppression method
CN111220955B (en) Airborne weather radar ground clutter suppression method based on vertical array covariance matrix characteristic decomposition
CN102788970B (en) High-efficiency self-adapting suppression method for short-range clutters of airborne digital array radar
CN103744067A (en) Non-adaptive airborne non-side-looking radar short-range clutter suppression method
CN104345300A (en) Onboard non-positive side view array radar STAP (Space Time Adaptive Processing) method for clutter space-time spectrum linear compensation
CN109765529B (en) Millimeter wave radar anti-interference method and system based on digital beam forming
CN104345299A (en) Airborne MIMO (Multiple Input Multiple Output) radar space-time self-adaptive processing method based on simplified EC
CN110082744A (en) The MIMO airborne bistatic radar clutter suppression method of Doppler's stepped multiplexing
CN111781603B (en) Ground clutter suppression method for airborne weather radar
CN101907702A (en) Two-dimensional multi-pulse canceller for MIMO radar
CN107748364A (en) Low wind field speed estimation method based on contraction multistage wiener filter
CN110646765B (en) Riemann distance-based generalized sidelobe cancellation algorithm
CN103728596A (en) Method for inhibiting distributed interference of bi-static MIMO radar
CN108872947B (en) Sea clutter suppression method based on subspace technology
CN109143235A (en) A kind of biradical forward sight synthetic aperture radar Ground moving target detection method
CN109597034B (en) Space-time adaptive processing method based on Euclidean distance

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130925

Termination date: 20160614

CF01 Termination of patent right due to non-payment of annual fee