CN107785025A - Noise remove method and device based on room impulse response duplicate measurements - Google Patents

Noise remove method and device based on room impulse response duplicate measurements Download PDF

Info

Publication number
CN107785025A
CN107785025A CN201610726044.1A CN201610726044A CN107785025A CN 107785025 A CN107785025 A CN 107785025A CN 201610726044 A CN201610726044 A CN 201610726044A CN 107785025 A CN107785025 A CN 107785025A
Authority
CN
China
Prior art keywords
signal
denoising
module
frequency
impulse response
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
Application number
CN201610726044.1A
Other languages
Chinese (zh)
Other versions
CN107785025B (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 Sonic Acoustic Engineering Technology Ltd By Share Ltd
Original Assignee
Shanghai Sonic Acoustic Engineering Technology Ltd By Share Ltd
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 Sonic Acoustic Engineering Technology Ltd By Share Ltd filed Critical Shanghai Sonic Acoustic Engineering Technology Ltd By Share Ltd
Priority to CN201610726044.1A priority Critical patent/CN107785025B/en
Publication of CN107785025A publication Critical patent/CN107785025A/en
Application granted granted Critical
Publication of CN107785025B publication Critical patent/CN107785025B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • G10L21/0216—Noise filtering characterised by the method used for estimating noise
    • G10L21/0224—Processing in the time domain

Abstract

The present invention relates to a kind of noise remove method and device based on room impulse response duplicate measurements, including：1) signal emission module sends swept-frequency signal；2) signal acquisition module collection room acoustic pressure；3) denoising prepulse response computing module calculates noisy impulse response；4) signal to noise ratio determination module carries out signal to noise ratio judgement；5) pendulous frequency determining module determines pendulous frequency；6) signal pre-analysis module selection denoising method；7) superposition denoising module is overlapped to multiple measurement signal；8) time-frequency coefficients reconstruct denoising module carries out time-frequency coefficients reconstruct to multiple measurement signal；9) denoising afterpulse response computing module calculates the room impulse response after denoising；10) acoustical parameters computing module calculates acoustical parameters；11) result is shown and memory module, shows test signal feature and acoustical parameters before and after denoising, and preserve result of calculation.Compared with prior art, the present invention can effectively improve the anti-noise ability and precision of live sound quality measurement.

Description

Noise remove method and device based on room impulse response duplicate measurements
Technical field
The present invention relates to the e measurement technology of field of acoustics, more particularly, to a kind of based on room impulse response duplicate measurements The method and device of noise remove.
Background technology
Impulse response is to reflect the important physical amount of auditorium, and the reverberation time can be directly calculated from impulse response (RT60), early stage decay time (EDT), the definition factor (D50), the perspectivity factor (C50), the loudness factor (G), ears are related The auditorium parameters such as coefficient (IACC).Hall impulse response acoustical signal is accurately measured to check and accept auditorium acoustical design, tonequality Development and engineering practice with environmental noise abatement technology is significant.During live impulse response measurement often by To voice, traffic and construction equipment noise jamming, particularly there are scene or the in-site measurement of construction period of spectators, noise is to surveying The interference of examination is almost inevitable.For the measurement of impulse response, to ensure measurement accuracy, typically require that actual acoustic measured Signal to noise ratio higher Cheng Zhongyou.As provided in ISO3382, the measurement for the reverberation time, 35dB feelings are reached in signal to noise ratio Measurement is just effective under condition, but in addition to laboratory and the late into the night quiet environment, in-site measurement is extremely difficult to so high letter Make an uproar ratio.Test invalidation can generally be caused in the case where signal to noise ratio is less than 25dB, therefore anti-noise technology is for impulse response scene Measurement is very necessary.
The content of the invention
It is an object of the present invention to overcome the above-mentioned drawbacks of the prior art and provide one kind accurately and fast, certainly The high noise remove method and device based on room impulse response duplicate measurements of dynamicization degree, in the range of 100~8000Hz, The signal to noise ratio measured needed for impulse response 10dB or so be can be reduced to, the precision and reliability of in-site measurement substantially increased.
The purpose of the present invention can be achieved through the following technical solutions：
A kind of method of the noise remove based on room impulse response duplicate measurements, it is characterised in that comprise the following steps：
1) signal emission module by the data signal generation unit that is sequentially connected, digital-to-analogue conversion card, power amplifier and Loudspeaker sends linear frequency sweep acoustical signal into room；
2) signal acquisition module is connect by microphone, preamplifier, analog-to-digital conversion card and the data signal being sequentially connected Unit is received, gathers the sound pressure signal at room measuring point, and sends it to denoising prepulse response computing module；
3) signal and the data signal of swept-frequency signal that collect are made computing by denoising prepulse response computing module, are obtained The impulse response in room, and it is transferred to signal to noise ratio determination module；
4) signal to noise ratio under initial condition is calculated using the response of denoising prepulse for signal to noise ratio determination module, judges that it is The no snr threshold for reaching acoustical parameters calculating, if it is satisfied, then sending impulse response to acoustical parameters computing module, such as It is unsatisfactory for, then sends signal to noise ratio to pendulous frequency determining module；
5) pendulous frequency determining module determines duplicate measurements number according to signal to noise ratio condition；
6) signal pre-analysis module carries out time-frequency characteristic analysis to the signal collected, and selects noise-reduction method accordingly：It is folded Denoising or time-frequency coefficients are added to reconstruct denoising；
7) signal repeatedly measured is overlapped by superposition denoising module, is realized noise suppressed and useful signal enhancing, is obtained Signal after to denoising, and denoising afterpulse response computing module is sent to, and perform step 9)；
8) signal repeatedly measured is carried out time-frequency coefficients reconstruct by coefficient reconstruct denoising module, realizes that noise is gone in environment Remove, obtain the signal after denoising, and be sent to denoising afterpulse response computing module, and perform step 9)；
9) sound pressure signal after denoising and the data signal of swept-frequency signal are made computing by denoising afterpulse response computing module, The room impulse response after denoising is obtained, and is sent to acoustical parameters computing module；
10) acoustical parameters computing module calculates acoustical parameters using the forward and backward impulse response of denoising, and is sent to result and shows Show and memory module；
11) result is shown and memory module shows sound pressure signal curve before and after denoising, pulse respond, under sound pressure level Curve, acoustical parameters result of calculation are dropped, and stores its data signal and result of calculation.
A kind of method of the noise remove based on room impulse response duplicate measurements, it is characterised in that launch including signal Module, signal acquisition module, denoising prepulse response computing module, signal to noise ratio determination module, pendulous frequency determining module, superposition Denoising module, time-frequency coefficients reconstruct denoising module, denoising afterpulse response computing module, acoustical parameters computing module and result show Show and memory module.
Described signal emission module includes the data signal generation being sequentially connected and weighted units, digital-to-analogue conversion card, work( Rate amplifier and loudspeaker, for completing the task to hall internal emission line frequency sweep acoustical signal.
Described data signal is generated as a data signal generation unit, described data signal generation unit generation signal Calculating process it is as follows：
Wherein, A is signal amplitude；T is signal total duration；T is time variable；f1It is off lower-frequency limit；f2It is off frequency The rate upper limit.
Selection signal transmission channel and acquisition channel；Input transmission signal duration, collection signal duration, sample frequency, letter Number amplitude, the initial frequency of swept-frequency signal and cut-off frequency.
Described signal acquisition module includes microphone, preamplifier, analog-to-digital conversion card and the numeral letter being sequentially connected Number receiving unit, to complete the task that sound pressure signal gathers in room, and will collection signal be sent to signal pre-analysis module and Denoising prepulse responds computing module.
Described signal acquisition part uses triggering collection, ensures that repeatedly measure the signal collected corresponds point by point, And the signal collected is sent to denoising prepulse response computing module.
Described denoising prepulse response computing module is to utilize the digital linear frequency sweep letter modulated in signal emission module Number e0(t) the sound pressure signal x (t) collected with measuring point is calculated, and obtains noisy acoustic pressure impulse response p of the measuring point on room (t), and it is transferred to signal to noise ratio determination module.
Described signal to noise ratio determination module calculates signal to noise ratio using impulse response, determines whether to reach according to signal to noise ratio size Snr threshold：I.e. if signal to noise ratio is more than 25dB, then acoustical parameters computing module is sent to, is such as unsatisfactory for, be then sent to survey Measure number determining module.
Described pendulous frequency determining module, duplicate measurements number is determined according to signal to noise ratio size：Signal to noise ratio be 15- 25dB duplicate measurements 5 times, the duplicate measurements that signal to noise ratio is 5-15dB 10 times.
Described signal preanalysis is to carry out VARIANT FILTER OF DISCRETE GABOR TRANSFORM to the signal collected：
Wherein, W=exp (j2 π/N)；It is the amplitude of time frequency analysis coefficient；It is the acoustic pressure time-domain signal of collection； Δ M is time-domain analysis step-length；Δ N is frequency-domain analysis step-length；It is the analytic function of Gabor transformation；Pass through signal Preanalysis is obtained gathering signal time frequency analysis, and intuitive judgment is provided to noise characteristic.
Window function type, window function width, time sampling step-length, frequency sampling step-length are selected during signal preanalysis.
Described superposition denoising module, for the signal repeatedly measured to be carried out into time domain superposition, realize that useful signal strengthens Suppress with flucturation noise, the signal after processing is then sent to denoising afterpulse response computing module.
Described time-frequency coefficients reconstruct denoising module, for the signal repeatedly measured to be carried out into time-frequency coefficients reconstruct, then Signal after processing is sent to denoising afterpulse response computing module.
Described denoising afterpulse response computing module is to utilize the digital linear frequency sweep letter modulated in signal emission module Number e0(t) make correlation computations with the sound pressure signal y ' (t) after denoising, obtain the room impulse response after denoising, and by its Send acoustical parameters computing module to；
Described acoustical parameters computing module utilizes the impulse response before and after denoising to calculate the reverberation time (RT60), decline in early days Become time (EDT), the definition factor (D50), the perspectivity factor (C50), the loudness factor (G), the Room of subjective assessment (IACC) Hall acoustical parameters.
Described result display module includes the time-domain signal display unit before the gathered denoising of display, acoustic energy decay song Impulse response display unit before line display unit, denoising, time-domain signal display unit, acoustic energy decay curve after denoising are shown Impulse response display unit after unit, denoising；And the sound pressure signal memory cell before denoising, the impulse response before denoising are deposited The impulse response memory cell after sound pressure signal memory cell, denoising after storage unit, denoising；The acoustical parameters of also each frequency band Result of calculation；
Wherein, acoustic pressure time-domain signal is bent in sound receiving room when described time-domain signal display unit real-time display measures every time Line, described acoustic energy decay curve display unit with 63Hz, 80Hz, 100Hz, 125Hz, 160Hz, 200Hz, 250Hz, 315Hz, 400Hz、500Hz、630Hz、800Hz、1000Hz、1250Hz、1600Hz、2000Hz、2500Hz、3150Hz、4000Hz、 The 1/3 octave band acoustic energy decay curve of frequency centered on 5000Hz, 6300Hz, 8000Hz；Described impulse response display unit Show acoustic pressure and Energy impulse response curve map；Described sound pressure signal memory cell is stored in the acoustic pressure letter collected in room Number, described impulse response memory cell storage acoustic pressure pulse responds；When described acoustical parameters computing module is mainly reverberation Between；Described acoustical parameters result of calculation store with 63Hz, 80Hz, 100Hz, 125Hz, 160Hz, 200Hz, 250Hz, 315Hz, 400Hz、500Hz、630Hz、800Hz、1000Hz、1250Hz、1600Hz、2000Hz、2500Hz、3150Hz、4000Hz、 The 1/3 octave band reverberation time of frequency is stored with excel forms centered on 5000Hz, 6300Hz, 8000Hz.
Compared with prior art, the present invention has advantages below：
1st, sound field is encouraged using linear frequency sweep acoustical signal, arteries and veins is obtained by reception signal and Digital Sweep signal correlation computations Punching response；Linear chirp requires low to system time-varying, anti-distortion ability is strong, can ensure useful signal in multiple measurement It is almost consistent, and noise is because interference when its strong randomness to measure every time to signal is inconsistent, it is convenient to be surveyed using multiple The consistency of signal and the variability of noise carry out denoising to signal during amount.
2nd, signal to noise ratio determination module provides direct judgement to test condition, to whether needing repeatedly measurement, determining repeatedly to survey Amount number provides visual reference foundation.
3rd, the method that will repeatedly measure the signal progress time domain superposition collected, especially suitable for multiple measurement signal noise Than close to situation that is consistent, being disturbed compared with stationary noise.When repeatedly measurement signal to noise ratio is more steady close to consistent and noise characteristic When, meet between signal to noise ratio size SNR and duplicate measurements times N that denoising can be lifted：SNR=10log N, duplicate measurements time Number is more, and the signal to noise ratio that can be lifted is bigger.
4th, by multiple measurement signal carry out time-frequency coefficients reconstruct method, especially suitable for repeatedly measurement interference noise when The obvious situation of frequency property difference.If any 3 measurements, disturbed respectively by different types of noise every time, through time-frequency coefficients weight After structure, at least 10dB signal to noise ratio can be lifted.
By multiple measurement signal to be carried out to the denoising method of time domain superposition or time-frequency coefficients reconstruct, language can be effectively removed All kinds of stable states such as sound, traffic, equipment and construction, unstable state and impulsive noise, it can will measure necessary signal to noise ratio and be down to by 35dB 10dB, the precision and reliability of impulse response in-site measurement are significantly improved, reduce the difficulty of in-site measurement, survey crew need not Call powerful acoustical signal emitter (power of acoustical signal emitter can be reduced into 100 to obtain 35dB signal to noise ratio Times), it need not also wait until that the late into the night carries out measurement again.
The signals and associated noises feature provided when the 5th, for measurement every time in signal preanalysis, flexibly selects denoising method, 100 In the range of~8000Hz, the signal to noise ratio measured needed for impulse response can be decreased to 10dB, substantially increase in-site measurement has Effect property and reliability, reduce the requirement to acoustical signal emission system power, realize the portability of measuring system.
In the market both at home and abroad similar impulse response measurement system largely without the complete acoustical signal denoising system of use, It must be measured under conditions of signal to noise ratio is more than 35dB.This defect greatly limit during measuring system measures at the scene Using.
6th, there is data transmitting, collection, storage, Treatment Analysis integration, avoid artificial treatment in traditional measurement method Error caused by data and mistake, significantly improve the precision and efficiency of measurement.Denoising method can be according to actual test situation Flexibly selected, greatly improve measurement and treatment effeciency.
7th, denoising process and result visualization, tester can decline song by observing impulse response and each sound control The form of line can tentatively judge the validity of test and denoising result, it is found that wave chopping, signal to noise ratio deficiency, sound field diffusion level are poor The problem of causing false test result, in time adjustment test and denoising scheme, obtain correct impulse response test result.Original creation Display interface parameter setting, crucial denoising and measurement result are subjected to image conversion shown, make operator real according to image When judge whether denoising and measurement result correct, overcome both at home and abroad at present measuring system to signal carry out secret operation, cause The defects of tester has no way of judging measurement result reliability.
Brief description of the drawings
Fig. 1 is the structural representation of apparatus of the present invention；
Fig. 2 is that (A is ideal signal to signal preanalysis time-frequency characteristic analysis chart, and B is to be believed by the test after white noise acoustic jamming Number)；
Fig. 3 is signal of change acoustic energy decay curve figure after time domain superposition；
Fig. 4 is that (A is ideal signal to signal preanalysis time-frequency distributions performance plot, and B is voice noise interference signal, and C is pulse Noise interferences, D are traffic noise interference signal)；
Fig. 5 is signal of change acoustic energy decay curve figure after time-frequency coefficients reconstruct.
Embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
Embodiment
As shown in figure 1, the device of the noise remove based on room impulse response duplicate measurements, it is characterised in that including letter Number transmitter module 1, signal acquisition module 2, denoising prepulse response computing module 3, signal to noise ratio determination module profit 4, pendulous frequency Determining module 5, signal pre-analysis module 6, superposition denoising module 7, coefficient reconstruct denoising module 8, the response of denoising afterpulse calculate Module 9, acoustical parameters computing module 10, result are shown and memory module 11.
Described signal emission module 1 passes through data signal generation unit, digital-to-analogue conversion card, power amplifier and loudspeaker Cooperate, Linear chirp is sent to room；Described signal acquisition module 2 passes through the microphone that is sequentially connected, preposition Amplifier, analog-to-digital conversion card and digital signal receiving unit gather the sound pressure signal in room, and are passed along arteries and veins before denoising Punching response computing module 3；The signal collected and the data signal of swept-frequency signal are made phase by denoising prepulse response computing module Computing is closed, obtains the impulse response in room, and be transferred to signal to noise ratio determination module 4；Signal to noise ratio determination module utilizes denoising Prepulse response computing module calculates signal to noise ratio, judges whether that needs take multiple measurements according to signal to noise ratio, if signal to noise ratio is big It is small it is accurate enough calculate acoustical parameters, then send the response of denoising prepulse to acoustical parameters computing module 10, such as if not satisfied, Then send signal to noise ratio to pendulous frequency determining module 5；Pendulous frequency determining module carries out n times weight according to determination pendulous frequency N Repetition measurement amount, and the signal that collects will be measured every time and send signal pre-analysis module 6 to；Signal pre-analysis module is to measurement every time The signal collected carries out preanalysis, provides its time-frequency representation, and denoising method is determined according to the time-frequency characteristic of signal, if selection Time domain superposition carries out denoising to signal, then the signal by multiple measurement collection sends superposition denoising module 7 to；If select time-frequency The method of coefficient reconstruct carries out denoising to signal, then the signal by multiple measurement collection sends reconstruct denoising module 8 to；So Send the signal of denoising to denoising afterpulse response computing module 9 afterwards；Denoising afterpulse response computing module, which utilizes, to be gone The data signal of signal of change and swept-frequency signal after making an uproar makees related operation, obtains room impulse response after denoising, is then transported on To acoustical parameters computing module 10；Acoustical parameters computing module calculates acoustical parameters using the forward and backward impulse response of denoising, concurrently Result is given to show and memory module 11；As a result display and memory module show measure in a room sound pressure signal curve, It is acoustic pressure pulse response curve after sound pressure signal curve, original acoustic pressure pulse response curve, denoising after denoising, original Energy impulse response curve, denoising afterpulse response curve, acoustical parameters result of calculation, and store its data signal and tied with calculating Fruit.
The present invention use step be：
1) parameter setting：Selection signal transmission channel and acquisition channel；A length of 10s, collection signal when inputting transmission signal Shi Changwei 20s, sample frequency are 44100Hz (should be greater than at least 2 times of highest cut-off frequency), signal amplitude 1, swept-frequency signal Initial frequency is 20Hz and cut-off frequency is 20kHz, and signal pre-analysis module analysis window function is taken as Gaussian window, window function width 2048 are taken as, time-domain analysis step-length is 16, i.e. time-domain resolution ratio is 16/44100=3.9 × 10-3s；Frequency-domain analysis step-length is 1024, i.e. frequency domain resolution is 44100/1024=43Hz；
2) experiment is carried out in acoustic investigation institute of Tongji University reverberation chamber, by laying certain sound absorption in reverberation chamber Body comes analog functuion room, and wherein 1000Hz reverberation time is 2.4s.According to the noise characteristic in actual test condition, divide Stationary white noise and nonstationary noise (voice noise, traffic noise and impulsive noise) are disturbed sounding trial signal.
3) when noiseless is disturbed, signal to noise ratio size calculates acoustical parameters enough, therefore without repeatedly measurement, directly calculate Acoustical parameters simultaneously preserve result.
4) when interference noise is stationary white noise, signal to noise ratio is about 10dB, is comparing result, pendulous frequency is distinguished It is set to 5 times, 10 times, 20 times；Signal preanalysis characteristic shows that the size of useful signal time-frequency conversion coefficient has with time and frequency Obvious variation tendency, time-frequency characteristic is notable, and the conversion coefficient of stationary white noise is evenly distributed on time frequency analysis figure, such as Fig. 2 It is shown.The method of time domain superposition is selected to carry out denoising to signal, result of calculation is as shown in Figure 3 and Table 1.As a result show：When Domain stacking fold is more, and signal to noise ratio lifting is bigger, and substantially meets SNR=10log N rule.Generalized time and efficiency, will Pendulous frequency is set to 10.
Signal of change reverberation time and its relative error after the superposition of the time domain of table 1
4) when interference noise is that nonstationary noise (tests 3 times respectively with voice noise, traffic noise and impulsive noise to letter Number disturbed), signal to noise ratio substantially 10dB, pendulous frequency is set to 3, found by signal preanalysis, each measurement signal Time-frequency characteristic significant difference, so, using time-frequency coefficients reconstruct method to signal carry out denoising.
Be utilized respectively each signal of change acoustic energy decay curve and calculate this acoustical parameters of reverberation time, as a result such as Fig. 5 and Shown in table 2.
Signal of change reverberation time and its relative error after the reconstruct of the time-frequency coefficients of table 2
It is above-mentioned test result indicates that：The sounding trial signal either disturbed by stationary noise interference or nonstationary noise, Method by the way that multiple measurement signal is carried out into time domain superposition or time-frequency coefficients reconstruct, is attained by preferable denoising effect.Say This bright denoising method has stronger adaptability and practicality with device.

Claims (17)

  1. A kind of 1. method of the noise remove based on room impulse response duplicate measurements, it is characterised in that comprise the following steps：
    1) signal emission module by the data signal generation unit that is sequentially connected, digital-to-analogue conversion card, power amplifier and is raised one's voice Device sends linear frequency sweep acoustical signal into room；
    2) signal acquisition module passes through the microphone, preamplifier, analog-to-digital conversion card and the digital signal reception list that are sequentially connected Member, the sound pressure signal at room measuring point is gathered, and send it to denoising prepulse response computing module；
    3) signal and the data signal of swept-frequency signal that collect are made computing by denoising prepulse response computing module, obtain room Impulse response, and be transferred to signal to noise ratio determination module；
    4) signal to noise ratio under initial condition is calculated using the response of denoising prepulse for signal to noise ratio determination module, judges whether it reaches The snr threshold calculated to acoustical parameters, if it is satisfied, then send impulse response to acoustical parameters computing module, it is such as discontented Foot, then send signal to noise ratio to pendulous frequency determining module；
    5) pendulous frequency determining module determines duplicate measurements number according to signal to noise ratio condition；
    6) signal pre-analysis module carries out time-frequency characteristic analysis to the signal collected, and selects noise-reduction method accordingly：Superposition is gone Make an uproar or time-frequency coefficients reconstruct denoising；
    7) signal repeatedly measured is overlapped by superposition denoising module, is realized noise suppressed and useful signal enhancing, is gone Signal after making an uproar, and denoising afterpulse response computing module is sent to, and perform step 9)；
    8) signal repeatedly measured is carried out time-frequency coefficients reconstruct by coefficient reconstruct denoising module, is realized noise remove in environment, is obtained Signal after to denoising, and denoising afterpulse response computing module is sent to, and perform step 9)；
    9) sound pressure signal after denoising and the data signal of swept-frequency signal are made computing by denoising afterpulse response computing module, are obtained Room impulse response after denoising, and it is sent to acoustical parameters computing module；
    10) acoustical parameters computing module calculates acoustical parameters using the forward and backward impulse response of denoising, and be sent to result show and Memory module；
    11) result show and memory module to show that sound pressure signal curve before and after denoising, pulse respond, sound pressure level decline bent Line, acoustical parameters result of calculation, and store its data signal and result of calculation.
  2. 2. a kind of method of noise remove based on room impulse response duplicate measurements according to claim 1, its feature It is, described data signal is generated as a data signal generation unit, described data signal generation unit generation signal Calculating process is as follows：
    Wherein, A is signal amplitude；T is signal total duration；T is time variable；f1It is off lower-frequency limit；f2It is off in frequency Limit.
  3. 3. a kind of method of noise remove based on room impulse response duplicate measurements according to claim 1, its feature It is, described signal acquisition part uses triggering collection, ensures that repeatedly measure the signal collected corresponds point by point, and will The signal collected is sent to denoising prepulse response computing module.
  4. 4. a kind of method of noise remove based on room impulse response duplicate measurements according to claim 1, its feature It is, described denoising prepulse response computing module is to utilize the digital linear swept-frequency signal e modulated in signal emission module0 (t) the sound pressure signal x (t) collected with measuring point is calculated, and obtains noisy acoustic pressure impulse response p (t) of the measuring point on room, And it is transferred to signal to noise ratio determination module.
  5. 5. a kind of method of noise remove based on room impulse response duplicate measurements according to claim 1, its feature It is, described signal to noise ratio determination module calculates signal to noise ratio using impulse response, determines whether to reach letter according to signal to noise ratio size Make an uproar and compare threshold value：I.e. if signal to noise ratio is more than 25dB, then acoustical parameters computing module is sent to, is such as unsatisfactory for, be then sent to measurement Number determining module.
  6. 6. a kind of method of noise remove based on room impulse response duplicate measurements according to claim 1, its feature It is, described pendulous frequency determining module, duplicate measurements number is determined according to signal to noise ratio size：Signal to noise ratio is 15-25dB weights Repetition measurement amount 5 times, the duplicate measurements that signal to noise ratio is 5-15dB 10 times.
  7. 7. a kind of method of noise remove based on room impulse response duplicate measurements according to claim 1, its feature It is, described signal preanalysis is to carry out VARIANT FILTER OF DISCRETE GABOR TRANSFORM to the signal collected：
    Wherein, W=exp (j2 π/N)；It is the amplitude of time frequency analysis coefficient；It is the acoustic pressure time-domain signal of collection；ΔM It is time-domain analysis step-length；Δ N is frequency-domain analysis step-length；It is the analytic function of Gabor transformation；It is pre- by signal Analysis is obtained gathering signal time frequency analysis, and intuitive judgment is provided to noise characteristic.
  8. 8. a kind of method of noise remove based on room impulse response duplicate measurements according to claim 1, its feature It is, described superposition denoising module, for the signal repeatedly measured to be carried out into time domain superposition, realizes useful signal enhancing and nothing Noise suppressed is advised, the signal after processing is then sent to denoising afterpulse response computing module.
  9. 9. a kind of method of noise remove based on room impulse response duplicate measurements according to claim 1, its feature It is, described coefficient reconstruct denoising module, for the signal repeatedly measured to be compared, for identical time-frequency point selection Minimum time-frequency coefficients, the signal after then time-frequency coefficients are reconstructed returns to time domain, and is sent to denoising afterpulse response meter Calculate module.
  10. 10. a kind of method of noise remove based on room impulse response duplicate measurements according to claim 1, its feature It is, described denoising afterpulse response computing module is to utilize the digital linear swept-frequency signal e modulated in signal emission module0 (t) make correlation computations with the sound pressure signal y ' (t) after denoising, obtain the room impulse response after denoising, and transmitted Give acoustical parameters computing module.
  11. 11. a kind of method of noise remove based on room impulse response duplicate measurements according to claim 1, its feature It is, described acoustical parameters computing module utilizes the impulse response before and after denoising to calculate the reverberation time (RT60), early stage decay when Between (EDT), the definition factor (D50), the perspectivity factor (C50), the loudness factor (G), the hall sound of subjective assessment (IACC) Matter parameter.
  12. 12. the device of the noise remove according to claim 1 based on room impulse response duplicate measurements, its feature exist In, including signal emission module, signal acquisition module, denoising prepulse response computing module, signal to noise ratio determination module, measurement time Number determining module, superposition denoising module, time-frequency coefficients reconstruct denoising module, denoising afterpulse response computing module, acoustical parameters Computing module and result are shown and memory module.
  13. 13. device according to claim 12, it is characterised in that described signal emission module includes the number being sequentially connected Word signal generation and weighted units, digital-to-analogue conversion card, power amplifier and loudspeaker, for completing to sweep to hall internal emission line The task of frequency acoustical signal.
  14. 14. device according to claim 12, it is characterised in that described signal acquisition module includes the biography being sequentially connected Sound device, preamplifier, analog-to-digital conversion card and digital signal receiving unit, to complete the task that sound pressure signal gathers in room, And collection signal is sent to signal pre-analysis module and denoising prepulse response computing module.
  15. 15. device according to claim 12, it is characterised in that selection signal transmission channel and acquisition channel；Input hair Penetrate signal duration, collection signal duration, sample frequency, signal amplitude, the initial frequency of swept-frequency signal and cut-off frequency.
  16. 16. device according to claim 12, it is characterised in that window function type, window function are selected during signal preanalysis Width, time sampling step-length, frequency sampling step-length.
  17. 17. device according to claim 12, it is characterised in that described result display module includes what display was gathered The impulse response display unit before time-domain signal display unit, acoustic energy decay curve display unit, denoising before denoising, after denoising Time-domain signal display unit, acoustic energy decay curve display unit, the impulse response display unit after denoising；And before denoising The sound pressure signal memory cell after impulse response memory cell, denoising before sound pressure signal memory cell, denoising, the arteries and veins after denoising Punching response memory cell；The acoustical parameters result of calculation of also each frequency band；
    Wherein, acoustic pressure time-domain signal curve in room, described when described time-domain signal display unit real-time display measures every time Acoustic energy decay curve display unit with 63Hz, 80Hz, 100Hz, 125Hz, 160Hz, 200Hz, 250Hz, 315Hz, 400Hz, 500Hz、630Hz、800Hz、1000Hz、1250Hz、1600Hz、2000Hz、2500Hz、3150Hz、4000Hz、5000Hz、 The 1/3 octave band acoustic energy decay curve of frequency centered on 6300Hz, 8000Hz；Described impulse response display unit shows acoustic pressure And Energy impulse response curve map；The sound pressure signal that described sound pressure signal memory cell storage collects indoors, it is described The storage acoustic pressure pulse response of impulse response memory cell；Described acoustical parameters computing module is mainly the reverberation time；Described Acoustical parameters result of calculation storage with 63Hz, 80Hz, 100Hz, 125Hz, 160Hz, 200Hz, 250Hz, 315Hz, 400Hz, 500Hz、630Hz、800Hz、1000Hz、1250Hz、1600Hz、2000Hz、2500Hz、3150Hz、4000Hz、5000Hz、 The 1/3 octave band reverberation time of frequency is stored with excel forms centered on 6300Hz, 8000Hz.
