CN104459659B - Wave beam Doppler domain space-time two-dimensional self-adaptive processing method based on data fitting - Google Patents
Wave beam Doppler domain space-time two-dimensional self-adaptive processing method based on data fitting Download PDFInfo
- Publication number
- CN104459659B CN104459659B CN201410276594.9A CN201410276594A CN104459659B CN 104459659 B CN104459659 B CN 104459659B CN 201410276594 A CN201410276594 A CN 201410276594A CN 104459659 B CN104459659 B CN 104459659B
- Authority
- CN
- China
- Prior art keywords
- data
- detected
- represent
- 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 a wave beam Doppler domain space-time two-dimensional self-adaptive processing method based on data fitting, and relates to radar technology. The method comprises the steps that 1, wave beam Doppler domain 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 optimal 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 wave beam Doppler domain 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, improves movable target detection probability, reduces detection false alarm rate, and has the good application prospect in measured data processing and STAP performance improvement.
Description
Technical field
The invention belongs to communication technical field, it is related to Radar Technology, many particularly to a kind of wave beam based on data matching
General Le domain 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. cn 101819269 a) disclose a kind of in non-homogeneous clutter environment super-resolution estimate miscellaneous
The overcomplete sparse representation method of ripple space-time two-dimensional spectrum.The method achieve in the case of independent same distribution sample number deficiency,
Estimate clutter covariance matrix using single frames training sample, thus avoiding strong non-homogeneous clutter environment to self-adaptive processing effect
Impact.But, the method yet suffers from following main deficiency is: one is the big problem of operand, carries out sparse table to clutter spectrum
The super complete radix mesh showing 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 exists
In the covariance matrix restructuring procedure of each range cell sample, required operand is very big, is unfavorable for real-time processing, from
And have influence on the effect in practical engineering application.Two is base mismatch problems, and used in the method, actually one group of base is inserted
Value discrete fourier dft vector, and the echo data of reality because there is the non-ideal factors such as error it is impossible to basic vector
Amount rarefaction representation;In addition, the method needs to carry out discretization expression to space-time plane, such that not being located at discrete grid block point
The signal at center 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 wave beam Doppler domain space-time adaptive processing method based on data matching, energy
Enough detections realized to echo signal, and solve the problems, such as in above-mentioned sparse recoverys stap method operand huge with base mismatch,
Simultaneously with respect to traditional wave beam Doppler domain stap method, 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 wave beam Doppler domain space-time adaptive processing method based on data matching is it is characterised in that include
Following steps:
Step 1, receives the space-time two-dimensional echo data x of airborne early warning radar, to space-time using mechanical early warning radar antenna
Two-dimentional echo data x carries out joint domain positioning jdl dimensionality reduction, obtains dimensionality reduction echo data z and the space-time guiding of wave beam Doppler domain
Vector sz-jdl;
Step 2, constructs the data basic matrix φ of range cell to be detectedk;
Step 3, using the dimensionality reduction space-time steering vector s of wave beam Doppler domainz-jdlDerive blocking matrix b;Using obstruction square
Battle array b and data z of range cell to be detectedk, obtain target and block later auxiliary echo dataThis range cell to be detected
Data zkTake from the dimensionality reduction echo data z of the wave beam Doppler domain obtaining in step 1 successively;Using blocking matrix b and to be checked
Survey the data basic matrix φ of range cellkObtain target and block later data basic matrix
Step 4, blocks later data basic matrix using targetTo auxiliary echo dataCarry out data matching, obtain
Take auxiliary echo dataOptimal fitting coefficient
Step 5, using the data basic matrix φ of range cell to be detectedkWith auxiliary echo dataOptimal fitting coefficientTreat detecting distance cell data zkCarry out data matching, obtain minimum error of fitting y of the data of range cell to be detectedk;
Data z of this range cell to be detectedkTake from the dimensionality reduction echo data z of the wave beam Doppler domain obtaining 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, this 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) the dimensionality reduction transition matrix p that joint domain positions jdl dimensionality reduction is obtained according to following formulajdl:
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 the 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, ps(gj-1,gj,gj+1) represent be made up of three spatial filters
Spatial domain transition matrix, the wave filter in the middle of these three spatial filters is the spatial domain steering vector of target, gjReturning for target
One change spatial frequency, gj-1,gj+1Be neighbouring two normalization spatial frequency, j be less than or equal to m and more than or equal to 1 from
So count;Represent Kronecker product;Dimensionality reduction transition matrix pjdlDimension be mn × 9 tie up;
1c) according to dimensionality reduction transition matrix pjdlObtain the dimensionality reduction echo data z of wave beam Doppler domain, be expressed as following formula:
Wherein, ()hRepresent conjugate transposition operation, according to dimensionality reduction transition matrix pjdlDimension be mn × 9 dimension and space-time
The dimension that two-dimentional echo data x can obtain dimensionality reduction echo data z for mn dimension is 9 dimensions.
The dimensionality reduction space-time steering vector s of wave beam Doppler domainz-jdl, it is expressed as following formula:
Wherein, ()tRepresent transposition operation, sz-jdlFor the space-time steering vector after the dimensionality reduction of wave beam Doppler domain.
(2) 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 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 dimensionality reduction echo data of wave beam Doppler domain
The dimension of z, l represents the number of the data closing on range cell, cd×lRepresent the complex matrix space of d × l dimension, zl(l=1,2 ...,
L) represent and close on the data of range cell, close on range cell and represent range cell near range cell to be detected, do not comprise
Two are used for protecting the range cell of target.
(3) step 3 specifically includes following sub-step:
3a) the dimensionality reduction space-time steering vector s to wave beam Doppler domainz-jdlTransposed 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 that left singular matrix, v represent that right singular matrix, λ represent singular value matrix;(·)tRepresent transposition behaviour
Make, ()hRepresent conjugate transposition operation;
3b) blocking matrix b is arranged to make up by the 2nd to n row order of right singular matrix, 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 blocking matrix bBlock later data basic matrix with targetAs follows
Shown in 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
The dimensionality reduction echo data z of the wave beam Doppler domain obtaining in rapid 1, k represents the sequence number of the data of range cell to be detected, for not surpassing
Cross the natural number of all range cell sums to be detected;Represent that target blocks later auxiliary echo data,Represent target
Block later data basic matrix, 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) pass through to solve above formula optimization problem, obtain assisting echo dataOptimal fitting coefficient
(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;And with assisting echo dataOptimal fitting coefficientIt is calculated optimization restrained boundary ηk, especially by following formula
Least square fitting formula is realized:
Wherein, αkRepresent the fitting coefficient of the data of range cell to be detected, ηkRepresent and optimize restrained boundary,
5b) by solving the optimal fitting coefficient that above formula least square fitting 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
The dimensionality reduction echo data z of the wave beam Doppler domain obtaining in rapid 1, k represents the sequence number of the data of range cell to be detected, for not surpassing
Cross the natural number of all range cell sums to be detected;φkRepresent the data basic matrix of range cell to be detected;Represent to be checked
Survey the optimal fitting coefficient of the data of range cell;||·||2Represent 2 norms solving vector.
Compared with prior art, the present invention has prominent substantive distinguishing features and significantly improves.The present invention with existing
Stap method is compared, and has the advantage that
(1) it is directed to the big problem of operand in existing sparse recovery stap method, sparse restoration methods are carried out to clutter spectrum
The super complete radix mesh of rarefaction representation 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 required in the covariance matrix restructuring procedure of each range cell sample operand is very big, is unfavorable for locating in real time
Reason, thus have influence on its practical engineering application effect.Institute's extracting method is positioned (joint domain with combining domain by the present invention
Localized, jdl) method combines, and the little optimization problem that the optimization problem of original mn dimension is decomposed into m 9 dimension to solve,
Greatly reduce the computation complexity of algorithm, solve the problems, such as that in sparse recovery stap method, operand is big.
(2) it is directed to base mismatch problems in existing sparse recovery stap method, used in sparse restoration methods, represent that base is real
It is one group of interpolation dft vector on border, and real data is because have the non-ideal factors such as error it is impossible to sparse with given base vector
Represent, in addition, the method needs to carry out discretization expression to space-time plane, such that not being located at discrete grid block dot center
Signal is revealed to all mesh points, thus destroying the openness of data.The present invention is directly used real data sample as base
Matrix, rather than the ideal flowing pattern being determined using certain, thus evaded asking of base mismatch in above-mentioned sparse recovery stap method
Topic.
(3) traditional wave beam Doppler domain stap method, such as joint domain localization method (jdl), by choosing distance to be detected
Data sample around unit, to estimate the statistical property of range cell data to be detected, to carry out clutter and AF panel.So
And when not containing interference data to be suppressed in the data sample chosen, it is present in the interference letter in range cell data to be detected
Number cannot suppress, form false-alarm.The inventive method on the basis of blocking operation is carried out to echo signal, to pending number
According to carrying out abundant matching expression, so that main lobe echo signal is retained, and (formula is also referred to as cheated to secondary lobe echo signal
Interference) carry out abundant suppression.
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 required training sample number, in order to obtain better performance, the present invention is directed to above-mentioned asking
Topic, closes on the data sample of range cell come Optimal Fitting range cell to be detected clutter data using range cell to be detected,
Thus improving the clutter recognition performance of Adaptive Signal Processing, improve 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 the distance-Doppler result figure after general pulse doppler processing;
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 result figure of the data of partial distance-doppler cells;Fig. 5 (a) is conventional wave beam Doppler domain
Distance-Doppler result figure after the process of stap algorithm jdl;Fig. 5 (b) is the distance-Doppler after the inventive method process
Result figure;
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 wave beam Doppler domain 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 mechanical early warning radar antenna
Two-dimentional echo data x carries out joint domain positioning jdl dimensionality reduction, obtains dimensionality reduction echo data z and the space-time guiding of wave beam Doppler domain
Vector sz-jdl;
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, this 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) the dimensionality reduction transition matrix p that joint domain positions jdl dimensionality reduction is obtained according to following formulajdl:
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 the 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, ps(gj-1,gj,gj+1) represent be made up of three spatial filters
Spatial domain transition matrix, the wave filter in the middle of these three spatial filters is the spatial domain steering vector of target, gjReturning for target
One change spatial frequency, gj-1,gj+1Be neighbouring two normalization spatial frequency, j be less than or equal to m and more than or equal to 1 from
So count;Represent Kronecker product;Dimensionality reduction transition matrix pjdlDimension be mn × 9 tie up.
1c) according to dimensionality reduction transition matrix pjdlObtain the dimensionality reduction echo data z of wave beam Doppler domain, be expressed as following formula:
Wherein, ()hRepresent conjugate transposition operation, according to dimensionality reduction transition matrix pjdlDimension be mn × 9 dimension and space-time
The dimension that two-dimentional echo data x can obtain dimensionality reduction echo data z for mn dimension is 9 dimensions.
The dimensionality reduction space-time steering vector s of wave beam Doppler domainz-jdl, it is expressed as following formula:
Wherein, ()tRepresent transposition operation, sz-jdlFor the space-time steering vector after the dimensionality reduction of wave beam Doppler domain.
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 respectively.Time domain and spatial domain steering vector
It 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 dimensionality reduction echo data of wave beam Doppler domain
The dimension of z, l represents the number of the data closing on range cell, cd×lRepresent the complex matrix space of d × l dimension, zl(l=1,2 ...,
L) represent and close on the data of range cell, close on range cell and represent range cell near range cell to be detected, do not comprise
Two are used for protecting the range cell of target.
Step 3, using the dimensionality reduction space-time steering vector s of wave beam Doppler domainz-jdlDerive blocking matrix b;Using obstruction square
Battle array b and data z of range cell to be detectedk, obtain target and block later auxiliary echo dataThis range cell to be detected
Data zkTake from the dimensionality reduction echo data z of the wave beam Doppler domain obtaining in step 1 successively;Using blocking matrix b and to be checked
Survey the data basic matrix φ of range cellkObtain target and block later data basic matrix
Step 3 specifically includes following sub-step:
3a) the dimensionality reduction space-time steering vector s to wave beam Doppler domainz-jdlTransposed 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 that left singular matrix, v represent that right singular matrix, λ represent singular value matrix;(·)tRepresent transposition behaviour
Make, ()hRepresent conjugate transposition operation;
3b) blocking matrix b is arranged to make up by the 2nd to n row order of right singular matrix, 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 blocking matrix bBlock later data basic matrix with targetAs follows
Shown in 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
The dimensionality reduction echo data z of the wave beam Doppler domain obtaining in rapid 1, k represents the sequence number of the data of range cell to be detected, for not surpassing
Cross the natural number of all range cell sums to be detected;Represent that target blocks later auxiliary echo data,Represent target
Block later data basic matrix, b is blocking matrix.
Step 4, blocks later data basic matrix using targetTo auxiliary echo dataCarry out data matching, obtain
Take auxiliary echo dataOptimal fitting coefficient
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) pass through to solve above formula optimization problem, obtain assisting echo dataOptimal fitting coefficient
Step 5, using the data basic matrix φ of range cell to be detectedkWith auxiliary echo dataOptimal fitting coefficientTreat detecting distance cell data zkCarry out data matching, obtain minimum error of fitting y of the data of range cell to be detectedk;
Data z of this range cell to be detectedkTake from the dimensionality reduction echo data z of the wave beam Doppler domain obtaining 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;And with assisting echo dataOptimal fitting coefficientIt is calculated optimization restrained boundary ηk, especially by following formula
Least square fitting formula is realized:
Wherein, αkRepresent the fitting coefficient of the data of range cell to be detected, ηkRepresent and optimize restrained boundary,
5b) by solving the optimal fitting coefficient that above formula least square fitting 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
The dimensionality reduction echo data z of the wave beam Doppler domain obtaining in rapid 1, k represents the sequence number of the data of range cell to be detected, for not surpassing
Cross the natural number of all range cell sums to be detected;φkRepresent the data basic matrix of range cell to be detected;Represent to be checked
Survey the optimal fitting coefficient of the data of range cell;||·||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 matlab 7.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.Transverse axis represents Doppler's channel position, and the longitudinal axis represents that distance is single
First sequence number.Range cell is also referred to as range gate in the present invention.Although as can be seen from Figure 2 this radar is operated in positive side battle array situation
Under, but the strong clutter component in space-time two-dimensional echo data, by various in the different modulation of landform reflectance and actual environment
The presence of false target, 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.When boundary compensation refers to optimize restrained boundary in solution procedure 5, the optimization upper bound of gained adds a compensation dosage.From Fig. 3
In as can be seen that basic near 0 compensation dosage, detection probability can obtain maximum, and this also just illustrates step 5 in the present invention
Calculation optimization restrained boundary is 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 jdl, data
Dimension (i.e. degree of freedom) is reduced to 9.Figure 4, it is seen that detection probability basically reaches maximum when sample number is more than 36, this
Just illustrate, in the present invention, data basic matrix can nearby 4 times of number of degrees of freedom, purposes close on range cell by range cell to be detected
Data, to constitute, obtains the detection performance of optimum.
Accompanying drawing 5 is the result figure of the data of partial distance-doppler cells.In the present invention, partial distance-how general
The data strangling unit includes the data of multigroup range cell to be detected.Transverse axis represents Doppler's passage, from 17 to 32;The longitudinal axis represents
Range gate sequence number, from 140 to 240.Range cell is also referred to as range gate in the present invention.In this experiment, two main lobe target (figures
In represented with circle) and four secondary lobe targets (in figure is represented with rhombus) be artificially added to respectively in experimental data.We need
To be detected is main lobe target, and secondary lobe target is Deceiving interference, shows as false-alarm.Figure (a) is that how general conventional wave beam is
Strangle the result figure of domain stap algorithm jdl, it can be seen that main lobe target can not be evident from, and secondary lobe target is not
Can be inhibited, false-alarm occurs.Figure (b) is the output result after the inventive method process, it is seen that secondary lobe target letter
Number obtain fine suppression, main lobe target also can be detected well.
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), classical joint domain localization method is compared in Fig. 6
And joint domain localization method (df-jdl) based on data matching proposed by the present invention (jdl).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.Classical joint domain localization method (jdl) and the joint domain positioning side based on data matching proposed by the present invention is compared in Fig. 7
Method (df-jdl).It can be seen from figure 7 that the inventive method, with respect to traditional jdl method, has higher under certain false alarm rate
Detection probability, and there is under identical detection probability lower false alarm rate.
Claims (5)
1. a kind of wave beam Doppler domain 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 two-dimensional using airborne early warning radar antenna
Echo data x carries out joint domain positioning jdl dimensionality reduction, obtains dimensionality reduction echo data z and the space-time steering vector of wave beam Doppler domain
sz-jdl;
Step 2, constructs the data basic matrix φ of range cell to be detectedk;
Step 3, using the dimensionality reduction space-time steering vector s of wave beam Doppler domainz-jdlDerive blocking matrix b;Using blocking matrix b and
Data z of range cell to be detectedk, obtain target and block later auxiliary echo dataThe data of this range cell to be detected
zkTake from the dimensionality reduction echo data z of the wave beam Doppler domain obtaining in step 1 successively;Using blocking matrix b and distance to be detected
The data basic matrix φ of unitkObtain target and block later data basic matrix
Step 3 specifically includes following sub-step:
3a) the dimensionality reduction space-time steering vector s to wave beam Doppler domainz-jdlTransposed form carry out singular value decomposition, obtain left strange
Different matrix u, right singular matrix v and singular value matrix λ;Decomposition formula such as following formula:
Wherein, u represents that left singular matrix, v represent that right singular matrix, a represent singular value matrix;()tRepresent transposition operation,
()hRepresent conjugate transposition operation;
3b) blocking matrix b is arranged to make up by the 2nd to n row order of right singular matrix, and by right singular matrix v, and following formula is derived
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 blocking matrix bBlock later data basic matrix with targetIt is shown below:
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
The dimensionality reduction echo data z of the wave beam Doppler domain obtaining, k represents the sequence number of the data of range cell to be detected, is less than institute
Need the natural number of detecting distance unit sum;Represent that target blocks later auxiliary echo data,Represent that target is blocked
Later data basic matrix, b is blocking matrix;
Step 4, blocks later data basic matrix using targetTo auxiliary echo dataCarry out data matching, obtain auxiliary
Echo dataOptimal fitting coefficient
Step 5, using the data basic matrix φ of range cell to be detectedkWith auxiliary echo dataOptimal fitting coefficientRight
Range cell data z to be detectedkCarry out data matching, obtain minimum error of fitting y of the data of range cell to be detectedk;This is treated
Data z of detecting distance unitkTake from the dimensionality reduction echo data z of the wave beam Doppler domain obtaining 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 wave beam Doppler domain space-time adaptive processing method based on data matching according to claim 1, its
It is characterised by, 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,
This 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) the dimensionality reduction transition matrix p that joint domain positions jdl dimensionality reduction is obtained according to following formulajdl:
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 the two normalization Doppler frequencies adjacent with target;I is less than or equal to m
And the natural number more than or equal to 1, m represents umber of pulse, ps(gj-1, gj, gj+1) represent the spatial domain being made up of three spatial filters
Transition matrix, the wave filter in the middle of these three spatial filters is the spatial domain steering vector of target, gjEmpty for the normalization of target
Between frequency, gj-1, gj+1Be neighbouring two normalization spatial frequency, j be less than or equal to m and more than or equal to 1 natural number;
Represent Kronecker product;Dimensionality reduction transition matrix pjdlDimension be mn × 9 tie up;
1c) according to dimensionality reduction transition matrix pjdlObtain the dimensionality reduction echo data z of wave beam Doppler domain, be expressed as following formula:
Wherein, ()hRepresent conjugate transposition operation, according to dimensionality reduction transition matrix pjdlDimension be mn × 9 dimension and space-time two-dimensional
The dimension that echo data x can obtain dimensionality reduction echo data z for mn dimension is 9 dimensions;
The dimensionality reduction space-time steering vector s of wave beam Doppler domainz-jdl, it is expressed as following formula:
Wherein, () t represents that transposition operates, sz-jdlFor the space-time steering vector after the dimensionality reduction of wave beam Doppler domain.
3. the wave beam Doppler domain space-time adaptive processing method based on data matching according to claim 1, its
It is characterised by, 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 dimensionality reduction echo data z's of wave beam Doppler domain
Dimension, l represents the number of the data closing on range cell, cd×lRepresent the complex matrix space of d × l dimension, zl(l=1,2 ..., l)
Represent and close on the data of range cell, close on range cell and represent range cell near range cell to be detected, do not comprise two
The individual range cell for protecting target.
4. the wave beam Doppler domain space-time adaptive processing method based on data matching according to claim 1, its
It is characterised by, step 4 specifically includes following sub-step:
4a) block later data basic matrix with targetLater auxiliary echo data is blocked to targetCarry out least square
Matching represents, is expressed as following formula:
Wherein, | | | |2Represent 2 norms solving vector,Represent auxiliary echo dataFitting coefficient;
4b) pass through to solve above formula optimization problem, obtain assisting echo dataOptimal fitting coefficient
5. the wave beam Doppler domain space-time adaptive processing method based on data matching according to claim 1, its
It is characterised by, 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;And with assisting echo dataOptimal fitting coefficientIt is calculated optimization restrained boundary ηk, especially by a following formula young waiter in a wineshop or an inn
Fitting formula is taken advantage of to realize:
Wherein, αkRepresent the fitting coefficient of the data of range cell to be detected, ηkRepresent and optimize restrained boundary,
5b) by solving the optimal fitting coefficient that above formula least square fitting 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
The dimensionality reduction echo data z of the wave beam Doppler domain obtaining, k represents the sequence number of the data of range cell to be detected, is less than institute
Need the natural number of detecting distance unit sum;φkRepresent the data basic matrix of range cell to be detected;Represent to be detected away from
Optimal fitting coefficient from the data of unit;||·||2Represent 2 norms solving vector.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410276594.9A CN104459659B (en) | 2014-06-19 | 2014-06-19 | Wave beam Doppler domain space-time two-dimensional self-adaptive processing method based on data fitting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410276594.9A CN104459659B (en) | 2014-06-19 | 2014-06-19 | Wave beam Doppler domain space-time two-dimensional self-adaptive processing method based on data fitting |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104459659A CN104459659A (en) | 2015-03-25 |
CN104459659B true CN104459659B (en) | 2017-01-25 |
Family
ID=52906012
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410276594.9A Expired - Fee Related CN104459659B (en) | 2014-06-19 | 2014-06-19 | Wave beam Doppler domain space-time two-dimensional self-adaptive processing method based on data fitting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104459659B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018223285A1 (en) * | 2017-06-06 | 2018-12-13 | 深圳大学 | Stap method and device with beam-doppler pattern sparsity constraint |
CN109212502B (en) * | 2018-09-10 | 2022-12-27 | 中国航天科工集团八五一一研究所 | Method for realizing sampling matrix inversion non-uniform detector based on reweighting |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975940A (en) * | 2010-09-27 | 2011-02-16 | 北京理工大学 | Segmentation combination-based adaptive constant false alarm rate target detection method for SAR image |
CN102288949A (en) * | 2011-05-13 | 2011-12-21 | 中国民航大学 | Optimal processor based maneuvering target detection method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010054980A (en) * | 2008-08-29 | 2010-03-11 | Hitachi Displays Ltd | Liquid crystal display device |
-
2014
- 2014-06-19 CN CN201410276594.9A patent/CN104459659B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975940A (en) * | 2010-09-27 | 2011-02-16 | 北京理工大学 | Segmentation combination-based adaptive constant false alarm rate target detection method for SAR image |
CN102288949A (en) * | 2011-05-13 | 2011-12-21 | 中国民航大学 | Optimal processor based maneuvering target detection method |
Non-Patent Citations (2)
Title |
---|
《基于STAP的机载预警雷达信号处理技术研究》;杨芸;《中国优秀硕士学位论文全文数据库》;20111215(第12期);正文第12-13、33页 * |
《机载相控阵雷达STAP及目标参数估计方法研究》;吴建新;《中国博士学位论文全文数据库》;20091215(第12期);正文第17-19、40页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104459659A (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 | |
CN103439692B (en) | STAP method based on wide symmetrical characteristic of covariance matrix | |
CN104459657B (en) | Extension factorization space-time two-dimensional self-adaptive processing method based on data fitting | |
CN104459660B (en) | Array element Doppler domain space-time two-dimensional self-adaptive processing method based on data fitting | |
CN103454624B (en) | The direct data domain moving target detection method of spectrum time empty based on dimensionality reduction sparse reconstruct | |
CN103364764B (en) | Airborne radar non-stationary clutter suppression method | |
CN105929371B (en) | A kind of airborne radar clutter suppression method based on covariance matrix | |
CN104569949B (en) | Radar target detection method based on combined adaptive normalized matched filter | |
CN106443633A (en) | Shipborne high frequency ground wave radar sea clutter time domain suppression method | |
CN105403856A (en) | DOA (direction of arrival) estimation method based on nested minimum redundant array | |
CN106529428A (en) | Underwater target recognition method based on deep learning | |
CN106054153A (en) | Sea clutter zone target detection and adaptive clutter inhibition method based on fractional transform | |
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 | |
CN107085205A (en) | Self-adapting detecting method based on clutter covariance matrix structural information | |
CN105403864A (en) | Two-dimension high-frequency ground wave radar sea clutter suppression method | |
CN104535973A (en) | Target detection method by use of airborne early warning radar | |
CN107942324A (en) | Multi-frame joint Small object double check method based on Doppler's guiding | |
CN105204007A (en) | Range-spread target polarization adaptive detection method without auxiliary data | |
CN103176168A (en) | Short-range cluster cancellation method for airborne non-side-looking array radar | |
CN104459659B (en) | Wave beam Doppler domain space-time two-dimensional self-adaptive processing method based on data fitting | |
CN105223554A (en) | Based on the space-time adaptive Monopulse estimation method of Doppler's triple channel Combined Treatment | |
CN104316914A (en) | Radar target self-adaptation detection method depending on shape parameters | |
CN105372635A (en) | Improved dimension-reduction space-time adaptive processing-based ship-borne high-frequency ground wave radar sea clutter suppression method | |
CN105445703A (en) | Two-stage time space adaptive processing method for airborne radar time space echo data |
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 |