CN110244258A - For extending DOA matrix method in double parallel battle array two dimension direction finding - Google Patents

For extending DOA matrix method in double parallel battle array two dimension direction finding Download PDF

Info

Publication number
CN110244258A
CN110244258A CN201910504096.8A CN201910504096A CN110244258A CN 110244258 A CN110244258 A CN 110244258A CN 201910504096 A CN201910504096 A CN 201910504096A CN 110244258 A CN110244258 A CN 110244258A
Authority
CN
China
Prior art keywords
matrix
doa
estimation
battle array
double parallel
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
Application number
CN201910504096.8A
Other languages
Chinese (zh)
Other versions
CN110244258B (en
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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201910504096.8A priority Critical patent/CN110244258B/en
Publication of CN110244258A publication Critical patent/CN110244258A/en
Application granted granted Critical
Publication of CN110244258B publication Critical patent/CN110244258B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Direction-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/02Direction-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/14Systems for determining direction or deviation from predetermined direction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Abstract

The invention discloses for extending DOA matrix method in the two dimension direction finding of double parallel battle array, the autocorrelation matrix and cross-correlation matrix of the reception data of double parallel battle array has been fully utilized in this method, construct the DOA matrix of an extension, then pass through the feature decomposition to DOA matrix, sense vector sum sense element to be estimated can be directly obtained, and thus obtains the two-dimentional DOA angle estimation of signal to be estimated.Traditional DOA matrix method is compared, method proposed by the present invention because the autocorrelation matrix and cross-correlation matrix of the reception data of double parallel battle array has been fully utilized there is better DOA to estimate performance.Method proposed by the present invention does not need spatial spectrum search, has lower algorithm complexity, and obtained DOA estimation angle is able to achieve automatic matching.

Description

For extending DOA matrix method in double parallel battle array two dimension direction finding
Technical field
The present invention relates to for extending DOA matrix method in the two dimension direction finding of double parallel battle array, and in particular to far-field signal source Angle estimating method belongs to array signal processing technology.
Background technique
Traditional DOA matrix method constructs DOA matrix according to the property of covariance matrix.Pass through the feature to DOA matrix It decomposes, sense vector sum sense element to be estimated can be directly obtained, and can thus estimate signal parameter, thus Multinomial search is avoided completely, and operand is smaller, but the auto-correlation information since array received signal not being fully utilized With cross-correlation information, so two dimension DOA angle estimation performance is lower.
Summary of the invention
The technical problems to be solved by the present invention are: provide for extending DOA matrix method in the two dimension direction finding of double parallel battle array, The auto-correlation information and cross-correlation information of array received signal has been fully utilized, has constructed an extension DOA matrix, solves double flat Arrival direction estimation problem under row battle array.
The present invention uses following technical scheme to solve above-mentioned technical problem:
For extending DOA matrix method in double parallel battle array two dimension direction finding, include the following steps:
Step 1, double parallel battle array, the uniform line temper that the double parallel battle array is parallel to each other by two are arranged on two-dimensional space Array XaAnd YaIt constitutes, when there is K irrelevant narrowbands to be incident on double parallel battle array with carrier signal, solves linear sub-arrays Xa And YaReceive the estimation of the autocorrelation matrix of signalWithLinear sub-arrays XaAnd YaReceive the cross-correlation matrix of signal EstimationAnd linear sub-arrays YaAnd XaReceive the estimation of the cross-correlation matrix of signal
Step 2, according to the estimation of autocorrelation matrixWithAnd the estimation of cross-correlation matrixWithObtain square Battle arrayWith
Step 3, according to the estimation of cross-correlation matrixWithAnd matrixWithDefine matrix R1And R2, and Building extension DOA matrix
Step 4, to extension DOA matrix R ' carry out Eigenvalues Decomposition, obtain two-dimentional DOA's according to characteristic value and feature vector Estimation.
As a preferred solution of the present invention, the estimation of autocorrelation matrix described in step 1WithAnd cross-correlation The estimation of matrixWithIt is obtained by following formula:
Wherein, L is number of snapshots, and x (t) indicates t moment linear sub-arrays XaReception signal, y (t) indicate t moment it is linear Subarray YaReception signal, ()HRepresenting matrix conjugate transposition.
As a preferred solution of the present invention, matrix described in step 2WithFormula it is as follows:
Wherein,WithIndicate autocorrelation matrix RxxAnd RyyEstimation, σ2For the variance of additive white Gaussian noise, IMFor The unit matrix of M × M, M are linear sub-arrays XaOr YaNumber of sensors.
As a preferred solution of the present invention, the formula that DOA matrix R ' is extended described in step 3 is as follows:
Wherein, R1And R2It is matrix,WithIndicate cross-correlation matrix RxyAnd RyxEstimation,WithIt is Matrix, ()+Representing matrix pseudoinverse.
As a preferred solution of the present invention, the detailed process of the step 4 are as follows:
DefinitionWithK is number of the irrelevant narrowband with carrier signal, To extension DOA matrix R ' carry out Eigenvalues Decomposition, matrix A is obtainedEAnd Φ, and have R ' AE=AEΦ, then according to the feature in Φ Value, obtains vkEstimated valueAnd then obtain DOA angle betakEstimation
By AEIt is divided into A and A Φ-1Two parts, after feature decomposition, this two-part estimation is respectivelyWithRecycling side DOA angle estimation is carried out to direction matrix to the Vandermonde feature of matrix, to obtain DOA angle [alpha]kEstimation
Wherein, It serves as reasonsObtained uk1Estimated value,It serves as reasons Obtained uk2Estimated value.
The invention adopts the above technical scheme compared with prior art, has following technical effect that
1, extension DOA matrix method proposed by the present invention, maintain traditional DOA matrix method can completely avoid it is multinomial While operand smaller advantage, the auto-correlation information and cross-correlation letter of array received signal has been fully utilized in formula search Breath constructs the DOA matrix of an extension, improves two-dimentional DOA angle estimation performance.
2, method proposed by the present invention does not need spatial spectrum search, has lower algorithm complexity.
3, the present invention is by the feature decomposition to DOA matrix, the angle DOA for obtaining the DOA angle of source signal, and obtaining Degree is able to achieve automatic matching.
Detailed description of the invention
Fig. 1 is array structure topological diagram of the invention.
Fig. 2 is the emulation scatter diagram of the method for the present invention.
Fig. 3 is the angle RMSE performance comparison figure of the method for the present invention and tradition DOA matrix method under the conditions of different signal-to-noise ratio.
Fig. 4 is the angle RMSE performance comparison figure of the method for the present invention and tradition DOA matrix method under the conditions of different snaps.
Fig. 5 is angle RMSE performance comparison figure of the method for the present invention under various information source number under the conditions of different signal-to-noise ratio.
Specific embodiment
Embodiments of the present invention are described below in detail, the example of the embodiment is shown in the accompanying drawings.Below by The embodiment being described with reference to the drawings is exemplary, and for explaining only the invention, and is not construed as limiting the claims.
Symbol indicates: using () in the present inventionTRepresenting matrix transposition, ()HRepresenting matrix conjugate transposition, () * indicate square Battle array conjugation, capital X representing matrix, lowercase x () indicate vector, IMIt is expressed as the unit matrix of M × M,It indicates Hadamard product, diag (v) indicate that the diagonal matrix constituted with element in v, E [] indicate it is expected Matrix Calculating, angle () Expression takes phase angle to operate.
One, data model
The two linear sub-arrays compositions being parallel to each other of signal receiving array as shown in Figure 1, two subarrays claim respectively For XaAnd Ya.Each submatrix shows M sensor, and two adjacent sensors spacing in X direction and two submatrix column pitch are d.Assuming that Space has K irrelevant narrowbands with carrier signal sk(t) (1≤k≤K) is incident on this array, and the angle with x-axis is αk, with The angle of y-axis is βk.Submatrix XaAnd YaIt is respectively as follows: in the corresponding output signal of t moment
X (t)=As (t)+nx(t)
Y (t)=A Φ s (t)+ny(t)
Wherein, nx(t) and nyIt (t) is the additive white Gaussian noise vector of two submatrixs, it is mutually indepedent with signal s (t).He Expression formula be respectively x (t)=[x1(t),…,xM(t)]T, y (t)=[y1 (t) ..., yM (t)]T, s (t)=[s1 (t),…,sK(t)]T, A=[a1,…,aK], and have
ak=[1, exp (j2 π dcos αk/λ),…,exp(j2π(M-1)dcosαk/λ)]T
Φ=diag [exp (j2 π dcos β1/λ),…,exp(j2πdcosβK/λ)]
Wherein, λ is wavelength.
Two, method derives
The autocorrelation matrix of the data x (t) received is Rxx, expression formula are as follows:
Rxx=E [x (t) xH(t)]=APAH2IM
Wherein P=E [s (t) sHIt (t)] is the covariance matrix of signal source, σ2For the variance of additive white Gaussian noise.
The autocorrelation matrix of the data y (t) received is Ryy, expression formula are as follows:
Ryy=E [y (t) yH(t)]=A Φ P ΦHAH2IM
=AP Φ ΦHAH2IM
=APAH2IM
Consider the independence of noise itself, and independently of signal, it is assumed that the cross-correlation matrix of y (t) and x (t) is Ryx, then:
Ryx=E [y (t) xH(t)]=A Φ PAH
It can similarly obtain, the cross-correlation matrix of x (t) and y (t) are as follows:
Rxy=E [x (t) yH(t)]=A Φ-1PAH
To RxxIt carries out Eigenvalues Decomposition (EVD), enables ε1,…,εKFor matrix RxxK big characteristic values, in the vacation of white noise It sets, can be averaged to obtain noise variance σ by a small characteristic value of M-K2Estimation.It then, can by removing the influence of noise To obtain:
Cxx=APAH=Rxx2IM
It can similarly obtain
Cyy=APAH=Ryy2IM
Definition
Definition
Therefore
R1=AEPAH
R2=AEΦPAH
According to the thought of DOA matrix method, the DOA matrix that can be defined as follows
Wherein
If A and P full rank, Φ is without identical diagonal element, then the K nonzero eigenvalue of DOA matrix R ' is equal in Φ K Diagonal element, and these are worth corresponding feature vector and are equal to corresponding sense vector, it may be assumed that
R′AE=AEΦ
DefinitionWithTo DOA matrix R ', feature decomposition is carried out, so that it may To matrix AEAnd Φ.According to the characteristic value in Φ, available vkEstimated valueAnd then obtain DOA angle betakEstimation:
According to AEDefinition, we are classified as A and A Φ-1Two parts, after feature decomposition, this two-part estimation difference ForWith
It estimatesWithAfterwards, ifA certain column ai, the Vandermonde feature of utilization orientation matrix is to direction square Battle array carries out DOA angle estimation.First to direction vector aiNormalization, makes its first term 1.Take angle (ai), estimate between its array Phase difference, finally estimate its DOA angle using least square method.Because of ai=[1, exp (j2 π dcos αi/λ),…,exp (j2π(M-1)dcosαi/λ)]T, it is possible to it obtains:
T=angle (ai)
=[0,2 π dcos αk/λ,…,2(M-1)πdcosαk/λ]T
Least square fitting is Bc=T, wherein c1=[c01,uk1]T, and have
It obtains
WhereinIt is exactly cos αk1Estimation, then it is availableEstimation:
Similarly, we can be fromIt obtainsTherefore DOA angle [alpha]kEstimation:
The method of the present invention step:
[1] autocorrelation matrix and the estimation of cross-correlation matrix of data x (t) and the y (t) received are asked
[2] influence of noise is removed to autocorrelation matrix, obtainedWith
[3] R is defined1And R2And construct extension DOA matrix
[4] to R ' carry out Eigenvalues Decomposition, the estimation of two-dimentional DOA is obtained according to characteristic value and feature vector.
Three, method analysis and emulation
Analysis of complexity is carried out to DOA angle estimating method of the invention, obtains the complexity of auto-correlation and cross-correlation matrix Degree is O { 4LM2, wherein L indicates to receive signal number of snapshots;It calculatesComplexity be O { 5M3};It calculatesAnswer Miscellaneous degree is O { 4M3};Complexity to R ' carry out feature decomposition is O { 8M3}.The total complexity of the algorithm asked is O { 4LM2+17M3}。
Auto-correlation information and the cross-correlation information of array received data has been fully utilized to construct one in method of the invention The DOA matrix of a extension, and traditional DOA matrix method, be not fully utilized array received data auto-correlation information and mutually Information is closed, therefore method of the invention has higher DOA angle estimation performance than traditional DOA matrix method.
Simulation result:
Assuming that three narrow band signal (α of space far-field1150 ° of)=(, 55 °), (α2260 ° of)=(, 65 °) and (α33) =(70 °, 75 °) are incident on array shown in Fig. 1, irrelevant between signal.The present invention uses 1000 Monte Carlo simulations Performance is estimated to assess DOA, and it is as follows to define root-mean-square error (RMSE) expression formula:
WhereinWithIndicate k-th of information source parameter estimation result, α in n-th Monte Carlo simulationkAnd βkIt indicates The parameter true value of k-th of information source.
Fig. 2 gives the scatter diagram of the extension DOA matrix method of double parallel battle array, and simulation parameter is L=500 and SNR= 20dB.This it appears that the two-dimentional DOA of information source from figure.
Fig. 3 gives traditional DOA matrix algorithm of double parallel battle array under the same conditions and extension DOA matrix algorithm with noise Than (SNR) variation angle estimation performance and with CRB (Cramér-Rao lower bound) performance comparison curve graph.Simulation parameter is set as double The array number M=8, number of snapshots L=500 of parallel battle array neutron array 1 and submatrix 2.From figure 3, it can be seen that extension DOA matrix method, Angle estimation performance with higher.
Fig. 4 is given under the conditions of identical signal-to-noise ratio, the traditional DOA matrix algorithm and extension DOA matrix algorithm of double parallel battle array The performance chart that changes with snap of angle estimation performance, signal-to-noise ratio is set as SNR=10dB.Figure 4, it is seen that With the increase of number of snapshots, the angle estimation performance for extending DOA matrix method is substantially better than traditional DOA matrix method.
Fig. 5 gives under the same conditions, when information source number changes, extends the angle estimation performance of DOA matrix method with letter It makes an uproar than the curve graph of (SNR) variation.From fig. 5, it can be seen that increasing with information source number, extends the angle estimation of DOA matrix method Performance is decreased obviously.
The above examples only illustrate the technical idea of the present invention, and this does not limit the scope of protection of the present invention, all According to the technical idea provided by the invention, any changes made on the basis of the technical scheme each falls within the scope of the present invention Within.

Claims (5)

1. for extending DOA matrix method in the two dimension direction finding of double parallel battle array, which comprises the steps of:
Step 1, double parallel battle array, the homogenous linear subarray that the double parallel battle array is parallel to each other by two are arranged on two-dimensional space XaAnd YaIt constitutes, when there is K irrelevant narrowbands to be incident on double parallel battle array with carrier signal, solves linear sub-arrays XaAnd Ya Receive the estimation of the autocorrelation matrix of signalWithLinear sub-arrays XaAnd YaReceive the estimation of the cross-correlation matrix of signalAnd linear sub-arrays YaAnd XaReceive the estimation of the cross-correlation matrix of signal
Step 2, according to the estimation of autocorrelation matrixWithAnd the estimation of cross-correlation matrixWithObtain matrixWith
Step 3, according to the estimation of cross-correlation matrixWithAnd matrixWithDefine matrix R1And R2, and construct expansion Open up DOA matrix
Step 4, to extension DOA matrix R ' carry out Eigenvalues Decomposition, estimating for two-dimentional DOA is obtained according to characteristic value and feature vector Meter.
2. according to claim 1 for extending DOA matrix method in the two dimension direction finding of double parallel battle array, which is characterized in that step 1 The estimation of the autocorrelation matrixWithAnd the estimation of cross-correlation matrixWithIt is obtained by following formula:
Wherein, L is number of snapshots, and x (t) indicates t moment linear sub-arrays XaReception signal, y (t) indicate t moment linear sub-arrays YaReception signal, ()HRepresenting matrix conjugate transposition.
3. according to claim 1 for extending DOA matrix method in the two dimension direction finding of double parallel battle array, which is characterized in that step 2 The matrixWithFormula it is as follows:
Wherein,WithIndicate autocorrelation matrix RxxAnd RyyEstimation, σ2For the variance of additive white Gaussian noise, IMFor M × M Unit matrix, M be linear sub-arrays XaOr YaNumber of sensors.
4. according to claim 1 for extending DOA matrix method in the two dimension direction finding of double parallel battle array, which is characterized in that step 3 The formula of the extension DOA matrix R ' is as follows:
Wherein, R1And R2It is matrix,WithIndicate cross-correlation matrix RxyAnd RyxEstimation,WithIt is matrix, (·)+Representing matrix pseudoinverse.
5. according to claim 1 for extending DOA matrix method in the two dimension direction finding of double parallel battle array, which is characterized in that described The detailed process of step 4 are as follows:
DefinitionWithK=1,2 ..., K, K is number of the irrelevant narrowband with carrier signal, to expansion DOA matrix R ' carry out Eigenvalues Decomposition is opened up, matrix A is obtainedEAnd Φ, and have R ' AE=AEΦ is obtained then according to the characteristic value in Φ To vkEstimated valueAnd then obtain DOA angle betakEstimation
By AEIt is divided into A and A Φ-1Two parts, after feature decomposition, this two-part estimation is respectivelyWithRecycle direction square The Vandermonde feature of battle array carries out DOA angle estimation to direction matrix, to obtain DOA angle [alpha]kEstimation
Wherein, It serves as reasonsObtained uk1Estimated value,It serves as reasonsIt obtains uk2Estimated value.
CN201910504096.8A 2019-06-12 2019-06-12 Method for expanding DOA matrix in two-dimensional direction finding of double parallel arrays Active CN110244258B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910504096.8A CN110244258B (en) 2019-06-12 2019-06-12 Method for expanding DOA matrix in two-dimensional direction finding of double parallel arrays

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910504096.8A CN110244258B (en) 2019-06-12 2019-06-12 Method for expanding DOA matrix in two-dimensional direction finding of double parallel arrays

Publications (2)

Publication Number Publication Date
CN110244258A true CN110244258A (en) 2019-09-17
CN110244258B CN110244258B (en) 2022-10-04

Family

ID=67886634

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910504096.8A Active CN110244258B (en) 2019-06-12 2019-06-12 Method for expanding DOA matrix in two-dimensional direction finding of double parallel arrays

Country Status (1)

Country Link
CN (1) CN110244258B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110888106A (en) * 2019-11-21 2020-03-17 南京航空航天大学 Angle and frequency joint estimation augmented DOA matrix method
CN111175692A (en) * 2020-01-17 2020-05-19 西北工业大学 Discrete sparse Bayesian DOA estimation method based on layered synthesis Lasso prior model
CN111337872A (en) * 2020-02-21 2020-06-26 南京航空航天大学 Generalized DOA matrix method for coherent information source direction finding
CN113126021A (en) * 2021-04-19 2021-07-16 电子科技大学 Single-snapshot two-dimensional DOA estimation method based on three parallel linear arrays

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104515969A (en) * 2013-09-29 2015-04-15 电子科技大学 Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method
CN106483493A (en) * 2016-09-13 2017-03-08 电子科技大学 A kind of sparse double parallel linear array and estimating two-dimensional direction-of-arrival method
CN107703478A (en) * 2017-10-27 2018-02-16 天津大学 Extension aperture arrival direction estimation method based on cross-correlation matrix
CN108020812A (en) * 2017-11-28 2018-05-11 天津大学 Arrival direction estimation method based on special three parallel linear array structures
CN109490820A (en) * 2018-11-13 2019-03-19 电子科技大学 A kind of arrival direction estimation method based on parallel nested battle array
CN109782218A (en) * 2019-02-01 2019-05-21 中国空间技术研究院 A kind of non-circular signal DOA estimation method of relevant distribution based on double parallel antenna array

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104515969A (en) * 2013-09-29 2015-04-15 电子科技大学 Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method
CN106483493A (en) * 2016-09-13 2017-03-08 电子科技大学 A kind of sparse double parallel linear array and estimating two-dimensional direction-of-arrival method
CN107703478A (en) * 2017-10-27 2018-02-16 天津大学 Extension aperture arrival direction estimation method based on cross-correlation matrix
CN108020812A (en) * 2017-11-28 2018-05-11 天津大学 Arrival direction estimation method based on special three parallel linear array structures
CN109490820A (en) * 2018-11-13 2019-03-19 电子科技大学 A kind of arrival direction estimation method based on parallel nested battle array
CN109782218A (en) * 2019-02-01 2019-05-21 中国空间技术研究院 A kind of non-circular signal DOA estimation method of relevant distribution based on double parallel antenna array

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110888106A (en) * 2019-11-21 2020-03-17 南京航空航天大学 Angle and frequency joint estimation augmented DOA matrix method
CN111175692A (en) * 2020-01-17 2020-05-19 西北工业大学 Discrete sparse Bayesian DOA estimation method based on layered synthesis Lasso prior model
CN111175692B (en) * 2020-01-17 2022-09-02 西北工业大学 Discrete sparse Bayesian DOA estimation method based on layered synthesis Lasso prior model
CN111337872A (en) * 2020-02-21 2020-06-26 南京航空航天大学 Generalized DOA matrix method for coherent information source direction finding
CN111337872B (en) * 2020-02-21 2023-08-01 南京航空航天大学 Generalized DOA matrix method for coherent source direction finding
CN113126021A (en) * 2021-04-19 2021-07-16 电子科技大学 Single-snapshot two-dimensional DOA estimation method based on three parallel linear arrays
CN113126021B (en) * 2021-04-19 2022-03-29 电子科技大学 Single-snapshot two-dimensional DOA estimation method based on three parallel linear arrays

Also Published As

Publication number Publication date
CN110244258B (en) 2022-10-04

Similar Documents

Publication Publication Date Title
CN110244258A (en) For extending DOA matrix method in double parallel battle array two dimension direction finding
CN110927659B (en) Method and system for estimating arbitrary array manifold DOA (direction of arrival) under cross-coupling condition and cross-coupling calibration
Wen et al. Angle estimation and mutual coupling self-calibration for ULA-based bistatic MIMO radar
Wen et al. Angle estimation for bistatic MIMO radar in the presence of spatial colored noise
CN106526530B (en) 2-L type array arrival direction estimation algorithm based on propagation operator
Zhang et al. Off-grid DOA estimation using array covariance matrix and block-sparse Bayesian learning
CN110244272B (en) Direction-of-arrival estimation method based on rank-denoising model
Porozantzidou et al. Azimuth and elevation angles estimation using 2-D MUSIC algorithm with an L-shape antenna
CN107544051A (en) Wave arrival direction estimating method of the nested array based on K R subspaces
CN109582919A (en) Method for parameter estimation when a kind of sky based on uniform linear array
CN109959891A (en) The dimensionality reduction spectrum peak search method of Space Angle and polarization parameter in L gusts of electromagnetic vector
CN111239678A (en) Two-dimensional DOA estimation method based on L-shaped array
Mao et al. An Improved DOA Estimation Algorithm Based on Wavelet Operator.
Zhang et al. Constant modulus waveform design for MIMO radar transmit beampattern with residual network
CN110895325B (en) Arrival angle estimation method based on enhanced quaternion multiple signal classification
CN108872930B (en) Extended aperture two-dimensional joint diagonalization DOA estimation method
CN113655444B (en) MIMO radar DOA estimation method based on re-weighting priori under array element failure
Xiao et al. A weighted forward-backward spatial smoothing DOA estimation algorithm based on TLS-ESPRIT
Wang et al. Improved DOA estimation based on real-valued array covariance using sparse Bayesian learning
CN112485755A (en) Angle estimation method based on 2 q-order nested MIMO array structure
Yang et al. Joint angular-frequency distribution estimation via spatial-temporal sparse sampling and low-rank matrix recovery
Zhao et al. Robust adaptive beamforming via improved worst-case performance optimization algorithm based on FDA-MIMO
CN111337872B (en) Generalized DOA matrix method for coherent source direction finding
Wen et al. Direction-of-arrival estimation for co-located multiple-input multiple-output radar using structural sparsity Bayesian learning
Tang et al. Iteratively reweighted lp norm minimization for DOD and DOA estimation in bistatic MIMO radar under impulsive noise

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant