Abstract:
An active noise controller using an adaptive notch filter performs a coefficient update operation by using only a signal obtained by processing an error signal with ½ cycle of the frequency of noise to be reduced for a predetermined period. This eliminates the necessity of product operation for the coefficient update operation, and in turn significantly reduces an operation load.

Description:
TECHNICAL FIELD 
       [0001]    The present invention relates to an active noise reduction device for actively reducing vibratory noise generated by a rotating machine such as an engine of a vehicle. 
       BACKGROUND ART 
       [0002]    For a conventional active noise reduction device, patent document 1 discloses an adaptive control method using an adaptive notch filter. 
         [0003]      FIG. 7  is a block diagram showing a configuration of a conventional active noise reduction device. In  FIG. 7 , discrete operation processing unit  115  executes a discrete operation for implementing the active noise reduction device. Engine rotational speed detector  101  outputs a pulse train as engine pulse p, which has a frequency proportional to an engine rotational speed. For example, engine pulse p is generated by extracting an output of a crank angle sensor. Frequency detecting unit  102  calculates and outputs noise frequency f based on engine pulse p. Basic signal generating unit  116  has sine wave table  103  holding, in a memory, a value of each point defined by equally dividing one cycle of a sine wave into predetermined parts. Selection unit  117  selects data from sine wave table  103 , and generates and outputs basic sine wave signal x 1 [ n ] and basic cosine wave signal x 2 [ n ] having the same frequency as noise frequency f. Reference signal generating unit  118  generates and outputs reference sine wave signal r 1 [ n ] and reference cosine wave signal r 2 [ n ] by using basic sine wave signal correction value table  119  (a basic sine wave signal correction value in a frequency of f (Hz) is referred to as C 1 [ f ]) and basic cosine wave signal correction value table  120  (a basic cosine wave signal correction value in a frequency of f (Hz) is referred to as C 2 [ f ]) simulating a transfer characteristic value from loudspeaker  110  to microphone  111 . 
         [0004]    First one-tap digital filter  107  filters basic sine wave signal x 1 [ n ] based on filter coefficient W 1 [ n ] held therein and generates first control signal y 1 [ n ]. Second one-tap digital filter  108  filters basic cosine wave signal x 2 [ n ] based on filter coefficient W 2 [ n ] held therein and generates second control signal y 2 [ n ]. Power amplifier  109  amplifies a signal obtained by adding first control signal y 1 [ n ] and second control signal y 2 [ n ] to each other. Loudspeaker  110  outputs an output signal from power amplifier  109  as a noise cancellation sound. Microphone  111  detects a sound generated as a result of interference of a noise and a noise cancellation sound as error signal ε[n]. 
         [0005]    First adaptive control algorithm operating unit  112  sequentially updates filter coefficient W 1 [ n ] by, for example, an LMS (Least Mean Square) algorithm that is a kind of steepest descent method based on reference sine wave signal r 1 [ n ] and error signal ε[n]. Similarly, second adaptive control algorithm operating unit  113  sequentially updates filter coefficient W 2 [ n ] based on reference cosine wave signal r 2 [ n ] and error signal ε[n]. 
         [0006]    Sequential updating formulae of coefficients W 1  and W 2  are represented by the following equations. 
         [0000]        W 1 [n+ 1 ]=W 1 [n]−μ×r 1 [n]×ε[n]   (1) 
         [0000]        W 2 [n+ 1 ]=W 2 [n]−μ×r 2 [n]×ε[n]   (2) 
         [0007]    Herein, μ denotes a constant called a convergence coefficient and relates to a time period in which coefficients W 1  and W 2  converge on suitable values. 
         [0008]    By repeating the above-mentioned processing with a predetermined cycle, noises can be reduced. 
         [0009]    However, in generating reference sine wave signal r 1 [ n ] and reference cosine wave signal r 2 [ n ], the above-mentioned conventional configuration carries out a product-sum operation of basic sine wave signal x 1 [ n ] and basic sine wave signal correction value C 1 [ k ] as well as a product-sum operation of basic cosine wave signal x 2 [ n ] and basic cosine wave signal correction value C 2 [ k ]. It needs two product operations in order to generate the respective reference signals. Furthermore, in order to determine coefficients W 1  and W 2  of the one-tap digital filters, the above-obtained reference sine wave signal r 1 [ n ] and reference cosine wave signal r 2 [ n ] need to be multiplied by convergence coefficient μ and error signal ε[n], respectively. Therefore, two product operations are required (see, formulae (1) and (2)). That is to say, in order to determine coefficients W 1  and W 2  of the respective one-tap digital filters, four product operations are required. As a result, operation load increases. 
         [0010]    [Patent document 1] Unexamined Japanese Patent Publication No. 2004-361721 
       SUMMARY OF THE INVENTION 
       [0011]    The present invention provides an active noise controller in which an operation load necessary for noise-suppression control is reduced by minimizing the execution of product operation. 
         [0012]    An active noise controller in accordance with the present invention includes a frequency detecting unit for detecting a frequency of noise to be controlled from a noise source; a sine wave generating unit for generating a sine wave having the same frequency as the frequency of the noise detected by the frequency detecting unit; a cosine wave generating unit for generating a cosine wave having the same frequency as the frequency of the noise detected by the frequency detecting unit; a first one-tap digital filter to which a sine wave signal from the sine wave generating unit is input; a second one-tap digital filter to which a cosine wave signal from the cosine wave generating unit is input; an interference signal generating unit to which a noise control signal obtained by adding an output from the first one-tap digital filter and an output from the second one-tap digital filter to each other is input and which outputs an interference signal to cause interference with the noise to be controlled from the noise source; an error signal detecting unit for detecting an error signal generated as a result of the interference of the interference signal output from the interference signal generating unit and the noise to be controlled from the noise source; a first coefficient updating unit for updating a filter coefficient of the first one-tap digital filter; and a second coefficient updating unit for updating a filter coefficient of the second one-tap digital filter. The first coefficient updating unit and the second coefficient updating unit update the coefficients of the first one-tap digital filter and the second one-tap digital filter so that noise in the error signal detecting unit is reduced, based on a coefficient update signal obtained by processing the error signal from the error signal detecting unit with ½ cycle of the frequency of the noise to be controlled for a predetermined period. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0013]      FIG. 1  is a block diagram showing an active noise controller in accordance with a first exemplary embodiment of the present invention. 
           [0014]      FIG. 2  is a graph showing an example of a sine wave table in the active noise controller in accordance with the first exemplary embodiment of the present invention. 
           [0015]      FIG. 3  shows an example of a sine wave table in the active noise controller in accordance with the first exemplary embodiment of the present invention. 
           [0016]      FIG. 4  is a graph showing an example of a transfer characteristic from a loudspeaker to a microphone of the active noise controller in accordance with the first exemplary embodiment of the present invention. 
           [0017]      FIG. 5  shows an example of a characteristic table having a lower limit passing point and an upper limit passing point corresponding to a transfer characteristic from the speaker to the microphone shown in  FIG. 4  in the active noise controller in accordance with the first exemplary embodiment of the present invention. 
           [0018]      FIG. 6A  is a graph showing a time-base waveform of a square wave formed by the active noise controller in accordance with the first exemplary embodiment of the present invention. 
           [0019]      FIG. 6B  is a graph showing a harmonic analysis of a square wave processed by the active noise controller in accordance with the first exemplary embodiment of the present invention. 
           [0020]      FIG. 7  is a block diagram showing a configuration of a conventional active noise reduction device. 
       
    
    
     REFERENCE MARKS IN THE DRAWINGS 
       [0000]    
       
           1  engine rotational speed detector 
           2  frequency detecting unit 
           3  sine wave table 
           4  characteristic table 
           5  sine wave generating unit 
           6  cosine wave generating unit 
           7  first one-tap digital filter 
           8  second one-tap digital filter 
           9  power amplifier 
           10  loudspeaker (interference signal generating unit) 
           11  microphone (error signal detecting unit) 
           12  first adaptive control algorithm operating unit (first coefficient updating unit) 
           13  second adaptive control algorithm operating unit (second coefficient updating unit) 
           14  coefficient update signal generating unit 
           15  discrete operation processing unit 
       
     
       DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
     First Exemplary Embodiment 
       [0036]    Hereinafter, an active noise controller in accordance with a first exemplary embodiment of the present invention is described with reference to drawings. 
         [0037]      FIG. 1  is a block diagram showing an active noise controller in accordance with the first exemplary embodiment of the present invention. In  FIG. 1 , engine rotational speed detector  1  outputs a pulse train as engine pulse p, which has a frequency proportional to a rotational speed of a vehicle-mounted engine as a noise source. Frequency detecting unit  2  calculates and outputs frequency f [Hz] of noise to be controlled from engine pulse p. Sine wave table  3  including discrete sine wave data holds, in a memory, a sine value at each point defined by equally dividing one cycle of the sine wave into N parts. Sine wave generating unit  5  reads out data in a position at readout point P[n] from the sine wave table in every sampling cycle, and generates basic sine wave signal x 1 [ n ]. At this time, the difference (P[n+1]−P[n]) between P[n+1] that is a value subsequent to the readout point and P[n] that is a value of the current readout point is N×f×T when the sampling cycle is denoted by T and the frequency of noise to be controlled is denoted by f. 
         [0038]    Similarly, cosine wave generating unit  6  reads out a point preceding the readout point of sine wave generating unit  5  by only N/4, that is, data in a position of P[n]+N/4 and generates basic cosine wave signal x 2 [ n ]. At this time, when each of readout points P[n] and P[n]+N/4 exceeds N, a point obtained by subtracting N from the readout point should be a new readout point. 
         [0039]    Characteristic table  4  holds lower limit point PP 1 [ f ] and upper limit point PP 2 [ f ] of the readout point for processing error signal ε[n] in every frequency based on the phase characteristic as the transfer characteristic from loudspeaker  10  to microphone  11 . 
         [0040]    Coefficient update error signal generating unit  14  reads lower limit point PP 1 [ f ] and upper limit point PP 2 [ f ] in frequency f of noise to be controlled from characteristic table  4  based on frequency f of noise to be controlled; based on them, processes error signal ε[n] detected by microphone  11 ; and generates coefficient update signals ε 1 [ n ] and ε 2   [n].    
         [0041]    Herein, ε 1 [ n ] satisfies the following equations: 
         [0000]      ε1[n]=ε[n] when PP1[f]≦P[n]≦PP2[f] is satisfied, and 
         [0000]      ε1 [n]=−ε[n]  when  PP 1 [f]+N/ 2 ≦P[n]≦PP 2 [f]+N/ 2 is satisfied. 
         [0042]    In the cases other than the above, the following equation is satisfied. 
         [0000]      ε1[n]=0  (3) 
         [0043]    Furthermore, ε[n] satisfies the following equations: 
         [0000]      ε2 [n]=ε[n ] when  PP 1 [f]+N/ 4 ≦P[n]≦PP 2 +N/ 4 [f ] is satisfied, and 
         [0000]      ε2 [n]=−ε[n ] when  NP 1 [f]+N× ¾ ≦P[n]≦NP 2 +N× ¾ [f ] is satisfied. 
         [0044]    In the cases other than the above, the following equation is satisfied. 
         [0000]      ε2[n]=0  (4) 
         [0045]    Next, first one-tap digital filter  7  holds first filter coefficient W 1 [ n ] therein, and outputs first control signal y 1 [ n ] based on basic sine wave signal x 1 [ n ] and first filter coefficient W 1 [ n ]. Second one-tap digital filter  8  holds second filter coefficient W 2 [ n ] therein, and outputs second control signal y 2 [ n ] based on basic cosine wave signal x 2 [ n ] and second filter coefficient W 2 [ n].    
         [0046]    Power amplifier  9  amplifies a noise control signal obtained by adding first control signal y 1 [ n ] and second control signal y 2 [ n ] to each other. Loudspeaker  10  outputs an output signal from power amplifier  9  as a noise cancellation sound. Microphone  11  detects a sound, which is generated as a result of interference of noise caused by engine vibration to be controlled and the noise cancellation sound, as error signal ε[n]. 
         [0047]    First adaptive control algorithm operating unit  12  as a first coefficient updating unit sequentially updates filter coefficient W 1 [ n ] of first one-tap digital filter  7  by using coefficient update signal ε 1 [ n ]. Second adaptive control algorithm operating unit  13  as a second coefficient updating unit sequentially updates filter coefficient W 2 [ n ] of second one-tap digital filter  8  by using coefficient update signal ε[n]. Thus, discrete operation processing unit  15  includes software. 
         [0048]    Next, a specific operation of the device is described. 
         [0049]    The generation of basic sine wave signal x 1 [ n ], generation of basic cosine wave signal x 2 [ n ], generation of first control signal y 1 [ n ], generation of second control signal y 2 [ n ], detection of error signal ε[n], update of first filter coefficient W 1 [ n ], and update of second filter coefficient W 2 [ n ] are all executed in the same cycle. In the following description, this cycle is defined as T (second). 
         [0050]    For example, frequency detecting unit  2  generates an interrupt in every rising edge of engine pulse p, measures a time period between the rising edges, and calculates frequency f of noise to be controlled based on the measurement results. Sine wave table  3  holds, in the memory, discrete data of the sine value at each point defined by equally dividing one cycle of the sine wave into N parts. When an arrangement that stores sine values from 0th point to (N−1)th point is represented by z[m] (0≦m&lt;N), the following relational expression (5) is satisfied. 
         [0000]        z[m ]=sin(360 °×m/N )  (5) 
         [0051]      FIG. 2  is a graph showing an example of a sine wave table in the active noise controller in accordance with the first exemplary embodiment of the present invention, which is visually shown for easy understanding.  FIG. 3  shows an example of the sine wave table in the active noise controller by using an example of numeric values as stored in a microprocessor for operation in accordance with the first exemplary embodiment of the present invention.  FIGS. 2 and 3  show an example in which N is 3000. 
         [0052]    Characteristic table  4  holds, in the memory, lower limit point PP 1 [ f ] and upper limit point PP 2 [ f ] for processing error signal ε[n] (f denotes a frequency (Hz)) based on the phase characteristic as a transfer characteristic from loudspeaker  10  to microphone  11 . 
         [0053]    When the phase characteristic at f(Hz) is θ[f] (°), the following relational expressions (6) are satisfied. 
         [0000]        PP 1 [f]=N×θ[f]/ 360+α 
         [0000]        PP 2 [f]=PP 1 [f]+β   (6) 
         [0000]    Herein, α and β may be any positive constants but they need to satisfy α&lt;N/4 and β+α&lt;N/4. 
         [0054]      FIG. 4  is a graph showing an example of a transfer characteristic from a loudspeaker to a microphone of the active noise controller in accordance with the first exemplary embodiment of the present invention.  FIG. 5  shows an example of a characteristic table having a lower limit passing point and an upper limit passing point corresponding to the transfer characteristic from the loudspeaker to the microphone shown in  FIG. 4  in the active noise controller in accordance with the first exemplary embodiment of the present invention.  FIG. 5  shows lower limit PP 1 [ f ] of the readout point and upper limit PP 2  [ f ] of the readout point in which N is 3000 and frequency f of noise to be controlled ranges from 30 Hz to 100 Hz. 
         [0055]    Sine wave generating unit  5  stores the current readout position P[n] of sine wave table  3  in the memory, and moves the current readout position at every cycle by formula (7) based on frequency f of noise to be controlled. 
         [0000]        P[n+ 1 ]=P[n]+N×f×T   (7) 
         [0000]    When the calculation result in the right side of the formula (7) is not less than N, a value obtained by subtracting N from the calculation result in the right side of the formula (7) is defined as P[n+1]. 
         [0056]    At the same time, sine wave generating unit  5  generates basic sine wave signal x 1 [ n ] having the same frequency as frequency f of noise to be controlled by the following formulae (8) and (9). 
         [0000]        i× 1 =P[n]   (8) 
         [0000]        x 1 [n]=z[i× 1]  (9) 
         [0000]    When the calculation result in the right side of the formula (8) is not less than N, a value obtained by subtracting N from the calculation result in the right side of the formula (8) is defined as i×1. 
         [0057]    Furthermore, cosine wave generating unit  6  generates basic cosine wave signal x 2 [ n ] having the same frequency as frequency f of noise to be controlled and preceding from basic sine wave signal x 1 [ n ] by ¼ by the following formulae (10) and (11). 
         [0000]        i× 2 =P[n]+N/ 4  (10) 
         [0000]        x 2 [n]=z[i× 2]  (11) 
         [0000]    When the calculation result in the right side of the formula (10) is not less than N, a value obtained by subtracting N from the calculation result in the right side of the formula (10) is defined as i×2. 
         [0058]    At the same time, coefficient update signal generating unit  14  reads lower limit point PP 1 [ f ] and upper limit point PP 2 [ f ] in frequency f of noise to be controlled from characteristic table  4  based on frequency f of noise to be controlled; based on them, processes error signal ε[n] detected by microphone  11  by formulae (3) and (4); and generates coefficient update signals ε 1 [ n ] and ε 2 [ n ], respectively. 
         [0059]    First and second one-tap digital filters  7  and  8  generate first and second control signals y 1 [ n ] and y 2 [ n ] by formulae (12) and (13), respectively. 
         [0000]        y 1 [n]=W 1 [n]×x 1 [n]   (12) 
         [0000]        y 2 [n]=W 2 [n]×x 2 [n]   (13) 
         [0060]    First and second adaptive control algorithm operating units  12  and  13  update filter coefficients W 1 [ n ] and W 2 [ n ] held by first and second one-tap digital filters  7  and  8  by formulae (14) and (15), respectively. 
         [0000]        W 1 [n+ 1 ]=W 1 [n]−ε 1 [n]   (14) 
         [0000]        W 2 [n+ 1 ]=W 2 [n]−ε 2 [n]   (15) 
         [0061]    Filter coefficient W 1 [ n ] and filter coefficient W 2 [ n ] are allowed to converge by using the above-mentioned procedure. Thus, noise to be controlled can be reduced. 
         [0062]    Herein, a mechanism for reducing noise having a frequency to be controlled by coefficient updating formulae (14) and (15) is described. 
         [0063]    A noise controller described in a conventional example sequentially updates filter coefficients W 1 [ n ] and W 2 [ n ] based on an LMS (Least Mean Square) algorithm. The followings are the updating formulae. 
         [0000]        W 1 [n+ 1 ]=W 1 [n]−μ×r 1 [n]×ε[n]   (1) 
         [0000]        W 2 [n+ 1 ]=W 2 [n]−μ×r 2 [n]×ε[n]   (2) 
         [0064]    Thus, in general, the product of error signal ε[n] and a sine wave signal and a cosine wave signal having a frequency of noise to be controlled as reference sine wave signal r 1 [ n ] and reference cosine wave signal r 2 [ n ] is used. This uses the orthogonal property between the sine wave and the cosine wave. In the sequential update for a long period of time (that is, n→∞), the product of frequency components having the same frequency as that of reference sine wave signal r 1  and reference cosine wave signal r 2  in error signal ε is accumulated, and the accumulated values of the product of the other frequency components becomes 0. Thus, W 1 [ n ] and W 2 [ n ] are updated so that frequency components having the same frequency as that of reference sine signal r 1  and reference cosine signal r 2  in error signal ε are reduced. Finally, when the frequency components having the same frequency as the frequency of the reference sine wave signal and the reference cosine wave signal in the error signal ε become 0, and average coefficient update of W 1 [ n ] and W 2 [ n ] becomes 0, and W 1 [ n ] and W 2 [ n ] converge. 
         [0065]    On the other hand, the present invention does not use so-called reference signals (r 1 [ n ] and r 2 [ n ]). The present invention updates coefficients by using only coefficient update signals ε 1 [ n ] and ε 2 [ n ] generated from error signal ε[n] by formulae (3) and (4). 
         [0066]    The ε 1 [ n ] and ε 2 [ n ] can be also represented as follows: 
         [0000]      ε1 [n]= 1 ×ε[n ] when  PP 1 [f]≦P[n]≦PP 2 [f ] is satisfied, and 
         [0000]      ε1 [n]=− 1 ×ε[n ] when  PP 1 [f]+N/ 2 ≦P[n]≦PP 2 [f]+N/ 2 is satisfied. 
         [0067]    In the cases other than the above, the following equation is satisfied. 
         [0000]      ε1 [n]= 0 ×ε[n]   (16) 
         [0000]      Similarly, 
         [0000]      ε2 [n]= 1 ×ε[n ] when  PP 1 [f]+N/ 4 ≦P[n]≦PP 2 [f]+N/ 4 is satisfied, and 
         [0000]      ε2 [n]=− 1 ×ε[n ] when  PP 1 [f]+N× ¾≦ P[n]≦PP 2 [f]+N× ¾ is satisfied. 
         [0068]    In the cases other than the above, the following equation is satisfied. 
         [0000]      ε1 [n]= 0 ×ε[n]   (17) 
         [0069]    In other words, ε 1 [ n ] and ε 2 [ n ] have frequency f of noise to be controlled that is the same as that of ε[n], which are equivalent to the product of square waves having one amplitude around 0. When the square wave signal at the side of ε 1 [ n ] is represented by H 1 [ n ] and the square wave signal at the side of ε 2 [ n ] is represented by H 2 [ n ], they can be represented as follows. 
         [0000]      ε1 [n]=H 1 [n]×ε[n]   (18) 
         [0000]      ε2 [n]=H 2 [n]×ε[n]   (19) 
         [0070]    Herein, formulae (16) and (17) show that H 1 [ n ] and H 2 [ n ] are different from each other by ¼ cycle. 
         [0071]      FIG. 6A  is a graph showing a time-base waveform of square wave signal H 1 [ n ] (H 2 [ n ]) formed for processing an error signal in the active noise controller in accordance with the first exemplary embodiment of the present invention.  FIG. 6B  is a graph showing a harmonic analysis of square wave signal H 1 [ n ] (H 2 [ n ]) formed for processing an error signal in the active noise controller in accordance with the first exemplary embodiment of the present invention.  FIGS. 6A and 6B  show that each of square wave signals H 1 [ n ] and H 2 [ n ] includes a fundamental frequency component and odd-order higher harmonic wave. These are generally represented by the following formulae. 
         [0000]        H 1 [n]=A 1 Sin(2 πfn/T )+ A 2 Sin(2π×3 fn/T )+ A 3 Sin(2π×5 fn/T )  (20) 
         [0000]        H 2 [n]=A 1 Cos(2 πfn/T )+ A 2 Cos(2π×3 fn/T )+ A 3 Cos(2π×5 fn/T )  (21) 
         [0072]    On the other hand, when coefficient updating formulae (16) and (17) of the digital filter are modified and the relations represented by the formulae (20) and (21) are substituted, the following equations are obtained: 
         [0000]      Δ W 1 =W 1 [n+ 1 ]−W 1 [n]=−ε[n]×H 1 [n]   
         [0000]      Δ W 2 =W 2 [n+ 1 ]−W 2 [n]=−ε[n]×H 2 [n]   
         [0000]        W 1 =ΣΔW 1=Σ(−ε[ n]×H 1 [n] )  (22) 
         [0000]        W 2 =ΣΔW 2=Σ(−ε[ n]×H 2 [n] )  (23) 
         [0000]    wherein W 1  and W 2  are proportional to the accumulated values of (−ε[ n ]×H 1 [ n ]) and (−ε[n]×H 2 [ n ]), respectively. 
         [0073]    Herein, when ε[n] is sine wave Sin(2πfn/T) with frequency f, W 1  is represented as follows from formulae (20) and (22). 
         [0000]    
       
         
           
             
               
                 
                   
                     W 
                      
                     
                         
                     
                      
                     1 
                   
                   = 
                     
                    
                   
                     Σ 
                      
                     
                       ( 
                       
                         
                           - 
                           
                             ɛ 
                              
                             
                               [ 
                               n 
                               ] 
                             
                           
                         
                         × 
                         H 
                          
                         
                             
                         
                          
                         
                           1 
                            
                           
                             [ 
                             
                               n 
                                
                               
                                   
                               
                                
                               1 
                             
                             ] 
                           
                         
                       
                       ) 
                     
                   
                 
               
             
             
               
                 
                   = 
                     
                    
                   
                     Σ 
                      
                     
                       { 
                       
                         
                           - 
                           
                             Sin 
                              
                             
                               ( 
                               
                                 2 
                                  
                                 π 
                                  
                                 
                                     
                                 
                                  
                                 
                                   fn 
                                   / 
                                   T 
                                 
                               
                               ) 
                             
                           
                         
                         × 
                         
                           ( 
                           
                             
                               
                                 
                                   
                                     A 
                                      
                                     
                                         
                                     
                                      
                                     1 
                                      
                                     
                                       Sin 
                                        
                                       
                                         ( 
                                         
                                           2 
                                            
                                           π 
                                            
                                           
                                               
                                           
                                            
                                           
                                             fn 
                                             / 
                                             T 
                                           
                                         
                                         ) 
                                       
                                     
                                   
                                   + 
                                 
                               
                             
                             
                               
                                 
                                   
                                     A 
                                      
                                     
                                         
                                     
                                      
                                     2 
                                      
                                     
                                       Sin 
                                        
                                       
                                         ( 
                                         
                                           2 
                                            
                                           π 
                                           × 
                                           3 
                                            
                                           
                                             fn 
                                             / 
                                             T 
                                           
                                         
                                         ) 
                                       
                                     
                                   
                                   + 
                                 
                               
                             
                             
                               
                                 
                                   
                                     A 
                                      
                                     
                                         
                                     
                                      
                                     3 
                                      
                                     
                                       Sin 
                                        
                                       
                                         ( 
                                         
                                           2 
                                            
                                           π 
                                           × 
                                           5 
                                            
                                           
                                             fn 
                                             / 
                                             T 
                                           
                                         
                                         ) 
                                       
                                     
                                   
                                   + 
                                   … 
                                 
                               
                             
                           
                           ) 
                         
                       
                       } 
                     
                   
                 
               
             
           
         
       
     
         [0074]    Furthermore, since the accumulated value of components having different frequencies becomes 0 due to the orthogonal property of sine waves, the following formula is satisfied. 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       
                         
                           W 
                            
                           
                               
                           
                            
                           1 
                         
                         = 
                           
                          
                         
                           Σ 
                            
                           
                             ( 
                             
                               
                                 - 
                                 
                                   ɛ 
                                    
                                   
                                     [ 
                                     n 
                                     ] 
                                   
                                 
                               
                               × 
                               H 
                                
                               
                                   
                               
                                
                               
                                 1 
                                  
                                 
                                   [ 
                                   n 
                                   ] 
                                 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                   
                     
                       
                         = 
                           
                          
                         
                           Σ 
                            
                           
                             ( 
                             
                               
                                 - 
                                 
                                   Sin 
                                    
                                   
                                     ( 
                                     
                                       2 
                                        
                                       π 
                                        
                                       
                                           
                                       
                                        
                                       
                                         fn 
                                         / 
                                         T 
                                       
                                     
                                     ) 
                                   
                                 
                               
                               × 
                               A 
                                
                               
                                   
                               
                                
                               1 
                                
                               
                                   
                               
                                
                               
                                 
                                   Sin 
                                    
                                   
                                     ( 
                                     
                                       2 
                                        
                                       π 
                                        
                                       
                                           
                                       
                                        
                                       
                                         fn 
                                         / 
                                         T 
                                       
                                     
                                     ) 
                                   
                                 
                                  
                                 
                                   [ 
                                   n 
                                   ] 
                                 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   24 
                   ) 
                 
               
             
           
         
       
     
         [0075]    The same is true to W 2 . That is to say, in both W 1  and W 2 , the product of only components with frequency f is accumulated, which is equal to the conventional ones using sine waves for reference signal. Coefficients W 1  and W 2  converge so that the noise with frequency f is reduced. Thus, similar to the conventional example using the sine wave for reference signal, the present invention can also reduce noise having the intended frequency f. 
         [0076]    Furthermore, the size of the difference between the upper limit point and the lower limit point (PP 2 [ f ]−PP 1 [ f ]) in the present invention can be selected arbitrarily. In other words, actually, they can be treated the same as μ (convergence coefficient) in a conventional example. That is to say, the larger PP 2 [ f ]−PP 1 [ f ] is, the higher the convergence rate is. Meanwhile, when PP 2 [ f ]−PP 1 [ f ] is small, the convergence rate becomes slow. Thus, the convergence rate can be adjusted by the size of the PP 2 [ f ]−PP 1 [ f].    
         [0077]    Herein, the method of the present invention and the method described in patent document 1 are compared with each other in terms of an operation load. The method described in patent document 1 uses basic sine wave signal correction value table  19  (a basic sine wave signal correction value in a frequency of f (Hz) is referred to as C 1 [ f ]) and reference cosine wave signal correction value table  20  (a basic cosine wave signal correction value in a frequency of f (Hz) is referred to as C 2 [ f ]) simulating a transfer characteristic value from loudspeaker  10  to microphone  11  so as to generate reference sine wave signal r 1 [ n ] and reference cosine wave signal r 2 [ n ] by formulae (25) and (26), respectively. 
         [0000]        r 1 [n]=C 1 [f]×x 1 [n]+C 2 [f]×x 2 [n]   (25) 
         [0000]        r 2 [n]=C 1 [f]×x 2 [n]−C 2 [f]×x 1 [n]   (26) 
         [0078]    Firstly, formulae (25) and (26) include two multiplications respectively while the present invention does not need multiplication because it does not use a reference signal. Furthermore, also in the update of coefficients, the method described in patent document 1 needs two multiplications as shown in formulae (1) and (2). Meanwhile, the present invention does not need multiplication as is apparent by formulae (14) and (15). 
         [0079]    Thus, the method described in patent document 1 needs four multiplications in every sampling cycle in order to obtain coefficients W 1  and W 2 , respectively. However, the present invention does not need multiplication. Therefore, according to the active noise controller in accordance with the present invention, an operation load can be reduced as compared with the method described in patent document 1. 
         [0080]    Furthermore, in the present invention, by preparing a plurality of first and second one-tap digital filters  7  and  8  and a plurality of first and second adaptive control algorithm operating units  12  and  13 , respectively, multiple-order components of noise to be controlled can be suppressed. 
       INDUSTRIAL APPLICABILITY 
       [0081]    An active noise controller in accordance with the present invention can reduce an operation load by minimizing the execution of product-sum operation, and thus it is useful as a practical and low-cost active noise controller.