CN104515969A - Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method - Google Patents
Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method Download PDFInfo
- Publication number
- CN104515969A CN104515969A CN201310454248.0A CN201310454248A CN104515969A CN 104515969 A CN104515969 A CN 104515969A CN 201310454248 A CN201310454248 A CN 201310454248A CN 104515969 A CN104515969 A CN 104515969A
- Authority
- CN
- China
- Prior art keywords
- array
- matrix
- signal
- hexagonal
- centerdot
- 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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention belongs to the array signal processing field and relates to a hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method, in particular to a high-precision, de-coherence, two-dimensional signal DOA subspace estimation method. According to the estimation method of the invention, an array is virtualized and extended through performing conjugate matrix solution on array receiving data, so that a hexagonal array including three sub arrays are obtained, and then, a correlation matrix of a subspace algorithm is extended and reconstructed according to an autocorrelation matrix and a cross correlation matrix of the three sub arrays, and therefore, de-coherence can be realized; and singular value decomposition is performed on the new correlation matrix, so that a signal subspace and a noise subspace can be obtained, and then, DOA estimation is performed on coherent receiving signals through utilizing a two-dimensional MUSIC (Multiple Signal Classification) estimation method. With the estimation method of the invention adopted, coherence between source signals can be eliminated, and the estimation accuracy of two-dimensional DOA can be further improved. The estimation method is mainly applied to two-dimensional high-precision DOA estimation of coherent signals.
Description
Technical field
The invention belongs to antenna array signals process field, it is in particular to the Subspace Estimation Method of a kind of high precision, decoherence, 2D signal arrival bearing (DOA).
Background technology
Super-resolution Mutual coupling (DOA) technology is an important research point of Estimation of Spatial Spectrum, is widely used in high-resolution array signal transacting.High precision and its application prospect of its angle estimation cause people to pay close attention to greatly.Multiple signal classification algorithm (MUSIC, multiple signal classification) and be two kinds of classical super-resolution subspace class algorithms based on the Signal parameter estimation algorithm (ESPRIT, estimation of signal parameter via rotational inviance techniques) of ESPRIT.Wherein, MUSIC algorithm can be applied to the array of arbitrary shape, and can estimate the source signal of multiple parameter.But in actual applications, when snap number is limited for noise, and therefore angle measurement accuracy and the resolving power of MUSIC algorithm are restricted, and array number can be less than by direction-finding signal number.In actual environment, the impact due to multipath effect makes the performance of classical super-resolution DOA algorithm for estimating decline rapidly, even complete failure.The Mutual coupling of coherent signal is the practical problems that Estimation of Spatial Spectrum needs solution badly.
At present, the method for coherent signal decorrelation LMS comprises: the frequency domain smoothing method in space smoothing class algorithm, matrix reconstruction algorithm and non-reducing dimension algorithm in reducing dimension algorithm and virtual array method of changing etc.There is the shortcoming of array aperture loss in reducing dimension algorithm, but not reducing dimension algorithm is often for specific environment.For the estimation problem in coherent signal source, propose the approach of many solutions.The MUSIC decorrelation LMS algorithm improved based on Search Space Smoothing has good performance, but need array to be divided into multiple submatrix, not only reduce array number and array aperture, also pair array structure has certain requirement, good decorrelation performance cannot be obtained in actual applications, for arbitrary shape array, as hexagonal array does not have the algorithm of good decoherence.
The present invention is under hexagonal array structural model, goes the arrival direction estimation problem solving coherent signal.Sexangle planar array antenna be by an antenna element being positioned at center and several concentric hexagonal rings battle arrays that are the center of circle with it (each containing element number of array for 6n, be the hexagonal rings battle array near the center of circle during n=1) composition, it also can be seen as the modified circle ring array that a center exists array element.It is narrower that this kind of array makes main lobe become, and secondary lobe is lower, is effectively suppressed the interference beyond main lobe.The present invention studies based on the simplest hexagonal array structure, and due to its array structure, the reducing dimension algorithm of decorrelation is no longer applicable.The present invention is according to the design feature of hexagonal battle array self, utilize the characteristic that actual array output signal is complex signal, the thought of conjugate matrices is utilized to carry out virtual extended array, and utilize based on merging the thought that extended matrix increases Correlation Moment rank of matrix in the algorithm of matrix decomposition, to realize decorrelation LMS object, the basis of not reducing array element number and array aperture proposes a kind of arrival direction estimation algorithm of the decoherence based on expansion reconstruct correlation matrix.
Summary of the invention
Technical matters to be solved by this invention is to provide a kind of under the prerequisite not reducing array element number and array aperture, utilize the arrival direction estimation method of expansion reconstruct correlation matrix decoherence, the method for non-uniform linear arrays, can be estimated the arrival bearing in coherent signal source accurately.
The thought that present invention utilizes conjugate matrices carrys out array extending, pair array receive data ask based on its conjugate matrices, expansion reconstruct correlation matrix, provide a kind of based on hexagonal array, decoherence, two dimension DOA estimate.The raw data matrix that snap receives and conjugate matrices thereof, being equivalent to available array length is original 2 times, takes full advantage of the auto-correlation between raw data matrix and Qi Xin matrix and cross-correlation information, improves the performance that DOA estimates.
Based on a coherent signal arrival direction estimation method for hexagonal array, comprise the steps:
Step one: sexangle antenna array elements number is N, exporting the fast umber of beats of data is K, then each array element exports the data matrix X complex matrix that (k)=AS (k)+N (k) is N × K, asks its conjugate matrices X
*(k), and the correlation matrix obtaining former array output data is R
xX=E{X (k) X (k)
h, H represents its conjugate transpose, and A is the steering vector matrix of array, and S (k) represents source signal vector matrix; N (k) represents that noise average that array exports be zero variance is σ
2additive white Gaussian noise, and uncorrelated with source signal;
Step 2: conjugation data matrix is for aerial array, the data being equivalent to the symmetrical virtual array of its conjugation export, for the array after conjugation virtual extended, three sub-hexagonal arrays can be divided into, first sub-hexagonal array is real hexagonal array, and second sub-hexagonal array is centered by reference array element, in submatrix array element wherein half belong to the part of true array, second half is the part of conjugation virtual array, and the 3rd sub-hexagonal array is classified as conjugate radical matroid row;
Step 3: according to the array element contained by every sub-hexagonal array of array after expansion, draw the output data matrix X of sub-hexagonal array
1(k), X
2(k), X
3(k);
Step 4: autocorrelation matrix and the cross-correlation matrix of obtaining three sub-hexagonal arrays respectively:
Step 5: construct the correlation matrix made new advances
Step 6: carry out svd to R, decompositing N number of eigenwert is λ
1>=λ
2>=...>=λ
p>=λ
p+1=...=λ
n=σ
2, by judging that the number of large eigenwert carrys out estimated signal source number, and obtain signal subspace Es and noise subspace matrix En respectively according to corresponding proper vector;
Step 7: utilize two-dimentional MUSIC algorithm to build spatial spectrum function,
θ and
be respectively the angle of pitch and the position angle of source signal,
be with θ and
the array steering vector to received signal of change.Make pitching angle theta at (0 °, 90 °) and position angle
in (0 °, 360 °) scope during change, find out spatial spectrum
angle corresponding to maximum point is the DOA of source signal.
Good effect of the present invention: the present invention receives data by pair array and asks conjugation to carry out virtual extended array, obtain the hexagonal battle array of three subarrays, then expand the correlation matrix of reconstruct Subspace algorithm according to the autocorrelation matrix of three subarrays and cross-correlation matrix, achieve the object of decoherence.Svd is carried out to new correlation matrix and obtains signal subspace and noise subspace, recycle two-dimentional MUSIC Power estimation method and DOA estimation is carried out to relevant incoming wave signal.The present invention further increases the estimated accuracy performance of two-dimentional DOA while coherence between source signal removing.
Accompanying drawing explanation
Fig. 1 is two kinds of array structures of hexagonal array: the sexangle planar array structural drawing of 7 array elements and 19 array elements;
Fig. 2 is the figure carrying out wave method of the simplest hexagonal array Received signal strength;
Fig. 3 is for the array of the minimum array number of hexagonal array, the array after the Conjugate extended drawn, and draws the structural drawing being divided into 3 subarrays;
Embodiment
Just further describe the course of work of this invention below by reference to the accompanying drawings.The arrival direction estimation process of decoherence can be divided into three parts.
Part I mainly receives data modeling, the virtual array extending of conjugation for hexagonal array, is divided into subarray form according to array structure, is mainly step one and step 2.
Suppose have narrow band signal source, D (D≤N) individual far field to incide in sexangle planar array as shown in Figure 2, element number of array is N, and array element distance is r, and the angle of pitch and the position angle of incoming signal are respectively
the noise that array exports is the additive white Gaussian noise of zero-mean, and variance is σ
2, and uncorrelated with source signal.
For being that example is to study the output data model of array containing the hexagonal battle array of minimum array number.
Get the leftmost array element of accompanying drawing 2 midplane array as true origin and reference array element, then the Received signal strength of this array is:
X(t)=AS(t)+N(t) (1)
Wherein A represents array manifold matrix; S (t)=[s
1(t), s
2(t) ... s
d(t)]
trepresent signal phasor matrix; N (t)=[n
1(t), n
2(t) ... n
n(t)]
trepresent the noise vector matrix of array, the transpose operation of T representing matrix.
Wherein:
c is the light velocity, c=3 × 10
8m/s, d=1,2 ... D;
In practical application, the output signal of array is complex signal, asks its conjugate matrices, carrys out array extending, as shown in Figure 3.For the data X (k) that array snap receives, then the reception data of its conjugation array are X
*(k), the output data of array can be expressed as:
X(k)=[x
1(k),…x
n(k),…x
N(k)]
T(3)
The correlation matrix of array is expressed as:
R
xX=E{X (k) X (k)
h(4) wherein, the conjugate transpose operation of H representing matrix.
The then data X of conjugate matrices
*(k) be:
X
*(k)=[x
1 *(k),…x
n *(k),…x
N *(k)]
T(5)
Namely as shown in Figure 2, be 3 crossing hexagonal arrays the array partition after virtual extended, array element spacing distance is respectively r to step 2.
Part II is the process for array after expansion, mainly comprises step 3 to step 5, comprises two parts content, is first for the array element contained by every sub-hexagonal array of array after expansion, draws the output data matrix X of 3 sub-hexagonal arrays
1(k), X
2(k), X
3k () is secondly autocorrelation matrix and the cross-correlation matrix of obtaining three sub-hexagonal arrays respectively, expansion reconstructs the correlation matrix made new advances, and reaches the object of decoherence.
(1) subarray exports the structure of data
The output data of each subarray can be expressed as:
Find out from the Array Model after expansion, submatrix 2 and submatrix 3 can be obtained to the distance of left r along x-axis successively by real hexagonal array, and according to formula (3) and (5), above formula can be expressed as again:
X
1(k)=AS(k)+N(k)
X
2(k)=A·Φ
0S(k)+N(k) (7)
X
3(k)=A·Φ
0 2S(k)+N(k)
In formula, Φ
0for diagonal matrix:
(2) structure of expansion reconstruct correlation matrix
According in (1), spatial autocorrelation matrix and the cross-correlation matrix of each subarray can be obtained:
In formula: E{ } represent statistical expection; H represents conjugate transpose operation; R
s=E{S (k) S
h(k) } represent signal covariance matrix;
R
n=σ
2i represents the covariance matrix of array noise, and I is unit matrix.
According to above-mentioned autocorrelation matrix and cross-correlation matrix, merging the correlation matrix be expanded is:
Can be expressed as again according to formula (8) and (9):
R=A[R
SΦ
0R
SΦ
0R
SΦ
0 H(Φ
0 2)R
S(Φ
0 2)R
SΦ
0 H(Φ
0 2)R
S(Φ
0 2)
H]A
H(10)
+[R
NR
NR
NR
NR
NR
N]
Order
G=[R
SΦ
0R
SΦ
0R
SΦ
0 H(Φ
0 2)R
S(Φ
0 2)R
SΦ
0 H(Φ
0 2)R
S(Φ
0 2)
H] ,
R
n=[R
nr
nr
nr
nr
nr
n], then:
R=AGA
H+R
n(11)
Theorem 1: establish Q to be the matrix (L≤M) of L × M dimension without zero row vector, P is that L × L ties up diagonal matrix, and its diagonal element is unequal mutually, if rank (Q)=r<L, then rank [Q PQ]=r+1.
Theorem 2: establish H to be the matrix (L≤M) of M × L dimension without zero row vector, P is that L × L ties up diagonal matrix, and its diagonal element is unequal mutually, if rank (H)=r<L, then rank [H HP]=r+1.
Suppose that D (D≤N) the individual signal reaching array is all relevant, then rank (R
s)=1, because each signal angle is not identical, when the angle of pitch sum at the position angle He another information source that there is not Arbitrary Information Sources be 90 ° (namely
) equal with any two angles of pitch time their position angle sum be 360 ° of (θ
i=θ
j,
) and the position angle of information source be not 90 ° and 270 ° time, Φ
0diagonal element be different, be non-singular matrix.Theory deduction from above: rank (R
s)+n≤rank (G)≤D, wherein n is the number of matrix in G.Then, can decomposite the signal subspace of D dimension according to formula (11) known R, recycling DOA algorithm for estimating, as MUSIC algorithm, obtains the angle information of D signal, namely achieves the object of decorrelation.The method is according to relevant rank of matrix new after its expansion reconstruct, and the maximum coherent signal number that can estimate is at least N-1, adds and can estimate the Number of Coherent Sources.
Part III, mainly on the basis of expansion reconstruct correlation matrix, utilizes Subspace algorithm to carry out the process of DOA estimation, comprises step 6 and step 7.Theoretical analysis according to Part II can be known, for coherent source signal, the new relevant rank of matrix increase of reconstruct can be removed relevant, and on this basis, step 6 make use of svd and decomposes matrix R, is σ because noise is zero-mean variance
2additive white Gaussian noise, matrix G is full rank, knownly can decomposite N number of eigenwert and is:
λ
1≥λ
2≥…≥λ
P≥λ
P+1=…=λ
N=σ
2(12)
According to above formula relation, can find out that the number of large eigenwert is signal source number.Its characteristic of correspondence vector is respectively e
1, e
2..., e
p, e
p+1... e
nthen
Wherein E=[e
1, e
2..., e
p, e
p+1... e
n], Σ=diag{ λ
1, λ
2..., λ
p, σ
2..., σ
2.Definition signal subspace is E
s=[e
1, e
2..., e
p], noise subspace is E
n=[e
p+1..., e
n].Thus, signal subspace Es and noise subspace matrix En is obtained according to eigenwert characteristic of correspondence vector; Step 7 have selected the wherein a kind of MUSIC algorithm in the class algorithm of subspace, utilizes noise subspace matrix En and steering vector
orthogonality, ask its two-dimensional space spectral function, shown in (14):
Wherein, θ and
be respectively the angle of pitch and the position angle of source signal,
be with θ and
the array steering vector to received signal of change.Make θ and
change, finds out spatial spectrum
the angle corresponding to maximum point be the DOA of source signal.
Claims (2)
1., based on a method for hexagonal array coherent signal arrival direction estimation, its feature is provided by following steps:
Step one: sexangle antenna array elements number is N, exporting the fast umber of beats of data is K, then each array element exports the data matrix X complex matrix that (k)=AS (k)+N (k) is N × K, asks its conjugate matrices X
*(k), and the correlation matrix obtaining former array output data is R
xX=E{X (k) X (k)
h, H represents its conjugate transpose, and A is the steering vector matrix of array, and S (k) represents source signal vector matrix; N (k) represents that noise average that array exports be zero variance is σ
2additive white Gaussian noise, and uncorrelated with source signal;
Step 2: conjugation data matrix is for aerial array, the data being equivalent to the symmetrical virtual array of its conjugation export, for the array after conjugation virtual extended, three sub-hexagonal arrays can be divided into, first sub-hexagonal array is real hexagonal array, and second sub-hexagonal array is centered by reference array element, in submatrix array element wherein half belong to the part of true array, second half is the part of conjugation virtual array, and the 3rd sub-hexagonal array is classified as conjugate radical matroid row;
Step 3: according to the array element contained by every sub-hexagonal array of array after expansion, draw the output data matrix X of sub-hexagonal array
1(k), X
2(k), X
3(k);
Step 4: autocorrelation matrix and the cross-correlation matrix of obtaining three sub-hexagonal arrays respectively:
Step 5: construct the correlation matrix made new advances
Step 6: carry out svd to R, decompositing N number of eigenwert is λ
1>=λ
2>=...>=λ
p>=λ
p+1=...=λ
n=σ
2, by judging that the number of large eigenwert carrys out estimated signal source number, and obtain signal subspace Es and noise subspace matrix En respectively according to corresponding proper vector;
Step 7: utilize two-dimentional MUSIC algorithm to build spatial spectrum function,
θ and
be respectively the angle of pitch and the position angle of source signal,
be with θ and
the array steering vector to received signal of change.Make pitching angle theta at (0 °, 90 °) and position angle
in (0 °, 360 °) scope during change, find out spatial spectrum
angle corresponding to maximum point is the DOA of source signal.
2., according to a kind of method based on hexagonal array coherent signal arrival direction estimation described in claim 1, it is characterized in that: in described step 4, for coherent signal source, its correlation matrix R
xXorder can be less than source signal number, correct signal subspace can not be obtained, and new correlation matrix R carries out expansion reconstruct on the basis of former correlation matrix, makes its order arrive source signal number, the coherence of releasing source signal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310454248.0A CN104515969B (en) | 2013-09-29 | 2013-09-29 | Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310454248.0A CN104515969B (en) | 2013-09-29 | 2013-09-29 | Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104515969A true CN104515969A (en) | 2015-04-15 |
CN104515969B CN104515969B (en) | 2017-01-18 |
Family
ID=52791563
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310454248.0A Active CN104515969B (en) | 2013-09-29 | 2013-09-29 | Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104515969B (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105403856A (en) * | 2015-10-30 | 2016-03-16 | 西安电子科技大学 | DOA (direction of arrival) estimation method based on nested minimum redundant array |
CN107728105A (en) * | 2017-10-12 | 2018-02-23 | 天津津航计算技术研究所 | A kind of hexagonal array DOA algorithm for estimating based on phased-array technique |
CN108181557A (en) * | 2017-12-29 | 2018-06-19 | 上海交通大学 | A kind of method in determining ultrahigh frequency partial discharge signal orientation |
CN109541526A (en) * | 2018-12-29 | 2019-03-29 | 西北工业大学 | A kind of ring array direction estimation method using matrixing |
CN109582919A (en) * | 2018-11-28 | 2019-04-05 | 四川九洲电器集团有限责任公司 | Method for parameter estimation when a kind of sky based on uniform linear array |
CN109672462A (en) * | 2019-01-15 | 2019-04-23 | 西安电子科技大学 | Two-dimensional antenna array airspace filter method based on fourth order cumulant |
CN110244258A (en) * | 2019-06-12 | 2019-09-17 | 南京航空航天大学 | For extending DOA matrix method in double parallel battle array two dimension direction finding |
CN110531312A (en) * | 2019-08-29 | 2019-12-03 | 深圳市远翰科技有限公司 | A kind of DOA estimation method and system based on sparse symmetric matrix column |
CN111142097A (en) * | 2019-07-05 | 2020-05-12 | 武汉工程大学 | Target direction estimation method in amplitude modulation broadcast external radiation source radar |
CN111175693A (en) * | 2020-01-19 | 2020-05-19 | 河北科技大学 | Direction-of-arrival estimation method and direction-of-arrival estimation device |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101592721A (en) * | 2009-06-19 | 2009-12-02 | 中国人民解放军理工大学 | Method for estimating angle of arrival of coherent signal based on eigenvalue reconstruction |
CN102253363A (en) * | 2011-03-29 | 2011-11-23 | 西安交通大学 | Device for estimating two-dimensional direction of arrival (DOA) of coherent signals based on L array and method thereof |
-
2013
- 2013-09-29 CN CN201310454248.0A patent/CN104515969B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101592721A (en) * | 2009-06-19 | 2009-12-02 | 中国人民解放军理工大学 | Method for estimating angle of arrival of coherent signal based on eigenvalue reconstruction |
CN102253363A (en) * | 2011-03-29 | 2011-11-23 | 西安交通大学 | Device for estimating two-dimensional direction of arrival (DOA) of coherent signals based on L array and method thereof |
Non-Patent Citations (3)
Title |
---|
SALEM AKKAR ET AL.: ""Concentric Circular Array for DOAs Estimation of Coherent Sources with ESPRIT Algorithm"", 《2010 INTERNATIONAL CONFERENCE ON DESIGN & TECHNOLOGY OF INTEGRATED SYSTEMS IN NANOSCALE ERA》 * |
XIAOPENG YANG ET AL.: ""Effect of Geometry of Planar Antenna Arrays on Cramer-Rao Bounds for DOA Estimation"", 《IEEE 10TH INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING PROCEEDINGS》 * |
ZHI TIAN ET AL.: ""DOA ESTIMATION WITH HEXAGONAL ARRAYS"", 《ACOUSTICS, SPEECH AND SIGNAL PROCESSING, 1998. PROCEEDINGS OF THE 1998 IEEE INTERNATIONAL CONFERENCE ON》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105403856A (en) * | 2015-10-30 | 2016-03-16 | 西安电子科技大学 | DOA (direction of arrival) estimation method based on nested minimum redundant array |
CN107728105A (en) * | 2017-10-12 | 2018-02-23 | 天津津航计算技术研究所 | A kind of hexagonal array DOA algorithm for estimating based on phased-array technique |
CN108181557B (en) * | 2017-12-29 | 2020-05-05 | 上海交通大学 | Method for determining ultrahigh frequency partial discharge signal direction |
CN108181557A (en) * | 2017-12-29 | 2018-06-19 | 上海交通大学 | A kind of method in determining ultrahigh frequency partial discharge signal orientation |
CN109582919A (en) * | 2018-11-28 | 2019-04-05 | 四川九洲电器集团有限责任公司 | Method for parameter estimation when a kind of sky based on uniform linear array |
CN109541526A (en) * | 2018-12-29 | 2019-03-29 | 西北工业大学 | A kind of ring array direction estimation method using matrixing |
CN109541526B (en) * | 2018-12-29 | 2022-09-02 | 西北工业大学 | Circular ring array azimuth estimation method by using matrix transformation |
CN109672462A (en) * | 2019-01-15 | 2019-04-23 | 西安电子科技大学 | Two-dimensional antenna array airspace filter method based on fourth order cumulant |
CN110244258A (en) * | 2019-06-12 | 2019-09-17 | 南京航空航天大学 | For extending DOA matrix method in double parallel battle array two dimension direction finding |
CN110244258B (en) * | 2019-06-12 | 2022-10-04 | 南京航空航天大学 | Method for expanding DOA matrix in two-dimensional direction finding of double parallel arrays |
CN111142097A (en) * | 2019-07-05 | 2020-05-12 | 武汉工程大学 | Target direction estimation method in amplitude modulation broadcast external radiation source radar |
CN110531312B (en) * | 2019-08-29 | 2021-09-17 | 深圳市远翰科技有限公司 | DOA estimation method and system based on sparse symmetric array |
CN110531312A (en) * | 2019-08-29 | 2019-12-03 | 深圳市远翰科技有限公司 | A kind of DOA estimation method and system based on sparse symmetric matrix column |
CN111175693A (en) * | 2020-01-19 | 2020-05-19 | 河北科技大学 | Direction-of-arrival estimation method and direction-of-arrival estimation device |
Also Published As
Publication number | Publication date |
---|---|
CN104515969B (en) | 2017-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104515969B (en) | Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method | |
CN109655799B (en) | IAA-based covariance matrix vectorization non-uniform sparse array direction finding method | |
Zheng et al. | DOA estimation for coprime linear arrays: An ambiguity-free method involving full DOFs | |
CN104698433B (en) | DOA Estimation in Coherent Signal method based on single snapshot data | |
CN107037392B (en) | Degree-of-freedom increased type co-prime array direction-of-arrival estimation method based on compressed sensing | |
CN101592721B (en) | Eigenvalue reconstruction based method for estimating angle of arrival of coherent signal | |
CN103941221B (en) | Space stretching Electromagnetic Vector Sensor Array method for parameter estimation | |
CN104898085B (en) | A kind of dimensionality reduction MUSIC algorithms of polarization sensitive array parameter Estimation | |
CN108120967B (en) | Plane array DOA estimation method and equipment | |
CN104991236B (en) | A kind of single base MIMO radar not rounded signal coherence source Wave arrival direction estimating method | |
Yang et al. | A unified array geometry composed of multiple identical subarrays with hole-free difference coarrays for underdetermined DOA estimation | |
CN105335615B (en) | A kind of two dimension angular and polarization parameter combined estimation method of low complex degree | |
Zheng et al. | Two-dimensional DOA estimation for coprime planar array: A coarray tensor-based solution | |
CN103353588B (en) | Two-dimensional DOA (direction of arrival) angle estimation method based on antenna uniform planar array | |
CN102662158B (en) | Quick processing method for sensor antenna array received signals | |
CN103344940A (en) | Low-complexity DOA estimation method and system | |
WO2021068495A1 (en) | Degree of freedom-enhanced spatial spectrum estimation method based on block sampled tensor signal construction using planar co-prime array | |
CN109507636B (en) | Direction-of-arrival estimation method based on virtual domain signal reconstruction | |
CN112269172A (en) | Nested MIMO radar angle estimation method and device based on tensor structure | |
CN109521393A (en) | A kind of DOA estimation algorithm based on signal subspace revolving property | |
CN111965591A (en) | Direction-finding estimation method based on fourth-order cumulant vectorization DFT | |
CN107315161A (en) | Not rounded signal wave based on compressed sensing reaches bearing estimate method | |
CN103399308A (en) | Rapid estimation method of radar target angle under main lobe and side lobe jamming backgrounds | |
Liu et al. | Two dimensional direction of arrival estimation for co-prime L-shaped array using sparse reconstruction | |
Shi et al. | Symmetric sum coarray based co-prime MIMO configuration for direction of arrival estimation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |