CN104931937A - Submatrix-level orthogonal projection (OP) wave beam forming method based on covariance matrix normalization - Google Patents
Submatrix-level orthogonal projection (OP) wave beam forming method based on covariance matrix normalization Download PDFInfo
- Publication number
- CN104931937A CN104931937A CN201510209368.3A CN201510209368A CN104931937A CN 104931937 A CN104931937 A CN 104931937A CN 201510209368 A CN201510209368 A CN 201510209368A CN 104931937 A CN104931937 A CN 104931937A
- Authority
- CN
- China
- Prior art keywords
- interference
- submatrix
- covariance matrix
- subarray
- self
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a submatrix-level OP wave beam forming method based on covariance matrix normalization. The method can be used to effectively inhibit interference, make a main lobe of a self-adaptive directional diagram conformal, reduce side lobes the self-adaptive directional diagram and obtain higher output SINR and convergence speed. The method comprises that submatrix-level reception signals are normalized, and a corresponding normalized sampling covariance matrix is calculated; the amount of interference signal sources is estimated by utilizing the MDL criterion, and further an interference subspace is obtained; and finally, the static weight vector is projected to an orthogonal complement space of the interference subspace to obtain a self-adaptive weight vector.
Description
Technical field
The present invention relates to array signal process technique field, be specifically related to a kind of based on covariance matrix normalized Subarray rectangular projection Beamforming Method.
Background technology
Array Signal Processing is an important branch in signal transacting field, and it is widely used in radar, sonar, communication, navigation, seismic monitoring, Speech processing and biomedical engineering etc.Adaptive beamformer is an important research content in Array Signal Processing, and its essence is exactly by adaptive weighted to each array element, carries out airspace filter, thus reaches and strengthen wanted signal, suppress undesired signal and weaken the object of noise signal.Response that minimum variance is undistorted (MVDR) is a kind of relatively conventional algorithm, and it is by being 1 at wanted signal direction constrain array gain, and makes array output power minimum, thus reaches the object suppressing to disturb.Covariance matrix is inverted (SMI), and algorithm is a kind of conventional method realizing MVDR algorithm, but when lower snap, the output SINR (Signal to Interference plus Noise Ratio) of this algorithm and the speed of convergence of adaptive direction figure slower.
In actual applications, consider hardware condition and environmental factor, the fast umber of beats of sampling calculating adaptive weight employing is less.In order in low snap situation, solve the problem that SMI algorithm brings, rectangular projection (OP) algorithm is widely used, and its essence is and is projected to by static weight vector in the orthogonal complement space (i.e. noise subspace) of interference space, and then obtain self-adaptation weight vector.In this algorithm, little eigenwert characteristic of correspondence vector does not participate in the calculating of self-adaptation weight vector, so under low snap condition, this algorithm can make output SINR and adaptive direction figure rapidly converge to optimal value.But when OP algorithm application is to Subarray, the uneven division of submatrix can cause each submatrix noise output power unequal, and then the accuracy of MDL criterion estimation can be affected, thus cause the interference space of estimation inaccurate, cause adaptive direction figure main lobe to be out of shape and secondary lobe rising, export SINR degradation.
Summary of the invention
In view of this, the invention provides a kind of based on covariance matrix normalized Subarray rectangular projection Beamforming Method, effectively can suppress interference, and the main lobe conformal of adaptive direction figure and secondary lobe can be made to reduce, and higher output SINR and speed of convergence faster can be obtained.
Of the present invention based on covariance matrix normalized Subarray rectangular projection Beamforming Method, comprise the steps:
Step 1, is normalized Subarray Received signal strength, calculates the interference-plus-noise covariance matrix R after normalization
sub_norm:
wherein, R
in_subfor the covariance matrix of Subarray; T
lfor normalization matrix:
wherein, L is submatrix number,
W
ibe the weighting coefficient of i-th array element, U=N
0+ N
1+ ... + N
l-1-J
0-J
1-...-J
l-1+ 1, Q=N
0+ N
1+ ... + N
l-J
0-J
1-...-J
l-1, N
ibe the array number that i-th (0≤i≤L-1) individual submatrix comprises, J
ibe i-th (0≤i≤L-2) individual submatrix with the overlapping array number of the i-th+1 submatrix; ()
hrepresent complex-conjugate transpose;
Step 2, utilizes MDL criterion to estimate interference space:
Step 2.1, to the interference-plus-noise covariance matrix R after normalization
sub_normcarry out Eigenvalues Decomposition, obtain eigenwert and characteristic of correspondence vector thereof, and eigenwert is carried out descending arrangement;
Step 2.2, utilizes MDL criterion to estimate the number of interference signal source for P, then in step 2.1, front P column vector of eigenvector forms interference space;
Step 3, utilizes the interference space that step 2 estimates, and adopts rectangular projection Adaptive beamformer method, solves self-adaptation weight vector;
Step 4, the self-adaptation weight vector utilizing step 3 to obtain, is weighted process to the echo data received, and obtains adaptive beam.
Beneficial effect:
The invention solves when the uneven and lower sampling snap of Subarray partition, the interference space that traditional Subarray orthogonal projection algorithm is estimated is inaccurate, interference is caused not to be effectively suppressed, and the main lobe distortion of adaptive direction figure, the problem that secondary lobe raises, can effectively complete Subarray Adaptive beamformer, fall in the adaptive formation zero of interference radiating way, and while effectively suppressing interference, make the main lobe conformal of adaptive direction figure and secondary lobe reduce, the present invention has higher output Signal to Interference plus Noise Ratio after Adaptive beamformer process, and output Signal to Interference plus Noise Ratio has speed of convergence faster.
Accompanying drawing explanation
Fig. 1 is process flow diagram of the present invention.
Fig. 2 is the inventive method and improves front method adaptive direction figure comparison diagram (when fast umber of beats is 2 times of submatrix numbers).
Fig. 3 is the inventive method and improves front method adaptive direction figure comparison diagram (when fast umber of beats is 10 times of submatrix numbers).
Fig. 4 is that the inventive method exports SINR with fast umber of beats change curve comparison diagram with the front method of improvement.
Fig. 5 is that the inventive method exports SINR with beam position change curve comparison diagram with the front method of improvement.
Embodiment
To develop simultaneously embodiment below in conjunction with accompanying drawing, describe the present invention.
The invention provides a kind of based on covariance matrix normalized Subarray rectangular projection Beamforming Method, first Subarray Received signal strength is normalized, and calculates corresponding normalization sample covariance matrix; Then utilize MDL criterion to estimate interference signal source number, and then obtain interference space; Finally static weight vector projected to the orthogonal complement space of interference space and obtain self-adaptation weight vector.When the uneven and lower sampling snap of Subarray partition, the present invention can suppress interference effectively, and the main lobe conformal of adaptive direction figure and secondary lobe can be made to reduce, and can obtain higher output SINR and speed of convergence faster.As shown in Figure 1, concrete steps are as follows for flow process of the present invention:
Step one, construct normalized Subarray covariance matrix
1. the foundation of signal model
Suppose an arrowband linear array, altogether N number of array element, array element is isotropy, P undesired signal, and undesired signal is far field narrow band signal, white noise when supposing that each array element noise is separate, power equal empty and undesired signal and noise are uncorrelated mutually.Then array received to signal model can be expressed as
X
in(t)=AS(t)+N(t) (1)
In formula, A=[a (θ
1), a (θ
2) ... a (θ
p)] be array manifold matrix, a (θ
i) (i=1,2 ... P) be the steering vector of undesired signal, if the spacing of the n-th array element and reference point is d
n(n=0,1,2 ..., N-1), usually with the 0th array element for reference point, now d
0=0, λ is wavelength, then
θ
i(i=1,2 ... P) be the incident angle of undesired signal, []
tfor matrix transpose, S (t)=[S
1(t), S
2(t) ... S
p(t)]
t, S
i(t) (i=1,2 ... P) be the complex envelope of i-th undesired signal, N (t)=[n
1(t), n
2(t) ..., n
n(t)] be background white noise.
Thus obtain array covariance matrix and be
R
in=E{X
in(t)X
in H(t)} (2)
In formula, E{} represents mathematical expectation, ()
hrepresent complex-conjugate transpose.
In practical application, according to maximal possibility estimation criterion, by limited fast beat of data X
in(t
i) estimate array covariance matrix,
In formula, X
in(t
i) for i (i=1,2 ..., K) and the sampled value of moment array, K is the fast umber of beats of sampling.
When Subarray process, be L submatrix by array partition, and (L >=P+1), can be non-overlapped submatrix or overlapping submatrix, submatrix transition matrix can be expressed as
T=φ
0WT
0(4)
Wherein
represent the effect of phase shifter, if beam position is identical with wanted signal direction; W=diag (w
n)
n=0,1 ..., N-1, wherein w
nbe the weighting coefficient of the n-th array element, for suppressing the sidelobe level of directional diagram; T
0for N × L submatrix formed matrix, its l (l=0,1 ..., L-1) in all elements that arranges, only have the element value corresponding with the array element sequence number of l submatrix to be 1, all the other be 0 (when non-overlapped submatrix, T
0column vector mutually orthogonal).
The sampling snap signal then Subarray received is
X
in_sub(t)=T
HX
in(t) (5)
Then the covariance matrix of Subarray is
2. the normalization of covariance matrix
First export each submatrix and be normalized, matrix T is passed through in normalization
lcomplete
Wherein
U=N
0+N
1+…+N
l-1-J
0-J
1-…-J
l-1+1
Q=N
0+N
1+…+N
l-J
0-J
1-…-J
l-1
N
ibe the array number that i-th (0≤i≤L-1) individual submatrix comprises, J
ibe i-th (0≤i≤L-2) individual submatrix with the overlapping array number of the i-th+1 submatrix.
Interference-plus-noise covariance matrix after normalization is
By normalized, make the noise power of each submatrix consistent, thus MDL criterion can be suitable for.
Step 2, estimation interference space
To normalized covariance matrix R
sub_normcarry out Eigenvalues Decomposition
In formula, λ
i(i=1,2 ..., L) and be covariance matrix R
sub_normeigenwert,
for with eigenvalue λ
icharacteristic of correspondence vector, λ
idescending arrangement
Utilize MDL criterion to estimate the number of interference signal source, and then estimate interference space.
The function of MDL criterion is
Wherein
From MDL criterion, when the numerical value change of d, the value of the d corresponding when formula (11) gets minimum value is the number P of interference signal source, selected characteristic vector
front P column vector composition interference space U
s, namely
by the known vector of mathematical knowledge
with vector a (θ
1), a (θ
2) ..., a (θ
p) open into same vector space, that is:
Wherein, span{} represents the space that vector is opened,
be the estimation of interference space.
Step 3, solve the self-adaptation weight vector of innovatory algorithm
Utilize the interference space that step 2 estimates, adopt rectangular projection Adaptive beamformer method, solve the self-adaptation weight vector of innovatory algorithm.
Conventional orthogonal is adopted to throw the thought of algorithm, by static weight vector w
q_subthe interference space U estimated in step 2
sthe orthogonal complement space projection, the self-adaptation weight vector of the algorithm that is improved is
In formula, I is that L × L ties up unit matrix, and η is a constant, w
q_subfor static weight vector, and each element is the L dimensional vector of 1,
effect be that the gain in antenna main beam direction is remained unchanged.
Step 4, carries out Adaptive beamformer to the echo received
After obtaining self-adaptation weight vector, process can be weighted to the echo data received:
Y=W
HX(t) (13)
In formula, X (t) is the echoed signal received, thus effectively removes interference and weaken noise, and strengthens wanted signal.
Since then, just complete a kind of based on the sub-Adaptive beamformer method of covariance matrix normalized Subarray rectangular projection to the process of echo data.
In order to verify that one that the present invention proposes is based on covariance matrix normalized Subarray rectangular projection Adaptive beamformer method, carry out the emulation of adaptive beam directional diagram and output Signal to Interference plus Noise Ratio (SINR), emulation adopts arrowband even linear array, and simulation parameter is as shown in table 1.The front algorithm of improvement is that the sampling snap signal of Subarray directly adopts rectangular projection (OP) algorithm to calculate self-adaptation weight vector.
Table 1 simulation parameter is arranged
Fig. 2 and Fig. 3 be innovatory algorithm of the present invention with improve before the comparing (emulating 1 time) of adaptive beam directional diagram of algorithm, fast umber of beats of sampling is respectively 20 and 100, beam position angle is 0 °, can find out, the adaptive beam major lobe of directional diagram that before improving, algorithm obtains is out of shape and sidelobe level seriously raises; The adaptive beam major lobe of directional diagram conformal that obtains of algorithm after improving and sidelobe level is lower, close to static beam pattern, performance improves compared with before-improvement widely.
Fig. 4 is that wanted signal angle is 0 °, and input signal-to-noise ratio is 0dB, the same Fig. 2 of other simulated conditions under difference sampling snap condition, the comparison of the output Signal to Interference plus Noise Ratio (SINR) of innovatory algorithm of the present invention and the front algorithm of improvement.From simulation result, after improving, the output SINR of algorithm is higher, and convergence is very fast; And the output SINR of algorithm is lower before improving, convergence is comparatively slow, and along with the increase of fast umber of beats, export SINR to decline to some extent, because fast umber of beats of sampling is higher, the degree of accuracy of the interference space of estimation reduces, interference can not be effectively suppressed, and causes the SINR meeting degradation exported.
Fig. 5 is when different beams orientation angle, the comparison of the output Signal to Interference plus Noise Ratio (SINR) of innovatory algorithm of the present invention and the front algorithm of improvement, input signal-to-noise ratio is 0dB, the same Fig. 2 of other simulated conditions, can find out that after improving, algorithm can suppress interference effectively, and the SINR exported is higher.
Can obtain from Fig. 2 ~ Fig. 5, innovatory algorithm of the present invention can strengthen wanted signal, has good interference free performance, is a kind of sane Subarray adaptive beam-forming algorithm.
In sum, these are only preferred embodiment of the present invention, be not intended to limit protection scope of the present invention.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (1)
1., based on a covariance matrix normalized Subarray rectangular projection Beamforming Method, it is characterized in that, comprise the steps:
Step 1, is normalized Subarray Received signal strength, calculates the interference-plus-noise covariance matrix R after normalization
sub_norm:
wherein, R
in_subfor the covariance matrix of Subarray; T
lfor normalization matrix:
wherein, L is submatrix number,
W
ibe the weighting coefficient of i-th array element, U=N
0+ N
1+ ... + N
l-1-J
0-J
1-...-J
l-1+ 1, Q=N
0+ N
1+ ... + N
l-J
0-J
1-...-J
l-1, N
ibe the array number that i-th (0≤i≤L-1) individual submatrix comprises, J
ibe i-th (0≤i≤L-2) individual submatrix with the overlapping array number of the i-th+1 submatrix; ()
hrepresent complex-conjugate transpose;
Step 2, utilizes MDL criterion to estimate interference space:
Step 2.1, to the interference-plus-noise covariance matrix R after normalization
sub_normcarry out Eigenvalues Decomposition, obtain eigenwert and characteristic of correspondence vector thereof, and eigenwert is carried out descending arrangement;
Step 2.2, utilizes MDL criterion to estimate the number of interference signal source for P, then in step 2.1, front P column vector of eigenvector forms interference space;
Step 3, utilizes the interference space that step 2 estimates, and adopts rectangular projection Adaptive beamformer method, solves self-adaptation weight vector;
Step 4, the self-adaptation weight vector utilizing step 3 to obtain, is weighted process to the echo data received, and obtains adaptive beam.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510209368.3A CN104931937B (en) | 2015-04-28 | 2015-04-28 | Based on the normalized Subarray rectangular projection Beamforming Method of covariance matrix |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510209368.3A CN104931937B (en) | 2015-04-28 | 2015-04-28 | Based on the normalized Subarray rectangular projection Beamforming Method of covariance matrix |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104931937A true CN104931937A (en) | 2015-09-23 |
CN104931937B CN104931937B (en) | 2017-09-29 |
Family
ID=54119190
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510209368.3A Active CN104931937B (en) | 2015-04-28 | 2015-04-28 | Based on the normalized Subarray rectangular projection Beamforming Method of covariance matrix |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104931937B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054142A (en) * | 2016-05-13 | 2016-10-26 | 西安电子科技大学 | Airborne multi-input-multi-output radar main lobe smart interference inhibition method and system |
CN106199547A (en) * | 2016-06-30 | 2016-12-07 | 西安电子科技大学 | Weak target detection method at a slow speed based on external illuminators-based radar |
CN107064884A (en) * | 2017-01-05 | 2017-08-18 | 西安电子科技大学 | The Adaptive beamformer method of rule-based overlapping submatrix |
CN114609651A (en) * | 2022-03-28 | 2022-06-10 | 电子科技大学 | Space domain anti-interference method of satellite navigation receiver based on small sample data |
CN116112323A (en) * | 2021-11-10 | 2023-05-12 | 大唐移动通信设备有限公司 | Interference suppression method, device, equipment and storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004257761A (en) * | 2003-02-24 | 2004-09-16 | Toshiba Corp | Radar signal processing device and method |
CN102830387A (en) * | 2012-08-23 | 2012-12-19 | 北京理工大学 | Data preprocessing based covariance matrix orthogonalization wave-beam forming method |
CN103885045A (en) * | 2014-04-09 | 2014-06-25 | 西安电子科技大学 | Sub-array division based circulation combined adaptive beam forming method |
-
2015
- 2015-04-28 CN CN201510209368.3A patent/CN104931937B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004257761A (en) * | 2003-02-24 | 2004-09-16 | Toshiba Corp | Radar signal processing device and method |
CN102830387A (en) * | 2012-08-23 | 2012-12-19 | 北京理工大学 | Data preprocessing based covariance matrix orthogonalization wave-beam forming method |
CN103885045A (en) * | 2014-04-09 | 2014-06-25 | 西安电子科技大学 | Sub-array division based circulation combined adaptive beam forming method |
Non-Patent Citations (2)
Title |
---|
杨志伟 等: ""子空间重构的一类自适应波束形成算法"", 《电子与信息学报》 * |
邓新红: ""子阵级自适应数字波束形成方法研究"", 《中国优秀博硕士学位论文全文数据库 (硕士) 信息科技辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054142A (en) * | 2016-05-13 | 2016-10-26 | 西安电子科技大学 | Airborne multi-input-multi-output radar main lobe smart interference inhibition method and system |
CN106054142B (en) * | 2016-05-13 | 2018-08-03 | 西安电子科技大学 | A kind of airborne MIMO radar main lobe smart munition suppressing method and system |
CN106199547A (en) * | 2016-06-30 | 2016-12-07 | 西安电子科技大学 | Weak target detection method at a slow speed based on external illuminators-based radar |
CN107064884A (en) * | 2017-01-05 | 2017-08-18 | 西安电子科技大学 | The Adaptive beamformer method of rule-based overlapping submatrix |
CN116112323A (en) * | 2021-11-10 | 2023-05-12 | 大唐移动通信设备有限公司 | Interference suppression method, device, equipment and storage medium |
CN116112323B (en) * | 2021-11-10 | 2024-06-07 | 大唐移动通信设备有限公司 | Interference suppression method, device, equipment and storage medium |
CN114609651A (en) * | 2022-03-28 | 2022-06-10 | 电子科技大学 | Space domain anti-interference method of satellite navigation receiver based on small sample data |
Also Published As
Publication number | Publication date |
---|---|
CN104931937B (en) | 2017-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103984676A (en) | Rectangular projection adaptive beamforming method based on covariance matrix reconstruction | |
CN106569181A (en) | Algorithm for reconstructing robust Capon beamforming based on covariance matrix | |
CN102830387B (en) | Data preprocessing based covariance matrix orthogonalization wave-beam forming method | |
CN101369014B (en) | Bilateral constraint self-adapting beam forming method used for MIMO radar | |
CN104270179A (en) | Self-adaptive beam forming method based on covariance reconstruction and guide vector compensation | |
CN103837861B (en) | The Subarray linear restriction Adaptive beamformer method of feature based subspace | |
CN103245956B (en) | A kind of GPS anti-multipath method based on robust ada-ptive beamformer algorithm | |
CN104931937A (en) | Submatrix-level orthogonal projection (OP) wave beam forming method based on covariance matrix normalization | |
CN105302936A (en) | Self-adaptive beam-forming method based on related calculation and clutter covariance matrix reconstruction | |
CN104408278A (en) | A method for forming steady beam based on interfering noise covariance matrix estimation | |
CN103885045B (en) | Based on the circulation associating Adaptive beamformer method of Subarray partition | |
CN109407055A (en) | The Beamforming Method utilized based on multipath | |
CN107462872A (en) | A kind of anti-major lobe suppression algorithm | |
CN109254261A (en) | Coherent signal null based on uniform circular array EPUMA deepens method | |
CN105354171B (en) | A kind of projection subspace estimation adaptive beam synthetic method for improving characteristic vector | |
CN106324625A (en) | Adaptive anti-interference method for satellite navigation system based on 2-norm multi-target optimization | |
CN109143190A (en) | A kind of broadband robust adaptive beamforming method of null broadening | |
CN105301580A (en) | Passive detection method based on split array cross-spectrum phase difference variance weighing | |
CN106019252A (en) | Sum-difference tracking angle measurement method based on Nested array | |
CN105158741A (en) | Matrix reconstruction based adaptive anti-interference multipath multi-beam forming method | |
CN109298395A (en) | A kind of thinned array Beamforming Method based on maximum Signal to Interference plus Noise Ratio | |
CN103293517A (en) | Diagonal-loading robust adaptive radar beam forming method based on ridge parameter estimation | |
CN104459635B (en) | Self adaptation air filter filtering method based on iterative shrinkage Weighted Fusion | |
CN104678362B (en) | MIMO sky-wave OTH radar waveform optimization method | |
Shirvani-Moghaddam et al. | A comprehensive performance study of narrowband DOA estimation algorithms |
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 |