CN103605108A - High-precision remote direction estimation method of acoustic vector array - Google Patents
High-precision remote direction estimation method of acoustic vector array Download PDFInfo
- Publication number
- CN103605108A CN103605108A CN201310322324.2A CN201310322324A CN103605108A CN 103605108 A CN103605108 A CN 103605108A CN 201310322324 A CN201310322324 A CN 201310322324A CN 103605108 A CN103605108 A CN 103605108A
- Authority
- CN
- China
- Prior art keywords
- vibration velocity
- acoustic
- direction estimation
- estimation method
- sensor array
- 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
Images
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
-
- 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
-
- 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/86—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 with means for eliminating undesired waves, e.g. disturbing noises
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 provides a high-precision remote direction estimation method of an acoustic vector array. (1) Narrow-band filtering is performed on a received signal of a two-dimensional vector hydrophone array, and a narrow-band output signal on a frequency point to be processed is obtained; (2) two vibration velocity components of the two-dimensional vector hydrophone array that are orthogonal to each other are linearly combined in a complex domain, and are converted into two new vibration velocity output components; (3) maximum likelihood estimation is used to determine sound pressure snapshoted by a certain number of samples and a cross covariance matrix pair of the double vibration velocity components in the complex domain; and (4) applying an ESPRIT algorithm of a matrix pencil to calculate a rotation invariant factor between the sound pressure vibration velocity cross covariance matrix pair based on the generalized acoustic energy flow, thereby performing direction estimation. The high-precision remote direction estimation method of the acoustic vector array can perform high-precision direction estimation under the circumstance of any array type and even an unknown array type, and the phenomenon that signals in certain directions are severely weakened and even completely shielded does not occur. In addition, the high-precision remote direction estimation method has a relatively low processable signal to noise ratio threshold.
Description
Technical field
What the present invention relates to is a kind of Underwater Detection and hydrolocation method, is specifically related to a kind of direction estimation method of acoustic vector sensor array.
Background technology
Sound field is with scalar field and vector field, and wherein acoustic pressure provides scalar information, and vibration velocity provides Vector Message, and this direction vector is consistent with sound direction of wave travel.Vector hydrophone is combined by pressure hydrophone and velocity hydrophone, can pick up acoustic pressure and vibration velocity information in sound field simultaneously, compared to traditional pressure hydrophone array, vector hydrophone arrays has obtained more information, thereby larger processing space can be provided.
Chinese scholars has been carried out a large amount of research to the orientation estimation application based on acoustic vector sensor array.Document [1], Malcolm Hawkes and Arye Nehorai (Acoustic Vector-sensor Corrections in Ambient Noise.IEEE J.Oceanic Eng, 2001,26 (3): 337-347.) correlation properties of the acoustic pressure in isotropic noise field and vibration velocity are inquired into, and provided the computing formula of acoustic pressure and vibration velocity related coefficient.Document [2] Bai Xingyu etc. (the long-range ESPRIT of a kind of new acoustic vector sensor array orientation algorithm for estimating. Harbin Engineering University's journal, 2006,27 (6): 891~895.), document [3] (acoustic vector sensor array information source number based on acoustic pressure vibration velocity Combined Treatment detects with orientation and estimates. acoustic journal, 2008,33 (1): 56~61.) propose the acoustic pressure vibration velocity Cross-covariance of acoustic vector sensor array to be applied in orientation estimation, the anti-isotropic noise ability of acoustic vector sensor array and high-resolution orientation algorithm for estimating are combined, obtained good effect.But the vibration velocity component of the method has spatial filtering effect, may cause interested signal conductively-closed to be fallen.Document [4] moon etc. (the non-space ESPRIT algorithm of vector array. Harbin Engineering University's journal, 2009,30 (4): 406~410.) utilize the acoustic pressure of acoustic vector sensor array and the rotational invariance between vibration velocity sensor array, can the in the situation that of unknown array manifold, to target, carry out orientation estimation, but the method is not utilized the uncorrelated characteristic of acoustic pressure and vibration velocity channel noise, estimated accuracy need to improve.
Summary of the invention
The object of the invention is to propose a kind ofly can be used in the even unknown formation acoustic vector of any formation array, and combine the anti-isotropic noise ability of acoustic vector sensor array, the acoustic vector sensor array High Precision Long-distance direction estimation method that precision is high.
The object of the present invention is achieved like this:
(1) the reception signal of two-dimensional vector hydrophone array is carried out to narrow-band filtering, obtain the arrowband output signal on pending frequency;
(2) two mutually orthogonal vibration velocity components of two-dimensional vector hydrophone array are carried out to linear combination in complex field, be converted to two new vibration velocity output components, construct the two vibration velocity components of complex field;
(3) adopt maximal possibility estimation to obtain acoustic pressure that some samplings take soon and the Cross-covariance pair of the two vibration velocity components of complex field, construct the acoustic pressure vibration velocity Cross-covariance pair based on broad sense acoustic energy flow;
(4) the ESPRIT algorithm of application matrix bundle calculate acoustic pressure vibration velocity Cross-covariance based on broad sense acoustic energy flow between the invariable rotary factor, thereby carry out orientation estimation.
The present invention can also comprise:
1. the two vibration velocity components of structure complex field described in specifically comprise: by the vibration velocity output component v of y direction
yj multiplies each other with complex unit, then with the vibration velocity output component v of x direction
xbe added and obtain new vibration velocity output component v
1; By the vibration velocity output component v of y direction
ymultiply each other with negative complex unit-j, then with the vibration velocity output component v of x direction
xbe added and obtain new vibration velocity output component v
2.V
1and v
2formed the two vibration velocity components of complex field, they have still carried the directional information of vibration velocity field, and have the receiving ability of omnidirectional, are a kind of improvement to conventional vibration velocity component.
2. described in, construct acoustic pressure vibration velocity Cross-covariance based on broad sense acoustic energy flow to specifically comprising: the output of the sound pressure sensor array of the snap quantity of necessarily sampling is multiplied each other with the conjugate transpose of two new complex field vibration velocity output components respectively, and ask its mean value, obtain Cross-covariance based on broad sense acoustic energy flow to R
1and R
2.
3. the ESPRIT algorithm based on broad sense acoustic energy flow described in, its thought is: solve the right non-zero generalized eigenvalue of acoustic pressure vibration velocity covariance matrix of above-mentioned gained, thereby provide the orientation estimated result of target.For fear of occurring potential numerical value difficulty in the process solving non-zero generalized eigenvalue, the ESPRIT algorithm of described application matrix bundle calculate acoustic pressure vibration velocity Cross-covariance based on broad sense acoustic energy flow between the invariable rotary factor in adopt svd (SVD) and total least square method (TLS) Ill-conditioned Least-squares Problems of a larger dimension to be converted into the non-ill total least square problem of a less dimension.
Advantage body of the present invention is mainly reflected in: the ESPRIT algorithm of (1) conventional acoustic vector sensor array adopts the autocorrelation matrix of data to carry out feature decomposition, even in the ideal case, still has noise in its signal subspace.For acoustic vector sensor array, when array element distance meets the noise correlation radius of acoustic pressure and vibration velocity, noise between acoustic pressure and vibration velocity component is incoherent, the present invention adopts the Cross-covariance of acoustic pressure and vibration velocity, take full advantage of the uncorrelated characteristic of acoustic pressure and vibration velocity noise field, there is good anti-isotropic noise ability; (2) conventional vibration velocity output component and the combination vibration velocity output based on Givess rotation all have space directivity, the incoming signal of different directions is had to different amplitude weightings, especially from the signal of the orthogonal directions incident of observed ray, can be masked completely.The present invention improves vibration velocity output component, and the incoming signal of all directions is all had to identical amplitude response, can completely receive from the signal of all directions incident; (3) between the vibration velocity output component that the present invention constructs, there is rotational invariance, its corresponding twiddle factor is only relevant with the incident direction of signal, and irrelevant with the locus of array, application ESPRIT algorithm can be realized in formation arbitrarily and even carry out the estimation of high-precision orientation unknown formation in the situation that.
Accompanying drawing explanation
Fig. 1 is the spatial correlation function curve of acoustic pressure vibration velocity in isotropic noise field;
Fig. 2 is the one group of irregular formation acoustic vector sensor array position view providing;
Fig. 3 is the block diagram of realizing of the present invention;
Fig. 4 (a)-Fig. 4 (e) is single goal signal two vibration velocity channel output spectrums of acoustic pressure, vibration velocity and complex field during from 30 ° of direction incident, wherein: Fig. 4 (a) sound pressure channel frequency spectrum, Fig. 4 (b) vibration velocity channel v
xfrequency spectrum, Fig. 4 (c) vibration velocity channel v
yfrequency spectrum, Fig. 4 (d) vibration velocity channel v
1frequency spectrum, Fig. 4 (e) vibration velocity channel v
2frequency spectrum;
Fig. 5 (a)-Fig. 5 (e) is that Bi-objective signal is respectively from 30 °, two vibration velocity channel output spectrums of acoustic pressure, vibration velocity and complex field during 40 ° of direction incident, wherein: Fig. 5 (a) sound pressure channel frequency spectrum, Fig. 5 (b) vibration velocity channel v
xfrequency spectrum, Fig. 5 (c) vibration velocity channel v
yfrequency spectrum, Fig. 5 (d) vibration velocity channel v
1frequency spectrum, Fig. 5 (e) vibration velocity channel v
2frequency spectrum;
Fig. 6 is that root-mean-square error curve is estimated in single goal orientation;
Fig. 7 is that root-mean-square error curve is estimated in Bi-objective orientation;
Fig. 8 (a)-Fig. 8 (e) is signal two vibration velocity channel output spectrums of acoustic pressure, vibration velocity and complex field during from the incident of y direction of principal axis, wherein: Fig. 8 (a) sound pressure channel frequency spectrum, Fig. 8 (b) vibration velocity channel v
xfrequency spectrum, Fig. 8 (c) vibration velocity channel v
yfrequency spectrum, Fig. 8 (d) vibration velocity channel v
1frequency spectrum, Fig. 8 (e) vibration velocity channel v
2frequency spectrum;
Fig. 9 (a)-Fig. 9 (b) is the orientation estimated result near the signal of direction incident y axle, Fig. 9 (a) is the orientation estimated result that adopts conventional vibration velocity output component, and Fig. 9 (b) is the orientation estimated result that adopts the two vibration velocity components of complex field.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is further described in detail:
1. calculation procedure
1.1. array output model
The polar coordinates of supposing i element position of the two-dimentional acoustic vector sensor array of N unit are [r
i, α
i], i=1,2 ..., N, far field, the arrowband plane wave that D centre frequency is ω incides on basic matrix, and wherein the incident direction of k signal is θ
k, the output signal of array can be designated as:
P in formula (t), v
x(t), v
y(t) be respectively the N * 1 dimension reception snap vector of sound pressure channel and x, y vibration velocity channel; S (t) is that D far-field signal source is at D * 1 of reference position dimension acoustic pressure vector; n
p(t), n
x(t), n
y(t) be the local noise vector of N * 1 dimension that each passage is corresponding.
A
p=[a
p1..., a
pD] be the sound pressure channel stream shape matrix that signal vector is corresponding, the stream shape vector of k signal can be written as:
a
pk=[exp(jωr
Ncos(θ
k-α
1)/c),…,exp(jωr
Ncos(θ
k-α
N)/c)]
T (2)
Operational symbol ()
trepresenting matrix transposition, A
x=[a
x1..., a
xD] be the x direction vibration velocity channel stream shape matrix that signal vector is corresponding, the stream shape vector of k signal can be written as:
a
xk=[exp(jωr
Ncos(θ
k-α
1)/c),…,exp(jωr
Ncos(θ
k-α
N)/c)]
Tcosθ
k (3)
A
y=[a
y1..., a
yD] be the y direction vibration velocity channel stream shape matrix that signal vector is corresponding, the stream shape vector of k signal can be written as:
a
yk=[exp(jωr
Ncos(θ
k-α
1)/c),…,exp(jωr
Ncos(θ
k-α
N)/c)]
Tsinθ
k (4)
From formula (2)-(4), A
p, A
x, A
ybetween there is following relation:
A
x=A
pΦ
x,A
y=A
pΦ
y,A
y=A
xΦ
xy (5)
In formula, Φ
x, Φ
y, Φ
xyare all diagonal matrix, and have:
From formula (5), can find out, for far field plane wave incoming signal, between the sound pressure sensor array of acoustic vector sensor array and each vibration velocity component sensor array, all there is invariable rotary relation, and this twiddle factor is only relevant with signal incident direction, and irrelevant with the locus of array, this is that acoustic vector sensor array can even be applied the physical basis that ESPRIT algorithm carries out orientation estimation unknown formation in the situation that in formation arbitrarily.
1.2. two vibration velocity output components of new complex field
In isotropic noise field, the spatial correlation function curve of acoustic pressure vibration velocity as shown in Figure 1.1. curve has shown the autocorrelation function curve of acoustic pressure, at d=0 (d is space length), locate correlativity the strongest, and be 0 at the auto-correlation function value that d=0.5 λ (λ is the wavelength of signal) locates acoustic pressure, therefore the algorithm based on acoustic pressure battle array adopts half-wavelength spacing conventionally.2. curve has shown the spatial cross correlation function of acoustic pressure and vibration velocity, and from figure, the acoustic pressure of space concurrent and vibration velocity are incoherent, and at d=0.715 λ place, its correlation function value is also 0, and the noise correlation radius of known acoustic pressure and vibration velocity is about 0.715 λ.Different from traditional acoustic pressure battle array, when if consider, the array element distance of acoustic vector sensor array meets the noise correlation radius of acoustic pressure and vibration velocity, now have:
In formula, operational symbol ()
hthe conjugate transpose of representing matrix.From formula (7), adopt the Cross-covariance of acoustic pressure and vibration velocity, will there is the ability of good anti-isotropic noise.It should be noted that the dipole directive property due to vibration velocity sensor self, conventional acoustic pressure vibration velocity Cross-covariance has the effect of spatial filtering, likely can cause interested parties to signal weakenedly even masked completely.The Cross-covariance of vibration velocity component of acoustic pressure and x direction of take is example, from the signal of y direction of principal axis incident, will be masked completely.In order to address this problem, the present invention proposes, in complex field, vibration velocity component is carried out to linear combination, the two vibration velocity output component v of complex field that structure makes new advances
1, v
2:
In formula, n
1, n
2it is corresponding ground unrest vector; Φ
x, Φ
yare all diagonal matrix, and have:
1.3. the acoustic pressure vibration velocity Cross-covariance based on broad sense acoustic energy flow
The two vibration velocity output components of new complex field are all identical to the amplitude response of the incoming signal of all directions, and its vector feature is only embodied in phase place, therefore can completely receive the signal of in any direction incident.Another feature is that the new two ground unrests of vibration velocity output component of complex field and the ground unrest of sound pressure channel is still incoherent, can find out from following formula:
Consider that acoustic pressure is as follows with the covariance matrix of new complex field pair vibration velocity components:
In formula: R
s=E (ss
h) be the incoming signal covariance matrix of D * D dimension.When incoming signal mutual when uncorrelated, R
sa diagonal matrix, the power that diagonal element is each signal.Because the two vibration velocity output components of new complex field are the expansions to former acoustic vector sensor array vibration velocity component, therefore new acoustic pressure vibration velocity Cross-covariance is called to the Cross-covariance based on broad sense acoustic energy flow.
R in practical application
1, R
2employing maximum likelihood estimator replaces:
In formula, M is the fast umber of beats of sampling.
1.4. the orientation based on pencil of matrix ESPRIT algorithm is estimated
Making γ is a constant, examination pencil of matrix (R
1, R
2):
Due to A
prow full ranks, and R
snonsingular, rank of matrix is asked in (13) formula both sides:
When
time, matrix
be nonsingular matrix, and work as
time, matrix
order wanes becomes singular matrix.
all pencil of matrix (R
1, R
2) non-zero generalized eigenvalue.
Generally, covariance matrix R
1, R
2be full rank not, in solving the process of non-zero generalized eigenvalue, likely cause potential numerical value difficulty.The present invention adopts svd (SVD) and total least square method (TLS), can avoid this problem.
To R
1svd:
As mentioned before, Σ in formula
1by R
1the diagonal matrix that forms of non-zero singular value; U
1, V
1respectively left and right singular vector corresponding to non-zero singular value; U
2, V
2zero left and right singular vector corresponding to singular value.So can obtain:
Visible, N * N dimension pencil of matrix (R
1, R
2) D non-zero generalized eigenvalue be that D * D ties up pencil of matrix
whole generalized eigenvalues.
This direction estimation method of the present invention, as shown in Figure 3, concrete calculation procedure is summarized as follows implementation method:
(1) the reception signal of pair array carries out narrow-band filtering, obtains the arrowband output signal on pending frequency.
(2) according to formula (8), calculate the two vibration velocity output components of new complex field.
(3) according to formula (12), estimate the acoustic pressure vibration velocity Cross-covariance based on broad sense acoustic energy flow.
(4) covariance matrix is carried out to svd, and estimate non-zero singular value matrix Σ
1left and right singular vector matrix U with correspondence
1, V
1.
(5) to pencil of matrix
carry out generalized eigenvalue decomposition, thereby obtain the azimuth estimation value of incoming signal.
2. the concrete calculated examples that orientation is estimated
As shown in Figure 2, the centre frequency of arrowband, incident far field plane wave is 1kHz to 5 yuan of acoustic vector sensor array element positions.
2.1. the root-mean-square error simulation example of target Bearing Estimation
Adopt the method in direction estimation method proposed by the invention and document [4] to carry out emulation experiment 500 times simultaneously, having carried out respectively single goal and Bi-objective orientation estimates, the fast umber of beats of sampling of each experiment is 100, and in single goal situation, target level position angle is 30 °; In Bi-objective situation, target level position angle is respectively 30 °, 40 °.
First, the reception data of pair array are carried out narrow-band filtering, and calculate the two vibration velocity output component v of new complex field according to formula (8)
1, v
2, when Fig. 4, Fig. 5 have shown respectively single goal and Bi-objective situation, the output spectrum of each passage when signal to noise ratio (S/N ratio) (SNR) is 0dB.As seen from the figure, due to the filtering characteristic of self, the spectrum peak of conventional vibration velocity output component is less than sound pressure channel; And the two vibration velocity components of complex field have the receiving ability of omnidirectional, its spectrum peak and sound pressure channel are basically identical.
Then, according to formula (12), estimate covariance matrix based on broad sense acoustic energy flow to R
1, R
2, last application matrix bundle ESPRIT method calculates matrix to R
1, R
2between the invariable rotary factor, obtain azimuth estimation value.Fig. 6, Fig. 7 have provided respectively root mean square (RMS) error of target level position angle estimation in two kinds of situations and the relation curve of SNR.In figure, horizontal ordinate is SNR, and the ordinate of Fig. 6 is the RMS error that single goal is estimated, the ordinate of Fig. 7 is the RMS error sum that Bi-objective is estimated.
By Fig. 6 and Fig. 7, can be found out, under logarithmic coordinate, the root-mean-square error of two kinds of methods of estimation and signal to noise ratio (S/N ratio) are roughly linear, and method proposed by the invention has lower evaluated error all the time compared to the method for document [4].In the single goal situation shown in Fig. 6, when SNR=5dB, the method root-mean-square error of document [4] is down to 1 ° of left and right, and method proposed by the invention is when SNR=0dB, is just down to 1 °; In the Bi-objective situation shown in Fig. 7, when SNR=17dB, the method root-mean-square error of document [4] is down to 1 ° of left and right, and method proposed by the invention is when SNR=10dB, just can be down to 1 °.Visible, method proposed by the invention effectively combines the ability of the anti-isotropic noise of vector hydrophone, thereby has the lower signal-noise ratio threshold processed.
2.2.y near axle, simulation example is estimated in the orientation of direction incoming signal
Consider the situation of single goal incident, the incident direction of signal changes to 93 ° with the step-length of 0.5 ° from 87 °, signal to noise ratio snr=0dB, and the fast umber of beats of sampling is 100.Adopt respectively the two vibration velocity output component structure of conventional vibration velocity output component and complex field acoustic pressure vibration velocity Cross-covariance, and carry out orientation estimation.When Fig. 8 has shown target from the incident of y direction of principal axis, the output signal spectrum of each passage, Fig. 9 has shown the orientation estimated result of this two schemes.
By Fig. 8 and Fig. 9, can be found out, when signal is during from the incident of y direction of principal axis, vibration velocity channel v
xthoroughly do not receive signal, therefore adopt the ESPRIT orientation estimated result of acoustic pressure and conventional vibration velocity component to lose efficacy, and the two vibration velocity channel v of complex field
1, v
2the complete signal that receives of energy, still can provide high-precision orientation estimated result.
Claims (5)
1. an acoustic vector sensor array High Precision Long-distance direction estimation method, is characterized in that:
(1) the reception signal of two-dimensional vector hydrophone array is carried out to narrow-band filtering, obtain the arrowband output signal on pending frequency;
(2) two mutually orthogonal vibration velocity components of two-dimensional vector hydrophone array are carried out to linear combination in complex field, be converted to two new vibration velocity output components, construct the two vibration velocity components of complex field;
(3) adopt maximal possibility estimation to obtain acoustic pressure that some samplings take soon and the Cross-covariance pair of the two vibration velocity components of complex field, construct the acoustic pressure vibration velocity Cross-covariance pair based on broad sense acoustic energy flow;
(4) the ESPRIT algorithm of application matrix bundle calculate acoustic pressure vibration velocity Cross-covariance based on broad sense acoustic energy flow between the invariable rotary factor, thereby carry out orientation estimation.
2. acoustic vector sensor array High Precision Long-distance direction estimation method according to claim 1, is characterized in that the two vibration velocity components of described structure complex field specifically comprise: by the vibration velocity output component v of y direction
yj multiplies each other with complex unit, then with the vibration velocity output component v of x direction
xbe added and obtain new vibration velocity output component v
1; By the vibration velocity output component v of y direction
ymultiply each other with negative complex unit-j, then with the vibration velocity output component v of x direction
xbe added and obtain new vibration velocity output component v
2.
3. acoustic vector sensor array High Precision Long-distance direction estimation method according to claim 1, it is characterized in that described constructing acoustic pressure vibration velocity Cross-covariance based on broad sense acoustic energy flow to specifically comprising: the output of the sound pressure sensor array of the snap quantity of necessarily sampling is multiplied each other with the conjugate transpose of two new complex field vibration velocity output components respectively, and ask its mean value, obtain Cross-covariance based on broad sense acoustic energy flow to R
1and R
2.
4. acoustic vector sensor array High Precision Long-distance direction estimation method according to claim 1 and 2, is characterized in that: the ESPRIT algorithm of described application matrix bundle calculate acoustic pressure vibration velocity Cross-covariance based on broad sense acoustic energy flow between the invariable rotary factor in adopt svd and total least square method the Ill-conditioned Least-squares Problems of a larger dimension to be converted into the non-ill total least square problem of a less dimension.
5. acoustic vector sensor array High Precision Long-distance direction estimation method according to claim 3, is characterized in that: the ESPRIT algorithm of described application matrix bundle calculate acoustic pressure vibration velocity Cross-covariance based on broad sense acoustic energy flow between the invariable rotary factor in adopt svd and total least square method the Ill-conditioned Least-squares Problems of a larger dimension to be converted into the non-ill total least square problem of a less dimension.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310322324.2A CN103605108B (en) | 2013-07-29 | 2013-07-29 | High-precision remote direction estimation method of acoustic vector array |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310322324.2A CN103605108B (en) | 2013-07-29 | 2013-07-29 | High-precision remote direction estimation method of acoustic vector array |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103605108A true CN103605108A (en) | 2014-02-26 |
CN103605108B CN103605108B (en) | 2015-09-30 |
Family
ID=50123348
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310322324.2A Active CN103605108B (en) | 2013-07-29 | 2013-07-29 | High-precision remote direction estimation method of acoustic vector array |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103605108B (en) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104180824A (en) * | 2014-08-18 | 2014-12-03 | 中国科学院上海应用物理研究所 | Method for improving measurement accuracy of probe based on principal component analysis algorithm |
CN104391268A (en) * | 2014-11-21 | 2015-03-04 | 哈尔滨工业大学深圳研究生院 | Direction-of-arrival angle detection method |
CN104392114A (en) * | 2014-11-11 | 2015-03-04 | 西北大学 | High-resolution target direction estimation method based on space-time data |
CN105223544A (en) * | 2015-08-26 | 2016-01-06 | 南京信息工程大学 | The constant Beamforming Method of the near field linear constraint adaptive weighted frequency of minimum variance |
CN105676168A (en) * | 2015-12-02 | 2016-06-15 | 江苏科技大学 | Acoustic vector array DOA estimation method |
CN105954709A (en) * | 2016-05-06 | 2016-09-21 | 哈尔滨工程大学 | Acoustic vector circular array source number detection method based on characteristic value multiple threshold correction |
CN106997038A (en) * | 2017-03-21 | 2017-08-01 | 陕西理工大学 | Any acoustic vector-sensor array row near field sources ESPRIT method for parameter estimation |
CN107728109A (en) * | 2017-09-18 | 2018-02-23 | 哈尔滨工程大学 | A kind of noncooperative target radiated noise measurement and positioning technology |
CN109443516A (en) * | 2018-12-25 | 2019-03-08 | 西北工业大学 | A kind of passive acquisition methods of Bottom sound speed based on the vertical vibration velocity signal of noise field |
CN109489799A (en) * | 2018-12-25 | 2019-03-19 | 西北工业大学 | A kind of Bottom sound speed substep inversion method based on double vector hydrophones |
CN109696657A (en) * | 2018-06-06 | 2019-04-30 | 南京信息工程大学 | A kind of coherent sound sources localization method based on vector hydrophone |
CN110082712A (en) * | 2019-03-14 | 2019-08-02 | 哈尔滨工程大学 | A kind of acoustic vector circle battle array Coherent Targets direction estimation method |
CN110824428A (en) * | 2019-11-06 | 2020-02-21 | 哈尔滨工程大学 | Vertical vector array underwater sound ray matching passive positioning method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4750147A (en) * | 1985-11-06 | 1988-06-07 | Stanford University | Method for estimating signal source locations and signal parameters using an array of signal sensor pairs |
US6545639B1 (en) * | 2001-10-09 | 2003-04-08 | Lockheed Martin Corporation | System and method for processing correlated contacts |
CN102226837A (en) * | 2011-04-08 | 2011-10-26 | 哈尔滨工程大学 | Vector circle array acoustic pressure and vibration velocity combined direction finding method on cylindrical form baffle condition |
CN102879764A (en) * | 2012-10-16 | 2013-01-16 | 浙江大学 | Underwater sound source direction estimating method |
-
2013
- 2013-07-29 CN CN201310322324.2A patent/CN103605108B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4750147A (en) * | 1985-11-06 | 1988-06-07 | Stanford University | Method for estimating signal source locations and signal parameters using an array of signal sensor pairs |
US6545639B1 (en) * | 2001-10-09 | 2003-04-08 | Lockheed Martin Corporation | System and method for processing correlated contacts |
CN102226837A (en) * | 2011-04-08 | 2011-10-26 | 哈尔滨工程大学 | Vector circle array acoustic pressure and vibration velocity combined direction finding method on cylindrical form baffle condition |
CN102879764A (en) * | 2012-10-16 | 2013-01-16 | 浙江大学 | Underwater sound source direction estimating method |
Non-Patent Citations (1)
Title |
---|
姚直象等: "基于声压振速联合处理的矢量阵旋转不变子空间方位估计方法", 《北京理工大学学报》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104180824A (en) * | 2014-08-18 | 2014-12-03 | 中国科学院上海应用物理研究所 | Method for improving measurement accuracy of probe based on principal component analysis algorithm |
CN104392114A (en) * | 2014-11-11 | 2015-03-04 | 西北大学 | High-resolution target direction estimation method based on space-time data |
CN104392114B (en) * | 2014-11-11 | 2018-04-17 | 西北大学 | A kind of high resolution target direction estimation method based on spatial-temporal data |
CN104391268B (en) * | 2014-11-21 | 2017-10-13 | 哈尔滨工业大学深圳研究生院 | A kind of ripple reaches azimuthal detection method |
CN104391268A (en) * | 2014-11-21 | 2015-03-04 | 哈尔滨工业大学深圳研究生院 | Direction-of-arrival angle detection method |
CN105223544A (en) * | 2015-08-26 | 2016-01-06 | 南京信息工程大学 | The constant Beamforming Method of the near field linear constraint adaptive weighted frequency of minimum variance |
CN105223544B (en) * | 2015-08-26 | 2018-01-12 | 南京信息工程大学 | Near field linear constrains the constant Beamforming Method of the adaptive weighted frequency of minimum variance |
CN105676168A (en) * | 2015-12-02 | 2016-06-15 | 江苏科技大学 | Acoustic vector array DOA estimation method |
CN105676168B (en) * | 2015-12-02 | 2018-03-06 | 江苏科技大学 | A kind of acoustic vector sensor array direction estimation method |
CN105954709A (en) * | 2016-05-06 | 2016-09-21 | 哈尔滨工程大学 | Acoustic vector circular array source number detection method based on characteristic value multiple threshold correction |
CN105954709B (en) * | 2016-05-06 | 2018-12-07 | 哈尔滨工程大学 | One kind is based on the modified acoustic vector circle battle array information source number detection method of characteristic value multi-threshold |
CN106997038A (en) * | 2017-03-21 | 2017-08-01 | 陕西理工大学 | Any acoustic vector-sensor array row near field sources ESPRIT method for parameter estimation |
CN106997038B (en) * | 2017-03-21 | 2019-10-18 | 陕西理工大学 | Any acoustic vector-sensor array column near field sources ESPRIT method for parameter estimation |
CN107728109A (en) * | 2017-09-18 | 2018-02-23 | 哈尔滨工程大学 | A kind of noncooperative target radiated noise measurement and positioning technology |
CN109696657A (en) * | 2018-06-06 | 2019-04-30 | 南京信息工程大学 | A kind of coherent sound sources localization method based on vector hydrophone |
CN109443516A (en) * | 2018-12-25 | 2019-03-08 | 西北工业大学 | A kind of passive acquisition methods of Bottom sound speed based on the vertical vibration velocity signal of noise field |
CN109489799A (en) * | 2018-12-25 | 2019-03-19 | 西北工业大学 | A kind of Bottom sound speed substep inversion method based on double vector hydrophones |
CN110082712A (en) * | 2019-03-14 | 2019-08-02 | 哈尔滨工程大学 | A kind of acoustic vector circle battle array Coherent Targets direction estimation method |
CN110824428A (en) * | 2019-11-06 | 2020-02-21 | 哈尔滨工程大学 | Vertical vector array underwater sound ray matching passive positioning method |
Also Published As
Publication number | Publication date |
---|---|
CN103605108B (en) | 2015-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103605108B (en) | High-precision remote direction estimation method of acoustic vector array | |
CN103901395B (en) | Coherent signal direction of arrival Dynamic Tracking under a kind of impulsive noise environment | |
CN103353596B (en) | Wave beam space domain meter wave radar height measurement method based on compressed sensing | |
CN102279387B (en) | Method for estimating target arrival angle of multiple input multiple output (MIMO) radar | |
CN106526530B (en) | 2-L type array arrival direction estimation algorithm based on propagation operator | |
CN105589056B (en) | A kind of multiple target distance field mixing source localization method | |
CN103353595B (en) | Meter wave radar height measurement method based on array interpolation compression perception | |
CN104155648B (en) | The high-frequency ground wave radar single snap MUSIC direction-finding method reset based on array data | |
CN103926572B (en) | A kind of clutter suppression method of battle array radar self adaptation subspace, airborne anon-normal side | |
CN109471082A (en) | Array element defect MIMO radar angle estimating method based on signal subspace reconstruct | |
CN106802403B (en) | Acoustic vector sensors two-dimensional array MUSIC decorrelation LMS method for parameter estimation | |
CN103744076B (en) | MIMO radar moving target detection method based on non-convex optimization | |
CN103344940B (en) | The DOA estimation method of low complex degree and system | |
CN104407335B (en) | DOA estimation method of 3-axis cross array | |
CN104515969B (en) | Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method | |
CN103760527B (en) | Single base MIMO radar coherent source Wave arrival direction estimating method | |
CN106353738B (en) | A kind of robust adaptive beamforming method under new DOA mismatch condition | |
CN103091661A (en) | Broadband signal arriving direction estimation method based on iteration spectral reconfiguration | |
CN106019234A (en) | L-shaped antenna array low computation complexity two-dimensional DOA estimation method | |
CN102662158B (en) | Quick processing method for sensor antenna array received signals | |
CN102393525A (en) | Navigation interference suppression and signal amplification method for subspace projection | |
CN104035069B (en) | Arrowband based on partial correction linear array symmetrically and evenly near-field signals source location method | |
CN104020440B (en) | Interfere the two-dimentional direction of arrival estimation method of formula linear array based on L-type | |
CN103116162B (en) | High-resolution sonar location method based on sparsity of objective space | |
CN103364772A (en) | Target low elevation estimation method based on real number field generalized multiple-signal sorting algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |