# Patent

## Title
Active noise control system and method

## Abstract

An Active Noise Control (ANC) for controlling a noise produced by a noise source may include an acoustic sensor (212) to sense a noise pattern and to produce a noise signal corresponding to the sensed noise pattern, an estimator (202) to produce a predicted noise signal by applying an estimation function to the noise signal, and an acoustic transducer (216) to produce a noise destructive pattern based on the predicted noise signal.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
The present application is a US National Phase of PCT Application No. PCT/IL2004/000863, filed on Sep. 19, 2004, which claims the benefit under 35 U.S.C. 119(e) of US Provisional Application No. 60/503,471 filed Sep. 17, 2003 and is a continuation-in-part of U.S. application Ser. No. 09/120,973 filed Jul. 22, 1998 which claims benefit from Israeli Application 121555 filed Aug. 14, 1997, the disclosure of which is incorporated herein by reference.


FIELD OF THE INVENTION
The invention relates to the field of active noise control.
BACKGROUND
Conventional passive noise control systems may include “insulation” elements, silencers, vibration mounts, damping treatments, absorptive treatments, e.g., ceiling tiles, and/or conventional mufflers, e.g., mufflers as may be used in the automobile industry. The dimensions and/or mass of such passive noise control systems may usually depend on the acoustic pattern length of the noise intended to be reduced. Generally, passive noise control systems implemented to reduce noises of relatively low frequencies are bulky, large, heavy and/or expensive.
SUMMARY
According to embodiments of the invention, Active Noise Control (ANC) may be used to reduce noise energy and wave amplitude of a source noise pattern via an ANC sound system, which produces a noise-destructive pattern related to the source noise pattern such that a reduced noise zone may be created.
According to an exemplary embodiment of the invention, the ANC system may include an acoustic sensor, e.g., a microphone, to sense a noise pattern and to produce a noise signal corresponding to the sensed noise pattern; an estimator to produce a predicted noise signal by applying an estimation function to the noise signal; and an acoustic transducer, e.g., a speaker, to produce a noise destructive pattern based on the predicted noise signal.
According to some exemplary embodiments of the invention, the estimation function may include a non-linear estimation function, e.g., a radial basis function.
The estimator may be able to adapt one or more parameters of the estimation function based on a noise error at a predetermined location. For example, the ANC system may include an error evaluator to evaluate the noise error based on the noise signal and the predicted noise signal. Additionally or alternatively, the system may include an error sensing acoustic sensor to sense the noise error at the predetermined location.
The error evaluator may include a speaker transfer function module to produce an estimation of the noise destructive pattern, e.g., by applying a speaker transfer function to the predicted noise signal; a modulation transfer function module to produce an estimation of the noise pattern at the predetermined location, e.g., by applying a modulation transfer function to the noise signal; and a subtractor to subtract the estimation of the noise destructive pattern from the estimation of the noise pattern.
According to some exemplary embodiments, the estimator may be able to adapt the one or more parameters based on a predetermined criterion. For example, the estimator may be able to reduce, e.g., minimize, the error value by adapting the one or more parameters.
According to another exemplary embodiment of the invention, the ANC system may include a primary acoustic sensor, e.g., a microphone, to sense a noise pattern and to produce a corresponding primary noise signal; at least one secondary acoustic sensor, e.g., microphone, to sense a residual noise pattern and to produce at least one secondary noise signal corresponding to the residual noise pattern sensed by the at least one secondary microphone, respectively, wherein the at least one secondary acoustic sensor is separated from the noise source by a distance larger than a distance between the primary acoustic sensor and the noise source; and a controller to control an acoustic transducer to produce a noise destructive pattern based on the primary noise signal and the at least one secondary noise signal.
The controller may include, for example, a primary estimator to produce a predicted primary signal, e.g., by applying a primary estimation function to the primary noise signal; and at least one secondary estimator to produce at least one predicted secondary signal by applying at least one secondary estimation function to the at least one secondary noise signal, respectively.
The primary estimator may be able, for example, to iteratively adapt one or more parameters of the primary estimation function based on a noise error. The at least one secondary estimator may be able, for example, to iteratively adapt one or more parameters of the at least one secondary estimation function, respectively, based on the noise error.
The controller may control the acoustic transducer based on a combination of the predicted primary signal and the at least one predicted secondary signal.



BREIF DESCRIPTION OF THE DRAWINGS
The subject matter regarded as the invention is particularly pointed out and distinctly claimed in the concluding portion of the specification. The invention, however, both as to organization and method of operation, together with objects, features and advantages thereof, may best be understood by reference to the following detailed description when read with the accompanied drawings in which:
FIG. 1 is a schematic illustration of an active noise control system according to an exemplary embodiment of the invention;
FIG. 2 is a schematic illustration of a controller according to some exemplary embodiments of the invention that may be used, for example, in conjunction with the system of FIG. 1;
FIG. 3 is a schematic illustration of an active noise control system according to another exemplary embodiment of the invention; and
FIG. 4 is a schematic illustration of a controller according to other exemplary embodiments of the invention that may be used, for example, in conjunction with the system of FIG. 3.