CN201610726044.1A 2016-08-25 2016-08-25 Noise removal method and device based on repeated measurement of room impulse response Active CN107785025B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610726044.1A CN107785025B (en) 2016-08-25 2016-08-25 Noise removal method and device based on repeated measurement of room impulse response

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610726044.1A CN107785025B (en) 2016-08-25 2016-08-25 Noise removal method and device based on repeated measurement of room impulse response

Publications (2)

Publication Number Publication Date
CN107785025A true CN107785025A (en) 2018-03-09
CN107785025B CN107785025B (en) 2021-06-22

Family

ID=61439660

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610726044.1A Active CN107785025B (en) 2016-08-25 2016-08-25 Noise removal method and device based on repeated measurement of room impulse response

Country Status (1)

Country Link
CN (1) CN107785025B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109545242A (en) * 2018-12-07 2019-03-29 广州势必可赢网络科技有限公司 A kind of audio data processing method, system, device and readable storage medium storing program for executing
CN110398338A (en) * 2018-04-24 2019-11-01 广州汽车集团股份有限公司 Wind is obtained in wind tunnel test to make an uproar the method and system of speech intelligibility contribution amount
CN111535470A (en) * 2020-05-09 2020-08-14 广州大学 Method and system for improving indoor language definition
CN115862665A (en) * 2023-02-27 2023-03-28 广州市迪声音响有限公司 Visual curve interface system of echo reverberation effect parameters

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101109731A (en) * 2007-08-08 2008-01-23 哈尔滨工业大学 Gabor translating self-adapting window width selecting method represented by ultrasound signal
CN101510426A (en) * 2009-03-23 2009-08-19 北京中星微电子有限公司 Method and system for eliminating noise
CN101620004A (en) * 2008-07-04 2010-01-06 中国科学院声学研究所 Method for measuring direct wave signal sound pressure of sound wave transmitted in limited regional medium
US20100322437A1 (en) * 2009-06-23 2010-12-23 Fujitsu Limited Signal processing apparatus and signal processing method
CN102801861A (en) * 2012-08-07 2012-11-28 歌尔声学股份有限公司 Voice enhancing method and device applied to cell phone
CN102918588A (en) * 2010-03-29 2013-02-06 弗兰霍菲尔运输应用研究公司 A spatial audio processor and a method for providing spatial parameters based on an acoustic input signal
CN102928073A (en) * 2012-11-26 2013-02-13 昆山北极光电子科技有限公司 Phase-congruency-based vibration noise filtering method
CN103414845A (en) * 2013-07-24 2013-11-27 中国航天科工集团第三研究院第八三五七研究所 Self-adaptive video image noise reducing method and noise reducing system
CN103929706A (en) * 2013-01-11 2014-07-16 克里佩尔有限公司 Arrangement and method for measuring the direct sound radiated by acoustical sources
CN103994820A (en) * 2014-04-21 2014-08-20 中国科学院上海微系统与信息技术研究所 Moving target identification method based on micro-aperture microphone array
CN104363545A (en) * 2014-10-17 2015-02-18 深圳市时尚德源文化传播有限公司 Audio denoising circuit, smart terminal adopting same and teaching method
US20150211917A1 (en) * 2010-01-03 2015-07-30 Acousticeye Ltd Method and system for improving mls-based measurements
CN105371945A (en) * 2015-10-10 2016-03-02 同济大学 Method and device for removing noise in hall pulse response measurement sound signals

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101109731A (en) * 2007-08-08 2008-01-23 哈尔滨工业大学 Gabor translating self-adapting window width selecting method represented by ultrasound signal
CN101620004A (en) * 2008-07-04 2010-01-06 中国科学院声学研究所 Method for measuring direct wave signal sound pressure of sound wave transmitted in limited regional medium
CN101510426A (en) * 2009-03-23 2009-08-19 北京中星微电子有限公司 Method and system for eliminating noise
US20100322437A1 (en) * 2009-06-23 2010-12-23 Fujitsu Limited Signal processing apparatus and signal processing method
US20150211917A1 (en) * 2010-01-03 2015-07-30 Acousticeye Ltd Method and system for improving mls-based measurements
CN102918588A (en) * 2010-03-29 2013-02-06 弗兰霍菲尔运输应用研究公司 A spatial audio processor and a method for providing spatial parameters based on an acoustic input signal
CN102801861A (en) * 2012-08-07 2012-11-28 歌尔声学股份有限公司 Voice enhancing method and device applied to cell phone
CN102928073A (en) * 2012-11-26 2013-02-13 昆山北极光电子科技有限公司 Phase-congruency-based vibration noise filtering method
CN103929706A (en) * 2013-01-11 2014-07-16 克里佩尔有限公司 Arrangement and method for measuring the direct sound radiated by acoustical sources
CN103414845A (en) * 2013-07-24 2013-11-27 中国航天科工集团第三研究院第八三五七研究所 Self-adaptive video image noise reducing method and noise reducing system
CN103994820A (en) * 2014-04-21 2014-08-20 中国科学院上海微系统与信息技术研究所 Moving target identification method based on micro-aperture microphone array
CN104363545A (en) * 2014-10-17 2015-02-18 深圳市时尚德源文化传播有限公司 Audio denoising circuit, smart terminal adopting same and teaching method
CN105371945A (en) * 2015-10-10 2016-03-02 同济大学 Method and device for removing noise in hall pulse response measurement sound signals

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
HOUBIAO LIU: "Adaptive Iterative Truncated Arithmetic Mean Filter in Image Denoising", 《2012 IEEE 12TH INTERNATIONAL CONFERENCE ON COMPUTER AND INFORMATION TECHNOLOGY》 *
姜春香: "金属矿地震资料降噪处理技术研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 *
杨春庄等: "混响时间测量虚拟仪器研究", 《声学技术》 *
王海玲等: "基于时频分析的声信号去噪", 《声学技术》 *
蒋毅: "基于听感知特性的噪声分析与防护", 《电声技术》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110398338A (en) * 2018-04-24 2019-11-01 广州汽车集团股份有限公司 Wind is obtained in wind tunnel test to make an uproar the method and system of speech intelligibility contribution amount
CN109545242A (en) * 2018-12-07 2019-03-29 广州势必可赢网络科技有限公司 A kind of audio data processing method, system, device and readable storage medium storing program for executing
CN111535470A (en) * 2020-05-09 2020-08-14 广州大学 Method and system for improving indoor language definition
CN111535470B (en) * 2020-05-09 2021-11-19 广州大学 Method and system for improving indoor language definition
CN115862665A (en) * 2023-02-27 2023-03-28 广州市迪声音响有限公司 Visual curve interface system of echo reverberation effect parameters

