CN104459657B - Extension factorization space-time two-dimensional self-adaptive processing method based on data fitting - Google Patents
Extension factorization space-time two-dimensional self-adaptive processing method based on data fitting Download PDFInfo
- Publication number
- CN104459657B CN104459657B CN201410276078.6A CN201410276078A CN104459657B CN 104459657 B CN104459657 B CN 104459657B CN 201410276078 A CN201410276078 A CN 201410276078A CN 104459657 B CN104459657 B CN 104459657B
- Authority
- CN
- China
- Prior art keywords
- data
- represent
- detected
- matrix
- range cell
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses an extension factorization space-time two-dimensional self-adaptive processing method based on data fitting and relates to radar technology. The method comprises the specific steps that 1, dimensionality reduction echo data and space-time guide vectors are obtained; 2, a data basis matrix of a distance unit to be detected is constructed; 3, a data basis matrix with targets blocked is obtained; 4, the fitting coefficient of auxiliary echo data is obtained; 5, the minimum fitting error of data of the distance unit to be detected is obtained; 6, unit average constant false alarm rate detection is carried out on the minimum fitting error of the data of the distance unit to be detected. The extension factorization space-time two-dimensional self-adaptive processing method mainly solves the problem that a traditional STAP method is huge in computation burden and strict in sample requirement, and has the good application prospect in measured data processing and STAP performance improvement.
Description
Technical field
The invention belongs to communication technical field, be related to Radar Technology, particularly to a kind of extension based on data matching because
Sonization space-time adaptive processing method, for airborne early warning Radar Signal Processing.
Background technology
The main task of airborne early warning radar is detection target in complex clutter background, and carries out locating and tracking to it,
And it is the core means improving airborne early warning radar service behaviour that clutter is effectively suppressed.Space-time adaptive is processed
(space-time adaptive processing, stap) technology makes full use of spatially and temporally information, to echo signal
While carrying out coherent accumulation, processed by space-time adaptive and filter ground clutter, realizing airborne early warning radar has to target
Effect detection, such as the e2-d airborne early warning radar of the U.S. just adopt this technology.
In actual applications, stap technology is primarily present following two aspect problems: on the one hand, in clutter environment heterogeneous
In, enough independent same distribution (independent and identically for estimate covariance matrix will be obtained
Distributed, iid) training sample is extremely difficult;On the other hand, even if the demand of training sample is met, at full space-time
The excessive problem of reason amount of calculation can lead to real-time to be difficult to ensure that.
For solving the above problems, promote stap technology more practical, there has been proposed many corrective measures or method.Clearly
Patent of invention " the space-time adaptive processing method under non-homogeneous clutter environment " (number of patent application of Hua Da application
201010129723.3, publication No. cn101819269a) and disclose a kind of super-resolution estimation clutter in non-homogeneous clutter environment
The overcomplete sparse representation method of space-time two-dimensional spectrum.The method achieve in the case of independent same distribution sample number deficiency, profit
Estimate clutter covariance matrix with single frames training sample, thus avoiding the shadow to self-adaptive processing effect for the strong non-homogeneous clutter environment
Ring.But, the method yet suffers from following main deficiency is: one is the big problem of operand, carries out rarefaction representation to clutter spectrum
Super complete radix mesh uncertain, but be much larger than degree of freedom in system, and degree of freedom in system is generally thousands of in practice, so every
In the covariance matrix restructuring procedure of one range cell sample, required operand is very big, is unfavorable for real-time processing, thus
Have influence on the effect in practical engineering application.Two is base mismatch problems, and used in the method, base is actually one group of interpolation
Discrete fourier dft vector, and the echo data of reality because there is the non-ideal factors such as error it is impossible to base vector
Rarefaction representation;In addition, the method needs to carry out discretization expression to space-time plane, such that not being located in discrete grid block point
The signal of the heart is revealed to all mesh points, thus destroying the openness of echo data.
Content of the invention
Overcomplete sparse representation side for above-mentioned super-resolution estimation clutter space-time two-dimensional spectrum in non-homogeneous clutter environment
The deficiency of method, the present invention proposes a kind of spreading factor space-time adaptive processing based on data matching
(extendedfactored approach stap, efa-stap) method, is capable of the detection to echo signal, and solves
The huge problem with base mismatch of operand in prior art, simultaneously with respect to conventional Extension factorization efa-stap method, can
Improve the detection probability of moving-target, reduce the false alarm rate of detection.
For reaching above-mentioned purpose, the present invention employ the following technical solutions pre- with realize.
A kind of spreading factor space-time adaptive processing method based on data matching it is characterised in that include with
Lower step:
Step 1, receives the space-time two-dimensional echo data x of airborne early warning radar, to space-time using airborne early warning radar antenna
Two-dimentional echo data x is extended factorization efa dimensionality reduction, obtains dimensionality reduction echo data z and space-time steering vector sz-efa;
Step 2, constructs the data basic matrix φ of range cell to be detectedk;
Step 3, using space-time steering vector sz-efaDerive blocking matrix b;Using blocking matrix b and range cell to be detected
Data zk, obtain target and block later auxiliary echo dataData z of this range cell to be detectedkTake from step successively
Dimensionality reduction echo data z in rapid 1;Data basic matrix φ using blocking matrix b and range cell to be detectedkObtain target to block
Later data basic matrix
Step 4, blocks later data basic matrix using targetTo auxiliary echo dataCarry out data matching, obtain
Take auxiliary echo dataFitting coefficient
Step 5, using the data basic matrix φ of range cell to be detectedkWith auxiliary echo dataFitting coefficientRight
Data z of range cell to be detectedkCarry out data matching, obtain data z of range cell to be detectedkMinimum error of fitting yk;
Data z of this range cell to be detectedkTake from the dimensionality reduction echo data z in step 1 successively;
Step 6, treats minimum error of fitting y of the data of detecting distance unitkCarry out CA-CFAR detection, and
There is target or there is not target in output.
The feature of technique scheme and further improvement is that:
(1) step 1 specifically includes following sub-step:
1a) utilize airborne early warning radar antenna, receive the space-time two-dimensional number of echoes of ground return within the coherent accumulation time
According to x;Space-time two-dimensional echo data x ties up for mn, and wherein m represents umber of pulse, and n represents the array number of airborne early warning radar antenna;
1b) be expanded the dimensionality reduction transition matrix p of factorization efa dimensionality reduction according to following formulaefa:
Wherein, pt(fi-1,fi,fi+1) represent the time domain transition matrix being made up of three Doppler filters closing on, fiTable
Show the normalization Doppler frequency of target, fi-1,fi+1Represent two normalization Doppler frequencies adjacent with target, i be less than
Natural number equal to m and more than or equal to 1, m represents umber of pulse, inRepresent the unit matrix of n × n,Represent Kronecker product,
N represents the array number of airborne early warning radar antenna, dimensionality reduction transition matrix pefaDimension be mn × 3n tie up;
1c) according to dimensionality reduction transition matrix pefaBe expanded the dimensionality reduction echo data z after factorization efa dimensionality reduction, is expressed as down
Formula:
Wherein, ()hRepresent conjugate transposition operation, according to dimensionality reduction transition matrix pefaDimension be mn × 3n dimension and empty
The dimension that Shi Erwei echo data x can obtain dimensionality reduction echo data z for mn dimension is tieed up for 3n;
Space-time steering vector s after spreading factor efa dimensionality reductionz-efa, it is expressed as following formula:
ss(fs)=[1 exp (j2 π fs) exp(j2πfs2) … exp(j2πfs(n-1))]t
Wherein, ()tRepresent transposition operation, ss(fs) represent target spatial domain steering vector, fsRepresent the normalization of target
Spatial frequency, n represents the array number of airborne early warning radar antenna.
(2) step 2 specifically includes following sub-step:
Data basic matrix φkIt is expressed as following formula:
φk=[z1z2… zl]∈cd×l
Wherein, φkRepresent the data basic matrix of k-th range cell to be detected;K represents the data of range cell to be detected
Sequence number, be less than the natural number of all range cells to be detected sums;D represents the dimension of dimensionality reduction echo data z, and l represents
Close on the number of the data of range cell, cd×lRepresent the complex matrix space of d × l dimension, zl(l=1,2 ..., l) represent close on away from
From the data of unit, close on range cell and represent range cell near range cell to be detected, do not comprise two and be used for protecting
The range cell of target.
(3) step 3 specifically includes following sub-step:
3a) to space-time steering vector sz-efaTransposed form carry out singular value decomposition, obtain left singular matrix u, right unusual
Matrix v and singular value matrix λ;Decomposition formula such as following formula:
Wherein, u represents left singular matrix, and v represents right singular matrix, and λ represents singular value matrix, ()tRepresent transposition behaviour
Make, ()hRepresent conjugate transposition operation;
3b) blocking matrix b is arranged to make up by the 2nd to the n row order of right singular matrix v, by right singular matrix v and following formula
Derive blocking matrix b;
B=[v (:, 2:n)]t
Wherein, v (:, 2:n) represents all row of right singular matrix v, the 2nd arrives n column matrix element;(·)tRepresent transposition behaviour
Make;
3c) auxiliary echo data is obtained according to following formulaBlock later data basic matrix with target
Wherein, zkRepresent the data of k-th range cell to be detected, the data of this range cell to be detected takes from step successively
Dimensionality reduction echo data z in rapid 1, k represent the sequence number of the data of range cell to be detected, are single less than all distances to be detected
The natural number of unit's sum;Represent that target blocks later auxiliary echo data,Represent that target blocks later data group moment
Battle array, b is blocking matrix.
(4) step 4 specifically includes following sub-step:
4a) block later data basic matrix with targetLater auxiliary echo data is blocked to targetCarry out
A young waiter in a wineshop or an inn takes advantage of matching to represent, is expressed as following formula:
Wherein, | | | |2Represent 2 norms solving vector,Represent auxiliary echo dataFitting coefficient;
4b) according to step 4a) in least square fitting represent, solve auxiliary echo dataFitting coefficient?
A young waiter in a wineshop or an inn takes advantage of solutionIt is expressed as following formula:
Wherein, ()-1Representing matrix inversion operation, ()hRepresent conjugate transposition operation,After representing that target is blocked
Data basic matrix,Represent that target blocks later auxiliary echo data.
(5) step 5 specifically includes following sub-step:
5a) with the data basic matrix φ of range cell to be detectedk, treat data z of detecting distance unitkCarry out a young waiter in a wineshop or an inn
Take advantage of matching;With solution auxiliary echo dataFitting coefficientLeast square solutionIt is calculated optimization restrained boundary ηk, tool
Body is realized by following formula:
Wherein, αkRepresent the fitting coefficient of the data of range cell to be detected, ykThe 2 of the minimum error of fitting vector of expression
Norm, ηkRepresent and optimize restrained boundary,
5b) by solving the optimal fitting coefficient that above formula obtains the data of range cell to be detected
The optimal fitting coefficient of the data according to range cell to be detectedObtain the data of range cell to be detected
Little error of fitting yk, obtained by following formula:
Wherein, zkRepresent the data of k-th range cell to be detected, the data of this range cell to be detected takes from step successively
Dimensionality reduction echo data z in rapid 1, k represent the sequence number of the data of range cell to be detected, are single less than all distances to be detected
The natural number of unit's sum;φkRepresent the data basic matrix of range cell to be detected;Represent the data of range cell to be detected
Optimal fitting coefficient;||·||2Represent 2 norms solving vector.
Compared with prior art, the present invention has prominent substantive distinguishing features and significantly improves.The present invention is dilute with existing
Thin recovery stap method is compared and is had the advantage that
(1) it is directed to the problem of operand, existing sparse restoration methods carry out the super complete radix of rarefaction representation to clutter spectrum
Mesh is uncertain, but is much larger than degree of freedom in system, and degree of freedom in system is generally thousands of in practice, so in each range cell sample
In this covariance matrix restructuring procedure, required operand is very big, is unfavorable for real-time processing, thus having influence on its reality
Border engineer applied effect.Stap method is combined by the present invention with efa, and the optimization problem of original mn dimension is decomposed into m 3n dimension
Little optimization problem to solve, greatly reduce the computation complexity of algorithm, solve operand in sparse recovery stap method
Big problem.
(2) it is directed to base mismatch problems, used in existing sparse restoration methods, represent actually one group interpolation dft arrow of base
Amount, and real data because there is the non-ideal factors such as error it is impossible to use given base vector rarefaction representation, in addition, the method need
Discretization expression is carried out to space-time plane, such that the signal not being located at discrete grid block dot center is let out to all mesh points
Dew, thus destroy the openness of data.The present invention is directly used real data sample as basic matrix, rather than adopts certain
Plant the ideal flowing pattern determining, thus having evaded the problem of base mismatch in above-mentioned sparse recovery stap method.
It can be seen that, in prior art, stap carries out airborne radar signal processing generally in high-dimensional data space, thus can
Greatly increase computational complexity and the problem of required training sample number, in order to obtain better performance, using to be detected
The data sample that range cell closes on range cell carrys out Optimal Fitting range cell to be detected clutter data, thus improving self adaptation
The clutter recognition performance of signal processing, improves the detection probability of target.
Brief description
The present invention will be further described with reference to the accompanying drawings and detailed description.
Fig. 1 is flow chart of the present invention;
Fig. 2 is general pulse doppler processing result figure;
Fig. 3 is the change curve with border compensation dosage for the detection probability;
Fig. 4 is the change curve with sample number for the detection probability;
Fig. 5 is the output result figure of the minimum error of fitting of data of target proximity range cell;Fig. 5 (a) is through normal
The output result of the minimum error of fitting of a segment distance unit scope after rule efa method;Fig. 5 (c) is through transmission from one meridian to another side of the present invention
The output result of the minimum error of fitting of a segment distance unit scope after method process;Fig. 5 (b) be through conventional efa method it
The output result of the minimum error of fitting of another segment distance unit scope afterwards;Fig. 5 (d) is another section after the inventive method
The output result of the minimum error of fitting of range cell scope;
Fig. 6 is the change curve with target power for the detection probability;
Fig. 7 is the change curve graph of a relation with false alarm rate for the detection probability.
Specific embodiment
With reference to Fig. 1, a kind of spreading factor space-time adaptive processing side based on data matching of the present invention is described
Method, the present invention is used for airborne early warning Radar Signal Processing, and its specific implementation step is as follows:
Step 1, receives the space-time two-dimensional echo data x of airborne early warning radar, to space-time using airborne early warning radar antenna
Two-dimentional echo data x is extended factorization efa dimensionality reduction, obtains dimensionality reduction echo data z and space-time steering vector sz-efa;
Step 1 specifically includes following sub-step:
1a) utilize airborne early warning radar antenna, receive the space-time two-dimensional number of echoes of ground return within the coherent accumulation time
According to x;Space-time two-dimensional echo data x ties up for mn, and wherein m represents umber of pulse, and n represents the array number of airborne early warning radar antenna;
1b) be expanded the dimensionality reduction transition matrix p of factorization efa dimensionality reduction according to following formulaefa:
Wherein, pt(fi-1,fi,fi+1) represent the time domain transition matrix being made up of three Doppler filters closing on, fiTable
Show the normalization Doppler frequency of target, fi-1,fi+1Represent two normalization Doppler frequencies adjacent with target, i be less than
Natural number equal to m and more than or equal to 1, m represents umber of pulse, inRepresent the unit matrix of n × n,Represent Kronecker product,
N represents the array number of airborne early warning radar antenna, dimensionality reduction transition matrix pefaDimension be mn × 3n tie up;
1c) according to dimensionality reduction transition matrix pefaBe expanded the dimensionality reduction echo data z after factorization efa dimensionality reduction, is expressed as down
Formula:
Wherein, ()hRepresent conjugate transposition operation, according to dimensionality reduction transition matrix pefaDimension be mn × 3n dimension and empty
The dimension that Shi Erwei echo data x can obtain dimensionality reduction echo data z for mn dimension is tieed up for 3n;
Space-time steering vector s after spreading factor efa dimensionality reductionz-efa, it is expressed as following formula:
ss(fs)=[1 exp (j2 π fs) exp(j2πfs2) … exp(j2πfs(n-1))]t
Wherein, ()tRepresent transposition operation, ss(fs) represent target spatial domain steering vector, fsRepresent the normalization of target
Spatial frequency, n represents the array number of airborne early warning radar antenna.
General space-time steering vector s is represented by the Kronecker product form of time domain steering vector and spatial domain steering vector,
I.e.fdAnd fsRepresent normalized Doppler frequency and spatial frequency, time domain and spatial domain guiding arrow respectively
Amount is represented by respectively
Step 2, constructs the data basic matrix φ of range cell to be detectedk;
Data basic matrix φkIt is expressed as following formula:
φk=[z1z2… zl]∈cd×l
Wherein, φkRepresent the data basic matrix of k-th range cell to be detected;K represents the data of range cell to be detected
Sequence number, be less than the natural number of all range cells to be detected sums;D represents the dimension of dimensionality reduction echo data z, and l represents
Close on the number of the data of range cell, cd×lRepresent the complex matrix space of d × l dimension, zl(l=1,2 ..., l) represent close on away from
From the data of unit, close on range cell and represent range cell near range cell to be detected, do not comprise two and be used for protecting
The range cell of target.
Step 3, using space-time steering vector sz-efaDerive blocking matrix b;Using blocking matrix b and range cell to be detected
Data zk, obtain target and block later auxiliary echo dataData z of this range cell to be detectedkTake from step successively
Dimensionality reduction echo data z in 1;Data basic matrix φ using blocking matrix b and range cell to be detectedkObtain target block with
Data basic matrix afterwards
Step 3 specifically includes following sub-step:
3a) to space-time steering vector sz-efaTransposed form carry out singular value decomposition, obtain left singular matrix u, right unusual
Matrix v and singular value matrix λ;Decomposition formula such as following formula:
Wherein, u represents left singular matrix, and v represents right singular matrix, and λ represents singular value matrix, ()tRepresent transposition behaviour
Make, ()hRepresent conjugate transposition operation;
3b) blocking matrix b is arranged to make up by the 2nd to the n row order of right singular matrix v, by right singular matrix v and following formula
Derive blocking matrix b;
B=[v (:, 2:n)]t
Wherein, v (:, 2:n) represents all row of right singular matrix v, the 2nd arrives n column matrix element;(·)tRepresent transposition behaviour
Make;
3c) auxiliary echo data is obtained according to following formulaBlock later data basic matrix with target
Wherein, zkRepresent the data of k-th range cell to be detected, the data of this range cell to be detected takes from step successively
Dimensionality reduction echo data z in rapid 1, k represent the sequence number of the data of range cell to be detected, are single less than all distances to be detected
The natural number of unit's sum;Represent that target blocks later auxiliary echo data,Represent that target blocks later data group moment
Battle array, b is blocking matrix.
Step 4, blocks later data basic matrix using targetTo auxiliary echo dataCarry out data matching, obtain
Auxiliary echo dataFitting coefficient
Step 4 specifically includes following sub-step:
4a) block later data basic matrix with targetLater auxiliary echo data is blocked to targetCarry out
A young waiter in a wineshop or an inn takes advantage of matching to represent, is expressed as following formula:
Wherein, | | | |2Represent 2 norms solving vector,Represent auxiliary echo dataFitting coefficient;
4b) according to step 4a) in least square fitting represent, solve auxiliary echo dataFitting coefficient?
A young waiter in a wineshop or an inn takes advantage of solutionIt is expressed as following formula:
Wherein, ()-1Representing matrix inversion operation, ()hRepresent conjugate transposition operation,After representing that target is blocked
Data basic matrix,Represent that target blocks later auxiliary echo data.
Step 5, using the data basic matrix φ of range cell to be detectedkWith auxiliary echo dataFitting coefficientRight
Data z of range cell to be detectedkCarry out data matching, obtain data z of range cell to be detectedkMinimum error of fitting yk;
Data z of this range cell to be detectedkTake from the dimensionality reduction echo data z in step 1 successively;
Step 5 specifically includes following sub-step:
5a) with the data basic matrix φ of range cell to be detectedk, treat data z of detecting distance unitkCarry out a young waiter in a wineshop or an inn
Take advantage of matching;With solution auxiliary echo dataFitting coefficientLeast square solutionIt is calculated optimization restrained boundary ηk, tool
Body is realized by following formula:
Wherein, αkRepresent the fitting coefficient of the data of range cell to be detected, ykThe 2 of the minimum error of fitting vector of expression
Norm, ηkRepresent and optimize restrained boundary
5b) by solving the optimal fitting coefficient that above formula obtains the data of range cell to be detected
The optimal fitting coefficient of the data according to range cell to be detectedObtain the data of range cell to be detected
Little error of fitting yk, obtained by following formula:
Wherein, zkRepresent the data of k-th range cell to be detected, the data of this range cell to be detected takes from step successively
Dimensionality reduction echo data z in rapid 1, k represent the sequence number of the data of range cell to be detected, are single less than all distances to be detected
The natural number of unit's sum;φkRepresent the data basic matrix of range cell to be detected;Represent the data of range cell to be detected
Optimal fitting coefficient;||·||2Represent 2 norms solving vector.
Step 6, treats minimum error of fitting y of the data of detecting distance unitkCarry out CA-CFAR detection, and
There is target or there is not target in output.
Treat minimum error of fitting y of the data of detecting distance unitkCarry out CA-CFAR detection (cell-
Averaging constant false alarm rate, ca-cfar), will minimum error of fitting ykSingle with distance to be detected
The meansigma methodss of the minimum error of fitting of the data of the range cell of surrounding of unit compare, according to minimum error of fitting ykWith meansigma methodss
Ratio size, to determine whether there is target, finally there will be target or there is not target output.
With reference to emulation experiment, the effect of the present invention is described further.
(1) experiment condition:
The experiment of the present invention is carried out under matlab7.11 software.In the experimental design of the present invention, in order to check this
Inventive method effectiveness in actual applications, processing data is derived from the measured data mcarm data set of external admission, its radar
It is operated under positive side battle array pattern, front is made up of the spatial domain array element that can be used for Adaptive Signal Processing that 2 row 11 arrange, its part system
System parameter is with reference to table 1.In order to test conveniently, the present invention extracts 11 array elements of the first row in mcarm data, front 32 pulses
Echo data is verified.
Table 1
(2) Comparison of experiment results
Accompanying drawing 2 is general pulse doppler processing result figure, is also range Doppler figure.Transverse axis represents Doppler's passage sequence
Number, the longitudinal axis represents range cell sequence number.Range cell is also referred to as range gate in the present invention.Although as can be seen from Figure 2 this thunder
Reach in the case of being operated in positive side battle array, but the strong clutter component in space-time two-dimensional echo data, by different the adjusting of landform reflectance
In system and actual environment, the presence of various false targets, has certain heterogeneity.
Accompanying drawing 3 is the change curve with border compensation dosage for the detection probability.Transverse axis represents boundary compensation, and the longitudinal axis represents detection
Probability.Boundary compensation is that the optimization upper bound showing gained when optimizing restrained boundary in solution procedure 5 adds a compensation dosage.From
As can be seen that substantially near 0 compensation dosage, detection probability can obtain maximum in Fig. 3, this also just illustrates step in the present invention
Rapid 5 calculation optimization restrained boundaries are effective and optimum.
Accompanying drawing 4 is the change curve with sample number for the detection probability.Sample number refers to close on range cell in the present invention
The number of data.Transverse axis represents sample number, and the longitudinal axis represents detection probability.In the present invention after the dimension-reduction treatment of data efa, data
Dimension (i.e. degree of freedom) is reduced to 3n, and wherein n is array number, 3n=33 in this experiment.It can be seen that detection probability exists
Maximum is basically reached, this just illustrates, in the present invention, data basic matrix can be attached by range cell to be detected when sample number is more than 100
Nearly 3 times of number of degrees of freedom, purposes close on the data of range cell to constitute, and obtain the detection performance of optimum.
Fig. 5 is the output result figure of the minimum error of fitting of the data of target proximity range cell.In the present invention, target
Nearby the data of range cell includes the data of multigroup range cell to be detected.Transverse axis represents the sequence number of range cell, longitudinal axis table
Show power, represented with db.In Fig. 5, curve represents the minimum error of fitting of a segment distance unit scope in target Doppler passage
Output result, horizontal linear represents the average output value of the minimum error of fitting of this segment distance unit scope, and circle represents target
The minimum error of fitting output valve of range cell, in Fig. 5, " difference " represents the minimum error of fitting output valve of target range unit
The difference and average output value of minimum error of fitting of this segment distance unit scope between.From figure 5 it can be seen that through this
After the reason of daylight, target output becomes apparent from, and Fig. 5 (a) power raising when range cell is 150 compared with 5 (c) exceedes
15db, Fig. 5 (b) power when range cell is 200 compared with 5 (d) improves more than 25db, such that it is able to be more beneficial for target
The detection of signal.
Fig. 6 is the change curve with target power for the detection probability.Transverse axis represents the target power of artificial setting, longitudinal axis table
Show the detection probability of this target.The sparse restoration methods of direct data domain (d3sr), conventional Extension factorization method is compared in Fig. 6
And the spreading factor method (df-efa) based on data matching proposed by the present invention (efa).From fig. 6 it can be seen that this
The detection performance of bright method is much better than other two methods.
Fig. 7 is the change curve graph of a relation with false alarm rate for the detection probability.Transverse axis represents false alarm rate, and the longitudinal axis represents that detection is general
Rate.Conventional Extension factorization method (efa) and the spreading factor side based on data matching proposed by the present invention is compared in Fig. 7
Method (df-efa).It can be seen from figure 7 that the inventive method, with respect to traditional efa method, has higher under certain false alarm rate
Detection probability, and there is under identical detection probability lower false alarm rate.
Claims (6)
1. a kind of spreading factor space-time adaptive processing method based on data matching is it is characterised in that include following
Step:
Step 1, receives the space-time two-dimensional echo data x of airborne early warning radar, to space-time two-dimensional using airborne early warning radar antenna
Echo data x is extended factorization efa dimensionality reduction, obtains dimensionality reduction echo data z and space-time steering vector sz-efa;
Step 2, constructs the data basic matrix φ of range cell to be detectedk;
Step 3, using space-time steering vector sz-efaDerive blocking matrix b;Number using blocking matrix b and range cell to be detected
According to zk, obtain target and block later auxiliary echo dataData z of this range cell to be detectedkTake from successively in step 1
Dimensionality reduction echo data z;Data basic matrix φ using blocking matrix b and range cell to be detectedkAfter obtaining target obstruction
Data basic matrix
Step 4, blocks later data basic matrix using targetTo auxiliary echo dataCarry out data matching, obtain auxiliary
Echo dataFitting coefficient
Step 5, using the data basic matrix φ of range cell to be detectedkWith auxiliary echo dataFitting coefficientTo be checked
Survey data z of range cellkCarry out data matching, obtain data z of range cell to be detectedkMinimum error of fitting yk;This is treated
Data z of detecting distance unitkTake from the dimensionality reduction echo data z in step 1 successively;
Step 6, treats minimum error of fitting y of the data of detecting distance unitkCarry out CA-CFAR detection, and export
There is target or there is not target.
2. the spreading factor space-time adaptive processing method based on data matching according to claim 1, it is special
Levy and be, step 1 specifically includes following sub-step:
1a) utilize airborne early warning radar antenna, receive the space-time two-dimensional echo data x of ground return within the coherent accumulation time;
Space-time two-dimensional echo data x ties up for mn, and wherein m represents umber of pulse, and n represents the array number of airborne early warning radar antenna;
1b) be expanded the dimensionality reduction transition matrix p of factorization efa dimensionality reduction according to following formulaefa:
Wherein, pt(fi-1,fi,fi+1) represent the time domain transition matrix being made up of three Doppler filters closing on, fiRepresent mesh
Target normalization Doppler frequency, fi-1,fi+1Represent two normalization Doppler frequencies adjacent with target, i be less than or equal to
M and the natural number more than or equal to 1, m represents umber of pulse, inRepresent the unit matrix of n × n,Represent Kronecker product, n represents
The array number of airborne early warning radar antenna, dimensionality reduction transition matrix pefaDimension be mn × 3n tie up;
1c) according to dimensionality reduction transition matrix pefaBe expanded the dimensionality reduction echo data z after factorization efa dimensionality reduction, is expressed as following formula:
Wherein, ()hRepresent conjugate transposition operation, according to dimensionality reduction transition matrix pefaDimension be mn × 3n dimension and space-time two-dimensional
The dimension that echo data x can obtain dimensionality reduction echo data z for mn dimension is tieed up for 3n;
Space-time steering vector s after spreading factor efa dimensionality reductionz-efa, it is expressed as following formula:
ss(fs)=[1 exp (j2 π fs) exp(j2πfs2) … exp(j2πfs(n-1))]t
Wherein, ()tRepresent transposition operation, ss(fs) represent target spatial domain steering vector, fsRepresent the normalization space of target
Frequency, n represents the array number of airborne early warning radar antenna.
3. the spreading factor space-time adaptive processing method based on data matching according to claim 1, it is special
Levy and be, step 2 specifically includes:
Data basic matrix φkIt is expressed as following formula:
φk=[z1z2…zl]∈cd×l
Wherein, φkRepresent the data basic matrix of k-th range cell to be detected, k represents the sequence of the data of range cell to be detected
Number, it is less than the natural number of all range cell sums to be detected;D represents the dimension of dimensionality reduction echo data z, and l represents and closes on
The number of the data of range cell, cd×lRepresent the complex matrix space of d × l dimension, zl(l=1,2 ..., l) represent and close on distance list
The data of unit, closes on range cell and represents range cell near range cell to be detected, do not comprise two and be used for protecting target
Range cell.
4. the spreading factor space-time adaptive processing method based on data matching according to claim 1, it is special
Levy and be, step 3 specifically includes following sub-step:
3a) to space-time steering vector sz-efaTransposed form carry out singular value decomposition, obtain left singular matrix u, right singular matrix
V and singular value matrix λ;Decomposition formula such as following formula:
Wherein, u represents left singular matrix, and v represents right singular matrix, and λ represents singular value matrix, ()tRepresent transposition operation,
(·)hRepresent conjugate transposition operation;
3b) blocking matrix b is arranged to make up by the 2nd to the n row order of right singular matrix v, is derived by right singular matrix v and following formula
Blocking matrix b;
B=[v (:, 2:n)]t
Wherein, v (:, 2:n) represents all row of right singular matrix v, the 2nd arrives n column matrix element;(·)tRepresent transposition operation;
3c) auxiliary echo data is obtained according to following formulaBlock later data basic matrix with target
Wherein, zkRepresent the data of k-th range cell to be detected, the data of this range cell to be detected is taken from step 1 successively
Dimensionality reduction echo data z, k represents the sequence number of the data of range cell to be detected, is total less than all range cells to be detected
The natural number of number;Represent that target blocks later auxiliary echo data,Represent that target blocks later data basic matrix, b
For blocking matrix.
5. the spreading factor space-time adaptive processing method based on data matching according to claim 1, it is special
Levy and be, step 4 specifically includes following sub-step:
4a) block later data basic matrix with targetLater auxiliary echo data is blocked to targetCarry out a young waiter in a wineshop or an inn
Take advantage of matching to represent, be expressed as following formula:
Wherein, | | | |2Represent 2 norms solving vector,Represent auxiliary echo dataFitting coefficient;
4b) according to step 4a) in least square fitting represent, solve auxiliary echo dataFitting coefficientA young waiter in a wineshop or an inn
Take advantage of solutionIt is expressed as following formula:
Wherein, ()-1Representing matrix inversion operation, ()hRepresent conjugate transposition operationRepresent that target blocks later data
Basic matrix,Represent that target blocks later auxiliary echo data.
6. the spreading factor space-time adaptive processing method based on data matching according to claim 1, it is special
Levy and be, step 5 specifically includes following sub-step:
5a) with the data basic matrix φ of range cell to be detectedk, treat data z of detecting distance unitkCarry out least square plan
Close;With solution auxiliary echo dataFitting coefficientLeast square solutionIt is calculated optimization restrained boundary ηk, specifically lead to
Cross following formula to realize:
Wherein, αkRepresent the fitting coefficient of the data of range cell to be detected, ykRepresent 2 norms of minimum error of fitting vector,
ηkRepresent and optimize restrained boundary,
5b) by solving the optimal fitting coefficient that above formula obtains the data of range cell to be detected
The optimal fitting coefficient of the data according to range cell to be detectedObtain the minimum matching of the data of range cell to be detected
Error yk, obtained by following formula:
Wherein, zkRepresent the data of k-th range cell to be detected, the data of this range cell to be detected is taken from step 1 successively
Dimensionality reduction echo data z, k represents the sequence number of the data of range cell to be detected, is total less than all range cells to be detected
The natural number of number;φkRepresent the data basic matrix of range cell to be detected;Represent the optimum of the data of range cell to be detected
Fitting coefficient;||·||2Represent 2 norms solving vector.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410276078.6A CN104459657B (en) | 2014-06-19 | 2014-06-19 | Extension factorization space-time two-dimensional self-adaptive processing method based on data fitting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410276078.6A CN104459657B (en) | 2014-06-19 | 2014-06-19 | Extension factorization space-time two-dimensional self-adaptive processing method based on data fitting |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104459657A CN104459657A (en) | 2015-03-25 |
CN104459657B true CN104459657B (en) | 2017-01-25 |
Family
ID=52906010
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410276078.6A Expired - Fee Related CN104459657B (en) | 2014-06-19 | 2014-06-19 | Extension factorization space-time two-dimensional self-adaptive processing method based on data fitting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104459657B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104950296B (en) * | 2015-06-12 | 2017-06-20 | 西安电子科技大学 | Based on the remaining sane non-homogeneous detection method of weight Weighted adaptive power |
CN105093189B (en) * | 2015-09-08 | 2017-06-20 | 西安电子科技大学 | Airborne radar object detection method based on GCV |
CN105911527B (en) * | 2016-06-30 | 2018-03-06 | 西安电子科技大学 | Airborne radar space-time adaptive processing method based on EFA and MWF |
CN106546966B (en) * | 2016-10-31 | 2019-03-29 | 西安电子科技大学 | Radar noise power estimation method under clutter background based on fitting of a polynomial |
CN109212502B (en) * | 2018-09-10 | 2022-12-27 | 中国航天科工集团八五一一研究所 | Method for realizing sampling matrix inversion non-uniform detector based on reweighting |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5748143A (en) * | 1996-12-09 | 1998-05-05 | The United States Of America As Represented By The Secretary Of The Air Force | Adaptive post-doppler sequential beam processor |
CN101819269A (en) * | 2010-03-19 | 2010-09-01 | 清华大学 | Space-time adaptive processing method under non-homogeneous clutter environment |
CN103399303A (en) * | 2013-07-22 | 2013-11-20 | 西安电子科技大学 | Method and system for resisting dense forwarding type defraud interference of airborne radar |
CN103399309A (en) * | 2013-07-22 | 2013-11-20 | 西安电子科技大学 | Space-time two-dimensional clutter spectrum estimation method based on iterative weighted minimum variance |
CN103439692A (en) * | 2013-09-01 | 2013-12-11 | 西安电子科技大学 | STAP method based on wide symmetrical characteristic of covariance matrix |
CN103728606A (en) * | 2014-01-16 | 2014-04-16 | 西安电子科技大学 | Doppler channel correlation two-stage dimension reduction method for onboard multiple input multiple output (MIMO) radar |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6801161B2 (en) * | 2001-12-05 | 2004-10-05 | Lockheed Martin Corporation | System and method for auto calibrated reduced rank adaptive processor |
-
2014
- 2014-06-19 CN CN201410276078.6A patent/CN104459657B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5748143A (en) * | 1996-12-09 | 1998-05-05 | The United States Of America As Represented By The Secretary Of The Air Force | Adaptive post-doppler sequential beam processor |
CN101819269A (en) * | 2010-03-19 | 2010-09-01 | 清华大学 | Space-time adaptive processing method under non-homogeneous clutter environment |
CN103399303A (en) * | 2013-07-22 | 2013-11-20 | 西安电子科技大学 | Method and system for resisting dense forwarding type defraud interference of airborne radar |
CN103399309A (en) * | 2013-07-22 | 2013-11-20 | 西安电子科技大学 | Space-time two-dimensional clutter spectrum estimation method based on iterative weighted minimum variance |
CN103439692A (en) * | 2013-09-01 | 2013-12-11 | 西安电子科技大学 | STAP method based on wide symmetrical characteristic of covariance matrix |
CN103728606A (en) * | 2014-01-16 | 2014-04-16 | 西安电子科技大学 | Doppler channel correlation two-stage dimension reduction method for onboard multiple input multiple output (MIMO) radar |
Non-Patent Citations (6)
Title |
---|
"Extended Factored Space-Time Processing for Airborne Radar";Robert C. DiPietro;《Proceedings of the 26th Asilomar conference on signals,systems and computing》;19921231;第425-430页 * |
"SAR-MTI Processing of Multi-Channel Airborne Radar Measurement (MCARM)";Mehrdad Soumekh,等;《Radar conference,2002,Proceeding of the IEEE》;20021231;第24-48页 * |
"Space-Time Adaptive Radar";WILLIAM L. MELVIN;《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》;20000831;第36卷(第2期);第621-633页 * |
"基于非线性最小二乘的空中机动目标检测方法";李海,等;《电子与信息学报》;20120930;第34卷(第9期);第2143-2148页 * |
"机载雷达空时二维信号处理自适应权值";廖桂生,等;《信号处理》;19980930;第14卷(第3期);第233-238页 * |
"机载雷达简易STAP方法及其应用";王彤;《中国优秀博硕士学位论文全文数据库(博士) 信科技辑》;20020615(第01(2012)期);第I136-109页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104459657A (en) | 2015-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104459658B (en) | Unite-domain positioning space-time two-dimensional self-adaptive processing method based on data fitting | |
CN104459657B (en) | Extension factorization space-time two-dimensional self-adaptive processing method based on data fitting | |
CN103439692B (en) | STAP method based on wide symmetrical characteristic of covariance matrix | |
CN103364764B (en) | Airborne radar non-stationary clutter suppression method | |
CN102147464B (en) | Intelligent system and method for forecasting robust radar sea clutter | |
CN103454624B (en) | The direct data domain moving target detection method of spectrum time empty based on dimensionality reduction sparse reconstruct | |
CN104459660B (en) | Array element Doppler domain space-time two-dimensional self-adaptive processing method based on data fitting | |
CN103942449B (en) | Feature interference cancellation beam forming method based on estimation of number of information sources | |
CN106529428A (en) | Underwater target recognition method based on deep learning | |
CN106468770A (en) | Closely optimum radar target detection method under K Distribution Clutter plus noise | |
CN107390194B (en) | A kind of radar target detection method based on the graceful divergence of full Donald Bragg | |
CN104535973A (en) | Target detection method by use of airborne early warning radar | |
CN102866388B (en) | Iterative computation method for self-adaptive weight number in space time adaptive processing (STAP) | |
CN107085205A (en) | Self-adapting detecting method based on clutter covariance matrix structural information | |
CN103176168A (en) | Short-range cluster cancellation method for airborne non-side-looking array radar | |
CN105204007A (en) | Range-spread target polarization adaptive detection method without auxiliary data | |
CN105223554A (en) | Based on the space-time adaptive Monopulse estimation method of Doppler's triple channel Combined Treatment | |
CN105158749A (en) | High-frequency radar sea-clutter amplitude statistical distribution test method | |
CN108229750A (en) | A kind of stock yield Forecasting Methodology | |
CN106202491A (en) | A kind of based on seasonal effect in time series data digging method and system | |
CN103413296B (en) | To polarization SAR order object detection method in property field | |
CN104215939A (en) | Knowledge assisted space-time adaptive processing method integrating generalized symmetrical structure information | |
CN105372635A (en) | Improved dimension-reduction space-time adaptive processing-based ship-borne high-frequency ground wave radar sea clutter suppression method | |
CN102183745B (en) | Sea clutter forecasting system and method for intelligent radar | |
CN104459659B (en) | Wave beam Doppler domain space-time two-dimensional self-adaptive processing method based on data fitting |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
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: 20170125 Termination date: 20170619 |
|
CF01 | Termination of patent right due to non-payment of annual fee |