Patent Application
20170098454
This application is directed to a method and system for producing low-noise acoustical impulse responses and more particularly to a method and system for producing low-noise acoustical impulse responses at high sampling rates for use in acoustical systems.
A well-established approach for measuring impulse responses (IRs) of acoustical systems at sampling rates of either 44.1 kHz or 48 kHz involves the use of an exponential sine sweep (ESS) up to the Nyquist frequency. However, some acoustical systems have responses that extend beyond 24 kHz, and using this conventional ESS approach presents problems.
If the conventional ESS approach is used to measure a system whose response extends beyond the Nyquist frequency, the anti-aliasing filter of the A/D converter will low-pass filter the system response with a cut-off near the Nyquist frequency, eliminating the response that exists above this point. This would manifest itself as a spurious (non-causal) pre-response in the measured IR. Therefore, in order to accurately characterize the entire system response, it is necessary to adopt higher sampling rates. Doing so would eliminate the pre-response, but poses the risk of damaging the transducers when attempting a sweep to a higher Nyquist frequency. This is especially true when using long-duration sweeps as is done in existing methods. Additionally, as the system response may now fall below the noise floor prior to the Nyquist frequency, high-frequency noise can potentially contaminate the measurement and yield a less than maximal signal-to-noise ratio.
Abrupt termination of the sweep before the Nyquist frequency introduces artifacts such as an end-of-sweep “pop” in the transducer that corrupts the measurement and may damage the transducers. This issue can be prevented by applying a fade-out to the end of the sweep as is done in existing methods. However, this solution is only viable when employing the time-reversed sweep inversion approach derived by Angelo Farina and described in the article entitled “Simultaneous Measurement of Impulse Response and Distortion with a Swept-Sine Technique,” presented at the AES 108th Convention, February 2000, since the exact, frequency-domain inverse of the faded-out sweep will cause excessive high-frequency noise amplification. However, the time-reversed sweep inversion approach produces a pre-response if the sweep does not cover the entire system response.
A critical step in such impulse response interpolations is locating impulse response onsets for time-alignment, a task which becomes difficult in the presence of pre-responses. Furthermore, it has been found that pre-responses can be audible as described by Peter G. Craven in “Antialias Filters and System Transient Response at High Sample Rates,” J. Audio Eng. Soc., 52(3):216-242 (2004), and so their suppression is desirable.
It is therefore an object of the present invention to provide a method and system of measuring low-noise acoustical impulse responses at high sampling rates that enable low noise measurement without non-causal pre-response contamination.
The method and system for measuring low-noise acoustical impulse responses at high sampling rates of the present invention utilizes two exponential sine sweeps (ESSs) to measure the impulse responses. The first ESS is a quick sweep up to the Nyquist frequency to provide an estimate of the system response and sample the ambient noise. This measurement is used to algorithmically determine an appropriate pass-band of the system. A second, slower sweep through the pass-band alone is then executed and a corresponding band-pass filter is applied to the resulting output signal to suppress noise. The result is a measured impulse response with an improved signal-to-noise ratio and a much-reduced pre-response.
In general, the method and system for measuring low-noise acoustical impulse responses at high sampling rates of the present invention carries out a quick (approximately 1 second or shorter) initial ESS up to the Nyquist frequency to provide an estimate of the system response and sample the ambient noise. This measurement is used to algorithmically determine an appropriate pass-band of the system. A second, slower sweep through the pass-band alone is then executed and a corresponding band-pass filter is applied to the resulting output signal to suppress noise. The result is a measured impulse response with an improved signal-to-noise ratio and a much-reduced pre-response.
Referring to
Following the first sweep, in step 14, an “optimal-SNR” pass-band is determined by finding the widest possible frequency range over which the signal-to-noise ratio (SNR) is positive. For this pass-band, a measure of the effects of filtering the system response with a corresponding band-pass filter is estimated. The measure of interest is called the peak deviation amplitude (PDA), which quantifies the largest possible pre-response amplitude that will result from band-pass filtering over the previously found pass-band. The PDA is estimated based on the measured frequency response magnitude at the pass-band edge frequencies using models derived by Joseph G. Tylka et al. in the article “A New Approach to Impulse Response Measurements at High Sampling Rates,” presented at the AES 137th Convention, October 2014.
In step 16, the user specifies as to whether the estimated PDA value is acceptable. For example, the user (which can be a user of the method or any system that relies on the method) may decide that the PDA must be less than or equal to 1% of the impulse response's peak amplitude. If this upper-limit is not exceeded, then the user selects to accept the predicted PDA value. Alternatively, in an automatic version of the method, the minimum acceptable PDA is specified a priori. In step 18, the user decides to accept or not accept the predicted PDA depending on the requirements of the particular application the method is used for. If the user elects to accept the predicted PDA value, then, in step 24, a refined, band-limited ESS that sweeps through the optimal-SNR pass-band alone is designed and executed, including a fade-out to prevent an end-of-sweep “pop”. Ideally, this second sweep is as long as possible considering the environment in order to increase the signal-to-noise ratio. In one preferred embodiment, the second sweep is in the range of 5-10 seconds. The microphone signal is then band-pass filtered in step 26 and deconvolved by the input sweep in step 28 to get the final impulse response. At that point the process is concluded in step 30.
If the estimated PDA value is determined to be not acceptable in step 18, then the user is asked to specify a maximum tolerable PDA value in step 20. A new “constrained-PDA” pass band is now found in step 22 to satisfy the specified maximum. Using the same PDA models as in step 14, the lowest low-edge frequency that results in one half of the maximum PDA value is found. Similarly, the highest high-edge frequency that results in one half of the maximum PDA value is found. Therefore, the net effect of applying the corresponding “constrained-PDA” band-pass filter will not result in a PDA value that exceeds the maximum set by the user. In step 24, a refined, band-limited ESS that sweeps through the constrained-PDA pass-band alone is designed and executed, including a fade-out to prevent an end-of-sweep “pop”. The second sweep here is preferably as long as possible considering the environment in order to increase the signal-to-noise ratio. The microphone signal is then band-pass filtered in step 26 and deconvolved by the input sweep in step 28 to get the final impulse response. At that point the process is concluded in step 30.
Referring to
The method and system of the present invention described above can be used in many applications that require high-fidelity impulse response measurements at high sampling rates. Such applications include high-resolution 3D headphones processors, high-resolution audio components that rely on calibration though impulse response measurements, and instruments for characterizing acoustical systems at high resolution such as transducers, acoustical spaces and for measuring head-related impulse responses (HRIR). The method and system of the present invention could also be used in other commercial and industrial products that rely on, or require, low-noise high sample rate measurements such as virtual reality audio systems and other non-audio applications where low-noise, high sample rate impulse measurements are needed or required.
The method and system of the present invention also allow low-noise measurements, free of non-causal pre-response contamination, of acoustical impulse responses at high sampling rates. As mentioned above, an important application of this method and system is its use in 3D headphones processors where interpolation between two or more impulse responses is required to apply the appropriate digital filter as a function of the tracked head rotation coordinate (yaw angle) in order to fix the perceived audio image in 3D space so that that the image does not move with the listener's head as such motion severely degrades the ability of a listener to head externalize the sound.
By using high sampling rates of 96 kHz or more, the method of the present invention reduces the effects of time smearing, which may otherwise prevent the capture of the minimum perceptible inter-aural time difference, a vital spatialization cue for accurate binaural audio.
The only relative disadvantage the new method has over the prior art methods is the requirement of two sweeps. This disadvantage is a minor issue as the sweeps can be easily automated in a processor using pre-defined user preferences, and the benefits justify the additional sweep.
While the foregoing invention has been described with reference to its preferred embodiments, various alterations and modifications will occur to those skilled in the art. All such variations and modifications are intended to fall within the scope of the appended claims.
This application relates to and claims priority under 35 U.S.C. §119(e) to U.S. Provisional Patent Application No. 62/237,739, titled “METHOD AND SYSTEM FOR PRODUCING LOW-NOISE ACOUSTICAL IMPULSE RESPONSES AT HIGH SAMPLING RATE,” which was filed on Oct. 6, 2015 and is hereby incorporated by reference and is hereby incorporated by reference herein in its entirety.
| Number | Date | Country | |
|---|---|---|---|
| 62237739 | Oct 2015 | US |