Also Published As

Publication number Publication date
CN107785025B (en) 2021-06-22

Similar Documents

Publication Publication Date Title
Garai Measurement of the sound-absorption coefficient in situ: the reflection method using periodic pseudo-random sequences of maximum length
CN105371945B (en) Remove the method and device of noise in hall impulse response measurement acoustical signal
CN106768266B (en) A kind of measurement method of workshop reverberation time that rejecting background noise energy
CN101194535B (en) Method for correcting electroacoustic converter acoustic paramenter and device accomplishing the method
CN107785025A (en) Noise remove method and device based on room impulse response duplicate measurements
CN104569154B (en) The detection method and device of quick nondestructive fruit quality
Zhang et al. An improved AE detection method of rail defect based on multi-level ANC with VSS-LMS
CN110186546B (en) Hydrophone sensitivity free field broadband calibration method based on pink noise
CN103675104B (en) The measuring method of material random incidence acoustical absorption coefficient or absorption and measurement mechanism thereof
CN109238436B (en) Method for measuring transient sound source characteristics in reverberation pool
CN104535647A (en) Prediction apparatus for sound absorption and insulation performance of multilayer material and method
CN103698403A (en) Measuring method and measuring device for sound insulation quantities of material and members
CN106031197A (en) Acoustic processing device, acoustic processing method, and acoustic processing program
US9060222B2 (en) Method for determining an averaged frequency-dependent transmission function for a disturbed linear time-invariant system, evaluation device and computer program product
CN104869519B (en) A kind of method and system of test microphone background noise
CN108718438B (en) Method for judging time starting point of response signal of electroacoustic product detection system
Duroux et al. Extracting guided waves from cross-correlations of elastic diffuse fields: Applications to remote structural health monitoring
Simmons Noise interference with echo delay discrimination in bat biosonar
Nutter et al. Measurement of sound power and absorption in reverberation chambers using energy density
Vergara et al. Measurement of cement porosity by centroid frequency profiles of ultrasonic grain noise
Olsson et al. Incremental use of FFT as a solution for low BT-product reverberation time measurements
CN103050127A (en) Hall objective audio quality parameter measuring system and method
CN115549813A (en) Sound compatibility lake test method of communication underwater sound equipment
JP2019032242A (en) Attenuation time analytic method, device, and program
CN110988130A (en) Test system for rock slope indoor test damage identification

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