CN105223557B  Airborne early warning radar clutter suppression method based on accessory channel  Google Patents
Airborne early warning radar clutter suppression method based on accessory channel Download PDFInfo
 Publication number
 CN105223557B CN105223557B CN201510717941.1A CN201510717941A CN105223557B CN 105223557 B CN105223557 B CN 105223557B CN 201510717941 A CN201510717941 A CN 201510717941A CN 105223557 B CN105223557 B CN 105223557B
 Authority
 CN
 China
 Prior art keywords
 num
 early warning
 accessory channel
 warning radar
 dimensionality reduction
 Prior art date
Links
Abstract
Description
Technical field
The invention belongs to radar clutter suppression technology field, more particularly to a kind of airborne early warning radar based on accessory channel Clutter suppression method, suitable for solving the problems, such as the decline of airborne early warning radar clutter rejection nonhomogeneous clutter environment, change It is apt to its clutter recognition performance.
Background technology
Airborne early warning radar can flexibly, rapidly be deployed in required local and in widespread attention, its master with it It is that target is detected in clutter background to want task, and the target to detecting carries out locating and tracking, and clutter is effectively pressed down System is to improve the core of airborne early warning radar service behaviour.Therefore, before the target that locating and tracking detects, it is necessary first to suppress miscellaneous Clutter present in ripple background or interference.If when under airborne early warning radar depending on being radiated at relatively flat low scattering region, production Raw clutter or interference can be very weak, is handled using conventional method.But generated due to the motion of airborne early warning radar Clutter spectrum spreads in the main lobe broadening and sidelobe clutter of Doppler frequency domain so that shows very strong spacetime coupled characteristic, therefore Need to suppress caused clutter or interference using spacetime adaptive signal transacting (STAP) technology.Spacetime adaptive processing (STAP) technology can make full use of spatial information (si) and timedomain information, and can effective clutter reduction, but as a rule Enough independent same distributions (independent and identically distributed, IID) instruction can not almost be obtained Practice sample to estimate spacetime covariance matrix.Even if enough independent same distribution number of training are obtained, for high level matrix There is also the difficulty that amount of calculation and precision aspect are difficult to for the computing inverted.
In the eighties, German doctor R.Klemm is opened up to spacetime adaptive signal transacting (STAP) technology Property theoretical research, he, which passes through, carries out thoroughgoing and painstaking analysis to noise performance, finds the big characteristic value of spacetime covariance matrix Number be no more than N+M1, wherein N is the array number of airborne early warning radar, and M is airborne early warning radar in a Coherent processing The umber of pulse of transmitting, illustrates that the full spacetime adaptive signal transacting (STAP) for clutter reduction exists at dimensionality reduction really in interval The possibility of reason, on this basis, he proposes accessory channel method (Auxiliary Channel ReceiverACR), at dimensionality reduction Dimension after reason is by NM to N+M1.Research show this kind of dimensionreduction treatment in performance close to optimal full space time processing effect, but There is also following two problems in actual applications：First, dimensionreduction treatment is in performance close to optimal full space time processing effect Obtained in the case of without amplitude phase error, if it is considered that spatial domain error, clutter spectrum caused by airborne early warning radar can be along sky Domain Directional Extension so that its clutter dimension is significantly increased, and process performance is decreased obviously；Second, in the array element of airborne early warning radar When number N is bigger, required processor dimension is also bigger.Therefore, accessory channel method still need in actual applications into Onestep optimization or improvement, so that clutter recognition performance is improved.
The content of the invention
The problem of existing for above prior art, it is a kind of based on the airborne of accessory channel it is an object of the invention to propose Early warning radar clutter suppression method, this method can solve the problem that tolerance difference and array number of traditional accessory channel method to spatial domain error The problem of being difficult to application when more, and the utilization rate by improving clutter ridge carries out dimensionreduction treatment, so as to reduce required independence With distribution number of training, while recess can be also formed at clutter ridge, improve clutter recognition effect.
To reach abovementioned technical purpose, the present invention, which adopts the following technical scheme that, to be achieved.
A kind of airborne early warning radar clutter suppression method based on accessory channel, it is characterised in that comprise the following steps：
Step 1, the threedimensional echo data X that airborne early warning radar receives is obtained_{N×M×L}, and airborne early warning thunder is obtained accordingly Up to the threedimensional echo data X received_{N×M×L}The clutter ridge of formation and the main beam spatial domain frequency θ of airborne early warning radar emission_{s}, Then threedimensional echo data X airborne early warning radar received_{N×M×L}Rearranged in the way of row, obtain airborne early warning thunder Up to the twodimentional echo data X received_{NM×L}；Wherein, N is the element number of array of airborne early warning radar, and M is airborne early warning radar one The umber of pulse of transmitting in individual coherent processing inteval, L are the threedimensional echo data X that airborne early warning radar receives_{N×M×L}Range gate Number；
Step 2, the threedimensional echo data X received according to airborne early warning radar_{N×M×L}The clutter ridge of formation, is obtained respectively Num accessory channel and num search passage, and then respectively obtain by each selfcorresponding spatial domain frequency structure of num accessory channel Into spatial domain frequency vector θ, the temporal frequency vector that is made up of each selfcorresponding temporal frequency of num accessory channelAnd Each selfcorresponding time domain steering vector of num search passage, is then calculated dimensionality reduction matrix corresponding to num accessory channel T_{b}；
Step 3, according to the main beam spatial domain frequency θ of airborne early warning radar emission_{s}, airborne early warning radar emission is calculated Main beam spatial domain steering vector G_{s}, further according to each selfcorresponding time domain steering vector of num search passage, appoint to take kth and search Suo Tongdao time domain steering vector F_{k}, then according to the main beam spatial domain steering vector G of airborne early warning radar emission_{s}Searched with kth Suo Tongdao time domain steering vector F_{k}, search passage column vector S corresponding to kth of search passage is calculated_{k}, and then obtain num Search passage matrix S corresponding to individual search passage；Wherein, k ∈ { 1,2 ... num }；
Step 4, the dimensionality reduction square corresponding to search passage matrix S and num accessory channel according to corresponding to num search passage Battle array T_{b}, optimization dimensionality reduction matrix corresponding to num_s × num_t accessory channel is calculatedAnd then it is calculated based on auxiliary The transformation matrix T of passage；Wherein, k ∈ { 1,2 ... num }, num_s are the optimization dimensionality reduction matrixComprising spatial frequency Number, num_t are the optimization dimensionality reduction matrixComprising Doppler frequency number, num_s × num_t<<num；
Step 5, according to the transformation matrix T based on accessory channel, the twodimentional echo data received to airborne early warning radar X_{NM×L}With num search passage corresponding to search passage matrix S carry out dimensionreduction treatment respectively, respectively obtain the number of echoes after dimensionality reduction According to X_{T}With the search passage steering vector S after dimensionality reduction_{T}；
Step 6, according to the echo data X after dimensionality reduction_{T}With the search passage steering vector S after dimensionality reduction_{T}, dimensionality reduction is calculated Filtering weight vector afterwards；
Step 7, according to the filtering weight vector W after dimensionality reduction_{T}, to the echo data X after dimensionality reduction_{T}Clutter recognition processing is carried out, is obtained To range Doppler figure.
Compared with prior art, advantages of the present invention and improvement are：
First, the present invention utilizes clutter ridge information, different dimensionality reduction matrixes can be formed to different search passages, and lead to Cross the clutter that part accessory channel having identical Doppler frequency with search passage, being formed along clutter ridge offsets search passage Component, deficiency of the original auxiliary channel algorithm to spatial domain error tolerance difference is effectively improved, improves practical application.
Second, the present invention utilizes the dimensionality reduction transition matrix based on accessory channel method, the echo received to airborne early warning radar While data carry out doppler filtering processing, dimensionreduction treatment is also carried out, effectively improves original auxiliary channel algorithm due to can The problem of clutter covariance matrix estimation is inaccurate caused by independent same distribution number of training deficiency, and airborne early warning The free degree of radar it is big and caused by amount of calculation is excessive and the problem of equipment cost increase, so as to reduce estimate covariance square Sample number required for battle array so that the present invention can will not also reduce clutter recognition performance while number of training deficiency.
Brief description of the drawings
The present invention is described in further detail with reference to the accompanying drawings and detailed description.
Fig. 1 is a kind of schematic flow sheet of airborne early warning radar clutter suppression method based on accessory channel of the present invention；
Fig. 2 (a) be using the range Doppler figure obtained after pulse Doppler (PD) algorithm process,
Fig. 2 (b) is using the range Doppler figure obtained after the processing of original accessory channel method (ACR)；
Fig. 3 is the range Doppler figure obtained after being handled using the inventive method；
Fig. 4 (a) is the global twodimentional response diagram obtained using the inventive method,
Fig. 4 (b) is the partial enlargement two dimension response diagram obtained using the inventive method；
After Fig. 5 is handles using pulse Doppler (PD) algorithm, original accessory channel method (ACR) and the inventive method respectively Obtained clutter dump power comparison diagram.
Embodiment
Reference picture 1, the flow for a kind of airborne early warning radar clutter suppression method based on accessory channel of the present invention are shown It is intended to, airborne early warning radar clutter suppression method of this kind based on accessory channel, comprises the following steps：
Step 1, the threedimensional echo data X that airborne early warning radar receives is obtained_{N×M×L}, and airborne early warning thunder is obtained accordingly Up to the threedimensional echo data X received_{N×M×L}The clutter ridge of formation and the main beam spatial domain frequency θ of airborne early warning radar emission_{s}, Then threedimensional echo data X airborne early warning radar received_{N×M×L}Rearranged in the way of row, obtain airborne early warning thunder Up to the twodimentional echo data X received_{NM×L}；Wherein, N is the element number of array of airborne early warning radar, and M is airborne early warning radar one The umber of pulse of transmitting in individual coherent processing inteval, L are the threedimensional echo data X that airborne early warning radar receives_{N×M×L}Range gate Number.
Specifically, airborne early warning radar chooses positive side battle array airborne early warning radar, and positive side battle array airborne early warning radar antenna includes N Individual array element, N number of array element receive the threedimensional echo data X of ground scatter body reflection_{N×M×L}, as acquisition airborne early warning radar The threedimensional echo data X received_{N×M×L}；Wherein, N is the element number of array of airborne early warning radar, and M is airborne early warning radar one The umber of pulse of transmitting in individual coherent processing inteval, L are the threedimensional echo data X that airborne early warning radar receives_{N×M×L}Range gate Number.
Wherein, because airborne early warning radar uses positive side battle array airborne early warning radar, and the day of positive side battle array airborne early warning radar Linear array member axially it is consistent with the heading of carrier aircraft, the main beam direction of positive side battle array airborne early warning radar illumination ground scatter body and The cone angle cosine value that the bay of positive side battle array airborne early warning radar is axially formed, and the Doppler frequency of ground scatter body echo A kind of linear relationship, and the ground scatter body echo is exactly clutter, the clutter the cone angle cosine value with it is described how general Strangle in frequency space, clutter distribution is straight line, using this straight line as clutter ridge, and obtains airborne early warning radar emission Main beam spatial domain frequency θ_{s}For 0, while the threedimensional echo data X that airborne early warning radar is received_{N×M×L}The weight in the way of row New arrangement, obtains the twodimentional echo data X that airborne early warning radar receives_{NM×L}。
Step 2, the threedimensional echo data X received according to airborne early warning radar_{N×M×L}The clutter ridge of formation, is obtained respectively Num accessory channel and num search passage, and then respectively obtain by each selfcorresponding spatial domain frequency structure of num accessory channel Into spatial domain frequency vector θ, the temporal frequency vector that is made up of each selfcorresponding temporal frequency of num accessory channelAnd Each selfcorresponding time domain steering vector of num search passage, is then calculated dimensionality reduction matrix corresponding to num accessory channel T_{b}。
Specifically, threedimensional echo data X airborne early warning radar received_{N×M×L}Element number of array and airborne early warning thunder Up to the space for the umber of pulse composition launched in a coherent processing inteval, the threedimensional echo received as airborne early warning radar Data X_{N×M×L}Array elementpulse domain twodimensional space, it is and airborne pre according to the twodimensional space interior edge of the array elementpulse domain Alert radar receives threedimensional echo data X_{N×M×L}The clutter ridge of formation, obtains num accessory channel respectively and num search is logical Road, detailed process are：
If the threedimensional echo data X that airborne early warning radar receives_{N×M×L}The clutter ridge slope of formation is β, airborne early warning thunder Up to the threedimensional echo data X received_{N×M×L}Maximum of the clutter ridge of formation in the frequency of spatial domain is θ_{max}, airborne early warning radar The threedimensional echo data X received_{N×M×L}Minimum value of the clutter ridge of formation in the frequency of spatial domain is θ_{min}, airborne early warning radar connects The threedimensional echo data X received_{N×M×L}Maximum of the clutter ridge of formation in Doppler frequency beAirborne early warning radar The threedimensional echo data X received_{N×M×L}Minimum value of the clutter ridge of formation in Doppler frequency beIts expression formula point It is not：
Wherein, N is the element number of array of airborne early warning radar, and M is that airborne early warning radar is sent out in a coherent processing inteval The umber of pulse penetrated, L are the threedimensional echo data X that airborne early warning radar receives_{N×M×L}Range gate number, V be carrier aircraft flight speed Degree, λ are wavelength, f_{r}For pulse recurrence frequency, d represents the adjacent array element interval of airborne early warning radar.
Then, the threedimensional echo data X received to airborne early warning radar_{N×M×L}The clutter ridge spatial domain frequency and machine of formation Carry the threedimensional echo data X that early warning radar receives_{N×M×L}The clutter ridge Doppler frequency of formation is evenly dividing respectively, is obtained To num accessory channel and num search passage, and each selfcorresponding time domain steering vector of num search passage.
If the spatial domain frequency of ith of accessory channel is θ_{i}, then the temporal frequency of ith of accessory channel beIts Expression formula is respectively：
The spatial domain frequency vector θ being then made up of each selfcorresponding spatial domain frequency of num accessory channel, and it is auxiliary by num The temporal frequency vector for helping each selfcorresponding temporal frequency of passage to formExpression formula be respectively：
θ=[θ_{1},θ_{2},…θ_{num}]
Further according to the spatial domain frequency vector θ being made up of each selfcorresponding spatial domain frequency of num accessory channel, and by num The temporal frequency vector that individual each selfcorresponding temporal frequency of accessory channel is formedIt is calculated respectively by num accessory channel The spatial domain steering vector column vector G that each selfcorresponding spatial domain steering vector is formed_{b}With by num accessory channel it is each selfcorresponding when The time domain steering vector column vector F that domain steering vector is formed_{b}, its expression is respectively：
G_{b}=[1；e^{j2πθ}；…；e^{j(N1)2πθ}]
And then dimensionality reduction matrix T corresponding to num accessory channel is calculated_{b}, its expression formula is：
Wherein, N be airborne early warning radar element number of array, G_{b}To be oriented to by each selfcorresponding spatial domain of num accessory channel The spatial domain steering vector column vector that vector is formed, F_{b}It is made up of each selfcorresponding time domain steering vector of num accessory channel Time domain steering vector column vector, T_{b}For dimensionality reduction matrix corresponding to num accessory channel,For Kronecker product operation symbol.
Step 3, according to the main beam spatial domain frequency θ of airborne early warning radar emission_{s}, airborne early warning radar emission is calculated Main beam spatial domain steering vector G_{s}, further according to each selfcorresponding time domain steering vector of num search passage, appoint to take kth and search Suo Tongdao time domain steering vector F_{k}, then according to the main beam spatial domain steering vector G of airborne early warning radar emission_{s}Searched with kth Suo Tongdao time domain steering vector F_{k}, search passage column vector S corresponding to kth of search passage is calculated_{k}, and then obtain num Search passage matrix S corresponding to individual search passage；Wherein, k ∈ { 1,2 ... num }.
Specifically, search passage column vector S corresponding to kth of search passage in num search passage_{k}It is according to airborne pre The main beam spatial domain steering vector G of alert radar emission_{s}With the time domain steering vector F of kth of search passage_{k}It is calculated, it is expressed Formula is：
Wherein, G_{s}For the main beam spatial domain steering vector of airborne early warning radar emission, and num search passage is each right The main beam spatial domain steering vector for the airborne early warning radar emission answered is the same；F_{k}It is oriented to for the time domain of kth of search passage Vector, S_{k}For search passage column vector corresponding to kth of search passage,For Kronecker product operation symbol.
The main beam spatial domain steering vector G of airborne early warning radar emission_{s}With the time domain steering vector F of kth of search passage_{k} Expression formula be respectively：
θ_{s}=0
Wherein, θ_{s}For the main beam spatial domain frequency of airborne early warning radar emission,For the time domain frequency of kth of search passage Rate, N are the element number of array of airborne early warning radar, and M is the umber of pulse that airborne early warning radar is launched in a coherent processing inteval.
According to search passage column vector S corresponding to kth of search passage_{k}, it is logical to obtain search corresponding to num search passage Road matrix S.
Step 4, the dimensionality reduction square corresponding to search passage matrix S and num accessory channel according to corresponding to num search passage Battle array T_{b}, optimization dimensionality reduction matrix corresponding to num_s × num_t accessory channel is calculatedAnd then it is calculated based on auxiliary The transformation matrix T of passage；Wherein, k ∈ { 1,2 ... num }, num_s are the optimization dimensionality reduction matrixComprising spatial frequency Number, num_t are the optimization dimensionality reduction matrixComprising Doppler frequency number, num_s ∈ { 1,2 ... num } num_t ∈ 1, 2 ... num }, num_s × num_t<<num.
Specifically, if the Doppler frequency of (i ∈ { 1,2 ... num }) accessory channel is ith in num accessory channelThen choose num accessory channel in (i (num_t1)/2) it is individual to (i+ (num_t1)/2) individual accessory channel each Corresponding Doppler frequency, as the Doppler frequency of optimization accessory channel, choose (i (num_s in num accessory channel 1) it is/2) individual to (i+ (num_s1)/2) individual each selfcorresponding spatial frequency of accessory channel, the space as optimization accessory channel Frequency, then optimize accessory channel and include num_s spatial frequency and num_t Doppler frequency, and obtain optimization auxiliary accordingly The spatial domain frequency vector θ of passage_{s'}With the temporal frequency vector of optimization accessory channelIts expression formula is respectively：
θ_{s'}=[θ_{i(num_s1)/2},θ_{i[(num_s1)/2]+1},…θ_{i},…,θ_{i+[(num_s1)/2]1},θ_{i+(num_s1)/2}]
Typically, num_s, num_t and num relation meet num_s × num_t<<Num, and num_s and num_t difference Take odd number value.
According to the spatial domain frequency vector θ of optimization accessory channel_{s'}With the temporal frequency vector of optimization accessory channelRespectively The spatial domain steering vector G of optimization accessory channel is calculated_{b'}With the time domain steering vector F of optimization accessory channel_{b'}, its expression formula Respectively：
And then optimization dimensionality reduction matrix corresponding to num_s × num_t accessory channel is calculatedIts expression formula is：
Wherein, N is the element number of array of airborne early warning radar, and M is that airborne early warning radar is sent out in a coherent processing inteval The umber of pulse penetrated.
Therefore, the transformation matrix T based on accessory channel be according to corresponding to num search passage search passage matrix S and Optimization dimensionality reduction matrix corresponding to num_s*num_t accessory channelIt is calculated, its expression formula is：
Wherein, S is search passage matrix corresponding to num search passage,For num_s × num_t accessory channel pair The optimization dimensionality reduction matrix answered, it is dimensionality reduction matrix T corresponding to num accessory channel_{b}Subset, num_s for optimization accessory channel bag The spatial frequency number contained, num_t are the Doppler frequency number that optimization accessory channel includes, and T is the change based on accessory channel Matrix is changed, H represents conjugate transposition.
Step 5, according to the transformation matrix T based on accessory channel, the twodimentional echo data received to airborne early warning radar X_{NM×L}With num search passage corresponding to search passage matrix S carry out dimensionreduction treatment respectively, respectively obtain the number of echoes after dimensionality reduction According to X_{T}With the search passage steering vector S after dimensionality reduction_{T}。
Specifically, according to the transformation matrix T based on accessory channel, the twodimentional echo data received to airborne early warning radar X_{NM×L}With num search passage corresponding to search passage matrix S carry out dimensionreduction treatment respectively, obtain the echo data X after dimensionality reduction_{T} With the search passage steering vector S after dimensionality reduction_{T}, its expression formula is respectively：
X_{T}=T^{H}X_{NM×L}
S_{T}=T^{H}S_{k}
Wherein, X_{T}For the echo data after dimensionality reduction, T is the transformation matrix based on accessory channel, X_{NM×L}For airborne early warning thunder Up to the twodimentional echo data received, S_{T}The search passage steering vector after dimensionality reduction is represented, S is corresponding to num search passage Search passage matrix S, H represent conjugate transposition.
By dimensionreduction treatment back echo data X_{T}Dimension by NM dimensionality reductions to ((num_s × num_t)+1), so as to realize machine Carry the dimensionreduction treatment for the twodimentional echo data that early warning radar receives.
Step 6, according to the echo data X after dimensionality reduction_{T}With the search passage steering vector S after dimensionality reduction_{T}, utilize maximum likelihood The filtering weight vector after dimensionality reduction is calculated in method.
Specifically, according to the echo data X after dimensionality reduction_{T}With the search passage steering vector S after dimensionality reduction_{T}, utilize maximum likelihood Echo data X after method estimation dimensionality reduction_{T}Covariance matrix, that is, utilize the echo data X after dimensionality reduction_{T}As independent same distribution Training sample, estimate to obtain the echo data X after dimensionality reduction by maximum likelihood method_{T}Covariance matrix R_{T}, its expression formula is：
L1=2 × num_s × num_t
Wherein, R_{T}Represent the echo data X after dimensionality reduction_{T}Covariance matrix, L1 represents the number of echoes of pending search passage According to required independent same distribution number of samples, X_{Tj}After representing that the echo data after the dimensionality reduction of jth of call number, j represent dimensionality reduction Echo data X_{T}Jth of call number, H represent conjugate transposition.
Adaptively weight vector W calculation formula is：
Wherein, R_{Per}For the original covariance matrix obtained using traditional accessory channel method, S_{Per}To be led to using tradition auxiliary The initial search passage steering vector that Dow process obtains.
Then, by the echo data X after dimensionality reduction_{T}Covariance matrix R_{T}The original obtained instead of using traditional accessory channel method Beginning covariance matrix R_{Per}, with the search passage steering vector S after dimensionality reduction_{T}Obtained instead of using traditional accessory channel method original Search passage steering vector S_{Per}, and the filtering weight vector W after dimensionality reduction is calculated_{T}, its expression formula is：
Wherein, W_{T}Represent the filtering weight vector after dimensionality reduction, R_{T}For the echo data X after dimensionality reduction_{T}Covariance matrix, S_{T}For Search passage steering vector after dimensionality reduction.
Step 7, according to the filtering weight vector W after dimensionality reduction_{T}, and using spacetime adaptive processing method to the echo after dimensionality reduction Data X_{T}Clutter recognition processing is carried out, obtains range Doppler figure.
Specifically, kth of search passage forms kth of transformation matrix based on accessory channel method, and num search passage is each The selfcorresponding transformation matrix based on accessory channel method, the num different transformation matrixs based on accessory channel method are formed, So as to form num different filtering weight vectors, i.e., the filtering weight vector W after dimensionality reduction_{T}.Therefore, according to the filter after dimensionality reduction Ripple weight vector W_{T}, and (STAP) method is handled to the echo data X after dimensionality reduction using spacetime adaptive_{T}Clutter recognition processing is carried out, is obtained Echo data Y to after clutter recognition processing, and the echo data Y outputs after clutter recognition is handled, obtain range Doppler Figure.
The expression formula of echo data Y after clutter recognition processing is：
Y=W_{T} ^{H}X_{T}
Wherein, W_{T}For the filtering weight vector after dimensionality reduction, X_{T}For the echo data after dimensionality reduction, H represents conjugate transposition.
Make further checking explanation to effect of the present invention with reference to emulation experiment.
(1) echo data emulation and experiment condition：
What the emulation experiment of the present invention was carried out under the softwares of MATLAB 7.11, it is airborne pre in the emulation experiment of the present invention The antenna of alert radar is using the evenly distributed linear array of 128 array element, and adjacent array element spacing and wavelength d/ λ ratio are 0.5, main ripple The antenna normal direction of Shu Zhixiang airborne early warning radars, i.e., it is 0 ° with front normal direction angle, that is, positive side array antenna. In the emulation experiment of the present invention, the threedimensional echo data that the airborne early warning radar that uses receives is according to Lincoln laboratory The clutter model simulation that J.Ward is proposed produces, and adds white Gaussian noise, and specific simulation parameter is as shown in table 1：
Table 1
(2) emulation content
In order to illustrate to improve the superiority of accessory channel algorithm, Fig. 2 (a) and Fig. 2 (b) give the place of other several algorithms Result is managed, wherein Fig. 2 (a) is using the range Doppler figure obtained after existing pulse Doppler (PD) algorithm process, Fig. 2 (b) To use the range Doppler figure obtained after the processing of original accessory channel method (ACR).
Fig. 2 (a) and Fig. 2 (b) transverse axis represents Doppler's channel position respectively, and the longitudinal axis represents range gate sequence number respectively；Fig. 2 (a) white portion in is to be distributed using the main clutter dump power obtained after existing pulse Doppler (PD) algorithm process, face The shallower subregion of color is to be distributed using the sidelobe clutter dump power obtained after existing pulse Doppler (PD) algorithm process, black Color region is to be distributed using the noise dump power obtained after existing pulse Doppler (PD) algorithm process；White in Fig. 2 (b) Region is to be distributed using the clutter dump power after the processing of original accessory channel method (ACR), and black region is at existing algorithm ACR Noise dump power distribution after reason.
From Fig. 2 (a) it can be seen that stronger residual spur occupies more distanceDoppler unit area, direct shadow The Radar Targets'Detection in the distancedoppler cells region is rung, from Fig. 2 (b) it can be seen that white portion and color are shallower Region has greatly reduced, and shows that the clutter component in echo data is suppressed.
Fig. 3 is the range Doppler figure obtained after being handled using the inventive method.White portion and color in Fig. 3 is shallower Region has greatly reduced, and shows that the clutter component in the threedimensional echo data that echo data airborne early warning radar receives obtains Preferably suppress, and independent same distribution number of training is as (2* (num+1) 3) required for original accessory channel method (ACR) It is individual be reduced to the present invention required for (2* (num_s*num_t+1) 3) it is individual, reduce required independent same distribution training sample Number.
Fig. 4 (a) is the global twodimentional response diagram obtained using the inventive method, and Fig. 4 (b) is to be obtained using the inventive method Partial enlargement two dimension response diagram；Fig. 4 (a) and Fig. 4 (b) transverse axis represents to normalize Doppler frequency, longitudinal axis difference table respectively Show normalization spatial frequency.It can be seen that can be received using the inventive method in airborne early warning radar three from Fig. 4 (b) Tie up and recess is formed at the clutter ridge that echo data is formed.
Fig. 5 is to be contrasted respectively using the clutter dump power obtained after existing algorithm PD, ACR and the inventive method processing Figure；Fig. 5 transverse axis represents Doppler's passage, and the longitudinal axis represents clutter dump power；Solid line in Fig. 5 represents more using existing pulse The clutter dump power obtained after Pu Le (PD) algorithm process, dotted line are to be obtained using after the processing of original accessory channel method (ACR) Clutter dump power, pecked line is obtained clutter dump power after being handled using the inventive method.
The inventive method is lacked relative to what the main lobe gain of original auxiliary channel algorithm reduced as can see from Figure 5, right Spatial domain error tolerance is good, and can utilize less independent same distribution number of training suppressed sidelobes clutter, and main lobe Clutter area narrows, and obtains more preferable clutter recognition effect.
The simulation result shows that the present invention can obtain preferably miscellaneous in the case where reducing training sample and array element error Ripple inhibition.
In summary, emulation experiment demonstrates the correctness of the present invention, validity and reliability.
Obviously, those skilled in the art can carry out the essence of various changes and modification without departing from the present invention to the present invention God and scope；So, if these modifications and variations of the present invention belong to the scope of the claims in the present invention and its equivalent technologies Within, then the present invention is also intended to comprising including these changes and modification.
Claims (10)
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201510717941.1A CN105223557B (en)  20151029  20151029  Airborne early warning radar clutter suppression method based on accessory channel 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201510717941.1A CN105223557B (en)  20151029  20151029  Airborne early warning radar clutter suppression method based on accessory channel 
Publications (2)
Publication Number  Publication Date 

CN105223557A CN105223557A (en)  20160106 
CN105223557B true CN105223557B (en)  20171121 
Family
ID=54992615
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201510717941.1A CN105223557B (en)  20151029  20151029  Airborne early warning radar clutter suppression method based on accessory channel 
Country Status (1)
Country  Link 

CN (1)  CN105223557B (en) 
Families Citing this family (2)
Publication number  Priority date  Publication date  Assignee  Title 

CN105929371B (en) *  20160422  20180504  西安电子科技大学  A kind of airborne radar clutter suppression method based on covariance matrix 
CN106772303B (en) *  20161222  20190201  西安电子工程研究所  The channel level clutter suppression method of MTD radar 
Citations (7)
Publication number  Priority date  Publication date  Assignee  Title 

DE3639500A1 (en) *  19861120  19880601  Forschungsgesellschaft Fuer An  Radar receiver for mobile radar apparatuses having an antenna array with clutter suppression which acts in a twodimensional manner 
US6297772B1 (en) *  19740923  20011002  The United States Of America As Represented By The Secretary Of The Navy  Predicting coherent sidelobe canceller 
CN101226236A (en) *  20071218  20080723  西安电子科技大学  Sky wave overthehorizon radar selfadaption interference rejection method based on sidelobe constraint 
CN101561497A (en) *  20090522  20091021  西安电子科技大学  Airborne radar clutter suppression method 
CN103954942A (en) *  20140425  20140730  西安电子科技大学  Method for partial combination clutter suppression in airborne MIMO radar threedimensional beam space 
CN104360325A (en) *  20141126  20150218  西安电子科技大学  Spacetime adaptive processing method for airborne forwardlooking array radar 
CN104459658A (en) *  20140619  20150325  西安电子科技大学  Unitedomain positioning spacetime twodimensional selfadaptive processing method based on data fitting 
Family Cites Families (2)
Publication number  Priority date  Publication date  Assignee  Title 

JP5193455B2 (en) *  20061106  20130508  株式会社東芝  Radar signal processing device 
JP2008157679A (en) *  20061221  20080710  Toshiba Corp  Radar signal processor 

2015
 20151029 CN CN201510717941.1A patent/CN105223557B/en active IP Right Grant
Patent Citations (7)
Publication number  Priority date  Publication date  Assignee  Title 

US6297772B1 (en) *  19740923  20011002  The United States Of America As Represented By The Secretary Of The Navy  Predicting coherent sidelobe canceller 
DE3639500A1 (en) *  19861120  19880601  Forschungsgesellschaft Fuer An  Radar receiver for mobile radar apparatuses having an antenna array with clutter suppression which acts in a twodimensional manner 
CN101226236A (en) *  20071218  20080723  西安电子科技大学  Sky wave overthehorizon radar selfadaption interference rejection method based on sidelobe constraint 
CN101561497A (en) *  20090522  20091021  西安电子科技大学  Airborne radar clutter suppression method 
CN103954942A (en) *  20140425  20140730  西安电子科技大学  Method for partial combination clutter suppression in airborne MIMO radar threedimensional beam space 
CN104459658A (en) *  20140619  20150325  西安电子科技大学  Unitedomain positioning spacetime twodimensional selfadaptive processing method based on data fitting 
CN104360325A (en) *  20141126  20150218  西安电子科技大学  Spacetime adaptive processing method for airborne forwardlooking array radar 
NonPatent Citations (4)
Title 

"A modified dimensionreduced spacetime adaptive processing method";C Luo等;《IEEE Radar Conference》;20141231;第31卷(第5期);第724728页 * 
"Adaptive airborne MTI: an auxiliary channel approach";R Klemm;《Communications Radar and Signal Processing Iee Proceedings F》;19871231;第134卷(第3期);第269276页 * 
"在非均匀环境下辅助通道法STAP处理的性能改善";王万林等;《西安电子科技大学学报（自然科学版）》;20041031;第31卷(第5期);第761764页 * 
"空时二维处理辅助通道法应用研究";杜瑞;《科技信息》;20111231;第767768页 * 
Also Published As
Publication number  Publication date 

CN105223557A (en)  20160106 
Similar Documents
Publication  Publication Date  Title 

Zhang et al.  Imaging of a moving target with rotating parts based on the Hough transform  
Melvin et al.  Assessment of multichannel airborne radar measurements for analysis and design of spacetime processing architectures and algorithms  
US7369083B2 (en)  Method for performing bistatic radar functions  
US7154433B1 (en)  Method and device for the detection and track of targets in high clutter  
Sun et al.  Direct data domain STAP using sparse representation of clutter spectrum  
US20070200750A1 (en)  Unique space time adaptive system (uss)  
Gierull et al.  Twostep detector for RADARSAT2's experimental GMTI mode  
He et al.  Novel rangespread target detectors in nonGaussian clutter  
CN101561497B (en)  Airborne radar clutter suppression method  
Folster et al.  An automotive radar network based on 77 GHz FMCW sensors  
Rohling  Some radar topics: waveform design, range CFAR and target recognition  
CN103257346A (en)  Automotive anticollision radar multitarget detecting method and system  
Smith et al.  Robust throughthewall radar image classification using a targetmodel alignment procedure  
Zhang et al.  Geometryinformationaided efficient radial velocity estimation for moving target imaging and location based on radon transform  
CN103399303B (en)  Airborne radar resists intensive deceiving jamming method and system  
EP1946145B1 (en)  Pulse doppler coherent method and system for snr enhancement  
US8232908B2 (en)  Inverse synthetic aperture radar image processing  
CN103744068B (en)  The movingtarget detection formation method of dual pathways Continuous Wave with frequency modulation SAR system  
Liang  Waveform design and diversity in radar sensor networks: theoretical analysis and application to automatic target recognition  
CN102156279A (en)  Method for detecting moving target on ground by utilizing bistatic radar based on MIMO (Multiple Input Multiple Output)  
Chen et al.  Radonfractional ambiguity functionbased detection method of lowobservable maneuvering target  
CN102721947A (en)  Efficient threedimensional spacetime adaptive clutter suppression method for airborne digital array radar  
EP2533069A1 (en)  Signal processing unit and method  
Hyun et al.  A pedestrian detection scheme using a coherent phase difference method based on 2D rangeDoppler FMCW radar  
US9594159B2 (en)  2D object detection in radar applications 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
C06  Publication  
SE01  Entry into force of request for substantive examination  
C10  Entry into substantive examination  
GR01  Patent grant  
GR01  Patent grant 