CN104730491A - Virtual array DOA estimation method based on L type array - Google Patents
Virtual array DOA estimation method based on L type array Download PDFInfo
- Publication number
- CN104730491A CN104730491A CN201510100318.1A CN201510100318A CN104730491A CN 104730491 A CN104730491 A CN 104730491A CN 201510100318 A CN201510100318 A CN 201510100318A CN 104730491 A CN104730491 A CN 104730491A
- Authority
- CN
- China
- Prior art keywords
- array
- signal
- matrix
- submatrix
- virtual
- 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/80—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 ultrasonic, sonic or infrasonic waves
- G01S3/802—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)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a virtual array DOA estimation method based on an L type array. The method comprises the steps that 1, based on the shift invariant property, a subarray Zx and a subarray Zy of the L type array horizontally shift to obtain a virtual array Zx' and a virtual array Zy', rotation invariance of two sub signals is formed due to the shift invariant property of the subarrays, and the virtual subarray signals are equal to the L type subarray Zx and L type subarray Zy input signals multiplied by a twiddle factor respectively; 2, the output of the four subarrays are combined to form a virtual array output signal Z (t); 3, the signal subspace and the noise subspace can be described by decomposing the features of covariance matrixes output by the array, mutual correlation processing is carried out on the array output signal Z(t) to obtain Rzz, and eigenvalue decomposition is carried out to obtain signal subspaces; 4, the twiddle factor is solved through linear operation, and the signal wave arrival direction can be obtained through the diagonal element of the twiddle factor. According to the virtual array DOA estimation method, no spectral function needs to be calculated, the phenomenon that the wave arrival direction is indirectly calculated without searching for the peak value is avoided, the complexity is lowered, the equipment complexity and cost are reduced, and the positioning precision is high.
Description
Technical field
The invention belongs to Array Signal Processing arrival direction estimation technical field, in particular a kind of Virtual array DOA estimation method based on L-type array.
Background technology
Target Bearing Estimation (DOA) estimation is an important branch in the signal transacting fields such as sonar, radar, radio communication, medical imaging, microphone array column processing, simultaneously the basic problem that DOA solves determines to be in the locus of multiple interested echo signal in a certain space field, and namely each echo signal arrives the deflection of sensor array.Direction-finding method based on conventional beamformer scanning has intrinsic restriction, and be subject to the restriction of " Rayleigh criterion ", the resolution of estimation depends on array length, only when the separation angle in space between two information sources is greater than array aperture reciprocal, just can be resolved.In order to improve the resolution of basic matrix, when signal wavelength lambda must be, generally only has the aperture length increasing basic matrix, namely increase array element m or increase battle array spacing d, but increasing element number of array can improve equipment complexity and cost, increases array element distance and will cause secondary maximum again, for the restriction of actual conditions, the yardstick of array also can not do very large, therefore only relies on increase array aperture to reach and puies forward high-resolution way, be difficult to be suitable in Practical Project.In order to overcome this limitation, devise a kind of virtual array DOA estimation method based on L-type battle array, its basic thought is under the array case of limited dimensions, is obtained the virtual expansion of array by optimized algorithm, thus improves resolution.
Summary of the invention
The present invention is directed to the deficiencies in the prior art, propose a kind of virtual array DOA estimation method based on L-type battle array.
Technical scheme of the present invention is as follows:
Based on a virtual array DOA estimation method for L-type array, wherein, comprise the following steps:
Step 1: structure L-type array, determines the signal model of array received;
Receive acoustic pressure time-domain signal by 2M-1 array element in L-type array, this L-type array is the even linear array Z of M by array number in x-axis
xwith the even linear array Z that array number in y-axis is M
yform, wherein, 2M-1 is L-type array elements number, M be not less than 2 integer, d is array element distance.Hypothesis space has K information source to incide on array, and its 2-d direction finding is
be respectively the elevation angle and the position angle of a kth signal source.
Suppose that the signal number incided on this array is K, then signal x-axis, y-axis received by M array element is respectively respectively as shown in the formula (1) and formula (2):
x(t)=A
xs(t)+n(t) (1)
y(t)=A
ys(t)+n(t) (2)
In formula, s (t) is signal source matrix, and n (t) is noise matrix, A
x, A
y∈ C
m × K, be respectively the direction matrix in L-type array x-axis, y-axis, can be expressed as:
Step 2: obtain output signal matrix Z by ESPRIT constructing virtual array;
By the submatrix Z of ESPRIT by L battle array
x, Z
ycarrying out virtual expansion is submatrix Z
x', Z
y', the motion immovability due to submatrix defines the rotational invariance of two submatrix signals, i.e. Z
x'submatrix signal be actual submatrix Z
xinput signal be multiplied by twiddle factor φ
xobtain, Z
y'submatrix signal be actual submatrix Z
yinput signal be multiplied by twiddle factor φ
yobtain, first obtained the output signal of virtual submatrix by formula (1) and formula (2), then four submatrixs are exported and merged, form output signal matrix z (t) of whole array as shown in the formula (3):
Wherein
suppose that the direction of arrival of each information source is different, then
column vector between Line independent,
And
Wherein, matrix φ
x, φ
yfor the diagonal matrix of K × K, its diagonal element is that signal is respectively at Z
x, Z
yphase delay on array between any array element, diag represents diagonal matrix, and the element namely except principal diagonal is the square formation of zero.
X-axis direction homogenous linear submatrix Z is comprised such as formula (3) z (t)
xoutput signal x (t)
,y-axis direction homogenous linear submatrix Z
youtput signal y (t), Z
xthe virtual submatrix Z that translation obtains
x'output signal x'(t), Z
ythe virtual submatrix Z that translation obtains
y'output signal y'(t), the noise that each array received arrives is identical, virtual submatrix Z
x', Z
y'the even linear array of to be all array number be M.
Step 3: obtain correlation matrix R from array output signal matrix Z
z,
Signal subspace and noise subspace can obtain, shown in (4) by the feature decomposition that array exports the covariance matrix of Z:
R in formula
sfor the autocorrelation matrix of signal, σ
2for noise variance, I is unit matrix, the E [.] in formula (4), (.)
hbe expressed as mathematical expectation, conjugate transpose operation.
Step 4: by correlation matrix R
zdo feature decomposition, estimated signal number;
Array correlation matrix R
ztwo spaces can be divided into, namely
the eigenwert characteristic of correspondence vector E of K
s=[s
1, s
2... s
k] composition signal subspace, the non-singular matrix T that there is a K × K meets
and move invariant feature E due to array
scan be analyzed to 4 parts, E
x, E
y, E
x', E
y'∈ C
m × K, corresponding subarray is respectively Z
x, Z
y, Z
x', Z
y', shown in (5),
Step 5: structure φ
x, φ
ysimilar matrix F, H;
Formula (6) can be derived by formula (5):
E
x'=E
xT
-1φ
xT=E
xT E
y'=E
yT
-1φ
yT=E
yH (6)
Wherein, F=T
-1φ
xt, H=T
-1φ
yt, T are non-singular matrix, therefore F and φ
x, H and φ
yfor similar matrix, have identical eigenwert, and its eigenwert is twiddle factor φ
x, φ
ydiagonal element.
Step 6: least square method solves twiddle factor φ
x, φ
y, calculate direction of arrival
Twiddle factor φ is solved by least square method
x, φ
yshown in (7), just therefrom can draw the direction of arrival of signal.
for E
xpseudoinverse,
for E
ypseudoinverse, Eigenvalues Decomposition is carried out to F and obtains
obtain simultaneously
the u of estimated value
k, Eigenvalues Decomposition is carried out to H and obtains
obtain simultaneously
the v of estimated value
k.
can be estimated by formula (8):
Tool of the present invention has the following advantages: (1), compared with traditional coherent signal subspace algorithm, the present invention does not need to calculate spectral function, indirectly solves direction of arrival, reduce complexity without the need to search peak; (2) when number of sensors is determined, increase array aperture by virtual array, decrease equipment complexity and cost; (3) use the estimated value of linear operation direct solution sound source two dimension DOA, there is higher positioning precision.
Accompanying drawing explanation
Fig. 1 is L-type array structure schematic diagram;
Fig. 2 is the virtual array Z of translation along the y-axis direction
x';
Fig. 3 is the virtual array Z of translation along the x-axis direction
y'.
Fig. 4 is the actual direction of arrival of Fig. 4 sound source (/ °).
Fig. 5 is MUSIC algorithm estimated performance (SNR=20dB) in three information source situations
Fig. 6 is two-dimensional virtual array algorithm estimated performance (SNR=20dB) in three information source situations
Embodiment
Below in conjunction with drawings and Examples, the invention will be further described:
Step 1: structure L-type array, determines the signal model of array received;
As shown in Figure 1, receive acoustic pressure time-domain signal by 2M-1 array element in L-type array, this L-type array is the even linear array Z of M by array number in x-axis
xwith the even linear array Z that array number in y-axis is M
yform, wherein, 2M-1 is L-type array elements number, M be not less than 2 integer, d is array element distance.Hypothesis space has K information source to incide on array, and its 2-d direction finding is
be respectively the elevation angle and the position angle of a kth signal source.
Suppose that the signal number incided on this array is K, then signal x-axis, y-axis received by M array element is respectively respectively as shown in the formula (1) and formula (2):
x(t)=A
xs(t)+n(t) (1)
y(t)=A
ys(t)+n(t) (2)
In formula, s (t) is signal source matrix, and n (t) is noise matrix, A
x, A
y∈ C
m × K, be respectively the direction matrix in L-type array x-axis, y-axis, can be expressed as:
Step 2: obtain output signal matrix Z by ESPRIT constructing virtual array;
As shown in Figure 2 and Figure 3, by the submatrix Z of ESPRIT by L battle array
x, Z
ycarrying out virtual expansion is submatrix Z
x', Z
y', the motion immovability due to submatrix defines the rotational invariance of two submatrix signals, i.e. Z
x'submatrix signal be actual submatrix Z
xinput signal be multiplied by twiddle factor φ
xobtain, Z
y'submatrix signal be actual submatrix Z
yinput signal be multiplied by twiddle factor φ
yobtain, first obtained the output signal of virtual submatrix by formula (1) and formula (2), then four submatrixs are exported and merged, form output signal matrix z (t) of whole array as shown in the formula (3):
Its
suppose that the direction of arrival of each information source is different, then
column vector between Line independent,
And
Wherein, matrix φ
x, φ
yfor the diagonal matrix of K × K, its diagonal element is that signal is respectively at Z
x, Z
yphase delay on array between any array element, diag represents diagonal matrix, and the element namely except principal diagonal is the square formation of zero.
X-axis direction homogenous linear submatrix Z is comprised such as formula (3) z (t)
xoutput signal x (t)
,y-axis direction homogenous linear submatrix Z
youtput signal y (t), Z
xthe virtual submatrix Z that translation obtains
x'output signal x'(t), Z
ythe virtual submatrix Z that translation obtains
y'output signal y'(t), the noise that each array received arrives is identical, virtual submatrix Z
x', Z
y'the even linear array of to be all array number be M.
Step 3: obtain correlation matrix R from array output signal matrix Z
z,
Signal subspace and noise subspace can obtain, shown in (4) by the feature decomposition that array exports the covariance matrix of Z:
R in formula
sfor the autocorrelation matrix of signal, σ
2for noise variance, I is unit matrix, the E [.] in formula (4), (.)
hbe expressed as mathematical expectation, conjugate transpose operation.
Step 4: by correlation matrix R
zdo feature decomposition, estimated signal number;
Array correlation matrix R
ztwo spaces can be divided into, namely
the eigenwert characteristic of correspondence vector E of K
s=[s
1, s
2... s
k] composition signal subspace, the non-singular matrix T that there is a K × K meets
and move invariant feature E due to array
scan be analyzed to 4 parts, E
x, E
y, E
x', E
y'∈ C
m × K, corresponding subarray is respectively Z
x, Z
y, Z
x', Z
y', shown in (5),
Step 5: structure φ
x, φ
ysimilar matrix F, H;
Formula (6) can be derived by formula (5):
E
x'=E
xT
-1φ
xT=E
xT E
y'=E
yT
-1φ
yT=E
yH (6)
Wherein, F=T
-1φ
xt, H=T
-1φ
yt, T are non-singular matrix, therefore F and φ
x, H and φ
yfor similar matrix, have identical eigenwert, and its eigenwert is twiddle factor φ
x, φ
ydiagonal element.Step 6: least square method solves twiddle factor φ
x, φ
y, calculate direction of arrival
Twiddle factor φ is solved by least square method
x, φ
yshown in (7), just therefrom can draw the direction of arrival of signal.
for E
xpseudoinverse,
for E
ypseudoinverse, Eigenvalues Decomposition is carried out to F and obtains
obtain simultaneously
the u of estimated value
k, Eigenvalues Decomposition is carried out to H and obtains
obtain simultaneously
the v of estimated value
k.
can be estimated by formula (8):
Step 7: based on the Dynamic simulation result of the virtual array DOA estimation method of L-type battle array;
Simulated conditions is: the L-type array longitudinal axis of employing, transverse axis are all the even linear array of 8 omnidirectional's array element compositions, the aerial velocity of propagation of sound is c=340 m/s, and frequency of source is f=3000 Hz, and array element distance gets d=λ/2, i.e. d=5.7 cm, signal to noise ratio snr=20 dB.Contrast the DOA estimated performance based on MUSIC algorithm and virtual array algorithm, analogue simulation three sound sources, actual direction of arrival is (10 °, 15 °), (30 °, 25 °), (50 °, 35 °), as shown in Figure 4, in figure, stain is the real incident angle of sound source.Use MUSIC algorithm to estimate the direction of arrival of sound source by search peak-to-peak value, as shown in Figure 5, the direction of arrival of 1 sound source can only be estimated, DOA is (28 °, 24 °), and secondary lobe is more, positioning precision is fuzzyyer, accurately can not identify the accurate location of sound source.By the relation between sound source position and microphone, the ripple using the virtual array DOA estimation method based on L-type battle array directly to calculate three sound sources reaches estimated value and is (10.01 °, 15.07 °), (30 °, 25.01 °), (49.98 °, 34.99 °), as shown in Figure 6, in figure, stain is the DOA of estimation, compared with real direction of arrival, the DOA calculated changes between ± 1 °, error is less, and the DOA estimation method positioning precision based on L-type battle array virtual array algorithm is relatively high.
Can find out in result, compared with MUSIC algorithm, virtual array algorithm does not need to calculate spectral function, indirectly solves direction of arrival, reduce complexity without the need to search peak; When number of sensors is determined, increase array aperture by virtual array, decrease equipment complexity and cost; Use the estimated value of linear operation direct solution sound source two dimension DOA, there is higher positioning precision.
Claims (1)
1., based on a virtual array DOA estimation method for L-type array, it is characterized in that: the method comprises the following steps:
Step 1: structure L-type array, determines the signal model of array received;
Receive acoustic pressure time-domain signal by 2M-1 array element in L-type array, this L-type array is the even linear array Z of M by array number in x-axis
xwith the even linear array Z that array number in y-axis is M
yform, wherein, 2M-1 is L-type array elements number, M be not less than 2 integer, d is array element distance; Hypothesis space has K information source to incide on array, and its 2-d direction finding is
be respectively the elevation angle and the position angle of a kth signal source;
Suppose that the signal number incided on this array is K, then signal x-axis, y-axis received by M array element is respectively respectively as shown in the formula (1) and formula (2):
x(t)=A
xs(t)+n(t) (1)
y(t)=A
ys(t)+n(t) (2)
In formula, s (t) is signal source matrix, and n (t) is noise matrix, A
x, A
y∈ C
m × K, be respectively the direction matrix in L-type array x-axis, y-axis, can be expressed as:
Step 2: obtain output signal matrix Z by ESPRIT constructing virtual array;
By the submatrix Z of ESPRIT by L battle array
x, Z
ycarrying out virtual expansion is submatrix Z
x', Z
y', the motion immovability due to submatrix defines the rotational invariance of two submatrix signals, i.e. Z
x'submatrix signal be actual submatrix Z
xinput signal be multiplied by twiddle factor φ
xobtain, Z
y'submatrix signal be actual submatrix Z
yinput signal be multiplied by twiddle factor φ
yobtain, first obtained the output signal of virtual submatrix by formula (1) and formula (2), then four submatrixs are exported and merged, form output signal matrix z (t) of whole array as shown in the formula (3):
Wherein
suppose that the direction of arrival of each information source is different, then
column vector between Line independent, and
Wherein, matrix φ
x,φ
yfor the diagonal matrix of K × K, its diagonal element is that signal is respectively at Z
x, Z
yphase delay on array between any array element, diag represents diagonal matrix, and the element namely except principal diagonal is the square formation of zero;
X-axis direction homogenous linear submatrix Z is comprised such as formula (3) z (t)
xoutput signal x (t)
,y-axis direction homogenous linear submatrix Z
youtput signal y (t), Z
xthe virtual submatrix Z that translation obtains
x'output signal x'(t), Z
ythe virtual submatrix Z that translation obtains
y'output signal y'(t), the noise that each array received arrives is identical, virtual submatrix Z
x',z
y'the even linear array of to be all array number be M;
Step 3: obtain correlation matrix R from array output signal matrix Z
z,
Signal subspace and noise subspace can obtain, shown in (4) by the feature decomposition that array exports the covariance matrix of Z:
R in formula
sfor the autocorrelation matrix of signal, σ
2for noise variance, I is unit matrix, the E [.] in formula (4), (.)
hbe expressed as mathematical expectation, conjugate transpose operation;
Step 4: by correlation matrix R
zdo feature decomposition, estimated signal number;
Array correlation matrix R
ztwo spaces can be divided into, namely
the eigenwert characteristic of correspondence vector E of K
s=[s
1, s
2... s
k] composition signal subspace, the non-singular matrix T that there is a K × K meets
, and move invariant feature E due to array
scan be analyzed to 4 parts, E
x, E
y, E
x', E
y'∈ C
m × K, corresponding subarray is respectively Z
x, Z
y, Z
x', Z
y', shown in (5),
Step 5: structure φ
x, φ
ysimilar matrix F, H;
Formula (6) can be derived by formula (5):
E
x'=E
xT
-1φ
xT=E
xT E
y'=E
yT
-1φ
yT=E
yH (6)
Wherein, F=T
-1φ
xt, H=T
-1φ
yt, T are non-singular matrix, therefore F and φ
x, H and φ
yfor similar matrix, have identical eigenwert, and its eigenwert is twiddle factor φ
x, φ
ydiagonal element;
Step 6: least square method solves twiddle factor φ
x, φ
y, calculate direction of arrival
Twiddle factor φ is solved by least square method
x, φ
yshown in (7), just therefrom can draw the direction of arrival of signal;
for E
xpseudoinverse
, for E
ypseudoinverse, Eigenvalues Decomposition is carried out to F and obtains
obtain simultaneously
sin θ
kthe u of estimated value
k, Eigenvalues Decomposition is carried out to H and obtains
obtain simultaneously
sin θ
kthe v of estimated value
k;
can be estimated by formula (8):
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510100318.1A CN104730491B (en) | 2015-03-06 | 2015-03-06 | A kind of virtual array DOA estimation method based on L-type battle array |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510100318.1A CN104730491B (en) | 2015-03-06 | 2015-03-06 | A kind of virtual array DOA estimation method based on L-type battle array |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104730491A true CN104730491A (en) | 2015-06-24 |
CN104730491B CN104730491B (en) | 2017-05-31 |
Family
ID=53454570
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510100318.1A Expired - Fee Related CN104730491B (en) | 2015-03-06 | 2015-03-06 | A kind of virtual array DOA estimation method based on L-type battle array |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104730491B (en) |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105242237A (en) * | 2015-09-09 | 2016-01-13 | 电子科技大学 | Electromagnetic vector array parameter estimation method based on compressed sensing |
CN106154220A (en) * | 2016-06-20 | 2016-11-23 | 陕西理工学院 | L-type simplifies acoustic vector-sensor array row multiparameter Combined estimator quaternary counting method |
CN106443573A (en) * | 2016-09-14 | 2017-02-22 | 清华大学 | Angle-of-arrival single snapshot quick estimation method and device for massive antenna array signals |
CN106526530A (en) * | 2016-09-30 | 2017-03-22 | 天津大学 | Propagation operator-based 2-L type array two-dimensional DOA estimation algorithm |
CN106772224A (en) * | 2016-12-28 | 2017-05-31 | 武汉大学 | A kind of L-type array estimating two-dimensional direction-of-arrival algorithm of use time frequency analysis |
CN106950529A (en) * | 2017-03-21 | 2017-07-14 | 陕西理工大学 | Acoustic vector near field sources ESPRIT and MUSIC method for parameter estimation |
CN107422317A (en) * | 2017-08-03 | 2017-12-01 | 中国人民解放军空军工程大学 | Low angle target arrival direction estimation method based on smoothing matrix collection |
CN108008348A (en) * | 2017-11-16 | 2018-05-08 | 华南理工大学 | Underwater Wave arrival direction estimating method and device based on adjustable angle even linear array |
CN108181557A (en) * | 2017-12-29 | 2018-06-19 | 上海交通大学 | A kind of method in determining ultrahigh frequency partial discharge signal orientation |
CN108414967A (en) * | 2018-04-11 | 2018-08-17 | 华南理工大学 | Based on L gusts of underwater two-dimension Wave arrival direction estimating method and device of angle adjustable double |
CN108957391A (en) * | 2018-07-24 | 2018-12-07 | 北京理工大学 | A kind of estimating two-dimensional direction-of-arrival method of the inverted-L antenna battle array based on nested array |
CN109239657A (en) * | 2018-09-05 | 2019-01-18 | 南京航空航天大学 | Load the radiation source high-precision locating method under nested battle array unmanned aerial vehicle platform |
CN109375152A (en) * | 2018-09-05 | 2019-02-22 | 南京航空航天大学 | The DOA and polarization combined estimation method of L gusts of electromagnetic vector nesting lower low complex degrees |
CN109521392A (en) * | 2018-10-24 | 2019-03-26 | 华南理工大学 | Underwater one-dimensional DOA estimation method and device based on non-circular signal and L-type linear array |
CN109633600A (en) * | 2019-02-28 | 2019-04-16 | 电子科技大学 | A kind of DOA estimation method of minimum-redundancy linear arrays MIMO-OTHR |
CN109725292A (en) * | 2019-03-05 | 2019-05-07 | 中国电子科技集团公司第三研究所 | Underwater operation multiple target high-precision ultra-short baseline localization method and device |
CN109932680A (en) * | 2019-04-04 | 2019-06-25 | 哈尔滨工程大学 | A kind of non-circular method for estimating signal wave direction based on the relatively prime array of translation |
CN110018438A (en) * | 2019-04-23 | 2019-07-16 | 北京邮电大学 | A kind of Wave arrival direction estimating method and device |
CN110208738A (en) * | 2019-06-13 | 2019-09-06 | 哈尔滨工业大学 | Signal frequency and two dimension DOA combined estimation method based on array modulation wide-band transducer |
CN110286351A (en) * | 2019-07-12 | 2019-09-27 | 电子科技大学 | A kind of arrival direction estimation method and device based on L-type nesting battle array |
CN110673085A (en) * | 2019-09-25 | 2020-01-10 | 南京航空航天大学 | Coherent information source direction finding method based on fast convergence parallel factor under uniform area array |
CN111413667A (en) * | 2020-03-25 | 2020-07-14 | 北京邮电大学 | Method and device for determining direction of arrival of signal and electronic equipment |
CN111965599A (en) * | 2020-07-03 | 2020-11-20 | 重庆大学 | Sound source identification method for two-dimensional dynamic grid compressed beam forming |
CN113253213A (en) * | 2021-04-07 | 2021-08-13 | 哈尔滨工业大学(威海) | Method and system for measuring ionosphere echo space distribution characteristics of high-frequency ground wave radar |
CN113296056A (en) * | 2021-05-10 | 2021-08-24 | 华中科技大学 | Sound array configuration and sound source positioning method and system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4965732A (en) * | 1985-11-06 | 1990-10-23 | The Board Of Trustees Of The Leland Stanford Junior University | Methods and arrangements for signal reception and parameter estimation |
CN103364762A (en) * | 2013-07-13 | 2013-10-23 | 西安电子科技大学 | Estimation method for arriving direction of monostatic MIMO radar based on random array manifolds |
CN103941222A (en) * | 2014-03-12 | 2014-07-23 | 西安电子科技大学 | Parameter estimation method based on Rayleigh-Ritz entropy dimensionality reduction MUSIC algorithm |
CN104035069A (en) * | 2014-06-18 | 2014-09-10 | 西安交通大学 | Partly-corrected symmetrical uniform linear array based narrow-band near-field signal source positioning method |
-
2015
- 2015-03-06 CN CN201510100318.1A patent/CN104730491B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4965732A (en) * | 1985-11-06 | 1990-10-23 | The Board Of Trustees Of The Leland Stanford Junior University | Methods and arrangements for signal reception and parameter estimation |
CN103364762A (en) * | 2013-07-13 | 2013-10-23 | 西安电子科技大学 | Estimation method for arriving direction of monostatic MIMO radar based on random array manifolds |
CN103941222A (en) * | 2014-03-12 | 2014-07-23 | 西安电子科技大学 | Parameter estimation method based on Rayleigh-Ritz entropy dimensionality reduction MUSIC algorithm |
CN104035069A (en) * | 2014-06-18 | 2014-09-10 | 西安交通大学 | Partly-corrected symmetrical uniform linear array based narrow-band near-field signal source positioning method |
Non-Patent Citations (3)
Title |
---|
J.K.OH ET.AL: "《Optimized DOA estimation for an array antenna》", 《SPRINGER LINK》 * |
唐建红 等: "《快速四阶累积量旋转不变子空间算法》", 《西安交通大学学报》 * |
顾陈: "《高分辨阵列信号处理方法研究》", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105242237A (en) * | 2015-09-09 | 2016-01-13 | 电子科技大学 | Electromagnetic vector array parameter estimation method based on compressed sensing |
CN105242237B (en) * | 2015-09-09 | 2017-08-08 | 电子科技大学 | A kind of electromagnetic vector array parameter method of estimation based on compressed sensing |
CN106154220A (en) * | 2016-06-20 | 2016-11-23 | 陕西理工学院 | L-type simplifies acoustic vector-sensor array row multiparameter Combined estimator quaternary counting method |
CN106154220B (en) * | 2016-06-20 | 2019-10-18 | 陕西理工大学 | L-type simplifies acoustic vector-sensor array column multi-parameter Combined estimator quaternary counting method |
CN106443573B (en) * | 2016-09-14 | 2018-12-28 | 清华大学 | The direction of arrival list snap method for quick estimating and device of large-scale antenna array signal |
CN106443573A (en) * | 2016-09-14 | 2017-02-22 | 清华大学 | Angle-of-arrival single snapshot quick estimation method and device for massive antenna array signals |
CN106526530A (en) * | 2016-09-30 | 2017-03-22 | 天津大学 | Propagation operator-based 2-L type array two-dimensional DOA estimation algorithm |
CN106772224A (en) * | 2016-12-28 | 2017-05-31 | 武汉大学 | A kind of L-type array estimating two-dimensional direction-of-arrival algorithm of use time frequency analysis |
CN106950529A (en) * | 2017-03-21 | 2017-07-14 | 陕西理工大学 | Acoustic vector near field sources ESPRIT and MUSIC method for parameter estimation |
CN106950529B (en) * | 2017-03-21 | 2019-07-05 | 陕西理工大学 | Acoustic vector near field sources ESPRIT and MUSIC method for parameter estimation |
CN107422317A (en) * | 2017-08-03 | 2017-12-01 | 中国人民解放军空军工程大学 | Low angle target arrival direction estimation method based on smoothing matrix collection |
CN108008348A (en) * | 2017-11-16 | 2018-05-08 | 华南理工大学 | Underwater Wave arrival direction estimating method and device based on adjustable angle even linear array |
CN108181557A (en) * | 2017-12-29 | 2018-06-19 | 上海交通大学 | A kind of method in determining ultrahigh frequency partial discharge signal orientation |
CN108181557B (en) * | 2017-12-29 | 2020-05-05 | 上海交通大学 | Method for determining ultrahigh frequency partial discharge signal direction |
CN108414967A (en) * | 2018-04-11 | 2018-08-17 | 华南理工大学 | Based on L gusts of underwater two-dimension Wave arrival direction estimating method and device of angle adjustable double |
CN108957391A (en) * | 2018-07-24 | 2018-12-07 | 北京理工大学 | A kind of estimating two-dimensional direction-of-arrival method of the inverted-L antenna battle array based on nested array |
CN109375152B (en) * | 2018-09-05 | 2020-08-07 | 南京航空航天大学 | Low-complexity DOA and polarization joint estimation method under electromagnetic vector nested L array |
CN109239657A (en) * | 2018-09-05 | 2019-01-18 | 南京航空航天大学 | Load the radiation source high-precision locating method under nested battle array unmanned aerial vehicle platform |
CN109375152A (en) * | 2018-09-05 | 2019-02-22 | 南京航空航天大学 | The DOA and polarization combined estimation method of L gusts of electromagnetic vector nesting lower low complex degrees |
CN109521392A (en) * | 2018-10-24 | 2019-03-26 | 华南理工大学 | Underwater one-dimensional DOA estimation method and device based on non-circular signal and L-type linear array |
CN109521392B (en) * | 2018-10-24 | 2022-12-16 | 华南理工大学 | Underwater one-dimensional DOA estimation method and device based on non-circular signal and L-shaped linear array |
CN109633600A (en) * | 2019-02-28 | 2019-04-16 | 电子科技大学 | A kind of DOA estimation method of minimum-redundancy linear arrays MIMO-OTHR |
CN109633600B (en) * | 2019-02-28 | 2020-09-08 | 电子科技大学 | DOA estimation method of minimum redundant linear array MIMO-OTHR |
CN109725292B (en) * | 2019-03-05 | 2021-04-23 | 中国电子科技集团公司第三研究所 | Multi-target high-precision ultra-short baseline positioning method and device for underwater operation |
CN109725292A (en) * | 2019-03-05 | 2019-05-07 | 中国电子科技集团公司第三研究所 | Underwater operation multiple target high-precision ultra-short baseline localization method and device |
CN109932680A (en) * | 2019-04-04 | 2019-06-25 | 哈尔滨工程大学 | A kind of non-circular method for estimating signal wave direction based on the relatively prime array of translation |
CN109932680B (en) * | 2019-04-04 | 2022-05-20 | 哈尔滨工程大学 | Non-circular signal direction of arrival estimation method based on translational co-prime array |
CN110018438A (en) * | 2019-04-23 | 2019-07-16 | 北京邮电大学 | A kind of Wave arrival direction estimating method and device |
CN110018438B (en) * | 2019-04-23 | 2020-09-25 | 北京邮电大学 | Direction-of-arrival estimation method and device |
CN110208738B (en) * | 2019-06-13 | 2022-12-06 | 哈尔滨工业大学 | Signal frequency and two-dimensional DOA joint estimation method based on array modulation broadband converter |
CN110208738A (en) * | 2019-06-13 | 2019-09-06 | 哈尔滨工业大学 | Signal frequency and two dimension DOA combined estimation method based on array modulation wide-band transducer |
CN110286351A (en) * | 2019-07-12 | 2019-09-27 | 电子科技大学 | A kind of arrival direction estimation method and device based on L-type nesting battle array |
CN110673085A (en) * | 2019-09-25 | 2020-01-10 | 南京航空航天大学 | Coherent information source direction finding method based on fast convergence parallel factor under uniform area array |
CN111413667A (en) * | 2020-03-25 | 2020-07-14 | 北京邮电大学 | Method and device for determining direction of arrival of signal and electronic equipment |
CN111965599A (en) * | 2020-07-03 | 2020-11-20 | 重庆大学 | Sound source identification method for two-dimensional dynamic grid compressed beam forming |
CN113253213A (en) * | 2021-04-07 | 2021-08-13 | 哈尔滨工业大学(威海) | Method and system for measuring ionosphere echo space distribution characteristics of high-frequency ground wave radar |
CN113253213B (en) * | 2021-04-07 | 2024-01-30 | 哈尔滨工业大学(威海) | Method and system for measuring ionosphere echo spatial distribution characteristics of high-frequency ground wave radar |
CN113296056A (en) * | 2021-05-10 | 2021-08-24 | 华中科技大学 | Sound array configuration and sound source positioning method and system |
Also Published As
Publication number | Publication date |
---|---|
CN104730491B (en) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104730491A (en) | Virtual array DOA estimation method based on L type array | |
CN103353596B (en) | Wave beam space domain meter wave radar height measurement method based on compressed sensing | |
CN111123192B (en) | Two-dimensional DOA positioning method based on circular array and virtual extension | |
CN103353595B (en) | Meter wave radar height measurement method based on array interpolation compression perception | |
CN103383452B (en) | Distributive array target angle-of-arrival estimation method | |
CN104020469B (en) | A kind of MIMO radar distance-angle two-dimensional super-resolution rate imaging algorithm | |
CN110927661A (en) | Single-basis expansion co-prime array MIMO radar DOA estimation method based on MUSIC algorithm | |
CN106950529B (en) | Acoustic vector near field sources ESPRIT and MUSIC method for parameter estimation | |
Baig et al. | Comparison of direction of arrival (DOA) estimation techniques for closely spaced targets | |
Zhang et al. | Two-dimensional direction of arrival estimation for coprime planar arrays via polynomial root finding technique | |
CN109143153A (en) | A kind of super nested array Wave arrival direction estimating method based on sparse reconstruct | |
Zhang et al. | A rank-reduction based 2-D DOA estimation algorithm for three parallel uniform linear arrays | |
CN104730513A (en) | Multistage sub-array focusing MVDR wave beam forming method | |
CN106997037A (en) | Acoustic vector-sensor array column space rotates decorrelation LMS angle-of- arrival estimation method | |
CN106405487A (en) | General spatial spectrum estimation method based on extended ESPRIT | |
CN104931923A (en) | Grid iterative estimation of signal parameters via rotational invariance techniques (ESPRIT), namely, extensible rapid estimation algorithm capable of being used for uniform circular array 2-dimensional direction of arrival (2D DOA) | |
CN112462363A (en) | Coherent target parameter estimation method for non-uniform sparse polarization array | |
CN108614234B (en) | Direction-of-arrival estimation method based on multi-sampling snapshot co-prime array received signal fast Fourier inverse transformation | |
CN106908754A (en) | L-type acoustic vector-sensor array row ESPRIT decorrelation LMS method for parameter estimation | |
Reaz et al. | A comprehensive analysis and performance evaluation of different direction of arrival estimation algorithms | |
Bhargav et al. | Compressive sensing based doa estimation for multi-path environment | |
CN114563760B (en) | Second-order super-beam forming method, equipment and medium based on SCA array | |
CN105703841B (en) | A kind of separation method of multipath propagation broadband active acoustical signal | |
Chen et al. | Two-dimensional DOA estimation for planar array using a successive propagator method | |
CN112666558B (en) | Low-complexity MUSIC direction finding method and device suitable for automobile FMCW radar |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170531 Termination date: 20190306 |