CN102064892A - Two-level interference suppression method for sub-array level adaptive single pulse - Google Patents

Two-level interference suppression method for sub-array level adaptive single pulse Download PDF

Info

Publication number
CN102064892A
CN102064892A CN2010106089947A CN201010608994A CN102064892A CN 102064892 A CN102064892 A CN 102064892A CN 2010106089947 A CN2010106089947 A CN 2010106089947A CN 201010608994 A CN201010608994 A CN 201010608994A CN 102064892 A CN102064892 A CN 102064892A
Authority
CN
China
Prior art keywords
delta
sub
main lobe
level
interference suppression
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.)
Pending
Application number
CN2010106089947A
Other languages
Chinese (zh)
Inventor
胡航
阎晓莉
徐颖
肖勇
邱朝阳
吴群
张皓
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN2010106089947A priority Critical patent/CN102064892A/en
Publication of CN102064892A publication Critical patent/CN102064892A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a two-level interference suppression method for a sub-array level adaptive single pulse, and relates to an interference suppression method for a sub-array level adaptive single pulse. The method solves the problems of poorer single pulse performance and higher calculation cost of the conventional interference suppression method for the sub-array level adaptive single pulse. The method comprises two levels of interference suppression: the first level of interference suppression is side lobe interference suppression with main lobe shape keeping, and the interference suppression method comprises the following steps of: combining sub-array level adaptive digital beam forming (ADBF) with directional diagram control performance and the main lobe shape keeping, and suppressing the side lobe interference at the same time of keeping main lobe shape; the second level of interference suppression is main lobe interference suppression, and the interference suppression method comprises the following step of: performing main lobe interference suppression by adopting a 4-channel single pulse system; and the 4 channels of the 4-channel single pulse system mean 4 receiving channels for forming 4 beams respectively, namely a sum beam, a pitching difference beam, an azimuth difference beam and a double difference beam. The method is suitable for interference suppression of the sub-array level adaptive single pulse.

Description

The two-stage disturbance restraining method of submatrix level adaptation pulse
Technical field
The present invention relates to a kind of disturbance restraining method of submatrix level adaptation pulse.
Background technology
Array antenna is made of by certain regularly arranged antenna a plurality of.Phase array is a kind of crucial array antenna form, has important application in fields such as communication, direction-finding station, smart antenna, seismic prospectings.It is compared tool with traditional mechanical scanning method and has great advantage by the beam position of phase shifter control directional diagram.It need not antenna is carried out mechanical rotation, has improved beam scanning speed greatly, has the advantage of quick discovery target; Even under target-rich environment, also can follow the tracks of a plurality of targets fast by changing phase shifter.Phase array utilizes the phase difference between the array element to carry out direction finding, can obtain narrow beam by increasing array aperture, thereby improve direction finding precision greatly.
After connecing a digital receive path (comprising mixing, amplification, filtering and A/D conversion etc.) after each array element of phase array, just obtain the phase array of digital array.The phase array of digital array can adopt array signal process technique, and many advanced persons' array processing method is used, thereby the performance of phased array system is greatly improved.The phase array of digital array has the incomparable superiority of traditional phase array, and the information that has improved phased array system greatly detects and handling property.
Array number is all a lot of in a lot of phased array systems, usually reaches hundreds of to thousands of, and all can't bear cost and operational capability on if still adopt receive path with the array number equivalent amount this moment.Usually adopt subarray configuration for this reason, some adjacent array elements are combined as a submatrix, after each submatrix, connect a receive path and carry out digitlization, thereby can reduce the array cost greatly and reduce the dimension of signal processing.
Yet, behind the employing subarray configuration, can't obtain the digitlization output of each array element, can only obtain the digitlization output of each submatrix, thereby conventional array element level array processing method is no longer suitable, must employing submatrix level array processing method.It is significantly different that submatrix level ARRAY PROCESSING and conventional array element level ARRAY PROCESSING have, and exists many submatrix levels to handle specific question to be solved.
A lot of phased array systems adopt monopulse technologies to carry out angle measurement, it based on the output of, difference beam, receive echo-signal simultaneously with the several separate receiving branch, the signal by relatively each branch road reception is to obtain angle information; Only utilize an echo impulse just can determine target direction.
In the practical application, electronic system may be subjected to various interference; Accuracy for target detection and angle in the bonding pulse technique are estimated must suppress it.Can adopt adaptive technique for this reason; The monopulse technology of having used adaptive technique is called as the self adaptation pulse, and it is the key technology that pulse is estimated.Requirement to the self adaptation pulse is, when effectively suppressing to disturb, should keep good pulse performance (the pulse ratio: promptly difference beam with the ratio of the directivity function of wave beam), with the accuracy that guarantees that angle is estimated.
Under the interference environment, ADBF (the adaptive digital wave beam forms, i.e. Adaptive Digital Beamforming) can effectively improve Signal to Interference plus Noise Ratio, recovers the good detection performance.For monopulse technology, secondary lobe can adopt ADBF under disturbing; Because zero the falling into of self adaptation this moment is positioned at the secondary lobe zone, the distortion that the directional diagram main beam produces is very little.But under disturbing, main lobe can not adopt ADBF, be positioned near the beam position because disturb this moment, self adaptation zero falls into the main beam that will make directional diagram and produces very big distortion, thereby make the self adaptation pulse than and static state answer pulse to compare to produce very mistake, the serious deviation that causes angle to be estimated than when noiseless (pulse than).Main lobe disturbs down, how to keep good pulse performance when suppressing interference, is the key issue that the self adaptation pulse should solve.
Self adaptation single pulse method under abroad main lobe being disturbed has been carried out a lot of researchs, and its overwhelming majority provides certain special pulse processing method, and itself and adaptive direction figure are taken into account.Even the main beam characteristic of adaptive direction figure matches with pulse, avoid the pulse performance near the decline of main beam zero point.These methods comprise: 1, maximum likelihood (Maximumlikelihood is abbreviated as ML) method, at linear array 3 kinds of different maximum likelihood pulse expression formulas are proposed as R.C.Davies etc., and on statistical significance, solved the angle estimation problem; U.Nickel estimates to promote to the maximum likelihood of self adaptation and wave beam, obtains two-dimentional monopulse technology, and proposes 1 rank Taylor expansion of self adaptation pulse; 2. least variance method.Promptly carry out the target direction search, to determine optimum difference beam weights adaptively based on least mean-square error.Yet said method all need utilize and, difference passage output, be the pulse ratio that obtains proofreading and correct, need very high computing cost.R.L.Fante proposed constraint self adaptation single pulse method in 1999, suppress the angular error of pulse by the self adaptation difference beam is increased constraints, but reduced degree of freedom in system, and output signal-to-noise ratio is descended to some extent.
K.B.Yu equals calendar year 2001 and proposes a kind of secondary lobe to be disturbed and main lobe disturbs the self adaptation single pulse method that suppresses respectively: interference suppresses to secondary lobe at first to utilize the ADBF technology, carries out the main lobe conformal simultaneously; And then main lobe disturbed suppress, keep the pulse performance simultaneously.Yet there are three deficiencies in this method: the applied environment that 1, is not suitable for submatrix level phase array; 2, after the 1st level adaptation was handled, the secondary lobe of adaptive direction figure increased greatly; 3, the pulse performance is undesirable, and the self adaptation pulse is than having bigger distortion departing from the beam position place, and its reason is that the main lobe conformal effect of the 1st level adaptation in handling is bad.
Summary of the invention
The present invention is for the pulse poor-performing of the disturbance restraining method that solves existing submatrix level adaptation pulse, computing cost problem of higher, thereby a kind of two-stage disturbance restraining method of submatrix level adaptation pulse is provided.
The two-stage disturbance restraining method of submatrix level adaptation pulse is characterized in that: it comprises that the two-stage interference suppresses:
The first order disturbs inhibition to suppress for the secondary lobe with main lobe conformal disturbs, and its method is: combine with the main lobe conformal by the submatrix level ADBF that will have the directional diagram control performance, when keeping the main lobe shape, interference suppresses to secondary lobe;
The second level is disturbed and is suppressed to suppress for main lobe disturbs, and its method is: adopt 4 passage monopulse systems to carry out the main lobe interference and suppress; The implication of 4 passages of described 4 passage monopulse systems is: each passage forms 1 wave beam in 4 receive paths, totally 4 wave beams, that is: and wave beam, trim wave beam, gun parallax wave beam and two difference beam.
The concrete grammar of first order disturbance restraining method is: when not having interference, introduce 1 each element and be 1 static dominant vector, construct the mismatch guiding vector simultaneously, and the guiding vector that the submatrix level optimal beam of main lobe conformal forms device is revised with the mismatch guiding vector, promptly be the product of 1 static dominant vector with the covariance matrix of the interference plus noise of submatrix level output under the noiseless situation and each element, replace the guiding vector of submatrix level phase array at the beam position place, thereby the submatrix level adaptation of the main lobe conformal that obtains revising power is finished the first order and is disturbed inhibition.
In the disturbance restraining method of the second level, if the orientation to pitching to directional diagram separate, in 4 passage monopulse systems, suppress interference and form self adaptation zero to fall into along a direction, along another orthogonal direction keep non-self-adapting and, difference beam, thereby obtain undistorted pulse ratio, promptly with static pulse than identical.
In 4 passage monopulse systems, suppress to disturb and form self adaptation zero along pitch orientation and fall into, and along azimuth direction keep non-self-adapting and, difference beam, this moment, described difference beam was meant the gun parallax wave beam.
Suppress to disturb and form self adaptation zero along azimuth direction and fall into, and along pitch orientation keep non-self-adapting and, difference beam, this moment, described difference beam was meant the trim wave beam.
Beneficial effect: the present invention carries out combination with multiple technologies such as submatrix level ADBF, main lobe conformal, directional diagram control and main lobe interference inhibition, given full play to the advantage of distinct methods, effectively suppressing can to keep good pulse performance in the interference, its self adaptation pulse is than more approaching than very with static pulse; Simultaneously, the present invention adopts the two-stage adaptation scheme, secondary lobe interference and main lobe interference are suppressed respectively, need not to design complicated pulse program cooperates with the main lobe of adaptive direction figure, compare with methods such as maximum likelihood, minimum variances, need not to self adaptation and, the output of difference beam proofreaies and correct, the computing cost is reduced significantly.
Description of drawings
Fig. 1 be the four-way pulse described in the specific embodiment of the invention one constitute schematic diagram with passage; Fig. 2 is that the trim passage of the four-way pulse described in the specific embodiment of the invention one constitutes schematic diagram; Fig. 3 is that the gun parallax passage of the four-way pulse described in the specific embodiment of the invention one constitutes schematic diagram; Fig. 4 is that two difference passages of the four-way pulse described in the specific embodiment of the invention one constitute schematic diagram; Fig. 5 is the planar array coordinate system schematic diagram described in the specific embodiment of the invention one.
Embodiment
The two-stage disturbance restraining method of embodiment one, submatrix level adaptation pulse, it comprises that the two-stage interference suppresses:
The first order disturbs inhibition to suppress for the secondary lobe with main lobe conformal disturbs, and its method is: combine with the main lobe conformal by the submatrix level ADBF that will have the directional diagram control performance, when keeping the main lobe shape, interference suppresses to secondary lobe;
The second level is disturbed and is suppressed to suppress for main lobe disturbs, and its method is: adopt 4 passage monopulse systems to carry out the main lobe interference and suppress; The implication of 4 passages of described 4 passage monopulse systems is: each passage forms 1 wave beam in 4 receive paths, totally 4 wave beams, that is: and wave beam, trim wave beam, gun parallax wave beam and two difference beam.
The concrete grammar of first order disturbance restraining method is: introduce 1 each element and be 1 static dominant vector, construct the mismatch guiding vector simultaneously, and the guiding vector that the submatrix level optimal beam of main lobe conformal forms device is revised with the mismatch guiding vector, promptly be the product of 1 static dominant vector with the covariance matrix of the interference plus noise of submatrix level output under the noiseless situation and each element, replace the guiding vector of submatrix level phase array at the beam position place, thereby the submatrix level adaptation of the main lobe conformal that obtains revising power is finished the first order and is disturbed inhibition.
The main lobe conformal: main lobe disturbs down, the self adaptation of array pattern zero falls into the main beam that will make directional diagram and produces very big distortion, thereby make the pulse angle estimate to produce bigger deviation. the main lobe conformal is meant and keeps the main lobe shape invariance, promptly not forming self adaptation zero at main lobe falls into. and specific implementation is: at first estimate the main lobe interference radiating way, from the covariance matrix of the interference plus noise of submatrix output, deduct then and comprise that part that main lobe disturbs, thereby being about to main lobe disturbs the influence that produces to remove from the covariance matrix of the interference plus noise of submatrix output. do not suppressed the self adaptation power that main lobe disturbs, thereby it is sunken not form self adaptation zero in main lobe, thereby kept the main lobe shape.
Go to disturb: the signal that arrives by array received, obtain submatrix level adaptation power,, make array pattern form zero adaptively and fall into (submatrix level adaptation power is adjusted adaptively along with the difference of disturbing), thereby interference is curbed from the effect that should weigh by the submatrix level at interference radiating way.
In the disturbance restraining method of the second level, if the orientation to pitching to directional diagram separate, in 4 passage monopulse systems, suppress interference and form self adaptation zero to fall into along a direction, along another orthogonal direction keep non-self-adapting and, difference beam, thereby obtain undistorted pulse ratio, promptly with static pulse than identical.Thereby, suppressing pitching when main lobe disturbs, kept the orientation to and, the difference beam shape, promptly kept the orientation to the pulse performance; And suppressing the orientation when main lobe disturbs, kept pitching to and, the difference beam shape, promptly kept pitching to the pulse performance.Like this, when suppressing the main lobe interference, kept good pulse performance.
Operation principle: 1. four-way monopulse system and signal model
If array is positioned on the rectangular mesh on xoy plane, by N 1* N 2Individual omnidirectional array element is formed, and is d in the array element distance of x and y direction.With the formation schematic diagram of, trim, gun parallax and two difference passages respectively shown in Fig. 1 to 4, as seen from Figure 1 and the output of passage be all delayed output signals sums; As seen from Figure 2, the output of trim passage is the difference of top two quadrant array elements output sum and following two quadrant array elements output sum; As seen from Figure 3, the output of gun parallax passage is the difference of two quadrant array element output of the output sum and the left side sum of two the quadrant array elements in the right side; As seen from Figure 4, the output of two poor passages is differences of diagonal quadrant array output sum.
Coordinate system as shown in Figure 5.With θ and Represent the angle of pitch and azimuth respectively,
Figure BDA0000041153980000052
V=sin θ represents respectively
Figure BDA0000041153980000053
The unit vector of direction is in the projection of x and y direction.
With array partition is L 1* L 2Individual submatrix.If T xN for the x direction 1* L 1The dimension submatrix forms matrix, at its l 1(l 1=0 ..., L 1-1) in Lie all elements, has only and l 1The element value of the array element sequence number correspondence of individual submatrix is 1, and all the other are 0.Similarly, establish T yN for the y direction 2* L 2The dimension submatrix forms matrix.Use T 0The submatrix of representing whole array forms matrix, then has
T 0 = T x ⊗ T y - - - ( 1 )
In the formula
Figure BDA0000041153980000055
Expression Kronecker is long-pending.
For different passages, adopt different weights on the array element.To use w respectively with the array element level weighting of, trim, gun parallax and two difference beams ,
Figure BDA0000041153980000056
And
Figure BDA0000041153980000057
Expression.Definition w Tay_uFor adopting Taylor, array x direction adds temporary N 1The dimension weighing vector, w Tay_vFor adopting Taylor to add temporary N on the y direction 2The right-safeguarding vector, w Bay_uFor adopting Bayliss, the x direction adds temporary N1 right-safeguarding vector, w Bay_vFor adopting Bayliss, the y direction adds temporary N 2The right-safeguarding vector; Then have:
w Σ = w Tay _ v ⊗ w Tay _ u w Δ E = w Bay _ v ⊗ w Tay _ u w Δ A = w Tay _ v ⊗ w Bay _ u w Δ Δ = w Bay _ v ⊗ w Bay _ u - - - ( 2 )
With w ,
Figure BDA0000041153980000059
And
Figure BDA00000411539800000510
In n (n=0 ..., N 1N 2-1) individual element is used w respectively ∑ _ n,
Figure BDA00000411539800000511
And
Figure BDA00000411539800000512
Expression, and order:
W Σ = diag ( w Σ _ n ) n = 0 , . . . , N 1 · N 2 - 1 W Δ E = diag ( w Δ E - n ) n = 0 , . . . , N 1 · N 2 - 1 W Δ A = diag ( w Δ A _ n ) n = 0 , . . . , N 1 · N 2 - 1 W Δ Δ = diag ( w Δ Δ _ n ) n = 0 , . . . , N 1 · N 2 - 1 - - - ( 3 )
Order:
Figure BDA0000041153980000061
And establish:
Figure BDA0000041153980000062
Use T ,
Figure BDA0000041153980000063
With
Figure BDA0000041153980000064
Respectively expression and, the submatrix transition matrix of trim, gun parallax and two difference beams, then:
T Σ = Φ W Σ T 0 T Δ E = Φ W Δ E T 0 T Δ A = Φ W Δ A T 0 T Δ Δ = Φ W Δ Δ T 0 - - - ( 6 )
If the interference plus noise of array element output is x (t), then its covariance matrix is R=E[x (t) x H(t)].For with, trim, gun parallax and two difference beams, the interference plus noise of submatrix level output is respectively:
x sub _ Σ ( t ) = ( T Σ ) H x ( t ) x sub _ Δ E ( t ) = ( T Δ E ) H x ( t ) x sub _ Δ A ( t ) = ( T Δ A ) H x ( t ) x sub _ Δ Δ ( t ) = ( T Δ Δ ) H x ( t ) - - - ( 7 )
The covariance matrix of the submatrix level interference plus noise of each passage output is:
R sub _ Σ = ( T Σ ) H RT Σ R sub _ Δ E = ( T Δ E ) H RT Δ E R sub _ Δ A = ( T Δ A ) H RT Δ A R sub _ Δ Δ = ( T Δ Δ ) H RT Δ Δ - - - ( 8 )
2. the 1st level adaptation
Order
α ( u ) = [ 1 , . . . exp ( j 2 πd n 1 u / λ ) , . . . exp ( j 2 πd ( N 1 - 1 ) u / λ ) ] T α ( v ) = [ 1 , . . . exp ( j 2 πd n 2 v / λ ) , . . . exp ( j 2 πd ( N 2 - 1 ) v / λ ) ] T - - - ( 9 )
If
Figure BDA0000041153980000071
And order
Figure BDA0000041153980000072
° expression Hadamard product in the formula.If
Figure BDA0000041153980000073
v 0=sin θ 0, (u then 0, v 0) the array element level guiding vector of each wave beam is on the direction:
a Σ ( u 0 , v 0 ) = α Σ ( u 0 ) ⊗ a Σ ( v 0 ) a Δ E ( u 0 , v 0 ) = α Σ ( u 0 ) ⊗ a Δ ( v 0 ) a Δ A ( u 0 , v 0 ) = a Δ ( u 0 ) ⊗ a Σ ( v 0 ) a Δ Δ ( u 0 , u 0 ) a Δ ( u 0 ) ⊗ a Δ ( v 0 ) - - - ( 11 )
Thereby the submatrix level guiding vector of each wave beam is:
a sub _ Σ ( u 0 , v 0 ) = ( T Σ ) H a Σ ( u 0 , v 0 ) a sub _ Δ E ( u 0 , v 0 ) = ( T Δ E ) H a Δ E ( u 0 , v 0 ) a sub _ Δ A ( u 0 , v 0 ) = ( T Δ A ) H a Δ A ( u 0 , v 0 ) a sub _ Δ Δ ( u 0 , v 0 ) = ( T Δ Δ ) H a Δ Δ ( u 0 , v 0 ) - - - ( 12 )
For each wave beam, the self adaptation power that forms device based on the submatrix level optimal beam of LCMV criterion (Linearly Constrained Minimum Variance, linear restriction minimum variance) is:
w sub _ Σ = μ ( R sub _ Σ ) - 1 a sub _ Σ ( u 0 , v 0 ) w sub _ Δ E = μ ( R sub _ Δ E ) - 1 a sub _ Δ E ( u 0 , v 0 ) w sub _ Δ A = μ ( R sub _ Δ A ) - 1 a sub _ Δ A ( u 0 , v 0 ) w sub _ Δ Δ = μ ( R sub _ Δ Δ ) - 1 a sub _ Δ Δ ( u 0 , v 0 ) - - - ( 13 )
In the formula, μ is a constant.
Below submatrix level optimal beam being formed device combines with the main lobe conformal.If only exist a main lobe to disturb, and direction is (u MJ, v MJ).At first, utilize Estimation of Spatial Spectrum method such as submatrix level MUSIC to estimate interference radiating way, the covariance matrix of the submatrix level interference plus noise of each passage output after obtaining removing main lobe then and disturbing:
R ^ sub _ Σ = R sub _ Σ - P MJ [ a sub _ Σ ( u MJ , v MJ ) · a sub _ Σ H ( u MJ , v MJ ) ] R ^ sub _ Δ E = R sub _ Δ A - P MJ [ a sub _ Δ A ( u MJ , v MJ ) · a sub _ Δ E H ( u MJ , v MJ ) ] R ^ sub _ Δ A = R sub _ Δ A - P MJ [ a sub _ Δ A ( u MJ , v MJ ) · a sub _ Δ A H ( u MJ , v MJ ) ] R ^ sub _ Δ Δ = R sub _ Δ Δ - P MJ [ a sub _ Δ Δ ( u MJ , v MJ ) · a sub _ Δ Δ H ( u MJ , v MJ ) ] - - - ( 14 )
In the formula, P MJBe the main lobe interference power.For each wave beam, the weight vector that the submatrix level optimal beam of main lobe conformal forms device is:
w ^ sub _ Σ = μ ( R ^ sub _ Σ ) - 1 a sub _ Σ ( u 0 , v 0 ) w ^ sub _ Δ E = μ ( R ^ sub _ Δ E ) - 1 a sub _ Δ E ( u 0 , v 0 ) w ^ sub _ Δ A = μ ( R ^ sub _ Δ A ) - 1 a sub _ Δ A ( u 0 , v 0 ) w ^ sub _ Δ Δ = μ ( R ^ sub _ Δ Δ ) - 1 a sub _ Δ Δ ( u 0 , v 0 ) - - - ( 15 )
Below, submatrix level optimal beam formation device to the main lobe conformal improves, the submatrix level optimal beam that replaces the main lobe conformal with the submatrix level adaptation Beam-former of the main lobe conformal with directional diagram control performance forms device, promptly adopts the submatrix level adaptation Beam-former with directional diagram control performance of main lobe conformal.Its purpose has two: (1), suppress self adaptation and, the secondary lobe of difference directional diagram; (2), improve main lobe conformal effect.For this reason, introduce the L of 1 each element equal 1 1* L 2Dimensional vector q, and the guiding vector in the formula (12) is modified to:
a sub _ Σ ( mod i ) ( u 0 , v 0 ) = R ^ sub _ Σ ( NJ ) · q a sub _ Δ E ( mod i ) ( u 0 , v 0 ) = R ^ sub _ Δ E ( NJ ) · q a sub _ Δ E ( mod i ) ( u 0 , v 0 ) = R ^ sub _ Δ A ( NJ ) · q a sub _ Δ Δ ( mod i ) ( u 0 , v 0 ) = R ^ sub _ Δ Δ ( NJ ) · q - - - ( 16 )
In the formula,
Figure BDA0000041153980000084
Represent when noiseless (subscript (NJ)Represent noiseless), the rest may be inferred by analogy.With the guiding vector in the following formula alternate form (15), the weight vector that obtains the submatrix level adaptation Beam-former with directional diagram control performance of main lobe conformal is:
w ^ sub _ Σ ( mod i ) = μ ( R ^ sub _ Σ ) - 1 R ^ sub _ Σ ( NJ ) q w ^ sub _ Δ E ( mod i ) = μ ( R ^ sub _ Δ E ) - 1 R ^ sub _ Δ E ( NJ ) q w ^ sub _ Δ A ( mod i ) = μ ( R ^ sub _ Δ A ) - 1 R ^ sub _ Δ A ( NJ ) q w ^ sub _ Δ Δ ( mod i ) = μ ( R ^ sub _ Δ Δ ) - 1 R ^ sub _ Δ Δ ( NJ ) q - - - ( 17 )
Do not having under the situation about disturbing, following formula becomes:
( w ^ sub _ Σ ( mod i ) ) ( NJ ) = μ ( R ^ sub _ Σ ( NJ ) ) - 1 R ^ sub _ Σ ( NJ ) q = μq ( w ^ sub _ Δ E ( mod i ) ) ( NJ ) = μ ( R ^ sub _ Δ E ( NJ ) ) - 1 R ^ sub _ Δ E ( NJ ) q = μq ( w ^ sub _ Δ A ( mod i ) ) ( NJ ) = μ ( R ^ sub _ Δ A ( NJ ) ) - 1 R ^ sub _ Δ A ( NJ ) q = μq ( w ^ sub _ Δ Δ ( mod i ) ) ( NJ ) = μ ( R ^ sub _ Δ Δ ( NJ ) ) - 1 R ^ sub _ Δ Δ ( NJ ) q = μq - - - ( 18 )
At this moment, for each wave beam, each submatrix level adaptation weight coefficient is identical, thereby adaptive direction figure is identical with desired static directional diagram, thereby does not have the secondary lobe loss.And under situation about exist disturbing, this method forms device with the submatrix level optimal beam of main lobe conformal and compares, also effective suppressed sidelobes.
3. the 2nd level adaptation
If
a ( u , v ) = a ( u ) ⊗ a ( v ) - - - ( 19 )
After then the 1st level adaptation was handled, the directivity function of 4 wave beams was:
F sub _ Σ ( u , v ) = ( w ^ sub _ Σ ( mod i ) ) H [ ( T Σ ) H a ( u , v ) ] = f sub _ Σ A ( u ) f sub _ Σ E ( v ) F sub _ Δ E ( u , v ) = ( w ^ sub _ Δ A ( mod i ) ) H [ ( T Δ E ) H a ( u , v ) ] = f sub _ Δ A ( u ) f sub _ Σ E ( v ) F sub _ Δ AΔ ( u , v ) = ( w ^ sub _ Δ A ( mod i ) ) H [ ( T Δ A ) H a ( u , v ) ] = f sub _ Δ A ( u ) f sub _ Σ E ( v ) F sub _ Δ Δ ( u , v ) = ( w ^ sub _ Δ Δ ( mod i ) ) H [ ( T Δ Δ ) H a ( u , v ) ] = f sub _ Δ A ( u ) f sub _ Δ E ( v ) - - - ( 20 )
In the formula With
Figure BDA0000041153980000096
Be respectively the orientation to wave beam, pitching to wave beam, orientation to difference beam and pitching adaptive direction figure to difference beam.
Consideration suppresses the situation that main lobe disturbs along azimuth direction.If y The sub_ ∑,
Figure BDA0000041153980000097
With
Figure BDA0000041153980000098
Be respectively the 1st level adaptation and handle the output of back 4 passages, then:
y sub _ Σ = ( w ^ sub _ Σ ( mod i ) ) H x sub _ Σ ( t ) y sub _ Δ E = ( w ^ sub _ Δ E ( mod i ) ) H x sub _ Δ E ( t ) y sub _ Δ A = ( w ^ sub _ Δ A ( mod i ) ) H x sub _ Δ A ( t ) y sub _ Δ Δ = ( w ^ sub _ Δ Δ ( mod i ) ) H x sub _ Δ Δ ( t ) - - - ( 21 )
After suppressing main lobe and disturb along azimuth direction, pitch orientation and, the difference passage is output as:
y ^ sub _ E Σ = y sun _ Σ - w adap _ E Σ y sub _ Δ A y ^ sub _ E Δ = y sub _ Δ E - w adap _ E Δ y sub _ Δ Δ - - - ( 22 )
In the formula, With
Figure BDA0000041153980000104
For being used to suppress the self adaptation power that the azimuth direction main lobe disturbs, be used to make
Figure BDA0000041153980000105
With
Figure BDA0000041153980000106
The power output minimum; And
Figure BDA0000041153980000107
With
Figure BDA0000041153980000108
Be the auxiliary signal in the self-adaptive processing.
If:
R ΣΔ A = E [ y sub _ Σ · ( y sub _ Δ A ) * ] P MJ F sub _ Σ ( u MJ , v MJ ) F sub _ Δ A * ( u MJ , v MJ ) R Δ A Δ A = E [ y sub _ Δ A · ( y sub _ Δ A ) * ] = P MJ F sub _ Δ A ( u MJ , v MJ ) F sub _ Δ A * ( u MJ , v MJ ) + P N R Δ E Δ Δ = E [ y sub _ Δ E · ( y sub _ Δ Δ ) * ] = P MJ F sub _ Δ E ( u MJ , v MJ ) F sub _ Δ Δ * ( u MJ , v MJ ) R Δ Δ Δ Δ = E [ y sub _ Δ Δ · ( y sub _ Δ Δ ) * ] = P MJ F sub _ Δ Δ ( u MJ , v MJ ) F sub _ Δ Δ * ( u MJ , v Mj ) + P N - - - ( 23 )
* represents conjugation in the formula, P NBe noise power, then:
w adap _ E Σ = R Σ Δ A / R Δ A Δ A w adap _ E Δ = R Δ E Δ Δ / R Δ Δ Δ Δ - - - ( 24 )
With formula (23) substitution formula (24), then make an uproar than down working energetically, With
Figure BDA00000411539800001012
Approximate identical.Establish for this reason:
w adap _ E = ( w adap _ E Σ + w adap _ E Δ ) / 2 - - - ( 25 )
And order:
w adap _ E Σ = w adap _ E w adap _ E Δ = w adap _ E - - - ( 26 )
After then suppressing main lobe and disturbing, pitching face and, the adaptive direction function of difference beam is:
F sub _ E Σ ( MJS ) ( u , v ) = F sub _ Σ ( u , v ) - w adap _ E F sub _ Δ A ( u , v ) F sub _ E Δ ( MJS ) ( u , v ) = F sub _ Δ E ( u , v ) - w adap _ E F sub _ Δ Δ ( u , v ) - - - ( 27 )
Thereby be along the self adaptation pulse ratio of pitch orientation:
K E ( u , v ) = F sub _ E Δ ( NJS ) ( u , v ) F sub _ E Σ ( MJS ) ( u , v ) = F sub _ Δ E ( u , v ) - w adap _ E F sub _ Δ Δ ( u , v ) F sub _ Σ ( u , v ) - w adap _ E F sub _ Δ A ( u , v )
= f sub _ Δ E ( v ) [ f sub _ Σ A ( u ) - w adap _ E f sub _ Δ A ( u ) ] f sub _ Σ E ( v ) [ f sub _ Σ A ( u ) - w adap _ E f sub _ Δ A ( u ) ]
= f sub _ Δ E ( v ) f sub _ Σ E ( v ) - - - ( 28 )
As seen, the pulse of pitch orientation is than the both direction ratio of functions that only depends on this direction
Figure BDA0000041153980000114
On the other hand, after the interference of azimuth direction inhibition main lobe, can be not right
Figure BDA0000041153980000115
Exert an influence.Thereby by formula (28) as can be known, after the interference of azimuth direction inhibition main lobe, the self adaptation pulse ratio of pitch orientation remains unchanged.
Consider to suppress the situation that the pitch orientation main lobe disturbs below.Similar with above-mentioned situation, the self adaptation pulse ratio that can obtain azimuth direction is:
K A ( u , v ) = f sub _ Δ A ( u ) f sub _ Σ A ( u ) - - - ( 29 )
As seen, after the interference of pitch orientation inhibition main lobe, the self adaptation pulse of azimuth direction is than the both direction ratio of functions that only depends on this direction; Thereby its pulse ratio is preceding identical with the interference of inhibition main lobe, thereby kept original pulse performance.
The two-stage disturbance restraining method that contrasts existing self adaptation pulse illustrates principle of the present invention and advantage: the pulse performance of existing method is undesirable: compare with static pulse ratio, its self adaptation pulse is than bigger distortion being arranged departing from the beam position place.Its reason is in existing methods the 1st level adaptation that its submatrix level ADBF has adopted submatrix level optimal beam to form device, combines with the main lobe conformal simultaneously.The main lobe conformal effect of this method is undesirable, the pulse poor-performing; Another limitation is, because submatrix level ADBF has adopted submatrix level optimal beam to form device, thereby makes the sidelobe level of adaptive direction figure compare increase greatly with static directional diagram.
The present invention has introduced the directional diagram control method in the 1st level adaptation; Its submatrix level ADBF adopts submatrix level optimal beam to form device, but adopts the submatrix level adaptation Beam-former with directional diagram control performance.By this improvement, compare with existing method, not only suppressed self adaptation and, the secondary lobe of difference directional diagram, the more important thing is main lobe conformal effect is significantly improved, thereby effectively improved the pulse performance.
For the planar rectangular battle array that comprises 56 * 42 array elements, to establish array and be divided into 6 * 6 submatrixs, each submatrix is the rectangle battle array; Beam position is (0 a °, 0 °); Have a main lobe to disturb and a secondary lobe interference, wherein the main lobe interference radiating way is (1 a °, 2 °), and the secondary lobe interference radiating way is (10 °, 15 °), dried making an uproar than being 30dB.Adopt present technique, when the angle of pitch was-2.7 °, resulting self adaptation pulse had only 3.49% than with the error of static pulse ratio, and adopts existing two-stage disturbance restraining method, and error is up to 69.70%.Error of the present invention has fallen 66.21% on a year-on-year basis.

Claims (5)

1. the two-stage disturbance restraining method of submatrix level adaptation pulse is characterized in that: it comprises that two-stage is disturbed and suppresses:
The first order disturbs inhibition to suppress for the secondary lobe with main lobe conformal disturbs, and its method is: combine with the main lobe conformal by the submatrix level ADBF that will have the directional diagram control performance, when keeping the main lobe shape, interference suppresses to secondary lobe;
The second level is disturbed and is suppressed to suppress for main lobe disturbs, and its method is: adopt 4 passage monopulse systems to carry out the main lobe interference and suppress; The implication of 4 passages of described 4 passage monopulse systems is: each passage forms 1 wave beam in 4 receive paths, totally 4 wave beams, that is: and wave beam, trim wave beam, gun parallax wave beam and two difference beam.
2. the two-stage disturbance restraining method of submatrix level adaptation pulse according to claim 1, the concrete grammar that it is characterized in that first order disturbance restraining method is: introduce 1 each element and be 1 static dominant vector, construct the mismatch guiding vector simultaneously, and the guiding vector that the submatrix level optimal beam of main lobe conformal forms device is revised with the mismatch guiding vector, promptly be the product of 1 static dominant vector with the covariance matrix of the interference plus noise of submatrix level output under the noiseless situation and each element, replace the guiding vector of submatrix level phase array at the beam position place, thereby the submatrix level adaptation of the main lobe conformal that obtains revising power is finished the first order and is disturbed inhibition.
3. the two-stage disturbance restraining method of submatrix level adaptation pulse according to claim 1, it is characterized in that, in the disturbance restraining method of the second level, if the orientation to pitching to directional diagram separate, in 4 passage monopulse systems, suppress to disturb and form self adaptation zero along direction and fall into, along another orthogonal direction keep non-self-adapting and, difference beam, thereby obtain undistorted pulse ratio, promptly with static pulse than identical.
4. the two-stage disturbance restraining method of submatrix level adaptation pulse according to claim 3, it is characterized in that, in 4 passage monopulse systems, suppress interference and form self adaptation zero to fall into along pitch orientation, and along azimuth direction keep non-self-adapting and, difference beam, this moment, described difference beam was meant the gun parallax wave beam.
5. the two-stage disturbance restraining method of submatrix level adaptation pulse according to claim 3, it is characterized in that, in 4 passage monopulse systems, suppress interference and form self adaptation zero to fall into along azimuth direction, and along pitch orientation keep non-self-adapting and, difference beam, this moment, described difference beam was meant the trim wave beam.
CN2010106089947A 2010-12-28 2010-12-28 Two-level interference suppression method for sub-array level adaptive single pulse Pending CN102064892A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010106089947A CN102064892A (en) 2010-12-28 2010-12-28 Two-level interference suppression method for sub-array level adaptive single pulse

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010106089947A CN102064892A (en) 2010-12-28 2010-12-28 Two-level interference suppression method for sub-array level adaptive single pulse

Publications (1)

Publication Number Publication Date
CN102064892A true CN102064892A (en) 2011-05-18

Family

ID=44000013

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010106089947A Pending CN102064892A (en) 2010-12-28 2010-12-28 Two-level interference suppression method for sub-array level adaptive single pulse

Country Status (1)

Country Link
CN (1) CN102064892A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103178881A (en) * 2011-12-23 2013-06-26 中兴通讯股份有限公司 Main lobe interference suppression method and device
CN104076336A (en) * 2014-06-25 2014-10-01 西安电子科技大学 Sum/difference beam formation method based on radar folded area array rule digital subarray
CN104702379A (en) * 2013-12-10 2015-06-10 北京信威通信技术股份有限公司 Dynamic zero-trap method in wireless communication system
CN104777468A (en) * 2015-05-12 2015-07-15 上海无线电设备研究所 Method of using sum and difference channels of radar antenna for restraining radar side lobe interception
CN107888237A (en) * 2016-09-30 2018-04-06 北京三星通信技术研究有限公司 It is initially accessed the method, base station equipment and user equipment with Stochastic accessing
CN108141268A (en) * 2015-10-12 2018-06-08 Sk电信有限公司 The wireless communications method and device shaped using mixed-beam
CN113156380A (en) * 2021-03-19 2021-07-23 哈尔滨工业大学 Main lobe noise suppression interference self-adaptive suppression method based on transceiving mismatch processing
CN113406620A (en) * 2021-04-28 2021-09-17 北京理工大学 Distributed array angle measurement method for array decomposition

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101291165A (en) * 2007-04-17 2008-10-22 大唐移动通信设备有限公司 Sequence detecting method and apparatus for multi-antenna system
CN101702017A (en) * 2009-11-30 2010-05-05 中国人民解放军空军雷达学院 Multi-input multi-output radar waveform design and processing method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101291165A (en) * 2007-04-17 2008-10-22 大唐移动通信设备有限公司 Sequence detecting method and apparatus for multi-antenna system
CN101702017A (en) * 2009-11-30 2010-05-05 中国人民解放军空军雷达学院 Multi-input multi-output radar waveform design and processing method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《电子学报》 20090930 胡航等 一种改进的两级子阵级自适应单脉冲方法 第2、3节,图1 3-5 , *
胡航等: "一种改进的两级子阵级自适应单脉冲方法", 《电子学报》, 30 September 2009 (2009-09-30) *
胡航等: "子阵级自适应单脉冲的四通道主瓣干扰抑制", 《电波科学学报》, 31 October 2009 (2009-10-31), pages 2 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103178881B (en) * 2011-12-23 2017-08-25 南京中兴新软件有限责任公司 Main lobe interference suppression method and device
CN103178881A (en) * 2011-12-23 2013-06-26 中兴通讯股份有限公司 Main lobe interference suppression method and device
CN104702379A (en) * 2013-12-10 2015-06-10 北京信威通信技术股份有限公司 Dynamic zero-trap method in wireless communication system
CN104702379B (en) * 2013-12-10 2018-03-09 北京信威通信技术股份有限公司 A kind of dynamic null method in wireless communication system
CN104076336A (en) * 2014-06-25 2014-10-01 西安电子科技大学 Sum/difference beam formation method based on radar folded area array rule digital subarray
CN104076336B (en) * 2014-06-25 2016-09-07 西安电子科技大学 Based on radar folding face battle array rule numeral submatrix and the forming method of/difference beam
CN104777468A (en) * 2015-05-12 2015-07-15 上海无线电设备研究所 Method of using sum and difference channels of radar antenna for restraining radar side lobe interception
CN104777468B (en) * 2015-05-12 2017-08-25 上海无线电设备研究所 Suppress the method that radar secondary lobe is intercepted and captured using radar antenna and poor passage
CN108141268A (en) * 2015-10-12 2018-06-08 Sk电信有限公司 The wireless communications method and device shaped using mixed-beam
CN108141268B (en) * 2015-10-12 2021-07-06 Sk电信有限公司 Wireless communication method and apparatus using hybrid beamforming
CN107888237A (en) * 2016-09-30 2018-04-06 北京三星通信技术研究有限公司 It is initially accessed the method, base station equipment and user equipment with Stochastic accessing
CN107888237B (en) * 2016-09-30 2022-06-21 北京三星通信技术研究有限公司 Initial access and random access method, base station equipment and user equipment
CN113156380A (en) * 2021-03-19 2021-07-23 哈尔滨工业大学 Main lobe noise suppression interference self-adaptive suppression method based on transceiving mismatch processing
CN113156380B (en) * 2021-03-19 2022-05-17 哈尔滨工业大学 Main lobe noise suppression interference self-adaptive suppression method based on transceiving mismatch processing
CN113406620A (en) * 2021-04-28 2021-09-17 北京理工大学 Distributed array angle measurement method for array decomposition

Similar Documents

Publication Publication Date Title
CN102064892A (en) Two-level interference suppression method for sub-array level adaptive single pulse
CN101369014B (en) Bilateral constraint self-adapting beam forming method used for MIMO radar
CN105137399B (en) The radar self-adaption Beamforming Method filtered based on oblique projection
US9279884B2 (en) Method and device for estimating direction of arrival
CN103399303B (en) Airborne radar resists intensive deceiving jamming method and system
CN108508423A (en) Submatrix number based on special-shaped battle array and poor Monopulse estimation method
CN109946664B (en) Array radar seeker monopulse angle measurement method under main lobe interference
CN102195701A (en) Method for suppressing side lobes of sum beams and difference beams of planar phased array only by utilizing one kind of analogue weighting
CN102608580B (en) Ultra-low side lobe adaptive digital beam forming (ADBF) method for digital array
CN108051782A (en) Extensive phased array difference beam based on Subarray partition forms system
CN102142609A (en) Sub-array-class adaptive digital beam forming device with low side-lobe characteristics
CN103116162A (en) High-resolution sonar location method based on sparsity of objective space
CN107703478A (en) Extension aperture arrival direction estimation method based on cross-correlation matrix
CN105372633A (en) Phased array radar dimension reduction four-channel mainlobe sidelobe interference-resisting method
CN104931937B (en) Based on the normalized Subarray rectangular projection Beamforming Method of covariance matrix
CN110646765B (en) Riemann distance-based generalized sidelobe cancellation algorithm
CN110208757B (en) Steady self-adaptive beam forming method and device for inhibiting main lobe interference
CN103323810B (en) L-array azimuthal angle and pitch angle paired signal processing method
CN107064904A (en) A kind of carrier-borne high-frequency ground wave radar orientation High Resolution Method based on virtual aperture
CN100585429C (en) Passive channel correcting method based on non-linear antenna array
CN112858994B (en) Amplitude comparison direction finding method based on uniform circular array
CN103248412B (en) A kind of method of satellite navigation Multibeam synthesis
CN111551908B (en) Method for reducing complexity of phased array system array element activation algorithm
CN104868946A (en) Adaptive weighted interference suppression method of subarray level mixed MIMO-phased array system
CN103235959B (en) The method that aerial array output forms digital beam is made in read write line

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20110518