Room acoustic response modeling and equalization with linear predictive coding and parametric filters
First Claim
1. A method for equalizing audio signals provided to speakers, the method comprising:
- measuring loudspeaker-room acoustics to obtain time domain room response data;
processing the time domain room response data with a Linear Predictive Coding (LPC) model to obtain smoothed time domain room response data;
performing an FFT on the smoothed time domain room response data to obtain smoothed frequency domain room response data;
computing center frequency Fc, a gain G, a bandwidth BW, and a term Q factor for a plurality of parametric Infinite-duration Impulse Response (IIR) filters from the smoothed frequency domain room response data includes a method for determining a number of peaks num_peaks comprising;
computing a frequency response HH2 from a linear predication coefficient q wherein elements in the frequency response HH2 correspond to an array of frequencies FF;
given an interested frequency range between LO_FREQ and HI_FREQ in Hz, and a low bin bin_lo and a high bin, bin_hi, updating the array of frequencies FF such that the elements of the updated the array of frequencies FF are the frequencies in the interested frequency range;
computing HH2_abs, the magnitude of HH2, based on the bin_lo and bin_hi in the interested frequency range;
determining peak locations peak_loc and valley_locations valley_loc while ensuring that a first peak occurs before a first valley;
saving a number of peak locations as the num_peaks; and
determining the center frequency Fc and gain G of each peak based on the peak location peak_loc and the magnitude HH2_abs at the peak location;
cascading the plurality of parametric IIR filters to form an equalizing filter;
equalizing a signal with the equalization filter to obtain an improved loudspeaker-room response; and
providing the equalized signal to a speaker to provide improved sound production.
6 Assignments
0 Petitions
Accused Products
Abstract
A method for determining coefficients of a family of cascaded second order Infinite Impulse Response (IIR) parametric filters used for equalizing a room response. The method includes determining parameters of each IIR parametric filter from poles or roots of a reasonably high-order Linear Predictive Coding (LPC) model. The LPC model is able to accurately model the low-frequency room response modes providing better equalization of loudspeaker and room acoustics, particularly at the low frequencies. Advantages of the method include fast and efficient computation of the LPC model using a Levinson-Durbin recursion to solve the normal equations that arise from the least squares formulation. Due to possible band interactions between the cascaded IIR parametric filters, the method further includes optimizing the Q value of each filter to better equalize the room response.
-
Citations
15 Claims
-
1. A method for equalizing audio signals provided to speakers, the method comprising:
-
measuring loudspeaker-room acoustics to obtain time domain room response data; processing the time domain room response data with a Linear Predictive Coding (LPC) model to obtain smoothed time domain room response data; performing an FFT on the smoothed time domain room response data to obtain smoothed frequency domain room response data; computing center frequency Fc, a gain G, a bandwidth BW, and a term Q factor for a plurality of parametric Infinite-duration Impulse Response (IIR) filters from the smoothed frequency domain room response data includes a method for determining a number of peaks num_peaks comprising; computing a frequency response HH2 from a linear predication coefficient q wherein elements in the frequency response HH2 correspond to an array of frequencies FF; given an interested frequency range between LO_FREQ and HI_FREQ in Hz, and a low bin bin_lo and a high bin, bin_hi, updating the array of frequencies FF such that the elements of the updated the array of frequencies FF are the frequencies in the interested frequency range; computing HH2_abs, the magnitude of HH2, based on the bin_lo and bin_hi in the interested frequency range; determining peak locations peak_loc and valley_locations valley_loc while ensuring that a first peak occurs before a first valley; saving a number of peak locations as the num_peaks; and determining the center frequency Fc and gain G of each peak based on the peak location peak_loc and the magnitude HH2_abs at the peak location; cascading the plurality of parametric IIR filters to form an equalizing filter; equalizing a signal with the equalization filter to obtain an improved loudspeaker-room response; and providing the equalized signal to a speaker to provide improved sound production. - View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9)
-
-
10. A method for computing coefficients of a family of cascaded parametric IIR filters and using the filters to filter signals provided to speakers, the method comprising:
-
collecting unprocessed time domain room response data; performing an FFT on the time domain room response data to obtain a frequency domain room response; normalizing the frequency domain room response in a frequency range of interest to obtain a normalized frequency domain room response; performing an inverse FFT on the normalized frequency domain room response to obtain a normalized time domain room response data; representing the normalized time domain room response data using an LPC model to obtain smoothed time domain room response data; performing an FFT on the smoothed time domain room response data to obtain smoothed frequency domain room response data; inverting the smoothed frequency domain room response data to obtain an equalization frequency response; computing a magnitude of the equalization frequency response; detecting peaks and valleys of the magnitude of the equalization frequency response; computing gains G, center frequencies Fc, bandwidths BW and bandwidth term Q factors of each of the peaks for a plurality of the parametric Infinite-duration Impulse Response (IIR) filters using a step for determining a number of peaks num_peaks, the step comprising; computing a frequency response HH2 from a linear predication coefficient q wherein elements in the frequency response HH2 correspond to an array of frequencies FF; given an interested frequency range between LO_FREQ and HI_FREQ in Hz, and a low bin bin_lo and a high bin, bin_hi, updating the array of frequencies FF such that the elements of the updated array of frequencies FF are frequencies in the interested frequency range; computing HH2_abs, a magnitude of the frequency response HH2, based on the bin_lo and the bin_hi in the interested frequency range; determining peak locations peak_loc and valley_locations valley_loc while ensuring that a first peak occurs before a first valley; saving a number of peak locations as the num_peaks; and determining the center frequency Fc and gain G of each peak based on the peak location peak_loc and the magnitude of the HH2_abs at the peak location; optimizing the gains and the Q factors; computing parametric filter coefficients from the optimized gains and the optimized Q factors; using the parametric filter coefficients to equalize signals; and providing the equalized signals to speakers to provide improved sound production. - View Dependent Claims (11, 12, 13)
-
-
14. A method for computing the coefficients of a family of cascaded parametric IIR filters and using the filters to filter signals provided to speakers, the method comprising:
-
collecting unprocessed time domain room response data; performing an FFT on the time domain room response data to obtain a frequency domain room response; normalizing the frequency domain room response in a frequency range of interest to obtain a normalized frequency domain room response; performing an inverse FFT on the normalized frequency domain room response to obtain a normalized time domain room response data; representing the normalized time domain room response data using an LPC model to obtain smoothed time domain room response data; performing an FFT on the smoothed time domain room response data to obtain smoothed frequency domain room response data; computing the magnitude of the smoothed frequency domain room response; detecting peaks and valleys of the magnitude of the smoothed frequency domain room response; computing gains G, center frequencies Fc, bandwidths BW and bandwidth term Q factors of each of the peaks for a plurality of the parametric Infinite-duration Impulse Response (IIR) filters using a step for determining a number of peaks num_peaks, the step comprising; computing a frequency response HH2 from a linear predication coefficient q wherein elements in the frequency response HH2 correspond to an array of frequencies FF; given an interested frequency range between LO_FREQ and HI_FREQ in Hz, and a low bin bin_lo and a high bin, bin_hi, updating the array of frequencies FF such that the elements of the updated array of frequencies FF are frequencies in the interested frequency range; computing HH2_abs, a magnitude of the frequency response HH2, based on the bin_lo and the bin_hi in the interested frequency range; determining peak locations peak_loc and valley_locations valley_loc while ensuring that a first peak occurs before a first valley; saving a number of peak locations as the num_peaks; and determining the center frequency Fc and gain G of each peak based on the peak location peak_loc and the magnitude of the HH2_abs at the peak location; optimizing the gains and the Q factors; computing parametric filter coefficients from the center frequencies and the optimized gains and the optimized Q factors; determining poles and zeros of each of the parametric IIR filters in a cascade based on the parametric filter coefficients; computing minimum-phase zeros from the zeros of each of the parametric filters in the cascade; reflecting each minimum-phase zero as a reflected pole and reflecting each pole as a reflected zero for each parametric filter in the cascade; expanding each reflected zero and its complex conjugate into a real second order numerator polynomial and expanding each reflected pole and its complex conjugate into a real second order denominator polynomial for each parametric filter in the cascade; using the parametric filters to equalize signals; and providing the equalized signals to speakers to provide improved sound production.
-
-
15. A method for computing the coefficients of a family of cascaded parametric IIR filters and using the filters to filter signals provided to speakers, the method comprising:
-
collecting unprocessed time domain room response data; performing an FFT on the time domain room response data to obtain a frequency domain room response; normalizing the frequency domain room response in a frequency range of interest to obtain a normalized frequency domain room response; performing an inverse FFT on the normalized frequency domain room response to obtain a normalized time domain room response data; representing the normalized time domain room response data using an LPC model to obtain smoothed time domain room response data; performing an FFT on the smoothed time domain room response data to obtain smoothed frequency domain room response data; computing the magnitude of the smoothed frequency domain room response to obtain a magnitude response; inverting the magnitude response; detecting peaks and valleys of the inverted magnitude response; computing gains G, center frequencies Fc, bandwidths BW and bandwidth term Q factors of each of the peaks for a plurality of the parametric Infinite-duration Impulse Response (IIR) filters using a step for determining a number of peaks num_peaks, the step comprising; computing a frequency response HH2 from a linear predication coefficient q wherein elements in the frequency response HH2 correspond to an array of frequencies FF; given an interested frequency range between LO_FREQ and HI_FREQ in Hz, and a low bin bin_lo and a high bin, bin_hi, updating the array of frequencies FF such that the elements of the updated array of frequencies FF are frequencies in the interested frequency range; computing HH2_abs, a magnitude of the frequency response HH2, based on the bin_lo and the bin_hi in the interested frequency range; determining peak locations peak_loc and valley_locations valley_loc while ensuring that a first peak occurs before a first valley; saving a number of peak locations as the num_peaks; and determining the center frequency Fc and gain G of each peak based on the peak location peak_loc and the magnitude of the HH2_abs at the peak location; optimizing the gains and the Q factors; computing parametric filter coefficients from the optimized center frequencies, the optimized gains, and the optimized Q factors; using the parametric filters to equalize signals; and providing the equalized signals to speakers to provide improved sound production.
-
Specification