CN104459658A - Unite-domain positioning space-time two-dimensional self-adaptive processing method based on data fitting - Google Patents
Unite-domain positioning space-time two-dimensional self-adaptive processing method based on data fitting Download PDFInfo
- Publication number
- CN104459658A CN104459658A CN201410276307.4A CN201410276307A CN104459658A CN 104459658 A CN104459658 A CN 104459658A CN 201410276307 A CN201410276307 A CN 201410276307A CN 104459658 A CN104459658 A CN 104459658A
- Authority
- CN
- China
- Prior art keywords
- data
- matrix
- represent
- detected
- range unit
- 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.)
- Granted
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 unite-domain positioning space-time two-dimensional self-adaptive processing method based on data fitting, and relates to radar technology. The method comprises the steps that 1, JDL 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 unite-domain positioning 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, reduces the false alarm rate, improves movable target detection probability, and has the good application prospect in measured data processing and STAP performance improvement.
Description
Technical field
The invention belongs to communication technical field, relate to Radar Technology, particularly location, a kind of associating territory based on data fitting space-time adaptive processing method, for airborne early warn ing radar signal transacting.
Background technology
The main task of airborne early warn ing radar is the detection of a target in complex clutter background, and positions tracking to it, and carries out effectively suppressing being the core means improving airborne early warn ing radar serviceability to clutter.Space-time adaptive process (space-time adaptive processing, STAP) technology makes full use of spatial domain and time-domain information, while coherent accumulation is carried out to echo signal, by space-time adaptive process filtering ground clutter, realize the effective detection of airborne early warn ing radar to target, the E2-D airborne early warn ing radar as the U.S. just adopts this technology.
In actual applications, mainly there are following two aspect problems in STAP technology: on the one hand, in clutter environment heterogeneous, obtain the abundant independent same distribution for estimate covariance matrix (independent and identically distributed, IID) training sample very difficult; On the other hand, even if the demand of training sample is met, the excessive problem of full space time processing calculated amount can cause real-time to be difficult to ensure.
For solving the problem, promoting STAP technology more practical, there has been proposed many innovative approachs or method.The patent of invention " space-time adaptive processing method under non-homogeneous clutter environment " (number of patent application 201010129723.3, publication No. CN101819269A) of Tsing-Hua University's application discloses a kind of overcomplete sparse representation method of super-resolution estimation clutter space-time two-dimensional spectrum in non-homogeneous clutter environment.The method achieve when independent same distribution sample number deficiency, utilize single frames training sample to estimate clutter covariance matrix, thus avoid strong non-homogeneous clutter environment on the impact of self-adaptive processing effect.But, still there is following main deficiency in the method: one is the problem that operand is large, super complete radix order clutter spectrum being carried out to rarefaction representation is uncertain, but much larger than degree of freedom in system, and degree of freedom in system is usually thousands of in reality, operand required in the covariance matrix restructuring procedure of each range unit sample is like this very large, is unfavorable for real-time process, thus has influence on the effect in practical engineering application.Two is base mismatch problems, and the base used in the method is actually one group of interpolation discrete fourier DFT vector, and the echo data of reality is because exist the non-ideal factors such as error, can not use given base vector rarefaction representation; In addition, the method needs to carry out discretize expression to plane during sky, the signal not being positioned at discrete grid block dot center can be made like this to reveal to all net points, thus destroy the openness of echo data.
Summary of the invention
For the deficiency of the overcomplete sparse representation method of above-mentioned super-resolution estimation clutter space-time two-dimensional spectrum in non-homogeneous clutter environment, the present invention proposes location, a kind of associating territory based on data fitting space-time adaptive processing (Joint-Domain Localized STAP, JDL-STAP) method, the detection to echo signal can be realized, and solve the problem of the huge and base mismatch of operand in above-mentioned sparse recovery STAP method, simultaneously relative to location, classical joint territory JDL-STAP method, improve the detection probability of moving-target, reduce the false alarm rate of detection.
For achieving the above object, the present invention is pre-to realize by the following technical solutions.
Based on location, an associating territory space-time adaptive processing method for data fitting, it is characterized in that, comprise the following steps:
Step 1, utilizes airborne early warn ing radar aerial receiver to carry the space-time two-dimensional echo data x of early warning radar, carries out location, associating territory JDL dimensionality reduction to space-time two-dimensional echo data x, obtains dimensionality reduction echo data z and steering vector s time empty
z-JDL;
Step 2, constructs the Data Base matrix Φ of range unit to be detected
k;
Step 3, steering vector s when utilizing empty
z-JDLderive blocking matrix B; Utilize the data z of blocking matrix B and range unit to be detected
k, obtain target and block later auxiliary echo data
the data z of this range unit to be detected
ktake from the dimensionality reduction echo data z in step 1 successively; Utilize the Data Base matrix Φ of blocking matrix B and range unit to be detected
kobtain target and block later Data Base matrix
Step 4, utilizes target to block later Data Base matrix
to auxiliary echo data
carry out data fitting, obtain auxiliary echo data
fitting coefficient
Step 5, utilizes the Data Base matrix Φ of range unit to be detected
kwith auxiliary echo data
fitting coefficient
treat detecting distance cell data z
kcarry out data fitting, obtain the data z of range unit to be detected
kminimum error of fitting y
k; The data z of this range unit to be detected
ktake from the dimensionality reduction echo data z in step 1 successively;
Step 6, treats the minimum error of fitting y of the data of detecting distance unit
kcarry out CA-CFAR detection, and there is target or there is not target in output.
Feature and further improvement of technique scheme are:
(1) step 1 specifically comprises following sub-step:
1a) utilize airborne early warn ing 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 is MN dimension, wherein M indicating impulse number, and N represents the array number of airborne early warn ing radar antenna;
The dimensionality reduction transition matrix P of location, associating territory JDL dimensionality reduction 1b) is obtained according to following formula
jDL:
Wherein, P
t(f
i-1, f
i, f
i+1) represent the time domain transition matrix become by three Doppler filter group closed on, f
irepresent the normalization Doppler frequency of target, f
i-1, f
i+1represent two the normalization Doppler frequencies adjacent with target; I be less than or equal to M and be more than or equal to 1 natural number, M indicating impulse number, P
s(g
j-1, g
j, g
j+1) representing the spatial domain transition matrix be made up of three spatial filters, the wave filter in the middle of these three spatial filters is the spatial domain steering vector of target, g
jfor the normalization spatial frequency of target, g
j-1, g
j+1for two contiguous normalization spatial frequencys, j be less than or equal to M and be more than or equal to 1 natural number;
represent Kronecker product; Dimensionality reduction transition matrix P
jDLdimension be MN × 9 tie up;
1c) according to dimensionality reduction transition matrix P
jDLobtain the dimensionality reduction echo data z after the JDL dimensionality reduction of location, associating territory, be expressed as following formula:
Wherein, ()
hrepresent conjugate transposition operation, according to dimensionality reduction transition matrix P
jDLdimension be that MN × 9 are tieed up and space-time two-dimensional echo data x be the dimension that MN dimension can obtain dimensionality reduction echo data z is 9 dimensions;
Steering vector s during empty after the JDL dimensionality reduction of location, associating territory
z-JDL, be expressed as following formula:
Wherein, ()
trepresent matrix transpose operation, s
z-JDLfor steering vector during empty after dimensionality reduction.
(2) step 2 specifically comprises following sub-step:
Data Base matrix Φ
kbe expressed as following formula:
Φ
k=[z
1z
2…z
L]∈C
D×L
Wherein, Φ
krepresent the Data Base matrix of a kth range unit to be detected, k represents the sequence number of the data of range unit to be detected, for being no more than the natural number needing all range unit sums to be detected to be processed; D represents the dimension of dimensionality reduction echo data z, and L represents the number of the data of closing on range unit, C
d × Lrepresent the complex matrix space of D × L dimension, z
l(l=1,2 ..., L) represent and close on the data of range unit, close on range unit and represent and range unit near range unit to be detected do not comprise two range units for the protection of target.
(3) step 3 specifically comprises following sub-step:
3a) to steering vector s during sky
z-JDLtransposed form carry out svd, obtain left singular matrix U, right singular matrix V and singular value matrix Λ; Decomposition formula as shown in the formula:
Wherein, U represents left singular matrix, V represents right singular matrix, Λ represents singular value matrix; ()
trepresent matrix transpose operation, ()
hrepresent conjugate transposition operation;
3b) blocking matrix B is made up of to the arrangement of N row order the 2nd of right singular matrix, derives blocking matrix B by right singular matrix V and following formula;
B=[V(:,2:N)]
T
Wherein, V (:, 2:N) represent that all row, the 2nd of right singular matrix V are to N column matrix element; ()
trepresent matrix transpose operation;
3c) obtain auxiliary echo data according to blocking matrix B
later Data Base matrix is blocked with target
be shown below:
Wherein, z
krepresent the data of a kth range unit to be detected, the data of this range unit to be detected take from the dimensionality reduction echo data z in step 1 successively, k represents the sequence number of the data of range unit to be detected, for being no more than the natural number needing all range unit sums to be detected to be processed;
represent that target blocks later auxiliary echo data,
represent that target blocks later Data Base matrix, B is blocking matrix.
(4) step 4 specifically comprises following sub-step:
4a) block later Data Base matrix by target
later auxiliary echo data is blocked to target
carry out least square fitting expression, be expressed as following formula:
Wherein, || ||
2represent 2 norms solving vector,
represent auxiliary echo data
fitting coefficient;
4b) according to step 4a) in least square fitting represent, solve auxiliary echo data
fitting coefficient
least square solution
be expressed as following formula:
Wherein, ()
-1representing matrix inversion operation, ()
hrepresent conjugate transposition operation,
represent that target blocks later Data Base matrix,
represent that target blocks later auxiliary echo data.
(5) step 5 specifically comprises following sub-step:
5a) with the Data Base matrix Φ of range unit to be detected
k, treat the data z of detecting distance unit
kcarry out least square fitting; And with auxiliary echo data
fitting coefficient
least square solution
calculate and optimize restrained boundary η
k, realize especially by following formula:
Wherein, α
krepresent the fitting coefficient of the data of range unit to be detected, y
krepresent 2 norms of minimum error of fitting vector, η
krepresent and optimize restrained boundary,
The optimal fitting coefficient of the data of range unit to be detected 5b) is obtained by solving above formula
According to the optimal fitting coefficient of the data of range unit to be detected
obtain the minimum error of fitting y of the data of range unit to be detected
k, obtained by following formula:
Wherein, z
krepresent the data of a kth range unit to be detected, the data of this range unit to be detected take from the dimensionality reduction echo data z in step 1 successively, k represents the sequence number of the data of range unit to be detected, for being no more than the natural number needing all range unit sums to be detected to be processed; Φ
krepresent the Data Base matrix of range unit to be detected;
represent the optimal fitting coefficient of the data of range unit to be detected; || ||
2represent 2 norms solving vector.
Compared with prior art, the present invention has outstanding substantive distinguishing features and significant progressive, and the present invention sparsely recovers compared with STAP method with existing, and the present invention has the following advantages:
(1) for the problem of operand, the super complete radix order that sparse restoration methods carries out rarefaction representation to clutter spectrum is uncertain, but much larger than degree of freedom in system, and degree of freedom in system is usually thousands of in reality, operand required in the covariance matrix restructuring procedure of each range unit sample is like this very large, be unfavorable for real-time process, thus have influence on its practical engineering application effect.The present invention is by combining with dimensionality reduction STAP method JDL, and the little optimization problem optimization problem that original MN ties up being decomposed into M 9 dimensions solves, and greatly reduces the computation complexity of algorithm, solves the problem that in sparse recovery STAP method, operand is large.
(2) for base mismatch problems, the expression base used in sparse restoration methods is actually one group of interpolated DFT vector, and because there is the non-ideal factors such as error in real data, can not with give base vector rarefaction representation, in addition, the method needs to carry out discretize expression to plane during sky, the signal not being positioned at discrete grid block dot center can be made like this to reveal to all net points, thus destroy the openness of data.The present invention directly uses real data sample as basis matrix, instead of the ideal flowing pattern adopting certain to determine, thus has evaded the problem of base mismatch in above-mentioned sparse recovery STAP method.
Visible, in prior art, STAP carries out airborne radar signal transacting usually in high-dimensional data space, greatly will increase computational complexity and required training sample number like this, in order to obtain better performance, the present invention is directed to the problems referred to above, the data sample utilizing range unit to be detected to close on range unit carrys out Optimal Fitting range unit clutter data to be detected, thus improves the clutter recognition performance of Adaptive Signal Processing, improves the detection probability of target.
Accompanying drawing explanation
Below in conjunction with the drawings and specific embodiments, the present invention will be further described.
Fig. 1 is process flow diagram of the present invention;
Fig. 2 is general pulse doppler processing result figure;
Fig. 3 is the change curve of detection probability with border compensation rate;
Fig. 4 is the change curve of detection probability with sample number;
Fig. 5 is the Output rusults figure of the minimum error of fitting of the data of target proximity range unit; Fig. 5 (a) is the Output rusults of the minimum error of fitting of a segment distance unit scope after conventional JDL method; Fig. 5 (c) is the Output rusults through the minimum error of fitting of a segment distance unit scope after the inventive method process; Fig. 5 (b) is the Output rusults of the minimum error of fitting of another segment distance unit scope after conventional JDL method; Fig. 5 (d) is the Output rusults of the minimum error of fitting of another segment distance unit scope after the inventive method;
Fig. 6 is the change curve of detection probability with target power;
Fig. 7 is the change curve graph of a relation of detection probability with false alarm rate.
Embodiment
Below in conjunction with accompanying drawing, the invention will be further described.
With reference to Fig. 1, location, a kind of associating territory based on data fitting of the present invention space-time adaptive processing method is described, the present invention is used for airborne early warn ing radar signal transacting, and its concrete implementation step is as follows:
Step 1, utilizes airborne early warn ing radar aerial receiver to carry the space-time two-dimensional echo data x of early warning radar, carries out location, associating territory JDL dimensionality reduction to space-time two-dimensional echo data x, obtains dimensionality reduction echo data z and steering vector s time empty
z-JDL;
Step 1 specifically comprises following sub-step:
1a) utilize airborne early warn ing 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 is MN dimension, wherein M indicating impulse number, and N represents the array number of airborne early warn ing radar antenna;
The dimensionality reduction transition matrix P of location, associating territory JDL dimensionality reduction 1b) is obtained according to following formula
jDL:
Wherein, P
t(f
i-1, f
i, f
i+1) represent the time domain transition matrix become by three Doppler filter group closed on, f
irepresent the normalization Doppler frequency of target, f
i-1, f
i+1represent two the normalization Doppler frequencies adjacent with target; I be less than or equal to M and be more than or equal to 1 natural number, M indicating impulse number, P
s(g
j-1, g
j, g
j+1) representing the spatial domain transition matrix be made up of three spatial filters, the wave filter in the middle of these three spatial filters is the spatial domain steering vector of target, g
jfor the normalization spatial frequency of target, g
j-1, g
j+1for two contiguous normalization spatial frequencys, j be less than or equal to M and be more than or equal to 1 natural number;
represent Kronecker product; Dimensionality reduction transition matrix P
jDLdimension be MN × 9 tie up;
1c) according to dimensionality reduction transition matrix P
jDLobtain the dimensionality reduction echo data z after the JDL dimensionality reduction of location, associating territory, be expressed as following formula:
Wherein, ()
hrepresent conjugate transposition operation, according to dimensionality reduction transition matrix P
jDLdimension be that MN × 9 are tieed up and space-time two-dimensional echo data x be the dimension that MN dimension can obtain dimensionality reduction echo data z is 9 dimensions;
Steering vector s during empty after the JDL dimensionality reduction of location, associating territory
z-JDL, be expressed as following formula:
Wherein, ()
trepresent matrix transpose operation, s
z-JDLfor steering vector during empty after dimensionality reduction.
During general sky, steering vector s can be expressed as the Kronecker product form of time domain steering vector and spatial domain steering vector, namely
f
dand f
srepresent normalized Doppler frequency and spatial frequency respectively; Time domain and spatial domain steering vector can be expressed as respectively
Step 2, constructs the Data Base matrix Φ of range unit to be detected
k;
Data Base matrix Φ
kbe expressed as following formula:
Φ
k=[z
1z
2… z
L]∈C
D×L
Wherein, Φ
krepresent the Data Base matrix of a kth range unit to be detected, k represents the sequence number of the data of range unit to be detected, for being no more than the natural number needing all range unit sums to be detected to be processed; D represents the dimension of dimensionality reduction echo data z, and L represents the number of the data of closing on range unit, C
d × Lrepresent the complex matrix space of D × L dimension, z
l(l=1,2 ..., L) represent and close on the data of range unit, close on range unit and represent and range unit near range unit to be detected do not comprise two range units for the protection of target.
Step 3, steering vector s when utilizing empty
z-JDLderive blocking matrix B; Utilize the data z of blocking matrix B and range unit to be detected
k, obtain target and block later auxiliary echo data
the data z of this range unit to be detected
ktake from the dimensionality reduction echo data z in step 1 successively; Utilize the Data Base matrix Φ of blocking matrix B and range unit to be detected
kobtain target and block later Data Base matrix
Step 3 specifically comprises following sub-step:
3a) to steering vector s during sky
z-JDLtransposed form carry out svd, obtain left singular matrix U, right singular matrix V and singular value matrix Λ; Decomposition formula as shown in the formula:
Wherein, U represents left singular matrix, V represents right singular matrix, Λ represents singular value matrix; ()
trepresent matrix transpose operation, ()
hrepresent conjugate transposition operation;
3b) blocking matrix B is made up of to the arrangement of N row order the 2nd of right singular matrix, derives blocking matrix B by right singular matrix V and following formula;
B=[V(:,2:N)]
T
Wherein, V (:, 2:N) represent that all row, the 2nd of right singular matrix V are to N column matrix element; ()
trepresent matrix transpose operation;
3c) obtain auxiliary echo data according to blocking matrix B
later Data Base matrix is blocked with target
be shown below:
Wherein, z
krepresent the data of a kth range unit to be detected, the data of this range unit to be detected take from the dimensionality reduction echo data z in step 1 successively, k represents the sequence number of the data of range unit to be detected, for being no more than the natural number needing all range unit sums to be detected to be processed;
represent that target blocks later auxiliary echo data,
represent that target blocks later Data Base matrix, B is blocking matrix.
Step 4, utilizes target to block later Data Base matrix
to auxiliary echo data
carry out data fitting, obtain auxiliary echo data
fitting coefficient
Step 4 specifically comprises following sub-step:
4a) block later Data Base matrix by target
later auxiliary echo data is blocked to target
carry out least square fitting expression, be expressed as following formula:
Wherein, || ||
2represent 2 norms solving vector,
represent auxiliary echo data
fitting coefficient;
4b) according to step 4a) in least square fitting represent, solve auxiliary echo data
fitting coefficient
least square solution
be expressed as following formula:
Wherein, ()
-1representing matrix inversion operation, ()
hrepresent conjugate transposition operation,
represent that target blocks later Data Base matrix,
represent that target blocks later auxiliary echo data.
Step 5, utilizes the Data Base matrix Φ of range unit to be detected
kwith auxiliary echo data
fitting coefficient
treat detecting distance cell data z
kcarry out data fitting, obtain the data z of range unit to be detected
kminimum error of fitting y
k; The data z of this range unit to be detected
ktake from the dimensionality reduction echo data z in step 1 successively;
Step 5 specifically comprises following sub-step:
5a) with the Data Base matrix Φ of range unit to be detected
k, treat the data z of detecting distance unit
kcarry out least square fitting; And with auxiliary echo data
fitting coefficient
least square solution
calculate and optimize restrained boundary η
k, realize especially by following formula:
Wherein, α
krepresent the fitting coefficient of the data of range unit to be detected, y
krepresent 2 norms of minimum error of fitting vector, η
krepresent and optimize restrained boundary,
The optimal fitting coefficient of the data of range unit to be detected 5b) is obtained by solving above formula
According to the optimal fitting coefficient of the data of range unit to be detected
obtain the minimum error of fitting y of the data of range unit to be detected
k, obtained by following formula:
Wherein, z
krepresent the data of a kth range unit to be detected, the data of this range unit to be detected take from the dimensionality reduction echo data z in step 1 successively, k represents the sequence number of the data of range unit to be detected, for being no more than the natural number needing all range unit sums to be detected to be processed; Φ
krepresent the Data Base matrix of range unit to be detected;
represent the optimal fitting coefficient of the data of range unit to be detected; || ||
2represent 2 norms solving vector.
Step 6, treats the minimum error of fitting y of the data of detecting distance unit
kcarry out CA-CFAR detection, and there is target or there is not target in output.
Treat the minimum error of fitting y of the data of detecting distance unit
kcarry out CA-CFAR detection (Cell-averaging constant false alarm rate, CA-CFAR), by minimum error of fitting y
kcompare with the mean value of the minimum error of fitting of the data of the range unit of the surrounding of range unit to be detected, according to minimum error of fitting y
kjudge whether to there is target with toaverage ratio size, finally will there is target or there is not target output.
Below in conjunction with emulation experiment, effect of the present invention is described further.
(1) experiment condition:
Experiment of the present invention is carried out under MATLAB7.11 software.In experimental design of the present invention, in order to check the inventive method validity in actual applications, process data are from the measured data MCARM data set of external admission, under its radar is operated in positive side battle array pattern, the spatial domain array element that can be used for Adaptive Signal Processing that front is arranged by 2 row 11 forms, and its part system parameter is with reference to table 1.In order to test conveniently, the present invention extracts the first row 11 array elements in MCARM data, the echo data of front 32 pulses 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 range unit sequence number.Range unit also claims range gate in the present invention.Under as can be seen from Figure 2 although this radar is operated in positive side battle array situation, the strong clutter component in space-time two-dimensional echo data, is subject to the existence of various false target in the different modulation of landform reflectivity and actual environment, has certain heterogeneity.
Accompanying drawing 3 is the change curve of detection probability with border compensation rate.Transverse axis represents boundary compensation, and the longitudinal axis represents detection probability.The additional compensation rate in the optimization upper bound of gained when boundary compensation refers in solution procedure 5 and optimizes restrained boundary.As can be seen from Figure 3, substantially near 0 compensation rate, detection probability can obtain maximal value, and this also just describes step 5 calculation optimization restrained boundary in the present invention is effective and optimum.
Accompanying drawing 4 is the change curve of detection probability with sample number.Sample number refers to the number of the data of closing on range unit in the present invention.Transverse axis represents sample number, and the longitudinal axis represents detection probability.In the present invention after data JDL dimension-reduction treatment, data dimension (i.e. degree of freedom) reduces to 9.As can be seen from Figure 4, detection probability reaches maximal value substantially when sample number is more than 36, this just illustrates, the data that in the present invention, Data Base matrix can close on range unit by 4 times of number of degrees of freedom, objects near range unit to be detected are formed, and obtain optimum detection perform.
Fig. 5 is the Output rusults figure of the minimum error of fitting of the data of target proximity range unit.In the present invention, the data of target proximity range unit comprise the data of many groups range unit to be detected.Transverse axis represents the sequence number of range unit, and the longitudinal axis represents power, represents with dB.In Fig. 5, curve represents the Output rusults of the minimum error of fitting of a segment distance unit scope in target Doppler passage, horizontal linear represents the average output value of the minimum error of fitting of this segment distance unit scope, circle represents the minimum error of fitting output valve of target range unit, and in Fig. 5, " difference " represents the difference between the average output value of the minimum error of fitting output valve of target range unit and the minimum error of fitting of this segment distance unit scope.As can be seen from Figure 5, after process of the present invention, target exports more obvious, Fig. 5 (a) compared with Fig. 5 (c) when range unit is 150 power improve 11.9dB, Fig. 5 (b) compared with Fig. 5 (d) when range unit is 200 power improve 11.1dB, thus more can be conducive to the detection of echo signal.
Fig. 6 is the change curve of detection probability with target power.Transverse axis represents the artificial target power arranged, and the longitudinal axis represents the detection probability of this target.Associating territory localization method (DF-JDL) based on data fitting that the sparse restoration methods of direct data domain (D3SR), classical joint territory localization method (JDL) and the present invention propose is compared in Fig. 6.As can be seen from Figure 6, the detection perform of the inventive method is much better than other two kinds of methods.
Fig. 7 is the change curve graph of a relation of detection probability with false alarm rate.Transverse axis represents false alarm rate, and the longitudinal axis represents detection probability.Associating territory localization method (DF-JDL) based on data fitting of classical joint territory localization method (JDL) and the present invention's proposition is compared in Fig. 7.As can be seen from Figure 7, the inventive method, relative to traditional JDL method, has higher detection probability, and have lower false alarm rate under identical detection probability under certain false alarm rate.
Claims (6)
1., based on location, an associating territory space-time adaptive processing method for data fitting, it is characterized in that, comprise the following steps:
Step 1, utilizes airborne early warn ing radar aerial receiver to carry the space-time two-dimensional echo data x of early warning radar, carries out location, associating territory JDL dimensionality reduction to space-time two-dimensional echo data x, obtains dimensionality reduction echo data z and steering vector s time empty
z-JDL;
Step 2, constructs the Data Base matrix Φ of range unit to be detected
k;
Step 3, steering vector s when utilizing empty
z-JDLderive blocking matrix B; Utilize the data z of blocking matrix B and range unit to be detected
k, obtain target and block later auxiliary echo data
the data z of this range unit to be detected
ktake from the dimensionality reduction echo data z in step 1 successively; Utilize the Data Base matrix Φ of blocking matrix B and range unit to be detected
kobtain target and block later Data Base matrix
Step 4, utilizes target to block later Data Base matrix
to auxiliary echo data
carry out data fitting, obtain auxiliary echo data
fitting coefficient
Step 5, utilizes the Data Base matrix Φ of range unit to be detected
kwith auxiliary echo data
fitting coefficient
treat detecting distance cell data z
kcarry out data fitting, obtain the data z of range unit to be detected
kminimum error of fitting y
k; The data z of this range unit to be detected
ktake from the dimensionality reduction echo data z in step 1 successively;
Step 6, treats the minimum error of fitting y of the data of detecting distance unit
kcarry out CA-CFAR detection, and there is target or there is not target in output.
2. location, the associating territory based on data fitting according to claim 1 space-time adaptive processing method, it is characterized in that, step 1 specifically comprises following sub-step:
1a) utilize airborne early warn ing 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 is MN dimension, wherein M indicating impulse number, and N represents the array number of airborne early warn ing radar antenna;
The dimensionality reduction transition matrix P of location, associating territory JDL dimensionality reduction 1b) is obtained according to following formula
jDL:
Wherein, P
t(f
i-1, f
i, f
i+1) represent the time domain transition matrix become by three Doppler filter group closed on, f
irepresent the normalization Doppler frequency of target, f
i-1, f
i+1represent two the normalization Doppler frequencies adjacent with target; I be less than or equal to M and be more than or equal to 1 natural number, M indicating impulse number, P
s(g
j-1, g
j, g
j+1) representing the spatial domain transition matrix be made up of three spatial filters, the wave filter in the middle of these three spatial filters is the spatial domain steering vector of target, g
jfor the normalization spatial frequency of target, g
j-1, g
j+1for two contiguous normalization spatial frequencys, j be less than or equal to M and be more than or equal to 1 natural number;
represent Kronecker product; Dimensionality reduction transition matrix P
jDLdimension be MN × 9 tie up;
1c) according to dimensionality reduction transition matrix P
jDLobtain the dimensionality reduction echo data z after the JDL dimensionality reduction of location, associating territory, be expressed as following formula:
Wherein, ()
hrepresent conjugate transposition operation, according to dimensionality reduction transition matrix P
jDLdimension be that MN × 9 are tieed up and space-time two-dimensional echo data x be the dimension that MN dimension can obtain dimensionality reduction echo data z is 9 dimensions;
Steering vector s during empty after the JDL dimensionality reduction of location, associating territory
z-JDL, be expressed as following formula:
Wherein, ()
trepresent matrix transpose operation, s
z-JDLfor steering vector during empty after dimensionality reduction.
3. location, the associating territory based on data fitting according to claim 1 space-time adaptive processing method, it is characterized in that, step 2 specifically comprises:
Data Base matrix Φ
kbe expressed as following formula:
Φ
k=[z
1z
2…z
L]∈C
D×L
Wherein, Φ
krepresent the Data Base matrix of a kth range unit to be detected, k represents the sequence number of the data of range unit to be detected, for being no more than the natural number needing all range unit sums to be detected to be processed; D represents the dimension of dimensionality reduction echo data z, and L represents the number of the data of closing on range unit, C
d × Lrepresent the complex matrix space of D × L dimension, z
l(l=1,2 ..., L) represent and close on the data of range unit, close on range unit and represent and range unit near range unit to be detected do not comprise two range units for the protection of target.
4. location, the associating territory based on data fitting according to claim 1 space-time adaptive processing method, it is characterized in that, step 3 specifically comprises following sub-step:
3a) to steering vector s during sky
z-JDLtransposed form carry out svd, obtain left singular matrix U, right singular matrix V and singular value matrix Λ; Decomposition formula as shown in the formula:
Wherein, U represents left singular matrix, V represents right singular matrix, Λ represents singular value matrix; ()
trepresent matrix transpose operation, ()
hrepresent conjugate transposition operation;
3b) blocking matrix B is made up of to the arrangement of N row order the 2nd of right singular matrix, derives blocking matrix B by right singular matrix V and following formula;
B=[V(:,2:N)]
T
Wherein, V (:, 2:N) represent that all row, the 2nd of right singular matrix V are to N column matrix element; ()
trepresent matrix transpose operation;
3c) obtain auxiliary echo data according to blocking matrix B
later Data Base matrix is blocked with target
be shown below:
Wherein, z
krepresent the data of a kth range unit to be detected, the data of this range unit to be detected take from the dimensionality reduction echo data z in step 1 successively, k represents the sequence number of the data of range unit to be detected, for being no more than the natural number needing all range unit sums to be detected to be processed;
represent that target blocks later auxiliary echo data,
represent that target blocks later Data Base matrix, B is blocking matrix.
5. location, the associating territory based on data fitting according to claim 1 space-time adaptive processing method, it is characterized in that, step 4 specifically comprises following sub-step:
4a) block later Data Base matrix by target
later auxiliary echo data is blocked to target
carry out least square fitting expression, be expressed as following formula:
Wherein, || ||
2represent 2 norms solving vector,
represent auxiliary echo data
fitting coefficient;
4b) according to step 4a) in least square fitting represent, solve auxiliary echo data
fitting coefficient
least square solution
be expressed as following formula:
Wherein, ()
-1representing matrix inversion operation, ()
hrepresent conjugate transposition operation,
represent that target blocks later Data Base matrix,
represent that target blocks later auxiliary echo data.
6. location, the associating territory based on data fitting according to claim 1 space-time adaptive processing method, it is characterized in that, step 5 specifically comprises following sub-step:
5a) with the Data Base matrix Φ of range unit to be detected
k, treat the data z of detecting distance unit
kcarry out least square fitting; And with auxiliary echo data
fitting coefficient
least square solution
calculate and optimize restrained boundary η
k, realize especially by following formula:
Wherein, α
krepresent the fitting coefficient of the data of range unit to be detected, y
krepresent 2 norms of minimum error of fitting vector, η
krepresent and optimize restrained boundary,
The optimal fitting coefficient of the data of range unit to be detected 5b) is obtained by solving above formula
According to the optimal fitting coefficient of the data of range unit to be detected
obtain the minimum error of fitting y of the data of range unit to be detected
k, obtained by following formula:
Wherein, z
krepresent the data of a kth range unit to be detected, the data of this range unit to be detected take from the dimensionality reduction echo data z in step 1 successively, k represents the sequence number of the data of range unit to be detected, for being no more than the natural number needing all range unit sums to be detected to be processed; Φ
krepresent the Data Base matrix of range unit to be detected;
represent the optimal fitting coefficient of the data of range unit to be detected; || ||
2represent 2 norms solving vector.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410276307.4A CN104459658B (en) | 2014-06-19 | 2014-06-19 | Unite-domain positioning space-time two-dimensional self-adaptive processing method based on data fitting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410276307.4A CN104459658B (en) | 2014-06-19 | 2014-06-19 | Unite-domain positioning space-time two-dimensional self-adaptive processing method based on data fitting |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104459658A true CN104459658A (en) | 2015-03-25 |
CN104459658B CN104459658B (en) | 2017-01-25 |
Family
ID=52906011
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410276307.4A Expired - Fee Related CN104459658B (en) | 2014-06-19 | 2014-06-19 | Unite-domain positioning space-time two-dimensional self-adaptive processing method based on data fitting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104459658B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105022040A (en) * | 2015-07-08 | 2015-11-04 | 西安电子科技大学 | Array element error estimation method based on clutter data combined fitting |
CN105044688A (en) * | 2015-08-24 | 2015-11-11 | 西安电子科技大学 | Radar robust space-time adaption processing method based on iterative subspace tracking algorithm |
CN105093189A (en) * | 2015-09-08 | 2015-11-25 | 西安电子科技大学 | Aircraft radar target detection method based on GCV (generalized cross validation) |
CN105223557A (en) * | 2015-10-29 | 2016-01-06 | 西安电子科技大学 | Based on the airborne early warn ing radar clutter suppression method of accessory channel |
CN105319538A (en) * | 2015-11-24 | 2016-02-10 | 西安电子科技大学 | Accessory channel-based space-time self-adaptive clutter suppression method |
CN107121665A (en) * | 2017-05-03 | 2017-09-01 | 电子科技大学 | A kind of passive location method of the near field coherent source based on Sparse Array |
CN109238547A (en) * | 2018-08-28 | 2019-01-18 | 西北工业大学 | A kind of edge filling method of aircraft airborne equipment monitoring data |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6292592B1 (en) * | 1998-10-19 | 2001-09-18 | Raytheon Company | Efficient multi-resolution space-time adaptive processor |
CN1088199C (en) * | 1998-12-14 | 2002-07-24 | 中国人民解放军空军雷达学院 | Method for processing space-time multibeam adaptive signals |
CN101561497B (en) * | 2009-05-22 | 2011-09-21 | 西安电子科技大学 | Airborne radar clutter suppression method |
CN101819269A (en) * | 2010-03-19 | 2010-09-01 | 清华大学 | Space-time adaptive processing method under non-homogeneous clutter environment |
-
2014
- 2014-06-19 CN CN201410276307.4A patent/CN104459658B/en not_active Expired - Fee Related
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105022040A (en) * | 2015-07-08 | 2015-11-04 | 西安电子科技大学 | Array element error estimation method based on clutter data combined fitting |
CN105044688A (en) * | 2015-08-24 | 2015-11-11 | 西安电子科技大学 | Radar robust space-time adaption processing method based on iterative subspace tracking algorithm |
CN105044688B (en) * | 2015-08-24 | 2017-07-18 | 西安电子科技大学 | The sane space-time adaptive processing method of radar based on iteration subspace tracking algorithm |
CN105093189A (en) * | 2015-09-08 | 2015-11-25 | 西安电子科技大学 | Aircraft radar target detection method based on GCV (generalized cross validation) |
CN105093189B (en) * | 2015-09-08 | 2017-06-20 | 西安电子科技大学 | Airborne radar object detection method based on GCV |
CN105223557A (en) * | 2015-10-29 | 2016-01-06 | 西安电子科技大学 | Based on the airborne early warn ing radar clutter suppression method of accessory channel |
CN105223557B (en) * | 2015-10-29 | 2017-11-21 | 西安电子科技大学 | Airborne early warning radar clutter suppression method based on accessory channel |
CN105319538A (en) * | 2015-11-24 | 2016-02-10 | 西安电子科技大学 | Accessory channel-based space-time self-adaptive clutter suppression method |
CN107121665A (en) * | 2017-05-03 | 2017-09-01 | 电子科技大学 | A kind of passive location method of the near field coherent source based on Sparse Array |
CN107121665B (en) * | 2017-05-03 | 2019-11-05 | 电子科技大学 | A kind of passive location method of the near field coherent source based on Sparse Array |
CN109238547A (en) * | 2018-08-28 | 2019-01-18 | 西北工业大学 | A kind of edge filling method of aircraft airborne equipment monitoring data |
Also Published As
Publication number | Publication date |
---|---|
CN104459658B (en) | 2017-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104459658A (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 | |
CN103364764B (en) | Airborne radar non-stationary clutter suppression method | |
CN103983944B (en) | Based on arrowband, the far field DOA estimation method of covariance matrix rarefaction representation | |
CN103353596B (en) | Wave beam space domain meter wave radar height measurement method based on compressed sensing | |
CN107589399B (en) | Estimation method of direction of arrival of co-prime array based on singular value decomposition of multi-sampling virtual signal | |
CN103454624B (en) | The direct data domain moving target detection method of spectrum time empty based on dimensionality reduction sparse reconstruct | |
CN105403856A (en) | DOA (direction of arrival) estimation method based on nested minimum redundant array | |
CN104535973A (en) | Target detection method by use of airborne early warning radar | |
CN104459660B (en) | Array element Doppler domain space-time two-dimensional self-adaptive processing method based on data fitting | |
CN102866388B (en) | Iterative computation method for self-adaptive weight number in space time adaptive processing (STAP) | |
CN106872935A (en) | A kind of Electromagnetic Vector Sensor Array Wave arrival direction estimating method based on quaternary number | |
CN103760527B (en) | Single base MIMO radar coherent source Wave arrival direction estimating method | |
CN106443633A (en) | Shipborne high frequency ground wave radar sea clutter time domain suppression method | |
CN104991236A (en) | Monostatic MIMO radar non-circular signal coherent source DOA (Direction Of Arrival) estimation method | |
CN103116162B (en) | High-resolution sonar location method based on sparsity of objective space | |
CN104392034A (en) | L1/2-norm-based sparse linear array optimization method | |
CN103605121B (en) | Wideband radar data fusion method based on rapid sparse Bayesian learning algorithm | |
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 | |
CN105372635A (en) | Improved dimension-reduction space-time adaptive processing-based ship-borne high-frequency ground wave radar sea clutter suppression method | |
CN103792523B (en) | UHF wave band Multichannel radar radial velocity detection method based on tensor product | |
CN105445703A (en) | Two-stage time space adaptive processing method for airborne radar time space echo data | |
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 |