Active room compensation in loudspeaker system

• US 10,448,187 B2
• Filed: 12/16/2015
• Issued: 10/15/2019
• Est. Priority Date: 10/08/2015
• Status: Active Grant

• Current Assignee: Bang & Olufsen a/s
• Original Assignee: Bang & Olufsen a/s
• Inventors: Dyreby, Jakob
• Primary Examiner(s): Paul, Disler
• Application Number:
  US15/757,927
• Publication Number:
  US 20180249272A1
• Time in Patent Office: 1,399 Days
• CPC Class Codes: H04R 3/04 for correcting frequency response
• CPC Class Codes: H04R 3/12 for distributing signals to two or more loudspeakers specially adapted for hearing aids H04R25/407
• CPC Class Codes: H04R 5/02 Spatial or constructional arrangements of loudspeakers
• CPC Class Codes: H04R 5/04 Circuit arrangements, e.g. for selective connection of amplifier inputs/outputs to loudspeakers, for loudspeaker detection, or for adaptation of settings to personal preferences or hearing impairments combinations of amplifiers H03F3/68; stereophonic systems H04S
• CPC Class Codes: H04S 7/301 Automatic calibration of stereophonic sound system, e.g. with test microphone
• CPC Class Codes: H04S 7/303 Tracking of listener position or orientation

Abstract

A method for compensating for acoustic influence of a listening room on an acoustic output from an audio system including at least a left and a right loudspeaker, the method comprising determining a left frequency response and a right frequency response, designing a left compensation filter F_L, and a right compensation filter F_R, and during playback applying the left and right filters to left and right channel inputs. According to the invention, a target response in the listening position is simulated, and the left and right compensation filters are designed to filter transfer functions based on the simulated target function multiplied by an inverse of the left/right frequency responses. By relying on a simulated target instead of relying on an empirical approach, the general impact of a room can be more accurately captured by the target functions.

Claims

1. A method for compensating for acoustic influence of a listening room on an acoustic output from an audio system including at least a left and a right loudspeaker, the method comprising:

determining a left frequency response LP_L between a signal applied to the left speaker and a resulting power average in a listening position,

determining a right frequency response LP_R between a signal applied to the right speaker and a resulting power average in the listening position,

designing a left compensation filter F_L,

designing a right compensation filter F_R,

during playback, applying the left compensation filter to a left input signal, and applying the right compensation filter to a right input signal,

providing a simulated target function H_T representing a simulated target response in the listening position,

designing the left compensation filter F_L to have a left filter transfer function based on the simulated target function H_T multiplied by an inverse of the left frequency response, and

designing the right compensation filter F_R to have a right filter transfer function based on the simulated target function H_T multiplied by an inverse of the right frequency response,

wherein the simulated target function is obtained by simulating the power emitted by a point source in a corner defined by three orthogonal walls into a one eighth sphere limited by the three walls, and defining the simulated target function as the transfer function between the point source and the emitted power,

wherein the simulated emitted power is a power average based on simulations in a plurality of points distributed on the one eighth sphere.

2. The method according to claim 1, wherein a radius of the one eighth sphere is based on size of listening room.

3. The method according to claim 2, wherein the radius of the one eighth sphere is in the range 2-8 m.

4. The method according to claim 1, wherein:

determining the left and right frequency responses involves measuring sound pressure in the listening position and in two complementary positions located in opposite corners of a rectangular cuboid having a center point in the listening position, said rectangular cuboid being aligned with a line of symmetry between the left and right speakers, and

forming an average sound pressure from the measured sound pressures.

5. The method according to claim 1, further comprising removing dips in at least one response, by:

providing a reference by smoothing the response with a reference smoothing width,

comparing the response and the reference, and

for each frequency, selecting the maximum of the response and the reference as dip removed response.

6. The method according to claim 5, wherein the reference smoothing width is at least two octaves.

7. The method according to claim 5, wherein the response is smoothed before the comparing step using a smoothing width narrower than the reference smoothing width.

8. A method according to claim 1, further comprising:

smoothing a frequency response between a signal applied to a speaker and a resulting power average in a listening position using a method comprising;

determining a number of peaks per octave in the response,

for a portion of the response where the number of peaks per octave is below a first threshold, smoothing the response with a first smoothing width,

for a portion of the response where the number of peaks per octave is above a second threshold, smoothing the response with a second smoothing width,

wherein said second threshold is greater than said first threshold and said second smoothing width is wider than said first smoothing width, and

for a portion of the response where the number of peaks per octave is between the first and second thresholds, smoothing with an intermediate smoothing width.

9. The method according to claim 8, wherein the intermediate smoothing width is frequency dependent as an interpolation of the first and second smoothing width.

10. The method according to claim 8, wherein the first, narrow smoothing width is less than ¼

octave, preferable 1/12 octave, and the second, wide smoothing width is at least one octave.

11. The method according to claim 8, wherein the first, smaller threshold is less than eight peaks per octave, preferably five peaks per octave, and the second, greater threshold is greater than eight peaks per octave, preferably ten peaks per octave.

12. The method according to claim 8, further comprising:

providing a reference by smoothing the response with a reference smoothing width, wherein the reference smoothing width is wider than the second, wide smoothing width,

comparing the smoothed response and the reference, and

for each frequency, selecting the maximum of the smoothed response and the reference as dip removed response.

13. The method according to claim 12, wherein the reference smoothing width is at least two octaves.

14. The method of claim 1, wherein the simulated emitted power is the power average based on simulations in the plurality of points, including more than 12 points, distributed on the one eighth sphere.

15. A method for compensating for acoustic influence of a listening room on an acoustic output from an audio system including at least a left and a right loudspeaker, the method comprising:

determining a left frequency response LP_L between a signal applied to the left speaker and a resulting power average in a listening position,

determining a right frequency response LP_R between a signal applied to the right speaker and a resulting power average in the listening position,

designing a left compensation filter F_L,

designing a right compensation filter F_R,

during playback, applying the left compensation filter to a left input signal, and applying the right compensation filter to a right input signal,

providing a simulated target function H_T representing a simulated target response in the listening position,

designing the left compensation filter F_L to have a left filter transfer function based on the simulated target function H_T multiplied by an inverse of the left frequency response, and

designing the right compensation filter F_R to have a right filter transfer function based on the simulated target function H_T multiplied by an inverse of the right frequency response,

the method further comprising;

determining a left roll-off frequency at which the left target function exceeds the left frequency response by a given threshold,

determining a right roll-off frequency at which the left target function exceeds the right frequency response by a given threshold,

calculating an average roll-off frequency based on the left and right roll-off frequencies,

estimating a roll-off function as a high pass filter with a cut-off frequency based on the average roll-off frequency, and

dividing the left and right frequency responses with the roll-off function before designing the left and right compensation filters.

16. The method according to claim 15, where the high pass filter is a Bessel filter.

17. The method according to claim 15, where the cut-off frequency is equal to the average roll-off frequency multiplied by a factor, and where the factor is in the range 1.2-1.5.

18. The method according to claim 15, wherein the given threshold is in the range 10-30 dB.

19. The method according to claim 15, further comprising:

setting the left filter transfer function below the left roll-off frequency to be equal to the left filter transfer function at the left roll-off frequency, and

setting the right filter transfer function below the right roll-off frequency to be equal to the right filter transfer function at the right roll-off frequency.

20. The method according to claim 15, wherein the left and right filter transfer functions are set equal to unity gain above 500 Hz.

21. The method according to claim 20, further comprising cross fading the transfer function to unity gain over a suitable frequency range, such as 200 Hz to 500 Hz.

22. The method according to claim 15, wherein the simulated target function is obtained by simulating the power emitted by a point source in a corner defined by three orthogonal walls into a one eighth sphere limited by the three walls, and defining the simulated target function as the transfer function between the point source and the emitted power.

23. The method according to claim 22, wherein the simulated emitted power is a power average based on simulations in a plurality of points distributed on the one eighth sphere.

24. The method according to claim 23, wherein the simulated emitted power is the power average based on simulations in the plurality of points, including more than 12 points, distributed on the one eighth sphere.

25. A method for compensating for acoustic influence of a listening room on an acoustic output from an audio system including at least a left and a right loudspeaker, the method comprising:

determining a left frequency response LP_L between a signal applied to the left speaker and a resulting power average in a listening position,

determining a right frequency response LP_R between a signal applied to the right speaker and a resulting power average in the listening position,

designing a left compensation filter F_L,

designing a right compensation filter F_R,

during playback, applying the left compensation filter to a left input signal, and applying the right compensation filter to a right input signal,

providing a simulated target function H_T representing a simulated target response in the listening position,

designing the left compensation filter F_L to have a left filter transfer function based on the simulated target function H_T multiplied by an inverse of the left frequency response, and

designing the right compensation filter F_R to have a right filter transfer function based on the simulated target function H_T multiplied by an inverse of the right frequency response,

the method further comprising;

determining a filtered mono frequency response LP_M according to LP_L F_L+LP_R F_R,

determining a filtered side frequency response LP_s according to LP_L F_L−

LP_R F_R,

