Abstract:
A hearing aid comprises at least one microphone (M) for converting input sound into an input signal, a subtraction node for subtracting a feedback cancellation signal from the input signal thereby generating a processor input signal, a hearing aid processor (G) for producing a processor output signal by applying an amplification gain to the processor input signal, a receiver (R) for converting the processor output signal into output sound, an adaptive feedback cancellation filter for adaptively deriving the feedback cancellation signal from the processor output signal by applying filter coefficients, calculation means for calculating the autocorrelation of a reference signal, and an adaptation means for adjusting the filter coefficients with an adaptation rate, wherein the adaptation rate is controlled in dependency of the autocorrelation of the reference signal.

Description:
RELATED APPLICATIONS 
       [0001]    The present application is a continuation-in-part of application no. PCT/EP2007053175 filed on Apr. 2, 2007 and published as WO-A1-2007113282, the contents of which are incorporated herein by reference. 
         [0002]    The present application is based on and claims priority from PA200600467 filed on Apr. 1, 2006 in Denmark, the contents of which are incorporated herein by reference. 
     
    
     BACKGROUND OF THE INVENTION 
       [0003]    1. Field of the Invention 
         [0004]    The present invention relates to hearing aids. The invention, more particularly relates to hearing aids that rely on adaptive feedback cancellation in order to reduce the problems caused by acoustic and mechanical feedback. The invention further relates to methods for control of the adaptation rate in feedback cancelling systems and hearing aids and to hearing aids and systems that incorporate such methods. 
         [0005]    2. Description of the Related Art 
         [0006]    Acoustic and mechanical feedback from a receiver to one or more microphones will limit the maximum amplification that can be applied in a hearing aid. Due to the feedback, the amplification in the hearing aid can cause resonances, which shape the spectrum of the output of the hearing aid in undesired ways and even worse, it can cause the hearing aid to become unstable, resulting in whistling or howling. The hearing aid usually employs compression to compensate hearing loss; that is, the amplification gain is reduced with increasing sound pressures. Moreover, an automatic gain control is commonly used on the output to limit the output level, thereby avoiding clipping of the signal. In case of instability, these compression effects will eventually make the system marginally stable, thus producing a howl or whistle of nearly constant sound level. 
         [0007]    Feedback cancellation is often used in hearing aids to compensate the acoustic and mechanical feedback. The acoustic feedback path can change dramatically over time as a consequence of, for example, amount of earwax, the user wearing a hat or holding a telephone to the ear or the user chewing or yawning. For this reason it is customary to apply an adaptation mechanism on the feedback cancellation to account for the time-variations. 
         [0008]    An adaptive feedback cancellation filter can be implemented in a hearing aid in several different ways. For example, it can be IIR, FIR or a combination of the two. It can be composed of a combination of a fixed filter and an adaptive filter. The adaptation mechanism can be implemented in several different ways, for example algorithms based on Least Mean Squares (LMS) or Recursive Least Squares (RLS). 
         [0009]    In A. Spriet, I. Proudler, M. Moonen, J. Wouters: “ Adaptive Feedback Cancellation in Hearing Aids With Linear Prediction of the Desired Signal ”, IEEE Trans. On Signal Processing, Vol. 53, No. 10, October 2005 it is described that the accuracy of the estimated feedback cancelling filter is degraded when the incoming signal is spectrally coloured. This is also mentioned in patent application WO 01/06812, “Feedback Cancellation with Low Frequency Input”. This patent describes a scheme in which an adaptive resonator filter is used for detecting if a dominating tone is present in the signal, in which case the adaptation rate is significantly increased. This allows for a rapid and efficient cancellation of feedback howl. The drawback is that if the tone is not due to feedback but is present in the environment, the adaptive feedback cancelling may react strongly on this signal, with the risk of noticeable audible artefacts. 
         [0010]    In Moonen et al. and WO 01/06812 it is further mentioned that it will lead to bias errors in the model of the acoustic feedback if the microphone signal is spectrally coloured. 
         [0011]    The patent application WO 99/26453, “Feedback Cancellation Apparatus and Methods” describes a feedback cancellation system in which separate cancellation filters are used for compensating the acoustic feedback to each microphone in a two-microphone hearing aid. In contrast to prior art in the field, this has the advantage that an adaptive directional system for spatial noise filtering is not treated as an integral part of the acoustic feedback path. 
         [0012]    The patent application WO 02/25996 describes a scheme for an adaptive feedback cancellation filter as well as a scheme for stabilization of the hearing aid by using a procedure for estimation of the current stability limit. 
         [0013]    LMS and other adaptation algorithms are derived and discussed in the book: S. Haykin:  Adaptive Filter Theory,  3rd Edition, Prentice-Hall, NJ, USA, 1996. 
         [0014]    Further details on convergence and behaviour of the LMS and Normalized LMS algorithms are provided in D. T. M Slock:  On the Convergence Behavior of the LMS and the Normalized LMS Algorithms , IEEE Trans. Signal Processing, Vol. 41, No. 9, September 1993, pp. 2811-2824. 
         [0015]    Even though many recommendations has been given in the prior art as to how the adaptation rate in such systems should be decided on, there still exists a need for improvements in this area. In particular, there exists a need for hearing aids in which methods for automatic adjustment of this rate, in dependency of the acoustic environment, have been implemented. 
       SUMMARY OF THE INVENTION 
       [0016]    On the background described herein, it is an object of some embodiments of the present invention to provide a method and a hearing aid of the kind defined, in which the deficiencies of the prior art methods and hearing aids are remedied by automatically adjusting the adaptation rate of feedback cancellation in dependency of the acoustic environment. 
         [0017]    Particularly, it is an object of some embodiments of the present invention to provide a method and a hearing aid allowing to implement specific procedures for selecting an appropriate adaptation step size in feedback cancellation. 
         [0018]    It is a further object of some embodiments of the present invention to provide a method and a hearing aid allowing to reduce the error in the estimate of the feedback path of the hearing aid. 
         [0019]    It is yet a further object of some embodiments of the present invention to provide a method and a hearing aid allowing to cope with the sensitivity of adaptive feedback cancelling systems to tonal input signals. 
         [0020]    It is still another object of some embodiments of the present invention provide a method and a hearing aid allowing to cope with the sensitivity of adaptive feedback cancelling systems to tonal input signals by preventing the onset of feedback initiated oscillation. 
         [0021]    It is yet another object of some embodiments of the present invention to provide a method and a hearing aid allowing to cope with the impact of the gain size onto the error in the estimate of the feedback path of the hearing aid. 
         [0022]    It is further an object of some embodiments of the present invention to provide a method and a hearing aid allowing to cope with the impact of non-continuous sound in the environment of the hearing aid onto the error in the estimate of the feedback path of the hearing aid. 
         [0023]    It is further an object of some embodiments of the present invention provide a method and a hearing aid allowing to cope with the impact of the presence of an adaptive microphone array, and hence the total gain of the hearing aid, onto the error in the estimate of the feedback path of the hearing aid. 
         [0024]    It is further an object of some embodiments of the present invention to provide a method and a hearing aid allowing to control the step size in the adaptive algorithm of a feedback cancelling system taking multiple aspects of the acoustic environment into account. 
         [0025]    According to the invention several suggestions as to how the adaptation rate should be controlled are given. In particular, it is suggested how the adaptation rate may be automatically adjusted in dependency of the acoustic environment. 
         [0026]    The invention, in a first aspect, provides a hearing aid comprising at least one microphone for converting input sound into an input signal; a subtraction node for subtracting a feedback cancellation signal from the input signal thereby generating a processor input signal; a hearing aid processor for producing a processor output signal by applying an amplification gain to the processor input signal; a receiver for converting the processor output signal into output sound; an adaptive feedback cancellation filter for adaptively deriving the feedback cancellation signal from the processor output signal by applying filter coefficients; calculation means for calculating an autocorrelation value of a reference signal; and an adaptation means for adjusting the filter coefficients with an adaptation rate, wherein the adaptation rate is set in dependency of the autocorrelation value calculated for the reference signal. This arrangement allows an improved adjustment of the adaptation rate taking the sensitivity of adaptive feedback systems like adaptive feedback cancellation filters to tonal input signals into account. 
         [0027]    A method for control of the adaptation rate in a hearing aid comprising: converting input sound into an input signal; subtracting a feedback cancellation signal from the input signal thereby generating a processor input signal; producing a processor output signal by applying an amplification gain to the processor input signal; converting the processor output signal into output sound; adaptively deriving the feedback cancellation signal from the processor output signal by applying filter coefficients; calculating an autocorrelation value of a reference signal; and adjusting the filter coefficients with an adaptation rate, wherein the adaptation rate is set in dependency of the autocorrelation value of the reference signal. This arrangement allows an improved adjustment of the adaptation rate taking the importance of gain size to the error in the filter coefficients and, hence, the error in the estimate of the feedback path of the hearing aid into account. 
         [0028]    According to an embodiment the hearing aid comprises detection means for detecting if the input signal represents a sudden increase in sound pressure of the input sound, and wherein the adaptation means is adapted to temporarily suspend the adjustment of the filter coefficients. This arrangement allows an improved adjustment of the adaptation rate taking the importance of non-continuous sound in the environment of the feedback path of the hearing aid into account. 
         [0029]    According to another embodiment, there is provided a hearing aid comprising at least two microphones converting the input sound in at least a first and a second spatial input signal providing a directional characteristic, at least two subtraction nodes for subtracting a first feedback cancellation signal from the first input signal and a second feedback cancellation signal from the second input signal thereby generating a resulting directional processor input signal, at least a first and a second adaptive feedback cancellation filter for adaptively deriving the first and second feedback cancellation signals, and wherein said adaptation means is adapted to further control the adaptation rate in dependency of the directional characteristic. This arrangement allows an improved adjustment of the adaptation rate taking the importance of the contribution of a directional microphone system providing momentary gain or attenuation to the overall system gain into account. 
         [0030]    The invention, in a third aspect, provides a computer program product comprising program code for performing, when run on a computer, a method for control of the adaptation rate in a hearing aid comprising: converting input sound into an input signal; subtracting a feedback cancellation signal from the input signal thereby generating a processor input signal; producing a processor output signal by applying an amplification gain to the processor input signal; converting the processor output signal into output sound; adaptively deriving the feedback cancellation signal from the processor output signal by applying filter coefficients; calculating an autocorrelation value of a reference signal; and adjusting the filter coefficients with an adaptation rate, wherein the adaptation rate is set in dependency of the autocorrelation value of the reference signal. 
         [0031]    The present invention lays out a number of schemes for adaptively setting the adaptation rate in an algorithm used for adjusting the coefficients in a feedback cancelling filter in a hearing aid. The adaptation rate is varied in accordance with the characteristics of the microphone signal(s) and the various internal parameters and signals inside the hearing aid. According to the present invention, specific ways are provided for adjusting the adaptation rate based on observations of the current microphone signal(s), the present state and/or the behaviour of the hearing aid. 
         [0032]    Further aspects, embodiments, and specific variations of the invention are defined by the further dependent claims. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0033]    The invention will now be described in greater detail based on non-limiting examples of preferred embodiments and with reference to the appended drawings. On the drawings: 
           [0034]      FIG. 1  shows a hearing aid with an adaptive feedback cancellation filter, according to the prior art; 
           [0035]      FIG. 2  shows a hearing aid with a feedback adaptation mechanism, according to the prior art; 
           [0036]      FIG. 3  shows a hearing aid with two microphones and two adaptive feedback cancellation filters, according to the prior art; 
           [0037]      FIG. 4  shows a schematic block diagram of a hearing aid according to an embodiment of the present invention; 
           [0038]      FIG. 5  shows a schematic block diagram of the hearing aid of  FIG. 4 , with schematic illustrations of the effect of signals with high autocorrelation; 
           [0039]      FIG. 6  shows a schematic block diagram of a hearing aid according to an embodiment of the present invention with means for detecting a sudden sound; 
           [0040]      FIG. 7  shows a schematic block diagram of a prior art hearing aid with directional characteristics; 
           [0041]      FIG. 8  shows a hearing aid with an adaptive feedback cancelling filter and with directional characteristic, according to an embodiment of the invention; 
           [0042]      FIG. 9  shows a hearing aid with an adaptive feedback cancelling filter and with a step-size control block, according to an embodiment of the invention; 
           [0043]      FIG. 10  shows a hearing aid with two microphones and with two adaptive feedback cancelling filters, according to an embodiment of the invention; 
           [0044]      FIG. 11  shows a hearing aid with two microphones and with one adaptive feedback cancelling filter, according to an embodiment of the invention; and 
           [0045]      FIG. 12  shows a hearing aid with two microphones and with a step-size control, according to an embodiment of the invention. 
       
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
       [0046]      FIGS. 1-3  show schematic block diagrams of prior art hearing aids implementing some basic feedback cancellation schemes. 
         [0047]    In  FIG. 1 , the microphone signal  1  from the microphone M is compensated by subtraction of the feedback cancelling signal  4 . The resulting signal  2  is used as input to the hearing aid processor  100 , and it is used as adaptation error in the adaptive feedback cancelling filter  101 . The output of the hearing aid processor is transmitted to the receiver R. The hearing aid processor  100  may comprise time-varying and frequency dependent filters to account for the hearing loss, suppression of noise, automatic gain control for handling large signals, and time-delays. The block  101  represents an adaptive feedback cancellation filter and embraces a simultaneous filtering and adaptation of filter coefficients. 
         [0048]    The diagram in  FIG. 2  shows a system like the one depicted in  FIG. 1  except that the adaptation mechanism implemented in block  103  is separated from the filtering function implemented in block  102 . The connection  5  symbolizes the filter coefficients. The advantage of this scheme over the one shown in  FIG. 1  is that a frequency shaping of the signals  2  and  3  can be made without disturbing the filtering performance. 
         [0049]    The diagram in  FIG. 3  shows how multiple feedback cancellation filters  202   a ,  202   b  can be used in the case of hearing aids with multiple microphones M 1 , M 2 . In this case two sets of filter coefficients  38   a ,  38   b  are passed on from the adaptation block  203 . In the example shown here, the two cancellation signals  35 ,  36  compensate the signals  30 ,  31 , which are created employing two spatial filters of the sound  206 ,  207 , each filter with its own fixed directional pattern (e.g., such than one is omnidirectional and one is bipolar). The compensated signals  32 ,  33  are subsequently weighted in order to achieve a resulting directional signal. This weighting can be time-varying as this will allow adaptation of the resulting directional pattern to the current sound environment. A band-split into several frequency bands is possible in e.g.,  205  as this will make it possible to vary the directional pattern over frequency, thus allowing improved noise reduction. The signal  34  will in this case be a multi-band signal. 
         [0050]    Further terms and prerequisites useful for understanding the present invention will be explained when describing particular embodiments of the present invention in the following. 
       Autocorrelation Dependency 
       [0051]    The extent to which a signal, x k , is spectrally coloured is often measured by the autocorrelation of the signal: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       R 
                       x 
                     
                      
                     
                       ( 
                       τ 
                       ) 
                     
                   
                   = 
                   
                     
                       ∑ 
                       
                         k 
                         = 
                         τ 
                       
                       N 
                     
                      
                     
                         
                     
                      
                     
                       
                         x 
                         k 
                       
                        
                       
                         x 
                         
                           k 
                           - 
                           τ 
                         
                       
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     1 
                   
                   ] 
                 
               
             
           
         
       
     
         [0052]    where τ is the time lag. For white noise, R x (τ)≈0 for all τ≠0. For periodic signals or other signals with a certain amount of predictability, the autocorrelation will be significantly larger than 0 for one or more time lags. 
         [0053]    To better allow comparison, the autocorrelation is often normalized with the window size or with the autocorrelation at lag 0: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         R 
                         x 
                         N 
                       
                        
                       
                         ( 
                         τ 
                         ) 
                       
                     
                     = 
                     
                       
                         1 
                         N 
                       
                        
                       
                         
                           ∑ 
                           
                             k 
                             = 
                             τ 
                           
                           N 
                         
                          
                         
                             
                         
                          
                         
                           
                             x 
                             k 
                           
                            
                           
                             x 
                             
                               k 
                               - 
                               τ 
                             
                           
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   or 
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     2 
                   
                   ] 
                 
               
             
             
               
                 
                   
                     
                       r 
                       x 
                     
                      
                     
                       ( 
                       τ 
                       ) 
                     
                   
                   = 
                   
                     
                       
                         ∑ 
                         
                           k 
                           = 
                           τ 
                         
                         N 
                       
                        
                       
                           
                       
                        
                       
                         
                           x 
                           k 
                         
                          
                         
                           x 
                           
                             k 
                             - 
                             τ 
                           
                         
                       
                     
                     
                       
                         ∑ 
                         
                           k 
                           = 
                           0 
                         
                         N 
                       
                        
                       
                           
                       
                        
                       
                         
                           x 
                           k 
                         
                          
                         
                           x 
                           k 
                         
                       
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     3 
                   
                   ] 
                 
               
             
           
         
       
     
         [0054]    The autocorrelation coefficients given by the last equation have the property that the values are limited to [−1; 1]. 
         [0055]    In a practical non-stationary setting, the autocorrelation must be calculated over a sliding window or according to some kind of recursive update. An embodiment of this is to use a sliding average in place of the sum in [Eq. 2]: 
         [0000]        R   x (τ, k )= R   x (τ, k− 1)+α·( x   k   x   k-τ   −R   x (τ, k− 1))  [Eq. 4] 
         [0000]    where αε]0;1[ controls the weighting between historic and current signal values. 
         [0056]    In a hearing aid context, this update can be quite costly to calculate because many multiplications are required. Particularly if many different lags, τ, are considered or if the calculation is carried out in several frequency bands. Instead, it might be relevant to consider updates that do not approximate the autocorrelation but something, which in a similar sense measures how systematic or predictable a signal is. Two embodiments, both quite simple to compute, as they do not depend on multiplications, are 
         [0000]        R   x (τ, k )= R   x ( r,k− 1)+α·( z (τ, k )− R   x (τ, k− 1)) 
         [0000]        z (τ, k )= x   k sign( x   k-τ ) 
         [0000]        z (τ, k )=sign( x   k )sign( x   k-τ )  [Eq. 5] 
         [0057]    The co-pending patent application DK 2006 00479 “Method for controlling signal processing in a hearing aid and a hearing aid implementing this method”, filed on Apr. 3, 2006, in Denmark, published as WO2007113283, and which is hereby incorporated by reference, describes these along with other signal characterization quantities related to the autocorrelation that can often be used instead of the true autocorrelation. 
         [0058]    The autocorrelation can be calculated for a wide-band signal or it can be calculated for a number of band-limited signals. In order to detect if a pure tone is present in the signal, it can be relevant to calculate the autocorrelation coefficients in a number of bands and subsequently look for the maximum of absolute values of the autocorrelation for several time lags and for all frequency bands. 
         [0059]    For several reasons, adaptive anti-feedback systems are often based on the adaptive scheme outlined by a variation of the Least Mean Square (LMS) algorithm. As a simple example, we can consider an adaptive FIR filter: 
         [0000]        {circumflex over (f)}   k   =w (0) x   k   +w (1) x   k-1   + . . . +w ( M ) x   k-M   [Eq. 6] 
         [0060]    Provided that y k  is the observed signal, which contains information about the underlying system we wish to model, the filter coefficients are adjusted according to e.g., 
       LMS: 
       [0061]        w   k ( i )= w   k-1 ( i )+μ x   k-i ( y   k   −{circumflex over (f)}   k )  [Eq. 7] 
       Normalized LMS, NLMS: 
       [0062]    
       
         
           
             
               
                 
                   
                     
                       w 
                       k 
                     
                      
                     
                       ( 
                       i 
                       ) 
                     
                   
                   = 
                   
                     
                       
                         w 
                         
                           k 
                           - 
                           1 
                         
                       
                        
                       
                         ( 
                         i 
                         ) 
                       
                     
                     + 
                     
                       
                         
                           μ 
                            
                           
                               
                           
                            
                           
                             x 
                             
                               k 
                               - 
                               i 
                             
                           
                         
                         
                           
                             ∑ 
                             
                               j 
                               = 
                               k 
                             
                             
                               j 
                               = 
                               
                                 k 
                                 + 
                                 M 
                               
                             
                           
                            
                           
                             x 
                             k 
                             2 
                           
                         
                       
                        
                       
                         ( 
                         
                           
                             y 
                             k 
                           
                           - 
                           
                             
                               f 
                               ^ 
                             
                             k 
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     8 
                   
                   ] 
                 
               
             
           
         
       
     
         [0000]    LMS with Variance Normalization: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         w 
                         k 
                       
                        
                       
                         ( 
                         i 
                         ) 
                       
                     
                     = 
                     
                       
                         
                           w 
                           
                             k 
                             - 
                             1 
                           
                         
                          
                         
                           ( 
                           i 
                           ) 
                         
                       
                       + 
                       
                         
                           μ 
                           
                             
                               σ 
                               ^ 
                             
                             k 
                             2 
                           
                         
                          
                         
                           
                             x 
                             
                               k 
                               - 
                               i 
                             
                           
                            
                           
                             ( 
                             
                               
                                 y 
                                 k 
                               
                               - 
                               
                                 
                                   f 
                                   ^ 
                                 
                                 k 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       
                         
                           σ 
                           ^ 
                         
                         k 
                         2 
                       
                       = 
                       
                         
                           ρ 
                            
                           
                               
                           
                            
                           
                             
                               σ 
                               ^ 
                             
                             
                               k 
                               - 
                               1 
                             
                             2 
                           
                         
                         + 
                         
                           
                             ( 
                             
                               1 
                               - 
                               ρ 
                             
                             ) 
                           
                            
                           
                             x 
                             k 
                             2 
                           
                         
                       
                     
                     , 
                     
                       0 
                       ≤ 
                       ρ 
                       &lt; 
                       1 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     9 
                   
                   ] 
                 
               
             
           
         
       
     
       Sign-Sign LMS: 
       [0063]        w   k ( i )= w   k-1 ( i )+μsign( x   k-i )sign( y   k   −{circumflex over (f)}   k )  [Eq. 10] 
         [0064]    A person skilled in the art however will appreciate that calling the latter an LMS-type algorithm is in a literal sense slightly misleading. 
         [0065]    The person skilled in the art will further appreciate that many variations can be made on both filter and algorithm. The adaptive FIR filter can be substituted by a warped delay line, a fixed pre-filter or post-filter can be used, or the filter can be an adaptive IIR-filter. There is a plethora of possible adaptation algorithms in addition to the ones shown. 
         [0066]    To accommodate the non-stationary nature of sound environments that a hearing aid user can be exposed to and the highly time-varying signal processing occurring in modern hearing aids, it is beneficial to let the step size, μ, be time-varying. The present invention deals with specific procedures for selecting an appropriate step size or adaptation speed or rate as will be described in detail below. 
         [0067]    The invention is particularly useful in relation to the NLMS algorithm as described in Eq. 8, or algorithms exhibiting a similar behaviour, such as the LMS with variance normalization, as described in Eq. 9. The principles are, however, relevant regardless of the implemented adaptation algorithm and may be implemented in various embodiments according to the present invention. 
         [0068]    With reference to  FIGS. 4 and 5 , an embodiment of the present invention will be discussed in connection with the presence of a spectrally coloured microphone signal. The hearing aid basically comprises microphone M, processor G, receiver R, and feedback cancellation filter {circumflex over (F)}. Considering  FIG. 5  but disregarding initially the adaptive feedback cancellation branch expressed by the filter {circumflex over (F)}, it is assumed that the incoming sound, v, is a pure tone (sinusoid). The microphone output y will then be a sinusoid, and if the hearing aid processing is assumed linear, the processor output x will be a sinusoid. The acoustic feedback signal, f will be a sinusoid. The incoming sound, v, and the acoustic feedback will be blended (summed), which yield another sinusoid (amplitude and phase altered), etc. 
         [0069]    The adaptive feedback cancellation filter {circumflex over (F)} relies on the processor output y as reference signal and produces output signal {circumflex over (f)}. The cancellation filter output signal {circumflex over (f)} is subtracted from the microphone output y to yield processor input signal e. 
         [0070]    If, in this case, one of the filter adaptation algorithms shown in Eqs. 7-10 is used to adjust the coefficients in the feedback cancellation filter {circumflex over (F)}, the cancellation filter will attempt to cancel y as this signal can be described as x with a simple change in amplitude and phase. The problem is that this is not the goal. The goal is to achieve that {circumflex over (f)}=f; not to remove tonal components in the environment. This example illustrates that if the external sound, v, is somehow “predictable”, one can expect large errors in the coefficients of the adaptive feedback cancellation filter. The present invention suggest to cope with this problem by providing a method according to which the adaptation will be halted if it is detected that an external tone is played as will be described in more detail below. 
         [0071]    It has been further observed in relation to the example above that a gain in the hearing aid processor, H, plays an important role for the accuracy of the feedback cancellation. If H represents a small amplification gain, the amplitude of the sinusoid, x, is small compared to the sinusoid, y, because only the amplitude of the feedback signal, f, is affected by the gain; not the incoming sinusoid, v. The reverse is the case when the gain is large. If the cancelling filter adaptation runs, the coefficients in {circumflex over (F)} are adjusted to make {circumflex over (f)} cancel the signal y. The error in the coefficients will consequently increase with a decreasing gain in the hearing aid processor. This is well in line with the result derived below with reference to Eq. 17. 
         [0072]    Generally, it has been observed that the more the signal x resembles a sinusoid with the less accuracy will the cancellation filter model the acoustic feedback (and instead attempt to attenuate the tone). This is a challenge because instability in the hearing aid will typically manifest itself as howling; a periodic signal resembling a tone. According to the present invention, there are at least two approaches provided which, at a first glance, seem to be completely contradictory: If an external tone is played, it is suggested to stop adaptation (μ=0) as otherwise the filter will be misadjusted; if a tone is generated internally due to feedback, it is to adapt fast in order to quickly compensate the tone. 
         [0073]    In the patent application WO 01/06812, a procedure is described, where an adaptive resonator filter is used for detecting whether a dominating tone is present. If it is, fast adaptation is used for attenuating the tone. This is an efficient procedure for eliminating feedback howling, but it will obviously produce severe artefacts when tones or whistling sounds are present in the environment. 
         [0074]    According to an embodiment of the present invention, another approach to cope with this problem is implemented by reducing the adaptation rate when the sound is spectrally coloured. This will reduce the ability to cancel feedback howling, so, according to a particular embodiment, the reduction of the adaptation rate is used along with a system for stabilizing the closed-loop system by limiting the amplification, thereby stopping the howling. 
         [0075]    Generally, modern hearing aids use compression for compensating the hearing-loss. Thus, the amplification in the hearing aid processor is decreased with increasing input sound levels. Without an anti-feedback system, the hearing aid processor will thus in worst case make the closed-loop system marginally stable; i.e., the level of the feedback howling will eventually be constant. To cope with this problem, according to an embodiment of the present invention, if feedback howling is observed then a small decrease in the amplification gain is applied which will stabilize the closed-loop system, resulting in removal of the howling. When the howling is removed, it is again safe to adapt the cancelling filter and eventually the filter will model the acoustic feedback better. This will in turn allow headroom for an increase in the amplification gain. 
         [0076]    Further approaches suggesting to stabilize the closed-loop system are disclosed in WO 02/25996, which provides a method for suppressing the time varying acoustic feedback with an adaptive filter, and co-pending patent application, filed on Mar. 31, 2006 with the title “Hearing aid and method of estimating dynamic gain limitation in a hearing aid”, PCT/EP2006/061215, published as WO2007112777, which provides an acoustic loop gain estimator for determining a dynamic gain limit, and which is herewith incorporated by reference. 
         [0077]    Rather than using a tone detector as described in WO 01/06812, according to an embodiment of the present invention, there is provided a method and a hearing aid using measures of either autocorrelation of the signal or one of the similar quantities as described in the previously mentioned co-pending patent application WO2007113283, “Method for controlling signal processing in a Hearing aid and a Hearing aid implementing this method” to detect whether an external tone is present. 
         [0078]    According to further embodiments of the present invention, the mentioned problems with spectral colouring can to some extent be further alleviated by the use of either adaptive notch filters to attenuate tones and/or by adaptive whitening filters to produce a spectral flattening of the signals. 
         [0079]    Since it is a complex issue to decide how the adaptation step size should optimally depend on the measure of signal autocorrelation, the present invention provides several methods and hearing aids, which at a first glance might be seen as following to some extent different and contradictory approaches, and which will be described now in more detail. 
         [0080]    According to an embodiment of the present invention, the step size of the feedback cancelling filter in a hearing aid is set in dependency of the autocorrelation value of the compensated signal e in  FIG. 5 . According to an embodiment, the cancelling filter is an FIR filter adjusted according to Eq. 8 or Eq. 9. According to a particular embodiment, an adaptive whitening filter is applied on the reference signal (and a similar filter is applied to the adaptation error). The step size is set according to the following formula resulting in a fast cancellation of tones for which the autocorrelation calculation gives a maximum correlation coefficient value &gt;0.98 so that a fast adaptation rate is applied. 
         [0000]    μ fast : A large step-size (fast adaptation rate).
 
μ slow : A small step-size (slow adaptation rate).
 
         [0000]    
       
         
           
             
               
                 r 
                 e 
               
                
               
                 ( 
                 τ 
                 ) 
               
             
             = 
             
               
                 
                   
                     1 
                     N 
                   
                    
                   
                     
                       ∑ 
                       
                         k 
                         = 
                         τ 
                       
                       N 
                     
                      
                     
                       
                         e 
                         k 
                       
                        
                       
                         e 
                         
                           k 
                           - 
                           τ 
                         
                       
                     
                   
                 
                 
                   
                     σ 
                     ^ 
                   
                   e 
                   2 
                 
               
                
               
                 : 
               
             
           
         
       
     
         [0000]    Autocorrelation coefficients based on the compensated signal. 
         [0000]    
       
         
           
             
               r 
               max 
             
             = 
             
               
                 
                   max 
                   τ 
                 
                  
                 
                   { 
                   
                      
                     
                       
                         r 
                         e 
                       
                        
                       
                         ( 
                         τ 
                         ) 
                       
                     
                      
                   
                   } 
                 
               
               : 
             
           
         
       
     
         [0000]    Maximum correlation coefficient.
 
A procedure for adjustment of the step size is:
 
         [0000]      If r max &gt;0.98 Then 
         [0000]      μ k =μ fast    
         [0000]      Else 
         [0000]      μ k =μ slow . 
         [0081]    According to another embodiment, the step size is decreased according to a monotonous function with increased autocorrelation of the reference signal. This embodiment allows to reduce the step size with increasing spectral colouring. 
         [0082]    According to an embodiment, the cancelling filter is an FIR filter adjusted according to Eq. 8 or Eq. 9. According to a particular embodiment, an adaptive whitening filter is applied on the reference signal (and a similar filter is applied to the adaptation error). The step size is decreased according to the following procedure for increasing maximum correlation coefficients in order to prevent the onset of undesired oscillation due to a distortion of the model of the feedback path modelled by the feedback cancelling filter coefficients. According to particular embodiments, an initiated feedback oscillation will be handled by further measures. The procedure is as follows: 
         [0000]    μ 1 , μ 2 , μ max : step-sizes of increasing magnitude, 0&lt;μ 1 &lt;μ 2 &lt;μ max &lt;2
 
T max , T 1 , T 2 : Autocorrelation thresholds of decreasing magnitude, 1&gt;T max &gt;T 1 &gt;T 2 &gt;0.
 
         [0000]    
       
         
           
             
               
                 r 
                 e 
               
                
               
                 ( 
                 τ 
                 ) 
               
             
             = 
             
               
                 
                   
                     1 
                     N 
                   
                    
                   
                     
                       ∑ 
                       
                         k 
                         = 
                         τ 
                       
                       N 
                     
                      
                     
                       
                         e 
                         k 
                       
                        
                       
                         e 
                         
                           k 
                           - 
                           τ 
                         
                       
                     
                   
                 
                 
                   
                     σ 
                     ^ 
                   
                   e 
                   2 
                 
               
                
               
                 : 
               
             
           
         
       
     
         [0000]    Autocorrelation coefficients. 
         [0000]    
       
         
           
             
               r 
               max 
             
             = 
             
               
                 max 
                 τ 
               
                
               
                 
                   { 
                   
                      
                     
                       
                         r 
                         e 
                       
                        
                       
                         ( 
                         τ 
                         ) 
                       
                     
                      
                   
                   } 
                 
                  
                 
                   : 
                 
               
             
           
         
       
     
         [0000]    Maximum correlation coefficient. 
         [0083]    According to the procedure, the step size is adjusted as follows: 
         [0000]      If r max &gt;T max  Then μ k =0 
         [0000]      Else If r max &gt;T 1  Then μ k =μ 1    
         [0000]      Else If r max &gt;T 2  Then μ k =μ 2    
         [0000]      Else μ k =μ max    
         [0084]    The embodiments described above can be varied in numerous ways. As most hearing aids operate in a number of frequency bands, the autocorrelation coefficients are calculated in several bands separately according a particular embodiment. In this way it is often easier to detect if spectral colouring occurs locally. The procedure is as follows: 
         [0000]    
       
         
           
             
               
                 r 
                 e 
                 
                   ( 
                   i 
                   ) 
                 
               
                
               
                 ( 
                 τ 
                 ) 
               
             
             = 
             
               
                 
                   
                     1 
                     N 
                   
                    
                   
                     
                       ∑ 
                       
                         k 
                         = 
                         τ 
                       
                       N 
                     
                      
                     
                       
                         e 
                         k 
                         
                           ( 
                           i 
                           ) 
                         
                       
                        
                       
                         e 
                         
                           k 
                           - 
                           τ 
                         
                         
                           ( 
                           i 
                           ) 
                         
                       
                     
                   
                 
                 
                   
                     ( 
                     
                       
                         σ 
                         ^ 
                       
                       e 
                       
                         ( 
                         i 
                         ) 
                       
                     
                     ) 
                   
                   2 
                 
               
                
               
                 : 
               
             
           
         
       
     
         [0000]    Autocorrelation coefficients. (i) is an index over bands, i={1, . . . , B}
 
and redefine
 
         [0000]    
       
         
           
             
               
                 r 
                 max 
               
                
               
                 
                   max 
                   i 
                 
                  
                 
                   
                     max 
                     τ 
                   
                    
                   
                     { 
                     
                        
                       
                         
                           r 
                           x 
                           
                             ( 
                             i 
                             ) 
                           
                         
                          
                         
                           ( 
                           τ 
                           ) 
                         
                       
                        
                     
                     } 
                   
                 
               
             
             : 
           
         
       
     
         [0000]    B. The coefficient over the bands is then used to adjust the step size as explained above. 
       Gain Dependency 
       [0085]    The description of embodiments of the present invention taking gain dependency into account is based on the derivations in Section 9.4 in S. Haykin:  Adaptive Filter Theory,  3 rd  Edition, Prentice-Hall, NJ, USA, 1996. It is advised to consult this book for intermediate results and further description of assumptions. 
         [0000]    First the following quantities are introduced:
 
ŵ k : Estimated weight vector at sample k.
   w : Optimum Wiener solution for coefficients in the cancelling filter (i.e., the true coefficients provided that the filter structure is sufficiently flexible to describe the acoustic feedback).
 
J k ≡E{e k   2 }: The mean squared error at sample k.
 
J min ≡E{ē k   2 }: The mean squared error evaluated in the Wiener solution. Assuming as above that the Wiener solution for the coefficients corresponds to the true acoustic feedback path then J min =E{v k   2 }.
 
ε k ≡  w −ŵ k : Coefficient error vector; the error between estimated and “true” coefficients.
 
K k ≡E{ε k ε k   T }: Correlation matrix for the coefficient error vector.
 
         [0086]    Furthermore, the assumption is made that the reference signal, x k , is white. In most practical sound environments this is not a valid assumption, but it can be achieved through the use of an adaptive whitening filter. According to an embodiment, the output signal x of the hearing aid processor H is input to the adaptive whitening filter (not shown in  FIGS. 4 and 5 ) and the output of the adaptive whitening filter is input to the adaptive cancelling filter. 
         [0087]    Consider first the setup shown in  FIG. 4  in which the compensated microphone input is multiplied by a simple gain, G, to produce x k . If x k  is assumed white then the environmental signal, v k , is also white. As mentioned; the whitening occurs as a consequence of an adaptive whitening filtering according to a particular embodiment. Further, the following definitions are made: 
         [0000]    R x =E{x k x k   T }=σ 2 I: is the correlation matrix for the reference signal.
 
R v =E{v k v k   T }=σ v   2 I: is the correlation matrix for the incoming signal. This equals J min  under the assumption that the cancelling filter length is sufficient.
 
         [0088]    According to S. Haykin:  Adaptive Filter Theory,  3 rd  Edition, Prentice-Hall, NJ, USA, 1996, the correlation matrix for the coefficient error vector in an LMS-algorithm develops according to 
         [0000]        K   k =( I−μR   x ) K   k-1 ( I−μR   x )+μ 2   J   min   R   x   [Eq. 11] 
         [0000]    Specializing this to white noise reference signals, R x =σ 2 I, gives 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           K 
                           k 
                         
                         = 
                           
                          
                         
                           
                             
                               ( 
                               
                                 I 
                                 - 
                                 
                                   μ 
                                    
                                   
                                       
                                   
                                    
                                   
                                     σ 
                                     2 
                                   
                                    
                                   I 
                                 
                               
                               ) 
                             
                              
                             
                               
                                 K 
                                 
                                   k 
                                   - 
                                   1 
                                 
                               
                                
                               
                                 ( 
                                 
                                   I 
                                   - 
                                   
                                     μ 
                                      
                                     
                                         
                                     
                                      
                                     
                                       σ 
                                       2 
                                     
                                      
                                     I 
                                   
                                 
                                 ) 
                               
                             
                           
                           + 
                           
                             
                               μ 
                               2 
                             
                              
                             
                               J 
                               min 
                             
                              
                             
                               σ 
                               2 
                             
                              
                             I 
                           
                         
                       
                     
                   
                   
                     
                       
                         = 
                           
                          
                         
                           
                             
                               
                                 ( 
                                 
                                   1 
                                   - 
                                   
                                     μ 
                                      
                                     
                                         
                                     
                                      
                                     
                                       σ 
                                       2 
                                     
                                   
                                 
                                 ) 
                               
                               2 
                             
                              
                             
                               K 
                               
                                 k 
                                 - 
                                 1 
                               
                             
                           
                           + 
                           
                             
                               μ 
                               2 
                             
                              
                             
                               J 
                               min 
                             
                              
                             
                               σ 
                               2 
                             
                              
                             I 
                           
                         
                       
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     12 
                   
                   ] 
                 
               
             
           
         
       
     
         [0000]    or in steady state 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           
                             ( 
                             
                               1 
                               - 
                               
                                 
                                   ( 
                                   
                                     1 
                                     - 
                                     
                                       μ 
                                        
                                       
                                           
                                       
                                        
                                       
                                         σ 
                                         2 
                                       
                                     
                                   
                                   ) 
                                 
                                 2 
                               
                             
                             ) 
                           
                            
                           
                             K 
                             ∞ 
                           
                         
                         = 
                           
                          
                         
                           
                             μ 
                             2 
                           
                            
                           
                             J 
                             min 
                           
                            
                           
                             σ 
                             2 
                           
                            
                           I 
                         
                       
                     
                   
                   
                     
                       
                         
                             
                         
                         ⇓ 
                                                                           
                       
                     
                   
                   
                     
                       
                         
                           K 
                           ∞ 
                         
                         = 
                           
                          
                         
                           
                             
                               
                                 μ 
                                 2 
                               
                                
                               
                                 J 
                                 min 
                               
                                
                               
                                 σ 
                                 2 
                               
                             
                             
                               
                                 2 
                                  
                                 μ 
                                  
                                 
                                     
                                 
                                  
                                 
                                   σ 
                                   2 
                                 
                               
                               - 
                               
                                 
                                   μ 
                                   2 
                                 
                                  
                                 
                                   σ 
                                   4 
                                 
                               
                             
                           
                            
                           I 
                         
                       
                     
                   
                   
                     
                       
                         = 
                           
                          
                         
                           
                             
                               
                                 J 
                                 min 
                               
                                
                               μ 
                                
                               
                                   
                               
                                
                               
                                 σ 
                                 2 
                               
                             
                             
                               
                                 2 
                                  
                                 
                                   σ 
                                   2 
                                 
                               
                               - 
                               
                                 μ 
                                  
                                 
                                     
                                 
                                  
                                 
                                   σ 
                                   4 
                                 
                               
                             
                           
                            
                           I 
                         
                       
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     13 
                   
                   ] 
                 
               
             
           
         
       
     
         [0089]    To simplify this, the LMS with variance normalization, which has a behaviour similar to that of the NLMS-algorithm, is used according to an embodiment. A more formal treatment relating to NLMS can be found in D. T. M Slock:  On the Convergence Behavior of the LMS and the Normalized LMS Algorithms , IEEE Trans. Signal Processing, Vol. 41, No. 9, September 1993, pp. 2811-2824. According to the embodiment, the step size is normalized with the exact variance of the reference signal; that is, the step size 
         [0000]    
       
         
           
             
               
                 
                   
                     μ 
                     = 
                     
                       
                         μ 
                         _ 
                       
                       
                         σ 
                         2 
                       
                     
                   
                   , 
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     14 
                   
                   ] 
                 
               
             
           
         
       
     
         [0000]    is inserted in the above: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           K 
                           ∞ 
                         
                         = 
                           
                          
                         
                           
                             
                               
                                 J 
                                 min 
                               
                                
                               
                                 μ 
                                 _ 
                               
                             
                             
                               
                                 2 
                                  
                                 
                                   σ 
                                   2 
                                 
                               
                               - 
                               
                                 
                                   μ 
                                   _ 
                                 
                                  
                                 
                                   σ 
                                   2 
                                 
                               
                             
                           
                            
                           I 
                         
                       
                     
                   
                   
                     
                       
                         
                           ≈ 
                             
                            
                           
                             
                               
                                 
                                   J 
                                   min 
                                 
                                  
                                 
                                   μ 
                                   _ 
                                 
                               
                               
                                 2 
                                  
                                 
                                   σ 
                                   2 
                                 
                               
                             
                              
                             I 
                           
                         
                         , 
                       
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     15 
                   
                   ] 
                 
               
             
           
         
       
     
         [0000]    J min  is not available, but instead an estimate of it is 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       used 
                        
                       
                         : 
                       
                        
                       
                           
                       
                        
                       
                         σ 
                         e 
                         2 
                       
                     
                     = 
                     
                       
                         E 
                          
                         
                           { 
                           
                             e 
                             k 
                             2 
                           
                           } 
                         
                       
                       = 
                       
                         
                           
                             σ 
                             2 
                           
                           
                             G 
                             2 
                           
                         
                         . 
                         
                           
 
                         
                          
                         Thus 
                       
                     
                   
                   , 
                   
                     
 
                   
                    
                   
                     
                       K 
                       ∞ 
                     
                     ≈ 
                     
                       
                         
                           μ 
                           _ 
                         
                         
                           2 
                            
                           
                             G 
                             2 
                           
                         
                       
                        
                       I 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     16 
                   
                   ] 
                 
               
             
           
         
       
     
         [0000]    or, if the uncertainty on the individual filter coefficients is considered: 
         [0000]    
       
         
           
             
               
                 
                   
                     δ 
                      
                     
                         
                     
                      
                     
                       w 
                       i 
                     
                   
                   = 
                   
                     
                       
                         K 
                         ∞ 
                         
                           ( 
                           
                             i 
                             , 
                             i 
                           
                           ) 
                         
                       
                     
                     ≈ 
                     
                       
                         
                           
                             
                               μ 
                               _ 
                             
                             / 
                             2 
                           
                         
                         G 
                       
                       . 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     17 
                   
                   ] 
                 
               
             
           
         
       
     
         [0090]    This result shows that if it is desired to maintain a specific uncertainty on the filter coefficients, the step size should be reduced by Δ 2  every time the gain is reduced by a factor Δ. 
         [0000]    In an embodiment, which is more relevant for a modern hearing aid, a bandsplit filter on the signal e in  FIG. 4  is used to generate a number of overlapping frequency bands, {e k   (1) , e k   (2) , . . . , e k   (B) }. On each of these bands, a separate amplification gain {G (1) , G (2) , . . . , G (B) } is used before the bands are added together to produce the signal x k . In order to ensure a certain maximum uncertainty on the filter coefficients, a safe approach is to scale the step size in accordance with changes in the smallest of the gains {G (1) , G (2) , . . . , G (B) }. 
       Amplification in the Hearing Aid Processor 
       [0091]    In the following, embodiments will be described which deal with amplification in the hearing aid processor. The resulting amplification in the hearing aid processor is usually composed of the output of various subsystems, such as a compression unit for compensating the hearing-loss, a temporal noise reduction system for attenuating unwanted noise, automatic gain control and more. Most often, these various systems operate in a number of frequency bands and separate gains are assigned to each band. In some hearing aids, the hearing aid processor is an adaptive wide-band filter and a mechanism is incorporated for adjusting the filter so that the amplitude response varies in accordance with the current sound pressure levels in a number of frequency bands. 
         [0092]    According to an embodiment, it is assumed that one of the algorithms NLMS in Eq. 8 or LMS with variance normalization in Eq. 9 is employed for adapting coefficients in the feedback cancelling filter and that the step size is constant. An important lesson learned from Eq. 17 is that if the amplification gain of the hearing aid processor is varied slowly compared to the adaptation rate, the stability margin will be more or less constant. If the amplification gain is increased, the cancelling filter becomes equally more accurate and vice versa. In most hearing aids, the amplification gain is, however, adjusted rapidly in comparison to the possible adaptation rate in the cancelling filter. Thus, if there has been a period of time with a small amplification gain, the accuracy of the cancelling filter is decreased. If suddenly the amplification goes up, the closed-loop system can become unstable. 
         [0093]    According to an embodiment, this problem is solved by providing higher accuracy when the hearing aid amplification is small. Thus, when the amplification goes down, the step size, μ, is reduced and vice versa. Following Eq. 17, a nominal step size is selected, which provides the desired accuracy at the maximum amplification gain, and then the step size is reduced proportional to the square of reductions in the amplification gain. 
         [0094]    According to another embodiment, the hearing aid processor corresponds to a simple amplification gain. The cancelling filter is an FIR filter adjusted according to Eq. 8 or Eq. 9 and an adaptive whitening filter is applied on the reference signal. According to a particular embodiment, a similar filter is applied to the adaptation error. It is: 
         [0000]    μ max : The maximum step-size (fastest adaptation rate).
 
G max : The maximum amplification gain used in the hearing aid processor. The maximum gain can be set according to the hearing-loss or according to an estimate of the stability limit (over which the hearing aid will howl).
 
G k : Current amplification gain.
 
With reference to Eq. 17, the step-size at sample number k is calculated as
 
         [0000]    
       
         
           
             
               
                 
                   
                     μ 
                     k 
                   
                   = 
                   
                     
                       
                         ( 
                         
                           
                             G 
                             k 
                           
                           
                             G 
                             max 
                           
                         
                         ) 
                       
                       2 
                     
                      
                     
                       μ 
                       max 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     18 
                   
                   ] 
                 
               
             
           
         
       
     
         [0000]    This step size is then used in a method or hearing aid providing a wide band solution. 
         [0095]    According to an embodiment providing a multi-band solution, in a multi-band hearing aid the signal is split into a number of frequency bands and an amplification gain is applied to each band before summing the bands. A conservative step-size control for this application is given below. 
         [0000]    G max,i : The maximum amplification gain used in the hearing aid processor for band i. The maximum can be set according to the hearing-loss or according to an estimate of the stability limit (over which the hearing aid will howl).
 
G i,k : Current amplification gain used in band i.
 
With reference to Eq. 17 and assuming we are operating with B frequency bands, the step-size at sample number k is calculated as
 
         [0000]    
       
         
           
             
               
                 
                   
                     μ 
                     k 
                   
                   = 
                   
                     
                       
                         ( 
                         
                           Min 
                            
                           
                             { 
                             
                               
                                 
                                   G 
                                   
                                     1 
                                     , 
                                     k 
                                   
                                 
                                 
                                   G 
                                   
                                     max 
                                     , 
                                     1 
                                   
                                 
                               
                               , 
                               
                                 
                                   G 
                                   
                                     2 
                                     , 
                                     k 
                                   
                                 
                                 
                                   G 
                                   
                                     max 
                                     , 
                                     2 
                                   
                                 
                               
                               , 
                               … 
                                
                               
                                   
                               
                               , 
                               
                                 
                                   G 
                                   
                                     B 
                                     , 
                                     k 
                                   
                                 
                                 
                                   G 
                                   
                                     max 
                                     , 
                                     B 
                                   
                                 
                               
                             
                             } 
                           
                         
                         ) 
                       
                       2 
                     
                      
                     
                       μ 
                       max 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     19 
                   
                   ] 
                 
               
             
           
         
       
     
       Adaptation Halt 
       [0096]    Sudden loud sounds, such as a door slamming or a hammer like sound, impose special risks when the cancelling filter is updated with an NLMS-like algorithm. The hearing aid processor will typically delay the signal, as most often it includes a filter bank, an FFT and/or other types of filters. This means that a sudden loud sound will quickly manifest itself in the adaptation error (e) in  FIG. 5 , but not until later on the reference for the cancellation filter (x). Therefore, the NLMS update as described in Eq. 8 will take very large adaptation steps right after the loud sound occurs because the denominator in Eq. 8 is small and the error signal is large. Moreover, it is adaptation steps, which are not governed by discrepancies between cancellation filter and acoustic feedback path. 
         [0097]    According to the invention, methods and hearing aids are provided to detect if a sudden increase in sound pressure occurs and temporarily suspend the adaptation afterwards. An embodiment of this is depicted in  FIG. 6  and will now be described. 
         [0098]    The input to the mechanism, which is part of a hearing aid, is for example the microphone signal  601  or an omnidirectional signal of the hearing aid. According to a particular embodiment, this signal is filtered. If, e.g., the feedback cancellation filter is implemented according to an embodiment so that it works in the high-frequency range only, it is not of much relevance what happens at lower frequencies. Thus, in order to detect sudden loud sounds with high-frequency components, the frequency weighting filter  602  could be a high-pass filter. The absolute value of the signal X is then taken by Abs-block  603  and this operation is then followed by a sliding averaging in averager  604  or some other type of magnitude calculation. The average of absolute values, Z, reflects the current sound pressure. The time-constant or window size in the average should at least correspond to the delay in the hearing aid processor and the length of the feedback cancelling filter. To detect if a loud sound occurs, the average signal Z is increased by a great amount, which is defined by a constant Threshold to get a signal A, which is then compared in block  606  to the momentary signal magnitude. If the momentary signal magnitude exceeds the signal A, the sound is classified as “a sudden loud sound”. In order to suspend the adaptation for a while after this happens, one solution is to use a peak holding block  605  applied on Y, which can store information about the signal maximum for a while after it occurred as signal B. If by the comparison of signals A and B in comparator  606  it is detected that A&lt;B, the adaptation is suspended by sending an adapt_disable signal  607 . 
         [0099]    Loud sounds (not necessarily sudden) can also cause a nonlinear behavior in one or more components of the hearing aid. The acoustic feedback path as it is seen from the cancelling filter&#39;s perspective embraces microphone(s), receiver and input- and output converters. Saturation or overload in one of these units thus corresponds to a non-linearity in the acoustic feedback path. Assuming a linear filter is used for feedback cancellation (such as an FIR filter), the filter is inadequate for modelling the highly nonlinear saturation function, thus leading to errors in the adaptation. Therefore, according to an embodiment, a detector (not shown) for recognition of these circumstances is included in the adaptation mechanism, and adaptation of the cancellation filter is temporarily suspended when the non-linearity occurs. The adaptation may, according to a particular embodiment, be suspended for a short while after one circumstance of that kind has been detected. 
       Dependency on Directional System—Calculating the Efficiency of a Spatial Filter 
       [0100]    The most advanced hearing aids today are supplied with directional microphones, with two or more omnidirectional microphones, or with a combination of omnidirectional and directional microphones. A directional microphone is a special microphone, which has two inlets and works according to the “delay-and-subtract” principle. Such a microphone will provide a signal, which has a fixed directional pattern. A directional system based on two or more omnidirectional microphones allows for an adaptive directional pattern and can also be extended to work in several frequency bands to enable a frequency dependent directional pattern. See for example patent application WO 01/01731 A1. In any case, spatial filtering is a highly efficient means of increasing the signal-to-noise ratio in many typical listening situations. An example of such a system is shown in  FIG. 7 . 
         [0101]    To determine the efficiency of a directional system at a given point in time it is useful to compare an estimated norm of the signals before and after the directional system. One can use the wide-band signal to get an estimate of the overall efficiency or number of band-pass filtered signals to get an estimate of the efficiency over frequency. 
         [0102]    Many norms can be considered and for practical use one will employ an approximation to reflect the value relevant in a window around the current point in time. The general p-norm definition along with some special cases of it is shown in [Eq. 20] and Table 1. 
         [0000]    The p-norm of a signal over some window is defined as: 
         [0000]    
       
         
           
             
               
                 
                   
                     N 
                     x 
                   
                   = 
                   
                     
                       
                          
                         x 
                          
                       
                       p 
                     
                     = 
                     
                       
                         ( 
                         
                           
                             ∑ 
                             
                               k 
                               = 
                               0 
                             
                             M 
                           
                            
                           
                             
                               F 
                               k 
                             
                              
                             
                               
                                  
                                 
                                   x 
                                   k 
                                 
                                  
                               
                               p 
                             
                           
                         
                         ) 
                       
                       
                         1 
                         / 
                         p 
                       
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     20 
                   
                   ] 
                 
               
             
           
         
       
     
         [0000]    {F k } represents a window or filter function. Various applicable norms are shown in Table 1 (shown with a rectangular window function of size M): 
         [0000]    
       
         
               
             
               
               
               
             
           
               
                 TABLE 1 
               
               
                   
               
               
                 Norm computation 
               
               
                   
               
             
             
               
                   
               
             
          
           
               
                   
                 1-norm 
                 ∥x∥ 1  = |x 1 |+ . . . + |x M | 
               
               
                   
                 Euclidean 
                 ∥x∥ 2  = {square root over (x 1   2  + .  .  .  + x M   2 )} 
               
               
                   
                 General 
                 ∥x∥ p  = (|x 1 | p  + . . . + |x M | p ) 1/p  for 1 ≦ p ≦ ∞ 
               
               
                   
                 Infinity 
                 ∥x∥ ∞ = max{|x 1 |, . . . , |x M |} 
               
               
                   
                 −Infinity 
                 ∥x∥ −∞ = min{|x 1 |, . . . , |x M |} 
               
               
                   
                   
               
             
          
         
       
     
         [0103]    A commonly used norm calculation within this category is based on the 1-norm. At sampling instant k, the norm is calculated by the recursive update with exponential forgetting: 
         [0000]        N   x ( k )=φ| x   k |+(1−φ)· N   x ( k− 1)  [Eq. 21] 
         [0000]    where φ is a constant, φε]0;1] (by this update the norm is also normalized to make it independent of window length). 
         [0104]    If N x  is the norm of an input signal, x, and N y  is the norm of an output signal, y, then the efficiency of the directional system in the frequency band to which x and y belongs can be calculated as 
         [0000]    
       
         
           
             
               
                 
                   G 
                   = 
                   
                     
                       N 
                       y 
                     
                     
                       N 
                       x 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     22 
                   
                   ] 
                 
               
             
           
         
       
     
         [0105]    If G is near 0, the directional system is highly efficient and is most likely removing a significant amount of noise or irrelevant signal components. 
         [0000]    Interaction with Multi-Microphone or Directional Microphone Systems 
         [0106]    A directional system for spatially filtering of the sound can be considered as a gain applied to the sound. Depending on the directional pattern selected and the location of the individual sound sources, this “gain” will take different values. Under fortunate circumstances a directional system can reduce the feedback problems, but generally one will not have exact knowledge of the sound source locations. When considering the directional system as a gain, it has been observed that in multi-microphone implementations like those depicted in  FIG. 10  and  FIG. 8 , the formula Eq. 17 plays a role for the accuracy of the feedback cancelling filter. 
         [0107]    The overall change of amplification gain due to the directional system can be calculated according to Eq. 21 and Eq. 22. 
         [0108]    According to an embodiment, Eq. 17 is used to govern the step size control. An implementation according to this embodiment will be described in the following with reference to  FIG. 8 . 
         [0109]      FIG. 8  shows a hearing aid with directional characteristics. The cancelling filters are FIR filters adjusted according to Eq. 8 or Eq. 9 and an adaptive whitening filter is applied on the reference signal. According to a particular embodiment, a similar filter is applied to the adaptation errors. The following definitions are made: 
         [0000]    N 1,k : The norm of the first spatial signal  32 . The norm is estimated according to Eq. 21.
 
N 2,k : The norm of the second spatial signal  33 . The norm is estimated according to Eq. 21.
 
P k : The norm of the resulting directional signal  34 . The norm is estimated according to Eq. 21.
 
         [0000]    
       
         
           
             
               G 
               
                 1 
                 , 
                 k 
               
             
             = 
             
               
                 
                   P 
                   k 
                 
                 
                   N 
                   
                     1 
                     , 
                     k 
                   
                 
               
               : 
             
           
         
       
     
         [0000]    Reduction of the first spatial signal  32  occurring in the directional weighting system  205 . 
         [0000]    
       
         
           
             
               G 
               
                 2 
                 , 
                 k 
               
             
             = 
             
               
                 
                   P 
                   k 
                 
                 
                   N 
                   
                     2 
                     , 
                     k 
                   
                 
               
               : 
             
           
         
       
     
         [0000]    Reduction of the second spatial signal  33  occurring in the directional weighting system  205 .
 
μ max : The maximum step-size (fastest adaptation rate).
 
         [0110]    To keep an upper limit on the accuracy of the cancelling filter, according to an embodiment changes of the step size are made by using Eq. 17. For sample k the step sizes used in the two feedback cancelling filters are then calculated as 
         [0000]      μ 1,k =G 1,k   2 μ max   [Eq. 23] 
         [0000]      μ 2,k =G 2,k   2 μ max   [Eq. 24] 
         [0000]    According to another embodiment, a multi-band directional system is used. If the signals  32  and  33  in  FIG. 8  are split into several frequency bands before being weighted together to achieve a further noise reduction compared to what is possible using a weighting of the broad-band signals, the gain reductions defined above must be calculated for each frequency band. A step size parameter can then be calculated for each band. The safest approach is then to take the minimum step size for each of the two branches and use these in the feedback cancelling filters: 
         [0000]      μ 1,k =Min{μ 1,k   (1) ,μ 1,k   (2) , . . . , μ 1,k   (B) }  [Eq. 25] 
         [0000]      μ 2,k =Min{μ 2,k   (1) ,μ 2,k   (2) , . . . , μ 2,k   (B) }  [Eq. 26] 
       FURTHER EMBODIMENTS 
       [0111]      FIGS. 8-12  show embodiments of hearing aid configurations including a subsystem for step size (adaptation rate) adjustment depicted as step size control block  104 ,  304  and  404 , which will be described in the following. 
         [0112]      FIG. 9  shows a hearing aid with one microphone like the one shown in  FIG. 2  except that the step size control block  104  has been introduced. The connection  7  symbolizes such information as amplification gains, state of automatic gain controller and noise reduction performance. The output  6  of block  104  is a step size parameter to be used in the adaptation block  103 . As it will appear in the following, the step size is set according to the output of the hearing aid processor  3 , the microphone signal  1  and the feedback cancelling signal  4 . 
         [0113]      FIG. 10  shows a hearing aid with two microphones and a separate feedback cancelling to each microphone signal. The compensated input signals  40 ,  41  are used as input to a spatial filtering system, which might be adaptive and work in multiple frequency bands. The resulting directional signal(s)  42  is (are) used as input to the hearing aid processor  100 . The filters  302   a ,  302   b  produce cancelling signals  43 ,  44  for each of the microphone signals  20 ,  21 . The adaptation of the cancelling filters takes place in adaptation block  303 , and outcome of this block is two sets of filter coefficients  46   a ,  46   b . The step size control block  304  works on parameters from the hearing aid processor  100 , one or both microphone signals, both cancelling filter outputs and the output of the hearing aid processor  100 . The step size control block  304  outputs one or two step size parameters  45   a ,  45   b . If both microphones are omnidirectional, the same step size parameter can be typically be used for adapting both cancelling filters. 
         [0114]      FIG. 11  shows a hearing aid with two omnidirectional microphones, a directional system for spatial noise filtering but only one feedback cancelling filter. This configuration is simpler than the one shown in  FIG. 10 , but the directional system becomes part of the acoustic feedback loop as it is seen from the perspective of the feedback cancelling filter. Thus, time-variations in the directional pattern require adaptation of the feedback cancelling filter coefficients. 
         [0115]      FIG. 12  shows a configuration similar to the one depicted in  FIG. 3 , but with the addition of a step size control Block  404 . This block provides two separate step size parameters  37   a ,  37   b  to be used for adaptation in block  403  of the coefficients  38   a ,  38   b  for each of the feedback cancelling filters  302   a ,  302   b . A consequence of using this concept as opposed to the one depicted in  FIG. 10 , is a highly different weighting of the adaptation error. Due to this difference, it is often easier to ensure stability of the hearing aid under the user of large amplification gains. 
         [0116]    In the following, further embodiments will be described which aim at providing an appropriate adaptation rate adjustment to remedy different adjustment problems. 
       Anti-Feedback Systems for Hearing Aids 
       [0117]    If one of the adaptation algorithms as defined in Eq. 7-Eq. 10 is used in a hearing aid like one of those depicted in  FIGS. 1-3  &amp;  8 - 12 , and the sound input represents a typical everyday sound environment, one will never achieve that the cancellation filter is an exact model of the acoustic feedback path. If an LMS-type adaptation algorithm is used with a constant step size, μ, the accuracy of the estimated feedback path will depend on several factors:
       1) The magnitude of the adaptation rate   2) The function and amplification in the hearing aid processor block  100 .   3) The “condition” of the microphone signal or signals; is the signal spectrally coloured or is it “noise-like”?   4) The performance of the multi-microphone directional system if such a system is integrated in the hearing aid.   5) The acoustic feedback path       
 
         [0123]    In order to make an accurate anti-feedback filter, the adaptation step size according to an embodiment is controlled in accordance with the items 2)-5). Further comments on each of the items mentioned will be given in the following along with a suggested adjustment of the step size parameter in each case. 
       Combining the Individual Effects 
       [0124]    Various observations about the signals entering the hearing aid and the state and behaviour of the hearing aid have been discussed above along with suggestions for adjusting the step size parameter accordingly. In the following, further embodiments will be described for how to combine the various effects into a single step size parameter for each feedback cancelling filter. 
         [0125]    At first, an embodiment of a hearing aid with directional system and a two-path feedback cancelling filter will be described with reference to  FIG. 12  depicting a hearing aid with a two-microphone implementation. According to a particular embodiment, the two feedback cancelling filters  302   a  and  302   b  are FIR-type filters, where the coefficients are adjusted using an adaptation block  403  such as LMS with variance normalization, as defined in Eq. 9, or an NLMS as defined in Eq. 8. The adaptation block  403 , according to an embodiment, contains an adaptive whitening filter which is applied on the reference signal  3  and the same filter is used on the adaptation errors, or, according to further embodiments, in a similar manner on signals  30 ,  31 ,  32 , and  33 . According to a particular embodiment, the hearing aid has B frequency bands and each band has a separate amplification gain and a separate directional pattern. The adaptation step size control unit  404  receives information about amplification gains from the hearing aid processor and band-splitted adaptation errors from either signals  51 ,  52  or, for simplicity, from signal  53 . The latter is used for calculating normalized autocorrelation or another type of self-similarity function for each band. It is further defined: 
         [0000]    N 1,k   (i) : The norm of the i&#39;th frequency band of the first spatial signal  51 .
 
The norm is estimated according to Eq. 21.
 
N 2,k   (i) : The norm of the i&#39;th frequency band of the second spatial signal  52 .
 
The norm is estimated according to Eq. 21.
 
P k   (i) : The norm of the i&#39;th frequency band of the resulting directional signal  53 . The norm is estimated according to Eq. 21.
 
         [0000]    
       
         
           
             
               G 
               
                 1 
                 , 
                 k 
               
               
                 ( 
                 i 
                 ) 
               
             
             = 
             
               
                 
                   P 
                   k 
                   
                     ( 
                     i 
                     ) 
                   
                 
                 
                   N 
                   
                     1 
                     , 
                     k 
                   
                   
                     ( 
                     i 
                     ) 
                   
                 
               
               : 
             
           
         
       
     
         [0000]    Reduction of the first spatial signal  51  occurring in the i&#39;th frequency band of the directional weighting system  205 . 
         [0000]    
       
         
           
             
               G 
               
                 2 
                 , 
                 k 
               
               
                 ( 
                 i 
                 ) 
               
             
             = 
             
               
                 
                   P 
                   k 
                   
                     ( 
                     i 
                     ) 
                   
                 
                 
                   N 
                   
                     2 
                     , 
                     k 
                   
                   
                     ( 
                     i 
                     ) 
                   
                 
               
               : 
             
           
         
       
     
         [0000]    Reduction of the second spatial signal  52  occurring in the i&#39;th frequency band of the directional weighting system  205 .
   G   k   (i) : The current amplification gain for band (i) as calculated in the hearing aid processor.
   G   max   (i) : The maximum amplification gain that can be used in the hearing aid processor. The maximum can be set according to the hearing-loss or according to an estimate of the stability limit (over which the hearing aid will howl).
 
         [0000]    
       
         
           
             
               
                 r 
                 e 
                 
                   ( 
                   i 
                   ) 
                 
               
                
               
                 ( 
                 τ 
                 ) 
               
             
             = 
             
               
                 
                   
                     1 
                     N 
                   
                    
                   
                     
                       ∑ 
                       
                         k 
                         = 
                         τ 
                       
                       N 
                     
                      
                     
                       
                         e 
                         k 
                         
                           ( 
                           i 
                           ) 
                         
                       
                        
                       
                         e 
                         
                           k 
                           - 
                           τ 
                         
                         
                           ( 
                           i 
                           ) 
                         
                       
                     
                   
                 
                 
                   
                     ( 
                     
                       
                         σ 
                         ^ 
                       
                       e 
                       
                         ( 
                         i 
                         ) 
                       
                     
                     ) 
                   
                   2 
                 
               
               : 
             
           
         
       
     
         [0000]    Autocorrelation coefficients for the i&#39;th band of the feedback compensated signal. τ 0 &lt;τ≦N. τ 0  is the standard transportation delay from the sound is send to the receiver until it is picked up by the microphone. N is the length of the tapped delay line used in the cancelling filters.
 
μ max : The maximum step-size (fastest adaptation rate).
 
For band i, calculate a step size decrement factor due to the amplification gain
 
         [0000]    
       
         
           
             
               
                 
                   
                     Δ 
                      
                     
                       
                         μ 
                         _ 
                       
                       k 
                       
                         ( 
                         i 
                         ) 
                       
                     
                   
                   = 
                   
                     
                       ( 
                       
                         
                           
                             G 
                             _ 
                           
                           k 
                           
                             ( 
                             i 
                             ) 
                           
                         
                         
                           
                             G 
                             _ 
                           
                           max 
                           
                             ( 
                             i 
                             ) 
                           
                         
                       
                       ) 
                     
                     2 
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     27 
                   
                   ] 
                 
               
             
           
         
       
     
         [0000]    and for each cancelling branch also a set of decrement factors due to the spatial filtering: 
         [0000]      Δμ 1,k   (i) =(G 1,k   (i) ) 2   [Eq. 28] 
         [0000]      Δμ 2,k   (i) =(G 2,k   (i) ) 2   [Eq. 29] 
         [0000]    Thus, a large decrement factor is equivalent to a small value Δμ. 
         [0126]    According to an embodiment, the autocorrelation coefficients in each frequency band are calculated from the feedback compensated inputs to the hearing aid processor. Then, a decrement factor is calculated in accordance with the maximum magnitude of the autocorrelation coefficients for each band (assuming the amplification gain is maximum): 
         [0000]    Δμ 1 , Δμ 2 : Decrement factors of decreasing magnitude, 0&lt;Δμ 1 &lt;Δμ 2 &lt;1
 
T max , T 1 , T 2 : Autocorrelation thresholds of decreasing magnitude, 1&gt;T max &gt;T 1 &gt;T 2 &gt;0.
 
         [0000]    
       
         
           
             
               
                 If 
                  
                 
                     
                 
                  
                 
                   
                     max 
                     τ 
                   
                    
                   
                     ( 
                     
                       
                         r 
                         k 
                         
                           ( 
                           i 
                           ) 
                         
                       
                        
                       
                         ( 
                         τ 
                         ) 
                       
                     
                     ) 
                   
                 
               
               &gt; 
               
                 
                   T 
                   1 
                 
                  
                 
                     
                 
                  
                 Then 
                  
                 
                     
                 
                  
                 Δ 
                  
                 
                     
                 
                  
                 
                   
                     μ 
                     ~ 
                   
                   k 
                   
                     ( 
                     i 
                     ) 
                   
                 
               
             
             = 
             
               Δ 
                
               
                   
               
                
               
                 μ 
                 1 
               
             
           
         
       
       
         
           
             
               
                 Else 
                  
                 
                     
                 
                  
                 If 
                  
                 
                     
                 
                  
                 
                   
                     max 
                     τ 
                   
                    
                   
                     ( 
                     
                       
                         r 
                         k 
                         
                           ( 
                           i 
                           ) 
                         
                       
                        
                       
                         ( 
                         τ 
                         ) 
                       
                     
                     ) 
                   
                 
               
               &gt; 
               
                 
                   T 
                   2 
                 
                  
                 
                     
                 
                  
                 Then 
                  
                 
                     
                 
                  
                 Δ 
                  
                 
                   
                     μ 
                     ~ 
                   
                   k 
                   
                     ( 
                     i 
                     ) 
                   
                 
               
             
             = 
             
               Δ 
                
               
                   
               
                
               
                 μ 
                 2 
               
             
           
         
       
     
         [0127]    The various decrement factors can be combined in different ways. According to a preferred embodiment, the step size decrement factors are compared within each band due to amplification gain and efficiency of the directional system, Δ  μ   k   (i) ·Δμ 1,k   (i) , to the step size decrement factors due to the colouring of the adaptation error: 
         [0000]    
       
         
           
             
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       μ 
                       
                         1 
                         , 
                         k 
                       
                     
                   
                   = 
                   
                     
                       min 
                       i 
                     
                      
                     
                       ( 
                       
                         min 
                          
                         
                           ( 
                           
                             
                               Δ 
                                
                               
                                   
                               
                                
                               
                                 
                                   
                                     μ 
                                     _ 
                                   
                                   k 
                                   
                                     ( 
                                     i 
                                     ) 
                                   
                                 
                                 · 
                                 Δ 
                               
                                
                               
                                   
                               
                                
                               
                                 μ 
                                 
                                   1 
                                   , 
                                   k 
                                 
                                 
                                   ( 
                                   i 
                                   ) 
                                 
                               
                             
                             , 
                             
                               
                                 
                                   
                                     Δ 
                                      
                                     
                                       
                                         
                                           μ 
                                           _ 
                                         
                                         k 
                                         
                                           ( 
                                           i 
                                           ) 
                                         
                                       
                                       · 
                                       Δ 
                                     
                                      
                                     
                                         
                                     
                                      
                                     
                                       μ 
                                       
                                         1 
                                         , 
                                         k 
                                       
                                       
                                         ( 
                                         i 
                                         ) 
                                       
                                     
                                   
                                 
                                 · 
                                 Δ 
                               
                                
                               
                                 
                                   μ 
                                   ~ 
                                 
                                 k 
                                 
                                   ( 
                                   i 
                                   ) 
                                 
                               
                             
                           
                           ) 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     30 
                   
                   ] 
                 
               
             
             
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       μ 
                       
                         2 
                         , 
                         k 
                       
                     
                   
                   = 
                   
                     
                       min 
                       i 
                     
                      
                     
                       ( 
                       
                         min 
                          
                         
                           ( 
                           
                             
                               Δ 
                                
                               
                                   
                               
                                
                               
                                 
                                   
                                     μ 
                                     _ 
                                   
                                   k 
                                   
                                     ( 
                                     i 
                                     ) 
                                   
                                 
                                 · 
                                 Δ 
                               
                                
                               
                                   
                               
                                
                               
                                 μ 
                                 
                                   2 
                                   , 
                                   k 
                                 
                                 
                                   ( 
                                   i 
                                   ) 
                                 
                               
                             
                             , 
                             
                               
                                 
                                   
                                     Δ 
                                      
                                     
                                       
                                         
                                           μ 
                                           _ 
                                         
                                         k 
                                         
                                           ( 
                                           i 
                                           ) 
                                         
                                       
                                       · 
                                       Δ 
                                     
                                      
                                     
                                         
                                     
                                      
                                     
                                       μ 
                                       
                                         2 
                                         , 
                                         k 
                                       
                                       
                                         ( 
                                         i 
                                         ) 
                                       
                                     
                                   
                                 
                                 · 
                                 Δ 
                               
                                
                               
                                 
                                   μ 
                                   ~ 
                                 
                                 k 
                                 
                                   ( 
                                   i 
                                   ) 
                                 
                               
                             
                           
                           ) 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     31 
                   
                   ] 
                 
               
             
           
         
       
     
         [0128]    As described previously, the error in the feedback cancelling filter will (in open-loop and for a fixed step size) be inverse proportional to the gain in the hearing aid processor. This dependency can be expressed by multiplying the decrement factors due to the colouring to the square root of the product of the two other types of decrement factor, as this square root is proportional to the decrement of the maximum amplification gain. Subsequent to these calculations, the largest decrement factor (smallest value) over bands is taken. The resulting step size for each branch is then 
         [0000]      μ 1,k =Δμ 1,k ·μ max   [Eq. 32] 
         [0000]      μ 2,k =Δμ 2,k ·μ max   [Eq. 33] 
         [0129]    According to an embodiment following a simpler, but quite conservative strategy, the decrements are multiplied within each band and subsequently take the factor leading to the largest decrement: 
         [0000]    
       
         
           
             
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       μ 
                       
                         1 
                         , 
                         k 
                       
                     
                   
                   = 
                   
                     
                       min 
                       i 
                     
                      
                     
                       ( 
                       
                         Δ 
                          
                         
                           
                             
                               μ 
                               _ 
                             
                             k 
                             
                               ( 
                               i 
                               ) 
                             
                           
                           · 
                           Δ 
                         
                          
                         
                           
                             
                               μ 
                               ~ 
                             
                             k 
                             
                               ( 
                               i 
                               ) 
                             
                           
                           · 
                           Δ 
                         
                          
                         
                             
                         
                          
                         
                           μ 
                           
                             1 
                             , 
                             k 
                           
                           
                             ( 
                             i 
                             ) 
                           
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     34 
                   
                   ] 
                 
               
             
             
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       μ 
                       
                         2 
                         , 
                         k 
                       
                     
                   
                   = 
                   
                     
                       min 
                       i 
                     
                      
                     
                       ( 
                       
                         Δ 
                          
                         
                           
                             
                               μ 
                               _ 
                             
                             k 
                             
                               ( 
                               i 
                               ) 
                             
                           
                           · 
                           Δ 
                         
                          
                         
                           
                             
                               μ 
                               ~ 
                             
                             k 
                             
                               ( 
                               i 
                               ) 
                             
                           
                           · 
                           Δ 
                         
                          
                         
                             
                         
                          
                         
                           μ 
                           
                             2 
                             , 
                             k 
                           
                           
                             ( 
                             i 
                             ) 
                           
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     35 
                   
                   ] 
                 
               
             
           
         
       
     
         [0000]    According to another embodiment also following a simple strategy, the autocorrelation-based decrements are treated separate from the other two types of decrements (gain-based and spectral colouring based). In this case, the Δ{tilde over (μ)} k   (i)  should not be correspond to the maximum gain but rather be appropriate for a typical gain: 
         [0000]    
       
         
           
             
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       μ 
                       
                         1 
                         , 
                         k 
                       
                     
                   
                   = 
                   
                     
                       min 
                       i 
                     
                      
                     
                       ( 
                       
                         min 
                          
                         
                           ( 
                           
                             
                               Δ 
                                
                               
                                 
                                   
                                     μ 
                                     _ 
                                   
                                   k 
                                   
                                     ( 
                                     i 
                                     ) 
                                   
                                 
                                 · 
                                 Δ 
                               
                                
                               
                                   
                               
                                
                               
                                 μ 
                                 
                                   1 
                                   , 
                                   k 
                                 
                                 
                                   ( 
                                   i 
                                   ) 
                                 
                               
                             
                             , 
                             
                               Δ 
                                
                               
                                 
                                   μ 
                                   ~ 
                                 
                                 k 
                                 
                                   ( 
                                   i 
                                   ) 
                                 
                               
                             
                           
                           ) 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     36 
                   
                   ] 
                 
               
             
             
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       μ 
                       
                         2 
                         , 
                         k 
                       
                     
                   
                   = 
                   
                     
                       min 
                       i 
                     
                      
                     
                       ( 
                       
                         min 
                          
                         
                           ( 
                           
                             
                               Δ 
                                
                               
                                 
                                   
                                     μ 
                                     _ 
                                   
                                   k 
                                   
                                     ( 
                                     i 
                                     ) 
                                   
                                 
                                 · 
                                 Δ 
                               
                                
                               
                                   
                               
                                
                               
                                 μ 
                                 
                                   2 
                                   , 
                                   k 
                                 
                                 
                                   ( 
                                   i 
                                   ) 
                                 
                               
                             
                             , 
                             
                               Δ 
                                
                               
                                 
                                   μ 
                                   ~ 
                                 
                                 k 
                                 
                                   ( 
                                   i 
                                   ) 
                                 
                               
                             
                           
                           ) 
                         
                       
                       ) 
                     
                   
                 
               
               
                 
                   [ 
                   
                     Eq 
                     . 
                     
                         
                     
                      
                     37 
                   
                   ] 
                 
               
             
           
         
       
     
         [0130]    According to particular embodiments, the calculated value of the step size parameter is overruled if either a large correlation is detected or a loud sound suddenly occurs. Under these circumstances, the adaptation of the cancelling filter coefficients is suspended. That is, 
         [0000]    
       
         
           
             
               
                 If 
                  
                 
                     
                 
                  
                 
                   
                     max 
                     i 
                   
                    
                   
                     ( 
                     
                       
                         max 
                         τ 
                       
                        
                       
                         ( 
                         
                           
                             r 
                             k 
                             
                               ( 
                               i 
                               ) 
                             
                           
                            
                           
                             ( 
                             τ 
                             ) 
                           
                         
                         ) 
                       
                     
                     ) 
                   
                 
               
               &gt; 
               
                 T 
                 max 
               
             
             , 
           
         
       
     
         [0000]    or if a sudden loud sound is detected according to the circuit shown in  FIG. 6 , Then μ 1,k =μ 2,k =0. 
         [0131]    In the following, measures according to embodiments of the pre-sent invention of how to adjust the adaptation rate of a feedback cancellation filter in a hearing aid in dependency of the acoustic environment of the hearing aid are summarised. 
         [0132]    When the amplification gain is increased (decreased) by a factor Δ compared to a nominal gain, the step size should be increased (decreased) by Δ 2  compared to the nominal step size. 
         [0133]    When operating with multiple frequency bands, the lowest amplification gain is decisive; if the lowest gain is increased (decreased) by a factor Δ compared to a nominal gain, the step size should be increased (decreased) by Δ 2  compared to the nominal step size. 
         [0134]    If the autocorrelation is high as measured by e.g., Eq. 2, Eq. 3, Eq. 4, or Eq. 5 the step size is increased substantially. 
         [0135]    A monotonic correspondence between the autocorrelation or a similar measure of a signals self-similarity and the step size is implemented such that the step size is reduced for increasing correlation or “self-similarity”. 
         [0136]    When the autocorrelation or similar measure of a signals self-similarity indicates that a pure tone is present in the signal, the adaptation is deactivated (step size=0). 
         [0137]    In a multi-band hearing aid, the autocorrelation or similar measure of a signals self-similarity can be calculated within each band. It is suggested to take the maximum of absolute values of the autocorrelation over bands and let this be decisive for the step size. 
         [0138]    If a sudden increase in sound pressure occurs in the incoming signal, the adaptation should be deactivated. This deactivation is maintained for a while after the incident. 
         [0139]    In a directional system working on wide-band signals, the efficiency of the system is defined by the ratio between the feedback compensated signal(s) and the directional output signal. If the norm is reduced by a factor Δ, the step size should be decreased by Δ 2  compared to the nominal step size. 
         [0140]    For a multi-band directional system the efficiency is calculated within in each band. The step size is reduced according to the largest factor Δ i   2  calculated over bands. 
         [0141]    In the multi-band case, combine amplification gain and efficiency of directional system for each band and then select step size as the maximum reduction of the nominal value. 
         [0142]    When operating with a multi-band system: combine “gain control”, “correlation control” and “directional filter control” in bands to find a set of equivalent step sizes. Next, take the minimum of these and use this as the resulting step size. 
         [0143]    According to further embodiments, these principles may well be applied to hearing aids with more than two microphones. 
         [0144]    All appropriate combinations of features described above are to be considered as belonging to the invention, even if they have not been explicitly described in their combination. 
         [0145]    According to embodiments of the present invention, hearing aids described herein may be implemented on signal processing devices suitable for the same, such as, e.g., digital signal processors, analogue/digital signal processing systems including field programmable gate arrays (FPGA), standard processors, or application specific signal processors (ASSP or ASIC). Obviously, it is preferred that the whole system is implemented in a single digital component even though some parts could be implemented in other ways—all known to the skilled person. 
         [0146]    Hearing aids, methods and devices according to embodiments of the present invention may be implemented in any suitable digital signal processing system. The hearing aids, methods and devices may also be used by, e.g., the audiologist in a fitting session. Methods according to the present invention may also be implemented in a computer program containing executable program code executing methods according to embodiments described herein. If a client-server-environment is used, an embodiment of the present invention comprises a remote server computer that embodies a system according to the present invention and hosts the computer program executing methods according to the present invention. According to another embodiment, a computer program product like a computer readable storage medium, for example, a floppy disk, a memory stick, a CD-ROM, a DVD, a flash memory, or any other suitable storage medium, is provided for storing the computer program according to the present invention. 
         [0147]    According to a further embodiment, the program code may be stored in a memory of a digital hearing device or a computer memory and executed by the hearing aid device itself or a processing unit like a CPU thereof or by any other suitable processor or a computer executing a method according to the described embodiments. 
         [0148]    Having described and illustrated the principles of the present invention in embodiments thereof, it should be apparent to those skilled in the art that the present invention may be modified in arrangement and detail without departing from such principles. Changes and modifications within the scope of the present invention may be made without departing from the spirit thereof, and the present invention includes all such changes and modifications.