It will be appreciated that for simplicity and clarity of illustration, elements shown in the drawings have not necessarily been drawn accurately or to scale. For example, the dimensions of some of the elements may be exaggerated relative to other elements for clarity or several physical components included in one functional block or element. Further, where considered appropriate, reference numerals may be repeated among the drawings to indicate corresponding or analogous elements. Moreover, some of the blocks depicted in the drawing may be combined into a single function.
DETAILED DESCRIPTION
In the following detailed description, numerous specific details are set forth in order to provide a thorough understanding of the invention. However, it will be understood by those of ordinary skill in the art that the present invention may be practiced without these specific details. In other instances, well-known methods, procedures, components and circuits may not have been described in detail so as not to obscure the present invention.
According to embodiments of the invention, Active Noise Control (ANC) may be used to reduce noise energy and wave amplitude of a source noise pattern, e.g., including one or more acoustic waves, via an ANC sound system, which produces a noise-destructive pattern, e.g., including one or more acoustic waves, related to the source noise pattern such that a reduced noise zone may be created.
Embodiments of the invention include ANC systems and methods, which may be efficiently implemented for reducing undesirable noises, e.g., at least noises of generally low frequencies, as described below.
Certain aspects of ANC methods and systems, in accordance with some exemplary embodiments of the invention, are described in U.S. patent application Ser. No. 09/120,973, filed Jul. 22, 1998, entitled “ACTIVE ACOUSTIC NOISE REDUCTION SYSTEM”; and in European Patent Application 02023483.7, filed Oct. 21, 2002, entitled “ACTIVE ACOUSTIC NOISE REDUCTION SYSTEM”, and published Apr. 28, 2004 as publication number 1414021. The entire disclosure of both of these applications is incorporated herein by reference.
Reference is made to FIG. 1, which schematically illustrates an ANC system 100 according to an exemplary embodiment of the invention.
ANC system 100 may include, for example, a acoustic sensor, e.g., a microphone 102, denoted MIC1, to sense the noise energy and/or wave amplitude of a noise pattern produced by a noise source 104. Microphone 102 may include any suitable microphone able to generate an output noise signal 103, corresponding to the noise pattern sensed by microphone 112. For example, microphone 102 may include microphone Part No. ECM6AP, available from ARIO Electronics Co. Ltd., Taoyuan, Taiwan. Noise signal 103 may include, for example, a sequence of N samples per second. For example, N may be 1000 samples per second, e.g., if microphone 103 operates at a sampling rate of about 10 KHz.
ANC system 100 may also include an acoustic transducer, e.g., a speaker 108, and a controller 106 to control speaker 108 to produce a noise destructive pattern to reduce or cancel the noise energy and/or wave amplitude of the noise pattern, e.g., within a reduced-noise zone 110, as described in detail below. Speaker 108 may include any suitable speaker, e.g., as is known in the art. For example, speaker 108 may include speaker Part No. AI 4.0, available from Cerwin-Vega Inc., Chatsworth, Calif.
According to some exemplary embodiments of the invention, controller 106 may be able to evaluate a noise error corresponding to an anticipated destructive interference between the noise pattern and the noise destructive pattern at a predetermined location 112 within zone 110, as described below. The noise error may be evaluated, for example, by controller 106, e.g., based on noise signal 103, as described below. Additionally or alternatively, the noise error may be sensed by an error-sampling microphone positioned at the predetermined location, as described below. Controller 106 may control speaker 108 to produce the noise destructive pattern, e.g., based on noise signal 103 and/or on the evaluated noise error, as described below.
According to some exemplary embodiments of the invention, it may be desired to control the timing at which the noise destructive pattern is produced, e.g., in order to efficiently control, e.g., reduce, the noise within zone 110. For example, it may be desired to controllably time the noise destructive pattern corresponding to a sample of the noise pattern such that the destructive noise pattern reaches a location within zone 110, e.g., location 112, at substantially the same time the sampled noise pattern reaches the same location.
According to embodiments of the invention, there may be a time delay between the time at which a currently sampled noise pattern reaches location 112 and the time at which the noise destructive pattern corresponding to the current sample of the noise pattern reaches location 112. This time delay may result, for example, from the time required for microphone 102 to sense the noise pattern, the time required for controller to process noise signal 103, and/or the time required for speaker 108 to produce the noise destructive pattern.
Thus, according to some exemplary embodiments of the invention, controller 106 may estimate a sample of the noise pattern succeeding the current sample (“the succeeding sample”) based on the current sample and/or one or more previous samples of the noise pattern. Controller 118 may provide an input to speaker 113, such that speaker 113 produces the noise destructive pattern based on the estimated succeeding sample, e.g., such that the noise destructive pattern may reach location 112 substantially at the same time the noise pattern reaches location 112.
An acoustic pattern, e.g., the noise pattern, may be characterized by a generally non-linear function. Thus, according to exemplary embodiments of the invention, controller 106 may use non-linear estimation to estimate the succeeding sample. Such non-linear estimation may provide, according to exemplary embodiments of the invention, a better estimation of the succeeding sample compared to a corresponding linear estimation. However, according to other embodiments of the invention, controller 106 may use any other suitable estimation, e.g., a linear estimation, to estimate the succeeding sample.
According to exemplary embodiments of the invention, controller 106 may include an estimator 121 to produce a predicted noise signal 114 by applying an estimation function to one or more samples of noise signal 103. Speaker 113 may produce the noise destructive pattern based on predicted noise signal 114, as described below.
Reference is made to FIG. 2, which schematically illustrates a controller 200 according to some exemplary embodiments of the invention. Although the invention is not limited in this respect, controller 200 may be implemented by ANC system 100 (FIG. 1).
According to exemplary embodiments of the invention, controller 200 may include an estimator 202 to receive from an acoustic sensor, e.g., a microphone 212, a noise signal 210, e.g., including a plurality of samples of a sensed noise pattern. Estimator 202 may generate a predicted noise signal 230 having a value, y(n), corresponding to an n-th sample, denoted MIC(n), received from microphone 212, by applying an estimation function F to the sample MIC(n) and to one or more other samples previously received from microphone 212, as described below. Controller 202 may control an acoustic transducer, e.g., a speaker 216, to generate a noise destructive pattern 218, e.g., based on output 230.
According to some exemplary embodiments of the invention, estimator 202 may implement a non-linear estimation algorithm, as described below.
According to some exemplary embodiments of the invention, estimator 202 may implement a Radial Basis Function (RBF) algorithm, as described below.
Estimator 202 may implement the RBF algorithm to estimate the value of a succeeding sample of the noise signal based on the values of one or more samples of the noise signal received from microphone 212. For example, the RBF algorithm may correspond to a combination of a set of K radial n-dimension functions, wherein each function may differ in one or more parameters, e.g., a center of the function parameter, denoted ck, an effective radius parameter, denoted vk, and/or and intensity of the function, denoted wk, as are known in the art. For example, estimator 202 may implement a RBF algorithm analogous to the one described by S. Haykin, “Adaptive Filter Theory”, 3rd edition, Prentice Hall, pp. 863-565.
According to some exemplary embodiments of the invention, estimator 202 may generate predicted noise 230 according to the following equation:




  
    
      
        
          y
          ⁡
          
            [
            n
            ]
          
        
        =
        
          
            ∑
            
              k
              =
              1
            
            K
          
          ⁢
          
            
              w
              k
            
            ⁢
            
              exp
              ⁡
              
                (
                
                  
                    -
                    
                      1
                      
                        2
                        ⁢
                        
                          v
                          k
                        
                      
                    
                  
                  ⁢
                  
                    
                      ∑
                      
                        i
                        =
                        0
                      
                      
                        L
                        -
                        1
                      
                    
                    ⁢
                    
                      
                        (
                        
                          
                            MIC
                            ⁡
                            
                              [
                              
                                n
                                -
                                i
                              
                              ]
                            
                          
                          -
                          
                            
                              c
                              k
                            
                            ⁡
                            
                              [
                              i
                              ]
                            
                          
                        
                        )
                      
                      2
                    
                  
                
                )
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




wherein L denotes a determined number of samples of the noise signal to be implemented for the estimation of y(n).
According to some exemplary embodiments of the invention, estimator 202 may iteratively adapt one or more parameters, e.g., one or more of the parameters ck, vk, and wk, of the estimation function F, e.g., based on a predetermined criterion, as described below.
According to some exemplary embodiments of the invention, estimator 202 may iteratively adapt one or more of the parameters ck, vk, and wk based on the evaluated noise error at a predetermined location, e.g., location 112 (FIG. 1), as described below.
According to some exemplary embodiments of the invention, controller 200 may also include an error evaluation module 203 to evaluate the noise error, e.g., based on noise signal 210 and predicted noise signal 230, as described below.
According to some exemplary embodiments of the invention, module 203 may include, for example, a Modulation Transfer Function (MTF) module 204 to apply to noise signal 210 a predetermined MTF, thereby to generate an output 241 having a value corresponding to an estimation, denoted d(n), of the n-th sample of the noise pattern at the predetermined location. The MTF may be determined, for example, based on characteristics of microphone 212 and/or based on geometrical and/or physical characteristics of a path and/or a medium, e.g., air, between microphone 212 and the predetermined location, e.g., as known in the art. MTF module 204 may include any suitable hardware and/or software, e.g., as are known in the art, to apply a predetermined MTF to noise signal 210.
According to exemplary embodiments of the invention, module 203 may also include a Speaker Transfer Function (STF) module 206 to apply a STF to predicted noise signal 230, thereby to generate an output 249 having a value corresponding to an estimation of noise destructive pattern 218 produced in response to predicted noise signal 230. The STF may be determined, for example, based on characteristics of speaker 216, e.g., as known in the art. STF module 206 may include any suitable hardware and/or software, e.g., as are known in the art, to apply a predetermined STF to predicted noise signal 230. For example, the value, denoted z(n), of output 249 may be calculated using the following equation:




  
    
      
        
          z
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            ∑
            
              s
              =
              0
            
            
              S
              -
              1
            
          
          ⁢
          
            
              STF
              ⁡
              
                (
                s
                )
              
            
            ⁢
            
              y
              ⁡
              
                (
                
                  n
                  -
                  s
                
                )
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




wherein S denotes a predetermined STF frequency parameter vector, as is known in the art.
Substituting Equation 1 in Equation 2 may yield the following equation:




  
    
      
        
          z
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            ∑
            
              s
              =
              0
            
            
              S
              -
              1
            
          
          ⁢
          
            
              STF
              ⁡
              
                (
                s
                )
              
            
            ⁢
            
              
                ∑
                
                  k
                  =
                  1
                
                K
              
              ⁢
              
                
                  w
                  k
                
                ⁢
                
                  exp
                  ⁡
                  
                    (
                    
                      
                        -
                        
                          1
                          
                            2
                            ⁢
                            
                              
                            
                            ⁢
                            
                              υ
                              k
                            
                          
                        
                      
                      ⁢
                      
                        
                          ∑
                          
                            i
                            =
                            0
                          
                          
                            L
                            -
                            1
                          
                        
                        ⁢
                        
                          
                            (
                            
                              
                                x
                                ⁡
                                
                                  (
                                  
                                    n
                                    -
                                    s
                                    -
                                    i
                                  
                                  )
                                
                              
                              -
                              
                                
                                  c
                                  k
                                
                                ⁡
                                
                                  (
                                  i
                                  )
                                
                              
                            
                            )
                          
                          2
                        
                      
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




According to exemplary embodiments of the invention, module 203 may also include a subtractor 208, which may be implemented by any suitable hardware and/or software as are known in the art. Subtractor 208 may subtract the value of the estimated noise destructive pattern, e.g., of output STF 249, from the estimated value of the noise pattern, e.g., of output 241, to produce an output 245 including the evaluated noise error, denoted e(n), corresponding to sample MIC(n).
According to exemplary embodiments of the invention, estimator 202 may implement an adaptive algorithm to iteratively adapt the values of one or more of the parameters vk, ck, and wk, e.g., based on the value of the noise error, as described below.
According to exemplary embodiments of the invention, the value of the noise error e(n), corresponding to the n-th sample of noise signal 210 may be estimated using the following equation:

e(n)=d(n)−z(n)  (4)

Substituting Equation 3 in Equation 4 may yield the following equation:




  
    
      
        
          e
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            d
            ⁡
            
              (
              n
              )
            
          
          -
          
            
              ∑
              
                s
                =
                0
              
              
                S
                -
                1
              
            
            ⁢
            
              
                STF
                ⁡
                
                  (
                  s
                  )
                
              
              ⁢
              
                
                  ∑
                  
                    k
                    =
                    1
                  
                  K
                
                ⁢
                
                  
                    w
                    k
                  
                  ⁢
                  
                    exp
                    ⁡
                    
                      (
                      
                        
                          -
                          
                            1
                            
                              2
                              ⁢
                              
                                
                              
                              ⁢
                              
                                υ
                                k
                              
                            
                          
                        
                        ⁢
                        
                          
                            ∑
                            
                              i
                              =
                              0
                            
                            
                              L
                              -
                              1
                            
                          
                          ⁢
                          
                            
                              (
                              
                                
                                  x
                                  ⁡
                                  
                                    (
                                    
                                      n
                                      -
                                      s
                                      -
                                      i
                                    
                                    )
                                  
                                
                                -
                                
                                  
                                    c
                                    k
                                  
                                  ⁡
                                  
                                    (
                                    i
                                    )
                                  
                                
                              
                              )
                            
                            2
                          
                        
                      
                      )
                    
                  
                
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




According to some exemplary embodiments of the invention, estimator 202 may iteratively adapt one or more of the parameters vk, ck, and wk, to reduce, e.g., minimize, the arithmetic mean, denoted E[(e(n))2], of the square of the noise error. For example, estimator 202 may be able to iteratively adapt one or more of the parameters of the estimation function such that the partial derivative of E[(e(n))2] with respect to one or more of the parameters, respectively, is equal to zero, as described below.
According to some exemplary embodiments of the invention, the arithmetic mean of the square of the estimated noise error may be calculated using the following equation:




  
    
      
        
          E
          ⁡
          
            [
            
              
                (
                
                  e
                  ⁡
                  
                    (
                    n
                    )
                  
                
                )
              
              2
            
            ]
          
        
        =
        
          E
          ⁡
          
            [
            
              
                (
                
                  
                    d
                    ⁡
                    
                      (
                      n
                      )
                    
                  
                  -
                  
                    
                      ∑
                      
                        s
                        =
                        0
                      
                      
                        S
                        -
                        1
                      
                    
                    ⁢
                    
                      
                        STF
                        ⁡
                        
                          (
                          s
                          )
                        
                      
                      ⁢
                      
                        
                          ∑
                          
                            k
                            =
                            1
                          
                          K
                        
                        ⁢
                        
                          
                            w
                            k
                          
                          ⁢
                          
                            
                              f
                              k
                            
                            ⁡
                            
                              [
                              
                                n
                                -
                                s
                              
                              ]
                            
                          
                        
                      
                    
                  
                
                )
              
              2
            
            ]
          
        
      
    
    
      
        (
        6
        )
      
    
  




wherein





  
    
      
        
          
            f
            k
          
          ⁡
          
            [
            
              n
              -
              s
            
            ]
          
        
        =
        
          exp
          (
          
            
              -
              
                1
                
                  2
                  ⁢
                  
                    
                  
                  ⁢
                  
                    υ
                    k
                  
                
              
            
            ⁢
            
              
                ∑
                
                  i
                  =
                  0
                
                
                  L
                  -
                  1
                
              
              ⁢
              
                
                  (
                  
                    
                      x
                      ⁡
                      
                        (
                        
                          n
                          -
                          s
                          -
                          i
                        
                        )
                      
                    
                    -
                    
                      
                        c
                        k
                      
                      ⁡
                      
                        (
                        i
                        )
                      
                    
                  
                  )
                
                2
              
            
          
          )
        
      
    
    
      
        (
        7
        )
      
    
  




The partial derivatives of Equation 6 with respect to the parameters ck, vk, and wk, respectively, may be calculated using the following equations:




  
    
      
        
          
            ∂
            
              E
              ⁡
              
                [
                
                  
                    (
                    
                      e
                      ⁡
                      
                        (
                        n
                        )
                      
                    
                    )
                  
                  2
                
                ]
              
            
          
          
            ∂
            
              w
              k
            
          
        
        =
        
          E
          ⁡
          
            [
            
              
                -
                2
              
              ⁢
              
                e
                ⁡
                
                  (
                  n
                  )
                
              
              ⁢
              
                
                  ∑
                  
                    s
                    =
                    0
                  
                  
                    S
                    -
                    1
                  
                
                ⁢
                
                  
                    STF
                    ⁡
                    
                      (
                      s
                      )
                    
                  
                  ⁢
                  
                    
                      f
                      k
                    
                    ⁡
                    
                      [
                      
                        n
                        -
                        s
                      
                      ]
                    
                  
                
              
            
            ]
          
        
      
    
    
      
        (
        8
        )
      
    
  
  
    
      
        
          
            ∂
            
              E
              ⁡
              
                [
                
                  
                    (
                    
                      e
                      ⁡
                      
                        (
                        n
                        )
                      
                    
                    )
                  
                  2
                
                ]
              
            
          
          
            ∂
            
              c
              k
            
          
        
        =
        
          -
          
            E
            ⁡
            
              [
              
                2
                ⁢
                
                  e
                  ⁡
                  
                    (
                    n
                    )
                  
                
                ⁢
                
                  w
                  k
                
                ⁢
                
                  
                    ∑
                    
                      s
                      =
                      0
                    
                    
                      S
                      -
                      1
                    
                  
                  ⁢
                  
                    STF
                    ⁢
                    
                      (
                      s
                      )
                    
                    ⁢
                    
                      
                        f
                        k
                      
                      ⁡
                      
                        [
                        
                          n
                          -
                          s
                        
                        ]
                      
                    
                    ⁢
                    
                      (
                      
                        
                          1
                          
                            υ
                            k
                          
                        
                        ⁢
                        
                          
                            ∑
                            
                              i
                              =
                              0
                            
                            
                              L
                              -
                              1
                            
                          
                          ⁢
                          
                            (
                            
                              
                                x
                                ⁡
                                
                                  (
                                  
                                    n
                                    -
                                    i
                                  
                                  )
                                
                              
                              -
                              
                                
                                  c
                                  k
                                
                                ⁡
                                
                                  (
                                  i
                                  )
                                
                              
                            
                            )
                          
                        
                      
                      )
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        9
        )
      
    
  
  
    
      
        
          
            ∂
            
              E
              ⁡
              
                [
                
                  
                    (
                    
                      e
                      ⁡
                      
                        (
                        n
                        )
                      
                    
                    )
                  
                  2
                
                ]
              
            
          
          
            ∂
            
              υ
              k
            
          
        
        =
        
          E
          ⁡
          
            [
            
              
                e
                ⁡
                
                  (
                  n
                  )
                
              
              ⁢
              
                w
                k
              
              ⁢
              
                
                  ∑
                  
                    s
                    =
                    0
                  
                  
                    S
                    -
                    1
                  
                
                ⁢
                
                  
                    STF
                    ⁡
                    
                      (
                      s
                      )
                    
                  
                  ⁢
                  
                    
                      f
                      k
                    
                    ⁡
                    
                      [
                      
                        n
                        -
                        s
                      
                      ]
                    
                  
                  ⁢
                  
                    1
                    
                      
                        (
                        
                          υ
                          k
                        
                        )
                      
                      2
                    
                  
                  ⁢
                  
                    
                      ∑
                      
                        i
                        =
                        0
                      
                      
                        L
                        -
                        1
                      
                    
                    ⁢
                    
                      
                        (
                        
                          
                            x
                            ⁡
                            
                              (
                              
                                n
                                -
                                i
                              
                              )
                            
                          
                          -
                          
                            
                              c
                              k
                            
                            ⁡
                            
                              (
                              i
                              )
                            
                          
                        
                        )
                      
                      2
                    
                  
                
              
            
            ]
          
        
      
    
    
      
        (
        10
        )
      
    
  




A minimum value of E[(e(n))2]may be determined by from the following equations:




  
    
      
        
          
            ∂
            
              E
              ⁡
              
                [
                
                  
                    (
                    
                      e
                      ⁡
                      
                        (
                        n
                        )
                      
                    
                    )
                  
                  2
                
                ]
              
            
          
          
            ∂
            
              w
              k
            
          
        
        =
        0
      
    
    
      
        (
        11
        )
      
    
  
  
    
      
        
          
            ∂
            
              E
              ⁡
              
                [
                
                  
                    (
                    
                      e
                      ⁡
                      
                        (
                        n
                        )
                      
                    
                    )
                  
                  2
                
                ]
              
            
          
          
            ∂
            
              c
              k
            
          
        
        =
        0
      
    
    
      
        (
        12
        )
      
    
  
  
    
      
        
          
            ∂
            
              E
              ⁡
              
                [
                
                  
                    (
                    
                      e
                      ⁡
                      
                        (
                        n
                        )
                      
                    
                    )
                  
                  2
                
                ]
              
            
          
          
            ∂
            
              υ
              k
            
          
        
        =
        0
      
    
    
      
        (
        13
        )
      
    
  




Applying the condition of Equation 11 to Equation 8 may result in the following relation between an adapted value, denoted wk(n+1), and the current value, wk(n), of the parameter wk:




  
    
      
        
          
            w
            k
          
          ⁡
          
            (
            
              n
              +
              1
            
            )
          
        
        =
        
          
            
              w
              k
            
            ⁡
            
              (
              n
              )
            
          
          -
          
            
              μ
              w
            
            ⁢
            
              e
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              
                ∑
                
                  s
                  =
                  0
                
                
                  S
                  -
                  1
                
              
              ⁢
              
                
                  STF
                  ⁡
                  
                    (
                    s
                    )
                  
                
                ⁢
                
                  
                    f
                    k
                  
                  ⁡
                  
                    [
                    
                      n
                      -
                      s
                    
                    ]
                  
                
              
            
          
        
      
    
    
      
        (
        14
        )
      
    
  




wherein μk is a determined convergence parameter corresponding to wk.

Applying the condition of Equation 12 to Equation 9 may result in the following relation between an adapted value, denoted ck(n+1), and the current value, ck(n), of the parameter ck:




  
    
      
        
          
            c
            k
          
          ⁡
          
            (
            
              n
              +
              1
            
            )
          
        
        =
        
          
            
              c
              k
            
            ⁡
            
              (
              n
              )
            
          
          -
          
            
              μ
              c
            
            ⁢
            
              e
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              w
              k
            
            ⁢
            
              
                ∑
                
                  s
                  =
                  0
                
                
                  S
                  -
                  1
                
              
              ⁢
              
                
                  STF
                  ⁡
                  
                    (
                    s
                    )
                  
                
                ⁢
                
                  
                    f
                    k
                  
                  ⁡
                  
                    [
                    
                      n
                      -
                      s
                    
                    ]
                  
                
                ⁢
                
                  (
                  
                    
                      1
                      
                        υ
                        k
                      
                    
                    ⁢
                    
                      
                        ∑
                        
                          i
                          =
                          0
                        
                        
                          L
                          -
                          1
                        
                      
                      ⁢
                      
                        (
                        
                          
                            x
                            ⁡
                            
                              (
                              
                                n
                                -
                                i
                              
                              )
                            
                          
                          -
                          
                            
                              c
                              k
                            
                            ⁡
                            
                              (
                              i
                              )
                            
                          
                        
                        )
                      
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




wherein μc is a determined convergence parameter corresponding to ck.

Applying the condition of Equation 13 to Equation 10 may result in the following relation between an adapted value, denoted vk(n+1), and the current value, vk(n), of the parameter vk:




  
    
      
        
          
            υ
            k
          
          ⁡
          
            (
            
              n
              +
              1
            
            )
          
        
        =
        
          
            
              υ
              k
            
            ⁡
            
              (
              n
              )
            
          
          -
          
            
              μ
              υ
            
            ⁢
            
              e
              ⁡
              
                (
                n
                )
              
            
            ⁢
            
              w
              k
            
            ⁢
            
              
                ∑
                
                  s
                  =
                  0
                
                
                  S
                  -
                  1
                
              
              ⁢
              
                
                  STF
                  ⁡
                  
                    (
                    s
                    )
                  
                
                ⁢
                
                  
                    f
                    k
                  
                  ⁡
                  
                    [
                    
                      n
                      -
                      s
                    
                    ]
                  
                
                ⁢
                
                  1
                  
                    
                      (
                      
                        υ
                        k
                      
                      )
                    
                    2
                  
                
                ⁢
                
                  
                    ∑
                    
                      i
                      =
                      0
                    
                    
                      L
                      -
                      1
                    
                  
                  ⁢
                  
                    
                      (
                      
                        
                          x
                          ⁡
                          
                            (
                            
                              n
                              -
                              i
                            
                            )
                          
                        
                        -
                        
                          
                            c
                            k
                          
                          ⁡
                          
                            (
                            i
                            )
                          
                        
                      
                      )
                    
                    2
                  
                
              
            
          
        
      
    
    
      
        (
        16
        )
      
    
  




wherein μv is a determined convergence parameter corresponding to vk.

According to some exemplary embodiments of the invention, adaptive estimator 202 may implement one or more of Equations 14-16 to iteratively adapt one or more of the parameters wk, ck, and vk, respectively.
Some exemplary embodiments of the invention relate to an ANC system, e.g., system 100 (FIG. 1), implementing an error evaluation module, e.g., module 203, to evaluate the noise error at a predetermined location, e.g., location 112 (FIG. 1). However, it will be appreciated by those skilled in the art, that according to other embodiments of the invention, any other one or more suitable modules may be implemented to evaluate the noise error. For example, an error-sensing microphone 239 may be located at the predetermined location, and an output 240 of error-sensing microphone 239 corresponding to the sensed noise error at the predetermined location may be provided to estimator 202.
Some exemplary embodiments of the invention relate to an ANC system, e.g., ANC system 100 (FIG. 1), including a controller, e.g., controller 106 (FIG. 1), to control an acoustic transducer, e.g., speaker 108 (FIG. 1), based on a noise signal of a noise pattern received from an acoustic sensor, e.g., microphone 102 (FIG. 1). However, other embodiments of the invention may refer to an ANC system including a controller able to control an acoustic transducer based on one or more noise signals of a noise pattern received from more than one acoustic sensor, e.g., as described below.
Reference is made to FIG. 3, which schematically illustrates an ANC system 300 according to another exemplary embodiment of the invention.
ANC system 300 may include, for example, a primary acoustic sensor, e.g., a microphone 302, denoted MIC1, to sample the noise energy and/or wave amplitude of a noise pattern produced by a noise source 304. Microphone 302 may include any suitable microphone, e.g., as described above with reference to microphone 102 (FIG. 1).
ANC system 300 may also include an acoustic transducer, e.g., a speaker 308, and a controller 306 able to control speaker 308 to produce a noise destructive pattern to reduce or cancel the noise energy and/or wave amplitude of the noise pattern, e.g., within a reduced-noise zone 310, as described in detail below. Speaker 308 may include any suitable speaker, e.g., as described above with reference to speaker 108 (FIG. 1).
According to some exemplary embodiments of the invention, controller 306 may be able to evaluate a noise error corresponding to a combination of, e.g., a difference between, the noise pattern and the noise destructive pattern, e.g., at a predetermined location 312 within zone 310, as described below. Controller 306 may control speaker 308 to produce the noise destructive pattern, for example, such that the noise error is reduced, e.g., minimized, as described below.
According to exemplary embodiments of the invention, a relatively good coherence between primary microphone 302 and the evaluation of the noise error, e.g., at the relevant frequencies of the noise pattern, may be required in order for ANC 300 to achieve an efficient degree of noise reduction, as described below. For example, the higher correlation between the noise pattern sampled by microphone 302 and the noise error, the higher the level of noise control, e.g., noise reduction, which may be achieved by ANC system 300. The coherence between the noise pattern sampled by microphone 302 and the noise error may depend, for example, on the geometric structure of the path between microphone 302 and location 312. Additionally or alternatively, the coherence between the noise pattern received by microphone 302 and the noise error may depend, for example, on the aerodynamic attributes, e.g., surface roughness, of the path. For example, no “eye contact” between microphone 302 and location 312 and/or a path having relatively rough surfaces may result in a reduced coherence between the signal received by microphone 302 and the evaluated noise error. Furthermore, the operation of ANC 300 to reduce the noise may be disturbed by formation of acoustic signals along the path between the microphone 302 and location 312, e.g., due to turbulent airflow and/or friction between the air and path materials, for example, if a structure of a device implementing one or more elements of ANC 300 does not have an aerodynamically optimized design, e.g., due to price and size constraints. Turbulent airflow may be characterized by stochastic formation of eddies which produce significant rustles, and friction between the air and the relatively rough surfaces may be characterized by conversion of kinetic energy into heat and noise energy.
According to exemplary embodiments of the invention, the noise error may be evaluated using a MTF, e.g., as described below with reference to FIG. 4. The MTF may be predetermined, e.g., based on one or more characteristics of the path between microphone 302 and location 312, and/or one or more expected characteristics of the noise-pattern. However, one or more of the characteristics of the path and/or the expected characteristics of the noise pattern may be different than the expected characteristics. As a result, the correlation between the noise error, e.g., evaluated based on the predetermined MTF, and the actual noise at location 312 may not be sufficiently accurate.
According to some exemplary embodiments of the invention, ANC system 300 may also include at least one secondary acoustic sensor, e.g., at least one secondary microphone 392, denoted MIC21, to sample the noise energy and/or wave amplitude of the noise pattern produced by noise source 304. Secondary microphone 392 may be separated from noise source 304 by a distance, d1, bigger than the distance, d2, between primary microphone 302 and noise source 304. For example, microphone 392 may be located along the path between microphone 302 and location 312. The distance d1-d2 between microphone 392 and microphone 302 may be large enough to allow microphone 392 to sample a residual noise pattern, e.g., a noise pattern formed by the path, which may not be received by microphone 302. Microphone 392 may include any suitable microphone, e.g., as described above with reference to microphone 102 (FIG. 1).
According to some exemplary embodiments of the invention, controller 306 may control speaker 308 to produce the noise destructive pattern based on the noise pattern sensed by microphone 302 and/or the residual noise pattern sensed by microphone 392, as described below.
Reference is made to FIG. 4, which schematically illustrates a controller 400 according to another exemplary embodiment of the invention. Although the invention is not limited in this respect, controller 400 may be implemented by ANC system 300 (FIG. 3).
According to exemplary embodiments of the invention, controller 400 may include a reference estimator 408 to receive from a primary microphone 402 a primary noise signal 412, e.g., including a plurality of samples. Estimator 408 may generate a predicted primary signal 414 having a value, y1(n), corresponding to an n-th sample, denoted MIC1(n), received from microphone 402, by applying a primary estimation function F1 to the sample MIC1(n) and to one or more other samples previously received from microphone 402, as described below.
According to exemplary embodiments of the invention, controller 400 may also include at least one secondary estimator 410 to receive from at least one secondary microphone 404 at least one secondary noise signal, respectively, e.g., including a plurality of samples. Estimator 410 may generate a predicted secondary signal 422 having a value, y2(n), corresponding to an n-th sample, denoted MIC21(n), received from microphone 404, by applying a secondary estimation function F2 to the sample MIC21(n) and to one or more other samples previously received from microphone 404, as described below.
Controller 400 may control an acoustic transducer, e.g., a speaker 406, to generate a noise destructive pattern 418, e.g., based on a combination of signal 414 and signal 422. For example, controller 400 may also include an adder 424, e.g., as is known in the art, to provide speaker 406 with an input 426 corresponding to the sum of signals 422 and 414.
According to some exemplary embodiments of the invention, estimator 408 may generate signal 414 according to the following equation:




  
    
      
        
          
            y
            1
          
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            ∑
            
              s
              =
              0
            
            
              L
              1
            
          
          ⁢
          
            
              
                W
                1
              
              ⁡
              
                (
                s
                )
              
            
            ⁢
            MIC
            ⁢
            
              
            
            ⁢
            1
            ⁢
            
              (
              
                n
                -
                s
              
              )
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




wherein W1 denotes a predetermined prediction filter (PF) vector of length L1 corresponding to estimation function F1.
According to some exemplary embodiments of the invention, estimator 410 may generate signal 422 according to the following equation:




  
    
      
        
          
            y
            2
          
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            ∑
            
              s
              =
              0
            
            
              L
              2
            
          
          ⁢
          
            
              
                W
                2
              
              ⁡
              
                (
                s
                )
              
            
            ⁢
            MIC
            ⁢
            
              
            
            ⁢
            21
            ⁢
            
              (
              
                n
                -
                s
              
              )
            
          
        
      
    
    
      
        (
        18
        )
      
    
  




wherein W2 denotes a predetermined PF vector of length L2 corresponding to estimation function F2.
According to some exemplary embodiments of the invention, estimator 408 may iteratively adapt the vector W1, and/or estimator 410 may iteratively adapt the vector W2, e.g., based on a predetermined criterion, as described below.
According to some exemplary embodiments of the invention, estimator 408 may iteratively adapt vector W1, based on the noise error corresponding to the combination of, e.g., the difference between, the noise pattern at the predetermined location, e.g., location 312 (FIG. 3), and an estimation of the contribution of signal y1(n) to noise destructive pattern 418, as described below.
According to some exemplary embodiments of the invention, controller 400 may also include a first evaluation module 430 to evaluate the noise error, e.g., based on signal 412 and signal 414, as described below.
According to some exemplary embodiments of the invention, module 430 may include, for example, a combiner 434 to combine signals 412 and 420. For example, combiner 434 may include a first MTF module 436 to apply a first predetermined MTF, denoted MTF1, to signal 412 and to divide the result by two. Combiner 434 may also include a second MTF module 438 to apply a second predetermined MTF, denoted MTF2, to signal 420 and to divide the result by two. For example, MTF1, may be determined, e.g., as known in the art, based on characteristics of microphone 402 and/or based on geometrical and/or physical characteristics of a path between microphone 412 and the certain location. MTF2, may be determined, for example, based on characteristics of microphone 404 and/or based on geometrical and/or physical characteristics of a path between microphone 404 and the predetermined location. Combiner 434 may also include an adder 440 to generate an output 442, denoted d(n), corresponding to an average between an estimation the n-th sample of the noise pattern at the certain location using MTF1, and an estimation the n-th sample of the noise pattern at the certain location using MTF2.
For example, d(n) may be calculated using the following equation:




  
    
      
        
          d
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            1
            2
          
          ⁢
          
            (
            
              
                
                  ∑
                  
                    x
                    =
                    0
                  
                  
                    M
                    1
                  
                
                ⁢
                
                  (
                  
                    
                      
                        MTF
                        1
                      
                      ⁡
                      
                        (
                        s
                        )
                      
                    
                    ⁢
                    Mic
                    ⁢
                    
                      
                    
                    ⁢
                    1
                    ⁢
                    
                      (
                      
                        n
                        -
                        s
                      
                      )
                    
                  
                  )
                
              
              +
              
                
                  ∑
                  
                    s
                    =
                    0
                  
                  
                    M
                    2
                  
                
                ⁢
                
                  (
                  
                    
                      
                        MTF
                        2
                      
                      ⁡
                      
                        (
                        s
                        )
                      
                    
                    ⁢
                    Mic
                    ⁢
                    
                      
                    
                    ⁢
                    21
                    ⁢
                    
                      (
                      
                        n
                        -
                        s
                      
                      )
                    
                  
                  )
                
              
            
            )
          
        
      
    
    
      
        (
        19
        )
      
    
  




wherein M1 denotes a predetermined number of samples of MTF1, and M2 denotes a predetermined number of samples of MTF2.

According to exemplary embodiments of the invention, module 430 may also include a STF module 450 to apply a STF to signal 414 to generate an output 452 representing an estimation of a primary part of the noise destructive pattern corresponding to predicted primary signal 414. The STF may be determined, for example, based on characteristics of speaker 406, e.g., as known in the art. STF module 450 may include any suitable hardware and/or software, e.g., as known in the art, to apply a predetermined STF to signal 414. For example, the value, denoted z1(n), of output 452 may be calculated using the following equation:




  
    
      
        
          
            z
            1
          
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            ∑
            
              p
              =
              0
            
            
              S
              -
              1
            
          
          ⁢
          
            
              STF
              ⁡
              
                (
                p
                )
              
            
            ⁢
            
              
                y
                1
              
              ⁡
              
                (
                
                  n
                  -
                  p
                
                )
              
            
          
        
      
    
    
      
        (
        20
        )
      
    
  




Substituting Equation 17 in Equation 20 may yield the following equation:




  
    
      
        
          
            z
            1
          
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            ∑
            
              p
              =
              0
            
            
              S
              -
              1
            
          
          ⁢
          
            
              STF
              ⁡
              
                (
                p
                )
              
            
            ⁢
            
              
                ∑
                
                  s
                  =
                  0
                
                
                  L
                  1
                
              
              ⁢
              
                
                  
                    W
                    1
                  
                  ⁡
                  
                    (
                    s
                    )
                  
                
                ⁢
                MIC
                ⁢
                
                  
                
                ⁢
                1
                ⁢
                
                  (
                  
                    n
                    -
                    s
                    -
                    p
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




According to exemplary embodiments of the invention, module 430 may also include a subtractor 454, e.g., implemented by any suitable hardware and/or software as known in the art. Subtractor 454 may subtract the value of output 452, from the value of output 442, to produce an output 455 including the evaluated noise error, denoted e1(n), corresponding to samples MIC1(n) and MIC21(n).
According to exemplary embodiments of the invention, estimator 408 may implement an adaptive algorithm to iteratively adapt the value of vector W1, e.g., based on the value of e1(n), as described below.
According to exemplary embodiments of the invention, the noise error, e1(n), corresponding to the n-th samples received from microphones 402 and 404 may be evaluated using the following equation:

e1(n)=d(n)−z1(n)  (22)

Substituting Equation 21 in Equation 22 may yield the following equation:




  
    
      
        
          
            e
            1
          
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            d
            ⁡
            
              (
              n
              )
            
          
          -
          
            
              ∑
              
                p
                =
                0
              
              
                S
                -
                1
              
            
            ⁢
            
              
                STF
                ⁡
                
                  (
                  p
                  )
                
              
              ⁢
              
                
                  ∑
                  
                    s
                    =
                    0
                  
                  
                    L
                    1
                  
                
                ⁢
                
                  
                    
                      W
                      1
                    
                    ⁡
                    
                      (
                      s
                      )
                    
                  
                  ⁢
                  MIC
                  ⁢
                  
                    
                  
                  ⁢
                  1
                  ⁢
                  
                    (
                    
                      n
                      -
                      s
                      -
                      p
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        23
        )
      
    
  




According to some exemplary embodiments of the invention, estimator 408 may iteratively adapt the value of vector W1, to reduce, e.g., minimize, the evaluated noise error e1(n). For example, estimator 408 may be able to iteratively adapt the value of vector W1 using the following equation:




  
    
      
        
          
            W
            1
          
          ⁡
          
            (
            
              n
              +
              1
            
            )
          
        
        =
        
          
            
              W
              1
            
            ⁡
            
              (
              n
              )
            
          
          -
          
            
              μ
              1
            
            ⁢
            
              
                ∑
                
                  s
                  =
                  0
                
                
                  S
                  -
                  1
                
              
              ⁢
              
                
                  STF
                  ⁡
                  
                    (
                    s
                    )
                  
                
                ⁢
                MIC
                ⁢
                
                  
                
                ⁢
                1
                ⁢
                
                  (
                  
                    n
                    -
                    s
                  
                  )
                
                ⁢
                
                  
                    e
                    1
                  
                  ⁡
                  
                    (
                    n
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        24
        )
      
    
  




wherein W1(n+1) denotes an adapted value of W1, W1(n) denotes the current value of W1, and μ1 denotes a predetermined convergence parameter corresponding to W1. For example, μ1 may be determined according the following condition:





  
    
      
        
          μ
          1
        
        <
        
          1
          
            2
            ⁢
            
              L
              1
            
          
        
      
    
    
      
        (
        25
        )
      
    
  




According to some exemplary embodiments of the invention, estimator 410 may iteratively adapt the value of vector W2 of the estimation function F2, based on an evaluated residual noise error corresponding to the combination of, e.g., the difference between, the evaluated noise error e1(n), and an estimation of the contribution of y2(n) to noise destructive pattern 418, as described below.
According to some exemplary embodiments of the invention, controller 400 may also include at least one secondary evaluation module 432 to evaluate the residual noise error, e.g., based on signal 422 and the evaluated noise error e1(n), as described below.
According to exemplary embodiments of the invention, module 432 may include a STF module 460 to apply a STF to signal 422 to generate an output 462 representing an estimation of a secondary part of th noise destructive pattern corresponding to signal 422. STF module 460 may include any suitable hardware and/or software, e.g., as known in the art, to apply a predetermined STF to signal 422. The STF may be predetermined, for example, based on characteristics of speaker 406, e.g., as known in the art. For example, the value, denoted z2(n), of output 462 may be calculated using the following equation:




  
    
      
        
          
            z
            2
          
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            ∑
            
              p
              =
              0
            
            
              S
              -
              1
            
          
          ⁢
          
            
              STF
              ⁡
              
                (
                p
                )
              
            
            ⁢
            
              
                y
                2
              
              ⁡
              
                (
                
                  n
                  -
                  p
                
                )
              
            
          
        
      
    
    
      
        (
        26
        )
      
    
  




Substituting Equation 18 in Equation 26 may yield the following equation:




  
    
      
        
          
            z
            2
          
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            ∑
            
              p
              =
              0
            
            
              S
              -
              1
            
          
          ⁢
          
            
              STF
              ⁡
              
                (
                p
                )
              
            
            ⁢
            
              
                ∑
                
                  x
                  =
                  0
                
                
                  L
                  2
                
              
              ⁢
              
                
                  
                    W
                    2
                  
                  ⁡
                  
                    (
                    s
                    )
                  
                
                ⁢
                MIC
                ⁢
                
                  
                
                ⁢
                21
                ⁢
                
                  (
                  
                    n
                    -
                    s
                    -
                    p
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        27
        )
      
    
  




According to exemplary embodiments of the invention, module 432 may also include a subtractor 464, e.g., implemented by any suitable hardware and/or software as known in the art. Subtractor 464 may subtract the value of output 462, from the value of output 452, to produce an output 466 including the evaluated residual noise error, denoted e2(n), corresponding to samples MIC1(n) and MIC21(n).
According to exemplary embodiments of the invention, estimator 410 may implement an adaptive algorithm to iteratively adapt the value of vector W2, e.g., based on the value of e2(n), as described below.
According to exemplary embodiments of the invention, the residual noise error, e2(n), corresponding to the n-th samples received from microphones 402 and 404 may be evaluated using the following equation:

e2(n)=e1(n)−z2(n)  (28)

Substituting Equations 23 and 27 in Equation 28 may yield the following equation:




  
    
      
        
          
            e
            2
          
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            d
            ⁡
            
              (
              n
              )
            
          
          -
          
            
              ∑
              
                p
                =
                0
              
              
                S
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                STF
                ⁡
                
                  (
                  p
                  )
                
              
              ⁢
              
                
                  ∑
                  
                    s
                    =
                    0
                  
                  
                    L
                    1
                  
                
                ⁢
                
                  
                
                ⁢
                
                  
                    
                      W
                      1
                    
                    ⁡
                    
                      (
                      s
                      )
                    
                  
                  ⁢
                  MIC
                  ⁢
                  
                    
                  
                  ⁢
                  1
                  ⁢
                  
                    (
                    
                      n
                      -
                      s
                      -
                      p
                    
                    )
                  
                
              
            
          
          -
          
            
              ∑
              
                p
                =
                0
              
              
                S
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                STF
                ⁡
                
                  (
                  p
                  )
                
              
              ⁢
              
                
                  ∑
                  
                    s
                    =
                    0
                  
                  
                    L
                    2
                  
                
                ⁢
                
                  
                
                ⁢
                
                  
                    
                      W
                      2
                    
                    ⁡
                    
                      (
                      s
                      )
                    
                  
                  ⁢
                  MIC
                  ⁢
                  
                    
                  
                  ⁢
                  21
                  ⁢
                  
                    (
                    
                      n
                      -
                      s
                      -
                      p
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        29
        )
      
    
  




According to some exemplary embodiments of the invention, estimator 410 may iteratively adapt the value of vector W2, to reduce, e.g., minimize, the evaluated residual noise error e2(n). For example, estimator 410 may be able to iteratively adapt one or more elements of vector W1 using the following equation:




  
    
      
        
          
            W
            2
          
          ⁡
          
            (
            
              n
              +
              1
            
            )
          
        
        =
        
          
            
              W
              2
            
            ⁡
            
              (
              n
              )
            
          
          -
          
            
              μ
              2
            
            ⁢
            
              
                ∑
                
                  p
                  =
                  0
                
                
                  S
                  -
                  1
                
              
              ⁢
              
                
              
              ⁢
              
                
                  STF
                  ⁡
                  
                    (
                    p
                    )
                  
                
                ⁢
                MIC
                ⁢
                
                  
                
                ⁢
                21
                ⁢
                
                  (
                  
                    n
                    -
                    s
                    -
                    p
                  
                  )
                
                ⁢
                
                  
                    e
                    2
                  
                  ⁡
                  
                    (
                    n
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        30
        )
      
    
  




wherein W2(n+1) denotes an adapted value of W2, W2(n) denotes the current value of W2, and μ2 denotes a predetermined convergence parameter corresponding to W2. For example, μ2 may be determined according the following condition:





  
    
      
        
          μ
          2
        
        <
        
          1
          
            2
            ⁢
            
              L
              2
            
          
        
      
    
    
      
        (
        31
        )
      
    
  




Some of the embodiments described above may refer to ANC systems implementing a controller, e.g., controller 400, able to control an acoustic transducer, e.g., speaker 406, to generate a noise destructive pattern based on a combination of an a primary noise signal of a primary acoustic sensor, e.g., microphone 402, and a secondary noise signal of a secondary acoustic sensor, e.g., microphone 404. However, it will be appreciated by those skilled in the art that according to other embodiments of the invention, these systems may be modified to implement one or more additional secondary acoustic sensors. For example, controller 400 may be modified to include an additional plurality of secondary estimators to receive one or more primary noise signals of the one or more additional secondary microphones, respectively. For example, an i-th estimator of the additional secondary estimators may generate, for example, an output, denoted yi(n), corresponding to the following equation:




  
    
      
        
          
            y
            i
          
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            ∑
            
              s
              =
              0
            
            
              L
              i
            
          
          ⁢
          
            
          
          ⁢
          
            
              
                W
                i
              
              ⁡
              
                (
                s
                )
              
            
            ⁢
            
              MICi
              ⁡
              
                (
                
                  n
                  -
                  s
                
                )
              
            
          
        
      
    
    
      
        (
        32
        )
      
    
  




wherein Wi denotes a predetermined prediction filter (PF) vector of length Li corresponding to the i-th estimator, and MICi denotes the output of the i-th additional secondary microphone.

Controller 400 may also be modified to include one or more additional residual noise error evaluators to evaluate a residual noise error, e.g., in analogy to evaluator 410. For example, an i-th residual error evaluator may evaluate the i-th residual noise error, ei(n), using the following equation:




  
    
      
        
          
            e
            i
          
          ⁡
          
            (
            n
            )
          
        
        =
        
          
            
              e
              
                i
                -
                1
              
            
            ⁡
            
              (
              n
              )
            
          
          -
          
            
              ∑
              
                p
                =
                0
              
              
                S
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                STF
                ⁡
                
                  (
                  p
                  )
                
              
              ⁢
              
                
                  ∑
                  
                    s
                    =
                    0
                  
                  
                    L
                    i
                  
                
                ⁢
                
                  
                
                ⁢
                
                  
                    
                      W
                      i
                    
                    ⁡
                    
                      (
                      s
                      )
                    
                  
                  ⁢
                  
                    MICi
                    ⁡
                    
                      (
                      
                        n
                        -
                        s
                        -
                        p
                      
                      )
                    
                  
                
              
            
          
        
      
    
    
      
        (
        33
        )
      
    
  




According to some exemplary embodiments, an i-th estimator of the additional estimators may iteratively adapt the value of the vector Wi, e.g., using the following equation:




  
    
      
        
          
            W
            i
          
          ⁡
          
            (
            
              n
              +
              1
            
            )
          
        
        =
        
          
            
              W
              i
            
            ⁡
            
              (
              n
              )
            
          
          -
          
            
              μ
              i
            
            ⁢
            
              
                ∑
                
                  s
                  =
                  0
                
                
                  S
                  -
                  1
                
              
              ⁢
              
                
              
              ⁢
              
                
                  STF
                  ⁡
                  
                    (
                    s
                    )
                  
                
                ⁢
                
                  
                    MICi
                    ⁡
                    
                      (
                      
                        n
                        -
                        s
                      
                      )
                    
                  
                  i
                
              
            
          
        
      
    
    
      
        (
        34
        )
      
    
  




wherein Wi(n+1) denotes an adapted value of Wi, Wi(in) denotes the current value of Wi, and μi denotes a predetermined convergence parameter corresponding to Wi. For example, μi may be determined according the following condition:





  
    
      
        
          μ
          i
        
        <
        
          1
          
            2
            ⁢
            
              L
              i
            
          
        
      
    
    
      
        (
        35
        )
      
    
  




Some of the embodiments described above may refer to ANC systems implementing a controller, e.g., controller 400, including one or more estimators, e.g., estimators 408 and/or 410, to apply an adaptive linear estimation algorithm to one or more respective noise signals, e.g., outputs 412 and/or 420. However, it will be appreciated by those skilled in the art that according to other embodiments of the invention, these systems may be modified to implement one or more estimators to apply an adaptive non-linear estimation algorithm to one or more respective noise signals. For example, controller 400 may be modified to implement one or more RBF estimation algorithms, e.g., in analogy to controller 200 (FIG. 2).
Embodiments of the present invention may be implemented by software, by hardware, or by any combination of software and/or hardware as may be suitable for specific applications or in accordance with specific design requirements. Embodiments of the present invention may include modules, units and sub-units, which may be separate of each other or combined together, in whole or in part, and may be implemented using specific, multi-purpose or general processors, or devices as are known in the art. Some embodiments of the present invention may include buffers, registers, storage units and/or memory units, for temporary or long-term storage of data and/or in order to facilitate the operation of a specific embodiment.
While certain features of the invention have been illustrated and described herein, many modifications, substitutions, changes, and equivalents may occur to those of ordinary skill in the art. It is, therefore, to be understood that the appended claims are intended to cover all such modifications and changes as fall within the true spirit of the invention.

## Claims

1. An active noise control system for controlling noise produced by a noise source, said system comprising:
an acoustic sensor to sense a noise pattern and to produce a noise signal corresponding to the sensed noise pattern;
an estimator to produce a predicted noise signal by applying a non-liner estimation function to said noise signal, wherein the predicted noise signal includes an estimation of a predicted sample of the noise signal, which is successive to a current sample of the noise signal, and wherein the estimator is to estimate the predicted sample by applying the estimation function to the current sample and to one or more samples preceding the current sample of the noise signal; and
an acoustic transducer to produce a noise destructive pattern based on said predicted noise signal,
wherein the noise destructive pattern has a non-linear relationship to the noise pattern sensed by the acoustic sensor.

2. The system of claim 1, wherein said estimator is able to adapt one or more parameters of said estimation function based on a noise error at a predetermined location.

3. The system of claim 2, wherein said noise error comprises an anticipated destructive interference between said noise pattern and said noise destructive pattern at said predetermined location.

4. The system of claim 2 comprising an error-sensing microphone to sense said noise error at said predetermined location.

5. The system of claim 2 comprising an error evaluator to evaluate said noise error based on said noise signal and said predicted noise signal.

6. The system of claim 5, wherein said error evaluator comprises:
a speaker transfer function module to produce an estimation of said noise destructive pattern by applying a speaker transfer function to said predicted noise signal;
a modulation transfer function module to produce an estimation of said noise pattern at said predetermined location by applying a modulation transfer function to said noise signal; and
a subtractor to subtract the estimation of said noise destructive pattern from the estimation of said noise pattern.

7. The system of claim 2, wherein said estimator is able to adapt said one or more parameters based on a predetermined criterion.

8. The system of any one of claim 7, wherein said estimator is able to reduce said error value by adapting said one or more parameters.

9. The system of claim 8, wherein said adaptive estimator is able to minimize said error value by adapting said one or more parameters.

10. The system of claim 2, wherein said one or more parameters comprise at least one parameter selected from the group consisting of a center parameter, an effective radius parameter, and an intensity parameter.

11. The system of claim 10, wherein said estimator is able to adapt said center parameter based on the following equation:




  
    
      c
      k
    
    ⁡
    
      (
      
        n
        +
        1
      
      )
    
  
  =
  
    
      
        c
        k
      
      ⁡
      
        (
        n
        )
      
    
    -
    
      
        μ
        c
      
      ⁢
      
        e
        ⁡
        
          (
          n
          )
        
      
      ⁢
      
        w
        k
      
      ⁢
      
        
          ∑
          
            s
            =
            0
          
          
            S
            -
            1
          
        
        ⁢
        
          
        
        ⁢
        
          
            STF
            ⁡
            
              (
              s
              )
            
          
          ⁢
          
            
              f
              k
            
            ⁡
            
              [
              
                n
                -
                s
              
              ]
            
          
          ⁢
          
            (
            
              
                1
                
                  υ
                  k
                
              
              ⁢
              
                
                  ∑
                  
                    i
                    =
                    0
                  
                  
                    L
                    -
                    1
                  
                
                ⁢
                
                  
                
                ⁢
                
                  (
                  
                    
                      x
                      ⁡
                      
                        (
                        
                          n
                          -
                          i
                        
                        )
                      
                    
                    -
                    
                      
                        c
                        k
                      
                      ⁡
                      
                        (
                        i
                        )
                      
                    
                  
                  )
                
              
            
            )
          
        
      
    
  




wherein ck(n+1) denotes an adapted value of said center parameter,
ck(n) denotes a current value of said center parameter,
wk denotes said intensity parameter,
L denotes a predetermined number of samples of said noise signal,
STF denotes a predetermined speaker transfer function,
S denotes a predetermined speaker transfer function frequency parameter,
μc denotes a predetermined convergence parameter corresponding to said center parameter,
vk denotes said effective radius parameter,
e(n) denotes said noise error,
fk denotes a predetermined function, and
x(n) denotes an n-th sample of said noise signal.

12. The system of claim 10, wherein said estimator is able to adapt said effective radius parameter based on the following equation:




  
    
      υ
      k
    
    ⁡
    
      (
      
        n
        +
        1
      
      )
    
  
  =
  
    
      
        υ
        k
      
      ⁡
      
        (
        n
        )
      
    
    -
    
      
        μ
        υ
      
      ⁢
      
        e
        ⁡
        
          (
          n
          )
        
      
      ⁢
      
        w
        k
      
      ⁢
      
        
          ∑
          
            s
            =
            0
          
          
            S
            -
            1
          
        
        ⁢
        
          
        
        ⁢
        
          
            STF
            ⁡
            
              (
              s
              )
            
          
          ⁢
          
            
              f
              k
            
            ⁡
            
              [
              
                n
                -
                s
              
              ]
            
          
          ⁢
          
            1
            
              
                (
                
                  υ
                  k
                
                )
              
              2
            
          
          ⁢
          
            
              ∑
              
                i
                =
                0
              
              
                L
                -
                1
              
            
            ⁢
            
              
            
            ⁢
            
              
                (
                
                  
                    x
                    ⁡
                    
                      (
                      
                        n
                        -
                        i
                      
                      )
                    
                  
                  -
                  
                    
                      c
                      k
                    
                    ⁡
                    
                      (
                      i
                      )
                    
                  
                
                )
              
              2
            
          
        
      
    
  




wherein vk(n+1) denotes an adapted value of said effective radius parameter,
vk(n) denotes a current value of said effective radius parameter,
wk denotes said intensity parameter,
L denotes a predetermined number of samples of said noise signal,
STF denotes a predetermined speaker transfer function,
S denotes a predetermined speaker transfer function frequency parameter,
μv, denotes a predetermined convergence parameter corresponding to said effective radius parameter,
ck denotes said center parameter,
e(n) denotes said noise error,
fk denotes a predetermined function, and
x(n) denotes an n-th sample of said noise signal.

13. The system of claim 10, wherein said estimator is able to adapt said intensity parameter based on the following equation:




  
    
      w
      k
    
    ⁡
    
      (
      
        n
        +
        1
      
      )
    
  
  =
  
    
      
        w
        k
      
      ⁡
      
        (
        n
        )
      
    
    -
    
      
        μ
        w
      
      ⁢
      
        e
        ⁡
        
          (
          n
          )
        
      
      ⁢
      
        
          ∑
          
            s
            =
            0
          
          
            S
            -
            1
          
        
        ⁢
        
          
        
        ⁢
        
          
            STF
            ⁡
            
              (
              s
              )
            
          
          ⁢
          
            
              f
              k
            
            ⁡
            
              [
              
                n
                -
                s
              
              ]
            
          
        
      
    
  




wherein wk(n+1) denotes an adapted value of said intensity parameter,
wk(n) denotes a current value of said intensity parameter,
wk denotes said intensity parameter,
L denotes a predetermined number of samples of said noise signal,
STF denotes a predetermined speaker transfer function,
S denotes a predetermined speaker transfer function frequency parameter,
μw, denotes a predetermined convergence parameter corresponding to said intensity parameter,
fk denotes a predetermined function, and
x(n) denotes an n-th sample of said noise signal.

14. The system of claim 1, wherein said estimation function comprises a non-linear estimation function,
wherein the estimator is able to estimate a noise error corresponding to an anticipated destructive interference between a pattern of the noise and the noise destructive pattern at a predetermined location, wherein said predetermined location is distinct from a location of said acoustic sensor.

15. The system of claim 14, wherein said non-linear function comprises a radial basis function.

16. The system of claim 1, wherein said acoustic sensor comprises a microphone, and wherein the noise destructive pattern produced by the acoustic transducer has an exponential relationship to the noise pattern sensed by the acoustic sensor.

17. The system of claim 1, wherein said acoustic transducer comprises a speaker,
wherein said acoustic sensor comprises an array of two or more microphones,
wherein the two or more microphones are located in two or more, respective, locations,
wherein the two or more microphones are adapted to achieve coherence between the sensed noise pattern and the noise produced by the noise source, by taking into account at least one or more of:
geometric structure of a path between said microphones and the noise source;
aerodynamic attributes of the path between said microphones and the noise source;
surface roughness along the path between said microphones and the noise source;
turbulent airflow along the path between said microphones and the noise source;
formation of acoustic signals along the path between said microphones and the noise source.

18. An active noise control system for controlling a noise produced by a noise source, said system comprising:
a primary acoustic sensor to sense a noise pattern and to produce a corresponding primary noise signal;
at least one secondary acoustic sensor to sense a residual noise pattern and to produce at least one secondary noise signal corresponding to the residual noise pattern sensed by said at least one secondary acoustic sensor, respectively,
wherein said at least one secondary acoustic sensor is separated from said noise source by a distance larger than a distance between said primary acoustic sensor and said noise source; and
a controller functionally associated with an acoustic transducer and at least one estimator to produce a predicted noise signal,
wherein the predicted noise signal includes an estimation of a predicted sample of at least one sampled signal of the primary noise signal and the secondary noise signal, which is successive to a current sample of the sampled signal, and wherein the estimator is to estimate the predicted sample by applying at least one non-linear estimation function to the current sample and to one or more samples preceding the current sample of the sampled signal,
wherein said controller is adapted to produce a noise destructive pattern based on said primary noise signal and said at least one secondary noise signal and said predicted noise signal,
and wherein the noise destructive pattern produced by the controller has a non-linear relationship to the noise pattern sensed by the primary acoustic sensor.

19. The system of claim 18, wherein said at least one estimator includes a primary estimator adapted to produce a predicted primary signal by applying a primary estimation function to said primary noise signal and at least one secondary estimator to produce at least one predicted secondary signal by applying at least one secondary estimation function to said at least one secondary noise signal, respectively.

20. The system of claim 19, wherein said primary estimator is able to iteratively adapt one or more parameters of said primary estimation function based on a noise error.

21. The system of claim 19, wherein said at least one secondary estimator is able to iteratively adapt one or more parameters of said at least one secondary estimation function, respectively, based on a noise error.

22. The system claim 19, wherein said controller is able to control said acoustic transducer based on a combination of said predicted primary signal and said at least one predicted secondary signal.

23. The system of claim 22, wherein said controller is able to control said acoustic transducer based on the sum of said predicted primary signal and said at least one predicted secondary signal.

24. The system claim 20, wherein said controller comprises a noise error evaluator to evaluate a noise error corresponding to an anticipated destructive interference between a pattern of the noise and the noise destructive pattern at a predetermined location, wherein said predetermined location is distinct from locations of said primary and secondary acoustic sensors.

25. The system of claim 24, wherein said noise error evaluator is able to evaluate said noise error based on said primary noise signal, said at least one secondary noise signal and said predicted primary signal.

26. The system of claim 25, wherein said noise error evaluator comprises:
a speaker transfer function module to produce an estimation of a primary part of said noise destructive pattern corresponding to said predicted primary signal by applying a speaker transfer function to said predicted primary signal;
a modulation transfer function module to produce an estimation of said noise pattern by applying a modulation transfer function to a combination of said primary noise signal and said at least one secondary noise signal; and
a subtractor to subtract the estimation of the primary part of said noise destructive pattern from the estimation of said noise pattern.

27. The system of claim 24, wherein said controller comprises at least one residual noise evaluator to evaluate at least one residual noise.

28. The system of claim 27, wherein said at least one residual noise evaluator is able to evaluate said residual noise based on said noise error and said at least one predicted secondary signal, respectively.

29. The system of claim 28, wherein said residual error evaluator comprises:
a speaker transfer function module to produce an estimation of a secondary part of said noise destructive pattern corresponding to said predicted secondary signal by applying a speaker transfer function to said predicted secondary signal; and
a subtractor to subtract the estimation of the secondary part of said noise destructive pattern from said noise error.

30. The system of claim 18, wherein at least one of said primary acoustic sensor and said at least one secondary acoustic sensor comprises a microphone, and wherein the noise destructive pattern produced by the acoustic transducer has an exponential relationship to the noise pattern sensed by the primary acoustic sensor.

31. The system of claim 18, wherein said acoustic transducer comprises a speaker,
wherein said primary acoustic sensor comprises an array of two or more microphones,
wherein the two or more microphones are located in two or more, respective, locations,
wherein the two or more microphones are adapted to achieve coherence between the sensed noise pattern and the noise produced by the noise source, by taking into account at least one or more of:
geometric structure of a path between said microphones and the noise source;
aerodynamic attributes of the path between said microphones and the noise source;
surface roughness along the path between said microphones and the noise source;
turbulent airflow along the path between said microphones and the noise source;
formation of acoustic signals along the path between said microphones and the noise source.