CN103809164B - A kind of rear Doppler's optimum clutter suppression method of reconstructed reference passage - Google Patents

Abstract

The invention belongs to airborne early warning radar Ground moving targets detection technical field, disclose rear Doppler's optimum clutter suppression method of a kind of reconstructed reference passage. Rear Doppler's optimum clutter suppression method of this reconstructed reference passage comprises the following steps: for each distance unit, the pulse data of each reception passage is made in time domain slide window processing; Data time domain sliding window after corresponding to each reception passage make orientation to Fourier transformation; Receive passage for the 1st and carry out time domain combined, then combine other reception channel datas and construct the data vector that each distance unit is corresponding; Ask for covariance matrix and try to achieve clutter recognition weight according to non-optimal coupling guiding; Finally utilize the clutter recognition weight tried to achieve that each distance unit is carried out clutter recognition.

Description

A kind of rear Doppler's optimum clutter suppression method of reconstructed reference passage

Technical field

The invention belongs to airborne early warning radar (AEW) Ground moving targets detection technical field, particularly to rear Doppler's optimum clutter suppression method of a kind of reconstructed reference passage.

Background technology

In modern war, grasping control of the air is important guarantee triumphantly of winning the war, and prior-warning plane plays very important effect wherein. In all previous war, prior-warning plane shows powerful power, thus being paid attention to by more people, the core of this prior-warning plane system is exactly airborne early warning (AEW) radar.

At AEW radar signal processing field, from 1973, after Brennan proposes space-time adaptive process (STAP) framework, the method of this space-time adaptive processing receives the concern of people, but owing to the dimension of this " optimal processor " is too high, making it be difficult to, the method that therefore the various quasi-optimals based on dimensionality reduction process becomes the focus that scholars study. It is that data are first carried out time-domain filtering and then combine the method that spatial domain carries out self-adaptive processing again by a class that rear Doppler (Post-Doppler) processes, such method can effectively reduce Clutter Degrees of Freedom while reducing processor dimension, is capable of quasi-optimal when goal orientation carries out accurately coupling and processes.

Time large area actual measurement contextual data is processed, if realize the suppression of scene clutter when goal orientation is mated, the amount of calculation so brought will cause that target detection is difficult to real-time implementation.Therefore, in actual treatment, people often adopt the mode that non-optimal coupling asks clutter recognition to weigh, and so can be substantially reduced operand, the detection performance of moving-target are not subject to too big loss simultaneously. Time-space processing (Filter-then-Adapt after the first sliding window that Brennan proposes, F $ A) method, when goal orientation is accurately mated, it is capable of the process performance of quasi-optimal, but in the middle of the non-optimal process of guiding inexact matching, the process performance of F $ A method is unsatisfactory.

Summary of the invention

It is an object of the invention to propose rear Doppler's optimum clutter suppression method of a kind of reconstructed reference passage. Rear Doppler's optimum clutter suppression method of this reconstructed reference passage use non-optimal coupling steering vector when carry out after doppler processing, realize Ground moving targets detection, clutter in main lobe can be carried out effective, quickly suppress, realize moving-target detection in main lobe, while reducing processor dimension, improve Ground moving targets detection performance.

For realizing above-mentioned technical purpose, the present invention adopts the following technical scheme that and is achieved.

Rear Doppler's optimum clutter suppression method of a kind of reconstructed reference passage comprises the following steps:

S1: utilize airborne early warning radar to be received back to wave datum, the echo data that described airborne early warning radar receives has the 1st distance unit to Q distance unit, and Q is the natural number more than 1; The reception passage of airborne early warning radar includes the 1st reception passage and receives passage to N, and N is the natural number more than 1; N-th receives passage receives K pulse data altogether, and n takes 1 to N; For the q distance unit of airborne early warning radar, the n-th pulse data receiving passage is made slide window processing in time domain, obtains the 1st Sliding window data of the n-th reception passageThe 3rd Sliding window data of passage is received to n-thQ takes 1 to Q;

S2: for q distance unit, receive the L time Sliding window data of passage to n-thK-2 pulse data make orientation to Fourier transformation, obtain the 1st orientation Doppler domain dataTo K-2 orientation Doppler domain dataL takes 1 to 3;

