CN112130112A - Information source number estimation method based on acoustic vector array joint information processing - Google Patents
Information source number estimation method based on acoustic vector array joint information processing Download PDFInfo
- Publication number
- CN112130112A CN112130112A CN202010991187.1A CN202010991187A CN112130112A CN 112130112 A CN112130112 A CN 112130112A CN 202010991187 A CN202010991187 A CN 202010991187A CN 112130112 A CN112130112 A CN 112130112A
- Authority
- CN
- China
- Prior art keywords
- channel
- array
- vibration velocity
- acoustic
- acoustic vector
- 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
- G01S3/8027—By vectorial composition of signals received by plural, differently-oriented transducers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
The invention belongs to the technical field of array signal processing, and particularly relates to an information source number estimation method based on acoustic vector array joint information processing. The method is based on acoustic vector array sound pressure and vibration velocity combined information processing, and is suitable for the problem that the information source number estimation algorithm is inaccurate under the condition of low signal to noise ratio in the underwater acoustic environment. The invention adopts acoustic vector array sound pressure and vibration velocity combined information processing to estimate the number of information sources, and aims at solving the problem that signals of certain directions are possibly filtered or weakened due to the space filtering action when the observation direction is fixed. The invention utilizes the advantage of acoustic vector array azimuth scanning, avoids the weakening of signals in a certain direction due to spatial filtering, and improves the performance of the traditional information source number estimation algorithm.
Description
Technical Field
The invention belongs to the technical field of array signal processing, and particularly relates to an information source number estimation method based on acoustic vector array joint information processing.
Background
Underwater acoustic array signal processing is an important technology for underwater detection, identification, communication and the like, and one of important tasks is to estimate the target orientation. The subspace type high-resolution spectrum estimation algorithm represented by the multiple signal classification algorithm is proposed successively, and the subspace type high-resolution spectrum estimation algorithm has good prospects in the field of underwater acoustic engineering for estimating the direction of arrival. The source number estimation is a prerequisite for the high-resolution subspace-like spectrum estimation algorithm to perform the azimuth estimation, so that the research of the source number estimation is one of the important research directions in the array signal processing.
Based on radar and sonar scalar arrays, the method is mainly divided into two categories of algorithms based on characteristic value information and characteristic vector information, the performance of the source number estimation algorithm is poor under low signal-to-noise ratio, and in underwater acoustic environment, the problem to be solved is how to accurately estimate the source number under low signal-to-noise ratio.
In the article, the acoustic vector array source number detection and orientation estimation (Chinese edition) based on the combined processing of sound pressure and vibration velocity is disclosed in 'Baixingyu, Jiangyu, Zhao Chunhui', 2008(01):56-61. The method adopts the method of acoustic vector array PVc combined information processing, cannot well inhibit noise, has performance decline under the condition of low signal-to-noise ratio, has narrow application range, and cannot adopt acoustic vector array sound pressure and vibration velocity combined processing to estimate the information source number with various criteria.
The Yao straight image is a vector array rotation invariant subspace orientation estimation method [ J ] based on sound pressure and vibration velocity combined processing, university of Beijing science and technology, 2012,32(05):513 +521 ], wherein a sound vector array is utilized to carry out cover-type circle criterion information source number estimation by adopting sound pressure and vibration velocity combined processing.
In summary, accurate estimation of underwater information source number under low signal-to-noise ratio is lacked in the current stage, so that the method is popularized to an acoustic vector array sound pressure and vibration velocity combined information processing feature projection method for information source number estimation algorithm.
Disclosure of Invention
The invention aims to provide an information source number estimation method based on acoustic vector array joint information processing, which can accurately estimate the underwater information source number under the underwater acoustic low signal-to-noise ratio.
The purpose of the invention is realized by the following technical scheme: the method comprises the following steps:
step 1: receiving far-field plane wave data of a plurality of information sources by utilizing a two-dimensional acoustic vector hydrophone uniform linear array to obtain an output component p of a sound pressure channel of an array element mm(t), output component v of vibration velocity x channelxm(t), output component v of the vibration velocity y channelym(t):
Wherein M is the number of acoustic vector arrays of the vector hydrophones; k is the number of space signals of incident sound to the vector array; thetakFor the k-th space signal sk(t) a two-dimensional spatial angle of arrival; n ispm(t),Respectively corresponding to the noise of the m-th array element sound pressure channel, the vibration velocity x channel and the vibration velocity y channel; a (theta)k) The vector is a guide vector of a uniform linear array;
step 2: arranging and combining the received data of each array element to obtain the output P (t) of the sound pressure channel of the acoustic vector hydrophone array and the output V of the vibration velocity x channelx(t) and output V of the Y channely(t);
And step 3: the vibration velocity x channel data Vx(t) and velocity y channel data Vy(t) combining them to obtain data V of combined channelc(t);
and 4, step 4: guiding azimuth angles to scan according to uniform intervals within an acoustic vector array observation angle range;
and 5: constructing a covariance matrix R (in a certain range) of sound pressure and vibration velocity joint processing under each azimuth guiding azimuth anglep+v)v;
Step 6: the covariance matrix constructed for each azimuth angle is arranged from large to small according to the row and norm thereof;
and 7: the covariance matrixes corresponding to the rows and the first 25 percent of the norm are added to obtain a new covariance matrix Rpv;
And 8: r is to bepvProjecting the image into an eigenvector matrix U to obtain a projection P ═ UHRpv(ii) a Wherein the projection P ═ UHRpvIs P in the ith column (i ═ 1, 2.., M)i=[p1i,p2i,...,pmi,...,pMi]T;
And step 9: calculate and makeThe m value with the maximum value is used as the estimated value of the number of the information sourcesWherein the content of the first and second substances,being the average of all the element modes of the m-th row of the projection P,
the invention has the beneficial effects that:
the invention provides the acoustic vector array joint information processing-based information source number estimation method which can accurately estimate the underwater information source number under the underwater acoustic low signal-to-noise ratio and popularize the information source number estimation algorithm to the acoustic vector array sound pressure and vibration velocity joint information processing. The method is based on acoustic vector array sound pressure and vibration velocity combined information processing, and is suitable for the problem that the information source number estimation algorithm is inaccurate under the condition of low signal to noise ratio in the underwater acoustic environment. The invention adopts acoustic vector array sound pressure and vibration velocity combined information processing to estimate the number of information sources, and aims at solving the problem that signals of certain directions are possibly filtered or weakened due to the space filtering action when the observation direction is fixed.
Drawings
FIG. 1 is a flow chart of the present invention.
Fig. 2 is a schematic diagram of the received signals of a uniform linear array of acoustic vector hydrophones.
FIG. 3 is a diagram of the accuracy of source estimation under the eigenspace Projection (ESP) source estimation and the SNR relationship between the method of the present invention and the prior art method.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
The invention provides the acoustic vector array joint information processing-based information source number estimation method which can accurately estimate the underwater information source number under the underwater acoustic low signal-to-noise ratio and popularize the information source number estimation algorithm to the acoustic vector array sound pressure and vibration velocity joint information processing.
Step 1: firstly, a two-dimensional acoustic vector hydrophone uniform linear array is utilized to receive far-field plane wave data of a plurality of information sources, an acoustic vector array of M vector hydrophones is utilized, K space signals are incident on the vector array, and the kth space signal sk(t) has a two-dimensional angle of arrival of θkTo obtain the output component p of the acoustic pressure channel of the array element mm(t), output component v of vibration velocity x channelxm(t), output component v of the vibration velocity y channelym(t):
npm(t),Respectively corresponding to the m-th array element sound pressure channel, the vibration velocity x channel and the vibration velocity y channel. a (theta)i) Steering vector for uniform linear array:
step 2: arranging and combining the received data of each array element to obtain an acoustic vector hydrophone array sound pressure channel P (t) and a vibration velocity x channel Vx(t) vibration velocity y channel VyOutput of (t):
and step 3: the vibration velocity x channel data Vx(t) and velocity y channel data Vy(t) combining with each other to obtain data of a combined channel:whereinIs the leading azimuth;
and 4, step 4: guiding azimuth angles to scan according to uniform intervals within an acoustic vector array observation angle range;
and 5: constructing a covariance matrix of sound pressure and vibration velocity combined processing under each azimuth guiding azimuth angle:
step 6: the covariance matrix constructed for each azimuth angle is arranged from large to small according to the row and norm thereof;
and 7: the covariance matrixes corresponding to the rows and the first 25 percent of the norm are added to obtain a new covariance matrix Rpv。
And 8: r is to bepvProjecting the image into an eigenvector matrix U to obtain a projection P ═ UHRpvThe ith column (i ═ 1, 2.., M) indicates: pi=[p1i,p2i,...,pmi,...,pMi]TThe source number estimation is performed by processing the projection matrix P. The mean of all the element modes in the m-th row of P is recorded asHas a value ofpmL can be used as criterion for estimating the number of information sources, pmI is divided into two groups, the larger group corresponds to the incident sound signals, and the value of m in the larger group is the number of the information sources;
and step 9: the feature Space Projection (ESP) algorithm performs discrimination by using the following relational expressionSource number estimationValue ofTaking the value of m which maximizes the value of the formula;
compared with the existing information source number estimation method, the method has the following beneficial effects:
firstly, the method is based on acoustic vector array sound pressure and vibration velocity combined information processing, and is suitable for the problem that the information source number estimation algorithm is inaccurate under the condition of low signal to noise ratio in the underwater acoustic environment.
Secondly, the invention adopts acoustic vector array sound pressure and vibration velocity combined information processing to estimate the number of the information sources, and aims at solving the problem that signals of certain directions are possibly filtered or weakened due to the space filtering action when the observation direction is fixed.
Example 1:
the method provided by the invention and the existing method are adopted for simulation calculation, the simulation conditions are calculated by using MATLAB on a computer, and the specific parameters of the simulation environment are as follows:
array properties: a two-dimensional acoustic vector hydrophone array;
the number of array elements: 8;
array element spacing: half of the incident wavelength of the signal;
number of signal sources: 3, the number of the cells is 3;
signal source properties: mutually independent signals;
the incident angle of the signal source is as follows: (15 DEG, 30 DEG, 60 DEG)
Fast beat number: 1000, parts by weight;
signal-to-noise ratio: -20-10 dB;
experimental contents and results:
experiment 1: the method adopts the characteristic space projection to carry out the information source number estimation, and the performance is obviously improved
(1) Two mutually independent equal-power narrow-band signal sources are arranged and respectively marked as a signal source 1 and a signal source 2, the center frequency of the two independent equal-power narrow-band signal sources is 150Hz, and the incident angles are respectively 15 degrees, 30 degrees and 60 degrees; the array is an 8-element acoustic vector uniform linear array, the distance between array elements is half wavelength, the sampling frequency of the array is 10kHz, and the receiving bandwidth is 20 Hz-1500 Hz; the observation direction of the acoustic vector hydrophone array is set to scan between the range of minus 60 degrees and 60 degrees. The background noise is band-limited white Gaussian noise, the range of the signal-to-noise ratio is-20 dB to 10dB, the variation step length is 1dB, the snapshot number is 1000, and 500 times of Monte Carlo simulation is carried out on each signal-to-noise ratio.
Fig. 3 shows the relationship between the accuracy of the source number estimation using eigenspace Projection (ESP) and the signal-to-noise ratio in the method of the present invention and the method of the prior art.
The prior direction adopts a method proposed in the article of Yao Zheng Xiang, Jiangcao, Guohui and Chenjiang, a vector array rotation invariant subspace direction estimation method based on sound pressure and vibration velocity combined processing [ J ], university of Beijing university of science and technology, 2012,32(05): 513-.
Referring to fig. 3, the method provided by the present invention and the method in the prior art adopt the relationship between the accuracy of the characteristic Space Projection (ESP) information source estimation and the signal-to-noise ratio, and it can be seen from the figure that under the same signal-to-noise ratio, the success probability of the information source detection of the method of the present invention is far higher than that of the other two methods, and when the signal-to-noise ratio is-13 dB, the success probability of the detection of the method of the present invention completely approaches 1, and the method of the present invention has excellent performance at low signal-to-noise ratio. The method provided by the invention is obviously superior to the existing method and the method adopting sound pressure information processing, and the success probability of information source number estimation is obviously improved.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (1)
1. A method for estimating the information source number based on acoustic vector array joint information processing is characterized by comprising the following steps:
step 1: receiving multiple information sources by using two-dimensional acoustic vector hydrophone uniform linear arrayObtaining the output component p of the acoustic pressure channel of the array element m from the far-field plane wave datam(t), output component v of vibration velocity x channelxm(t), output component v of the vibration velocity y channelym(t):
Wherein M is the number of acoustic vector arrays of the vector hydrophones; k is the number of space signals of incident sound to the vector array; thetakFor the k-th space signal sk(t) a two-dimensional spatial angle of arrival;respectively corresponding to the noise of the m-th array element sound pressure channel, the vibration velocity x channel and the vibration velocity y channel; a (theta)k) The vector is a guide vector of a uniform linear array;
step 2: arranging and combining the received data of each array element to obtain the output P (t) of the sound pressure channel of the acoustic vector hydrophone array and the output V of the vibration velocity x channelx(t) and output V of the Y channely(t);
And step 3: the vibration velocity x channel data Vx(t) and velocity y channel data Vy(t) combining them to obtain data V of combined channelc(t);
and 4, step 4: guiding azimuth angles to scan according to uniform intervals within an acoustic vector array observation angle range;
and 5: constructing a covariance matrix R of sound pressure and vibration velocity combined processing under each azimuth guiding azimuth angle(p+v)v;
Step 6: the covariance matrix constructed for each azimuth angle is arranged from large to small according to the row and norm thereof;
and 7: the covariance matrixes corresponding to the rows and the first 25 percent of the norm are added to obtain a new covariance matrix Rpv;
And 8: r is to bepvProjecting the image into an eigenvector matrix U to obtain a projection P ═ UHRpv(ii) a Wherein the projection P ═ UHRpvIs P in the ith column (i ═ 1, 2.., M)i=[p1i,p2i,...,pmi,...,pMi]T;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010991187.1A CN112130112B (en) | 2020-09-20 | 2020-09-20 | Information source number estimation method based on acoustic vector array joint information processing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010991187.1A CN112130112B (en) | 2020-09-20 | 2020-09-20 | Information source number estimation method based on acoustic vector array joint information processing |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112130112A true CN112130112A (en) | 2020-12-25 |
CN112130112B CN112130112B (en) | 2022-09-27 |
Family
ID=73842978
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010991187.1A Active CN112130112B (en) | 2020-09-20 | 2020-09-20 | Information source number estimation method based on acoustic vector array joint information processing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112130112B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113325365A (en) * | 2021-05-18 | 2021-08-31 | 哈尔滨工程大学 | Quaternion-based coherent signal two-dimensional DOA estimation method |
CN115469265A (en) * | 2022-09-02 | 2022-12-13 | 哈尔滨工程大学 | Acoustic vector array joint processing azimuth estimation method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103942449A (en) * | 2014-05-05 | 2014-07-23 | 北京理工大学 | Feature interference cancellation beam forming method based on estimation of number of information sources |
WO2016089300A1 (en) * | 2014-12-02 | 2016-06-09 | Thales Solutions Asia Pte Ltd. | Methods and systems for spectral analysis of sonar data |
CN105954709A (en) * | 2016-05-06 | 2016-09-21 | 哈尔滨工程大学 | Acoustic vector circular array source number detection method based on characteristic value multiple threshold correction |
CN107167776A (en) * | 2017-07-02 | 2017-09-15 | 中国航空工业集团公司雷华电子技术研究所 | The adaptive beam-forming algorithm compensated based on subspace |
CN107861114A (en) * | 2017-11-07 | 2018-03-30 | 哈尔滨工程大学 | A kind of noise suppressing method based on the reversion of underwater sound array spatial domain |
-
2020
- 2020-09-20 CN CN202010991187.1A patent/CN112130112B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103942449A (en) * | 2014-05-05 | 2014-07-23 | 北京理工大学 | Feature interference cancellation beam forming method based on estimation of number of information sources |
WO2016089300A1 (en) * | 2014-12-02 | 2016-06-09 | Thales Solutions Asia Pte Ltd. | Methods and systems for spectral analysis of sonar data |
CN105954709A (en) * | 2016-05-06 | 2016-09-21 | 哈尔滨工程大学 | Acoustic vector circular array source number detection method based on characteristic value multiple threshold correction |
CN107167776A (en) * | 2017-07-02 | 2017-09-15 | 中国航空工业集团公司雷华电子技术研究所 | The adaptive beam-forming algorithm compensated based on subspace |
CN107861114A (en) * | 2017-11-07 | 2018-03-30 | 哈尔滨工程大学 | A kind of noise suppressing method based on the reversion of underwater sound array spatial domain |
Non-Patent Citations (3)
Title |
---|
LEI HUANG: "Robust estimation of the number of sources using an MMSE-based MDL method", 《2009 IET INTERNATIONAL RADAR CONFERENCE》 * |
司伟建: "基于特征子空间投影的信源数估计方法", 《系统工程与电子技术》 * |
王琨: "估计白噪声中正弦信号个数的比率准则方法", 《西北工业大学学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113325365A (en) * | 2021-05-18 | 2021-08-31 | 哈尔滨工程大学 | Quaternion-based coherent signal two-dimensional DOA estimation method |
CN115469265A (en) * | 2022-09-02 | 2022-12-13 | 哈尔滨工程大学 | Acoustic vector array joint processing azimuth estimation method |
Also Published As
Publication number | Publication date |
---|---|
CN112130112B (en) | 2022-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sheng et al. | Maximum likelihood multiple-source localization using acoustic energy measurements with wireless sensor networks | |
CN112130112B (en) | Information source number estimation method based on acoustic vector array joint information processing | |
CN109375154B (en) | Coherent signal parameter estimation method based on uniform circular array in impact noise environment | |
CN111537982B (en) | Distortion drag array line spectrum feature enhancement method and system | |
CN113866718B (en) | Matching field passive positioning method based on mutual mass array | |
CN103323810B (en) | L-array azimuthal angle and pitch angle paired signal processing method | |
CN112731280B (en) | ESPRIT-DOA estimation method in inter-mass array mixed noise environment | |
CN113759303A (en) | Non-grid DOA (angle of arrival) estimation method based on particle swarm optimization | |
CN108088547A (en) | A kind of weak target passive detection method based on small-bore two-dimensional vector hydrophone battle array | |
CN112098938B (en) | Six-element cone vector array-based underwater acoustic target dimension reduction matching sound field positioning method | |
He et al. | Particle-velocity-field difference smoothing for coherent source localization in spatially nonuniform noise | |
CN111381212A (en) | Virtual ultrashort baseline positioning method based on subarray division | |
CN109541572B (en) | Subspace orientation estimation method based on linear environment noise model | |
CN113341371B (en) | DOA estimation method based on L array and two-dimensional ESPRIT algorithm | |
CN112415469B (en) | Rapid interference direction finding method for two-dimensional digital array radar | |
CN110632579B (en) | Iterative beam forming method using subarray beam domain characteristics | |
CN112800596A (en) | Robust dynamic direction finding method based on nested array under strong impact noise | |
US20030043695A1 (en) | System and method for stochastic characterization of sparse, four-dimensional, underwater-sound signals | |
CN112363106B (en) | Signal subspace direction of arrival detection method and system based on quantum particle swarm | |
CN112630724B (en) | High-resolution target azimuth estimation method suitable for UUV platform | |
CN113514808B (en) | Intelligent identification method for judging target number of small unmanned aerial vehicle | |
CN116879834A (en) | Extended covariance matrix two-dimensional DOA estimation method based on double parallel line arrays | |
CN114142900B (en) | Multi-channel self-adaptive beam forming method based on LCMV (liquid Crystal display television) criterion | |
CN112114287B (en) | Outlier real-time eliminating method for azimuth observation data | |
CN115825854B (en) | Underwater target azimuth estimation method, medium and system based on deep learning |
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 |