wherein LP_L is the left frequency response, LP_L is the right frequency response, F_L is the left filter and F_R is the right filter,

determining a mono target function based on the simulated target function H_T,

determining a side target function based on the simulated target function H_T,

designing a mono compensation filter F_M having a mono filter transfer function based on the mono target function multiplied by an inverse of the mono response,

designing a side compensation filter F_s having a side filter transfer function based on the side target function multiplied by an inverse of the side response, and

during playback, applying the mono compensation filter to a mono signal based on the left and right signal inputs, and applying the side compensation filter to a side signal based on the left and right input signals.

26. The method according to claim 25, wherein:

the mono signal is formed as the sum of a left input signal and a right input signal,

the side signal is formed as the difference between a left input signal and a right input signal,

the left filter input is formed as the sum of the filtered mono channel input and the filtered side channel input, and

the right filter input is formed as the difference between the filtered mono channel input and the side channel input.

27. The method according to claim 25, wherein the mono target function is determined as the simulated target function multiplied by a shelving filter with a center frequency in the order of 100 Hz and a gain in the order of one dB.

28. The method according to claim 25, wherein the side target function is determined as the mono target function reduced by a difference between a smoothed filtered mono response and a smoothed filtered side response.

29. The method according to claim 25, further comprising:

measuring a mono frequency response in the listening position,

applying the mono compensation filter to the measured mono response to form a filtered mono response,

forming a difference between the filtered mono frequency response and the mono target,

forming a peak removing component as portions of said difference smaller than zero, and

subtracting the peak removing component from the mono compensation filter and side compensation filter to form a peak cancelling mono compensation filter and a peak cancelling side compensation filter.

30. A method for compensating for acoustic influence of a listening room on an acoustic output from an audio system including at least a left and a right loudspeaker, the method comprising:

determining a left frequency response LP_L between a signal applied to the left speaker and a resulting power average in a listening position,

determining a right frequency response LP_R between a signal applied to the right speaker and a resulting power average in the listening position,

designing a left compensation filter F_L,

designing a right compensation filter F_R,

during playback, applying the left compensation filter to a left input signal, and applying the right compensation filter to a right input signal,

providing a simulated target function H_T representing a simulated target response in the listening position,

designing the left compensation filter F_L to have a left filter transfer function based on the simulated target function H_T multiplied by an inverse of the left frequency response,

designing the right compensation filter F_R to have a right filter transfer function based on the simulated target function H_T multiplied by an inverse of the right frequency response, and

smoothing at least one response using a method comprising;

determining a number of peaks per octave in the response,

for a portion of the response where the number of peaks per octave is below a first threshold, smoothing the response with a first smoothing width,

for a portion of the response where the number of peaks per octave is above a second threshold, smoothing the response with a second smoothing width,

wherein said second threshold is greater than said first threshold and said second smoothing width is wider than said first smoothing width, and

for a portion of the response where the number of peaks per octave is between the first and second thresholds, smoothing with an intermediate smoothing width.

31. An audio system including:

at least a left and a right loudspeaker arranged in a listening room;

at least one microphone arranged in a listening position;

a signal processing system for compensating for acoustic influence of the listening room on an acoustic output from the loudspeakers, said signal processing system being configured to;

apply a test signal to the left speaker, determine a power average based on a signal measured in the microphone, and determine a left frequency response LP_L between the test signal and the power average,

apply a test signal to the right speaker, determine a power average based on a signal measured in the microphone, and determine a right frequency response LP_L between the test signal and the power average,

design a left compensation filter F_L, and

design a right compensation filter F_R; and

a filtering system configured to;

during playback, apply the left compensation filter to a left channel input, and applying the right compensation filter to a right channel input,

wherein the signal processing system is provided with a simulated target function H_T representing a simulated target response in the listening position, and in that the signal processing system is configured to design the left compensation filter F_L to have a left filter transfer function based on the simulated target function H_T multiplied by an inverse of the left frequency response, and design the right compensation filter F_R to have a right filter transfer function based on the simulated target function H_T multiplied by an inverse of the right frequency response,

wherein the simulated target function is obtained by simulating the power emitted by a point source in a corner defined by three orthogonal walls into a one eighth sphere limited by the three walls, and defining the simulated target function as the transfer function between the point source and the emitted power,

wherein the simulated emitted power is a power average based on simulations in a plurality of points distributed on the one eighth sphere.

32. The system in claim 31, wherein the loudspeakers are directivity controlled loudspeakers.

33. The system of claim 31, wherein the simulated emitted power is the power average based on simulations in the plurality of points, including more than 12 points, distributed on the one eighth sphere.