CN106768266A - A kind of measuring method of the factory building reverberation time for rejecting background noise energy - Google Patents
A kind of measuring method of the factory building reverberation time for rejecting background noise energy Download PDFInfo
- Publication number
- CN106768266A CN106768266A CN201710131962.4A CN201710131962A CN106768266A CN 106768266 A CN106768266 A CN 106768266A CN 201710131962 A CN201710131962 A CN 201710131962A CN 106768266 A CN106768266 A CN 106768266A
- Authority
- CN
- China
- Prior art keywords
- signal
- reverberation
- factory building
- reverberation time
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H7/00—Measuring reverberation time ; room acoustic measurements
Abstract
The invention provides a kind of measuring method of the factory building reverberation time for rejecting background noise energy, comprise the following steps：Microphone is arranged in factory building, sending impulse response by sound source encourages, and reverberant sound signals are acquired using microphone；Selected octave center frequency；To gather signal corresponding to frequency range be filtered, obtain the sound pressure signal under this centre frequency；Choosing one section of sufficiently long later stage plateau ambient noise signal carries out energy integral, obtains unit interval noise energy average：Long enough reverberation decay period stage signal is chosen, the degradation period acoustic energy value after rejecting background noise energy is drawn：Linear fit is carried out to sound pressure level decay curve, slope after being fitted calculates the reverberation time by slope；The present invention can eliminate the influence that ambient interferences noise was measured the reverberation time, be applied to the factory building that there is ambient interferences noise, measure its reverberation time.
Description
Technical field
The present invention relates to acoustic measurement field, more particularly to a kind of survey of the factory building reverberation time for rejecting background noise energy
Amount method.
Background technology
Reverberation effect is produced through the multiple continuous reflection of surroundings wall within the enclosed space by sound.Description reverberation effect is strong
A weak important parameter is the reverberation time.Reverberation time is not only the physics of sound attenuating speed degree in description closing space
Amount, is also an important objective evaluation parameter in acoustics design.Common reverberation time measuring method has steady sound field
The process of chopping, impulse response integration method etc..
Some factory's later stages need to measure the reverberation time in its factory building when reconstructing, but some auxiliary equipments cannot be closed, and be caused
In the presence of stronger ambient interferences noise and cannot be eliminated during the measurement reverberation time.
The content of the invention
For Shortcomings in the prior art, the invention provides a kind of factory building reverberation time for rejecting background noise energy
Measuring method, the influence that ambient interferences noise was measured the reverberation time can be eliminated, can be applied to there are ambient interferences and make an uproar
The factory building of sound, measures its reverberation time.
The present invention is to realize above-mentioned technical purpose by following technological means.
A kind of measuring method of the factory building reverberation time for rejecting background noise energy, comprises the following steps：
S01：Microphone is arranged in factory building, sending impulse response by sound source encourages, using microphone to reverberation sound
Signal is acquired, and sample frequency is fc, the signal importing computer that will be gathered；
S02：According to the selected octave center frequency f of GBT 50076-2013m=125 × 2mHz, wherein 0≤m≤5, first
Choose as m=0, fm=f0=125Hz；
S03：To gather signal with fmCorresponding frequency range is filtered, and obtains the sound pressure signal p under this centre frequencym
(τ)；
S04：From pmOne section of sufficiently long later stage plateau ambient noise signal is chosen in (τ) carries out energy integral, obtains
To unit interval noise energy average
In formula：
t3- choose ambient noise signal initial time；
t4- choose the ambient noise signal termination time；
S05：In pmLong enough reverberation decay period stage signal is chosen in (τ), being used for pulse reverse integration method is carried out
Integration, subtracts the background noise energy value corresponding to integration time, draws the degradation period acoustic energy value after rejecting background noise energy
Em(t)：
In formula：
t1- voice signal decay period initial value is the time corresponding to sound pressure level peak；
t2Certain point corresponding time in-plateau ambient noise；
To EmT () does logarithm operation and obtains sound pressure level decay curve：Lpm(t)=10lg [Em(t)]；
S06：To sound pressure level decay curve LpmT () carries out linear fit, slope k after being fittedm, calculated by slope
Reverberation time Tm(unit：s)：
S07：Work as m<When 5, repeat step two to seven draws the reverberation time T corresponding to each centre frequencym；Work as m=5
When, i.e. fm=f5=4000Hz, draws the reverberation time corresponding to each centre frequency.
Further, sound source described in S01 is starting pistol or firecracker.
Beneficial effect of the present invention：
1. the measuring method of the factory building reverberation time for rejecting background noise energy of the present invention, can eliminate background and do
Disturb the influence that noise was measured the reverberation time.
2. the measuring method of the factory building reverberation time for rejecting background noise energy of the present invention, is applied to the presence of background
The factory building of interference noise, measures its reverberation time.
3. it is of the present invention reject background noise energy the factory building reverberation time measuring method, by calculating background
Noise energy average, carries out the calculating of reverberation time again after rejecting background noise energy in reverberation signal energy.
4. the measuring method of the factory building reverberation time for rejecting background noise energy of the present invention, does not have to ambient noise
It is required that, do not limited by ambient noise power, it is a kind of more accurate reverberation time measuring method.
Brief description of the drawings
Fig. 1 is example centre frequency f of the present invention0Sound pressure signal time domain after frequency range filtering corresponding to=125Hz
Figure.
Fig. 2 is example centre frequency f of the present invention5Sound pressure signal time domain after frequency range filtering corresponding to=4000Hz
Figure.
Fig. 3 is the flow chart of the measuring method of the factory building reverberation time of rejecting background noise energy of the invention.
Specific embodiment
Below in conjunction with the accompanying drawings and specific embodiment the present invention is further illustrated, but protection scope of the present invention is simultaneously
Not limited to this.
According to Fig. 3, a kind of measuring method of the factory building reverberation time of rejecting background noise energy comprises the following steps：
S01：Microphone is arranged in factory building, sending impulse response by sound source encourages, using microphone to reverberation sound
Signal is acquired, and sample frequency is fc, the signal importing computer that will be gathered；
S02：According to the selected octave center frequency f of GBT 50076-2013m=125 × 2mHz, wherein 0≤m≤5, first
Choose as m=0, fm=f0=125Hz；
S03：To gather signal with fmCorresponding frequency range is filtered, and obtains the sound pressure signal p under this centre frequencym
(τ)；
S04：From pmOne section of sufficiently long later stage plateau ambient noise signal is chosen in (τ) carries out energy integral, obtains
To unit interval noise energy average
In formula：
t3- choose ambient noise signal initial time；
t4- choose the ambient noise signal termination time；
S05：In pmLong enough reverberation decay period stage signal is chosen in (τ), being used for pulse reverse integration method is carried out
Integration, subtracts the background noise energy value corresponding to integration time, draws the degradation period acoustic energy value after rejecting background noise energy
Em(t)：
In formula：
t1- voice signal decay period initial value is the time corresponding to sound pressure level peak；
t2Corresponding time (the t of certain point in-plateau ambient noise2Generally relative to t1Choose)；
To EmT () does logarithm operation and obtains sound pressure level decay curve：Lpm(t)=10lg [Em(t)]；
S06：To sound pressure level decay curve LpmT () carries out linear fit, slope k after being fittedm, calculated by slope
Reverberation time Tm(unit：s)：
S07：Work as m<When 5, repeat step two to seven draws the reverberation time T corresponding to each centre frequencym；Work as m=5
When, i.e. fm=f5=4000Hz, draws the reverberation time corresponding to each centre frequency.
Specific embodiment：
S01：Three groups of microphones are arranged in factory building, the source of students of pulse excitation is produced by starting pistol or firecracker, used
Microphone is acquired to reverberant sound signals, and sample frequency is fc=51200, sampling number N=2048000, by what is gathered
Signal imports computer and is ready for analysis；
S02：According to the selected octave center frequency f of GBT 50076-2013m=125 × 2mHz, wherein 0≤m≤5, first
Choose as m=0, fm=f0=125Hz；
S03：To first group of signal of microphone pick with f0Corresponding frequency range (90Hz~180Hz) is filtered, and obtains
Sound pressure signal p under to this centre frequency0(τ), as shown in Figure 1；
S04：From p0One section of sufficiently long later stage plateau ambient noise signal is chosen in (τ) carries out energy integral, obtains
To unit interval noise energy average
In formula：t3The time point corresponding to N=800000；t4The time point corresponding to N=900000.
S05：In p0Long enough reverberation decay period stage signal is chosen in (τ), being used for pulse reverse integration method is carried out
Integration, subtracts the background noise energy value corresponding to integration time, draws the degradation period acoustic energy value after rejecting background noise energy
E0(t)：
In formula：
T1-voice signal decay period initial value is the time corresponding to sound pressure level peak；
Corresponding time (the t of certain point in t2-plateau ambient noise2Generally relative to t1After increase the regular hour
Constant, is designated as：t2=t1+ C, wherein C are constant)
It is simple and direct to calculate for convenience, 200000 points are taken herein as object is analyzed, i.e.,；C=199999/fc
To E0T () does logarithm operation and obtains sound pressure level decay curve：Lp0(t)=10lg [E0(t)]；
S06：To sound pressure level decay curve Lp0T () carries out linear fit, slope k after being fitted0=-4.461, by oblique
Rate calculates reverberation time T0(unit：s)：
S07：Work as m<When 5, m=m+1, repeat step S02-S06 is made to draw the reverberation time corresponding to each centre frequency
Tm；Calculate f1=250Hz, f2=500Hz, f3=1000Hz, f4=2000Hz, f5=4000Hz, corresponding T1,T2,
T3,T4,T5；Wherein Fig. 2 is f5Sound pressure signal time-domain diagram after frequency range filtering corresponding to=4000Hz.
As m=5, circulation terminates, that is, when drawing the reverberation corresponding to each centre frequency that first group of microphone is measured
Between, by data record.
S03-S07 is repeated, the data to second group of microphone and the 3rd group of microphone measurement are analyzed calculating, as a result see
Following table：
Preferred embodiment but the present invention is not limited to above-mentioned implementation method to the embodiment for of the invention, not
In the case of substance of the invention, any conspicuously improved, replacement that those skilled in the art can make
Or modification belongs to protection scope of the present invention.
Claims (3)
1. it is a kind of reject background noise energy the factory building reverberation time measuring method, it is characterised in that comprise the following steps：
S01：Microphone is arranged in factory building, sending impulse response by sound source encourages, using microphone to reverberant sound signals
It is acquired, sample frequency is fc, the signal importing computer that will be gathered；
S02：According to the selected octave center frequency f of GBT 50076-2013m=125 × 2mHz, wherein 0≤m≤5, choose first
As m=0, fm=f0=125Hz；
S03：To gather signal with fmCorresponding frequency range is filtered, and obtains the sound pressure signal p under this centre frequencym(τ)；
S04：From pmOne section of sufficiently long later stage plateau ambient noise signal is chosen in (τ) carries out energy integral, obtains unit
Temporal noise average energy value
In formula：
t3- choose ambient noise signal initial time；
t4- choose the ambient noise signal termination time；
S05：In pmLong enough reverberation decay period stage signal is chosen in (τ), pulse reverse integration method is used for and is integrated,
The background noise energy value corresponding to integration time is subtracted, the degradation period acoustic energy value E after rejecting background noise energy is drawnm(t)：
In formula：
t1- voice signal decay period initial value is the time corresponding to sound pressure level peak；
t2Certain point corresponding time in-plateau ambient noise；
To EmT () does logarithm operation and obtains sound pressure level decay curve：Lpm(t)=10lg [Em(t)]；
S06：To sound pressure level decay curve LpmT () carries out linear fit, slope k after being fittedm, reverberation is calculated by slope
Time Tm(unit：s)：
S07：Work as m<When 5, m=m+1 is made, repetition S02-S06 steps draw the reverberation time T corresponding to each centre frequencym；Work as m
When=5, i.e. fm=f5=4000Hz, draws the reverberation time corresponding to each centre frequency.
2. it is according to claim 1 reject background noise energy the factory building reverberation time measuring method, it is characterised in that
Sound source described in S01 is starting pistol or firecracker.
3. it is according to claim 1 reject background noise energy the factory building reverberation time measuring method, it is characterised in that
T described in S052T can be expressed as2=t1+ C, wherein：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710131962.4A CN106768266B (en) | 2017-03-07 | 2017-03-07 | A kind of measurement method of workshop reverberation time that rejecting background noise energy |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710131962.4A CN106768266B (en) | 2017-03-07 | 2017-03-07 | A kind of measurement method of workshop reverberation time that rejecting background noise energy |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106768266A true CN106768266A (en) | 2017-05-31 |
| CN106768266B CN106768266B (en) | 2019-04-30 |
Family
ID=58962848
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710131962.4A Active CN106768266B (en) | 2017-03-07 | 2017-03-07 | A kind of measurement method of workshop reverberation time that rejecting background noise energy |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106768266B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107621653A (en) * | 2017-07-17 | 2018-01-23 | 上海多克实业有限公司 | Whether a kind of quick detection container is empty instrument and method |
| CN110160640A (en) * | 2018-05-18 | 2019-08-23 | 李文杰 | Acoustic power level hemisphere measuring surface method detects the method that console mode toilet seat rinses noise |
| CN110160635A (en) * | 2018-05-18 | 2019-08-23 | 李文杰 | A weighted acoustic energy grade mensuration detects sanitary ware to the method for/draining noise |
| CN110160632A (en) * | 2018-05-18 | 2019-08-23 | 李文杰 | Acoustic energy grade cylindrical measurement face method detects the method that ceramic pedestal pan rinses noise |
| CN110432885A (en) * | 2019-09-11 | 2019-11-12 | 东北大学 | A kind of photoplethysmographic noise remove method |
| CN111711914A (en) * | 2020-06-15 | 2020-09-25 | 杭州艾力特数字科技有限公司 | Sound amplification system with reverberation time measuring function |
| CN113916357A (en) * | 2021-09-03 | 2022-01-11 | 中国电力科学研究院有限公司 | Method and system for denoising station boundary noise measured value of flow conversion station |
| CN115424636A (en) * | 2022-11-04 | 2022-12-02 | 广州声博士声学技术有限公司 | Voice quality measuring method and device based on voice transmission index |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010091897A (en) * | 2008-10-10 | 2010-04-22 | Kyushu Univ | Voice signal emphasis device |
| US20110188350A1 (en) * | 2010-02-02 | 2011-08-04 | Russo Donato M | System and method for depth determination of an impulse acoustic source |
| CN102928067A (en) * | 2012-10-16 | 2013-02-13 | 华南理工大学 | System and method for measuring room acoustic parameters |
| JP2015137931A (en) * | 2014-01-22 | 2015-07-30 | ヤマハ株式会社 | reverberation time analyzer |
| CN105115585A (en) * | 2015-07-23 | 2015-12-02 | 南京大学 | Measuring method for reverberation time of reverberation room |
-
2017
- 2017-03-07 CN CN201710131962.4A patent/CN106768266B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010091897A (en) * | 2008-10-10 | 2010-04-22 | Kyushu Univ | Voice signal emphasis device |
| US20110188350A1 (en) * | 2010-02-02 | 2011-08-04 | Russo Donato M | System and method for depth determination of an impulse acoustic source |
| CN102928067A (en) * | 2012-10-16 | 2013-02-13 | 华南理工大学 | System and method for measuring room acoustic parameters |
| JP2015137931A (en) * | 2014-01-22 | 2015-07-30 | ヤマハ株式会社 | reverberation time analyzer |
| CN105115585A (en) * | 2015-07-23 | 2015-12-02 | 南京大学 | Measuring method for reverberation time of reverberation room |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107621653A (en) * | 2017-07-17 | 2018-01-23 | 上海多克实业有限公司 | Whether a kind of quick detection container is empty instrument and method |
| CN110160640A (en) * | 2018-05-18 | 2019-08-23 | 李文杰 | Acoustic power level hemisphere measuring surface method detects the method that console mode toilet seat rinses noise |
| CN110160635A (en) * | 2018-05-18 | 2019-08-23 | 李文杰 | A weighted acoustic energy grade mensuration detects sanitary ware to the method for/draining noise |
| CN110160632A (en) * | 2018-05-18 | 2019-08-23 | 李文杰 | Acoustic energy grade cylindrical measurement face method detects the method that ceramic pedestal pan rinses noise |
| CN110160635B (en) * | 2018-05-18 | 2021-02-12 | 李文杰 | Method for detecting sanitary ware water supply/drainage noise by A weighting acoustic energy level measurement method |
| CN110160640B (en) * | 2018-05-18 | 2021-02-26 | 李文杰 | Method for detecting flushing noise of floor type toilet by sound power level hemisphere measuring surface method |
| CN110432885A (en) * | 2019-09-11 | 2019-11-12 | 东北大学 | A kind of photoplethysmographic noise remove method |
| CN111711914A (en) * | 2020-06-15 | 2020-09-25 | 杭州艾力特数字科技有限公司 | Sound amplification system with reverberation time measuring function |
| CN113916357A (en) * | 2021-09-03 | 2022-01-11 | 中国电力科学研究院有限公司 | Method and system for denoising station boundary noise measured value of flow conversion station |
| CN115424636A (en) * | 2022-11-04 | 2022-12-02 | 广州声博士声学技术有限公司 | Voice quality measuring method and device based on voice transmission index |
| CN115424636B (en) * | 2022-11-04 | 2023-02-03 | 广州声博士声学技术有限公司 | Voice quality measurement method and device based on voice transmission index |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106768266B (en) | 2019-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106768266A (en) | A kind of measuring method of the factory building reverberation time for rejecting background noise energy | |
| CN105371945B (en) | Remove the method and device of noise in hall impulse response measurement acoustical signal | |
| CN103675104B (en) | The measuring method of material random incidence acoustical absorption coefficient or absorption and measurement mechanism thereof | |
| RU2013123696A (en) | ECHO SUPPRESSION CONTAINING MODELING OF LATE REVERB COMPONENTS | |
| CN104215694A (en) | Acoustic insulation testing device for fabric | |
| CN104535647A (en) | Prediction apparatus for sound absorption and insulation performance of multilayer material and method | |
| CN107785025B (en) | Noise removal method and device based on repeated measurement of room impulse response | |
| CN101403666A (en) | Dynamic elastic modulus E and damping ratio zeta measurement method for timber wood and wooden composite material | |
| Katalin | Studying noise measurement and analysis | |
| CN204228665U (en) | Fabric sound-adsorbing sound-insulating hung proving installation | |
| CN104869519A (en) | Method and system for testing background noise of microphone | |
| CN108718438B (en) | Method for judging time starting point of response signal of electroacoustic product detection system | |
| US9088857B2 (en) | Audio apparatus, control method for the audio apparatus, and storage medium for determining sudden noise | |
| CN106710602B (en) | Acoustic reverberation time estimation method and device | |
| CN104034970B (en) | A kind of sound level measurement method based on computer sound card | |
| JP2009211021A (en) | Reverberation time estimating device and reverberation time estimating method | |
| Baggenstoss et al. | Comparing shift-autocorrelation with cepstrum for detection of burst pulses in impulsive noise | |
| CN103050127A (en) | Hall objective audio quality parameter measuring system and method | |
| CN103558410B (en) | Based on the jamproof device and method of ultrasonic wind velocity indicator that extrinsic frequency excites | |
| Trémolières et al. | A multi‐parametric study of impact noise‐induced TTS | |
| Qin et al. | Measurement of field complex noise using a novel acoustic detection system | |
| CN108872402A (en) | Ultrasonic wave Butterworth, Hanning window combination with hinder filtering method | |
| JP2019032242A (en) | Attenuation time analytic method, device, and program | |
| CN107582068A (en) | A kind of distortion product otoacoustic emissions detecting strategy based on the cycle from relevant principle | |
| Yu et al. | Measurement of acoustic attenuation coefficient of stored grain |
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 |