CN1202407C - Method for identifying non-stationary acoustic source characteristics by bias coherent technology - Google Patents

Method for identifying non-stationary acoustic source characteristics by bias coherent technology Download PDF

Info

Publication number
CN1202407C
CN1202407C CN 03129404 CN03129404A CN1202407C CN 1202407 C CN1202407 C CN 1202407C CN 03129404 CN03129404 CN 03129404 CN 03129404 A CN03129404 A CN 03129404A CN 1202407 C CN1202407 C CN 1202407C
Authority
CN
China
Prior art keywords
alpha
reference source
signal
source
stationary
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.)
Expired - Fee Related
Application number
CN 03129404
Other languages
Chinese (zh)
Other versions
CN1472515A (en
Inventor
蒋伟康
万泉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN 03129404 priority Critical patent/CN1202407C/en
Publication of CN1472515A publication Critical patent/CN1472515A/en
Application granted granted Critical
Publication of CN1202407C publication Critical patent/CN1202407C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Holo Graphy (AREA)

Abstract

The present invention relates to a method for identifying non-stationary acoustic source characteristics by bias coherent technology, which belongs to the field of the noise class in the physics. The present invention is used for the occasion of the complex circulation stationary acoustic field of a plurality of acoustic sources, a smooth circulation theory replaces Fourier transformation, bias coherent analysis technology which is suitable for the circulation stationary acoustic field is put forward and is used for the separation of the part of the acoustic field of the circulation stationary acoustic field. The present invention comprises the method as follows: firstly the number and the position of a reference source are determined, the microphone array of the reference source is arranged, the signal of the reference source is extracted, then the microphone array is designed again to scan and measure a holographic surface and collect the data of the holographic surface; then the part of the acoustic field of each acoustic source is separated by means of the bias coherent analysis technology of the stationary circulation. The present invention can analyze the complex circulation stationary acoustic field which is formed by the a plurality of circulation stationary acoustic sources, and can obtain the part of the acoustic field which is respectively formed by each acoustic source, three-dimensional acoustic field distribution which is formed by each acoustic source can be calculated by the signals of acoustic pressure which is measured from a holographic plane, and the visualization of the propagation path of each acoustic source can be realized.

Description

Adopt the method for inclined to one side coherent technique identification non-stationary sound source characteristic
Technical field
What the present invention relates to is a kind of method of identifying sound source characteristic, and particularly a kind of method of utilizing the inclined to one side phase stem portion sound field separation technique identification non-stationary sound source characteristic of cyclo-stationary sound field belongs to the noise field in the physics class.
Background technology
Understanding to the noise source characteristic is the prerequisite of control noise, therefore, in order to control noise effectively, before noise reduction measure is implemented, must at first carry out the noise source diagnosis, determines each position, overriding noise seedbed, and characteristic.Along with the development of modern signal processing technology ground, spectral analysis technology, relevant and partial coherence analysis technology, sound intensity analytical technology and sound near-field holography technology etc. have obtained developing by leaps and bounds.Find by literature search, M.A.Tomlinson is at " Applied Acoustics " (57 (1999): write articles " Partial source discrimination innear field acoustic holography " (" applied acoustics " 243-261), part source identification in the near field acoustic holography), this article has proposed the part sound source isolation technics of near field acoustic holography partial coherence analysis, can diagnose noise source effectively.But these technology can only be analyzed steady sound field, therefore, are necessary to propose new technology, are used for the analysis of non-stationary sound field.Yet, for general non-stationary sound field, become when the statistical property parameter of acoustical signal is, thereby also replace ensemble average with regard to unrenewable time average, make data acquisition very difficult, be difficult to analyze the characteristic of sound field.
The cyclo-stationary signal is the special non-stationary signal of a class because the cyclic stationary of self uniqueness, make single acquisition to record have the cycle ergodic property, this specific character provides non-stationary signal to analyze possibility.The cyclo-stationary signal has crucial realistic meaning in engineering is used, for example rotating machinery is because the physical arrangement of symmetry or near symmetrical and periodic working motion pattern, its sound field has the obvious periodic time varying characteristic, and acoustical signal has cyclostationarity.In further retrieving, find as yet and the identical or similar bibliographical information of theme of the present invention.
Summary of the invention
The objective of the invention is to overcome deficiency of the prior art, a kind of method of identification non-stationary sound source characteristic is provided, achieve and obtain the part sound field that each sound source forms respectively, near-field holography by the cyclo-stationary sound field is rebuild, can also obtain three-dimensional sound field by the measurement data on the plane and distribute, realize that each sound source route of transmission is visual.
The present invention is achieved by the following technical solutions, the present invention is in the occasion of the steady sound field of complex loops of a plurality of sound sources generations, adopt the theoretical alternative traditional Fourier transform of cyclo-stationary, proposition is applicable to the partial coherence analysis technology of cyclo-stationary sound field, and it part sound field that is used for the cyclo-stationary sound field separated, at first determine the quantity and the position of reference source, and layout reference source microphone array, extract reference source signal, design microphone array again holographic facet is carried out scanning survey, gather the holographic facet data, then, utilize circulation partial coherence analysis technology, separate the part sound field of each sound source, obtain the relevant Density Distribution of spectrum of part sound field.
Below the inventive method is done further to limit, method step is as follows:
1, determines the quantity and the position of reference source, and arrange and extract reference source signal by the reference source microphone array.For this reason, need to determine the number K and the position of reference source, can adopt the near field compbined test analysis (or claiming inclined to one side odd value analysis) of steady sound field to obtain.
2, arrange holographic facet measuring microphone array,, can on whole hologram plane, arrange microphone if port number is abundant; If port number is not enough, microphone can be arranged to linear array, measure in the enterprising line scanning of whole holographic facet; In the holographic facet data acquisition, utilize the reference source array to gather the reference source acoustical signal.During measurement, the acoustical signal of the whole microphone passages of synchronous acquisition, recording storage are in computing machine, magnetic tape recorder or other equipment, and establishing the holographic facet port number is Q.
3, analyze the acoustical signal of being gathered, choose the frequency and the cycle frequency that can reflect sound field characteristic, can analyze the spectral density function of each reference source signal and select.
4, utilize the holographic facet time domain acoustical signal data collect and the time domain acoustical signal data of reference source, calculate the circulation of the reference source signal on selected frequency f and the cycle frequency α and compose relevant density matrix (S Rr α(f)) K * K, and the relevant density matrix (S of circulation cross-spectrum of microphone signal on reference source signal and the holographic facet Rp α(x, y, z h, f)) K * QAnd (S Pr α(x, y, z h, f)) Q * KCan be by these three the relevant density matrix of circulation spectrum, the circulation that obtains acoustical signal on the holographic facet is from the relevant density matrix (S of spectrum Pp α(x, y, z h, f)) Q * 1, computing method are suc as formula (1):
( S pp α ( f ) ) Q × 1 = diag ( ( S pr α ( f ) ) Q × K inv ( ( S rr α ( f ) ) K × K ) ( S rp α ( f ) ) K × Q ) - - - - - - ( 1 )
Wherein, the diagonal element of matrix is extracted in diag () expression, and inv () represents that this relational expression is set up to matrix inversion in whole sound field, this formula has reflected the time phase relation of reference source and holographic facet information, can be used to scanning collection to holographic facet acoustic pressure information carry out phase-locking.
5, utilize holographic facet data and the reference source data that collect, the partial coherence analysis that circulates separates the part sound field that obtains each sound source, and method is as follows:
Set mark r and represent r reference source signal, the signal of measuring point p on the hologram plane that subscript p represents to measure, then after the influence of removing r-1 reference source, r reference source in the relevant density of spectrum of the output signal of holographic facet measuring point place generation is:
S p : r · ( r - 1 ) ! α ( f ) = γ rp · ( r - 1 ) ! α 2 ( f ) S pp · ( r - 1 ) ! α ( f ) - - - - - - ( 2 )
Wherein:
γ rp · ( r - 1 ) ! α 2 ( f ) = S rp · ( r - 1 ) ! α ( f ) S pr · ( r - 1 ) ! α ( f ) S rr · ( r - 1 ) ! α ( f ) S pp · ( r - 1 ) ! α ( f ) Be the inclined to one side coherence function that circulates.
The present invention has substantive distinguishing features and marked improvement, the present invention utilizes the periodicity of the uniqueness of cyclo-stationary sound field, on the basis of partial coherence analysis technology, circulation partial coherence analysis technology has been proposed, and on the basis of the inclined to one side coherent technique harmony near-field holography technology of circulation, invented circulation partial coherence analysis part sound field separation technique, can analyze, obtained the part sound field that each sound source forms respectively the steady sound field of complex loops that a plurality of cyclo-stationary sound sources form; Near-field holography by the cyclo-stationary sound field is rebuild, and can also extrapolate the three-dimensional sound field distribution that single sound source forms the sound pressure signal that records on hologram plane, realizes that each sound source route of transmission is visual.
Embodiment
Provide following examples in conjunction with the inventive method content:
1, supposes known sound source number and position, adopt two loudspeaker sounding to form two cyclo-stationary sound sources, arrange two microphones, extract reference source signal.
The driving source of loudspeaker is:
Vsourcel=Acos(2πf 1t)*noise(t)
Vsource2=B(1+Ccos(2πf bt))*cos(2πf at)
Wherein, A=10, B=C=1, f1=300, f a=300, f b=120, noise is connected with coloured noise for band.
2,25 microphones are arranged to linear array, measure, form 25 * 25 holographic facet array in the enterprising line scanning of whole holographic facet; In the holographic facet data acquisition, utilize two reference source microphones to gather the reference source acoustical signal.During measurement, the acoustical signal of the whole microphone passages of synchronous acquisition, recording storage is in computing machine or other equipment.
3, analyze the acoustical signal of being gathered in the laboratory playback, choose the frequency f=120Hz and the cycle frequency alpha=600Hz that can reflect sound field characteristic.
4, utilize holographic facet data and the reference source data collect, the circulation of analyzing the reference source signal on selected frequency f and the cycle frequency α is from composing relevant density matrix (S Rr α(f)) K * K, the relevant density matrix (S of circulation cross-spectrum of reference source signal and holographic facet upper sensor signal Rp α(x, y, z h, f)) K * QAnd (S Pr α(x, y, z h, f)) Q * K
5, separate the part sound field that obtains each sound source.
6, utilize the near-field holography reconstruction formula of cyclo-stationary sound field to carry out sound field rebuilding, to obtain the three-dimensional information of each sound source part sound field.
Find when two reference source signal coherences that obtain are smaller in this embodiment, separating effect is very good, can analyze significantly and obtain two sound sources sound field separately, when two reference source signal coherences that obtain are very strong, separating effect is still very good, and the inherent coherence that can weaken to a great extent when adopting circulation partial coherence analysis separating part sound field between the reference source array is described.

Claims (4)

1, a kind of method of identification non-stationary sound source characteristic, it is characterized in that the occasion of the steady sound field of complex loops that produces in a plurality of sound sources is at first determined the quantity and the position of reference source, and layout reference source microphone array, extract reference source signal, design microphone array again holographic facet is carried out scanning survey, gather the holographic facet data, then, utilize cyclo-stationary partial coherence analysis technology, separate the part sound field of each sound source, obtain the relevant Density Distribution of spectrum of part sound field.
2, the method for identification non-stationary sound source characteristic according to claim 1 is characterized in that, below the present invention is made further qualification, and concrete method step is as follows:
(1) adopts the near field compbined test analysis of steady sound field or claim quantity and the position that inclined to one side odd value analysis is determined reference source, and arrange and extract reference source signal by the reference source microphone array;
(2) arrange holographic facet measuring microphone array, this array can be plane or curved array, perhaps use the linear array scanning survey, in the holographic facet data acquisition, utilize the reference source array to gather the reference source acoustical signal, during measurement, the acoustical signal of the whole microphone passages of synchronous acquisition, recording storage is in registering instrument or computing machine, and establishing the holographic facet port number is Q;
(3) analyze the acoustical signal of being gathered, choose the frequency and the cycle frequency of reflection sound field characteristic by the spectral density function of analyzing each reference source signal;
(4) utilize the holographic facet time domain acoustical signal data collect and the time domain acoustical signal data of reference source, calculate selected frequency f and compose relevant density matrix with the reference source signal circulation on the cycle frequency α ( S rr α ( f ) ) K × K , And the relevant density matrix of circulation cross-spectrum of microphone signal on reference source signal and the holographic facet ( S rp α ( x , y , z h , f ) ) K × Q And ( S pr α ( x , y , z h , f ) ) Q × K , Can be by these three the relevant density matrix of circulation spectrum, the circulation that obtains acoustical signal on the holographic facet is from the relevant density matrix of spectrum ( S pp α ( x , y , z h , f ) ) Q × 1 , Computing method are as follows:
( S pp α ( f ) ) Q × 1 = diag ( ( S pr α ( f ) ) Q × K inv ( ( S rr α ( f ) ) K × K ) ( S rp α ( f ) ) K × Q )
3, the method for identification non-stationary sound source characteristic according to claim 2 is characterized in that, in the step (4), utilizes the holographic facet data and the reference source data that collect, and the partial coherence analysis that circulates separates the part sound field that obtains each sound source, and method is as follows:
If S represents the spectral density function that circulates, subscript r represents r reference source signal, the signal of measuring point p on the hologram plane that subscript p represents to measure, subscript α represents cycle frequency, f represents the frequency of acoustical signal, the same meaning of symbol in subscript ": " and " " and the common inclined to one side coherence function, wherein ": " represents the coherent relationships between two signals, r-1 is removed in following target mid point " r (r-1) " expression from r influence, then after the influence of removing r-1 reference source, r reference source computing, the relevant density of circulation spectrum of the output signal that produces at holographic facet measuring point place is:
S p : r · ( r - 1 ) ! α ( f ) = γ rp · ( r - 1 ) ! α 2 ( f ) S pp · ( r - 1 ) ! α ( f )
Wherein:
γ rp · ( r - 1 ) ! α 2 ( f ) = S rp · ( r - 1 ) ! α ( f ) S pr · ( r - 1 ) ! α ( f ) S rr · ( r - 1 ) ! α ( f ) S pp · ( r - 1 ) ! α ( f ) Be the inclined to one side coherence function that circulates.
4, the method for identification non-stationary sound source characteristic according to claim 1 and 2, it is characterized in that, near-field holography by the cyclo-stationary sound field is rebuild, and obtains the three-dimensional sound field distribution that single sound source forms the sound pressure signal that records on hologram plane, realizes that each sound source route of transmission is visual.
CN 03129404 2003-06-19 2003-06-19 Method for identifying non-stationary acoustic source characteristics by bias coherent technology Expired - Fee Related CN1202407C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 03129404 CN1202407C (en) 2003-06-19 2003-06-19 Method for identifying non-stationary acoustic source characteristics by bias coherent technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 03129404 CN1202407C (en) 2003-06-19 2003-06-19 Method for identifying non-stationary acoustic source characteristics by bias coherent technology

Publications (2)

Publication Number Publication Date
CN1472515A CN1472515A (en) 2004-02-04
CN1202407C true CN1202407C (en) 2005-05-18

Family

ID=34153519

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 03129404 Expired - Fee Related CN1202407C (en) 2003-06-19 2003-06-19 Method for identifying non-stationary acoustic source characteristics by bias coherent technology

Country Status (1)

Country Link
CN (1) CN1202407C (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100442030C (en) * 2006-10-27 2008-12-10 合肥工业大学 A separating method for sound field

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101936818B (en) * 2010-08-27 2012-09-05 上海交通大学 Diagnostic system of non-contact type rotary mechanical failure
CN102333265B (en) * 2011-05-20 2014-02-19 南京大学 Replay method of sound fields in three-dimensional local space based on continuous sound source concept
CN106052849B (en) * 2016-05-20 2020-02-18 西南交通大学 Method for identifying non-stationary abnormal noise source in automobile
CN106872019A (en) * 2017-01-20 2017-06-20 湖北文理学院 A kind of part based on particle vibration velocity decomposition method
CN109409341A (en) * 2018-12-10 2019-03-01 中国航发四川燃气涡轮研究院 A kind of aero-engine noise source discrimination method near field
CN112729528B (en) * 2020-12-07 2022-09-23 潍柴动力股份有限公司 Noise source identification method, device and equipment

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100442030C (en) * 2006-10-27 2008-12-10 合肥工业大学 A separating method for sound field

Also Published As

Publication number Publication date
CN1472515A (en) 2004-02-04

Similar Documents

Publication Publication Date Title
CN101251412B (en) Method for rebuilding circulation calm sound source by overlapping spherical wave
CN100582690C (en) Method for rebuilding circulated steady sound source by boundary element method
CN101865789B (en) Fault detecting device of near field acoustic holography sound image mode identification and detecting method thereof
Goumas et al. Classification of washing machines vibration signals using discrete wavelet analysis for feature extraction
CN1207700C (en) Near field acoustic holographic method of disting uishing non-stationary sound source
US20060190180A1 (en) Petroleum exploration and prediction apparatus and method
CN110058305A (en) A kind of DAS seismic data noise-reduction method based on convolutional neural networks
CN1202408C (en) Process for identifying non-stable acoustical source characteristic applying principal component analyzing technique
CN1202407C (en) Method for identifying non-stationary acoustic source characteristics by bias coherent technology
US6958950B2 (en) Method of visualizing sound fields of individual sound sources using acoustic holography
US20120039150A1 (en) Unique seismic source encoding
Zheng et al. Deep neural networks design and analysis for automatic phase pickers from three-component microseismic recordings
Bouguerriou et al. Novel cyclostationarity-based blind source separation algorithm using second order statistical properties: Theory and application to the bearing defect diagnosis
CN109765303A (en) The detection method for degree of coming to nothing behind for liner structure
CN116576956A (en) Multisource vibration signal separation method based on distributed optical fiber acoustic wave sensing
Fagerström et al. One-to-many conversion for percussive samples
CN114265118A (en) Method, device and system for extracting time difference of offshore acoustic logging while drilling
CN103239235A (en) Multifunctional comprehensive otoacoustic emission detector
Wu et al. Audio-based music visualization for music structure analysis
CA2378839A1 (en) Improved method for determining subscriber loop make-up
Hou et al. Bearing fault diagnosis based on spatial features of 2.5 dimensional sound field
CN109917459A (en) Method, device and system for suppressing earthquake noise
Guilment et al. Modeling the acoustic repertoire of Cuvier's beaked whale clicks
Dell'Aversana A bridge between geophysics and digital music–applications to hydrocarbon exploration
Wang et al. Identification of the defective transmission devices using the wavelet transform

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee