EP0759606A2 - Active adaptive selective control system - Google Patents
Active adaptive selective control system Download PDFInfo
- Publication number
- EP0759606A2 EP0759606A2 EP96306121A EP96306121A EP0759606A2 EP 0759606 A2 EP0759606 A2 EP 0759606A2 EP 96306121 A EP96306121 A EP 96306121A EP 96306121 A EP96306121 A EP 96306121A EP 0759606 A2 EP0759606 A2 EP 0759606A2
- Authority
- EP
- European Patent Office
- Prior art keywords
- signal
- model
- output
- input
- error
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
- G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
- G10K11/1787—General system configurations
- G10K11/17879—General system configurations using both a reference signal and an error signal
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
- G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
- G10K11/1781—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase characterised by the analysis of input or output signals, e.g. frequency range, modes, transfer functions
- G10K11/17813—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase characterised by the analysis of input or output signals, e.g. frequency range, modes, transfer functions characterised by the analysis of the acoustic paths, e.g. estimating, calibrating or testing of transfer functions or cross-terms
- G10K11/17817—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase characterised by the analysis of input or output signals, e.g. frequency range, modes, transfer functions characterised by the analysis of the acoustic paths, e.g. estimating, calibrating or testing of transfer functions or cross-terms between the output signals and the error signals, i.e. secondary path
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
- G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
- G10K11/1783—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase handling or detecting of non-standard events or conditions, e.g. changing operating modes under specific operating conditions
- G10K11/17833—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase handling or detecting of non-standard events or conditions, e.g. changing operating modes under specific operating conditions by using a self-diagnostic function or a malfunction prevention function, e.g. detecting abnormal output levels
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
- G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
- G10K11/1785—Methods, e.g. algorithms; Devices
- G10K11/17853—Methods, e.g. algorithms; Devices of the filter
- G10K11/17854—Methods, e.g. algorithms; Devices of the filter the filter being an adaptive filter
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
- G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
- G10K11/1787—General system configurations
- G10K11/17879—General system configurations using both a reference signal and an error signal
- G10K11/17881—General system configurations using both a reference signal and an error signal the reference signal being an acoustic signal, e.g. recorded with a microphone
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/10—Applications
- G10K2210/112—Ducts
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/30—Means
- G10K2210/301—Computational
- G10K2210/3012—Algorithms
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/30—Means
- G10K2210/301—Computational
- G10K2210/3017—Copy, i.e. whereby an estimated transfer function in one functional block is copied to another block
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/30—Means
- G10K2210/301—Computational
- G10K2210/3028—Filtering, e.g. Kalman filters or special analogue or digital filters
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/30—Means
- G10K2210/301—Computational
- G10K2210/3029—Fuzzy logic; Genetic algorithms
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/30—Means
- G10K2210/301—Computational
- G10K2210/3037—Monitoring various blocks in the flow chart
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
- G10K2210/30—Means
- G10K2210/301—Computational
- G10K2210/3049—Random noise used, e.g. in model identification
Definitions
- the invention relates to active adaptive control systems, and more particularly to improvements for selectively controlling performance of the active adaptive model.
- Active acoustic attenuation involves injecting a canceling acoustic wave to destructively interfere with and cancel an input acoustic wave.
- the output acoustic wave is sensed with an error transducer, such as a microphone or an accelerometer, which supplies an error signal to an adaptive filter control model which in turn supplies a correction signal to a canceling output transducer or actuator, such as a loudspeaker or a shaker, which injects an acoustic wave to destructively interfere with the input acoustic wave and cancel or reduce same such that the output acoustic wave at the error transducer is zero or some other desired value.
- an error transducer such as a microphone or an accelerometer
- an adaptive filter control model which in turn supplies a correction signal to a canceling output transducer or actuator, such as a loudspeaker or a shaker, which injects an acoustic wave to destructively interfere with the input acoustic wave and cancel or reduce same such that the
- An active adaptive control system minimizes an error signal by introducing a control signal from an output transducer to combine with the system input signal and yield a system output signal.
- the system output signal is sensed with an error transducer providing the error signal.
- An adaptive filter model has a model input from a reference signal correlated with the system input signal, an error input from the error signal, and outputs a correction signal to the output transducer to introduce a control signal matching the system input signal, to minimize the error signal.
- the filter coefficients are updated according to a weight update signal which is the product of the reference signal and the error signal.
- the present invention is applicable to active adaptive control systems, including active acoustic attenuation systems.
- the invention maximizes model performance by concentrating model adaptation in frequency ranges of interest, and protects the output transducer against overdriving of same.
- Performance of the model is spectrally controlled to maximize the correction signal sent to the output transducer such that at frequencies where maximum power from the output transducer reaches the error transducer, the correction signal supplied to the output transducer is maximized, and at frequencies where minimum power from the output transducer reaches the error transducer, the correction signal supplied to the output transducer is minimized.
- model adaptation is selectively controlled to provide desired performance.
- Model performance is controlled by fuzzy logic to provide self-designing control architecture using fuzzy rules and/or to control a filter transfer function and/or to control filter weights used in an update process for feedforward and/or feedback, including FIR (finite impulse response) and IIR (infinite impulse response) applications and/or to control magnitude and/or rate of change of a leak signal degrading performance of the model.
- fuzzy logic to provide self-designing control architecture using fuzzy rules and/or to control a filter transfer function and/or to control filter weights used in an update process for feedforward and/or feedback, including FIR (finite impulse response) and IIR (infinite impulse response) applications and/or to control magnitude and/or rate of change of a leak signal degrading performance of the model.
- Fig. 1 is a schematic illustration of an active adaptive control system known in the prior art.
- Fig. 2 is a schematic illustration of an active adaptive control system in accordance with co-pending U.S. Application S.N. 08/264,510, filed June 23, 1994.
- Fig. 3 is a graph showing performance of the system of Fig. 2.
- Fig. 4 is a graph further showing performance of the system of Fig. 2.
- Fig. 5 is a graph showing an alternate performance of the system of Fig. 2.
- Fig. 6 is a graph further showing alternate performance of the system of Fig. 2.
- Fig. 7 is a schematic illustration of an active adaptive control system in accordance with the present invention.
- Fig. 8 is a graph showing frequency versus output and illustrates performance of the system of Fig. 7.
- Fig. 9 is a schematic illustration of an active adaptive control system.
- Fig. 10 is a schematic illustration of an adaptive filter model, and illustrates a principle employed by a system in accordance with copending U.S. Application Serial No. 08/166,698, filed December 14, 1993.
- Fig. 11 is like Fig. 10 and shows another manner of implementing the principle thereof.
- Fig. 12 is a schematic illustration of an active adaptive control system in accordance with the system of the '698 application.
- Fig. 13 shows a further embodiment of the system of Fig. 12.
- Fig. 14 shows a further embodiment of the system of Fig. 12.
- Fig. 15 is a schematic illustration of an active adaptive control system in accordance with the present invention.
- Fig. 16 shows a further embodiment of the system of Fig. 15.
- Fig. 17 shows a further embodiment of the system of Fig. 15.
- Fig. 18 is a schematic illustration of an active adaptive control system in accordance with the present invention.
- Fig. 19 shows a further embodiment of the system of Fig. 18.
- Fig. 20 shows a further embodiment of the system of Fig. 18.
- Fig. 21 shows a further embodiment of the system of Fig. 18.
- Figs. 22-24 illustrate a further embodiment of the invention.
- Fig. 1 shows an active adaptive control system similar to that shown in U.S. Patent 4,677,676, incorporated herein by reference, and uses like reference numerals therefrom where appropriate to facilitate understanding.
- the system introduces a control signal from a secondary source or output transducer 14, such as a loudspeaker, shaker, or other actuator or controller, to combine with the system input signal 6 and yield a system output signal 8.
- An input transducer 10 such as a microphone, accelerometer, or other sensor, senses the system input signal and provides a reference signal 42.
- An error transducer 16 such as a microphone, accelerometer, or other sensor, senses the system output signal and provides an error signal 44.
- Adaptive filter model 40 adaptively models the system and has a model input from reference signal 42 correlated to system input signal 6, and an output outputting a correction signal 46 to output transducer 14 to introduce the control signal according to a weight update signal 74.
- Reference signal 42 and error signal 44 are combined at multiplier 72 to provide the weight update signal through delay element 73.
- the reference signal 42 may be provided by one or more error signals, in the case of a periodic system input signal, "Active Adaptive Sound Control In A Duct: A Computer Simulation", J.C. Burgess, Journal of Acoustic Society of America, 70(3), September 1981, pages 715-726, U.S. Patents 5,206,911, 5,216,722, incorporated herein by reference.
- one or more previous weights are added to the current product of reference signal 42 and error signal 44 at summer 75.
- Leakage factor ⁇ at 77 multiplies one or more previous weights, after passage through one or more delay elements 73, by an exponential decay factor less than one before adding same at summer 75 to the current product of reference signal 42 and error signal 44, Adaptive Signal Processing , Widrow and Stearns, Prentice-Hall, Inc., Engelwood Cliffs, NJ, 1985, pages 376-378, including equations 13.27 and 13.31.
- a deficiency of this method is that it reduces control effort and degrades performance across all power levels, regardless of whether such reduced effort is desired.
- leakage of the weight update signal is provided in response to a given condition of a given parameter, to control performance of the model on an as needed basis.
- leakage is varied as a function of correction signal 46.
- a variable leakage factor ⁇ is provided at 79 in Fig. 2, replacing fixed ⁇ 77 of Fig. 1.
- Leakage factor ⁇ at 79 is varied from a maximum value of 1.0 affording maximum control effort, to a minimum value such as zero providing minimum control effort.
- Leakage is varied as a function of the output power of correction signal 46 supplied from the output of model 40 to output transducer 14.
- the leakage is varied as a discontinuous step function of the output power of the correction signal.
- ⁇ is abruptly, nonlinearly changed as a step function from a first level 83 to a second level 85.
- the reduction at 85 reduces the weight update signal summed at summer 75 with the product of the reference signal 42 and error signal 44 from multiplier 72, and hence reduces the weight update signal supplied to model 40.
- the noted reduction of ⁇ at threshold 81 increases leakage of the weight update signal, Fig. 4, from level 87 to level 89.
- leakage is varied as a continuous function of the output power of the correction signal.
- ⁇ is maintained at level 83 until output power reaches threshold 81, and then is linearly decreased as shown at 91 as a continuous linearly changing value as a function of increasing output power above threshold 81.
- leakage is maintained at level 87 until output power reaches threshold 81, and then is linearly increased at 93 as a continuous linearly changing value as a function of increasing output power above threshold 81.
- leakage is adjustably varied to vary performance of the model by multiplying a previous weight update value by variable ⁇ 79 and adding the result at summer 75 to the product of reference signal 42 and error signal 44 from multiplier 72.
- ⁇ 79 is varied as a function of correction signal 46, preferably the output power of such correction signal.
- Fig. 7 illustrates the present invention and uses like reference numerals from Fig. 2 and from Figs. 19 and 20 of the incorporated '676 patent.
- the transfer function from output transducer 14 to error transducer 16 is modeled with an adaptive filter model C at 142, as in the incorporated '676 patent.
- Auxiliary random signal source 140 introduces a random signal into the output of model 40 at summer 152 and into the C model at 148.
- the auxiliary random signal from source 140 is random and uncorrelated with the system input signal 6 and in preferred form is provided by a Galois sequence, M.R. Schroeder, "Number Theory In Science And Communications", Berlin, Springer-Berlag, 1984, pages 252-261, though other random uncorrelated signal sources may be used.
- the Galois sequence is a pseudo random sequence that repeats after 2 M -1 points, where M is the number of stages in a shift register.
- the Galois sequence is preferred because it is easy to calculate and can easily have a period much longer than the response time of the system.
- the input 148 to C model 142 is multiplied with the error signal from error transducer 16 at multiplier 68, and the resultant product provided as weight update signal 67.
- Model 142 models the transfer function from output transducer 14 to error transducer 16, including the transfer function of each.
- the transfer function from output transducer 14 to error transducer 16 may be modeled without a random signal source, as in U.S. Patent 4,987,598, incorporated herein by reference.
- Auxiliary source 140 introduces an auxiliary random signal such that error transducer 16 also senses the auxiliary signal from the auxiliary source.
- the auxiliary signal may be introduced into the recursive loop of the A and B filters as in Fig. 19 of the incorporated '676 patent at summer 152, or alternatively the auxiliary signal may be introduced into the model after the recursive loop, i.e. introducing the auxiliary signal only to line 46, and not to line 47.
- a copy of model 142 is provided at 145 to compensate the noted transfer function, as in the incorporated '676 patent.
- the present system and method involves introducing a control signal from output transducer 14 to combine with system input signal 6 and yield system output signal 8, sensing the system output signal with error transducer 16 and providing an error signal 44, providing adaptive filter model 40 having a model input from reference signal 42 correlated to system input signal 6, and an output outputting a correction signal 46 to output transducer 14 to introduce the control signal according to weight update signal 74.
- performance of model 40 is spectrally controlled to maximize the signal sent to output transducer 14 at frequencies of interest or where it can be maximally effective, and minimize the signal sent to output transducer 14 at frequencies of non interest or where it is only minimally effective or is ineffective.
- a spectral leak signal is provided which degrades performance of the model.
- the leak signal is controlled according to frequency.
- the correction signal from the output of model 40 is monitored, and the leak signal is controlled in response thereto.
- the correction signal is filtered by filter 95 to provide the leak signal.
- the correction signal from the output of model 40 is spectrally processed by supplying the correction signal through the frequency responsive spectral transfer function provided by filter 95 to provide the leak signal to the error input of model 40.
- the correction signal is spectrally processed such that at frequencies where maximum power from output transducer 14 reaches error transducer 16, the correction signal is maximized, and at frequencies where minimum power from output transducer 14 reaches error transducer 16, the correction signal is minimized.
- correction signal 46 is spectrally processed by a function of model 142.
- such function is the inverse of C model 142 as provided at inverse C model, C -1 , at 95.
- the output of inverse C model 95 is supplied to an optional peak detector as provided by summer 97 comparing the output of inverse C model 95 with a desired peak value 99. When the output of model 95 rises above level 99, the positive output of summer 97 controls variable leakage factor ⁇ at 79, as above.
- the inverse C model includes the inverse of the transfer function of output transducer 14, inverse S, S -1 , and/or the inverse of the transfer function of the error path, inverse E, E -1 , between output transducer 14 and error transducer 16.
- control logic may be used to respond to the output of inverse C model 95 and control leakage factor ⁇ at 79 according to designated conditions or rules to generate or compute a leak value or control the leaking process, to be further described hereinafter.
- filter 95 is a displacement function of output transducer 14 such as a loudspeaker, to protect the latter against overdriving.
- an RMS (root mean square) to DC (direct current) conversion function is provided between filter 95 and peak detector 97 for average level detection to control convergence rate of the leak process.
- the transfer function of filter 95 may be linear or nonlinear.
- C model 142 has a maximum transfer characteristic
- the inverse of the C model at 95 has a minimum transfer characteristic, as shown at 105 and 107.
- the minimum transfer characteristic at 105 and 107 minimizes leakage of the update signal to model 40, to enable maximum output of model 40.
- C model 142 has a minimum transfer characteristic
- the inverse of the C model at 95 has a maximum transfer characteristic, as shown at 111.
- the maximum transfer characteristic at 111 maximizes leakage of the update signal, to minimize the output of model 40.
- Inverse C model 95 spectrally senses correction signal 46 and provides selective leakage of weight update signal 74 in response thereto, to control performance of model 40 according to frequency, to optimize performance of model 40 in frequency ranges such as 101 and 103 where the model can effectively control or cancel input signal 6. Outside of such ranges, e.g. at 218, inverse C model 95 minimizes performance of model 40, to avoid using computation or adaptation power where inefficient or unneeded, and to prevent the model from continually trying to generate an output in regions where it is ineffective to attempt to control or cancel input signal 6. Model 95 affords frequency weighting of the weight update signal.
- Fig. 9 is similar to Fig. 5 of U.S. Patent 4,677,676, incorporated herein by reference, and uses like reference numerals to facilitate understanding.
- the system introduces a control signal from an output transducer 14, such as a loudspeaker, shaker, or other actuator or controller, to combine with the system input signal 6 and yield a system output signal 8.
- An input transducer 10 such as a microphone, accelerometer, or other sensor, senses the system input signal.
- An error transducer 16 such as a microphone, accelerometer, or other sensor, senses the system output signal and provides an error signal.
- Adaptive filter model 40 adaptively models the system and has a model input 42 from input transducer 10, an error input 44 from error transducer 16, and a model output 46 outputting a correction signal to output transducer 14 to introduce the control signal.
- the input signal at 42 may be provided by one or more error signals, in the case of a periodic system input signal, "Active Adaptive Sound Control In A Duct: A Computer Simulation", J.C. Burgess, Journal of Acoustic Society of America, 70(3), September 1981, pages 715-726, U.S. Patents 5,206,911, 5,216,722, incorporated herein by reference.
- the system of the '698 application provides an active adaptive control system wherein the performance of model 40 is intentionally and selectively constrained by driving the output 46 of the model towards zero in response to a given condition of a given parameter.
- an active noise control system it may be desirable to cancel noise only in a given frequency band, and leave the noise uncanceled for frequencies outside the band.
- it may be desirable to selectively control the system output signal by selectively controlling introduction of the control signal from output transducer 14 to match or not match the system input signal.
- Fig. 10 One way of accomplishing this is shown in Fig. 10 wherein the output of model 40 is supplied to its error input, such that when the model adapts to drive its error input towards zero, the output of the model is necessarily also driven towards zero.
- Fig. 11 shows another manner of implementing this principle wherein a copy of the model is provided at 200, and the output of model copy 200 supplies the error signal to the error input of model 40.
- model 40 adapts to drive its error input towards zero, which in turn requires that the output of copy 200 be driven towards zero, which in turn means that the output of model 40 is driven towards zero because M copy 200 is a duplicate of model 40.
- Model 40 normally adapts to a converged condition wherein its output at 46 provides a correction signal to output transducer 14 which outputs a control signal matching the system input signal or a designated relative value correlated thereto. For example, in a noise cancellation system, the matching control signal from output transducer 14 cancels the input noise.
- the output of model 40 in response to a given condition of a given parameter, the output of model 40 is driven towards zero by driving the output of the model towards its error input, such that when the model adapts to drive the error signal towards zero, the output of the model is also driven towards zero.
- model 40 This is accomplished by providing a copy 200 of model 40, and supplying the output of the copy to an error input 202 of the model which is summed at summer 204 with the error signal at error input 44 from error transducer 16.
- the model adapts to drive the error input towards zero which in turn requires that the output of copy 200 and hence the output of model 40 are driven towards zero, to provide the noted constrained performance.
- the driving of model output 46 towards zero provides a zero or at least a reduced correction signal to output transducer 14 to constrain or reduce modification and/or cancellation of the system input signal 6.
- a random signal is provided at 206 from an auxiliary random signal source 208, preferably provided by a Galois sequence, M.R. Schroeder, "Number Theory In Science And Communications", Berlin, Springer-Berlag, 1984, pages 252-261, though other random signal sources may be used, uncorrelated with the system input signal 6.
- the Galois sequence is a pseudo random sequence that repeats after 2 M -1 points, where M is the number of stages in a shift register. The Galois sequence is preferred because it is easy to calculate and can easily have a period much longer than the response time of the system.
- the random signal is supplied through a stopband filter 210 to model copy 200 at 212. Stopband filter 210 blocks frequencies in the stopband, and passes frequencies outside the stopband.
- This arrangement provides a spectral leak signal at 202 in response to a given condition of a given parameter, for example a frequency outside the stopband.
- a given condition of a given parameter for example a frequency outside the stopband.
- the noted given parameter is frequency
- the given condition is a designated sub-optimum performance band outside the stopband.
- the spectral leak error signal at 202 drives the correction signal at model output 46 towards zero and provides sub-optimum performance of model 40. Outside of the sub-optimum performance band, i.e. within the stopband of filter 210, there is no signal at 212 and hence the output of copy 200 is undefined, and the error signal from error transducer 16 at error input 44 is maximally effective and model 40 optimally responds thereto and drives the correction signal at output 46 toward a value matching the system input signal 6.
- the spectral leak signal is present at error input 202, it constrains performance of model 40 by driving or at least attempting to drive the correction signal at model output 46 towards zero.
- the relative influence or amplitudes of the error signals at error inputs 44 and 202 are adjusted to provide the desired relative dominance.
- the noted relative amplitudes are set such that the error signal at error input 202 dominates the error signal at error input 44, and hence the correction signal at model output 46 is driven towards zero and away from a value matching the system input signal 6.
- the method of the '698 application involves driving error input 44 to drive the correction signal at model output 46 toward a value matching the system input signal, and selectively driving error input 202 to drive the correction signal at model output 46 away from the matching value by driving the correction signal towards zero.
- the arrangement provides a spectral leak signal to error input 202 in response to the noted given condition of the given parameter, e.g. a frequency outside the stopband, such that in the presence of the given condition, the spectral leak signal drives the correction signal at model output 46 towards zero, and in the absence of the given condition the error signal at error input 44 drives the correction signal at model output 46 towards a value matching the system input signal 6.
- Stopband filter 210 blocks frequencies in a given stopband at which modification or cancellation of the system input signal 6 by model 40 is desired. Filter 210 passes frequencies in a given passband at which modification or cancellation of the system input signal by model 40 is undesired or not possible.
- the control signal output by output transducer 14 is driven toward a value matching the system input signal 6 only for frequencies in the stopband. At frequencies in the stopband, the error signal at error input 44 is dominant, and the control signal output by output transducer 14 is driven toward a value matching the system input signal 6. At frequencies in the passband, the error signal at error input 202 is dominant, and the control signal output by output transducer 14 is driven away from a value matching the system input signal 6.
- model M at 40 is preferably an adaptive recursive filter having a transfer function with both poles and zeros.
- Model M is provided by an IIR, infinite impulse response, filter, e.g. a recursive least mean square, RLMS, filter having a first algorithm filter provided by an FIR, finite impulse response, filter, e.g. a least mean square, LMS, filter A at 12, and a second algorithm filter provided by an FIR filter, e.g. an LMS filter, B at 22.
- Filter A provides a direct transfer function
- filter B provides a recursive transfer function.
- the transfer function from output transducer 14 to error transducer 16 is modeled by a filter, e.g. an LMS or RLMS filter, C at 142, as in the incorporated '676 patent.
- Auxiliary random signal source 140 introduces a random signal into the output of model 40 at summer 152 and into the C model at 148.
- the auxiliary random signal from source 140 is random and uncorrelated with the system input signal 6, and is also uncorrelated with auxiliary random signal source 208, and in preferred form is provided by a Galois sequence, M.R. Schroeder, "Number Theory In Science And Communications", Berlin, Springer-Berlag, 1984, pages 252-261, though other random uncorrelated signal sources may be used.
- the Galois sequence is a pseudo random sequence that repeats after 2 M -1 points, where M is the number of stages in a shift register. The Galois sequence is preferred because it is easy to calculate and can easily have a period much longer than the response time of the system.
- Model 142 models the transfer function from output transducer 14 to error transducer 16, including the transfer function of each.
- the transfer function from output transducer 14 to error transducer 16 may be modeled without a random signal source, as in U.S. Patent 4,987,598, incorporated herein by reference.
- Auxiliary source 140 introduces an auxiliary random signal such that error transducer 16 also senses the auxiliary signal from the auxiliary source.
- the auxiliary signal may be introduced into the recursive loop of the A and B filters as in Fig.
- auxiliary signal may be introduced into the model after the recursive loop, i.e. introducing the auxiliary signal only to line 46, and not to line 47.
- copies of model 142 are provided at 144 and 146 to compensate the noted transfer function.
- the outputs of filters A and B are summed at summer 48, whose output is summed at summer 152 with the output of random signal source 140 to provide an output resultant sum which provides the model output at 46 supplying the noted correction signal to output transducer 14.
- the output of model 142 is summed at summer 64 with the output of error transducer 16, and the resultant sum supplied as the error input to model 142 and as an error input to model 40.
- the output of error transducer 16 may be supplied directly to an error input of model 40 without being supplied through summer 64.
- M copy 200, Fig. 12 is provided by a copy of A filter 12 at A copy 214, Fig. 13, and a copy of B filter 22 at B copy 216.
- Multiplier 218 multiplies the output 220 of A copy 214 and the input 222 of A copy 214, and supplies the output resultant product at 224 to summer 226.
- Multiplier 228 multiplies the output 230 of C copy 144 and the error input at 44, and supplies the output resultant product at 232 to summer 226.
- Summer 226 sums the inputs 232 and 224, and supplies the output resultant sum as weight update signal 74 to A filter 12.
- Multiplier 234 multiplies the output 236 of B copy 216 and the input 238 of B copy 216, and supplies the output resultant product at 240 to summer 242.
- Multiplier 244 multiplies the output 246 of C copy 146 and the error input at 44, and supplies the output resultant product at 248 to summer 242.
- the summer 242 sums the inputs 248 and 240, and supplies the output resultant sum as weight update signal 78 to B filter 22.
- the input to A copy 214 and to B copy 216 is provided by the output 212 of stopband filter 210 receiving the noted random input signal at 206 from random signal source 208.
- the error signals at error inputs 232 and 224 oppositely drive the model.
- the error signal at error input 232 of the direct transfer function filter A drives the correction signal at 46 towards a value matching the system input signal 6.
- the error signal at error input 224 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero.
- this is accomplished by using a copy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero.
- the signal at error input 224 is provided only in response to a given condition of a given parameter, e.g. when the frequency is outside the stopband of filter 210.
- the relative amplitudes of the input signals at error inputs 232 and 224 are adjusted such that the signal at error input 224 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree.
- the frequency is in the stopband of filter 210, there is no output at 212, and hence no input to A copy 214 and hence the latter's output is undefined, whereby error input 232 from error signal 44 from error transducer 16 dominates and hence drives correction signal 46 to a value which matches the system input signal 6 to provide modification and/or control of the latter.
- Model 40 has a first error input provided at 232 and 248 from error transducer 16 driving the output of the model towards a value matching the system input signal 6.
- Model 40 has a second error input at 224 and 240 selectively driving the output of model 40 away from such matching value and instead driving the correction signal 46 towards zero.
- Fig. 14 is similar to Fig. 13 and uses like reference numerals where appropriate to facilitate understanding.
- Summer 260 sums the output 230 of C copy 144 and the output 212 of stopband filter 210 which supplies the input to A copy 214, and supplies the output resultant sum at 262 to multiplier 264.
- Summer 266 sums the output 268 of A copy 214 and the error input at 44, and supplies the output resultant sum at 270 to multiplier 264.
- Multiplier 264 multiplies the inputs 262 and 270, and supplies the output resultant product as weight update signal 74 to A filter 12.
- Summer 272 sums the output 246 of C copy 146 and the output 212 of stopband filter 210 which supplies the input to B copy 216, and supplies the output resultant sum at 274 to multiplier 276.
- Summer 278 sums the output 280 of B copy 216 and the error input at 44, and supplies the output resultant sum at 282 to multiplier 276.
- Multiplier 276 multiplies inputs 282 and 274, and supplies the output resultant product as weight update signal 78 to B filter 22.
- the error signals at 44 and 268 oppositely drive the model.
- the error signal at error input 44 of the direct transfer function filter A drives the correction signal at 46 towards a value matching the system input signal 6.
- the error signal at error input 268 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero.
- this is accomplished by using a copy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero.
- the signal at error input 268 is provided only in response to a given condition of a given parameter, e.g., when the frequency is outside the stopband of filter 210.
- the relative amplitudes of the input signals at error inputs 44 and 268 to summer 266 are adjusted such that the signal at error input 268 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree.
- the frequency is in the stopband of filter 210, there is no output at 212, and hence no input to A copy 214 and hence the latter's output is undefined, whereby error input 44 to summer 266 from error transducer 16 dominates and hence the model drives correction signal 46 to a value which matches the system input signal 6 to provide modification and/or control of the latter.
- Model 40 in Fig. 14 has a first error input provided at 44 to summers 266 and 278 from error transducer 16 driving the output of the model towards a value matching the system input signal 6.
- Model 40 has a second error input at 268 and 280 selectively driving the output of model 40 away from such matching value and instead driving the correction signal 46 towards zero.
- Fig. 15 illustrates the present invention and uses like reference numerals from above and from Figs. 19 and 20 of the incorporated '676 patent.
- Random noise signal 206 from random noise source 208 is supplied through a first inverse C model copy 302 whose output is supplied to the filter input 304 of M copy 200.
- Correction signal 46 from the output of model 40 through summer 152 is supplied through a second inverse C model copy 306 to the convergence rate gain control input 308 of M copy 200.
- the spectral transfer function provided by filter 302 spectrally controls performance of model 40 to maximize the correction signal at frequencies of interest.
- the spectral transfer function is the inverse C model copy having the characteristics shown in Fig. 8, though other frequency responsive spectral transfer functions may be used.
- the spectral leak signal provided through transfer functions 302 and 200 degrades performance of model 40, as above described in conjunction with Fig. 12.
- Transfer function 302 controls the leak signal according to frequency. It is preferred that the leak signal also be controlled in response to the correction signal from the output of model 40.
- Transfer function 306 spectrally monitors the correction signal and filters same to provide a spectral control signal controlling the gain through M copy 200. Transfer function 306 and the gain control provided at control input 308 are optional.
- the spectral leak control provided through transfer function 302 to filter input 304 may be used with or without transfer function 306 and control input 308.
- Filter input 304 of M copy 200 receives a spectrally processed input signal from spectral transfer function 302.
- Control input 308 of M copy 200 controls gain of M copy 200 and receives a spectral control signal from transfer function 306.
- the spectral control signal is responsive to correction signal 46.
- Spectral transfer functions 302 and 306 may be the same, as shown, or may be different.
- model M at 40 is preferably an adaptive recursive filter having a transfer function with both poles and zeros.
- Model M is provided by an IIR, infinite impulse response, filter, e.g. a recursive least mean square, RLMS, filter having a first algorithm filter provided by an FIR, finite impulse response, filter, e.g. a least mean square, LMS, filter A at 12, and a second algorithm filter provided by an FIR filter, e.g. an LMS filter, B at 22.
- Filter A provides a direct transfer function
- filter B provides a recursive transfer function.
- the transfer function from output transducer 14 to error transducer 16 is modeled by a filter, e.g. an LMS or RLMS filter, C at 142, as in the incorporated '676 patent.
- Auxiliary random signal source 140 introduces a random signal into the output of model 40 at summer 152 and into the C model at 148.
- the auxiliary random signal from source 140 is random and uncorrelated with the system input signal 6, and is also uncorrelated with auxiliary random signal source 208, and in preferred form is provided by a Galois sequence, M.R. Schroeder, "Number Theory In Science And Communications", Berlin, Springer-Berlag, 1984, pages 252-261, though other random uncorrelated signal sources may be used.
- the Galois sequence is a pseudo random sequence that repeats after 2 M -1 points, where M is the number of stages in a shift register. The Galois sequence is preferred because it is easy to calculate and can easily have a period much longer than the response time of the system.
- Model 142 models the transfer function from output transducer 14 to error transducer 16, including the transfer function of each.
- the transfer function from output transducer 14 to error transducer 16 may be modeled without a random signal source, as in U.S. Patent 4,987,598, incorporated herein by reference.
- Auxiliary source 140 introduces an auxiliary random signal such that error transducer 16 also senses the auxiliary signal from the auxiliary source.
- the auxiliary signal may be introduced into the recursive loop of the A and B filters as in Fig.
- auxiliary signal may be introduced into the model after the recursive loop, i.e. introducing the auxiliary signal only to line 46, and not to line 47.
- copies of model 142 are provided at 144 and 146 to compensate the noted transfer function.
- the outputs of filters A and B are summed at summer 48, whose output is summed at summer 152 with the output of random signal source 140 to provide an output resultant sum which provides the model output at 46 supplying the noted correction signal to output transducer 14.
- the output of model 142 is summed at summer 64 with the output of error transducer 16, and the resultant sum supplied as the error input to model 142 and as an error input to model 40.
- the output of error transducer 16 may be supplied directly to an error input of model 40 without being supplied through summer 64.
- M copy 200, Fig. 15, is provided by a copy of A filter 12 at A copy 214, Fig. 16, and a copy of B filter 22 at B copy 216.
- Multiplier 218 multiplies the output 220 of A copy 214 and the input 222 of A copy 214, and supplies the output resultant product at 224 to summer 226.
- Multiplier 228 multiplies the output 230 of C copy 144 and the error input at 44, and supplies the output resultant product at 232 to summer 226.
- Summer 226 sums the inputs 232 and 224, and supplies the output resultant sum as weight update signal 74 to A filter 12.
- Multiplier 234 multiplies the output 236 of B copy 216 and the input 238 of B copy 216, and supplies the output resultant product at 240 to summer 242.
- Multiplier 244 multiplies the output 246 of C copy 146 and the error input at 44, and supplies the output resultant product at 248 to summer 242.
- the summer 242 sums the inputs 248 and 240, and supplies the output resultant sum as weight update signal 78 to B filter 22.
- the input to A copy 214 and to B copy 216 is provided by the output 304 of filter 302 receiving the noted random input signal at 206 from random signal source 208.
- the error signals at error inputs 232 and 224 oppositely drive the model.
- the error signal at error input 232 of the direct transfer function filter A drives the correction signal at 46 towards a value matching the system input signal 6.
- the error signal at error input 224 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero.
- this is accomplished by using a copy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero.
- the signal at error input 224 is provided only in response to a given condition of a given parameter, e.g. when the frequency is in the ranges passed by filter 302.
- the relative amplitudes of the input signals at error inputs 232 and 224 are adjusted such that the signal at error input 224 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree.
- the frequency is outside of the ranges passed by filter 302
- there is no output at 212, and hence no input to A copy 214 and hence the latter's output is undefined, whereby error input 232 from error signal 44 from error transducer 16 dominates and hence drives correction signal 46 to a value which matches the system input signal 6 to provide modification and/or control of the latter.
- the gain is controlled by the control signal at control input 308 which is the spectrally processed correction signal supplied through the transfer function filter provided by inverse C model copy 306.
- the error inputs 248 and 240 to the recursive transfer function filter B of model 40 function comparably to error inputs 232 and 224, respectively.
- Model 40 has a first error input provided at 232 and 248 from error transducer 16 driving the output of the model towards a value matching the system input signal 6.
- Model 40 has a second error input at 224 and 240 selectively driving the output of model 40 away from such matching value and instead driving the correction signal 46 towards zero.
- Fig. 17 uses like reference numerals from above and from Figs. 19 and 20 of the incorporated '676 patent.
- Summer 260 sums the output 230 of C copy 144 and the output 212 of stopband filter 210 which supplies the input to A copy 214, and supplies the output resultant sum at 262 to multiplier 264.
- Summer 266 sums the output 268 of A copy 214 and the error input at 44, and supplies the output resultant sum at 270 to multiplier 264.
- Multiplier 264 multiplies the inputs 262 and 270, and supplies the output resultant product as weight update signal 74 to A filter 12.
- Summer 272 sums the output 246 of C copy 146 and the output 304 of filter 302 which supplies the input to B copy 216, and supplies the output resultant sum at 274 to multiplier 276.
- Summer 278 sums the output 280 of B copy 216 and the error input at 44, and supplies the output resultant sum at 282 to multiplier 276.
- Multiplier 276 multiplies inputs 282 and 274, and supplies the output resultant product as weight update signal 78 to B filter 22.
- the error signals at 44 and 268 oppositely drive the model.
- the error signal at error input 44 of the direct transfer function filter A drives the correction signal at 46 towards a value matching the system input signal 6.
- the error signal at error input 268 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero.
- this is accomplished by using a copy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero.
- the signal at error input 268 is provided only in response to a given condition of a given parameter, e.g., when the frequency is in the ranges passed by transfer function filter 302.
- the relative amplitudes of the input signals at error inputs 44 and 268 to summer 266 are adjusted such that the signal at error input 268 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree.
- Model 40 has a second error input at 268 and 280 selectively driving the output of model 40 away from such matching value and instead driving the correction signal 46 towards zero.
- Fig. 18 illustrates the present invention and uses like reference numerals from above and from Figs. 19 and 20 of the incorporated '676 patent.
- spectral transfer function filter 302 has its input 310 supplied from the correction signal from the output of model 40, rather than from random noise signal source 208, Fig. 15.
- the spectral transfer function provided by filter 306 also includes a peak detector 312 similar to peak detector 97, Fig. 7.
- the output of inverse C model 306 is supplied to peak detector 312 comparing the output of inverse C model 306 with a desired peak value 314.
- model M at 40 is preferably an adaptive recursive filter having a transfer function with both poles and zeros.
- Model M is provided by an IIR, infinite impulse response, filter, e.g. a recursive least mean square, RLMS, filter having a first algorithm filter provided by an FIR, finite impulse response, filter, e.g. a least mean square, LMS, filter A at 12, and a second algorithm filter provided by an FIR filter, e.g. an LMS filter, B at 22.
- Filter A provides a direct transfer function
- filter B provides a recursive transfer function.
- the transfer function from output transducer 14 to error transducer 16 is modeled by a filter, e.g. an LMS or RLMS filter, C at 142, as in the incorporated '676 patent.
- Auxiliary random signal source 140 introduces a random signal into the output of model 40 at summer 152 and into the C model at 148.
- the auxiliary random signal from source 140 is random and uncorrelated with the system input signal 6, and is also uncorrelated with auxiliary random signal source 208, and in preferred form is provided by a Galois sequence, M.R. Schroeder, "Number Theory In Science And Communications", Berlin, Springer-Berlag, 1984, pages 252-261, though other random uncorrelated signal sources may be used.
- the Galois sequence is a pseudo random sequence that repeats after 2 M -1 points, where M is the number of stages in a shift register. The Galois sequence is preferred because it is easy to calculate and can easily have a period much longer than the response time of the system.
- Model 142 models the transfer function from output transducer 14 to error transducer 16, including the transfer function of each.
- the transfer function from output transducer 14 to error transducer 16 may be modeled without a random signal source, as in U.S. Patent 4,987,598, incorporated herein by reference.
- Auxiliary source 140 introduces an auxiliary random signal such that error transducer 16 also senses the auxiliary signal from the auxiliary source.
- the auxiliary signal may be introduced into the recursive loop of the A and B filters as in Fig.
- auxiliary signal may be introduced into the model after the recursive loop, i.e. introducing the auxiliary signal only to line 46, and not to line 47.
- copies of model 142 are provided at 144 and 146 to compensate the noted transfer function.
- the outputs of filters A and B are summed at summer 48, whose output is summed at summer 152 with the output of random signal source 140 to provide an output resultant sum which provides the model output at 46 supplying the noted correction signal to output transducer 14.
- the output of model 142 is summed at summer 64 with the output of error transducer 16, and the resultant sum supplied as the error input to model 142 and as an error input to model 40.
- the output of error transducer 16 may be supplied directly to an error input of model 40 without being supplied through summer 64.
- M copy 200, Fig. 18, is provided by a copy of A filter 12 at A copy 214, Fig. 19, and a copy of B filter 22 at B copy 216.
- Multiplier 218 multiplies the output 220 of A copy 214 and the input 222 of A copy 214, and supplies the output resultant product at 224 to summer 226.
- Multiplier 228 multiplies the output 230 of C copy 144 and the error input at 44, and supplies the output resultant product at 232 to summer 226.
- Summer 226 sums the inputs 232 and 224, and supplies the output resultant sum as weight update signal 74 to A filter 12.
- Multiplier 234 multiplies the output 236 of B copy 216 and the input 238 of B copy 216, and supplies the output resultant product at 240 to summer 242.
- Multiplier 244 multiplies the output 246 of C copy 146 and the error input at 44, and supplies the output resultant product at 248 to summer 242.
- the summer 242 sums the inputs 248 and 240, and supplies the output resultant sum as weight update signal 78 to B filter 22.
- the input to A copy 214 and to B copy 216 is provided by the output 304 of filter 302 receiving the model output at 47.
- the error signals at error inputs 232 and 224 oppositely drive the model.
- the error signal at error input 232 of the direct transfer function filter A drives the correction signal at 46 towards a value matching the system input signal 6.
- the error signal at error input 224 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero.
- this is accomplished by using a copy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero.
- the signal at error input 224 is provided only in response to a given condition of a given parameter, e.g. when the frequency is in the ranges passed by filter 302.
- the relative amplitudes of the input signals at error inputs 232 and 224 are adjusted such that the signal at error input 224 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree.
- the frequency is outside of the ranges passed by filter 302
- there is no output at 212, and hence no input to A copy 214 and hence the latter's output is undefined, whereby error input 232 from error signal 44 from error transducer 16 dominates and hence drives correction signal 46 to a value which matches the system input signal 6 to provide modification and/or control of the latter.
- the gain is controlled at control input 308 by the control signal supplied from peak detector 312 and filter 306 which supplies the spectrally processed correction signal based on the model output at 47.
- the error inputs 248 and 240 to the recursive transfer function filter B of model 40 function comparably to error inputs 232 and 224, respectively.
- Model 40 has a first error input provided at 232 and 248 from error transducer 16 driving the output of the model towards a value matching the system input signal 6.
- Model 40 has a second error input at 224 and 240 selectively driving the output of model 40 away from such matching value and instead driving the correction signal 46 towards zero.
- Fig. 20 uses like reference numerals from above and from Figs. 19 and 20 of the incorporated '676 patent.
- Summer 260 sums the output 230 of C copy 144 and the output of filter 302 which supplies the input to A copy 214, and supplies the output resultant sum at 262 to multiplier 264.
- Summer 266 sums the output 268 of A copy 214 and the error input at 44, and supplies the output resultant sum at 270 to multiplier 264.
- Multiplier 264 multiplies the inputs 262 and 270, and supplies the output resultant product as weight update signal 74 to A filter 12.
- Summer 272 sums the output 246 of C copy 146 and the output 304 of filter 302 which supplies the input to B copy 216, and supplies the output resultant sum at 274 to multiplier 276.
- Summer 278 sums the output 280 of B copy 216 and the error input at 44, and supplies the output resultant sum at 282 to multiplier 276.
- Multiplier 276 multiplies inputs 282 and 274, and supplies the output resultant product as weight update signal 78 to B filter 22.
- the error signals at 44 and 268 oppositely drive the model.
- the error signal at error input 44 of the direct transfer function filter A drives the correction signal at 46 towards a value matching the system input signal 6.
- the error signal at error input 268 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero.
- this is accomplished by using a copy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero.
- the signal at error input 268 is provided only in response to a given condition of a given parameter, e.g., when the frequency is in the ranges passed by filter 302.
- the relative amplitudes of the input signals at error inputs 44 and 268 to summer 266 are adjusted such that the signal at error input 268 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree.
- the frequency is outside of the ranges passed by filter 302
- there is no output at 212, and hence no input to A copy 214 and hence the latter's output is undefined, whereby error input 44 to summer 266 from error transducer 16 dominates and hence the model drives correction signal 46 to a value which matches the system input signal 6 to provide modification and/or control of the latter.
- the gain is controlled at control input 308 by the control signal provided from peak detector 312 and filter 306.
- the error inputs 44 and 280 at summer 278 to the recursive transfer function filter B of model 40 function comparably to error inputs 44 and 268 at summer 266, respectively.
- Model 40 in Fig. 6 has a first error input provided at 44 to summers 266 and 278 from error transducer 16 driving the output of the model towards a value matching the system input signal 6.
- Model 40 has a second error input at 268 and 280 selectively driving the output of model 40 away from such matching value and instead driving the correction signal 46 towards zero.
- Fig. 21 uses like reference numerals from Fig. 18 and illustrates a further embodiment.
- spectral transfer function filter 302 has its input 316 supplied from the reference signal from input transducer 10, rather than from the correction signal from the output of model 40, Fig. 18, or the random noise signal source 208, Fig. 15.
- the output of spectral transfer function filter 302 is supplied to filter input 304 of M copy 200, as above.
- the spectral transfer function provided by filter 306 includes a peak detector 312 similar to peak detector 97, Fig. 7.
- the output of inverse C model 306 is supplied to peak detector 312 comparing the output of inverse C model 306 with a desired peak value 314.
- the positive output of summer 312 at the control input 308 controls the gain of M copy 200 to in turn control the gain of the leak signal at filter input 304 supplied through the spectral transfer function 302 from the reference signal from input transducer 10.
- Filters 95, 302, 306 are each preferably an inverse C model copy provided by inverse S and/or inverse E, though other transfer functions may be used for any or all of such filters.
- performance of model 40 is controlled according to fuzzy logic to control the signal sent to output transducer 14.
- Fuzzy logic control is known in the prior art, for example "Adaptive Fuzzy Systems", E. Cox, IEEE Spectrum February 1993, pages 27-31, and the references noted therein at page 31, lower half of right column.
- fuzzy logic is used to provide self-designing control architecture using fuzzy rules and/or to control the filter weights which update the model and/or to control the leak signal which degrades performance of the model.
- the fuzzy logic controller having a given set of rules, is used for computing or setting the control architecture and/or the filter transfer function and/or the filter weights and/or the leak signal.
- Fuzzy logic enables control of model performance in a practical way for a multiplicity of factors which interact in a complex way.
- the fuzzy logic control is based upon relative values and qualitative trends, without requiring exact equation relationships for multiple variables. For example, referring to Figs. 22-24, which are similar to Fig. 4, page 29 of the above noted Cox article, if the magnitude of the weights is small, Fig. 22, and the magnitude of the correction signal output of model 40 is small, Fig. 23, then only a small leakage value is introduced, Fig. 24, to maximize the adaptive process. If the magnitude of the weights is large and the magnitude of the model output is large, then a large leak signal is introduced, to constrain model performance and minimize the adaptive process.
- the filter weights are monitored to provide a first input parameter, and the output of the model is monitored to provide a second input parameter. Fuzzy logic is applied to such first and second input parameters to provide an output parameter, Fig. 24, controlling the weights and/or leak and/or otherwise controlling performance of the model.
- the input parameters are fuzzified according to a fuzzy rule set, Figs. 22 and 23, a fuzzy leak is computed and then defuzzified using Fig. 24, to provide an output parameter to control performance of the model.
- the correction signal magnitude is nearing the capacity of output transducer 14 such as a loudspeaker, then leakage is increased. If the magnitude of the correction signal is not nearing the capacity of the loudspeaker, then leakage can be decreased.
- rate of change of the weights is monitored, and if the weights start increasing at a faster rate or at a rate above a given rate, then leakage is increased. If the weights are increasing rapidly and the magnitude of the correction signal is approaching the capacity of the loudspeaker, then leakage is increased very rapidly. If the weights start decreasing, then leakage is decreased.
- Inputs to the fuzzy logic include magnitude of the filter update weight, magnitude of the model output correction signal, rate of change of the filter weights, rate of change of the correction signal output by the model, magnitude of the reference signal 42 input to the model, rate of change of input signal 42, magnitude of the error signal from error transducer 16, rate of change of the error signal, spectral characteristics of the reference signal and/or correction signal and/or error signal, capability of output transducer 14, temperature, flow rates, environmental variables, fan speed in a duct application, whether or not such fan is running, thermostat settings, desired speed of adaptation, desired algorithm stability, system plant information, source information, etc. Fuzzy logic is a means of introducing human-like intuition into the controller adjustment process. Some factors may override others. For example, if the rate of increase of the weights is too large, then the algorithm may become unstable, and hence it is desirable to increase the leak or decrease the model update gain. If the capacity of a loudspeaker 14 is rapidly being reached, then leakage should be increased.
- the invention is applicable to multi-channel active acoustic attenuation systems, for example as shown in U.S. Patents 5,216,721 and 5,216,722, incorporated herein by reference.
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
- Feedback Control In General (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
Description
- The invention relates to active adaptive control systems, and more particularly to improvements for selectively controlling performance of the active adaptive model.
- The invention arose during continuing development efforts relating to the subject matter of U.S. Patents 4,837,834, 5,172,416, 5,278,913, 5,386,477, 5,390,255, 5,396,561, and co-pending U.S. Applications Serial No. 08/166,698, filed December 14, 1993, (EP-A-0 661 807), 08/247,561, filed May 23, 1994 (EP-A-0 684 594), 08/264,510, filed June 23, 1994, 08/340,613, filed November 16, 1994, 08/369,925, filed January 6, 1995, (EP-A-0 721 179) incorporated herein by reference.
- Active acoustic attenuation involves injecting a canceling acoustic wave to destructively interfere with and cancel an input acoustic wave. In an active acoustic attenuation system, the output acoustic wave is sensed with an error transducer, such as a microphone or an accelerometer, which supplies an error signal to an adaptive filter control model which in turn supplies a correction signal to a canceling output transducer or actuator, such as a loudspeaker or a shaker, which injects an acoustic wave to destructively interfere with the input acoustic wave and cancel or reduce same such that the output acoustic wave at the error transducer is zero or some other desired value.
- An active adaptive control system minimizes an error signal by introducing a control signal from an output transducer to combine with the system input signal and yield a system output signal. The system output signal is sensed with an error transducer providing the error signal. An adaptive filter model has a model input from a reference signal correlated with the system input signal, an error input from the error signal, and outputs a correction signal to the output transducer to introduce a control signal matching the system input signal, to minimize the error signal. The filter coefficients are updated according to a weight update signal which is the product of the reference signal and the error signal.
- The present invention is applicable to active adaptive control systems, including active acoustic attenuation systems. In one embodiment, the invention maximizes model performance by concentrating model adaptation in frequency ranges of interest, and protects the output transducer against overdriving of same. Performance of the model is spectrally controlled to maximize the correction signal sent to the output transducer such that at frequencies where maximum power from the output transducer reaches the error transducer, the correction signal supplied to the output transducer is maximized, and at frequencies where minimum power from the output transducer reaches the error transducer, the correction signal supplied to the output transducer is minimized. This maximizes model performance by concentrating model adaptation on the portion of the input signal which the model can control or cancel or where it is desired to do so, and constrains model adaptation as to those portions of the input signal which it cannot cancel or control or where it is not desired to do so. The latter is desired for stability of the model algorithm where active control solutions sometimes require more actuator power than is available or desirable. Actuators, amplifiers, etc. have limitations that adversely affect control algorithms. Pushed beyond capacity, the control output or power available from the secondary source or output transducer may exhibit saturation, clipping, or otherwise nonlinear behavior. Excessive control effort can result in damaged actuators, excessive power consumption, and instability in the control algorithm.
- It is known in the prior art to provide weight update signal leakage to counteract the adaptive process. This is done by implementing an exponential decay of the filter coefficients, intentionally defeating control effort, Widrow and Stearns, Adaptive Signal Processing, Prentice-Hall, Inc., Engelwood Cliffs, NJ, 1985, pages 376-378. The exponential decay is typically selected to be slow such that the adaptive process toward a control solution dominates. A deficiency of this method is that it unilaterally, across all power levels and frequencies, degrades performance. Such leakage is useful for limiting control effort and enhancing numerical stability, but performance suffers because of the lack of consideration for regions where the control effort is in an acceptable range.
- In other embodiments of the invention, model adaptation is selectively controlled to provide desired performance. Model performance is controlled by fuzzy logic to provide self-designing control architecture using fuzzy rules and/or to control a filter transfer function and/or to control filter weights used in an update process for feedforward and/or feedback, including FIR (finite impulse response) and IIR (infinite impulse response) applications and/or to control magnitude and/or rate of change of a leak signal degrading performance of the model.
- Fig. 1 is a schematic illustration of an active adaptive control system known in the prior art.
- Fig. 2 is a schematic illustration of an active adaptive control system in accordance with co-pending U.S. Application S.N. 08/264,510, filed June 23, 1994.
- Fig. 3 is a graph showing performance of the system of Fig. 2.
- Fig. 4 is a graph further showing performance of the system of Fig. 2.
- Fig. 5 is a graph showing an alternate performance of the system of Fig. 2.
- Fig. 6 is a graph further showing alternate performance of the system of Fig. 2.
- Fig. 7 is a schematic illustration of an active adaptive control system in accordance with the present invention.
- Fig. 8 is a graph showing frequency versus output and illustrates performance of the system of Fig. 7.
- Fig. 9 is a schematic illustration of an active adaptive control system.
- Fig. 10 is a schematic illustration of an adaptive filter model, and illustrates a principle employed by a system in accordance with copending U.S. Application Serial No. 08/166,698, filed December 14, 1993.
- Fig. 11 is like Fig. 10 and shows another manner of implementing the principle thereof.
- Fig. 12 is a schematic illustration of an active adaptive control system in accordance with the system of the '698 application.
- Fig. 13 shows a further embodiment of the system of Fig. 12.
- Fig. 14 shows a further embodiment of the system of Fig. 12.
- Fig. 15 is a schematic illustration of an active adaptive control system in accordance with the present invention.
- Fig. 16 shows a further embodiment of the system of Fig. 15.
- Fig. 17 shows a further embodiment of the system of Fig. 15.
- Fig. 18 is a schematic illustration of an active adaptive control system in accordance with the present invention.
- Fig. 19 shows a further embodiment of the system of Fig. 18.
- Fig. 20 shows a further embodiment of the system of Fig. 18.
- Fig. 21 shows a further embodiment of the system of Fig. 18.
- Figs. 22-24 illustrate a further embodiment of the invention.
- Fig. 1 shows an active adaptive control system similar to that shown in U.S. Patent 4,677,676, incorporated herein by reference, and uses like reference numerals therefrom where appropriate to facilitate understanding. The system introduces a control signal from a secondary source or
output transducer 14, such as a loudspeaker, shaker, or other actuator or controller, to combine with thesystem input signal 6 and yield asystem output signal 8. Aninput transducer 10, such as a microphone, accelerometer, or other sensor, senses the system input signal and provides areference signal 42. Anerror transducer 16, such as a microphone, accelerometer, or other sensor, senses the system output signal and provides anerror signal 44.Adaptive filter model 40 adaptively models the system and has a model input fromreference signal 42 correlated tosystem input signal 6, and an output outputting acorrection signal 46 tooutput transducer 14 to introduce the control signal according to aweight update signal 74.Reference signal 42 anderror signal 44 are combined atmultiplier 72 to provide the weight update signal throughdelay element 73. In a known alternative, thereference signal 42 may be provided by one or more error signals, in the case of a periodic system input signal, "Active Adaptive Sound Control In A Duct: A Computer Simulation", J.C. Burgess, Journal of Acoustic Society of America, 70(3), September 1981, pages 715-726, U.S. Patents 5,206,911, 5,216,722, incorporated herein by reference. - In updating the filter coefficients, and as is standard, one or more previous weights are added to the current product of
reference signal 42 anderror signal 44 atsummer 75. As noted above, it is known in the prior art to provide exponential decay of all of the filter coefficients in the system. Leakage factor γ at 77 multiplies one or more previous weights, after passage through one ormore delay elements 73, by an exponential decay factor less than one before adding same atsummer 75 to the current product ofreference signal 42 anderror signal 44, Adaptive Signal Processing, Widrow and Stearns, Prentice-Hall, Inc., Engelwood Cliffs, NJ, 1985, pages 376-378, including equations 13.27 and 13.31. As noted above, a deficiency of this method is that it reduces control effort and degrades performance across all power levels, regardless of whether such reduced effort is desired. - In the '510 application, selective leakage of the weight update signal is provided in response to a given condition of a given parameter, to control performance of the model on an as needed basis. In the preferred embodiment, leakage is varied as a function of
correction signal 46. A variable leakage factor γ is provided at 79 in Fig. 2, replacing fixedγ 77 of Fig. 1. Leakage factor γ at 79 is varied from a maximum value of 1.0 affording maximum control effort, to a minimum value such as zero providing minimum control effort. - Leakage is varied as a function of the output power of
correction signal 46 supplied from the output ofmodel 40 tooutput transducer 14. In the embodiment in Fig. 3, the leakage is varied as a discontinuous step function of the output power of the correction signal. When the output power exceeds a given threshold at 81, γ is abruptly, nonlinearly changed as a step function from afirst level 83 to asecond level 85. The reduction at 85 reduces the weight update signal summed atsummer 75 with the product of thereference signal 42 and error signal 44 frommultiplier 72, and hence reduces the weight update signal supplied tomodel 40. The noted reduction of γ atthreshold 81 increases leakage of the weight update signal, Fig. 4, fromlevel 87 tolevel 89. - In another embodiment as shown in Fig. 5, leakage is varied as a continuous function of the output power of the correction signal. In Fig. 5, γ is maintained at
level 83 until output power reachesthreshold 81, and then is linearly decreased as shown at 91 as a continuous linearly changing value as a function of increasing output power abovethreshold 81. As shown in Fig. 6, leakage is maintained atlevel 87 until output power reachesthreshold 81, and then is linearly increased at 93 as a continuous linearly changing value as a function of increasing output power abovethreshold 81. - Other variations of leakage are possible for providing selective leakage of the weight update signal to degrade performance of the model. The leakage is adjustably varied to vary performance of the model by multiplying a previous weight update value by
variable γ 79 and adding the result atsummer 75 to the product ofreference signal 42 and error signal 44 frommultiplier 72.γ 79 is varied as a function ofcorrection signal 46, preferably the output power of such correction signal. - Fig. 7 illustrates the present invention and uses like reference numerals from Fig. 2 and from Figs. 19 and 20 of the incorporated '676 patent. The transfer function from
output transducer 14 toerror transducer 16 is modeled with an adaptive filter model C at 142, as in the incorporated '676 patent. Auxiliaryrandom signal source 140 introduces a random signal into the output ofmodel 40 atsummer 152 and into the C model at 148. The auxiliary random signal fromsource 140 is random and uncorrelated with thesystem input signal 6 and in preferred form is provided by a Galois sequence, M.R. Schroeder, "Number Theory In Science And Communications", Berlin, Springer-Berlag, 1984, pages 252-261, though other random uncorrelated signal sources may be used. The Galois sequence is a pseudo random sequence that repeats after 2M-1 points, where M is the number of stages in a shift register. The Galois sequence is preferred because it is easy to calculate and can easily have a period much longer than the response time of the system. Theinput 148 toC model 142 is multiplied with the error signal fromerror transducer 16 atmultiplier 68, and the resultant product provided asweight update signal 67.Model 142 models the transfer function fromoutput transducer 14 toerror transducer 16, including the transfer function of each. Alternatively, the transfer function fromoutput transducer 14 toerror transducer 16 may be modeled without a random signal source, as in U.S. Patent 4,987,598, incorporated herein by reference.Auxiliary source 140 introduces an auxiliary random signal such thaterror transducer 16 also senses the auxiliary signal from the auxiliary source. The auxiliary signal may be introduced into the recursive loop of the A and B filters as in Fig. 19 of the incorporated '676 patent atsummer 152, or alternatively the auxiliary signal may be introduced into the model after the recursive loop, i.e. introducing the auxiliary signal only toline 46, and not to line 47. A copy ofmodel 142 is provided at 145 to compensate the noted transfer function, as in the incorporated '676 patent. - As in the '510 application, the present system and method involves introducing a control signal from
output transducer 14 to combine withsystem input signal 6 and yieldsystem output signal 8, sensing the system output signal witherror transducer 16 and providing anerror signal 44, providingadaptive filter model 40 having a model input fromreference signal 42 correlated tosystem input signal 6, and an output outputting acorrection signal 46 tooutput transducer 14 to introduce the control signal according toweight update signal 74. In the present invention, performance ofmodel 40 is spectrally controlled to maximize the signal sent tooutput transducer 14 at frequencies of interest or where it can be maximally effective, and minimize the signal sent tooutput transducer 14 at frequencies of non interest or where it is only minimally effective or is ineffective. A spectral leak signal is provided which degrades performance of the model. The leak signal is controlled according to frequency. In Fig. 7, the correction signal from the output ofmodel 40 is monitored, and the leak signal is controlled in response thereto. The correction signal is filtered byfilter 95 to provide the leak signal. The correction signal from the output ofmodel 40 is spectrally processed by supplying the correction signal through the frequency responsive spectral transfer function provided byfilter 95 to provide the leak signal to the error input ofmodel 40. In Fig. 7, the correction signal is spectrally processed such that at frequencies where maximum power fromoutput transducer 14 reacheserror transducer 16, the correction signal is maximized, and at frequencies where minimum power fromoutput transducer 14 reacheserror transducer 16, the correction signal is minimized. - The transfer function between
output transducer 14 anderror transducer 16 is modeled withC model 142.Correction signal 46 is spectrally processed by a function ofmodel 142. In Fig. 7, such function is the inverse ofC model 142 as provided at inverse C model, C-1, at 95. The output ofinverse C model 95 is supplied to an optional peak detector as provided bysummer 97 comparing the output ofinverse C model 95 with a desiredpeak value 99. When the output ofmodel 95 rises abovelevel 99, the positive output ofsummer 97 controls variable leakage factor γ at 79, as above. The inverse C model includes the inverse of the transfer function ofoutput transducer 14, inverse S, S-1, and/or the inverse of the transfer function of the error path, inverse E, E-1, betweenoutput transducer 14 anderror transducer 16. - In addition to or in place of
peak detector 97, control logic may be used to respond to the output ofinverse C model 95 and control leakage factor γ at 79 according to designated conditions or rules to generate or compute a leak value or control the leaking process, to be further described hereinafter. In other embodiments,filter 95 is a displacement function ofoutput transducer 14 such as a loudspeaker, to protect the latter against overdriving. In other embodiments, an RMS (root mean square) to DC (direct current) conversion function is provided betweenfilter 95 andpeak detector 97 for average level detection to control convergence rate of the leak process. The transfer function offilter 95 may be linear or nonlinear. - As seen in Fig. 8, at frequencies where maximum power from
output transducer 14 reacheserror transducer 16, as atfrequency regions C model 142 has a maximum transfer characteristic, and the inverse of the C model at 95 has a minimum transfer characteristic, as shown at 105 and 107. The minimum transfer characteristic at 105 and 107 minimizes leakage of the update signal to model 40, to enable maximum output ofmodel 40. At frequencies where minimum power fromoutput transducer 14 reacheserror transducer 16, as atregion 109,C model 142 has a minimum transfer characteristic, and the inverse of the C model at 95 has a maximum transfer characteristic, as shown at 111. The maximum transfer characteristic at 111 maximizes leakage of the update signal, to minimize the output ofmodel 40.Inverse C model 95 spectrallysenses correction signal 46 and provides selective leakage ofweight update signal 74 in response thereto, to control performance ofmodel 40 according to frequency, to optimize performance ofmodel 40 in frequency ranges such as 101 and 103 where the model can effectively control or cancelinput signal 6. Outside of such ranges, e.g. at 218,inverse C model 95 minimizes performance ofmodel 40, to avoid using computation or adaptation power where inefficient or unneeded, and to prevent the model from continually trying to generate an output in regions where it is ineffective to attempt to control or cancelinput signal 6.Model 95 affords frequency weighting of the weight update signal. - Fig. 9 is similar to Fig. 5 of U.S. Patent 4,677,676, incorporated herein by reference, and uses like reference numerals to facilitate understanding. The system introduces a control signal from an
output transducer 14, such as a loudspeaker, shaker, or other actuator or controller, to combine with thesystem input signal 6 and yield asystem output signal 8. Aninput transducer 10, such as a microphone, accelerometer, or other sensor, senses the system input signal. Anerror transducer 16, such as a microphone, accelerometer, or other sensor, senses the system output signal and provides an error signal.Adaptive filter model 40 adaptively models the system and has amodel input 42 frominput transducer 10, anerror input 44 fromerror transducer 16, and amodel output 46 outputting a correction signal tooutput transducer 14 to introduce the control signal. In a known alternative, the input signal at 42 may be provided by one or more error signals, in the case of a periodic system input signal, "Active Adaptive Sound Control In A Duct: A Computer Simulation", J.C. Burgess, Journal of Acoustic Society of America, 70(3), September 1981, pages 715-726, U.S. Patents 5,206,911, 5,216,722, incorporated herein by reference. - The system of the '698 application provides an active adaptive control system wherein the performance of
model 40 is intentionally and selectively constrained by driving theoutput 46 of the model towards zero in response to a given condition of a given parameter. For example, in an active noise control system, it may be desirable to cancel noise only in a given frequency band, and leave the noise uncanceled for frequencies outside the band. In other control applications, it may be desirable to selectively control the system output signal by selectively controlling introduction of the control signal fromoutput transducer 14 to match or not match the system input signal. - One manner of constraining system performance is to drive the output of
model 40 towards zero and away from a value matchingsystem input signal 6. One way of accomplishing this is shown in Fig. 10 wherein the output ofmodel 40 is supplied to its error input, such that when the model adapts to drive its error input towards zero, the output of the model is necessarily also driven towards zero. Fig. 11 shows another manner of implementing this principle wherein a copy of the model is provided at 200, and the output ofmodel copy 200 supplies the error signal to the error input ofmodel 40. In Fig. 11,model 40 adapts to drive its error input towards zero, which in turn requires that the output ofcopy 200 be driven towards zero, which in turn means that the output ofmodel 40 is driven towards zero becauseM copy 200 is a duplicate ofmodel 40. These principles are utilized in the present invention. -
Model 40, Fig. 12, normally adapts to a converged condition wherein its output at 46 provides a correction signal tooutput transducer 14 which outputs a control signal matching the system input signal or a designated relative value correlated thereto. For example, in a noise cancellation system, the matching control signal fromoutput transducer 14 cancels the input noise. In the system of the '698 application, in response to a given condition of a given parameter, the output ofmodel 40 is driven towards zero by driving the output of the model towards its error input, such that when the model adapts to drive the error signal towards zero, the output of the model is also driven towards zero. This is accomplished by providing acopy 200 ofmodel 40, and supplying the output of the copy to anerror input 202 of the model which is summed atsummer 204 with the error signal aterror input 44 fromerror transducer 16. The model adapts to drive the error input towards zero which in turn requires that the output ofcopy 200 and hence the output ofmodel 40 are driven towards zero, to provide the noted constrained performance. The driving ofmodel output 46 towards zero provides a zero or at least a reduced correction signal tooutput transducer 14 to constrain or reduce modification and/or cancellation of thesystem input signal 6. - A random signal is provided at 206 from an auxiliary
random signal source 208, preferably provided by a Galois sequence, M.R. Schroeder, "Number Theory In Science And Communications", Berlin, Springer-Berlag, 1984, pages 252-261, though other random signal sources may be used, uncorrelated with thesystem input signal 6. The Galois sequence is a pseudo random sequence that repeats after 2M-1 points, where M is the number of stages in a shift register. The Galois sequence is preferred because it is easy to calculate and can easily have a period much longer than the response time of the system. The random signal is supplied through astopband filter 210 tomodel copy 200 at 212.Stopband filter 210 blocks frequencies in the stopband, and passes frequencies outside the stopband. This arrangement provides a spectral leak signal at 202 in response to a given condition of a given parameter, for example a frequency outside the stopband. In such implementation, the noted given parameter is frequency, and the given condition is a designated sub-optimum performance band outside the stopband. - The spectral leak error signal at 202 drives the correction signal at
model output 46 towards zero and provides sub-optimum performance ofmodel 40. Outside of the sub-optimum performance band, i.e. within the stopband offilter 210, there is no signal at 212 and hence the output ofcopy 200 is undefined, and the error signal fromerror transducer 16 aterror input 44 is maximally effective andmodel 40 optimally responds thereto and drives the correction signal atoutput 46 toward a value matching thesystem input signal 6. When the spectral leak signal is present aterror input 202, it constrains performance ofmodel 40 by driving or at least attempting to drive the correction signal atmodel output 46 towards zero. The relative influence or amplitudes of the error signals aterror inputs filter 210, then the noted relative amplitudes are set such that the error signal aterror input 202 dominates the error signal aterror input 44, and hence the correction signal atmodel output 46 is driven towards zero and away from a value matching thesystem input signal 6. - The method of the '698 application involves driving
error input 44 to drive the correction signal atmodel output 46 toward a value matching the system input signal, and selectively drivingerror input 202 to drive the correction signal atmodel output 46 away from the matching value by driving the correction signal towards zero. The arrangement provides a spectral leak signal to errorinput 202 in response to the noted given condition of the given parameter, e.g. a frequency outside the stopband, such that in the presence of the given condition, the spectral leak signal drives the correction signal atmodel output 46 towards zero, and in the absence of the given condition the error signal aterror input 44 drives the correction signal atmodel output 46 towards a value matching thesystem input signal 6. -
Stopband filter 210 blocks frequencies in a given stopband at which modification or cancellation of thesystem input signal 6 bymodel 40 is desired. Filter 210 passes frequencies in a given passband at which modification or cancellation of the system input signal bymodel 40 is undesired or not possible. The control signal output byoutput transducer 14 is driven toward a value matching thesystem input signal 6 only for frequencies in the stopband. At frequencies in the stopband, the error signal aterror input 44 is dominant, and the control signal output byoutput transducer 14 is driven toward a value matching thesystem input signal 6. At frequencies in the passband, the error signal aterror input 202 is dominant, and the control signal output byoutput transducer 14 is driven away from a value matching thesystem input signal 6. - Fig. 13 is similar to Figs. 19 and 20 of the incorporated '676 patent, and uses like reference numerals where appropriate to facilitate understanding. As noted in the incorporated '676 patent, model M at 40 is preferably an adaptive recursive filter having a transfer function with both poles and zeros. Model M is provided by an IIR, infinite impulse response, filter, e.g. a recursive least mean square, RLMS, filter having a first algorithm filter provided by an FIR, finite impulse response, filter, e.g. a least mean square, LMS, filter A at 12, and a second algorithm filter provided by an FIR filter, e.g. an LMS filter, B at 22. Filter A provides a direct transfer function, and filter B provides a recursive transfer function. The transfer function from
output transducer 14 toerror transducer 16 is modeled by a filter, e.g. an LMS or RLMS filter, C at 142, as in the incorporated '676 patent. - Auxiliary
random signal source 140 introduces a random signal into the output ofmodel 40 atsummer 152 and into the C model at 148. The auxiliary random signal fromsource 140 is random and uncorrelated with thesystem input signal 6, and is also uncorrelated with auxiliaryrandom signal source 208, and in preferred form is provided by a Galois sequence, M.R. Schroeder, "Number Theory In Science And Communications", Berlin, Springer-Berlag, 1984, pages 252-261, though other random uncorrelated signal sources may be used. The Galois sequence is a pseudo random sequence that repeats after 2M-1 points, where M is the number of stages in a shift register. The Galois sequence is preferred because it is easy to calculate and can easily have a period much longer than the response time of the system. Theinput 148 toC model 142 is multiplied with the error signal fromerror transducer 16 atmultiplier 68, and the resultant product provided asweight update signal 67.Model 142 models the transfer function fromoutput transducer 14 toerror transducer 16, including the transfer function of each. Alternatively, the transfer function fromoutput transducer 14 toerror transducer 16 may be modeled without a random signal source, as in U.S. Patent 4,987,598, incorporated herein by reference.Auxiliary source 140 introduces an auxiliary random signal such thaterror transducer 16 also senses the auxiliary signal from the auxiliary source. The auxiliary signal may be introduced into the recursive loop of the A and B filters as in Fig. 19 of the incorporated '676 patent atsummer 152, or alternatively the auxiliary signal may be introduced into the model after the recursive loop, i.e. introducing the auxiliary signal only toline 46, and not to line 47. As in the incorporated '676 patent, copies ofmodel 142 are provided at 144 and 146 to compensate the noted transfer function. - The outputs of filters A and B are summed at
summer 48, whose output is summed atsummer 152 with the output ofrandom signal source 140 to provide an output resultant sum which provides the model output at 46 supplying the noted correction signal tooutput transducer 14. The output ofmodel 142 is summed atsummer 64 with the output oferror transducer 16, and the resultant sum supplied as the error input to model 142 and as an error input to model 40. Alternatively, the output oferror transducer 16 may be supplied directly to an error input ofmodel 40 without being supplied throughsummer 64. -
M copy 200, Fig. 12, is provided by a copy of Afilter 12 at Acopy 214, Fig. 13, and a copy ofB filter 22 atB copy 216.Multiplier 218 multiplies theoutput 220 of Acopy 214 and theinput 222 of Acopy 214, and supplies the output resultant product at 224 tosummer 226.Multiplier 228 multiplies theoutput 230 ofC copy 144 and the error input at 44, and supplies the output resultant product at 232 tosummer 226.Summer 226 sums theinputs weight update signal 74 to Afilter 12.Multiplier 234 multiplies theoutput 236 ofB copy 216 and theinput 238 ofB copy 216, and supplies the output resultant product at 240 tosummer 242.Multiplier 244 multiplies theoutput 246 ofC copy 146 and the error input at 44, and supplies the output resultant product at 248 tosummer 242. Thesummer 242 sums theinputs weight update signal 78 toB filter 22. The input to Acopy 214 and toB copy 216 is provided by theoutput 212 ofstopband filter 210 receiving the noted random input signal at 206 fromrandom signal source 208. - The error signals at
error inputs error input 232 of the direct transfer function filter A drives the correction signal at 46 towards a value matching thesystem input signal 6. The error signal aterror input 224 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero. As noted above, this is accomplished by using acopy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero. The signal aterror input 224 is provided only in response to a given condition of a given parameter, e.g. when the frequency is outside the stopband offilter 210. The relative amplitudes of the input signals aterror inputs error input 224 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree. When the frequency is in the stopband offilter 210, there is no output at 212, and hence no input to Acopy 214 and hence the latter's output is undefined, wherebyerror input 232 fromerror signal 44 fromerror transducer 16 dominates and hence drivescorrection signal 46 to a value which matches thesystem input signal 6 to provide modification and/or control of the latter. Theerror inputs model 40 function comparably to errorinputs Model 40 has a first error input provided at 232 and 248 fromerror transducer 16 driving the output of the model towards a value matching thesystem input signal 6.Model 40 has a second error input at 224 and 240 selectively driving the output ofmodel 40 away from such matching value and instead driving thecorrection signal 46 towards zero. - Fig. 14 is similar to Fig. 13 and uses like reference numerals where appropriate to facilitate understanding.
Summer 260 sums theoutput 230 ofC copy 144 and theoutput 212 ofstopband filter 210 which supplies the input to Acopy 214, and supplies the output resultant sum at 262 tomultiplier 264.Summer 266 sums theoutput 268 of Acopy 214 and the error input at 44, and supplies the output resultant sum at 270 tomultiplier 264.Multiplier 264 multiplies theinputs weight update signal 74 to Afilter 12.Summer 272 sums theoutput 246 ofC copy 146 and theoutput 212 ofstopband filter 210 which supplies the input toB copy 216, and supplies the output resultant sum at 274 tomultiplier 276.Summer 278 sums theoutput 280 ofB copy 216 and the error input at 44, and supplies the output resultant sum at 282 tomultiplier 276.Multiplier 276 multipliesinputs weight update signal 78 toB filter 22. - The error signals at 44 and 268 oppositely drive the model. The error signal at
error input 44 of the direct transfer function filter A drives the correction signal at 46 towards a value matching thesystem input signal 6. The error signal aterror input 268 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero. As noted above, this is accomplished by using acopy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero. The signal aterror input 268 is provided only in response to a given condition of a given parameter, e.g., when the frequency is outside the stopband offilter 210. The relative amplitudes of the input signals aterror inputs summer 266 are adjusted such that the signal aterror input 268 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree. When the frequency is in the stopband offilter 210, there is no output at 212, and hence no input to Acopy 214 and hence the latter's output is undefined, wherebyerror input 44 tosummer 266 fromerror transducer 16 dominates and hence the model drivescorrection signal 46 to a value which matches thesystem input signal 6 to provide modification and/or control of the latter. Theerror inputs summer 278 to the recursive transfer function filter B ofmodel 40 function comparably to errorinputs summer 266, respectively.Model 40 in Fig. 14 has a first error input provided at 44 tosummers error transducer 16 driving the output of the model towards a value matching thesystem input signal 6.Model 40 has a second error input at 268 and 280 selectively driving the output ofmodel 40 away from such matching value and instead driving thecorrection signal 46 towards zero. - Fig. 15 illustrates the present invention and uses like reference numerals from above and from Figs. 19 and 20 of the incorporated '676 patent. Random noise signal 206 from
random noise source 208 is supplied through a first inverseC model copy 302 whose output is supplied to thefilter input 304 ofM copy 200.Correction signal 46 from the output ofmodel 40 throughsummer 152 is supplied through a second inverseC model copy 306 to the convergence rategain control input 308 ofM copy 200. The spectral transfer function provided byfilter 302 spectrally controls performance ofmodel 40 to maximize the correction signal at frequencies of interest. In preferred form, the spectral transfer function is the inverse C model copy having the characteristics shown in Fig. 8, though other frequency responsive spectral transfer functions may be used. The spectral leak signal provided throughtransfer functions model 40, as above described in conjunction with Fig. 12.Transfer function 302 controls the leak signal according to frequency. It is preferred that the leak signal also be controlled in response to the correction signal from the output ofmodel 40.Transfer function 306 spectrally monitors the correction signal and filters same to provide a spectral control signal controlling the gain throughM copy 200.Transfer function 306 and the gain control provided atcontrol input 308 are optional. The spectral leak control provided throughtransfer function 302 to filterinput 304 may be used with or withouttransfer function 306 and controlinput 308.Filter input 304 ofM copy 200 receives a spectrally processed input signal fromspectral transfer function 302.Control input 308 ofM copy 200 controls gain ofM copy 200 and receives a spectral control signal fromtransfer function 306. The spectral control signal is responsive tocorrection signal 46.Spectral transfer functions - Fig. 16 uses like reference numerals from above and from Figs. 19 and 20 of the incorporated '676 patent. As noted in the incorporated '676 patent, model M at 40 is preferably an adaptive recursive filter having a transfer function with both poles and zeros. Model M is provided by an IIR, infinite impulse response, filter, e.g. a recursive least mean square, RLMS, filter having a first algorithm filter provided by an FIR, finite impulse response, filter, e.g. a least mean square, LMS, filter A at 12, and a second algorithm filter provided by an FIR filter, e.g. an LMS filter, B at 22. Filter A provides a direct transfer function, and filter B provides a recursive transfer function. The transfer function from
output transducer 14 toerror transducer 16 is modeled by a filter, e.g. an LMS or RLMS filter, C at 142, as in the incorporated '676 patent. - Auxiliary
random signal source 140 introduces a random signal into the output ofmodel 40 atsummer 152 and into the C model at 148. The auxiliary random signal fromsource 140 is random and uncorrelated with thesystem input signal 6, and is also uncorrelated with auxiliaryrandom signal source 208, and in preferred form is provided by a Galois sequence, M.R. Schroeder, "Number Theory In Science And Communications", Berlin, Springer-Berlag, 1984, pages 252-261, though other random uncorrelated signal sources may be used. The Galois sequence is a pseudo random sequence that repeats after 2M-1 points, where M is the number of stages in a shift register. The Galois sequence is preferred because it is easy to calculate and can easily have a period much longer than the response time of the system. Theinput 148 toC model 142 is multiplied with the error signal fromerror transducer 16 atmultiplier 68, and the resultant product provided asweight update signal 67.Model 142 models the transfer function fromoutput transducer 14 toerror transducer 16, including the transfer function of each. Alternatively, the transfer function fromoutput transducer 14 toerror transducer 16 may be modeled without a random signal source, as in U.S. Patent 4,987,598, incorporated herein by reference.Auxiliary source 140 introduces an auxiliary random signal such thaterror transducer 16 also senses the auxiliary signal from the auxiliary source. The auxiliary signal may be introduced into the recursive loop of the A and B filters as in Fig. 19 of the incorporated '676 patent atsummer 152, or alternatively the auxiliary signal may be introduced into the model after the recursive loop, i.e. introducing the auxiliary signal only toline 46, and not to line 47. As in the incorporated '676 patent, copies ofmodel 142 are provided at 144 and 146 to compensate the noted transfer function. - The outputs of filters A and B are summed at
summer 48, whose output is summed atsummer 152 with the output ofrandom signal source 140 to provide an output resultant sum which provides the model output at 46 supplying the noted correction signal tooutput transducer 14. The output ofmodel 142 is summed atsummer 64 with the output oferror transducer 16, and the resultant sum supplied as the error input to model 142 and as an error input to model 40. Alternatively, the output oferror transducer 16 may be supplied directly to an error input ofmodel 40 without being supplied throughsummer 64. -
M copy 200, Fig. 15, is provided by a copy of Afilter 12 at Acopy 214, Fig. 16, and a copy ofB filter 22 atB copy 216.Multiplier 218 multiplies theoutput 220 of Acopy 214 and theinput 222 of Acopy 214, and supplies the output resultant product at 224 tosummer 226.Multiplier 228 multiplies theoutput 230 ofC copy 144 and the error input at 44, and supplies the output resultant product at 232 tosummer 226.Summer 226 sums theinputs weight update signal 74 to Afilter 12.Multiplier 234 multiplies theoutput 236 ofB copy 216 and theinput 238 ofB copy 216, and supplies the output resultant product at 240 tosummer 242.Multiplier 244 multiplies theoutput 246 ofC copy 146 and the error input at 44, and supplies the output resultant product at 248 tosummer 242. Thesummer 242 sums theinputs weight update signal 78 toB filter 22. The input to Acopy 214 and toB copy 216 is provided by theoutput 304 offilter 302 receiving the noted random input signal at 206 fromrandom signal source 208. - The error signals at
error inputs error input 232 of the direct transfer function filter A drives the correction signal at 46 towards a value matching thesystem input signal 6. The error signal aterror input 224 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero. As noted above, this is accomplished by using acopy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero. The signal aterror input 224 is provided only in response to a given condition of a given parameter, e.g. when the frequency is in the ranges passed byfilter 302. The relative amplitudes of the input signals aterror inputs error input 224 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree. When the frequency is outside of the ranges passed byfilter 302, there is no output at 212, and hence no input to Acopy 214 and hence the latter's output is undefined, wherebyerror input 232 fromerror signal 44 fromerror transducer 16 dominates and hence drivescorrection signal 46 to a value which matches thesystem input signal 6 to provide modification and/or control of the latter. The gain is controlled by the control signal atcontrol input 308 which is the spectrally processed correction signal supplied through the transfer function filter provided by inverseC model copy 306. Theerror inputs model 40 function comparably to errorinputs Model 40 has a first error input provided at 232 and 248 fromerror transducer 16 driving the output of the model towards a value matching thesystem input signal 6.Model 40 has a second error input at 224 and 240 selectively driving the output ofmodel 40 away from such matching value and instead driving thecorrection signal 46 towards zero. - Fig. 17 uses like reference numerals from above and from Figs. 19 and 20 of the incorporated '676 patent.
Summer 260 sums theoutput 230 ofC copy 144 and theoutput 212 ofstopband filter 210 which supplies the input to Acopy 214, and supplies the output resultant sum at 262 tomultiplier 264.Summer 266 sums theoutput 268 of Acopy 214 and the error input at 44, and supplies the output resultant sum at 270 tomultiplier 264.Multiplier 264 multiplies theinputs weight update signal 74 to Afilter 12.Summer 272 sums theoutput 246 ofC copy 146 and theoutput 304 offilter 302 which supplies the input toB copy 216, and supplies the output resultant sum at 274 tomultiplier 276.Summer 278 sums theoutput 280 ofB copy 216 and the error input at 44, and supplies the output resultant sum at 282 tomultiplier 276.Multiplier 276 multipliesinputs weight update signal 78 toB filter 22. - The error signals at 44 and 268 oppositely drive the model. The error signal at
error input 44 of the direct transfer function filter A drives the correction signal at 46 towards a value matching thesystem input signal 6. The error signal aterror input 268 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero. As noted above, this is accomplished by using acopy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero. The signal aterror input 268 is provided only in response to a given condition of a given parameter, e.g., when the frequency is in the ranges passed bytransfer function filter 302. The relative amplitudes of the input signals aterror inputs summer 266 are adjusted such that the signal aterror input 268 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree. When the frequency is outside of the ranges passed byfilter 302, there is no output at 212, there is no output at 212, and hence no input to Acopy 214 and hence the latter's output is undefined, wherebyerror input 44 tosummer 266 fromerror transducer 16 dominates and hence the model drivescorrection signal 46 to a value which matches thesystem input signal 6 to provide modification and/or control of the latter. The gain is controlled by the control signal atcontrol input 308 provided by the spectrally processed correction signal supplied through the transfer function filter provided by inverseC model copy 306. Theerror inputs summer 278 to the recursive transfer function filter B ofmodel 40 function comparably to errorinputs summer 266, respectively.Model 40 in Fig. 6 has a first error input provided at 44 tosummers error transducer 16 driving the output of the model towards a value matching thesystem input signal 6.Model 40 has a second error input at 268 and 280 selectively driving the output ofmodel 40 away from such matching value and instead driving thecorrection signal 46 towards zero. - Fig. 18 illustrates the present invention and uses like reference numerals from above and from Figs. 19 and 20 of the incorporated '676 patent. In Fig. 18, spectral
transfer function filter 302 has itsinput 310 supplied from the correction signal from the output ofmodel 40, rather than from randomnoise signal source 208, Fig. 15. The spectral transfer function provided byfilter 306 also includes apeak detector 312 similar topeak detector 97, Fig. 7. The output ofinverse C model 306 is supplied to peakdetector 312 comparing the output ofinverse C model 306 with a desiredpeak value 314. When the output of inverseC model copy 306 rises above the level ofpeak value 314, the positive output ofsummer 312 at thecontrol input 308 controls the gain ofM copy 200 to in turn control the gain of the leak signal atfilter input 304 supplied through thespectral transfer function 302. - Fig. 19 uses like reference numerals from above and from Figs. 19 and 20 of the incorporated '676 patent. As noted in the incorporated '676 patent, model M at 40 is preferably an adaptive recursive filter having a transfer function with both poles and zeros. Model M is provided by an IIR, infinite impulse response, filter, e.g. a recursive least mean square, RLMS, filter having a first algorithm filter provided by an FIR, finite impulse response, filter, e.g. a least mean square, LMS, filter A at 12, and a second algorithm filter provided by an FIR filter, e.g. an LMS filter, B at 22. Filter A provides a direct transfer function, and filter B provides a recursive transfer function. The transfer function from
output transducer 14 toerror transducer 16 is modeled by a filter, e.g. an LMS or RLMS filter, C at 142, as in the incorporated '676 patent. - Auxiliary
random signal source 140 introduces a random signal into the output ofmodel 40 atsummer 152 and into the C model at 148. The auxiliary random signal fromsource 140 is random and uncorrelated with thesystem input signal 6, and is also uncorrelated with auxiliaryrandom signal source 208, and in preferred form is provided by a Galois sequence, M.R. Schroeder, "Number Theory In Science And Communications", Berlin, Springer-Berlag, 1984, pages 252-261, though other random uncorrelated signal sources may be used. The Galois sequence is a pseudo random sequence that repeats after 2M-1 points, where M is the number of stages in a shift register. The Galois sequence is preferred because it is easy to calculate and can easily have a period much longer than the response time of the system. Theinput 148 toC model 142 is multiplied with the error signal fromerror transducer 16 atmultiplier 68, and the resultant product provided asweight update signal 67.Model 142 models the transfer function fromoutput transducer 14 toerror transducer 16, including the transfer function of each. Alternatively, the transfer function fromoutput transducer 14 toerror transducer 16 may be modeled without a random signal source, as in U.S. Patent 4,987,598, incorporated herein by reference.Auxiliary source 140 introduces an auxiliary random signal such thaterror transducer 16 also senses the auxiliary signal from the auxiliary source. The auxiliary signal may be introduced into the recursive loop of the A and B filters as in Fig. 19 of the incorporated '676 patent atsummer 152, or alternatively the auxiliary signal may be introduced into the model after the recursive loop, i.e. introducing the auxiliary signal only toline 46, and not to line 47. As in the incorporated '676 patent, copies ofmodel 142 are provided at 144 and 146 to compensate the noted transfer function. - The outputs of filters A and B are summed at
summer 48, whose output is summed atsummer 152 with the output ofrandom signal source 140 to provide an output resultant sum which provides the model output at 46 supplying the noted correction signal tooutput transducer 14. The output ofmodel 142 is summed atsummer 64 with the output oferror transducer 16, and the resultant sum supplied as the error input to model 142 and as an error input to model 40. Alternatively, the output oferror transducer 16 may be supplied directly to an error input ofmodel 40 without being supplied throughsummer 64. -
M copy 200, Fig. 18, is provided by a copy of Afilter 12 at Acopy 214, Fig. 19, and a copy ofB filter 22 atB copy 216.Multiplier 218 multiplies theoutput 220 of Acopy 214 and theinput 222 of Acopy 214, and supplies the output resultant product at 224 tosummer 226.Multiplier 228 multiplies theoutput 230 ofC copy 144 and the error input at 44, and supplies the output resultant product at 232 tosummer 226.Summer 226 sums theinputs weight update signal 74 to Afilter 12.Multiplier 234 multiplies theoutput 236 ofB copy 216 and theinput 238 ofB copy 216, and supplies the output resultant product at 240 tosummer 242.Multiplier 244 multiplies theoutput 246 ofC copy 146 and the error input at 44, and supplies the output resultant product at 248 tosummer 242. Thesummer 242 sums theinputs weight update signal 78 toB filter 22. The input to Acopy 214 and toB copy 216 is provided by theoutput 304 offilter 302 receiving the model output at 47. - The error signals at
error inputs error input 232 of the direct transfer function filter A drives the correction signal at 46 towards a value matching thesystem input signal 6. The error signal aterror input 224 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero. As noted above, this is accomplished by using acopy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero. The signal aterror input 224 is provided only in response to a given condition of a given parameter, e.g. when the frequency is in the ranges passed byfilter 302. The relative amplitudes of the input signals aterror inputs error input 224 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree. When the frequency is outside of the ranges passed byfilter 302, there is no output at 212, and hence no input to Acopy 214 and hence the latter's output is undefined, wherebyerror input 232 fromerror signal 44 fromerror transducer 16 dominates and hence drivescorrection signal 46 to a value which matches thesystem input signal 6 to provide modification and/or control of the latter. The gain is controlled atcontrol input 308 by the control signal supplied frompeak detector 312 and filter 306 which supplies the spectrally processed correction signal based on the model output at 47. Theerror inputs model 40 function comparably to errorinputs Model 40 has a first error input provided at 232 and 248 fromerror transducer 16 driving the output of the model towards a value matching thesystem input signal 6.Model 40 has a second error input at 224 and 240 selectively driving the output ofmodel 40 away from such matching value and instead driving thecorrection signal 46 towards zero. - Fig. 20 uses like reference numerals from above and from Figs. 19 and 20 of the incorporated '676 patent.
Summer 260 sums theoutput 230 ofC copy 144 and the output offilter 302 which supplies the input to Acopy 214, and supplies the output resultant sum at 262 tomultiplier 264.Summer 266 sums theoutput 268 of Acopy 214 and the error input at 44, and supplies the output resultant sum at 270 tomultiplier 264.Multiplier 264 multiplies theinputs weight update signal 74 to Afilter 12.Summer 272 sums theoutput 246 ofC copy 146 and theoutput 304 offilter 302 which supplies the input toB copy 216, and supplies the output resultant sum at 274 tomultiplier 276.Summer 278 sums theoutput 280 ofB copy 216 and the error input at 44, and supplies the output resultant sum at 282 tomultiplier 276.Multiplier 276 multipliesinputs weight update signal 78 toB filter 22. - The error signals at 44 and 268 oppositely drive the model. The error signal at
error input 44 of the direct transfer function filter A drives the correction signal at 46 towards a value matching thesystem input signal 6. The error signal aterror input 268 of filter A drives the correction signal at 46 away from the noted matching value by driving the correction signal towards zero. As noted above, this is accomplished by using acopy 214 of the A filter and supplying the output of such copy as an error input to the adaptive model such that in attempting to drive the error input to zero, the model must drive its output to zero. The signal aterror input 268 is provided only in response to a given condition of a given parameter, e.g., when the frequency is in the ranges passed byfilter 302. The relative amplitudes of the input signals aterror inputs summer 266 are adjusted such that the signal aterror input 268 dominates when both are present, or the degree of dominance is adjusted to in turn adjust the amount of constrainment of performance of the model so that the correction signal at 46 is driven towards zero but never reaches zero, such that there is still some modification and/or cancellation of the system input signal, though to a reduced degree. When the frequency is outside of the ranges passed byfilter 302, there is no output at 212, and hence no input to Acopy 214 and hence the latter's output is undefined, wherebyerror input 44 tosummer 266 fromerror transducer 16 dominates and hence the model drivescorrection signal 46 to a value which matches thesystem input signal 6 to provide modification and/or control of the latter. The gain is controlled atcontrol input 308 by the control signal provided frompeak detector 312 andfilter 306. Theerror inputs summer 278 to the recursive transfer function filter B ofmodel 40 function comparably to errorinputs summer 266, respectively.Model 40 in Fig. 6 has a first error input provided at 44 tosummers error transducer 16 driving the output of the model towards a value matching thesystem input signal 6.Model 40 has a second error input at 268 and 280 selectively driving the output ofmodel 40 away from such matching value and instead driving thecorrection signal 46 towards zero. - Fig. 21 uses like reference numerals from Fig. 18 and illustrates a further embodiment. In Fig. 21, spectral
transfer function filter 302 has itsinput 316 supplied from the reference signal frominput transducer 10, rather than from the correction signal from the output ofmodel 40, Fig. 18, or the randomnoise signal source 208, Fig. 15. The output of spectraltransfer function filter 302 is supplied to filterinput 304 ofM copy 200, as above. The spectral transfer function provided byfilter 306 includes apeak detector 312 similar topeak detector 97, Fig. 7. The output ofinverse C model 306 is supplied to peakdetector 312 comparing the output ofinverse C model 306 with a desiredpeak value 314. When the output of inverseC model copy 306 rises above the level ofpeak value 314, the positive output ofsummer 312 at thecontrol input 308 controls the gain ofM copy 200 to in turn control the gain of the leak signal atfilter input 304 supplied through thespectral transfer function 302 from the reference signal frominput transducer 10. -
Filters - In further embodiments, performance of
model 40 is controlled according to fuzzy logic to control the signal sent tooutput transducer 14. Fuzzy logic control is known in the prior art, for example "Adaptive Fuzzy Systems", E. Cox, IEEE Spectrum February 1993, pages 27-31, and the references noted therein at page 31, lower half of right column. In the present invention, fuzzy logic is used to provide self-designing control architecture using fuzzy rules and/or to control the filter weights which update the model and/or to control the leak signal which degrades performance of the model. The fuzzy logic controller, having a given set of rules, is used for computing or setting the control architecture and/or the filter transfer function and/or the filter weights and/or the leak signal. - Fuzzy logic enables control of model performance in a practical way for a multiplicity of factors which interact in a complex way. The fuzzy logic control is based upon relative values and qualitative trends, without requiring exact equation relationships for multiple variables. For example, referring to Figs. 22-24, which are similar to Fig. 4, page 29 of the above noted Cox article, if the magnitude of the weights is small, Fig. 22, and the magnitude of the correction signal output of
model 40 is small, Fig. 23, then only a small leakage value is introduced, Fig. 24, to maximize the adaptive process. If the magnitude of the weights is large and the magnitude of the model output is large, then a large leak signal is introduced, to constrain model performance and minimize the adaptive process. In this example, the filter weights are monitored to provide a first input parameter, and the output of the model is monitored to provide a second input parameter. Fuzzy logic is applied to such first and second input parameters to provide an output parameter, Fig. 24, controlling the weights and/or leak and/or otherwise controlling performance of the model. The input parameters are fuzzified according to a fuzzy rule set, Figs. 22 and 23, a fuzzy leak is computed and then defuzzified using Fig. 24, to provide an output parameter to control performance of the model. - In another example, if the correction signal magnitude is nearing the capacity of
output transducer 14 such as a loudspeaker, then leakage is increased. If the magnitude of the correction signal is not nearing the capacity of the loudspeaker, then leakage can be decreased. - In another example, rate of change of the weights is monitored, and if the weights start increasing at a faster rate or at a rate above a given rate, then leakage is increased. If the weights are increasing rapidly and the magnitude of the correction signal is approaching the capacity of the loudspeaker, then leakage is increased very rapidly. If the weights start decreasing, then leakage is decreased.
- Inputs to the fuzzy logic include magnitude of the filter update weight, magnitude of the model output correction signal, rate of change of the filter weights, rate of change of the correction signal output by the model, magnitude of the
reference signal 42 input to the model, rate of change ofinput signal 42, magnitude of the error signal fromerror transducer 16, rate of change of the error signal, spectral characteristics of the reference signal and/or correction signal and/or error signal, capability ofoutput transducer 14, temperature, flow rates, environmental variables, fan speed in a duct application, whether or not such fan is running, thermostat settings, desired speed of adaptation, desired algorithm stability, system plant information, source information, etc. Fuzzy logic is a means of introducing human-like intuition into the controller adjustment process. Some factors may override others. For example, if the rate of increase of the weights is too large, then the algorithm may become unstable, and hence it is desirable to increase the leak or decrease the model update gain. If the capacity of aloudspeaker 14 is rapidly being reached, then leakage should be increased. - In further embodiments, the invention is applicable to multi-channel active acoustic attenuation systems, for example as shown in U.S. Patents 5,216,721 and 5,216,722, incorporated herein by reference.
- It is recognized that various equivalents, alternatives and modifications are possible within the scope of the appended claims.
Claims (57)
- An active adaptive control method comprising introducing a control signal from an output transducer to combine with a system input signal and yield a system output signal, sensing said system output signal with an error transducer providing an error signal, providing an adaptive filter model having a model input from a reference signal correlated to said system input signal, and an output outputting a correction signal to said output transducer to introduce said control signal, spectrally controlling performance of said model to maximize the signal sent to said output transducer at frequencies of interest.
- The method according to claim 1 comprising providing a spectral leak signal degrading performance of said model, and controlling said leak signal according to frequency.
- The method according to claim 2 comprising monitoring said correction signal and controlling said leak signal in response thereto.
- The method according to claim 3 comprising filtering said correction signal to provide said leak signal.
- The method according to claim 4 comprising spectrally processing said correction signal by supplying said correction signal through a frequency responsive spectral transfer function.
- The method according to claim 5 comprising modeling the transfer function between said output transducer and said error transducer with a second model, and wherein said spectral transfer function is a function of said second model.
- The method according to claim 5 comprising supplying said correction signal from said model output through said frequency responsive spectral transfer function to an error input of said model.
- The method according to claim 3 wherein said model output outputs said correction signal to said output transducer to introduce said control signal according to a weight update signal, and comprising combining said reference signal and said error signal to provide said weight update signal, and controlling said leak signal by controlling said weight update signal in response to said correction signal.
- The method according to claim 2 comprising providing a copy of said model, said copy having an output supplied to said error input of said model, said copy having a filter input receiving a spectrally processed input signal from a spectral transfer function, said copy having a control input controlling convergence gain of said copy and receiving a spectral control signal.
- The method according to claim 9 wherein said spectral control signal is responsive to said correction signal.
- The method according to claim 10 comprising providing a second spectral transfer function, and supplying said correction signal through said second spectral transfer function to said control input of said copy to provide said spectral control signal.
- The method according to claim 11 wherein said first and second spectral transfer functions are the same.
- The method according to claim 11 comprising modeling the transfer function between said output transducer and said error transducer with a second model, and wherein each of said first and second spectral transfer functions is a function of said second transfer function.
- The method according to claim 11 wherein said first and second spectral transfer functions are different.
- The method according to claim 2 comprising spectrally processing said leak signal such thatat frequencies where maximum power from said output transducer reaches said error transducer, the correction signal supplied to said output transducer is maximized, andat frequencies where minimum power from said output transducer reaches said error transducer, the correction signal supplied to said output transducer is minimized.
- The method according to claim 15 comprising modeling the transfer function between said output transducer and said error transducer with a second model, and spectrally processing said leak signal as a function of said second model.
- The method according to claim 16 wherein said function of said second model is the inverse of said second model.
- The method according to claim 17 wherein said first model is controlled by an update signal from said inverse of said second model such thatat frequencies where maximum power from said output transducer reaches said error transducer, said second model has a maximum transfer characteristic and said inverse of said second model has a minimum transfer characteristic minimizing leakage of said update signal, to enable maximum output of said first model, andat frequencies where minimum power from said output transducer reaches said error transducer, said second model has a minimum transfer characteristic and said inverse of said second model has a maximum transfer characteristic maximizing leakage of said update signal, to minimize the output of said first model.
- The method according to claim 2 comprising filtering said reference signal to provide said leak signal.
- The method according to claim 19 comprising spectrally processing said reference signal by supplying said reference signal through a frequency responsive spectral transfer function.
- The method according to claim 20 comprising modeling the transfer function between said output transducer and said error transducer with a second model, and wherein said spectral transfer function is a function of said second model.
- The method according to claim 19 comprising monitoring said correction signal and controlling said leak signal in response thereto.
- The method according to claim 1 comprising spectrally processing said leak signal by supplying said leak signal through a frequency responsive spectral transfer function.
- The method according to claim 23 wherein said leak signal is provided by said correction signal, and comprising supplying said correction signal through said frequency responsive spectral transfer function.
- The method according to claim 24 comprising supplying said correction signal from said model output through said frequency responsive spectral transfer function to said error input of said model.
- The method according to claim 25 comprising providing a copy of said model, and supplying said correction signal from said model output through said frequency responsive spectral transfer function and through said copy to said error input of said model.
- The method according to claim 23 comprising providing a copy of said model, and supplying said leak signal through said frequency responsive spectral transfer function and through said copy to said error input of said model.
- The method according to claim 27 wherein said copy has a filter input receiving the output of said frequency responsive spectral transfer function, said copy has a control input controlling convergence gain of said copy, and comprising providing a second frequency responsive spectral transfer function, and supplying said correction signal from said output of said model through said second frequency responsive spectral transfer function to said control input of said copy.
- The method according to claim 27 comprising supplying a random noise signal to the input of said spectral transfer function.
- The method according to claim 27 comprising supplying said correction signal from said output of said model to the input of said spectral transfer function.
- The method according to claim 1 comprising spectrally controlling performance of said model in response to said correction signal.
- The method according to claim 31 wherein said model outputs said correction signal to said output transducer to introduce said control signal according to a weight update signal, and comprising combining said reference signal and said error signal to provide said weight update signal, monitoring and spectrally sensing said correction signal and providing selective leakage of said weight update signal in response thereto to control performance of said model according to frequency, to optimize performance of said model in frequency ranges of interest.
- The method according to claim 31 comprising providing a copy of said model and supplying the output of said copy to said error input of said model and monitoring and spectrally sensing said correction signal and providing an input to said copy in response thereto to control performance of said model according to frequency, to optimize performance of said model in frequency ranges of interest.
- The method according to claim 2 comprising frequency weighting said leak signal to optimize performance of said model in frequency ranges of interest.
- The method according to claim 34 comprising frequency weighting said leak signal by spectrally processing said correction signal through a frequency responsive spectral transfer function.
- The method according to claim 34 comprising frequency weighting said leak signal by spectrally processing said correction signal through a frequency responsive spectral transfer function and through a copy of said model.
- The method according to claim 34 comprising frequency weighting said leak signal by spectrally processing a random noise signal through a frequency responsive spectral transfer function and through a copy of said model.
- The method according to claim 34 comprising frequency weighting said leak signal by spectrally processing a random noise signal through a first frequency responsive spectral transfer function and through a copy of said model, and spectrally processing said correction signal through a second frequency responsive spectral transfer function and through said copy of said model.
- The method according to claim 34 comprising frequency weighting said leak signal by spectrally processing said correction signal through a first frequency responsive spectral transfer function and through a copy of said model, and spectrally processing said correction signal through a second frequency responsive spectral transfer function and through said copy of said model.
- The method according to claim 1 comprising spectrally controlling performance of said model in response to said reference signal.
- The method according to claim 40 wherein said model outputs said correction signal to said output transducer to introduce said control signal according to a weight update signal, and comprising combining said reference signal and said error signal to provide said weight update signal, monitoring and spectrally sensing said reference signal and providing selective leakage of said weight update signal in response thereto to control performance of said model according to frequency, to optimize performance of said model in frequency ranges of interest.
- An active adaptive control system comprising an output transducer introducing a control signal to combine with a system input signal and yield a system output signal, an error transducer sensing said system output signal and providing an error signal, an adaptive filter model having a model input from a reference signal correlated to said system input signal, and an output outputting a correction signal to said output transducer to introduce said control signal, a spectral controller spectrally controlling performance of said model to maximize the signal sent to said output transducer at frequencies of interest.
- An active adaptive control method comprising introducing a control signal from an output transducer to combine with a system input signal and yield a system output signal, sensing said system output signal with an error transducer providing an error signal, providing an adaptive filter model having a model input from a reference signal correlated to said system input signal, and an output outputting a correction signal to said output transducer to introduce said control signal, controlling performance of said model according to fuzzy logic to control the signal sent to said output transducer.
- The method according to claim 43 comprising providing a leak signal degrading performance of said model, and controlling said leak signal according to said fuzzy logic.
- The method according to claim 43 comprising updating said model by a set of filter weights, and controlling said weights according to said fuzzy logic.
- The method according to claim 43 comprising fuzzifying a plurality of input parameters according to a fuzzy rule set, computing a fuzzy leak, and defuzzifying said leak to provide an output parameter to control performance of said model.
- The method according to claim 43 comprising updating said model by a set of filter weights, monitoring said weights to provide a first input parameter, monitoring the output of said model to provide a second input parameter, applying fuzzy logic to said first and second input parameters to provide an output parameter controlling said weights.
- The method according to claim 47 comprising providing a leak signal degrading performance of said model, and controlling said leak signal with said output parameter.
- An active adaptive control method comprising introducing a control signal from an output transducer to combine with a system input signal and yield a system output signal, fuzzifying a plurality of input parameters according to a fuzzy rule set and providing a control architecture therefrom outputting a correction signal to said output transducer to introduce said control signal.
- The method according to claim 49 comprising providing said control architecture by setting a filter transfer function.
- The method according to claim 50 comprising setting said filter transfer function by setting filter weights.
- The method according to claim 51 comprising providing a leak signal degrading performance of said filter transfer function, computing a fuzzy leak from said fuzzified input parameters and then defuzzifying said leak to provide an output parameter to control performance of said filter transfer function.
- An active adaptive control system comprising an output transducer introducing a control signal to combine with a system input signal and yield a system output signal, an error transducer sensing said system output signal and providing an error signal, an adaptive filter model having a model input from a reference signal correlated to said system input signal, and an output outputting a correction signal to said output transducer to introduce said control signal, fuzzy logic controlling performance of said model to control the signal sent to said output transducer.
- An active adaptive control method comprising introducing a control signal from an output transducer to combine with a system input signal and yield a system output signal, sensing said system output signal with an error transducer providing an error signal, providing an adaptive filter model having a model input from a reference signal correlated to said system input signal, and a model output outputting a correction signal to said output transducer to introduce said control signal according to a weight update signal, adaptively varying said weight update signal.
- The method according to claim 54 comprising adaptively varying said weight update signal as a function of said correction signal and supplying the adaptively varied weight update signal to said adaptive filter model to adapt said correction signal variably controlling leakage of said weight update signal.
- The method according to claim 55 comprising varying leakage of said weight update signal by multiplying a previous weight update value by a factor γ and adding the result to the product of said reference signal and said error signal, and varying γ as a function of said correction signal.
- The method according to claim 54 comprising adaptively varying said weight update signal as a function of said correction signal relative to a given threshold to alter convergence of said correction signal when the latter exceeds said threshold.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US517464 | 1995-08-21 | ||
US08/517,464 US5715320A (en) | 1995-08-21 | 1995-08-21 | Active adaptive selective control system |
Publications (3)
Publication Number | Publication Date |
---|---|
EP0759606A2 true EP0759606A2 (en) | 1997-02-26 |
EP0759606A3 EP0759606A3 (en) | 1998-12-30 |
EP0759606B1 EP0759606B1 (en) | 2003-03-26 |
Family
ID=24059910
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP96306121A Expired - Lifetime EP0759606B1 (en) | 1995-08-21 | 1996-08-21 | Active adaptive selective control system |
Country Status (4)
Country | Link |
---|---|
US (1) | US5715320A (en) |
EP (1) | EP0759606B1 (en) |
CA (1) | CA2183669A1 (en) |
DE (1) | DE69626902T2 (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999017275A1 (en) * | 1997-10-01 | 1999-04-08 | Digisonix, Inc. | Adaptive control system with efficiently constrained adaptation |
EP0973151A2 (en) * | 1998-07-16 | 2000-01-19 | Matsushita Electric Industrial Co., Ltd. | Noise control system |
WO2003015074A1 (en) * | 2001-08-08 | 2003-02-20 | Nanyang Technological University,Centre For Signal Processing. | Active noise control system with on-line secondary path modeling |
EP1414021A1 (en) * | 2002-10-21 | 2004-04-28 | Silentium Ltd. | Active acoustic noise reduction system |
EP1583075A1 (en) * | 2004-03-31 | 2005-10-05 | Nederlandse Organisatie voor toegepast-natuurwetenschappelijk Onderzoek TNO | System for actively reducing sound |
US7106866B2 (en) | 2000-04-06 | 2006-09-12 | Siemens Vdo Automotive, Inc. | Active noise cancellation stability solution |
US7317801B1 (en) | 1997-08-14 | 2008-01-08 | Silentium Ltd | Active acoustic noise reduction system |
US7853024B2 (en) | 1997-08-14 | 2010-12-14 | Silentium Ltd. | Active noise control system and method |
US8855329B2 (en) | 2007-01-22 | 2014-10-07 | Silentium Ltd. | Quiet fan incorporating active noise control (ANC) |
US9431001B2 (en) | 2011-05-11 | 2016-08-30 | Silentium Ltd. | Device, system and method of noise control |
US9928824B2 (en) | 2011-05-11 | 2018-03-27 | Silentium Ltd. | Apparatus, system and method of controlling noise within a noise-controlled volume |
Families Citing this family (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6850252B1 (en) | 1999-10-05 | 2005-02-01 | Steven M. Hoffberg | Intelligent electronic appliance system and method |
US5978489A (en) * | 1997-05-05 | 1999-11-02 | Oregon Graduate Institute Of Science And Technology | Multi-actuator system for active sound and vibration cancellation |
US6295364B1 (en) | 1998-03-30 | 2001-09-25 | Digisonix, Llc | Simplified communication system |
US20030040910A1 (en) * | 1999-12-09 | 2003-02-27 | Bruwer Frederick J. | Speech distribution system |
US8374875B2 (en) * | 2000-01-31 | 2013-02-12 | Intel Corporation | Providing programming information in response to spoken requests |
US20010046300A1 (en) * | 2000-04-17 | 2001-11-29 | Mclean Ian R. | Offline active control of automotive noise |
SG106582A1 (en) * | 2000-07-05 | 2004-10-29 | Univ Nanyang | Active noise control system with on-line secondary path modeling |
US20020039422A1 (en) * | 2000-09-20 | 2002-04-04 | Daly Paul D. | Driving mode for active noise cancellation |
US6559645B2 (en) | 2000-11-17 | 2003-05-06 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Detector apparatus and method |
US20020076058A1 (en) * | 2000-12-19 | 2002-06-20 | Astorino John Frank | Engine rotation reference signal for noise attenuation |
US6594368B2 (en) | 2001-02-21 | 2003-07-15 | Digisonix, Llc | DVE system with dynamic range processing |
US6549629B2 (en) | 2001-02-21 | 2003-04-15 | Digisonix Llc | DVE system with normalized selection |
US6665411B2 (en) | 2001-02-21 | 2003-12-16 | Digisonix Llc | DVE system with instability detection |
EP2261759B8 (en) * | 2001-04-20 | 2012-12-05 | Honda Giken Kogyo Kabushiki Kaisha | Plant control apparatus |
US20030112981A1 (en) * | 2001-12-17 | 2003-06-19 | Siemens Vdo Automotive, Inc. | Active noise control with on-line-filtered C modeling |
FR2845780A1 (en) * | 2002-10-11 | 2004-04-16 | St Microelectronics Sa | METHOD FOR CONTROLLING A DYNAMIC SYSTEM, USING A FOULD LOGIC MODEL OF AT LEAST ONE REVERSE TRANSFER FUNCTION OF THIS SYSTEM |
DE602005026178D1 (en) * | 2004-11-09 | 2011-03-10 | Truveon Corp | BUILDING CLIMATE CONTROL PROCESS AND SYSTEM |
EP1912468B1 (en) * | 2005-07-29 | 2013-08-14 | Panasonic Corporation | Loudspeaker device |
US8194873B2 (en) * | 2006-06-26 | 2012-06-05 | Davis Pan | Active noise reduction adaptive filter leakage adjusting |
US20070297619A1 (en) * | 2006-06-26 | 2007-12-27 | Bose Corporation*Ewc* | Active noise reduction engine speed determining |
US8204242B2 (en) * | 2008-02-29 | 2012-06-19 | Bose Corporation | Active noise reduction adaptive filter leakage adjusting |
US8306240B2 (en) * | 2008-10-20 | 2012-11-06 | Bose Corporation | Active noise reduction adaptive filter adaptation rate adjusting |
US8355512B2 (en) * | 2008-10-20 | 2013-01-15 | Bose Corporation | Active noise reduction adaptive filter leakage adjusting |
US8135140B2 (en) * | 2008-11-20 | 2012-03-13 | Harman International Industries, Incorporated | System for active noise control with audio signal compensation |
US9020158B2 (en) * | 2008-11-20 | 2015-04-28 | Harman International Industries, Incorporated | Quiet zone control system |
US8718289B2 (en) * | 2009-01-12 | 2014-05-06 | Harman International Industries, Incorporated | System for active noise control with parallel adaptive filter configuration |
US8019090B1 (en) * | 2009-02-12 | 2011-09-13 | United States Of America As Represented By The Secretary Of The Navy | Active feedforward disturbance control system |
US8189799B2 (en) | 2009-04-09 | 2012-05-29 | Harman International Industries, Incorporated | System for active noise control based on audio system output |
US8199924B2 (en) * | 2009-04-17 | 2012-06-12 | Harman International Industries, Incorporated | System for active noise control with an infinite impulse response filter |
WO2010127276A1 (en) * | 2009-05-01 | 2010-11-04 | Bose Corporation | Multi-element electroacoustical transducing |
US8077873B2 (en) * | 2009-05-14 | 2011-12-13 | Harman International Industries, Incorporated | System for active noise control with adaptive speaker selection |
DE102009055889B4 (en) * | 2009-11-26 | 2015-02-26 | Austriamicrosystems Ag | Signal processing circuit for generating a loudspeaker signal and signal processing method |
WO2013063659A1 (en) * | 2011-11-03 | 2013-05-10 | Minelab Electronics Pty Limited | Improved metal detection method |
US10891935B2 (en) * | 2019-05-07 | 2021-01-12 | Harman International Industries, Incorporated | In-vehicle noise cancellation adaptive filter divergence control |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4677676A (en) * | 1986-02-11 | 1987-06-30 | Nelson Industries, Inc. | Active attenuation system with on-line modeling of speaker, error path and feedback pack |
EP0340974A2 (en) * | 1988-05-04 | 1989-11-08 | Nelson Industries, Inc. | Active acoustic attenuation system with differential filtering |
US5018202A (en) * | 1988-09-05 | 1991-05-21 | Hitachi Plant Engineering & Construction Co., Ltd. | Electronic noise attenuation system |
EP0471290A2 (en) * | 1990-08-16 | 1992-02-19 | Hughes Aircraft Company | Active adaptive noise canceller without training mode |
EP0517525A2 (en) * | 1991-06-06 | 1992-12-09 | Matsushita Electric Industrial Co., Ltd. | Noise suppressor |
WO1993015501A1 (en) * | 1992-01-28 | 1993-08-05 | Active Noise And Vibration Technologies Inc. | Active noise cancellation |
EP0661807A2 (en) * | 1993-12-14 | 1995-07-05 | DIGISONIX, Inc. | Active adaptive control system with spectral leak |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4791390A (en) * | 1982-07-01 | 1988-12-13 | Sperry Corporation | MSE variable step adaptive filter |
US4677677A (en) * | 1985-09-19 | 1987-06-30 | Nelson Industries Inc. | Active sound attenuation system with on-line adaptive feedback cancellation |
US4689821A (en) * | 1985-09-23 | 1987-08-25 | Lockheed Corporation | Active noise control system |
US4736431A (en) * | 1986-10-23 | 1988-04-05 | Nelson Industries, Inc. | Active attenuation system with increased dynamic range |
US5033082A (en) * | 1989-07-31 | 1991-07-16 | Nelson Industries, Inc. | Communication system with active noise cancellation |
US4987598A (en) * | 1990-05-03 | 1991-01-22 | Nelson Industries | Active acoustic attenuation system with overall modeling |
US5386472A (en) * | 1990-08-10 | 1995-01-31 | General Motors Corporation | Active noise control system |
US5146505A (en) * | 1990-10-04 | 1992-09-08 | General Motors Corporation | Method for actively attenuating engine generated noise |
US5396561A (en) * | 1990-11-14 | 1995-03-07 | Nelson Industries, Inc. | Active acoustic attenuation and spectral shaping system |
US5172416A (en) * | 1990-11-14 | 1992-12-15 | Nelson Industries, Inc. | Active attenuation system with specified output acoustic wave |
US5216721A (en) * | 1991-04-25 | 1993-06-01 | Nelson Industries, Inc. | Multi-channel active acoustic attenuation system |
JP3089082B2 (en) * | 1991-07-10 | 2000-09-18 | シャープ株式会社 | Adaptive digital filter |
US5216722A (en) * | 1991-11-15 | 1993-06-01 | Nelson Industries, Inc. | Multi-channel active attenuation system with error signal inputs |
US5206911A (en) * | 1992-02-11 | 1993-04-27 | Nelson Industries, Inc. | Correlated active attenuation system with error and correction signal input |
US5216692A (en) * | 1992-03-31 | 1993-06-01 | Motorola, Inc. | Method and apparatus for adjusting a power control threshold in a communication system |
DE69328851T2 (en) * | 1992-07-07 | 2000-11-16 | Sharp K.K., Osaka | Active control device with an adaptive digital filter |
US5278913A (en) * | 1992-07-28 | 1994-01-11 | Nelson Industries, Inc. | Active acoustic attenuation system with power limiting |
US5390255A (en) * | 1992-09-29 | 1995-02-14 | Nelson Industries, Inc. | Active acoustic attenuation system with error and model copy input |
US5386477A (en) * | 1993-02-11 | 1995-01-31 | Digisonix, Inc. | Active acoustic control system matching model reference |
-
1995
- 1995-08-21 US US08/517,464 patent/US5715320A/en not_active Expired - Lifetime
-
1996
- 1996-08-19 CA CA002183669A patent/CA2183669A1/en not_active Abandoned
- 1996-08-21 EP EP96306121A patent/EP0759606B1/en not_active Expired - Lifetime
- 1996-08-21 DE DE69626902T patent/DE69626902T2/en not_active Expired - Lifetime
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4677676A (en) * | 1986-02-11 | 1987-06-30 | Nelson Industries, Inc. | Active attenuation system with on-line modeling of speaker, error path and feedback pack |
EP0340974A2 (en) * | 1988-05-04 | 1989-11-08 | Nelson Industries, Inc. | Active acoustic attenuation system with differential filtering |
US5018202A (en) * | 1988-09-05 | 1991-05-21 | Hitachi Plant Engineering & Construction Co., Ltd. | Electronic noise attenuation system |
EP0471290A2 (en) * | 1990-08-16 | 1992-02-19 | Hughes Aircraft Company | Active adaptive noise canceller without training mode |
EP0517525A2 (en) * | 1991-06-06 | 1992-12-09 | Matsushita Electric Industrial Co., Ltd. | Noise suppressor |
WO1993015501A1 (en) * | 1992-01-28 | 1993-08-05 | Active Noise And Vibration Technologies Inc. | Active noise cancellation |
EP0661807A2 (en) * | 1993-12-14 | 1995-07-05 | DIGISONIX, Inc. | Active adaptive control system with spectral leak |
Non-Patent Citations (1)
Title |
---|
KIPERSZTOK O: "ACTIVE CONTROL OF BROADBAND NOISE USING FUZZY LOGIC" PROCEEDINGS OF THE INTERNATIONAL CONFERENCE ON FUZZY SYSTEMS, SAN FRANCISCO, MAR. 28 - APR. 1, 1993, vol. 2, no. CONF. 2, 28 March 1993, pages 906-911, XP000371528 INSTITUTE OF ELECTRICAL AND ELECTRONICS ENGINEERS * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7317801B1 (en) | 1997-08-14 | 2008-01-08 | Silentium Ltd | Active acoustic noise reduction system |
US8630424B2 (en) | 1997-08-14 | 2014-01-14 | Silentium Ltd. | Active noise control system and method |
US7853024B2 (en) | 1997-08-14 | 2010-12-14 | Silentium Ltd. | Active noise control system and method |
US6094601A (en) * | 1997-10-01 | 2000-07-25 | Digisonix, Inc. | Adaptive control system with efficiently constrained adaptation |
AU740931B2 (en) * | 1997-10-01 | 2001-11-15 | Digisonix, Inc. | Adaptive control system with efficiently constrained adaptation |
WO1999017275A1 (en) * | 1997-10-01 | 1999-04-08 | Digisonix, Inc. | Adaptive control system with efficiently constrained adaptation |
EP0973151A3 (en) * | 1998-07-16 | 2003-01-02 | Matsushita Electric Industrial Co., Ltd. | Noise control system |
EP0973151A2 (en) * | 1998-07-16 | 2000-01-19 | Matsushita Electric Industrial Co., Ltd. | Noise control system |
US7106866B2 (en) | 2000-04-06 | 2006-09-12 | Siemens Vdo Automotive, Inc. | Active noise cancellation stability solution |
WO2003015074A1 (en) * | 2001-08-08 | 2003-02-20 | Nanyang Technological University,Centre For Signal Processing. | Active noise control system with on-line secondary path modeling |
EP1414021A1 (en) * | 2002-10-21 | 2004-04-28 | Silentium Ltd. | Active acoustic noise reduction system |
EP1583075A1 (en) * | 2004-03-31 | 2005-10-05 | Nederlandse Organisatie voor toegepast-natuurwetenschappelijk Onderzoek TNO | System for actively reducing sound |
WO2005096269A1 (en) * | 2004-03-31 | 2005-10-13 | Nederlandse Organisatie Voor Toegepast- Natuurwetenschappelijk Onderzoek Tno | System for actively reducing sound |
US8855329B2 (en) | 2007-01-22 | 2014-10-07 | Silentium Ltd. | Quiet fan incorporating active noise control (ANC) |
US9431001B2 (en) | 2011-05-11 | 2016-08-30 | Silentium Ltd. | Device, system and method of noise control |
US9928824B2 (en) | 2011-05-11 | 2018-03-27 | Silentium Ltd. | Apparatus, system and method of controlling noise within a noise-controlled volume |
Also Published As
Publication number | Publication date |
---|---|
CA2183669A1 (en) | 1997-02-22 |
DE69626902T2 (en) | 2003-08-21 |
DE69626902D1 (en) | 2003-04-30 |
US5715320A (en) | 1998-02-03 |
EP0759606B1 (en) | 2003-03-26 |
EP0759606A3 (en) | 1998-12-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0759606B1 (en) | Active adaptive selective control system | |
US5586190A (en) | Active adaptive control system with weight update selective leakage | |
EP0340974B1 (en) | Active acoustic attenuation system with differential filtering | |
US5627896A (en) | Active control of noise and vibration | |
US5337366A (en) | Active control apparatus using adaptive digital filter | |
US5710822A (en) | Frequency selective active adaptive control system | |
US5278913A (en) | Active acoustic attenuation system with power limiting | |
EP0455479B1 (en) | Active acoustic attenuation system with overall modeling | |
US5475761A (en) | Adaptive feedforward and feedback control system | |
US4677677A (en) | Active sound attenuation system with on-line adaptive feedback cancellation | |
US5410606A (en) | Noise canceling method | |
EP0555585B1 (en) | Correlated active attenuation system with error and correction signal input | |
WO2003015074A1 (en) | Active noise control system with on-line secondary path modeling | |
US5602929A (en) | Fast adapting control system and method | |
LJ | II. recursive algorithms for active noise control | |
EP0661807B1 (en) | Active adaptive control system with spectral leak | |
EP3844741B1 (en) | Systems and methods for noise-cancellation with shaping and weighting filters | |
EP0492680B1 (en) | Method and apparatus for attenuating noise | |
JP3397245B2 (en) | Electronic silencer | |
Kohno et al. | Direct adaptive active noise control algorithms in case of uncertain secondary path dynamics | |
Phooi et al. | Nonlinear active noise control using Lyapunov theory and RBF network | |
JP3697702B2 (en) | Vehicle vibration control apparatus and vibration control method | |
JPH06124092A (en) | Noise eliminating device for air conditioner | |
Nam et al. | Adaptive active attenuation of noise using multiple model approaches | |
Davari et al. | A New Feedback ANC System Approach |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
AK | Designated contracting states |
Kind code of ref document: A2 Designated state(s): DE FR GB IT NL SE |
|
PUAL | Search report despatched |
Free format text: ORIGINAL CODE: 0009013 |
|
AK | Designated contracting states |
Kind code of ref document: A3 Designated state(s): DE FR GB IT NL SE |
|
17P | Request for examination filed |
Effective date: 19990628 |
|
RIN1 | Information on inventor provided before grant (corrected) |
Inventor name: ERIKSSON, LARRY J. Inventor name: ALLIE, MARK C. |
|
17Q | First examination report despatched |
Effective date: 20010601 |
|
GRAG | Despatch of communication of intention to grant |
Free format text: ORIGINAL CODE: EPIDOS AGRA |
|
GRAG | Despatch of communication of intention to grant |
Free format text: ORIGINAL CODE: EPIDOS AGRA |
|
GRAG | Despatch of communication of intention to grant |
Free format text: ORIGINAL CODE: EPIDOS AGRA |
|
GRAH | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOS IGRA |
|
GRAH | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOS IGRA |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
AK | Designated contracting states |
Designated state(s): DE FR GB IT NL SE |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: IT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT;WARNING: LAPSES OF ITALIAN PATENTS WITH EFFECTIVE DATE BEFORE 2007 MAY HAVE OCCURRED AT ANY TIME BEFORE 2007. THE CORRECT EFFECTIVE DATE MAY BE DIFFERENT FROM THE ONE RECORDED. Effective date: 20030326 |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D |
|
REF | Corresponds to: |
Ref document number: 69626902 Country of ref document: DE Date of ref document: 20030430 Kind code of ref document: P |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: SE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20030626 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: GB Payment date: 20030807 Year of fee payment: 8 |
|
ET | Fr: translation filed | ||
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: NL Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20040301 |
|
26N | No opposition filed |
Effective date: 20031230 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: FR Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20040430 |
|
NLV4 | Nl: lapsed or anulled due to non-payment of the annual fee |
Effective date: 20040301 |
|
REG | Reference to a national code |
Ref country code: FR Ref legal event code: ST |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: GB Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20040821 |
|
GBPC | Gb: european patent ceased through non-payment of renewal fee |
Effective date: 20040821 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: DE Payment date: 20151029 Year of fee payment: 20 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R071 Ref document number: 69626902 Country of ref document: DE |