Method for reducing interference in acoustic signals using an adaptive filtering method involving spectral subtraction

• US 6,643,619 B1
• Filed: 06/20/2000
• Issued: 11/04/2003
• Est. Priority Date: 10/30/1997
• Status: Expired due to Term

• Current Assignee: Nuance Communications, Inc. (Microsoft Corporation)
• Original Assignee: DaimlerChrysler AG (Mercedes-Benz Group AG)
• Inventors: Haulick, Tim, Linhard, Klaus
• Primary Examiner(s): Chawan, Vijay
• Application Number:
  US09/530,527
• Time in Patent Office: 1,232 Days
• Field of Search: 704/233, 704/226, 704/227, 704/228, 704/210, 704/214, 704/219, 704/230, 704/234, 704/200, 704/207, 704/240, 704/253, 379/406.08, 381/94.1, 381/94.2, 381/94.3, 381/94.7
• US Class Current: 704/233
• CPC Class Codes: G10L 21/02 Speech enhancement, e.g. noise reduction or echo cancellation reducing echo effects in line transmission systems H04B3/20; echo suppression in hands-free telephones H04M9/08

Abstract

A method for reducing interference in acoustic signals by using of an adaptive filter method involving spectral subtraction. The inventive method enables a significant reduction of interference in acoustic signals, especially voice signals, without causing any substantial falsification of said signals such as echo or musical tones, and significantly reduces computational requirements in comparison with other methods known per se that are similarly designed to improve signal quality.

Claims

1. A method for reducing interference in disturbed acoustic signals using an adaptive filtering process including spectral subtraction, the method comprising:

filtering the signals in a plurality of respective time segments and a plurality of respective discrete frequencies i segmentwise using an adaptive filtering function;

determining a respective noise-input ratio for each of the plurality of respective time segments and respective discrete frequencies so that each respective noise-input ratio has a small respective value for signals having a relatively low disturbing noise component and a high respective value for signals having a relatively high disturbing noise component;

adapting the adaptive filtering function so that respective information on a respective a priori signal-to-noise ratio is used for a calculation of each of a plurality of characteristic values of the adaptive filtering function; and

using at least one of the plurality of characteristic values from a respective at least one preceding time segment as the respective information on each respective a priori signal-to-noise ratio;

the adaptive filtering function having a characteristic curve including two parts and having a break edge positioned such that the filtering for signals having a high respective noise-input ratio results in a signal-independent relatively strong damping and the filtering for signals having a low noise-input ratio results in a signal-dependent relatively low damping.

2. The method as recited in claim 1 wherein in the using step the characteristic value from only the immediately preceding time segment is used as the information on the a priori signal-to-noise ratio.

3. The method as recited in claim 1 wherein each of the plurality of characteristic values is calculated using a respective corrected noise-input ratio, the respective corrected noise-input ratio being calculated using the respective noise-input ratio so as to use the information on the respective a priori signal-to-noise ratio.

4. The method as recited in claim 3 wherein each respective noise-input ratio is calculated using ${\mathrm{NIR}}^{\prime }$



(k,i);

=NIR

(k,i)/∑

j=1N

 



wj

H

(k-j,i)

where NIR′

(k, i) is the corrected noise-input ratio, NIR(k, i) is the noise-input ratio, k is the respective time segment, i is the respective frequency, H(k−

j, i) is a respective one of the plurality of characteristic values, the weighting factors w_j are real numbers smaller than 1, and N is a natural number greater than or equal to 1.

5. The method as recited in claim 3, wherein the filtering function is calculated using at least one of:

6. The method as recited in claim 5 wherein a is an element of an interval from 1 to 4 and b is an element of an interval from 0.1 to 0.3.

7. The method as recited in claim 1 further comprising adapting the position of the break edge of the characteristic curve of the adaptive filtering function to the frequency of the signal being filtered.

8. The method as recited in claim 7 wherein each of the plurality of characteristic values is calculated using the respective noise-input ratio and wherein the adapting of the position of the break edge is performed by replacing each respective noise-input ratio with a respective corrected noise-input ratio for the calculating of the respective characteristic value.

9. The method as recited in claim 8 wherein the corrected noise-input ratio is calculated using:

${\mathrm{NIR}}^{\mathrm{\prime <br><br>}}<br><br>\left(k,i\right);<br><br>=\mathrm{NIR}<br><br>\left(k,i\right)/\left[c<br><br>\left(i\right)+\left(1-c<br><br>\left(i\right)\right)<br><br>\underset{j=1}{\overset{N}{\sum <br><br>}}<br><br>\text{\hspace{1em}<br><br>}<br><br>w_j<br><br>H<br><br>\left(k-j,i\right)\right]$

where NIR′

(k, i) is the corrected noise-input ratio, NIR(K, i) is the noise-input ratio, k is the respective time segment, i is the respective frequency, H(k−

j, i) is a respective one of the plurality of characteristic values, the weighting factors w_j are real numbers smaller than 1, and N is a natural number greater than or equal to 1.

10. The method as recited in claim 9 wherein the corrected noise-input ratio is calculated using:

11. The method as recited in claim 8 further comprising correcting the respective characteristic values from the at least one preceding time segment prior to calculating each respective corrected noise-input ratio.

12. The method as recited in claim 11 wherein the correcting of each of the respective at characteristic values is performed using

13. The method as recited in claim 1 further comprising adapting the position of the break edge as a function of a presence of a speech signal and a presence of a speech pause.

14. The method as recited in claim 13 wherein each of the plurality of characteristic values is calculated using the respective noise-input ratio and wherein the adapting of the position of the break edge is performed by replacing each respective noise-input ratio with a respective corrected noise-input ratio for the calculating of the respective characteristic value.

15. The method as recited in claim 14 wherein the corrected noise-input ratio is calculated using:

${\mathrm{NIR}}^{\mathrm{\prime <br><br>}}<br><br>\left(k,i\right);<br><br>=\mathrm{NIR}<br><br>\left(k,i\right)/\left[c<br><br>\left(i\right)+\left(1-c<br><br>\left(i\right)\right)<br><br>\underset{j=1}{\overset{N}{\sum <br><br>}}<br><br>\text{\hspace{1em}<br><br>}<br><br>w_j<br><br>H<br><br>\left(k-j,i\right)\right]$

where NIR′

(k, i) is the corrected noise-input ratio, NIR(k, i) is the noise-input ratio, k is the respective time segment, i is the respective frequency, H(k−

j, i) is a respective one of the plurality of characteristic values, the weighting factors w_j are real numbers smaller than 1, and N is a natural number greater than or equal to 1.

16. The method as recited in claim 15 wherein the corrected noise-input ratio is calculated using:

17. The method as recited in claim 14 further comprising correcting the respective characteristic values from the at least one preceding time segment prior to calculating each respective corrected noise-input ratio.

18. The method as recited in claim 17 wherein the correcting of each of the respective at characteristic values is performed using

19. The method as recited in claim 1 wherein the acoustic signals are speech signals.