CN103389495A - Correction method for improving resolution ratio of low frequency sound source - Google Patents
Correction method for improving resolution ratio of low frequency sound source Download PDFInfo
- Publication number
- CN103389495A CN103389495A CN2013103033515A CN201310303351A CN103389495A CN 103389495 A CN103389495 A CN 103389495A CN 2013103033515 A CN2013103033515 A CN 2013103033515A CN 201310303351 A CN201310303351 A CN 201310303351A CN 103389495 A CN103389495 A CN 103389495A
- Authority
- CN
- China
- Prior art keywords
- sound source
- low frequency
- resolution ratio
- matrix
- frequency sound
- 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
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The principle of a correction method for improving the resolution ratio of a low frequency sound source is that wave beam focusing of noise signals is accelerated through phase information of the noise signals which are connected through a control array element to improve the resolution ratio of low frequency noise. The calculation process of the correction method comprises calculating a self-spectrum and cross spectrum matrix through the noise signals which are corresponding to the analysis frequency and are connected through every measuring point; jointly determining a correction matrix according to measuring point parameters such as signal phases, target wave beam widths and current wave beam widths in the analysis frequency; calculating the sound source intensity of every grid point through the self-spectrum and cross spectrum matrix and the correction matrix and generating a wave beam cloud picture according to the calculated sound source intensity of the grid points to achieve correction of the resolution ratio of the low frequency sound source. The correction method for improving the resolution ratio of the low frequency sound source has the advantages of being capable of improving the resolution ratio of low frequency noise signals rapidly and efficiently and capable of significantly improving the resolution ratio of the low frequency sound source, greatly reducing the post-processing time of a phase array test and meeting the requirements of the phase array noise test.
Description
Technical field
The present invention relates to a kind of high-level efficiency modification method for improving phase array test low-frequency sound source resolution.
Background technology
For the test of phase array, classical beamforming algorithm is the most basic aerodynamic noise acoustic source identification research method, is also the comparison effective method.But the analysis frequency of classical beamforming algorithm is very limited, and especially the sound source resolution for lower frequency is very low.This is mainly to be subjected to the restriction of phase array array element characteristic distributions due to the resolution of classic algorithm, therefore need sometimes to take certain methods to improve the resolution of low-frequency noise signal, for example in noise fest, can be according to the array of the different designs different pore size size of analysis of experiments frequency.Sometimes, in order to guarantee the effective rate of utilization of sensor, also need the repeated disassembled and assembled sensor, the time that this has wasted test greatly, increased the cost of testing.
Along with the development of phase array experimental technique, some new phase array post-processing algorithm grow up gradually as CLEAN-PSF, CLEAN-SC, DAMAS, NNLS etc.The development and the progress that have gone out to represent the research of phase array technique of these algorithms.But because these algorithms need to iterate, counting yield is lower, and computing time is long, can not meet the needs that acoustic source identification shows in real time.
Summary of the invention
Technology of the present invention is dealt with problems: overcome the deficiencies in the prior art, a kind of modification method for improving phase array test low-frequency sound source resolution is provided, the method is by the phase information of the noise signal controlling array element and collect, the acceleration noise signal beam focuses on, thereby improves the resolution of Low Frequency Noise Generator.
Technical solution of the present invention:
A kind of modification method for improving low-frequency sound source resolution, step is as follows:
(1) according to formula
Calculate from the element G that composes on the cross-spectrum diagonal of a matrix
ii(f), according to formula
Calculate from the element G that composes on cross-spectrum matrix off-diagonal
ij(f); Wherein, p
i(f) represent that the frequency that i measuring point records is the sound-source signal of f,
For p
i(f) contrary, i, j represent to survey period, and i, j are all integer, and f represents analysis frequency.
(2) pass through formula
Determine correction matrix E, wherein,
(1) according to formula
Calculate the strength of sound source b of m analyzing spot
m, wherein, G is the cross-spectrum of the spectrum certainly matrix that obtains in step (1),
Be the control vector of m analyzing spot,
Expression
Contrary, m is positive integer;
(4) according to the strength of sound source b of all grid nodes of calculating
m, generate the wave beam cloud atlas, namely the resolution of low-frequency sound source has obtained correction.
Described measurement point refers to measure the sensing station of described low-frequency sound source signal, and the sum of sensor equals J.Described analyzing spot refers to the grid node after grid division on the sound source plane.
The invention has the advantages that:
(1) the present invention reaches by the phase place of the signal of control sound source correspondence analysis frequency the purpose that the acceleration noise signal beam focuses on, and method is simply effective, easily operation;
(2) modification method calculates and does not relate to the numerical value iteration, and computing time is short, and counting yield is higher.
Description of drawings
Fig. 1 (a) is the classical beamforming algorithm result of single sound source;
Fig. 1 (b) is the result of calculation of single sound source modification method;
Fig. 2 (a) is the classical beamforming algorithm results of two sound sources;
Fig. 2 (b) is the result of calculation of two sound source modification methods;
Fig. 3 is the inventive method process flow diagram.
Embodiment
As shown in Figure 3, the invention provides a kind of modification method for improving low-frequency sound source resolution, the method is controlled the noise signal phase information, and the acceleration noise signal beam focuses on and mainly by correction matrix, realizes., because the width size that wave beam focuses on is determined by the size of analysis frequency, the most directly be limited by the impact of wave length of sound corresponding to analysis frequency.The low frequency correction algorithm, by the difference of more current beam angle and object beam width, is determined modifying factor, and then the beam angle of control output is the object beam width.
Step is as follows:
(1) according to formula
Calculate from the element G that composes on the cross-spectrum diagonal of a matrix
ii(f), according to formula
Calculate from the element G that composes on cross-spectrum matrix off-diagonal
ij(f); Wherein, p
i(f) represent that the frequency that i measuring point records is the sound-source signal of f,
For p
i(f) contrary, i, j represent to survey period, and i, j are all integer, and f represents analysis frequency, and described measurement point refers to measure the sensing station of described low-frequency sound source signal, and the sum of sensor equals J;
(2) pass through formula
Determine correction matrix E, wherein,
(3) according to formula
Calculate the strength of sound source b of m analyzing spot
m, wherein, G is the cross-spectrum of the spectrum certainly matrix that obtains in step (1),
Be the control vector of m analyzing spot,
Expression
Contrary, m is positive integer;
(4) according to the strength of sound source b of all grid nodes of calculating
m, generate the wave beam cloud atlas, namely the resolution of low-frequency sound source has obtained correction.
Fig. 1, Fig. 2 are the result contrast of classic algorithm and modification method, and in figure, the concentric circles district represents the low-frequency sound source distribution range.Usually sound source resolution is all to weigh by the sound source distribution range, and the sound source distribution range is wider, and sound source resolution is lower.Therefore as can be seen from the figure, no matter for single low-frequency sound source or the situation of two low-frequency sound sources, modification method to the resolution of sound source all far away higher than the result of classic algorithm.The invention of modification method has made up the deficiency of classical beamforming algorithm, has expanded the analysis frequency scope of phase array test.Because the method does not relate to the folded band of numerical value, counting yield other optimized algorithm (CLEAN-PSF, CLEAN-SC, DAMAS, NNLS) relatively is high a lot, has effectively met the needs that acoustic source identification shows in real time.
Claims (3)
1. one kind for the modification method that improves low-frequency sound source resolution, it is characterized in that step is as follows:
(1) according to formula
Calculate from the element G that composes on the cross-spectrum diagonal of a matrix
ii(f), according to formula
Calculate from the element G that composes on cross-spectrum matrix off-diagonal
ij(f); Wherein, p
i(f) represent that the frequency that i measuring point records is the sound-source signal of f,
For p
i(f) contrary, i, j represent to survey period, and i, j are all integer, and f represents analysis frequency.
(2) pass through formula
Determine correction matrix E, wherein,
(3) according to formula
Calculate the strength of sound source b of m analyzing spot
m, wherein, G is the cross-spectrum of the spectrum certainly matrix that obtains in step (1),
Be the control vector of m analyzing spot,
Expression
Contrary, m is positive integer;
(4) according to the strength of sound source b of all grid nodes of calculating
m, generate the wave beam cloud atlas, namely the resolution of low-frequency sound source has obtained correction.
2. a kind of modification method for improving low-frequency sound source resolution according to claim 1, it is characterized in that: described measurement point refers to measure the sensing station of described low-frequency sound source signal, and the sum of sensor equals J.
3. a kind of modification method for improving low-frequency sound source resolution according to claim 1, it is characterized in that: described analyzing spot refers to the grid node after grid division on the sound source plane.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310303351.5A CN103389495B (en) | 2013-07-18 | 2013-07-18 | Correction method for improving resolution ratio of low frequency sound source |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310303351.5A CN103389495B (en) | 2013-07-18 | 2013-07-18 | Correction method for improving resolution ratio of low frequency sound source |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103389495A true CN103389495A (en) | 2013-11-13 |
| CN103389495B CN103389495B (en) | 2015-04-22 |
Family
ID=49533819
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310303351.5A Active CN103389495B (en) | 2013-07-18 | 2013-07-18 | Correction method for improving resolution ratio of low frequency sound source |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103389495B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113362808A (en) * | 2021-06-02 | 2021-09-07 | 云知声智能科技股份有限公司 | Target direction voice extraction method and device, electronic equipment and storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0394743A (en) * | 1989-02-08 | 1991-04-19 | Centre Natl De La Rech Scient | Method and apparatus for localizing reflective target acoustically |
| GB2386947A (en) * | 2002-03-27 | 2003-10-01 | Qinetiq Ltd | Calibration of a multichannel receiver |
| CN101149435A (en) * | 2007-10-23 | 2008-03-26 | 中国船舶重工集团公司第七一五研究所 | U-shaped array beam forming weighting method |
| EP2017388A2 (en) * | 2001-08-14 | 2009-01-21 | Maba Fertigteilindustrie GmbH | Concrete wall element |
-
2013
- 2013-07-18 CN CN201310303351.5A patent/CN103389495B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0394743A (en) * | 1989-02-08 | 1991-04-19 | Centre Natl De La Rech Scient | Method and apparatus for localizing reflective target acoustically |
| EP2017388A2 (en) * | 2001-08-14 | 2009-01-21 | Maba Fertigteilindustrie GmbH | Concrete wall element |
| GB2386947A (en) * | 2002-03-27 | 2003-10-01 | Qinetiq Ltd | Calibration of a multichannel receiver |
| CN101149435A (en) * | 2007-10-23 | 2008-03-26 | 中国船舶重工集团公司第七一五研究所 | U-shaped array beam forming weighting method |
Non-Patent Citations (1)
| Title |
|---|
| 赵小见等: "基于相干声源ＣＬＥＡＮ算法的常规风洞声源辨识研究", 《空气动力学学报》, vol. 31, no. 2, 30 April 2013 (2013-04-30), pages 239 - 243 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113362808A (en) * | 2021-06-02 | 2021-09-07 | 云知声智能科技股份有限公司 | Target direction voice extraction method and device, electronic equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103389495B (en) | 2015-04-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108802674B (en) | Joint search method and device for direct positioning | |
| JP6508336B2 (en) | Fourier transform mass spectrometry | |
| CN105403860B (en) | A kind of how sparse sound localization method related based on domination | |
| KR102341760B1 (en) | Systems and methods for parallelizing and pipelining a tunable blind source separation filter | |
| Herold et al. | A comparison of microphone array methods for the characterization of rotating sound sources | |
| Le Courtois et al. | Genetic optimisation of a plane array geometry for beamforming. Application to source localisation in a high speed train | |
| CN107462317B (en) | A kind of water surface target radiation noise simulation method under complex scene | |
| CN106205637B (en) | Noise detection method and device for audio signal | |
| CN109188362A (en) | A kind of microphone array auditory localization signal processing method | |
| JP6763332B2 (en) | Sound collectors, programs and methods | |
| CN109343003B (en) | Method for identifying sound source formed by fast iterative shrinking wave beams | |
| CN104811886A (en) | Phase difference measurement-based microphone array direction finding method | |
| CN106842112A (en) | Sound localization method based on parametrization Bayes's dictionary learning under strong reverberant ambiance | |
| CN103389495B (en) | Correction method for improving resolution ratio of low frequency sound source | |
| US10757519B2 (en) | Neural network-based parameter estimation of loudspeakers | |
| JP2006279740A (en) | Apparatus for measuring sound strength by sound sources and sound source location analyzing method | |
| CN104656055B (en) | A kind of single signal angle-of- arrival estimation method based on large-scale multi-antenna system | |
| JP2015079080A (en) | Sound source position estimation device, method, and program | |
| Pierzga et al. | Studies on infrasound noise generated by operation οf low-power wind turbine | |
| CN109471073B (en) | NLFM signal generation method and device based on augmented Lagrange particle swarm optimization | |
| CN106779091B (en) | A kind of periodic vibration signal localization method based on transfinite learning machine and arrival distance | |
| Sarradj | Three-dimensional acoustic source mapping | |
| CN116306359A (en) | Random wave simulation method for multidimensional space non-stationary non-uniform wind field | |
| JP5826502B2 (en) | Sound processor | |
| CN112162235A (en) | Smooth segmented stochastic resonance enhanced acoustic vector signal orientation method |
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 |