S3: rightWithMake time domain combined, draw time domain combined rear dataWherein, k is any integer value in 1 to K-2; UtilizeExtremelyConstruct the data vector x that q distance unit is corresponding^(q);

S4: utilize x⁽¹⁾To x^(Q)It is 3 (N-1)+1 that combination forms the line number of data vector matrix X, X, and columns is Q;

S5: draw clutter recognition weight W according to below equation:Wherein, S is column vector, and in S, the element of the first row is 1, and all the other elements are 0;

S6: respectively each distance unit is carried out clutter recognition according to below equation: Z^(q)=W^Hx^(q), Z^(q)It it is the result after the clutter recognition that q distance unit is corresponding.

The feature of the present invention and further improvement is that:

In step sl, n-th the 1st Sliding window data of passage is receivedReceiving passage K-2 pulse data including the n-th reception passage the 1st pulse data to n-th, n-th receives the 2nd Sliding window data of passageReceiving passage K-1 pulse data including the n-th reception passage the 2nd pulse data to n-th, n-th receives the 3rd Sliding window data of passagePassage K pulse data is received including the n-th reception passage the 3rd pulse data to n-th.

In step s3,For: $X_1^{\left(q\right)}\left(k\right)=W_t^T{[X_{\mathrm{1,1}}^{\left(q\right)}\left(k\right),X_{\mathrm{1,2}}^{\left(q\right)}\left(k\right),X_{\mathrm{1,3}}^{\left(q\right)}\left(k\right)]}^T,$ Wherein, W_t=[1,exp(j2πf_dk),exp(j4πf_dk)]^H, f_dk=η/(K-2), η are Frequency point;

In step s3, the data vector x that q distance unit is corresponding^(q)For:

$x^{\left(q\right)}={[X_1^{\left(q\right)}\left(k\right),X_{\mathrm{2,1}}^{\left(q\right)}\left(k\right),X_{\mathrm{2,2}}^{\left(q\right)}\left(k\right),X_{\mathrm{2,3}}^{\left(q\right)}\left(k\right),...,X_{N,1}^{\left(q\right)}\left(k\right),X_{N,2}^{\left(q\right)}\left(k\right){,X}_{N,3}^{\left(q\right)}\left(k\right)]}^T.$

In step sl, to each reception channel reception to echo impulse data carry out pretreatment;Described to each reception channel reception to echo impulse data carry out pretreatment and include: to each reception channel reception to echo impulse data be sequentially carried out Range compress, quadratic term compensates and doppler centroid correction; Then each reception passage is lost burst process.

The invention have the benefit that

1) moving-target detection performance in airborne early warning radar main lobe is improved.

Each passage is directly extracted identical Doppler's passage after time domain sliding window and carries out localization process by F $ A algorithm, and its reference channel is the identical Doppler's passage only differing a pulse delay in three time domains. The data that 1st receives passage are first carried out the associating in time domain by rear Doppler's optimum clutter suppression method of the reconstructed reference passage of the present invention, make to receive signal energy in passage to be added up, then pass through other reception passages and realize effective suppression of clutter in main lobe, improve moving-target detection performance.

2) computational complexity is reduced.

Compared with F $ A algorithm, in rear Doppler's optimum clutter suppression method of the reconstructed reference passage of the present invention, first carry out the associating in time domain by receiving the data of passage to the 1st, reduce computational complexity, when adopting F $ A algorithm, computational complexity is o [(3N)³], and the computational complexity of rear Doppler's optimum clutter suppression method of the reconstructed reference passage of the present invention is o [(3 (N-1)+1)³], thus the present invention improves the response speed of airborne early warning radar.

3) it is prone to Project Realization.

Rear Doppler's optimum clutter suppression method of the reconstructed reference passage of the present invention is primarily directed to F $ A algorithm hydraulic performance decline when using non-optimal coupling steering vector and proposes relatively greatly, rear Doppler's optimum clutter suppression method of the reconstructed reference passage of the present invention is more suitable for the process to large area contextual data, the suppression of scene clutter can be realized effectively and quickly, it is simple to Project Realization.

Accompanying drawing explanation

Fig. 1 is the schematic flow sheet of rear Doppler's optimum clutter suppression method of a kind of reconstructed reference passage of the present invention;

Fig. 2 is the position relationship schematic diagram of carrier aircraft and target in the embodiment of the present invention;

Fig. 3 is the first order mode schematic diagram receiving passage of the airborne early warning radar of the present invention;

Fig. 4 is the second order mode schematic diagram receiving passage of the airborne early warning radar of the present invention;

Fig. 5 is the distance-Doppler image schematic diagram of contextual data before clutter recognition;

Fig. 6 is the distance-Doppler territory residual plot after the clutter recognition after using F $ A algorithm;

Fig. 7 is the average energy value schematic diagram using F $ A algorithm to carry out each Doppler's passage of clutter recognition tailing edge;

Fig. 8 is the distance-Doppler territory residual plot using rear Doppler's optimum clutter suppression method of reconstructed reference passage of the present invention to obtain after clutter is suppressed;

Fig. 9 is the average energy value schematic diagram using rear Doppler's optimum clutter suppression method of reconstructed reference passage of the present invention to carry out each Doppler's passage of clutter recognition tailing edge;

Figure 10 is the schematic diagram of two the improvement factor curves obtained after clutter is suppressed by rear Doppler's optimum clutter suppression method of the reconstructed reference passage using F $ A algorithm and the present invention.

Detailed description of the invention

Below in conjunction with accompanying drawing, the invention will be further described:

With reference to Fig. 1, for the schematic flow sheet of rear Doppler's optimum clutter suppression method of a kind of reconstructed reference passage of the present invention. Rear Doppler's optimum clutter suppression method of this reconstructed reference passage comprises the following steps:

S1: utilize airborne early warning radar to be received back to wave datum, the echo data that described airborne early warning radar receives has the 1st distance unit to Q distance unit, and Q is the natural number more than 1.Airborne early warning radar by each distance echo impulse data (initial data) corresponding to unit of multiple reception channel reception, to each reception channel reception to echo impulse data carry out pretreatment. To each reception channel reception to echo impulse data carry out pretreatment and include: to each reception channel reception to echo impulse data be sequentially carried out Range compress, quadratic term compensates and doppler centroid correction. It is described as follows:

First to each reception channel reception to echo impulse data carry out Range compress: set the transmitting signal of airborne early warning radar as s (t), then through t₀The echo-signal received after time delay is represented by:

$s\left(t\right)=\mathrm{rect}\left(\frac{t-t_0}{T}\right)\exp \left\{\mathrm{j\π}{K}_r{(t-t_0)}^2\right\}$

Wherein, T is the pulse width of echo-signal, K_rFrequency modulation rate for the linear FM signal of airborne early warning radar. Generally, carry out matched filtering at frequency domain more convenient, signal spectrum s can be obtained hence with in facies principle (POSP)_rF () is approximately:

$S_r\left(f\right)=\mathrm{rect}\left(\frac{f}{\left|K_r\right|T}\right)\exp \{-\mathrm{j\π}\frac{f^2}{K_r}\}\exp \{-j2\mathrm{\πf}{t}_0\}$

Wherein, f represents frequency, thus can obtain Range compress frequency matching function H (f) and be:

$H\left(f\right)=\mathrm{rect}\left(\frac{f}{\left|K_r\right|T}\right)\exp \{+\mathrm{j\π}\frac{f^2}{K_r}\}$

With reference to Fig. 2, for the position relationship schematic diagram of carrier aircraft in the embodiment of the present invention Yu target. In embodiments of the present invention, set up three-dimensional cartesian coordinate system, the initial point of this three-dimensional cartesian coordinate system is the transmission channel projection on the ground of airborne early warning radar, the X-axis forward of three-dimensional cartesian coordinate system is parallel with the aircraft motion direction of airborne early warning radar, and the Z axis forward of three-dimensional cartesian coordinate system is direction straight up. The reception passage of airborne early warning radar includes the 1st reception passage and receives passage to N, and N is the natural number more than 1, and the 1st receives passage is arranged in order to N reception passage according to the direction of motion of carrier aircraft. The carrier aircraft height of airborne early warning radar is H, and for moving-target point, its azimuth is θ, and the angle of pitch isR₀Represent the shortest oblique distance. This moving-target o'clock receives the oblique distance R of passage relative to n-th_n(t) be:

$R_n\left(t\right)\≈R-V_{r}t+\frac{{(V_{\mathrm{rx}}t+d_n-x_0)}^2}{2R}$

$=R-V_{r}t+\frac{{(V_{\mathrm{rx}}t-x_0)}^2}{2R}+\frac{(V_{\mathrm{rx}}t-x_0)d_n}{R}+\frac{d_n^2}{2R}$

Wherein, n takes 1 to N, t express time, and R is the moving-target point oblique distance to the transmission channel of airborne early warning radar, d_nThe transmission channel of passage and airborne early warning radar interval on course (course of the carrier aircraft of airborne early warning radar), x is received for the n-th of airborne early warning radar₀Representing the X-coordinate of moving-target point, the carrier aircraft speed of airborne early warning radar is V, the speed of moving-target point component respectively V on X-axis, Y-axis, Z axis_x、V_yAnd V_z, the radial velocity (carrier aircraft relative to airborne early warning radar) of moving-target point is V_r, V_rx=V-V_x。

Space-time adaptive process in, when the coherent processing inteval of CPI(airborne early warning radar) shorter time, R_nT in (), the impact of the quadratic term of time can be ignored, but when CPI is longer, the impact of quadratic term can not be ignored, it is necessary to first quadratic term is compensated, by R_nT the computing formula of () can obtain corresponding quadratic term compensation reference function ref:

$\mathrm{ref}=\exp \left(j\frac{2\mathrm{\π}{V}^2{t}^2}{\mathrm{\λR}}\right)$

If airborne early warning radar is to receive echo impulse data when big stravismus, it is possible to cause doppler centroid not at zero-frequency. Doppler centroid can be estimated first with median method or averaging method, then the mid frequency estimated be compensated and go back.

To each reception channel reception to echo impulse data carry out pretreatment after, each reception passage is lost burst process. It is described as follows:

Formula below is utilized to calculate the umber of pulse that each reception passage needs to lose:

$\mathrm{num}=\frac{d}{2}\×\frac{f_r}{V}$

Wherein, d is the adjacent reception channel pitch of airborne early warning radar, f_rFor the pulse recurrence frequency of echo impulse data that airborne early warning radar receives, after drawing num, according to the position relationship between each passage and carrier aircraft course, calculate the umber of pulse that each passage needs to lose, and carry out losing burst process.

Be arranged in order according to the direction of motion of carrier aircraft owing to the 1st reception passage receives passage to N, then the 1st reception passage presents first order mode or second order mode to N reception passage; With reference to Fig. 3, for the first order mode schematic diagram receiving passage of the airborne early warning radar of the present invention, reference Fig. 4, for the second order mode schematic diagram receiving passage of the airborne early warning radar of the present invention. In first order mode, the 1st receives passage receives the passage front closer to the carrier aircraft of airborne early warning radar relative to N, and in second order mode, the 1st receives passage receives the passage rear closer to the carrier aircraft of airborne early warning radar relative to N. When the 1st reception passage to N reception passage presents first order mode, the n-th of airborne early warning radar is received passage, lose (n-1) × num the pulse data of its section start, and lose (N-n) × num the pulse data at its end place. When the 1st reception passage to N reception passage presents second order mode, the n-th of airborne early warning radar is received passage, lose (N-n) × num the pulse data of its section start, and lose (n-1) × num the pulse data at its end place. After losing burst process, it may be achieved each receives channel data spatial alignment, and improves clutter recognition performance.

After each reception passage is lost burst process, the pulse data of each reception passage is carried out slide window processing. It is described as follows:

N-th receives passage receives K pulse data altogether; For q distance unit, the n-th pulse data receiving passage being made in time domain slide window processing, q takes 1 to Q; Obtaining K-2 pulse data after slide window processing, during slide window processing, the time span of each slip is 1 pulse; Obtain after using slide window processing: n-th receives the 1st Sliding window data of passageN-th receives the 2nd Sliding window data of passageAnd n-th receive the 3rd Sliding window data of passageN-th receives the 1st Sliding window data of passagePassage K-2 pulse data (corresponding sliding window initial position) is received including the n-th reception passage the 1st pulse data to n-th. N-th receives the 2nd Sliding window data of passageReceiving passage K-1 pulse data (corresponding sliding window first time slide result) including the n-th reception passage the 2nd pulse data to n-th, n-th receives the 3rd Sliding window data of passagePassage K pulse data (corresponding sliding window second time slide result) is received including the n-th reception passage the 3rd pulse data to n-th.

S2: for q distance unit, receive the L time Sliding window data of passage to n-thK-2 pulse data make orientation to Fourier transformation, obtain the 1st orientation Doppler domain dataTo K-2 orientation Doppler domain dataL takes 1 to 3.

Specifically, for q distance unit, the 1st Sliding window data of passage is received to n-thK-2 pulse data (include the n-th reception passage the 1st pulse data to n-th receive passage K-2 pulse data) make orientation to Fourier transformation, obtainExtremelyFor q distance unit, receive the 2nd Sliding window data of passage to n-thK-2 pulse data (corresponding n-th receives passage the 2nd pulse data receives passage K-1 pulse data to n-th) make orientation to Fourier transformation, obtainExtremelyFor q distance unit, receive the 3rd Sliding window data of passage to n-thK-2 pulse data (corresponding n-th receives passage the 3rd pulse data receives passage K pulse data to n-th) make orientation to Fourier transformation, obtainExtremely

S3: the 1st Sliding window data receiving passage (i.e. reference channel) is made time domain combined, namely according to below equation pairWithMake time domain combined, draw time domain combined rear data

$X_1^{\left(q\right)}\left(k\right)=W_t^T{[X_{\mathrm{1,1}}^{\left(q\right)}\left(k\right),X_{\mathrm{1,2}}^{\left(q\right)}\left(k\right),X_{\mathrm{1,3}}^{\left(q\right)}\left(k\right)]}^T$

Wherein, W_t=[1,exp(j2πf_dk),exp(j4πf_dk)]^H, f_dk=η/(K-2);η is Frequency point; When K-2 is even number, η isArriveArbitrary integer; When K-2 is odd number, η isArriveArbitrary integer.

Then utilizeExtremelyConstruct the data vector x that q distance unit is corresponding^(q):

$x^{\left(q\right)}={[X_1^{\left(q\right)}\left(k\right),X_{\mathrm{2,1}}^{\left(q\right)}\left(k\right),X_{\mathrm{2,2}}^{\left(q\right)}\left(k\right),X_{\mathrm{2,3}}^{\left(q\right)}\left(k\right),...,X_{N,1}^{\left(q\right)}\left(k\right),X_{N,2}^{\left(q\right)}\left(k\right){,X}_{N,3}^{\left(q\right)}\left(k\right)]}^T.$

This shows the data vector x that q distance unit is corresponding^(q)For column vector, its element number is 3 (N-1)+1.

S4: owing to q takes 1 to Q, after step S3, is just obtaining x⁽¹⁾(i.e. the 1st distance data vector corresponding to unit of airborne early warning radar) is to x^(Q)(data vector that namely the Q distance unit of airborne early warning radar is corresponding). It is clear that at x⁽¹⁾To x^(Q)In, each data vector is column vector, and its element number is 3 (N-1)+1.

At this moment, x is utilized⁽¹⁾To x^(Q)It is 3 (N-1)+1 that combination forms the line number of data vector matrix X, X, and columns is Q. I.e. X=[x⁽¹⁾,...,x^(Q)]。

S5: draw clutter recognition weight W according to below equation:

$W=\frac{R^{-1}S}{S^H{R}^{-1}S},R=\frac{1}{Q}X{X}^H,$

Wherein, S is column vector, and in S, the element of the first row is 1, and all the other elements are 0; S is non-optimal coupling guiding.

S6: respectively each distance unit is carried out clutter recognition according to below equation: Z^(q)=W^Hx^(q), Z^(q)It it is the result after the clutter recognition that q distance unit is corresponding.

Below by emulation experiment, the invention will be further described:

1) experimental data

This emulation experiment adopts 4 reception passage 512 pulse airborne radar measured datas, carrier aircraft speed is 107.0570m/s, the spacing of adjacent reception passage is 0.27m, corresponding pulse recurrence frequency is 833.3333Hz, wavelength is 0.0313m, signal bandwidth 18MHz, pulse width 40 μ s, distance samples frequency 20MHz, linear FM signal frequency modulation rate is K_r=4.5×10¹¹, range gate width 7.5m, Range compress front distance unit adds up to 8192, the shortest oblique distance 20000m.

2) experiment content and result:

Experiment 1:

Initial data is carried out pretreatment; Draw the distance-Doppler image of contextual data before clutter recognition, with reference to Fig. 5, for the distance-Doppler image schematic diagram of contextual data before clutter recognition. In Fig. 5, gray level image represents clutter, and gray value is more big, then illustrate that energy is more big. When carrying out experiment 1,1100 distance unit (i.e. the 1st distance unit to the 1100th apart from unit), 34 pulse datas (the 230th pulse data to the 263rd pulse data) are intercepted, as seen from Figure 5, before carrying out clutter recognition, in airborne early warning radar main lobe, clutter is very strong, and target is submerged, and secondary lobe district clutter is more weak, the detection of moving-target can be made directly, if main-lobe clutter is effectively suppressed, it is possible to realize effective detection of main lobe internal object.

After initial data is carried out pretreatment, use rear Doppler's optimum clutter suppression method of the reconstructed reference passage of F $ A algorithm and the present invention that clutter is suppressed respectively. With reference to Fig. 6, for the residual plot of the distance-Doppler image after the clutter recognition after use F $ A algorithm. In figure 6, gray level image represents target, and gray value is more big, then illustrate that energy is more big. Although from fig. 6 it can be seen that F $ A algorithm is capable of the suppression of clutter, but signal energy corresponding with target after suppressing is more weak. It is use F $ A algorithm to carry out the average energy value schematic diagram of each Doppler's passage of clutter recognition tailing edge with reference to Fig. 7. In the figure 7, distance unit (abscissa) only takes 600 to 1100. From figure 7 it can be seen that compared with carrying out before clutter recognition, after using F $ A algorithm to carry out clutter recognition, it is possible to target and clutter are distinguished.

With reference to Fig. 8, for using the residual plot of distance-Doppler image that rear Doppler's optimum clutter suppression method of the reconstructed reference passage of the present invention obtains after clutter is suppressed. In Fig. 8, gray level image represents target, and gray value is more big, then illustrate that energy is more big. Compared with Fig. 6, after using rear Doppler's optimum clutter suppression method of reconstructed reference passage of the present invention that clutter is suppressed, the signal energy corresponding with target obtains enhancing. With reference to Fig. 9, for using rear Doppler's optimum clutter suppression method of the reconstructed reference passage of the present invention to carry out the average energy value schematic diagram of each Doppler's passage of clutter recognition tailing edge. In fig .9, distance unit (abscissa) only takes 600 to 1100. From fig. 9, it can be seen that compared with carrying out before clutter recognition, after using the present invention to carry out clutter recognition, peak value is more prominent, illustrates to be easier to distinguish target and clutter.

Experiment 2:

In this experiment, we adopt the mode adding target in scene, obtain the improvement factor curve of its correspondence by changing the speed of added target, thus verifying the effectiveness of rear Doppler's optimum clutter suppression method of the reconstructed reference passage of the present invention. Here moving-target is added to No. 17 Doppler's passage (be provided with in this experiment No. 1 Doppler's passage to No. 32 Doppler's passage) of the contextual data intercepted, change target velocity, then use rear Doppler's optimum clutter suppression method of the reconstructed reference passage of F $ A algorithm and the present invention that clutter is suppressed respectively. Draw two improvement factor curves of correspondence.

With reference to Figure 10, for using the schematic diagram of two improvement factor curves that rear Doppler's optimum clutter suppression method of the reconstructed reference passage of F $ A algorithm and the present invention obtains after clutter is suppressed. From fig. 10 it can be seen that the clutter recognition performance that the present invention is in main lobe district has greatly improved compared with F $ A, clutter recognition recess is narrower, and output signal energy is higher.

In sum, rear Doppler's optimum clutter suppression method of the reconstructed reference passage of the present invention is under the premise reducing computational complexity, the effective detection being capable of in main lobe moving-target, and the method to the accumulation degree of signal energy more than the accumulation to clutter energy, after making clutter recognition, signal energy is higher, output letter miscellaneous noise ratio (SCNR) is higher, it is simple to follow-up moving-target detects.

Obviously, the present invention can be carried out various change and modification without deviating from the spirit and scope of the present invention by those skilled in the art. So, if these amendments of the present invention and modification belong within the scope of the claims in the present invention and equivalent technologies thereof, then the present invention is also intended to comprise these change and modification.

Claims (3)

1. rear Doppler's optimum clutter suppression method of a reconstructed reference passage, it is characterised in that comprise the following steps:

S1: utilize airborne early warning radar to be received back to wave datum, the echo data that described airborne early warning radar receives has the 1st distance unit to Q distance unit, and Q is the natural number more than 1; The reception passage of airborne early warning radar includes the 1st reception passage and receives passage to N, and N is the natural number more than 1; N-th receives passage receives K pulse data altogether, and n takes 1 to N; For the q distance unit of airborne early warning radar, the n-th pulse data receiving passage is made slide window processing in time domain, obtains the 1st Sliding window data of the n-th reception passageThe 3rd Sliding window data of passage is received to n-thQ takes 1 to Q;

S2: for q distance unit, receive the L time Sliding window data of passage to n-thK-2 pulse data make orientation to Fourier transformation, obtain the 1st orientation Doppler domain dataTo K-2 orientation Doppler domain dataL takes 1 to 3;

S3: rightWithMake time domain combined, draw time domain combined rear dataWherein, k is any integer value in 1 to K-2;

$X_1^{\left(q\right)}\left(k\right)=W_t^T{\[X_{1,1}^{\left(q\right)}\left(k\right),X_{1,2}^{\left(q\right)}\left(k\right),X_{1,3}^{\left(q\right)}\left(k\right)\]}^T$

Wherein, W_t=[1, exp (j2 π f_dk),exp(j4πf_dk)]^H, f_dk=η/(K-2); η is Frequency point; When K-2 is even number,ArriveArbitrary integer; When K-2 is odd number, η isArriveArbitrary integer;

UtilizeExtremelyConstruct the data vector x that q distance unit is corresponding^(q), wherein, the data vector x that q distance unit is corresponding^(q)For:

$x^{\left(q\right)}={\[X_1^{\left(q\right)}\left(k\right),X_{2,1}^{\left(q\right)}\left(k\right),X_{2,2}^{\left(q\right)}\left(k\right),X_{2,3}^{\left(q\right)}\left(k\right),\mathrm{...},X_{N,1}^{\left(q\right)}\left(k\right),X_{N,2}^{\left(q\right)}\left(k\right),X_{N,3}^{\left(q\right)}\left(k\right)\]}^T;$

S4: utilize x⁽¹⁾To x^(Q)It is 3 (N-1)+1 that combination forms the line number of data vector matrix X, X, and columns is Q;

S5: draw clutter recognition weight W according to below equation:Wherein, S is column vector, and in S, the element of the first row is 1, and all the other elements are 0;

S6: respectively each distance unit is carried out clutter recognition according to below equation: Z^(q)=W^Hx^(q), Z^(q)It it is the result after the clutter recognition that q distance unit is corresponding.

2. rear Doppler's optimum clutter suppression method of a kind of reconstructed reference passage as claimed in claim 1, it is characterised in that in step sl, n-th receives the 1st Sliding window data of passageReceiving passage K-2 pulse data including the n-th reception passage the 1st pulse data to n-th, n-th receives the 2nd Sliding window data of passageReceiving passage K-1 pulse data including the n-th reception passage the 2nd pulse data to n-th, n-th receives the 3rd Sliding window data of passagePassage K pulse data is received including the n-th reception passage the 3rd pulse data to n-th.

3. rear Doppler's optimum clutter suppression method of a kind of reconstructed reference passage as claimed in claim 1, it is characterised in that in step sl, to each reception channel reception to echo impulse data carry out pretreatment; Described to each reception channel reception to echo impulse data carry out pretreatment and include: to each reception channel reception to echo impulse data be sequentially carried out Range compress, quadratic term compensates and doppler centroid correction; Then each reception passage is lost burst process.