CN105785326B - A kind of optimization method of non-working side battle array radar clutter spectrum registration - Google Patents
A kind of optimization method of non-working side battle array radar clutter spectrum registration Download PDFInfo
- Publication number
- CN105785326B CN105785326B CN201610156725.9A CN201610156725A CN105785326B CN 105785326 B CN105785326 B CN 105785326B CN 201610156725 A CN201610156725 A CN 201610156725A CN 105785326 B CN105785326 B CN 105785326B
- Authority
- CN
- China
- Prior art keywords
- battle array
- working side
- array radar
- clutter
- range gate
- 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.)
- Active
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/023—Interference mitigation, e.g. reducing or avoiding non-intentional interference with other HF-transmitters, base station transmitters for mobile communication or other radar systems, e.g. using electro-magnetic interference [EMI] reduction techniques
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a kind of optimization method of non-working side battle array radar clutter spectrum registration, its main thought is:The L range gate clutter data that non-working side battle array radar receives is obtained, calculates each self-corresponding covariance matrix of the L range gate clutter data respectively, chooses covariance matrix corresponding to k-th of range gateCovariance matrix is corresponded to as range gate to be detected, and covariance matrix is corresponded to as reference using the range gate to be detected, calculates the l' range gate clutter data Xl'Transformation matrix Tl', utilize the Xl'Transformation matrix Tl', to the l' range gate clutter data Xl'Processing is compensated, the l' range gate clutter data Y after being compensatedl';Choose k-th of range gate clutter data XkAs reference spurs, and according to the l' range gate clutter data Y after compensationl', covariance matrix R corresponding to k-th of range gate after compensation is calculatedk;According to covariance matrix R corresponding to k-th of range gate after compensationk, calculate the non-working side battle array radar clutter after space-time adaptive filtering process and suppress improvement factor.
Description
Technical field
The invention belongs to Radar Technology field, more particularly to a kind of optimization side of non-working side battle array radar clutter spectrum registration
Method, precision limitation problem is composed suitable for the radar clutter for solving to obtain using clutter registration method.
Background technology
Space-time adaptive handles (space-time adaptive processing, STAP) technology for strong clutter background
Under target detection at a slow speed it is significant, and carry out corresponding clutter recognition performance and its clutter during target detection at a slow speed
Covariance matrix precision is closely related.Traditional space-time adaptive processing (STAP) technology can obtain for positive side view battle array radar
To good clutter recognition purpose;And for non-working side battle array radar, because non-working side battle array radar clutter frequency spectrum has
Distance dependencies so that directly using the covariance matrix for closing on distance gate data estimation non-working side battle array radar clutter and truly
Value mismatches, so as to cause to suppress under the performance of non-working side battle array radar clutter using space-time adaptive processing (STAP) technology
Drop.
Therefore, the suppression radar clutter purpose that obtained to space-time adaptive processing (STAP) technology, it is necessary to carry out phase
It should pre-process;In general roadmap is non-homogeneous to data caused by the distance dependencies of non-working side battle array radar clutter frequency spectrum
Property compensate to mitigate variation characteristic of the radar clutter with distance, and using range gate clutter to be detected as reference sample, point
Each auxiliary range gate do not closed on by treating detecting distance door carries out Doppler domain or angle Doppler domain compensates so that faces
Closely the power space-time distribution of door clutter reaches unanimity with unit to be detected, so as to eliminate the distance dependencies of radar clutter.
Existing compensation method mainly has Doppler effect correction (Doppler compensation, DC) method, angle Doppler
(angle-Doppler compensation, the ADC) method of compensation, adaptive angle Doppler effect correction (adaptive angle-
Doppler compensation, A2DC) method and based on clutter spectrum registration (Registration-based compensation,
RBC) method;The radar clutter distance that these methods reduce under airborne radar non-working side battle array configuration to a certain extent relies on
Property so that the radar clutter inhibition than directly handling well can be obtained when using space-time adaptive processing (STAP) technology
Energy;But DC methods have only carried out main-lobe clutter spectrum center compensation in Doppler domain;ADC methods carry out miscellaneous in angle Doppler domain
Spectral center compensates, and needs the systematic parameter provided using inertial navigation;Although A2DC methods are not required to know configuration parameter, and resonable
Radar clutter rejection is also fine in the case of thinking, but the covariance matrix pole of each range cell radar clutter is not in practice
It is stable, cause radar clutter rejection to decline a lot;RBC methods can realize the full remuneration of radar clutter, and in preferable feelings
Good radar clutter rejection can be obtained under condition, but traditional RBC methods are using time domain that smoothly to carry out radar miscellaneous for sub- snap
The extraction of spectral peaks so that the estimated accuracy of radar clutter spectrum peak and the robustness of RBC methods all can be to using space-time certainly
Adaptation processing (STAP) technical performance impacts;Iteration self-adapting algorithm (the Iterative adaptive proposed in 2010
Algorithm, IAA), accurate clutter spectrum estimation can be also obtained in the case of single snap.
The content of the invention
Problem be present for above prior art, it is an object of the invention to propose a kind of non-working side battle array radar clutter spectrum
The optimization method of registration, the optimization method of this kind of non-working side battle array radar clutter spectrum registration are directed to using traditional based on clutter spectrum
Precision, which is restricted, when registering (Registration-based compensation, RBC) method carries out radar clutter Power estimation asks
Topic, estimate non-working side battle array radar clutter spectrum using iteration self-adapting algorithm (IAA), it is miscellaneous that non-working side battle array radar is improved with this
Wave spectrum estimated accuracy so that its clutter evenly, and then improves after the spectrum pretreatment of RBC method registration non-working side battle arrays radar clutter
Suppress the clutter performance of non-working side battle array radar using space-time adaptive treatment technology.
To reach above-mentioned technical purpose, the present invention, which adopts the following technical scheme that, to be achieved.
A kind of optimization method of non-working side battle array radar clutter spectrum registration, comprises the following steps:
Step 1, the L range gate clutter data that non-working side battle array radar receives is obtained, is designated as X1,X2,X3…XL;Its
In, L represents the range gate total number that the clutter that non-working side battle array radar receives includes;
Step 2, each self-corresponding covariance of L range gate clutter data that non-working side battle array radar receives is calculated respectively
Matrix, it is designated as
Step 3, covariance matrix corresponding to k-th of range gate is chosenCovariance square is corresponded to as range gate to be detected
Battle array, and covariance matrix is corresponded to as reference using the range gate to be detected, calculate the l' range gate clutter data Xl'Change
Change matrix Tl';Wherein, k ∈ { 1,2 ..., L }, l' ∈ { 1,2 ..., L } and l' ≠ k;
Step 4, the l' range gate clutter data X is utilizedl'Transformation matrix Tl', to the l' range gate clutter data
Xl'Registration process is carried out, obtains the l' range gate clutter data Y after registrationl';
Step 5, k-th of range gate clutter data X is chosenkAs reference spurs, and according to the l' distance after registration
Door clutter data Yl', covariance matrix R corresponding to k-th of range gate after registration is calculatedk;
Step 6, the covariance matrix R according to corresponding to k-th of range gate after registrationk, calculate at space-time adaptive filtering
Non-working side battle array radar clutter after reason suppresses improvement factor IF (fdt), the non-working side after the space-time adaptive filtering process
Battle array radar clutter suppresses improvement factor IF (fdt) be used to weigh the registration effect that non-working side battle array radar clutter is composed;Wherein, k ∈
{ 1,2 ..., L }, L represent the range gate total number that the clutter that non-working side battle array radar receives includes.
Beneficial effects of the present invention:The inventive method in the case of single snap using iteration self-adapting algorithm (IAA,
Iterative adaptive algorithm) estimation non-working side battle array radar clutter spectrum, it is possible to increase non-working side battle array radar
Clutter spectrum estimated accuracy, and covariance is composed according to non-working side battle array radar clutter Spectral structure characteristic and non-working side battle array radar clutter
The covariance matrix of relation reconstruct non-working side battle array radar clutter between matrix so that the pretreated non-working side battle array of registration
Radar clutter evenly, and then improves the clutter performance that space-time adaptive treatment technology suppresses non-working side battle array radar.
Brief description of the drawings
The present invention is described in further detail with reference to the accompanying drawings and detailed description.
Fig. 1 is the general flow chart of the inventive method;
Fig. 2 (a) is that lower No. 80 range cells of positive side view battle array radar do not compensate the clutter spectrum schematic diagram that clutter spectrum obtains,
Fig. 2 (b) is the clutter spectrum schematic diagram that lower No. 80 range cells of positive side view battle array radar are obtained using RBC methods,
Fig. 2 (c) is the clutter spectrum schematic diagram that lower No. 80 range cells of positive side view battle array radar are obtained using IAA-RBC methods,
Fig. 2 (d) is the lower real clutter spectrum schematic diagram of No. 80 range cells of positive side view battle array radar;
Fig. 3 is IAA-RBC method clutter recognition improvement factors and traditional RBC method clutter recognition improvement factor comparison diagrams.
Embodiment
Reference picture 1, it is the general flow chart of the inventive method;A kind of optimization side of non-working side battle array radar clutter spectrum registration
Method, comprise the following steps:
Step 1, the L range gate clutter data that non-working side battle array radar receives is obtained, is designated as X1,X2,X3…XL;Its
In, L represents the range gate total number that the clutter that non-working side battle array radar receives includes.
Specifically, the L range gate clutter data X that non-working side battle array radar receives is obtained1,X2,X3…XL, the L
Range gate clutter data X1,X2,X3…XLTie up column vector in respectively MN × 1;Wherein, M represents what non-working side battle array radar received
Pulse total number, N represent the receiving channel antenna total number of non-working side battle array radar.
Step 2, each self-corresponding covariance of L range gate clutter data that non-working side battle array radar receives is calculated respectively
Matrix, it is designated as
Specifically, l ∈ { 1,2 ..., L }, L represent that the range gate that the clutter that non-working side battle array radar receives includes is always individual
Number, calculate covariance matrix corresponding to l-th of range gateSub-step be:
2.1 set normalization spatial frequency section [- 1,1] and the non-working side battle array radar of non-working side battle array radar respectively
Normalize Doppler frequency section [- 1,1], and uniform sampling carried out to normalization spatial frequency section [- 1,1] respectively,
Obtain KsIndividual sampled point;Uniform sampling is carried out to normalization Doppler frequency section [- 1,1], obtains KtIndividual sampled point;With
fs,nThe normalization spatial frequency of n-th of sampled point on the non-working side battle array radar normalization spatial frequency section of setting is represented, is used
fd,mRepresent the normalization Doppler frequency of m-th of sampled point on the non-working side battle array radar normalization Doppler frequency section of setting
Rate;Wherein, n ∈ { 1,2 ..., Ks, m ∈ { 1,2 ..., Kt, the non-working side battle array radar of setting is then normalized into spatial frequency
The non-working side battle array radar of section and setting normalization Doppler frequency section is built into non-working side battle array radar initialization space-time
Plane, it is calculated in non-working side battle array radar initialization space-time plane at the beginning of the clutter power at any one sampled point (m, n) place
Initial value isIts expression formula is:
Wherein, v (fd,m,fs,n) represent any one sampled point (m, n) in non-working side battle array radar initialization space-time plane
The space-time steering vector at place, andfs,nExpression is set
The normalization spatial frequency of n-th of sampled point, f on fixed non-working side battle array radar normalization spatial frequency sectiond,mRepresent setting
Non-working side battle array radar normalization Doppler frequency section on m-th of sampled point normalization Doppler frequency,Represent
Kronecker is accumulated, ()HRepresent conjugate transposition, ()TTransposition is represented, | | represent the operation that takes absolute value, KsRepresent setting
The sampled point number that non-working side battle array radar normalization spatial frequency section includes, KtRepresent that the non-working side battle array radar of setting is returned
The sampled point number that one change Doppler frequency section includes.
2.2 initialization:β represents default iterations, and i represents iterations, and initial value is 1.
2.3 are calculated in the initialization space-time plane of the non-working side battle array radar after ith iteration any one successively adopts
The covariance matrix at sampling point (m, n) placeInitialized with the non-working side battle array radar after ith iteration any one in space-time plane
The clutter power value at individual sampled point (m, n) place
Wherein, u ∈ { 1,2 ..., Kt, v ∈ { 1,2 ..., Ks, n ∈ { 1,2 ..., Ks, m ∈ { 1,2 ..., Kt, u ≠ m,
V ≠ n, KsRepresent the sampled point number that the non-working side battle array radar normalization spatial frequency section of setting includes, KtRepresent setting
The sampled point number that non-working side battle array radar normalization Doppler frequency section includes.
Specifically, the non-working side battle array radar clutter power distribution matrix after ith iterationExpression formula be:
Wherein, diag () represents diagonal matrix, i.e. diagonal element isTo angular moment
Battle array,Any one sample point is miscellaneous in non-working side battle array radar initialization space-time plane after the i-th -1 time iteration of expression
Wave power value.
In order to calculate any one sampled point in the initialization space-time plane of the non-working side battle array radar after ith iteration (m,
N) the clutter power value at placeCalculate first any one in the non-working side battle array radar initialization space-time plane after ith iteration
The covariance matrix at individual sampled point (m, n) placeIts expression formula is:
Wherein, u ∈ { 1,2 ..., Kt, v ∈ { 1,2 ..., Ks, n ∈ { 1,2 ..., Ks, m ∈ { 1,2 ..., Kt, u ≠ m,
V ≠ n, KsRepresent the sampled point number that the non-working side battle array radar normalization spatial frequency section of setting includes, KtRepresent setting
The sampled point number that non-working side battle array radar normalization Doppler frequency section includes,Represent the non-positive side after ith iteration
Depending on the clutter power value at any one sampled point (u, v) place in battle array radar initialization space-time plane, v (fd,u,fs,v) represent anon-normal
The space-time steering vector at side view battle array radar initialization space-time plane up-sampling point (u, v) place, fd,uRepresent the non-working side battle array of setting
The normalization Doppler frequency of u-th of sampled point, f on radar normalization Doppler frequency sections,vRepresent the non-working side of setting
The normalization spatial frequency of v-th of sampled point, () on battle array radar normalization spatial frequency sectionHRepresent conjugate transposition.
Then, according to any one sampled point in the non-working side battle array radar initialization space-time plane after ith iteration (m,
N) covariance matrix at placeAny one is adopted in non-working side battle array radar initialization space-time plane after calculating ith iteration
The clutter power value at sampling point (m, n) placeIts expression formula is:
Wherein, v (fd,m,fs,n) represent any one sampled point (m, n) in non-working side battle array radar initialization space-time plane
The space-time steering vector at place,Represent that any one is adopted in the initialization space-time plane of the non-working side battle array radar after ith iteration
The covariance matrix at sampling point (m, n) place, fs,nRepresent on the non-working side battle array radar normalization spatial frequency section of setting n-th
The normalization spatial frequency of sampled point, fd,mRepresent on the non-working side battle array radar normalization Doppler frequency section of setting m-th
The normalization Doppler frequency of sampled point, n ∈ { 1,2 ..., Ks, m ∈ { 1,2 ..., Kt, KsRepresent the non-working side battle array of setting
The sampled point number that radar normalization spatial frequency section includes, KtRepresent the non-working side battle array radar normalization Doppler of setting
The sampled point number that frequency separation includes.
If 2.4Or i<β, make i add 1, return to sub-step 2.3;
IfOr i >=β, then iteration terminate, now obtain non-working side battle array radar initialization space-time
Plane up-samples the space-time steering vector v (f at point (1,1) placed,1,fs,1Adopted in)-non-working side battle array radar initialization space-time plane
Sampling point (Kt,Ks) place space-time steering vectorAnd covariance matrix corresponding to l-th of range gate is calculated
Wherein, WithRespectively Kt×KsTie up matrix.
Specifically, ifOr i<β, make i add 1, return to sub-step 2.3;
IfOr i >=β, then iteration terminate, now obtain non-working side battle array radar initialization space-time
Plane up-samples the space-time steering vector v (f at point (1,1) placed,1,fs,1Adopted in)-non-working side battle array radar initialization space-time plane
Sampling point (Kt,Ks) place space-time steering vectorAnd association side corresponding to l-th of range gate is calculated
Poor matrixWherein, set according to the convergence rate of the power spectrum l1 norm differences of adjacent iteration twice, WithRespectively Kt×KsTie up matrix, l ∈ { 1,2 ..., L }, L tables
Show the range gate total number that the clutter that non-working side battle array radar receives includes.
Covariance matrix corresponding to l-th of range gateExpression formula is:
Wherein, v (fd,m,fs,n) represent in non-working side battle array radar initialization space-time plane any one sampled point to (m,
N) the space-time steering vector at place,Represent the non-working side battle array radar clutter power distribution matrix obtained during iteration stopping;ε tables
Show the constant to go to zero of determination, β represents default iterations.
Step 3, covariance matrix corresponding to k-th of range gate is chosenCovariance square is corresponded to as range gate to be detected
Battle array, and covariance matrix is corresponded to as reference using the range gate to be detected, utilize minimum 2- normsMeter
Calculate the l' range gate clutter data Xl'Transformation matrix Tl';Wherein, k ∈ { 1,2 ..., L }, l' ∈ { 1,2 ..., L } and l' ≠
K, L represent the range gate total number that the clutter that non-working side battle array radar receives includes.
Specifically, the l' range gate clutter data Xl'Transformation matrix Tl', its expression formula is:
Wherein, VkRepresentCharacteristic vector, ΛkRepresentEigenvalue matrix, Vl'RepresentCharacteristic vector, Λl'
RepresentEigenvalue matrix,Covariance matrix corresponding to l-th of range gate is represented,Represent corresponding to the l' range gate
Covariance matrix,Represent covariance matrix corresponding to k-th of range gate, l ∈ { 1,2 ..., L }, k ∈ { 1,2 ..., L }, l' ∈
{ 1,2 ..., L } and l' ≠ k, L represent the range gate total number that the clutter that non-working side battle array radar receives includes.
Step 4, the l' range gate clutter data X is utilizedl'Transformation matrix Tl', to the l' range gate clutter data
Xl'Registration process is carried out, obtains the l' range gate clutter data Y after registrationl';Wherein, l' ∈ { 1,2 ..., L }, L are represented
The range gate total number that the clutter that non-working side battle array radar receives includes.
Specifically, the l' range gate clutter data Y after the registrationl', its expression formula is:Yl'=Tl'Xl';Wherein,
Tl'Represent the l' range gate clutter data Xl'Transformation matrix, l' ∈ { 1,2 ..., L }, L represent non-working side battle array radar connect
The range gate total number that the clutter received includes.
Step 5, k-th of range gate clutter data X is chosenkAs reference spurs, and according to the l' distance after registration
Door clutter data Yl', covariance matrix R corresponding to k-th of distance gate data after registration is calculatedk;Wherein, k ∈ 1,
2 ..., L }, L represents the range gate total number that the clutter that non-working side battle array radar receives includes.
Specifically, covariance matrix R corresponding to k-th of range gate after the registrationk, its expression formula is:
Wherein, Yl'Represent the l' range gate clutter data after registration, XkRepresent k-th of range gate clutter data, k ∈
{ 1,2 ..., L }, l' ∈ { 1,2 ..., L } and l' ≠ k, L represent the range gate that the clutter that non-working side battle array radar receives includes
Total number, ()HRepresent conjugate transposition.
Covariance matrix R corresponding to k-th of range gate after registration is calculatedk, can carry out weighing training and for sky
When the processing of adaption radar clutter recognition.
Step 6, the covariance matrix R according to corresponding to k-th of range gate after registrationk, calculate at space-time adaptive filtering
Non-working side battle array radar clutter after reason suppresses improvement factor IF (fdt), the non-working side after the space-time adaptive filtering process
Battle array radar clutter suppresses improvement factor IF (fdt) be used to weigh the registration effect that non-working side battle array radar clutter is composed;Wherein, k ∈
{ 1,2 ..., L }, L represent the range gate total number that the clutter that non-working side battle array radar receives includes.
Specifically, the covariance matrix R according to corresponding to k-th of range gate after the registrationk, calculate space-time adaptive filtering
Non-working side battle array radar clutter after processing suppresses improvement factor IF (fdt), the non-positive side after the space-time adaptive filtering process
Suppress improvement factor IF (f depending on battle array radar clutterdt) be used to weigh the registration effect that non-working side battle array radar clutter is composed.
Non-working side battle array clutter recognition improvement factor IF (f after the space-time adaptive filtering processdt), its expression formula
For:
Wherein, v (fdt,fst) the moving-target space-time that represents to include in the clutter that receives of non-working side battle array radar is oriented to arrow
Amount, fdtRepresent the moving-target normalization Doppler frequency included in the clutter that non-working side battle array radar receives, fstRepresent anon-normal
The moving-target normalization spatial frequency included in the clutter that side view battle array radar receives, PCRepresent the non-working side battle array radar set
Clutter input power, PNRepresent the non-working side battle array radar noise input power set.
The effect of the present invention can be further illustrated by following emulation experiment.
(1) emulation experiment data explanation
Contrasted in order to suppress non-working side battle array radar clutter performance with traditional RBC methods, method emulation of the present invention
Using even linear array;Appropriate repetition is selected, does not consider range ambiguity problem, the miscellaneous noise ratio of non-working side battle array radar array element level is
30dB, No. 80 range cell (19km) is handled.Time domain and the smooth skill of spatial domain sub-aperture is respectively adopted in traditional RBC methods
Art obtains multiple samples, and the spatially and temporally sub-aperture in method emulation of the present invention is respectively Ns=8 and Mt=6, what sub-aperture was brought
Lose as 10*log10 (MN/MtNs), non-working side battle array radar simulation parameter is as shown in table 1.
Table 1
(2) simulation result and analysis
The simulation result of the present invention is as shown in Fig. 2 (a)~Fig. 2 (d) and Fig. 3;Wherein, Fig. 2 (a) is positive side view battle array radar
Lower No. 80 range cells do not compensate the clutter spectrum schematic diagram that clutter spectrum obtains, and Fig. 2 (b) is that lower No. 80 distances of positive side view battle array radar are single
The clutter spectrum schematic diagram that member is obtained using RBC methods, Fig. 2 (c) are that lower No. 80 range cells of positive side view battle array radar use IAA-RBC methods
Obtained clutter spectrum schematic diagram, Fig. 2 (d) are the lower real clutter spectrum schematic diagram of No. 80 range cells of positive side view battle array radar, and Fig. 3 is
IAA-RBC method clutter recognition improvement factors and traditional based on clutter spectrum registration (Registration-based
Compensation, RBC) method clutter recognition improvement factor comparison diagram.
Due to sample number deficiency under non-homogeneous scene, uncompensated radar clutter spectrum secondary lobe is high, resolution ratio is poor;Use RBC
Method and the frequency spectrum that high-resolution can be obtained using IAA-RBC methods, but RBC methods are because time domain smoothly reduces dimension, from
And make it that the precision of spectrum estimation is restricted, such as clutter spectrum widening etc.;The IAA-RBC methods that the inventive method uses pass through repeatedly
Covariance matrix is reconstructed for method, smoothly caused time domain is lost in the absence of time domain, and the precision of Power estimation is higher.
In summary, emulation experiment demonstrates the correctness of the present invention, validity and reliability.
Obviously, those skilled in the art can carry out the essence of various changes and modification without departing from the present invention to the present invention
God and scope;So, if these modifications and variations of the present invention belong to the scope of the claims in the present invention and its equivalent technologies
Within, then the present invention is also intended to comprising including these changes and modification.
Claims (8)
1. a kind of optimization method of non-working side battle array radar clutter spectrum registration, it is characterised in that comprise the following steps:
Step 1, the L range gate clutter data that non-working side battle array radar receives is obtained, is designated as X1, X2, X3…XL;Wherein, L tables
Show the range gate total number that the clutter that non-working side battle array radar receives includes;
Step 2, to calculate the L range gate clutter data that non-working side battle array radar receives respectively using iteration self-adapting algorithm each
Self-corresponding covariance matrix, is designated as
Step 3, covariance matrix corresponding to k-th of range gate is chosenCovariance matrix is corresponded to as range gate to be detected,
And covariance matrix is corresponded to as reference using the range gate to be detected, calculate the individual range gate clutter data x of l 'l′Conversion square
Battle array Tl′;Wherein, k ∈ { 1,2 ..., L }, l ' ∈ { 1,2 ..., L } and l ' ≠ k;
Step 4, the individual range gate clutter data x of l ' are utilizedl′Transformation matrix Tl′, range gate clutter data x individual to l 'l′Enter
Row registration process, obtain the individual range gate clutter data Y of the l ' after registrationl′;
Step 5, k-th of range gate clutter data x is chosenkAs reference spurs, and according to the individual range gate clutters of l ' after registration
Data Yl′, covariance matrix R corresponding to k-th of range gate after registration is calculatedk;
Step 6, the covariance matrix R according to corresponding to k-th of range gate after registrationk, after calculating space-time adaptive filtering process
Non-working side battle array radar clutter suppresses improvement factor IF (fdt), the non-working side battle array radar after the space-time adaptive filtering process
Clutter recognition improvement factor IF (fdt) be used to weigh the registration effect that non-working side battle array radar clutter is composed;Wherein, k ∈ 1,2 ...,
L }, L represents the range gate total number that the clutter that non-working side battle array radar receives includes.
2. the optimization method of a kind of non-working side battle array radar clutter spectrum registration as claimed in claim 1, it is characterised in that in step
In rapid 1, the L range gate clutter data is respectively that column vector is tieed up in MN × 1;Wherein, M represents that non-working side battle array radar receives
Pulse total number, N represent non-working side battle array radar receiving channel antenna total number.
3. the optimization method of a kind of non-working side battle array radar clutter spectrum registration as claimed in claim 1, it is characterised in that in step
In rapid 2, each self-corresponding covariance matrix of L range gate clutter data that the non-working side battle array radar receives, in addition to;
Calculate covariance matrix corresponding to l-th of range gateIts sub-step is:
2.1 respectively set non-working side battle array radar normalization spatial frequency section and non-working side battle array radar normalization it is how general
Frequency separation is strangled, and uniform sampling is carried out to the normalization spatial frequency section respectively, obtains KsIndividual sampled point, to described
Normalize Doppler frequency section and carry out uniform sampling, obtain KtIndividual sampled point, then the non-working side battle array radar of setting is returned
One change spatial frequency section and the non-working side battle array radar normalization Doppler frequency section of setting are built into non-working side battle array thunder
Up to initialization space-time plane, any one sampled point (m, n) place in non-working side battle array radar initialization space-time plane is calculated
Clutter power initial value
Wherein, n ∈ { 1,2 ..., Ks, m ∈ { 1,2 ..., Kt, l ∈ { 1,2 ..., L }, L represent that non-working side battle array radar receives
The range gate total number that includes of clutter, KsWhat the non-working side battle array radar normalization spatial frequency section for representing to set included adopts
Sampling point number, KtRepresent the sampled point number that the non-working side battle array radar normalization Doppler frequency section of setting includes;
2.2 initialization:β represents default iterations, and i represents iterations, and initial value is 1;
2.3 are calculated any one sampled point in the initialization space-time plane of the non-working side battle array radar after ith iteration successively
The covariance matrix at (m, n) placeWith after ith iteration non-working side battle array radar initialization space-time plane on any one adopt
The clutter power value at sampling point (m, n) place
If 2.4Or i < β, make i add 1, return to sub-step 2.3;
IfOr i >=β, then iteration terminate, now obtain non-working side battle array radar initialization space-time plane
Up-sample the space-time steering vector v (f at point (1,1) placeD, 1, fS, 1)-non-working side battle array radar initialization space-time plane up-sampling point
(Kt, Ks) place space-time steering vectorAnd covariance matrix corresponding to l-th of range gate is calculatedIts
In, WithRespectively Kt×KsMatrix is tieed up, ε represents that what is determined becomes
In zero constant, β represents default iterations,RepresentNorm.
4. the optimization method of a kind of non-working side battle array radar clutter spectrum registration as claimed in claim 3, it is characterised in that described
The clutter power initial value at any one sampled point (m, n) place in non-working side battle array radar initialization space-time planeDescribed
The covariance matrix at any one sampled point (m, n) place in non-working side battle array radar initialization space-time plane after i iterationAny one sampled point (m, n) place is miscellaneous in non-working side battle array radar initialization space-time plane after the ith iteration
Wave power valueWith l-th of range gate corresponding to covariance matrixIts expression formula is respectively:
Wherein, v (fD, m, fS, n) represent any one sampled point (m, n) place in non-working side battle array radar initialization space-time plane
Space-time steering vector, fS, nRepresent the normalizing of n-th of sampled point on the non-working side battle array radar normalization spatial frequency section of setting
Change spatial frequency, fD, mRepresent the normalizing of m-th of sampled point on the non-working side battle array radar normalization Doppler frequency section of setting
Change Doppler frequency, ()HConjugate transposition is represented, | | the operation that takes absolute value is represented,Represent non-after ith iteration
The clutter power value at any one sampled point (u, v) place, v (f in positive side view battle array radar initialization space-time planeD, u, fS, v) represent
The space-time steering vector at non-working side battle array radar initialization space-time plane up-sampling point (u, v) place, fD, uRepresent the non-positive side of setting
Depending on the normalization Doppler frequency of u-th of sampled point on battle array radar normalization Doppler frequency section, fS, vRepresent the anon-normal of setting
The normalization spatial frequency of v-th of sampled point on side view battle array radar normalization spatial frequency section,When representing iteration stopping
Obtained non-working side battle array radar clutter power distribution matrix, ε represent the constant to go to zero determined, and β represents default iteration time
Number, u ∈ { 1,2 ..., Kt, v ∈ { 1,2 ..., Ks, n ∈ { 1,2 ..., Ks, m ∈ { 1,2 ..., Kt, u ≠ m, v ≠ n, KsRepresent
The sampled point number that the non-working side battle array radar normalization spatial frequency section of setting includes, KtRepresent the non-working side battle array of setting
The sampled point number that radar normalization Doppler frequency section includes.
5. the optimization method of a kind of non-working side battle array radar clutter spectrum registration as claimed in claim 1, it is characterised in that in step
In rapid 3, the individual range gate clutter data x of l 'l′Transformation matrix Tl′, its expression formula is:
Wherein, VkRepresentCharacteristic vector, ΛkRepresentEigenvalue matrix, Vl′RepresentCharacteristic vector, Λl′RepresentEigenvalue matrix,Covariance matrix corresponding to l-th of range gate is represented,Represent association side corresponding to the individual range gates of l '
Poor matrix,Represent covariance matrix corresponding to k-th of range gate, l ∈ { 1,2 ..., L }, k ∈ { 1,2 ..., L }, l ' ∈ 1,
2 ..., L } and l ' ≠ k, L represent the range gate total number that the clutter that receives of non-working side battle array radar includes.
A kind of 6. optimization method of non-working side battle array radar clutter spectrum registration as claimed in claim 1, it is characterised in that
In step 4, the individual range gate clutter data Y of l ' after the registrationl′, its expression formula is:Yl′=Tl′Xl′;Wherein, Tl′Represent
The individual range gate clutter data x of l 'l′Transformation matrix, l ' ∈ { 1,2 ..., L }, L represent non-working side battle array radar receive it is miscellaneous
The range gate total number that ripple includes.
7. the optimization method of a kind of non-working side battle array radar clutter spectrum registration as claimed in claim 1, it is characterised in that in step
In rapid 5, covariance matrix R corresponding to k-th of range gate after the registrationk, its expression formula is:
Wherein, Yl′Represent the individual range gate clutter datas of l ' after registration, XkRepresent k-th of range gate clutter data, k ∈ 1,
2 ..., L }, the range gate that the clutter that l ' ∈ { 1,2 ..., L } and l ' ≠ k, L expression non-working side battle array radar receive includes is always individual
Number, ()HRepresent conjugate transposition.
8. the optimization method of a kind of non-working side battle array radar clutter spectrum registration as claimed in claim 1, it is characterised in that in step
In rapid 6, the non-working side battle array radar clutter after the space-time adaptive filtering process suppresses improvement factor IF (fdt), its expression formula
For:
Wherein, v (fdt, fst) represent the moving-target space-time steering vector that is included in the clutter that receives of non-working side battle array radar, fdt
Represent the moving-target normalization Doppler frequency included in the clutter that non-working side battle array radar receives, fstRepresent non-working side battle array
The moving-target normalization spatial frequency included in the clutter that radar receives, PCRepresent that the non-working side battle array radar clutter of setting is defeated
Enter power, PNRepresent the non-working side battle array radar noise input power set.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610156725.9A CN105785326B (en) | 2016-03-18 | 2016-03-18 | A kind of optimization method of non-working side battle array radar clutter spectrum registration |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610156725.9A CN105785326B (en) | 2016-03-18 | 2016-03-18 | A kind of optimization method of non-working side battle array radar clutter spectrum registration |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105785326A CN105785326A (en) | 2016-07-20 |
CN105785326B true CN105785326B (en) | 2018-03-06 |
Family
ID=56394063
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610156725.9A Active CN105785326B (en) | 2016-03-18 | 2016-03-18 | A kind of optimization method of non-working side battle array radar clutter spectrum registration |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105785326B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106772253B (en) * | 2016-11-25 | 2019-12-24 | 西安电子科技大学 | Radar clutter suppression method under non-uniform clutter environment |
CN106970358B (en) * | 2017-04-28 | 2019-12-24 | 西安电子科技大学 | Optimization method for angular Doppler registration of clutter spectrum of non-normal side-looking array radar |
CN107092012A (en) * | 2017-05-19 | 2017-08-25 | 中国人民解放军空军工程大学 | Compensated distance method and apparatus |
CN108020817A (en) * | 2017-09-28 | 2018-05-11 | 西安电子科技大学 | Air-borne Forward-looking battle array radar clutter suppression method based on registration |
CN108761417B (en) * | 2018-05-31 | 2021-12-10 | 西安电子科技大学 | Airborne radar clutter suppression method based on knowledge-aided maximum likelihood |
CN112904289B (en) * | 2021-01-18 | 2022-12-02 | 西安电子科技大学 | Airborne non-forward side looking array radar steady clutter suppression method based on diagonal loading |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103969629A (en) * | 2014-04-21 | 2014-08-06 | 河海大学 | Airborne radar clutter self-adaption restraining method based on main-lobe clutter registering |
CN104345300A (en) * | 2014-10-30 | 2015-02-11 | 河海大学 | Onboard non-positive side view array radar STAP (Space Time Adaptive Processing) method for clutter space-time spectrum linear compensation |
-
2016
- 2016-03-18 CN CN201610156725.9A patent/CN105785326B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103969629A (en) * | 2014-04-21 | 2014-08-06 | 河海大学 | Airborne radar clutter self-adaption restraining method based on main-lobe clutter registering |
CN104345300A (en) * | 2014-10-30 | 2015-02-11 | 河海大学 | Onboard non-positive side view array radar STAP (Space Time Adaptive Processing) method for clutter space-time spectrum linear compensation |
Non-Patent Citations (4)
Title |
---|
"Clutter suppression approach for non-sidelooking airborne radar with medium pulse repetition frequency";Baoquan Dai et al.;《IET Radar, Sonar and Navigation》;20140930;第8卷(第9期);第999-1007页 * |
"Clutter suppression for airborne non-sidelooking radar using ERCB-STAP algorithm";J. Wu et al.;《IET Radar, Sonar and Navigation》;20100430;第4卷(第4期);第497-506页 * |
"基于主瓣杂波高效配准的机载非正侧视阵雷达STAP算法研究";王杰 等;《雷达学报》;20140430;第3卷(第2期);第235-240页 * |
"基于稀疏采样的双基地机载雷达杂波谱补偿方法";刘锦辉 等;《系统工程与电子技术》;20120630;第34卷(第6期);第1120-1125页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105785326A (en) | 2016-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105785326B (en) | A kind of optimization method of non-working side battle array radar clutter spectrum registration | |
CN103926572B (en) | A kind of clutter suppression method of battle array radar self adaptation subspace, airborne anon-normal side | |
CN105738879B (en) | Radar clutter space-time adaptive pre-filtering method based on sparse recovery | |
CN105929371B (en) | A kind of airborne radar clutter suppression method based on covariance matrix | |
CN105223560B (en) | Airborne radar object detection method based on the sparse recovery of clutter pitching azimuth spectrum | |
US7064702B1 (en) | System, method and computer program product for reducing quadratic phase errors in synthetic aperture radar signals | |
CN110113085B (en) | Wave beam forming method and system based on covariance matrix reconstruction | |
CN103616687B (en) | The fitting of a polynomial ISAR envelope alignment method that piecewise linearity is estimated | |
He et al. | Processing of ultrahigh-resolution spaceborne sliding spotlight SAR data on curved orbit | |
CN104345300B (en) | The airborne non-working side battle array radar STAP method of clutter space-time spectrum linear compensation | |
CN110221241A (en) | A kind of low elevation angle DOA estimation method based on RBF neural | |
CN106772253B (en) | Radar clutter suppression method under non-uniform clutter environment | |
CN109324322A (en) | A kind of direction finding and target identification method based on passive phased array antenna | |
CN105022040A (en) | Array element error estimation method based on clutter data combined fitting | |
CN103983959A (en) | SAR system movement target radial speed estimation method based on data reconstruction | |
CN101881822A (en) | Method for inhibiting same frequency interference of shared-spectrum radars | |
CN113050048A (en) | Orthogonal waveform optimization design method of LFM-PC composite modulation signal | |
CN104215939B (en) | Knowledge assisted space-time adaptive processing method integrating generalized symmetrical structure information | |
CN106970358B (en) | Optimization method for angular Doppler registration of clutter spectrum of non-normal side-looking array radar | |
Sun et al. | A postmatched-filtering image-domain subspace method for channel mismatch estimation of multiple azimuth channels SAR | |
CN106125039B (en) | Improvement space-time adaptive Monopulse estimation method based on local Combined Treatment | |
CN105044688A (en) | Radar robust space-time adaption processing method based on iterative subspace tracking algorithm | |
Wang et al. | Clutter suppression and ground moving target imaging approach for hypersonic vehicle borne multichannel radar based on two-step focusing method | |
CN113376606B (en) | Clutter suppression method for quickly converging sparse Bayes along clutter ridge | |
CN109581351A (en) | The method that joint pixel pait normalizes sample covariance matrix estimation target radial speed |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |