[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

EP0285632B1 - Active vibration control - Google Patents

Active vibration control Download PDF

Info

Publication number
EP0285632B1
EP0285632B1 EP87906587A EP87906587A EP0285632B1 EP 0285632 B1 EP0285632 B1 EP 0285632B1 EP 87906587 A EP87906587 A EP 87906587A EP 87906587 A EP87906587 A EP 87906587A EP 0285632 B1 EP0285632 B1 EP 0285632B1
Authority
EP
European Patent Office
Prior art keywords
filter
control system
vibration control
signal
active vibration
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.)
Expired - Lifetime
Application number
EP87906587A
Other languages
German (de)
French (fr)
Other versions
EP0285632A1 (en
Inventor
Stephen John Elliott
Philip Arthur Nelson
Ian Mcgregor Stothers
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Adaptive Control Ltd
Original Assignee
Adaptive Control Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=26291394&utm_source=google_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=EP0285632(B1) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Priority claimed from GB08624053A external-priority patent/GB2203016A/en
Application filed by Adaptive Control Ltd filed Critical Adaptive Control Ltd
Publication of EP0285632A1 publication Critical patent/EP0285632A1/en
Application granted granted Critical
Publication of EP0285632B1 publication Critical patent/EP0285632B1/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods 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/16Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178Methods 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
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/10Applications
    • G10K2210/128Vehicles
    • G10K2210/1282Automobiles
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/30Means
    • G10K2210/301Computational
    • G10K2210/3023Estimation of noise, e.g. on error signals
    • G10K2210/30232Transfer functions, e.g. impulse response
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/30Means
    • G10K2210/301Computational
    • G10K2210/3032Harmonics or sub-harmonics
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/30Means
    • G10K2210/301Computational
    • G10K2210/3046Multiple acoustic inputs, multiple acoustic outputs
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/30Means
    • G10K2210/301Computational
    • G10K2210/3048Pretraining, e.g. to identify transfer functions
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/30Means
    • G10K2210/301Computational
    • G10K2210/3049Random noise used, e.g. in model identification
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/30Means
    • G10K2210/321Physical
    • G10K2210/3212Actuator details, e.g. composition or microstructure

Definitions

  • the invention relates to active vibration control.
  • vibration includes sound or noise, and the invention is particularly concerned with active noise control.
  • an active vibration control system for reducing vibration generated by a primary source, is defined in Claim 1. It is characterised in that at least one reference signal containing selected harmonics of the said primary source vibration is supplied to means driving a plurality of secondary vibration sources, such that vibration energy detected by sensor means operable to sense the vibration field established by the primary and secondary sources is reduced.
  • control system is operable in accordance with an algorithm which adjusts the outputs from the secondary sources so as to substantially reduce a cost function on a time scale comparable with the delays associated with the propagation of vibration from the secondary sources to the sensor means.
  • the present invention is particularly concerned with an active noise reduction system which can control the sound throughout an enclosure of a car, or at one or a number of "quiet zones" within it, and which can quickly adapt to changes in the excitation of the sound field due to changes in, for example, engine load or speed.
  • a signal related to the engine crankshaft rotation rate for example a signal emitted by the engine ignition system, is used to generate a reference signal containing a number of sinusoids at harmonics (or subharmonics) of the engine crankshaft rotation frequency. These are known as engine order frequencies. These sinusoids may be obtained using a variety of methods outlined below. Rather than attempt to control all harmonics, only a selected set of engine order frequencies may, in accordance with the invention, be generated as the reference signal. For example, only the firing frequency (second engine order in a four cylinder car) and its second harmonic (fourth engine order) are used if the spectrum of the sound in the car is dominated by these components.
  • a signal containing all the engine order frequencies may be fed to a band pass filter which isolates only the particular frequency or frequencies exciting a particular resonance in the enclosure interior which could cause a "boom" to be excited.
  • the advantage of reducing the number of frequencies fed to the filter is that an adaptive filter having fewer coefficients than would otherwise be the case can be used. This makes implementation more efficient and allows faster adaption time. The faster adaption time is particularly important in automotive applications, for example, in which the active control system has to adapt sufficiently quickly to track changes in engine speed which may occur on a very short timescale.
  • reference signals may be obtained from transducers mounted on a road wheel hub or the suspension system of the car. Such reference signals would contain the harmonics of road wheel rotation or road noise. Transducers placed outside a vehicle may provide reference signals representative of wind noise. If the reference signals are periodic (deterministic), control at individual harmonics can be exercised as described hereinafter.
  • the secondary sources may be loudspeakers used as the low frequency drives of a car audio system.
  • a signal is obtained from the primary vibration source which contains components at all harmonics prevalent in the sound in a vehicle powered by an engine.
  • This signal is filtered so as only to leave the most important or dominant harmonics. Filtering is carried out by a filter whose centre frequency can be controlled by an external signal in such a manner that the critical filter frequencies have a constant ratio compared with the engine crankshaft rotation rate. This can be achieved by using, for example, charge coupled devices whose switching frequency is locked to the crankshaft rotation frequency, but can, also be implemented as a program running on a microprocessor, as described hereinafter.
  • a primary source signal rich in harmonics is filtered by a band pass filter, having a centre frequency fixed at that of a pronounced "boom" in the car enclosure and a characteristic such that the reference signal only contains the harmonic(s) which are particularly exciting the boom.
  • This may be extended such that the filter contains a number of resonances at a number of boom frequencies of the car, or even such that the filter has a frequency response which models the acoustic response of the car interior to the primary excitation.
  • the input signal to the filter may be a signal from the engine containing all important harmonics and may be in the form of a pulse train.
  • This method may be accomplished by using phase lock loops to generate sinusoidal signals, with frequencies bearing an integer relationship to a square wave signal from the engine, which frequencies are then added together to form the reference signal.
  • the signal derived from the engine can be used to control a number of tunable oscillators, each producing a sinusoid at a selected harmonic frequency.
  • ⁇ (n) is a unit sample sequence to initiate the oscillators formed by the difference equations above.
  • the second and fourth harmonic may be generated for example.
  • sample frequency (f s ) is derived from the counter frequency (f c ) by a frequency division circuit for example, so that the ratio f s /f c is exactly an integer number.
  • An alternative difference equation which may be used to implement the digital oscillator has the form of a series approximation to a trigonometric function, for example :-
  • the series approximation above can be used for y(n) in the range - ⁇ /2 ⁇ y(n) ⁇ /2 . For values of y(n) outside this range the symmetry properties of the cosine waveform are utilised, until y(n)> ⁇ .
  • the calculations of the coefficients used in the difference equations forming the digital oscillators, together with the difference equations themselves may be implemented on a dedicated processor, or may form part of the program which also implements a controller that generates the outputs used to drive the secondary sources from the reference signals described above.
  • the controller is designed to be adaptive so as to quickly track changes in engine speed and load.
  • the outputs of the secondary sources are adaptively controlled so that some measurable cost function is minimised.
  • This cost function would typically be the sum of the mean square outputs from a number of microphones in the enclosed space.
  • the controller can be implemented as a digital adaptive FIR filter, using the basic update algorithm described by S. Elliott and P. Nelson in “Electronics Letters", at pp. 979-981, 1985. A number of additions must be made to this basic algorithm, however, to enable it to work quickly and efficiently in this particular example.
  • the basic algorithm referred to hereinafter as the stochastic gradient algorithm, is presented below, in order to highlight the necessary alterations.
  • each of these coefficients should be adjusted at every sample according to the equation where ⁇ is a convergence coefficient, e l (n) is the sampled output from the l'th sensor and r l m (n) is a sequence formed by filtering the reference signal discussed above ( x(n), say) with a digital filter which models the response of the l'th sensor to excitation of the m'th secondary source.
  • is a convergence coefficient
  • e l (n) is the sampled output from the l'th sensor
  • r l m (n) is a sequence formed by filtering the reference signal discussed above ( x(n), say) with a digital filter which models the response of the l'th sensor to excitation of the m'th secondary source.
  • the digital filters must model the relevant response over a range of frequencies, governed by the frequency range which the active system is attempting to control. It has been found that under certain circumstances a digital filter need only model the overall delay in the response to ensure the stability of an adaptive filter. It is more common, however, to have the digital filters incorporate a delay and then some reverberant response. This may be implemented using either digital FIR or IIR filters whose coefficients are adjusted adaptively during an initialisation phase, so as to accurately match the desired responses. It is also possible to continue this initial adaption process during the operation of the active control system by feeding training signals to each secondary source which are suitably uncorrelated with each other and with the primary excitation.
  • this change may be detected with mechanical transducers and the information used to switch between a cariety of filters modelling the response of the enclosure under a variety of conditions.
  • a number of other adaptive algorithms may also be implemented to adjust the coefficients of the digital filters in the controller driving the secondary sources. These alternative algorithms are best described in matrix form.
  • y m (n) may be obtained by passing this reference signal through a digital filter whose i'th coefficient is w mi (n) at the n'th sample, so that
  • the sampled output from the l'th error sensor, e l (n) is equal to the sum of the contributions from the primary source, d l (n), and each of the secondary sources.
  • an enclosure 10 which is the interior of the passenger or driver compartment of an internal combustion engine driven vehicle, in this example, a motor car 100, is represented schematically together with an active sound control system 1 according to the invention.
  • the system 1 employs two secondary sound sources 11, comprising two low frequency loudspeakers of a stereo audio system fitted to the car, and three acoustic sensors, comprising microphones 12.
  • the loudspeakers 11 are driven by a controller circuit 13 which comprises a pair of adaptive filters 14.
  • Each adaptive filter 14 drives a respective one of the loudspeakers 11 with an output signal 3 which the filter 14 produces as a result of its action on a reference signal 4 supplied thereto by a reference signal generator 15.
  • the reference signal 4 is generated by the generator 15 from an input signal 16 which is periodic at the crankshaft rotation rate of the internal combustion engine 2.
  • the signal generator 15 may comprise a tracking filter.
  • the purpose of the outputs from the loudspeakers 11 driven by the controller 13 is to reduce the sound vibration field established by the primary and secondary sources, experienced within the enclosure 10. Since the primary source (engine 2) of the noise to be reduced is periodic, the reference signal 4 generated by the generator 15 is, in accordance with the invention, arranged to contain one or more sinusoidal components at harmonics (or subharmonics) of the crankshaft rotation rate of the engine 2.
  • the adaptive filters 14 are adjusted automatically by output signals 5 from the sensor-microphones 12, corresponding adjustment being made simultaneously to the outputs of the loudspeakers 11 so as to substantially minimise a cost function on a time scale comparable with the delays associated with the propagation of sound vibrations from the loudspeakers 11 to the microphones 12.
  • the cost function may comprise the sum of the mean square outputs of the microphones 12.
  • the control system 1 does not employ a stored solution. Instead, it makes use of a plurality of closed loops, each loop comprising a microphone 12, the controller 13, and a loudspeaker 11, whereby signals from the microphone 12 are used to adapt the filters 14 controlling the loudspeaker 11, which has an influence on the output of the microphone as a result of acoustic response within the enclosure 10.
  • the loudspeakers 11 and the microphones 12 are distributed in the enclosure 10 in spaced relationship. The distribution, which varies from car to car, is adjusted to get substantial sound reductions throughout the enclosure 10.
  • system 1 employs as many closed loops as the number of sensors (12) multiplied by the number of secondary sources (11).
  • system 1 employs more sensors (12) than secondary sources (11), whereby a controlled reduction of primary source vibration is achieved.
  • sensors (12) employ more sensors (12) than secondary sources (11), whereby a controlled reduction of primary source vibration is achieved.
  • the reference signal 4 contains harmonics of the input signal 16, which signal is periodic at the engine crankshaft rotation rate
  • the reference signal 4 contains engine order frequencies.
  • the signal generator 15 is arranged to select engine order frequencies that ensure that the sound produced by the loudspeakers 11 is of the same frequency or frequencies as the sound produced in the enclosure 10 by the engine 2, even during changes in engine conditions such as load or speed.
  • the number of engine order frequencies in the reference signal 4 is restricted so that the adaptive filters 14 have a relatively small number of coefficients and can therefore adapt quickly.
  • the input signal 16 can be obtained from another moving part of the engine or part of the ignition circuitry; for example.
  • Figure 2(a) illustrates graphically the response of the reference signal generator 15 when in the form of a tracking filter. That is to say, a filter having a centre frequency so controlled that the filter output frequencies have a constant ratio to the dominant input frequencies, so that in Figure 2(a), the frequency f o is N x (engine crankshaft rotation rate), and the reference signal 4 ( Figure 1) contains only the first N harmonics of the engine rotation rate, where N is an integer.
  • the filter input signal is a voltage pulse train as represented by Figure 2(b), where 8T is the periodic time of the engine rotation rate, the first eight harmonics of the engine rotation rate are present in the reference signal 4.
  • the spectrum, by Fourier analysis, of the reference signal 4 is then illustrated by Figure 2(c), in which A is amplitude. With the response of Figure 2(a), only the first six harmonics would be usable.
  • the tracing filter comprising the signal generator, can be in the form of charge coupled devices having a switching frequency locked to the engine crankshaft rotation rate.
  • FIG. 3 An alternative form of reference signal generator, which employs a plurality of tracking band pass filters is illustrated by Figure 3 in which the input signal 16, a square wave at, for example, 128 times the engine drive shaft rotation rate, is divided first by 32 and then by 2. Division by 32 is achieved by a divider 6 which produces a square wave signal 7 at four times the engine crankshaft rotation rate, which signal is supplied to a bandpass filter 17. The filter 17 has a centre frequency f4 which is arranged to track the fundamental frequency of the square wave signal 7 supplied thereto. The further division by 2 is achieved by a divider 8 which produces a square wave signal 9 at twice the engine drive shaft rotation rate.
  • the signal 9 is then supplied to a bandpass filter 18 having a centre frequency f2 and which is arranged to track the fundamental frequency of the square wave signal 9 supplied thereto.
  • the bandpass filters 17 and 18 produce respectively sinusoidal output signals 7 a , 9 a at f2 and f4 which are linearly summed at an adder 19 to produce the required reference signal 4.
  • reference signal generator may comprise a fixed frequency filtering circuit which selects harmonics and/or subharmonics from an input signal rich in the harmonics of the engine crankshaft rotation rate or firing rate.
  • the filtering circuit may comprise a bandpass filter having a centre frequency fixed at the frequency of a pronounced resonance excited in the enclosure 10 ( Figure 1) by the engine 2 or other primary source of vibration.
  • a bandpass filter may be arranged to have a response that models the acoustic response of a motor car passenger compartment to the engine.
  • a further form of reference signal generator may comprise a plurality of phase lock loops used to generate sinusoidal signals having respective frequencies with integer relationships to a square wave input signal from the engine 2 or other primary vibration source. The sinusoidal signals may then be added together to form the required reference signal. Thus a reference signal 4 comprising specific harmonics and/or harmonics locked to the primary source fundamental, such as engine crankshaft rotation rate, is generated.
  • FIG. 4 An alternative generator for such a reference signal is illustrated in Figure 4 in which a square wave 20 at the primary source fundamental is used to control a plurality of tunable digital oscillators 25, 26 each producing a sinusoidal signal at a chosen harmonic or subharmonic frequency to be added at the adder 19 which produce the reference signal 4 by simple addition of the sinusoidal signals.
  • the square wave signal 20 which is at an engine crankshaft rotation rate, is supplied to a bistable circuit 21 which divides the rate signal by two and thereby produces a pulse train signal 20 a in which the duration of each pulse is equal to the prevailing periodic time of the square wave signal 20.
  • This periodic time is then measured by a counter 22 which is enabled throughout the duration of each positive pulse from the bistable circuit 21 and counts clock pulses supplied by a clock pulse generator 23.
  • the clock pulses are generated at a fixed, suitably high rate f c .
  • the contents of the counter 22 are read at the end of each positive pulse from the bistable circuit 21 by a trigonometric function generator 24.
  • the digital sinusoidal outputs from the two oscillators 25 and 26 are superposed by a digital adder 191 which supplies the reference signal 4 as a digital signal.
  • the trigonometric function generator 24, oscillators 25 and 26 and adder 191 can be implemented by a microprocessor with a suitable program.
  • Figure 5 represents in block form an active sound control system 30 for reducing the level of engine generated noise in the passenger compartment of a motor car.
  • the car is provided with an ignition circuit including a low tension coil 31 from which a voltage signal 32 at the firing rate of the engine is taken and supplied to a waveform shaper 33 which in response thereto produces a pulse train at the engine firing rate.
  • the engine firing rate is twice the engine crankshaft rotation rate f o .
  • the shaper 33 provides a signal having a fundamental frequency which is a single harmonic, (2f0), of the crankshaft rate.
  • a reference signal generator is provided in the form of a proprietory tracking filter 34, manufactured by Bruel and Kjaer under type number 1623.
  • the tracking filter 34 receives the output of the shaper 33 as an input signal and as a trigger signal and produces a sinusoidal output signal at the selected harmonic 2f o .
  • This sinusoidal signal is sampled with an analog to digital converter 35 to produce a reference sequence x(n) of digitised samples which are supplied as data to a processor and memory unit 36.
  • the loudspeakers 371 and 372 are mounted within the motor car passenger compartment 10 (not shown in Figure 5) are two loudspeakers 371 and 372, which are in positions normally used for car stereo reproduction.
  • the loudspeakers 371, 372 are driven by a multiplexer 38 through respective low pass filters 39 and output amplifiers 40.
  • the filters 39 have a cut off frequency of 460 Hz and are provided to prevent aliasing.
  • the multiplexer 38 which contains sample and hold circuits for each output, is controlled by the processor and memory unit 36, through a control line 55, and receives a single input signal 57 from a digital to analog converter 41.
  • the purpose of the loudspeakers 371 and 372 is to generate, in the passenger compartment, audio waves that will cancel those set up directly by mechanical transmission from the engine to the compartment.
  • the digital to analog converter 41 is supplied by the processor and memory unit 36 with output data 58 which consists of two interleaved sequences of digitised samples y1(n) and y2(n).
  • the data 58 is converted by the converter 41 into interleaved sequences of analog samples and separated into respective sequences by the multiplexer 38 for application to the low pass filters 39.
  • the loudspeaker 371 is driven by the sequence y1(n) and the loudspeaker 372 is driven by the sequence y2(n).
  • each sequence of data 58 is represented by the expression y m (n), so that in this example m may be 1 or 2.
  • error signals are picked up from the passenger compartment and utilised by the processor and memory unit 36.
  • Acoustic error signals if present, are sensed by four microphones 421, 422, 423 and 424, which are placed respectively either side of a driver headrest and a passenger headrest, there being only two seats in the compartment in the present example.
  • the electrical outputs from the microphones 421 etc. are respectively amplified by amplifiers 43 and passed through low pass filters 44 to a four-input multiplexer 45 which supplies a single analog output to an analog to digital converter 46.
  • the filters 44 are provided to prevent aliasing and have a cut-off frequency of 460 Hz.
  • the multiplexer 45 is controlled by the processor unit 36 by way of control line 56.
  • the multiplexer 45 and the converter 46 convert the four filtered microphone outputs into a data stream 59 comprising four interleaved sequences of digitised samples e1(n), e2(n), e3(n) and e4(n), which correspond respectively to the filtered outputs of the microphones 421, 422, 423 and 424.
  • each sequence is represented by e (n) so that in this example l may be 1, 2, 3 or 4.
  • the processor and memory unit 36 receives a square wave signal 60 at 1.2 kilohertz from a sample rate oscillator 47 which determines the rate at which the converters 35, 41 and 46 convert samples and the frame duration of processing carried out by the unit 36.
  • the unit 36 completes its processing frame within 833 milliseconds.
  • a crystal clock oscillator 61 with a frequency of 10 Megahertz is included in the unit 36.
  • :- y1(n) w10x(n) + w11x(n - 1)
  • y2(n) w20x(n) + w21x(n - 1).
  • the filtered reference signal r lm is a sequence formed by filtering the reference signal x(n) with a filter that models the effect of the acoustic coupling between the m th loudspeaker and the l th microphone.
  • the unit 36 simulates this filtering as digital FIR (Finite Impulse Response) filtering. Coefficients for the digital FIR filtering are adjusted adaptively during an initialisation program in which a white noise generator 48 is energised.
  • a white noise signal is generated by the generator 48, which is then filtered by a low pass filter 49 to prevent aliasing, the filter 49 having a cut-off frequency of 460 Hz.
  • the signal is subsequently sampled and converted by an analog to digital converter 50.
  • the unit 36 performs a 35 coefficient FIR modelling of the impulse response between the m th loudspeaker and the l th microphone at the j th sample. Such modelling is described in "Adaptive Signal Processing" by B. Widrow and S.D. Stearns, published in 1985 by Prentice Hall.
  • the filtered reference sequence is then given by :-
  • the operation of the unit 36 is such that, having obtained the error samples e l (n) and the filtered reference signal r lm (n), each adaptive filter coefficient w mi for each output y m (n) is updated by a quantity proportional to the sum of the computed products of e l (n) and r lm (n - i) in accordance with the equation ;-
  • the new set of adaptive coefficients w mi is then stored and used to filter the next sample of the reference signal, x(n + 1).
  • the unit 36 includes RAM for temporary storage and computation, and EPROM for program storage. Calculated coefficients w mi and c lmj , and reference sequences r lm (n) are held in RAM. The convergence coefficient is entered at a set of manually operable switches (not shown).
  • the unit 36 includes a Texas Instruments TMS 32010 microprocessor.
  • the input signal rate from the ignition circuit, including the low tension coil 31, is 100 Hz to 200 Hz, and the waveform shaper 33 is a monostable circuit triggered by the leading edge of the input signal to produce pulses of a constant width which is small relative to the sample period set by the sample rate of 1.2 kilohertz.
  • the low pass filters 39, 44 and 49 are active filter modules supplied by Kemo Limited under No. 1431/L.
  • a family of algorithms which are alternatives to the SNM algorithm uses a bank of adaptive digital filters working in parallel for each secondary source.
  • Each individual filter is fed by a reference signal containing a subset of the harmonics or subharmonics to be controlled.
  • Figure 6 of the accompanying drawings shows two parallel FIR filters 70 each fed by pure tone reference signals 71, at the second and fourth engine order frequencies in this case. The outputs 72 of these filters are added together by an adder 73 to form an output 74 to the secondary source.
  • Each of the parallel filters 70 may be updated by any of the algorithms discussed above.
  • the said stochastic gradient algorithm may be modified so that :- where w Imi is the i'th coefficient of the FIR filter fed from the I'th harmonic of the engine, driving the m'th secondary source.
  • each harmonic frequency is controlled independently and the convergence of an harmonic does not couple with the convergence of any other harmonic, as is the case when a 2I coefficient filter is used to filter I harmonics simultaneously.
  • I filters each with 2 coefficients could be used to filter I harmonics individually, and their response combined afterwards.
  • the disadvantage of this algorithm is that a filtered reference signal needs to be generated for each source (m), sensor (l) and harmonic (I), to give each r lmI (n).
  • Another approach to controlling a number of harmonics is to take the Fourier transform of each of the error signals and to update a set of coefficients controlling each harmonic of each secondary output independently. The outputs of each of these filters are then combined together, for each secondary source, via an inverse Fourier transform, to generate the output waveform for this source, as indicated in Figure 7.
  • the complex value of the 'th error signal will be given by where A l is the value of E l with no active control, w m is the complex amplitude of the voltage to the m'th secondary source and c lm is the complex transfer function between the l'th sensor and m'th source at the frequency of the harmonic of interest.
  • the algorithm also reduces to the iterative matrix algorithm described in White and Cooper (1984, Applied Acoustics 17, 99-109. "An adaptive controller for multivariable active control”). See also U.K. Patent Specification No. 2,122,052A.
  • Figure 9 illustrates application of the invention to non-noise, i.e. mechanical, vibration control.
  • Figure 9 comprises a modification of the Figure 1 arrangement, wherein the microphone sensors have been replaced by accelerometers 90 and the loudspeakers by mechanical vibrators 91.
  • the accelerometers 90 and vibrators 91 are mounted on surface portions of the enclosure 10.
  • a sensing combination of microphones 12 and accelerometers 90 is used, and/or a source combination of loudspeakers 11 and vibrators 91.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)
  • Fittings On The Vehicle Exterior For Carrying Loads, And Devices For Holding Or Mounting Articles (AREA)

Abstract

To reduce noise inside a motor car passenger compartment, two loudspeakers (371, 372) are driven by signals derived from a reference signal x(n) by adaptive filtering carried out by a programmed microprocessor and memory unit (36) which adapts the filtering in dependence on error signals el(n) from four microphones (421, 422, 423 and 424) distributed in the passenger compartment. Reference filtering coefficients are initially determined by analysis of finite impulse responses when white noise is acoustically coupled from the loudspeakers (37) to the microphones (42), a white noise generator (48) being coupled to the unit (36). The reference signal x(n) is restricted to one or more selected harmonics or subharmonics of the fundamental noise frequency by a filter (34) which tracks the selected frequency. The selected frequency may be obtained from a coil (31) in the ignition circuit of the vehicle.

Description

    BACKGROUND TO THE INVENTION
  • The invention relates to active vibration control.
  • As used herein, the term "vibration" includes sound or noise, and the invention is particularly concerned with active noise control.
  • In passenger compartments of cars (automobiles), where a significant noise component is harmonically related to the rotation frequency of a reciprocating engine used to drive the car, the sound levels at low frequencies in such enclosures are difficult to attenuate using conventional passive methods and can give rise to subjectively annoying "boom". A method of actively attenuating a simple sound field by introducing a single secondary sound source driven so that its output is in antiphase with the original ambient noise is described in general terms by B. Chaplin in "The Chartered Mechanical Engineer" of January 1983, at pages 41 to 47. Other discussions are to be found in an article entitled "Active Attenuation of Noise - The State of the Art" by Glenn E. Warnaka at pages 100 to 110 in Noise Control Engineering, May-June 1982 and in Internoise 83 Proceedings, pages 457 to 458 and 461 to 464, and Internoise 84 Proceedings, pages 483 to 488. Particular methods and apparatus are also described in Bitish Patent Specification Nos. 1,577,322 and 2,149,614, as well as S.J. Elliot and I.M. Stothers: "A multichannel adaptive algorithm for the active control of start-up transients", Collogue Euromech 213, Sept. 1986, Marseille, France. The latter two documents, written by the inventors, relate closely to the invention. There is, however, a difference in the filter implementation which is used to take into account the reverberant response. This difference contributes to significantly improve the transient response of the system.
  • SUMMARY OF THE INVENTION
  • According to the present invention, an active vibration control system for reducing vibration generated by a primary source, is defined in Claim 1. It is characterised in that at least one reference signal containing selected harmonics of the said primary source vibration is supplied to means driving a plurality of secondary vibration sources, such that vibration energy detected by sensor means operable to sense the vibration field established by the primary and secondary sources is reduced.
  • As used herein, the term "harmonic" includes "subharmonics".
  • Preferably the control system is operable in accordance with an algorithm which adjusts the outputs from the secondary sources so as to substantially reduce a cost function on a time scale comparable with the delays associated with the propagation of vibration from the secondary sources to the sensor means.
  • The present invention is particularly concerned with an active noise reduction system which can control the sound throughout an enclosure of a car, or at one or a number of "quiet zones" within it, and which can quickly adapt to changes in the excitation of the sound field due to changes in, for example, engine load or speed.
  • In order to ensure that the sound produced by the secondary vibration sources is of the same frequency as that produced by the engine, a signal related to the engine crankshaft rotation rate, for example a signal emitted by the engine ignition system, is used to generate a reference signal containing a number of sinusoids at harmonics (or subharmonics) of the engine crankshaft rotation frequency. These are known as engine order frequencies. These sinusoids may be obtained using a variety of methods outlined below. Rather than attempt to control all harmonics, only a selected set of engine order frequencies may, in accordance with the invention, be generated as the reference signal. For example, only the firing frequency (second engine order in a four cylinder car) and its second harmonic (fourth engine order) are used if the spectrum of the sound in the car is dominated by these components.
  • Alternatively, a signal containing all the engine order frequencies may be fed to a band pass filter which isolates only the particular frequency or frequencies exciting a particular resonance in the enclosure interior which could cause a "boom" to be excited. The advantage of reducing the number of frequencies fed to the filter is that an adaptive filter having fewer coefficients than would otherwise be the case can be used. This makes implementation more efficient and allows faster adaption time. The faster adaption time is particularly important in automotive applications, for example, in which the active control system has to adapt sufficiently quickly to track changes in engine speed which may occur on a very short timescale.
  • Other reference signals may be obtained from transducers mounted on a road wheel hub or the suspension system of the car. Such reference signals would contain the harmonics of road wheel rotation or road noise. Transducers placed outside a vehicle may provide reference signals representative of wind noise. If the reference signals are periodic (deterministic), control at individual harmonics can be exercised as described hereinafter.
  • When the enclosed space is a motor vehicle interior, the secondary sources may be loudspeakers used as the low frequency drives of a car audio system.
  • Examples of methods in accordance with the invention for generating the reference signals from a signal from an internal combustion engine are now discussed.
  • 1. Selection of harmonics by filtering.
  • A signal is obtained from the primary vibration source which contains components at all harmonics prevalent in the sound in a vehicle powered by an engine. This signal is filtered so as only to leave the most important or dominant harmonics. Filtering is carried out by a filter whose centre frequency can be controlled by an external signal in such a manner that the critical filter frequencies have a constant ratio compared with the engine crankshaft rotation rate. This can be achieved by using, for example, charge coupled devices whose switching frequency is locked to the crankshaft rotation frequency, but can, also be implemented as a program running on a microprocessor, as described hereinafter.
  • 2. Selection of harmonics within a band by fixed filtering.
  • A primary source signal rich in harmonics is filtered by a band pass filter, having a centre frequency fixed at that of a pronounced "boom" in the car enclosure and a characteristic such that the reference signal only contains the harmonic(s) which are particularly exciting the boom. This may be extended such that the filter contains a number of resonances at a number of boom frequencies of the car, or even such that the filter has a frequency response which models the acoustic response of the car interior to the primary excitation. The input signal to the filter may be a signal from the engine containing all important harmonics and may be in the form of a pulse train.
  • 3. Generation of specific harmonics locked to the engine crankshaft rotation frequency.
  • This method may be accomplished by using phase lock loops to generate sinusoidal signals, with frequencies bearing an integer relationship to a square wave signal from the engine, which frequencies are then added together to form the reference signal. Alternatively, the signal derived from the engine can be used to control a number of tunable oscillators, each producing a sinusoid at a selected harmonic frequency. In one arrangement employing such a bank of tunable oscillators, the period of the square wave signal at the engine rotation rate is measured with a counter passed to a microprocessor which implements a number of digital oscillators using difference equations of the form

    x I (n) = δ(n) + 2 cos (Iω o )x I (n - 1) - x I (n - 2)
    Figure imgb0001


    where ω o = 2πf c /Nf s
    Figure imgb0002
    , I is the order of the harmonic or subharmonic to be generated, fc is the frequency o f the counter, which counts N pulses during a period, and fs is the sample rate used for the difference equation. δ(n) is a unit sample sequence to initiate the oscillators formed by the difference equations above. The second and fourth harmonic may be generated for example.
  • It is advantageous if the sample frequency (fs) is derived from the counter frequency (fc) by a frequency division circuit for example, so that the ratio fs/fc is exactly an integer number.
  • An alternative difference equation which may be used to implement the digital oscillator has the form of a series approximation to a trigonometric function, for example :-
    Figure imgb0003

    The variable y is the accumulated phase of the oscillator, given by
    Figure imgb0004

    ω i (n) = Iω o
    Figure imgb0005
    where ωo is calculated from the measured period of the reference signal, as above, for every new sample (n). The series approximation above can be used for y(n) in the range -π/2<y(n)<π/2
    Figure imgb0006
    . For values of y(n) outside this range the symmetry properties of the cosine waveform are utilised, until y(n)>π. When y(n)≧π, the natural overflow properties of the two's complement number representation is used to allow y(n) to "wrap-around" to y(n)≧-π and the series expansions and symmetry properties discussed above are again used. In this way the digital representation of y(n) is kept within the range -π<y(n)<π, and xI(n) is within the range ± 1.
  • The calculations of the coefficients used in the difference equations forming the digital oscillators, together with the difference equations themselves may be implemented on a dedicated processor, or may form part of the program which also implements a controller that generates the outputs used to drive the secondary sources from the reference signals described above.
  • The controller is designed to be adaptive so as to quickly track changes in engine speed and load. The outputs of the secondary sources are adaptively controlled so that some measurable cost function is minimised. This cost function would typically be the sum of the mean square outputs from a number of microphones in the enclosed space. The controller can be implemented as a digital adaptive FIR filter, using the basic update algorithm described by S. Elliott and P. Nelson in "Electronics Letters", at pp. 979-981, 1985. A number of additions must be made to this basic algorithm, however, to enable it to work quickly and efficiently in this particular example. The basic algorithm, referred to hereinafter as the stochastic gradient algorithm, is presented below, in order to highlight the necessary alterations. If the i'th coefficient of the adaptive filter driving the m'th secondary source at sample number n is wmi(n), then each of these coefficients should be adjusted at every sample according to the equation
    Figure imgb0007

    where α is a convergence coefficient, e(n) is the sampled output from the ℓ'th sensor and r
    Figure imgb0008
    m(n) is a sequence formed by filtering the reference signal discussed above ( x(n), say) with a digital filter which models the response of the ℓ'th sensor to excitation of the m'th secondary source. These digital filters, generating each rℓm(n), have only two coefficients in the implementation described in the "Electronics Letters" article, since control at only a single, fixed, frequency was being attempted. In the present example, however, the digital filters must model the relevant response over a range of frequencies, governed by the frequency range which the active system is attempting to control. It has been found that under certain circumstances a digital filter need only model the overall delay in the response to ensure the stability of an adaptive filter. It is more common, however, to have the digital filters incorporate a delay and then some reverberant response. This may be implemented using either digital FIR or IIR filters whose coefficients are adjusted adaptively during an initialisation phase, so as to accurately match the desired responses. It is also possible to continue this initial adaption process during the operation of the active control system by feeding training signals to each secondary source which are suitably uncorrelated with each other and with the primary excitation. This may be necessary to track changes in the acoustic response of the enclosure. Alternatively, if the change is due to some well-defined cause, such as a passenger sitting down or a car window being opened, this change may be detected with mechanical transducers and the information used to switch between a cariety of filters modelling the response of the enclosure under a variety of conditions.
  • Another important consideration concerning the use of adaptive algorithm in this example concerns the effect of unwanted, low level, harmonic or subharmonic components in the reference signal. Suppose that the method used to generate the reference signal, as described above, is designed to produce only I frequencies. Even though there may only be I desired harmonics in such a system, there will in practice also be a number of other harmonics or subharmonic frequencies at low level, because of the finite cut-off rate of the filters, for example. These components may also be being generated by the primary source and therefore present in the enclosure and hence at the outputs of sensors such as microphones and thus the adaptive algorithm will attempt to cancel them by enormously amplifying the low level, spurious harmonic reference signals. This can cause numerical overflow problems in the adaptive filter coefficients. This may be prevented in a number of ways:
    • (1) The use of only 2.I coefficients in the adaptive filter.
    • (2) The deliberate injection of random noise into the reference signal.
    • (3) The insertion of a "leak" into the algorithm so that, in the equation above, the past coefficient value, wmi(n), is multiplied by a factor close to but not equal to unity before being updated.
    • (4) The addition of an extra term in the update equation which minimises a cost function involving "effort" as well as "error", as described in ISVR Technical Report No. 136, 1985, published by the Institute of Sound and Vibration Research, University of Southampton.
  • A number of other adaptive algorithms may also be implemented to adjust the coefficients of the digital filters in the controller driving the secondary sources. These alternative algorithms are best described in matrix form.
  • Assuming the availability of a sampled reference signal, x(n), which is correlated with the output of the primary source, but is unaffected by the action of the secondary sources. The output to the m'th secondary sources, ym(n), may be obtained by passing this reference signal through a digital filter whose i'th coefficient is wmi(n) at the n'th sample, so that
    Figure imgb0009

       The sampled output from the ℓ'th error sensor, e(n), is equal to the sum of the contributions from the primary source, d(n), and each of the secondary sources. The response of the path between the m'th secondary source and ℓ'th error sensor is modelled as a J'th order FIR filter with coefficients cℓmj so that
    Figure imgb0010

    Therefore
    Figure imgb0011

       In order to obtain a matrix expression for the error surface we must now make the assumption that the filter coefficients in the controller are time invariant, i.e. the controller only adapts very slowly compared to the time scale of the response of the system to be controlled. Then w mi (n-j) = w mi
    Figure imgb0012
    and
    Figure imgb0013

    If we let
    Figure imgb0014

    where the sequences rℓm(n) for each ℓ and m are called the filtered reference signals, then
    Figure imgb0015

    or

    eℓ(n) = dℓ(n) + r ̲ T (n) w ̲
    Figure imgb0016


    where
    Figure imgb0017

    So if

    e ̲ T (n)=[e₁(n) e₂(n) ... e L (n)]
    Figure imgb0018

    d ̲ T (n) = [d₁(n) d₂(n) ... d L (n)]
    Figure imgb0019


    then

    e ̲ (n) = d ̲ (n) + R ̲ (n) w ̲
    Figure imgb0020


    where

    R T(n) = [r₁ (n) r₂ (n) ... r L (n) ]

    If the cost function is written as
    Figure imgb0021

    where E is the expectation operator, then

    J = E( d ̲ T (n) d ̲ (n) ) + 2 w ̲ T E( R ̲ T (n) d ̲ (n) ) + w ̲ T E( R ̲ T (n) R ̲ (n)) w ̲
    Figure imgb0022


       Using the standard theory of matrix quadratic forms, the minimum of J, Jmin, is obtained with

    w ̲ opt = - A ̲ ⁻¹ b ̲ = -E( R ̲ T (n) R ̲ (n))⁻¹E( R ̲ T (n) d ̲ (n))
    Figure imgb0023


    giving
    Figure imgb0024

       The true gradient may be written

    ∂J ∂W = 2E( R ̲ T (n) d ̲ (n) + R ̲ T (n) R ̲ (n) w ̲ (n))
    Figure imgb0025

    ∂J ∂W = 2E( R ̲ T (n) e ̲ (n))
    Figure imgb0026


    The time domain steepest descent algorithm may thus be written as

    w ̲ (n+1) = w ̲ (n) - 2µE( R ̲ T (n) e ̲ (n))
    Figure imgb0027


       In a practical implementation the true expectation could be approximated by an MA or AR averaging process. Alternatively, the instantaneous gradient could be used to update each filter coefficient every sample, as in the time domain "stochastic gradient" algorithm

    w ̲ (n+1) = w ̲ (n) - 2µ R ̲ T (n) e ̲ (n),
    Figure imgb0028


    which is the matrix representation of the algorithm described above. It is clear from this formulation that the stability and convergence properties of this algorithm, in the limit of slow adaption, are governed by the eigenvalue spread of the matrix E(R T(n) R(n)). This matrix depends only on the response of the system to be controlled, the positioning of the sources and sensors within that system and the spectral properties of the reference signal, x(n). It is possible that an unfortunate placing of these sources and sensors could cause this matrix to become ill conditioned, so that it has a large eigenvalue spread. This would create very slow "modes" in the convergence properties of such an algorithm.
  • This problem could be removed by using a Newton's method algorithm, the exact form of which may be written as

    w ̲ (n+1) = w ̲ (n) - 2µE([ R ̲ T (n) R ̲ (n)])⁻¹E( R ̲ T (n) e ̲ (n))
    Figure imgb0029


       Again, various types of averaging could be used to give a practical approximation to the expectation operator. It should be noted however that the time independent matrix E(R T(n) R(n)) depends only on the response of the system to be controlled and on the reference signal, and these are assumed to be known and stationary. This suggests a variety of stochastic Newton's method ("SNM") algorithms. The most obvious of these uses a modified or "normalised" set of reference signals Q(n), such that

    E([ R ̲ T (n) R ̲ (n)])⁻¹ R ̲ T (n) = Q ̲ T (n).
    Figure imgb0030


       The computation of each of these reference signals will take somewhat longer than for R(n) alone since none of the elements of Q(n) are necessarily time delayed versions of any other elements. The complete SNM algorithm, again using instantaneous versions of E(R T(n) e(n)), becomes

    w ̲ (n+1) = w ̲ (n) - 2µ Q ̲ T (n) e ̲ (n).
    Figure imgb0031
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • Embodiments of the invention will now be described, by way of example only, with reference to the accompanying drawings, wherein:
    • Figure 1 is a block schematic diagram of an active noise control system associated with an enclosed space;
    • Figure 2(a), (b) and (c) are graphical representations of the behaviour of an element of the system,
    • Figure 3 is a block diagram of one form of reference signal generator,
    • Figure 4 is a block diagram of another form of reference signal generator,
    • Figure 5 is a block diagram of a circuit which incorporates a microprocessor,
    • Figure 6 is a block diagram illustrating how two reference signals are combined,
    • Figure 7 is a flow chart illustrating processing in a particular embodiment of the invention,
    • Figure 8 is a schematic diagram illustrating a heterodyne and averaging method of obtaining inphase and quadrature components of an error sequence in another embodiment of the invention,
    • Figure 9 illustrates application of the invention to non-noise vibration control, and
    • Figure 10 illustrates a modification of the arrangement shown in Figure 9.
    DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
  • In Figure 1 an enclosure 10, which is the interior of the passenger or driver compartment of an internal combustion engine driven vehicle, in this example, a motor car 100, is represented schematically together with an active sound control system 1 according to the invention. In this example, the system 1 employs two secondary sound sources 11, comprising two low frequency loudspeakers of a stereo audio system fitted to the car, and three acoustic sensors, comprising microphones 12. The loudspeakers 11 are driven by a controller circuit 13 which comprises a pair of adaptive filters 14. Each adaptive filter 14 drives a respective one of the loudspeakers 11 with an output signal 3 which the filter 14 produces as a result of its action on a reference signal 4 supplied thereto by a reference signal generator 15. The reference signal 4 is generated by the generator 15 from an input signal 16 which is periodic at the crankshaft rotation rate of the internal combustion engine 2.
  • The signal generator 15 may comprise a tracking filter.
  • The purpose of the outputs from the loudspeakers 11 driven by the controller 13 is to reduce the sound vibration field established by the primary and secondary sources, experienced within the enclosure 10. Since the primary source (engine 2) of the noise to be reduced is periodic, the reference signal 4 generated by the generator 15 is, in accordance with the invention, arranged to contain one or more sinusoidal components at harmonics (or subharmonics) of the crankshaft rotation rate of the engine 2. The adaptive filters 14 are adjusted automatically by output signals 5 from the sensor-microphones 12, corresponding adjustment being made simultaneously to the outputs of the loudspeakers 11 so as to substantially minimise a cost function on a time scale comparable with the delays associated with the propagation of sound vibrations from the loudspeakers 11 to the microphones 12. The cost function may comprise the sum of the mean square outputs of the microphones 12.
  • A decision is made beforehand as to what harmonics are to be selected; a decision which may vary from car to car.
  • The control system 1 does not employ a stored solution. Instead, it makes use of a plurality of closed loops, each loop comprising a microphone 12, the controller 13, and a loudspeaker 11, whereby signals from the microphone 12 are used to adapt the filters 14 controlling the loudspeaker 11, which has an influence on the output of the microphone as a result of acoustic response within the enclosure 10.
  • Each loop is accounted for by part of an algorithm which adjusts the outputs of the loudspeakers 11 as aforesaid, the algorithm being of the form :-

    w ̲ (n+1) = w ̲ (n) -2µ R ̲ T (n) e ̲ (n).
    Figure imgb0032

  • It will be noted that the loudspeakers 11 and the microphones 12 are distributed in the enclosure 10 in spaced relationship. The distribution, which varies from car to car, is adjusted to get substantial sound reductions throughout the enclosure 10.
  • It will also be noted that the system 1 employs as many closed loops as the number of sensors (12) multiplied by the number of secondary sources (11).
  • The manner in which each secondary source 11 effects every one of the sensors 12 is reflected in the algorithms referred to herein.
  • Furthermore, the system 1 employs more sensors (12) than secondary sources (11), whereby a controlled reduction of primary source vibration is achieved. This contrasts with presently-known systems employing the same number of secondary sources as sensors, whereby near perfect cancellation can be achieved at the sensor locations but vibration levels away from these locations are increased.
  • Since the reference signal 4 contains harmonics of the input signal 16, which signal is periodic at the engine crankshaft rotation rate, the reference signal 4 contains engine order frequencies. The signal generator 15 is arranged to select engine order frequencies that ensure that the sound produced by the loudspeakers 11 is of the same frequency or frequencies as the sound produced in the enclosure 10 by the engine 2, even during changes in engine conditions such as load or speed. The number of engine order frequencies in the reference signal 4 is restricted so that the adaptive filters 14 have a relatively small number of coefficients and can therefore adapt quickly.
  • All the coefficients are constantly being updated by the system, on a sample by sample basis. Thus there is no waiting for a final response before making another adjustment. The sample time is only a small fraction of the fundamental frequency of the primary source 2.
  • The input signal 16 can be obtained from another moving part of the engine or part of the ignition circuitry; for example.
  • Figure 2(a) illustrates graphically the response of the reference signal generator 15 when in the form of a tracking filter. That is to say, a filter having a centre frequency so controlled that the filter output frequencies have a constant ratio to the dominant input frequencies, so that in Figure 2(a), the frequency fo is N x (engine crankshaft rotation rate), and the reference signal 4 (Figure 1) contains only the first N harmonics of the engine rotation rate, where N is an integer. If the filter input signal is a voltage pulse train as represented by Figure 2(b), where 8T is the periodic time of the engine rotation rate, the first eight harmonics of the engine rotation rate are present in the reference signal 4. The spectrum, by Fourier analysis, of the reference signal 4 is then illustrated by Figure 2(c), in which A is amplitude. With the response of Figure 2(a), only the first six harmonics would be usable.
  • The tracing filter comprising the signal generator, can be in the form of charge coupled devices having a switching frequency locked to the engine crankshaft rotation rate.
  • An alternative form of reference signal generator, which employs a plurality of tracking band pass filters is illustrated by Figure 3 in which the input signal 16, a square wave at, for example, 128 times the engine drive shaft rotation rate, is divided first by 32 and then by 2. Division by 32 is achieved by a divider 6 which produces a square wave signal 7 at four times the engine crankshaft rotation rate, which signal is supplied to a bandpass filter 17. The filter 17 has a centre frequency f₄ which is arranged to track the fundamental frequency of the square wave signal 7 supplied thereto. The further division by 2 is achieved by a divider 8 which produces a square wave signal 9 at twice the engine drive shaft rotation rate. The signal 9 is then supplied to a bandpass filter 18 having a centre frequency f₂ and which is arranged to track the fundamental frequency of the square wave signal 9 supplied thereto. The bandpass filters 17 and 18 produce respectively sinusoidal output signals 7a, 9a at f₂ and f₄ which are linearly summed at an adder 19 to produce the required reference signal 4.
  • Further dividers and tracking bandpass filters can, of course, be incorporated in the circuit of Figure 3 so that the reference signal 4 contains the desired set of engine order frequencies.
  • Another form of reference signal generator may comprise a fixed frequency filtering circuit which selects harmonics and/or subharmonics from an input signal rich in the harmonics of the engine crankshaft rotation rate or firing rate. The filtering circuit may comprise a bandpass filter having a centre frequency fixed at the frequency of a pronounced resonance excited in the enclosure 10 (Figure 1) by the engine 2 or other primary source of vibration. For example, a bandpass filter may be arranged to have a response that models the acoustic response of a motor car passenger compartment to the engine.
  • A further form of reference signal generator may comprise a plurality of phase lock loops used to generate sinusoidal signals having respective frequencies with integer relationships to a square wave input signal from the engine 2 or other primary vibration source. The sinusoidal signals may then be added together to form the required reference signal. Thus a reference signal 4 comprising specific harmonics and/or harmonics locked to the primary source fundamental, such as engine crankshaft rotation rate, is generated.
  • An alternative generator for such a reference signal is illustrated in Figure 4 in which a square wave 20 at the primary source fundamental is used to control a plurality of tunable digital oscillators 25, 26 each producing a sinusoidal signal at a chosen harmonic or subharmonic frequency to be added at the adder 19 which produce the reference signal 4 by simple addition of the sinusoidal signals.
  • In the signal generator 27 of Figure 4, the square wave signal 20, which is at an engine crankshaft rotation rate, is supplied to a bistable circuit 21 which divides the rate signal by two and thereby produces a pulse train signal 20a in which the duration of each pulse is equal to the prevailing periodic time of the square wave signal 20. This periodic time is then measured by a counter 22 which is enabled throughout the duration of each positive pulse from the bistable circuit 21 and counts clock pulses supplied by a clock pulse generator 23. The clock pulses are generated at a fixed, suitably high rate fc.
  • The contents of the counter 22 are read at the end of each positive pulse from the bistable circuit 21 by a trigonometric function generator 24. The generator 24 is triggered by the trailing pulse of each positive pulse from the circuit 21 and generates two digital outputs representing respectively cos(2 ₀) and cos(4 ₀) where ₀ is given by

    ω₀ = 2πfc/Nfs
    Figure imgb0033


    in which N is the number of clock pulses counted by the counter 22 in the duration of one positive pulse from the circuit 21, and fs is a sample rate used in the two digital oscillators 25 and 26 which receive respectively the digital outputs cos(2ω₀) and cos(4ω₀) from the function generator 24. The digital sinusoidal outputs from the two oscillators 25 and 26 are superposed by a digital adder 19¹ which supplies the reference signal 4 as a digital signal. The trigonometric function generator 24, oscillators 25 and 26 and adder 19¹ can be implemented by a microprocessor with a suitable program. The microprocessor can be used to produce the reference signal 4 in the form :-

    x I (n) = δ(n) + 2 cos(Iω₀)x I (n - 1) - x I (n - 2)
    Figure imgb0034


    where I is the order of the harmonic or subharmonic being generated, δ(n) is a unit sample sequence which initiates the simulation of the oscillator, and n is the sample number.
  • Figure 5 represents in block form an active sound control system 30 for reducing the level of engine generated noise in the passenger compartment of a motor car. The car is provided with an ignition circuit including a low tension coil 31 from which a voltage signal 32 at the firing rate of the engine is taken and supplied to a waveform shaper 33 which in response thereto produces a pulse train at the engine firing rate. It is assumed in the present example that the engine firing rate is twice the engine crankshaft rotation rate fo. Thus the shaper 33 provides a signal having a fundamental frequency which is a single harmonic, (2f₀), of the crankshaft rate. A reference signal generator is provided in the form of a proprietory tracking filter 34, manufactured by Bruel and Kjaer under type number 1623. The tracking filter 34 receives the output of the shaper 33 as an input signal and as a trigger signal and produces a sinusoidal output signal at the selected harmonic 2fo. This sinusoidal signal is sampled with an analog to digital converter 35 to produce a reference sequence x(n) of digitised samples which are supplied as data to a processor and memory unit 36.
  • Mounted within the motor car passenger compartment 10 (not shown in Figure 5) are two loudspeakers 37₁ and 37₂, which are in positions normally used for car stereo reproduction. The loudspeakers 37₁, 37₂ are driven by a multiplexer 38 through respective low pass filters 39 and output amplifiers 40. The filters 39 have a cut off frequency of 460 Hz and are provided to prevent aliasing. The multiplexer 38, which contains sample and hold circuits for each output, is controlled by the processor and memory unit 36, through a control line 55, and receives a single input signal 57 from a digital to analog converter 41. The purpose of the loudspeakers 37₁ and 37₂ is to generate, in the passenger compartment, audio waves that will cancel those set up directly by mechanical transmission from the engine to the compartment. The digital to analog converter 41 is supplied by the processor and memory unit 36 with output data 58 which consists of two interleaved sequences of digitised samples y₁(n) and y₂(n). The data 58 is converted by the converter 41 into interleaved sequences of analog samples and separated into respective sequences by the multiplexer 38 for application to the low pass filters 39. Thus, in effect, the loudspeaker 37₁ is driven by the sequence y₁(n) and the loudspeaker 37₂ is driven by the sequence y₂(n). In Figure 5 each sequence of data 58 is represented by the expression ym(n), so that in this example m may be 1 or 2.
  • In order to ensure that the acoustic outputs from the loudspeakers 37₁ and 37₂ have the correct phase and amplitude to effect cancellation of the engine noise, error signals are picked up from the passenger compartment and utilised by the processor and memory unit 36. Acoustic error signals, if present, are sensed by four microphones 42₁, 42₂, 42₃ and 42₄, which are placed respectively either side of a driver headrest and a passenger headrest, there being only two seats in the compartment in the present example. The electrical outputs from the microphones 42₁ etc. are respectively amplified by amplifiers 43 and passed through low pass filters 44 to a four-input multiplexer 45 which supplies a single analog output to an analog to digital converter 46. The filters 44 are provided to prevent aliasing and have a cut-off frequency of 460 Hz.
  • The multiplexer 45 is controlled by the processor unit 36 by way of control line 56.
  • The multiplexer 45 and the converter 46 convert the four filtered microphone outputs into a data stream 59 comprising four interleaved sequences of digitised samples e₁(n), e₂(n), e₃(n) and e₄(n), which correspond respectively to the filtered outputs of the microphones 42₁, 42₂, 42₃ and 42₄. In Figure 5, each sequence is represented by e (n) so that in this example ℓ may be 1, 2, 3 or 4.
  • The processor and memory unit 36 receives a square wave signal 60 at 1.2 kilohertz from a sample rate oscillator 47 which determines the rate at which the converters 35, 41 and 46 convert samples and the frame duration of processing carried out by the unit 36. Thus in the present example, the unit 36 completes its processing frame within 833 milliseconds. A crystal clock oscillator 61 with a frequency of 10 Megahertz is included in the unit 36.
  • The unit 36 simulates two adaptive filters, each having two coefficients, so that :-

    y m (n) = w m0 x(n) + w m1 x(n - 1)
    Figure imgb0035


    describes the relation between the output sequence ym(n) to a loudspeaker and the reference signal x(n), where the coefficients are wm0 and Wm1. Hence with the two loudspeakers 37₁ and 37₂ :-

    y₁(n) = w₁₀x(n) + w₁₁x(n - 1) and
    Figure imgb0036

    y₂(n) = w₂₀x(n) + w₂₁x(n - 1).
    Figure imgb0037

  • The values of the coefficients Wm0 and wm1 are calculated by the unit 36 from the relationship :-
    Figure imgb0038

    in which α is a fixed convergence coefficient, rℓm(n - i) is a value of a filtered reference signal rℓm, and i = 0 or 1.
  • The filtered reference signal rℓm is a sequence formed by filtering the reference signal x(n) with a filter that models the effect of the acoustic coupling between the mth loudspeaker and the ℓth microphone. The unit 36 simulates this filtering as digital FIR (Finite Impulse Response) filtering. Coefficients for the digital FIR filtering are adjusted adaptively during an initialisation program in which a white noise generator 48 is energised.
  • In the initialisation program, a white noise signal is generated by the generator 48, which is then filtered by a low pass filter 49 to prevent aliasing, the filter 49 having a cut-off frequency of 460 Hz. The signal is subsequently sampled and converted by an analog to digital converter 50. The digital output of the converter 50 is used to drive the loudspeakers 37₁ and 37₂, by way of the processor and memory unit 36, and the resulting digital input to the unit 36 from the microphones 42₁, 42₂, 42₃ and 42₄ is used to determine the values of reference filter coefficients cℓmj where j = 0,....,34. The unit 36 performs a 35 coefficient FIR modelling of the impulse response between the mth loudspeaker and the ℓth microphone at the jth sample. Such modelling is described in "Adaptive Signal Processing" by B. Widrow and S.D. Stearns, published in 1985 by Prentice Hall.
  • The filtered reference sequence is then given by :-
    Figure imgb0039

    The operation of the unit 36 is such that, having obtained the error samples e(n) and the filtered reference signal rℓm(n), each adaptive filter coefficient wmi for each output ym(n) is updated by a quantity proportional to the sum of the computed products of e(n) and rℓm(n - i) in accordance with the equation ;-
    Figure imgb0040

    The new set of adaptive coefficients wmi is then stored and used to filter the next sample of the reference signal, x(n + 1).
  • The unit 36 includes RAM for temporary storage and computation, and EPROM for program storage. Calculated coefficients wmi and cℓmj, and reference sequences rℓm(n) are held in RAM. The convergence coefficient is entered at a set of manually operable switches (not shown).
  • Preferably the unit 36 includes a Texas Instruments TMS 32010 microprocessor. The input signal rate from the ignition circuit, including the low tension coil 31, is 100 Hz to 200 Hz, and the waveform shaper 33 is a monostable circuit triggered by the leading edge of the input signal to produce pulses of a constant width which is small relative to the sample period set by the sample rate of 1.2 kilohertz. The low pass filters 39, 44 and 49 are active filter modules supplied by Kemo Limited under No. 1431/L.
  • Only a small amount of separate additional RAM is required with the above-mentioned TMS 32010 microprocessor, which has considerable internal RAM and operates as described in the TMS 32010 Users' Guide published in 1983 by Texas Instruments Inc. Data buses between the unit 36 and the converters 35, 41, 46 and 50 are 12 bit buses. Other buses and lines required for synchronisation and control are omitted for clarity. It will be noted that the values of the coefficients wmi and cℓmj can be initially set to zero.
  • A family of algorithms which are alternatives to the SNM algorithm uses a bank of adaptive digital filters working in parallel for each secondary source. Each individual filter is fed by a reference signal containing a subset of the harmonics or subharmonics to be controlled. For example, Figure 6 of the accompanying drawings shows two parallel FIR filters 70 each fed by pure tone reference signals 71, at the second and fourth engine order frequencies in this case. The outputs 72 of these filters are added together by an adder 73 to form an output 74 to the secondary source. Each of the parallel filters 70 may be updated by any of the algorithms discussed above. For example, the said stochastic gradient algorithm may be modified so that :-
    Figure imgb0041

    where wImi is the i'th coefficient of the FIR filter fed from the I'th harmonic of the engine, driving the m'th secondary source.
  • The advantage of such an algorithm is that each harmonic frequency is controlled independently and the convergence of an harmonic does not couple with the convergence of any other harmonic, as is the case when a 2I coefficient filter is used to filter I harmonics simultaneously. In this case I filters each with 2 coefficients could be used to filter I harmonics individually, and their response combined afterwards. The disadvantage of this algorithm is that a filtered reference signal needs to be generated for each source (m), sensor (ℓ) and harmonic (I), to give each rℓmI(n).
  • Another approach to controlling a number of harmonics is to take the Fourier transform of each of the error signals and to update a set of coefficients controlling each harmonic of each secondary output independently. The outputs of each of these filters are then combined together, for each secondary source, via an inverse Fourier transform, to generate the output waveform for this source, as indicated in Figure 7.
  • For a single harmonic in the frequency domain, the complex value of the 'th error signal will be given by
    Figure imgb0042

    where A is the value of E with no active control, wm is the complex amplitude of the voltage to the m'th secondary source and cℓm is the complex transfer function between the ℓ'th sensor and m'th source at the frequency of the harmonic of interest.
  • This may be expressed in matrix form as :-

    E ̲ = A ̲ + CW ̲
    Figure imgb0043


    where :-
    Figure imgb0044

    The cost function in this case may be written as J = E H E where the superscript H denotes the complex conjugate of the transposed vector or matrix. Therefore :-
    Figure imgb0045

    Note that
    Figure imgb0046

    so that the steepest descent algorithm may be written :-

    W ̲ k+1 = W ̲ k - 2µ C ̲ H E ̲ k
    Figure imgb0047


    where W k and E k are the complex values of the filter response and error output respectively at the k'th iteration.
  • This algorithm does not seem to appear in the active control literature. This is probably because the Newton's method algorithm, below, is no more difficult to implement after initialisation to form the matrix needed to premultiply E k. However, in some cases it may be that the matrix C changes with time and a separate "identification" algorithm is used in parallel with the adaptive control algorithm to track these changes. In such cases the steepest descent algorithm may be considerably more computationally efficient to implement than that below.
  • The frequency demain version of the Newton's method algorithm may be written :-

    W ̲ k+1 = W ̲ k - 2µ( C ̲ H C ̲ )⁻¹ C ̲ H E ̲ k
    Figure imgb0048


    A special case of this algorithm does appear in the literature, in which the number of error sensors is equal to the number of secondary sources (L = M), so that C is a square matrix, and the algorithm reduces to :-

    W ̲ k+1 = W ̲ k - 2 C ̲ ⁻¹ E ̲ k
    Figure imgb0049


       This algorithm has been presented by Pierce (1985, David W. Taylor Naval Ship Research and Development Center Report No. 85/047. An algorithm for active adaptive control of periodic interface). If the convergence coefficient, µ, is set equal to one half, the algorithm also reduces to the iterative matrix algorithm described in White and Cooper (1984, Applied Acoustics 17, 99-109. "An adaptive controller for multivariable active control"). See also U.K. Patent Specification No. 2,122,052A.
  • In practice the most advantageous algorithm will probably be derived from a judicious mixture of time and frequency domain concepts. For example, one implementation of a Fourier transformer operating on each error sequence e (n) to produce the inphase and quadrature frequency components of e(n) at Iωo is illustrated in Figure 8 wherein integrators 80 and multipliers 81 are used. The slowly varying outputs of this circuit represent the real and imaginary parts of the frequency domain signal E defined above, so these signals can be used with any of the frequency domain algorithms discussed above to update a bank of adaptive filters driven by this frequency component driving each secondary source, as in Figure 6.
  • Figure 9 illustrates application of the invention to non-noise, i.e. mechanical, vibration control.
  • The example illustrated by Figure 9 comprises a modification of the Figure 1 arrangement, wherein the microphone sensors have been replaced by accelerometers 90 and the loudspeakers by mechanical vibrators 91. The accelerometers 90 and vibrators 91 are mounted on surface portions of the enclosure 10.
  • In another modification, illustrated by Figure 10, a sensing combination of microphones 12 and accelerometers 90 is used, and/or a source combination of loudspeakers 11 and vibrators 91.
  • A combination of Figure 10 need not be coincident; the components thereof could instead be spaced from each other.

Claims (22)

  1. An active vibration control system for reducing vibration generated by a primary source of vibration possessing dominant harmonic frequencies which may quickly change, comprising a processing means (36); first low-pass filter means (39) provided at an output of said processing means (36) and having a fixed cut-off frequency; reference means (31,33,34,35) to supply at least one reference signal representing at least one selected harmonic of said primary source of vibration to said processing means (36); said processing means (36) being operative to generate at least one drive signal using said reference signal, and output the or each drive signal to a plurality of secondary vibration sources (37) through said first low-pass filter means (39); second low-pass filter means (44) provided at an input of said processing means (36), and having a fixed cut-off frequency; sensor means (42) provided at one or more locations, and operative to sense a vibration field established at at least one said location by said primary and secondary sources (37) and output at least one error signal through said second low-pass filter means (44) to said processing means (36); a sample rate oscillator (47) which provides a constant sample rate signal to said processor means such that said reference signal and the or each error signal is sampled at a constant rate; said processing means (36) comprising an adaptive response filter means (14,70) having first filter coefficients to model the delay and reverberant response of said sensor means (42) to at least one output of said secondary vibration sources (37) over a wide range of frequencies; said adaptive response filter means (14,70) being operative to adaptively determine second filter coefficients in response to the or each error signal and to adjust the or each drive signal using said first and second filter coefficients, to reduce the vibration field sensed by said sensor means (42).
  2. An active vibration control system as claimed in Claim 1, wherein said sensor means comprises a number of sensors, and said number of sensors is more than said number of secondary vibration sources.
  3. An active vibration control system as claimed in Claim 1 or Claim 2, wherein said reference means (31,33,34,35) is operative to supply a reference signal representing at least two harmonics of said primary source of vibration to said processing means (36).
  4. An active vibration control system as claimed in any preceding claim, wherein said adaptive response filter means comprises an array of filters each having 35 first coefficients (Cℓmj) which model the response of said sensor means (42) to at least one output of said secondary vibration sources (37).
  5. An active vibration control system as claimed in any of Claims 1 to 3, wherein said adaptive response filter means (14,70) comprises a number I of filters (70) each having two second coefficients (wmi), where I is the number of harmonics in said reference signal.
  6. An active vibration control system as claimed in any of Claims 1 or Claim 2, wherein said adaptive response filter means (14,70) comprises filters (70) which are provided with a plurality of reference signals (71) each representing a single harmonic, said filters (70) having their outputs (72) combined to form an output (74) to said secondary sources and being independently adjustable.
  7. An active vibration control system as claimed in any of Claims 1 to 3, wherein said adaptive response filter means (14,70) comprises an array of filters and is operative to take the Fourier transform of the or each error signal, update a set of complex second coefficients (W k) for said array of filters which control the or each harmonic of said drive signal, and combine outputs of said array of filters via an inverse Fourier transform to generate the or each drive signal.
  8. An active vibration control system as claimed in any preceding claim, wherein said reference means (31,33,34,35) includes a reference signal filter (34) to filter a periodic input signal having its fundamental frequency locked to a predominant frequency of said primary source of vibration.
  9. An adaptive vibration control system as claimed in Claim 8, wherein said reference signal filter (34) comprises a tracking filter.
  10. An active vibration control system claimed in any of Claims 1 to 7, wherein said reference means (31,33,34,35) includes at least one tunable oscillator (25,26), the frequency thereof being controlled by a signal (20) indicative of the fundamental frequency of the primary source of vibration.
  11. An active vibration control system as claimed in any preceding claim wherein said adaptive response filter means (14,70) comprises an array of filters each having a plurality of first coefficients (Cℓmj) which are adjusted adaptively during an initialisation phase of operation of said system.
  12. An active vibration control system as claimed in any preceding claim wherein said adaptive response filter means (14,70) comprises an array of filters each having a plurality of first coefficients (Cℓmj), which are adjusted adaptively during the operation of said system by feeding training signals over a wide frequency range to each secondary source, which signals are suitably uncorrelated with each other and with the primary source of vibration.
  13. An active vibration control system as claimed in any preceding claim wherein said adaptive response filter means is operative to adjust the or each drive signal in accordance with an algorithm so as to substantially minimise a cost function on a time scale comparable with the delays associated with the propagation of vibration from said secondary vibration sources (37) to said sensor means (42).
  14. An active vibration control system as claimed in Claim 13, wherein said adaptive response filter means is operative to adjust the or each drive signal in accordance with an algorithm of the form:

    w ̲ (n+1) = w ̲ (n) - 2µ R ̲ T (n) e ̲ (n).
    Figure imgb0050


    where
       w(n+1) represents a vector of values of the second filter coefficients for the (n+1)th sample;
       w(n) represents a vector of values of the second filter coefficients for the nth sample;
       µ represents a convergence factor;
       R T(n) represents a matrix of signals obtained by filtering said reference signal using said first filter coefficients; and
       e (n) represents a vector of values of the error signals for the nth sample.
  15. An active vibration control system as claimed in Claim 13, wherein said adaptive response filter means is operative to adjust the or each drive signal in accordance with an algorithm of the form:

    w ̲ (n+1) = w ̲ (n) - 2µ Q ̲ T ( n ̲ ) e(n).
    Figure imgb0051


    where
       w(n+1) represents a vector of values of the second filter coefficients for the (n+1)th sample;
       w(n) represents a vector of values for the filter coefficients for the nth sample;
       µ represents a convergence factor;
       Q T(n) represents a modified matrix of filtered reference signals; and
       e(n) represents the complex values for the error signals for the nth sample.
  16. An active vibration control system as claimed in Claim 13, wherein said adaptive response filter means is operative to adjust the or each drive signal at a single harmonic in accordance with an algorithm of the form:

    W ̲ k+1 = W ̲ k - 2µ C ̲ H E ̲ k
    Figure imgb0052


    where
       W k+1 represents a vector of complex values of the filter response at the (k+1)th iteration;
       W k represents a vector of complex values of the filter response at the kth iteration;
       µ represents a convergence factor;
       E k represents a vector of complex values of the Fourier transform of the error signals at the kth iteration;
       C represents the matrix of transfer functions; and
       H denotes the complex conjugate of the transposed vector or matrix.
  17. An active vibration control system as claimed in Claim 13, wherein said adaptive response filter means is operative to adjust the or each drive signal in accordance with an algorithm of the form:

    W ̲ k+1 = W ̲ k - 2µ( C ̲ H C ̲ )⁻¹ C ̲ H E ̲ k .
    Figure imgb0053


    where
       W k+1 represents a vector of complex values of the filter response at the (k+1)th iteration;
       W k represents a vector of complex values of the filter response at the kth iteration;
       µ represents a convergence factor;
       C represents the matrix of transfer functions;
       E k represents a vector of complex values of the Fourier transform of the error signals at the kth iteration; and
       H denotes the complex conjugate of the transposed vector or matrix.
  18. An active vibration control system as claimed in any preceding claim including said secondary vibration sources (37) which comprise loudspeakers, and wherein said sensor means (42) comprises at least one microphone.
  19. An active vibration control system as claimed in any of Claims 1 to 17 including said secondary vibration sources (37) which comprise vibrators, and wherein said sensor means (42) comprises accelerometers.
  20. An active vibration control system as claimed in Claim 18 and Claim 19, wherein said secondary vibration sources (37) comprise a mix of loudspeakers and vibrators and said sensor means (42) comprise a mix of microphones and accelerometers.
  21. An internal combustion engine driven vehicle including the active vibration control system of any of Claims 1 to 20.
  22. An internal combustion engine driven vehicle as claimed in Claim 21 as dependent on any of Claims 1 to 18 and 20, wherein said loudspeakers comprise loudspeakers of a stereo audio system fitted to said vehicle.
EP87906587A 1986-10-07 1987-10-06 Active vibration control Expired - Lifetime EP0285632B1 (en)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
GB08624053A GB2203016A (en) 1986-10-07 1986-10-07 Active sound control apparatus
GB8624053 1986-10-07
GB08630422A GB2201858A (en) 1986-10-07 1986-12-19 Active noise control
GB8630422 1986-12-19

Publications (2)

Publication Number Publication Date
EP0285632A1 EP0285632A1 (en) 1988-10-12
EP0285632B1 true EP0285632B1 (en) 1993-06-09

Family

ID=26291394

Family Applications (1)

Application Number Title Priority Date Filing Date
EP87906587A Expired - Lifetime EP0285632B1 (en) 1986-10-07 1987-10-06 Active vibration control

Country Status (4)

Country Link
EP (1) EP0285632B1 (en)
JP (1) JP2890196B2 (en)
DE (1) DE3786165T2 (en)
WO (1) WO1988002912A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8160264B2 (en) 2008-07-30 2012-04-17 Fujitsu Limited Transfer function estimating device, noise suppressing apparatus and transfer function estimating method

Families Citing this family (54)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8818087D0 (en) * 1988-07-29 1988-09-01 Univ Southampton Improvements in/relating to attenuation of electro-magnetic field components
GB8903201D0 (en) * 1989-02-13 1989-03-30 Lotus Group Plc Noise suppression in vehicles
US5371802A (en) * 1989-04-20 1994-12-06 Group Lotus Limited Sound synthesizer in a vehicle
GB8908980D0 (en) * 1989-04-20 1989-06-07 Lotus Group Plc Treatment of noise in vehicles
GB2230920B (en) * 1989-04-25 1993-12-22 Topexpress Ltd Active sound and/or vibration control
DE69030289T2 (en) * 1989-08-11 1997-09-04 Bridgestone Corp Active vibration control device
JP2685917B2 (en) * 1989-08-31 1997-12-08 株式会社東芝 Silencer
US5091953A (en) * 1990-02-13 1992-02-25 University Of Maryland At College Park Repetitive phenomena cancellation arrangement with multiple sensors and actuators
JP2573389B2 (en) * 1990-03-23 1997-01-22 晴夫 浜田 Electronic silencing method and device
US5233540A (en) * 1990-08-30 1993-08-03 The Boeing Company Method and apparatus for actively reducing repetitive vibrations
GB2252657A (en) * 1990-10-24 1992-08-12 Lotus Car Vibration and sound generation in vehicles
GB9023459D0 (en) * 1990-10-29 1990-12-12 Noise Cancellation Tech Active vibration control system
US5245552A (en) * 1990-10-31 1993-09-14 The Boeing Company Method and apparatus for actively reducing multiple-source repetitive vibrations
JP2533695B2 (en) * 1991-04-16 1996-09-11 株式会社日立製作所 Muffled sound reduction device
US5224168A (en) * 1991-05-08 1993-06-29 Sri International Method and apparatus for the active reduction of compression waves
JP3114074B2 (en) * 1991-06-21 2000-12-04 株式会社日立製作所 Medical diagnostic equipment
JP3471370B2 (en) * 1991-07-05 2003-12-02 本田技研工業株式会社 Active vibration control device
US5404409A (en) * 1991-07-31 1995-04-04 Fujitsu Ten Limited Adaptive filtering means for an automatic sound controlling apparatus
DE69330568T2 (en) * 1992-03-12 2001-11-22 Honda Giken Kogyo K.K., Tokio/Tokyo Vibration and noise control system for motor vehicles
CA2096926C (en) * 1992-05-26 1997-09-30 Masaaki Nagami Noise controller
FR2692709B1 (en) * 1992-06-18 1994-09-16 Renault Active noise control device in the passenger compartment of a motor vehicle.
GB9218465D0 (en) * 1992-08-29 1992-10-14 Adaptive Control Ltd Active sound control systems and sound reproduction systems
US5381485A (en) * 1992-08-29 1995-01-10 Adaptive Control Limited Active sound control systems and sound reproduction systems
JP3276214B2 (en) * 1992-08-31 2002-04-22 マツダ株式会社 Vehicle vibration reduction device
JP2924496B2 (en) * 1992-09-30 1999-07-26 松下電器産業株式会社 Noise control device
JP3506449B2 (en) * 1992-12-04 2004-03-15 富士通テン株式会社 Noise control device
DE4402086A1 (en) * 1993-01-25 1994-07-28 Mazda Motor Fastening device for a drive unit
DE69423531T2 (en) * 1993-02-02 2000-07-20 Honda Giken Kogyo K.K., Tokio/Tokyo Vibration / noise reduction device
US5530764A (en) * 1993-03-19 1996-06-25 Mazda Motor Corporation Vibration control system for an automotive vehicle
US5361303A (en) * 1993-04-01 1994-11-01 Noise Cancellation Technologies, Inc. Frequency domain adaptive control system
DE4432747C2 (en) * 1993-09-17 1997-03-27 Hitachi Medical Corp Device and method for noise reduction in an MRI scanner
JP3419865B2 (en) * 1993-12-28 2003-06-23 富士重工業株式会社 Noise reduction device
GB2287851A (en) * 1994-03-25 1995-09-27 Lotus Car Time domain adaptive control system for active noise cancellation
JPH0830278A (en) * 1994-07-14 1996-02-02 Honda Motor Co Ltd Active vibration control device
FR2723464B1 (en) * 1994-08-05 1996-10-11 Matra Cap Systems Sa METHOD AND DEVICE FOR ACTIVE DAMPING OF MECHANICAL WAVES
FR2724467B1 (en) * 1994-09-09 1996-11-22 Matra Cap Systems Sa METHOD AND DEVICE FOR ACTIVE DAMPING OF MECHANICAL WAVES WITH REMOTE SENSORS
GB2293898B (en) * 1994-10-03 1998-10-14 Lotus Car Adaptive control system for controlling repetitive phenomena
US5754662A (en) * 1994-11-30 1998-05-19 Lord Corporation Frequency-focused actuators for active vibrational energy control systems
US5692055A (en) * 1996-09-24 1997-11-25 Honda Giken Kogyo Kabushiki Kaisha Active noise-suppressive control method and apparatus
DE19746523B4 (en) * 1996-11-02 2006-06-08 Volkswagen Ag Method for generating sound in motor vehicles and a circuit arrangement for carrying out the method
US7062049B1 (en) 1999-03-09 2006-06-13 Honda Giken Kogyo Kabushiki Kaisha Active noise control system
WO2000068969A1 (en) 1999-05-10 2000-11-16 Nippon Electric Glass Co., Ltd. Glass bulb for cathode-ray tube
JP4581293B2 (en) * 2001-06-04 2010-11-17 株式会社デンソー Noise control device
JP3843082B2 (en) 2003-06-05 2006-11-08 本田技研工業株式会社 Active vibration noise control device
JP4077383B2 (en) 2003-09-10 2008-04-16 松下電器産業株式会社 Active vibration noise control device
GB0419512D0 (en) * 2004-09-02 2004-10-06 Avon Vibration Man Syst Ltd Controlling vibrations
JP4074612B2 (en) 2004-09-14 2008-04-09 本田技研工業株式会社 Active vibration noise control device
US8027484B2 (en) * 2005-07-27 2011-09-27 Panasonic Corporation Active vibration noise controller
DE102006045627A1 (en) * 2006-09-27 2008-04-03 Robert Bosch Gmbh Arrangement for active noise compensation for an electronically controlled windscreen wiper system of a motor vehicle
US8800302B2 (en) * 2012-07-16 2014-08-12 Sunpower, Inc. Driving an active vibration balancer to minimize vibrations at the fundamental and harmonic frequencies
US8860381B2 (en) 2012-07-16 2014-10-14 Sunpower, Inc. Balancing vibrations at harmonic frequencies by injecting harmonic balancing signals into the armature of a linear motor/alternator coupled to a Stirling machine
US11279265B2 (en) 2018-09-06 2022-03-22 Honda Motor Co., Ltd System and method for counteracting vibrations within a vehicle
DE102021120317A1 (en) 2021-08-04 2023-02-09 Sick Ag Sensor and method using a sensor
CN114368268B (en) * 2021-08-24 2023-11-03 华为数字能源技术有限公司 Power assembly, noise suppression method and electric automobile

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1983002031A1 (en) * 1981-11-26 1983-06-09 Chaplin, George, Brian, Barrie Improved method of an apparatus for cancelling vibrations from a source of repetitive vibrations
US4473906A (en) * 1980-12-05 1984-09-25 Lord Corporation Active acoustic attenuator

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2776020A (en) * 1955-02-09 1957-01-01 Gen Electric Noise reducing system for transformers
GB1577322A (en) * 1976-05-13 1980-10-22 Bearcroft R Active attenuation of recurring vibrations
US4232381A (en) * 1979-06-08 1980-11-04 Northrop Corporation Noise cancellation using tracking filter
GB2122052B (en) * 1982-06-09 1986-01-29 Plessey Co Plc Reducing noise or vibration
JPS599699A (en) * 1982-07-07 1984-01-19 日産自動車株式会社 Control of sound field in chamber of automobile
JPS59133595A (en) * 1982-11-26 1984-07-31 ロ−ド・コ−ポレ−シヨン Active sound attenuator
GB8328997D0 (en) * 1983-10-31 1983-11-30 Secr Defence Active noise reduction

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4473906A (en) * 1980-12-05 1984-09-25 Lord Corporation Active acoustic attenuator
WO1983002031A1 (en) * 1981-11-26 1983-06-09 Chaplin, George, Brian, Barrie Improved method of an apparatus for cancelling vibrations from a source of repetitive vibrations

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
B. Widrow und S.D. Stearns: "ADAPTIVE SIGNAL PROCESSING", Prentice-Hall, Inc., 1985 (pages 182-192, 142-143, 288-293) *
Black, "Modulation Theory", D. Van Nostrand Company Inc., 1953, Chapter 4 *
H.F. Olsen and E.G. May: "ELECTRONIC SOUND ABSORBER", J. Acoust. Soc. Am. Vol. 25, Nr. 6, November 1953 (pages 1130-1137) *
K. Kido and S. Onoda: "AUTOMATIC CONTROL OF ACOUSTIC NOISE EMITTED FROM POWER TRANSFORMER BY SYNTHESIZING DIRECTIVITY": Electrical Institute Repor 23 97HO, 1972, Tohoku University, Sendai, Japan (pages 97-110) *
Openheim and Schaefer, "Digital Signal Processing", Prentice-Hall, 1975, pages 26-30 *
S. J. Elliott and P. A. Nelson: "AN ADAPTIVE ALGORITHM FOR MULTICHANNEL ACTIVE CONTROL"; Proc. Inst. Acoustics Vol. 8(1), March 1986 (pages 135-147) *
S.J. Elliot and I.M. Stothers: "A MULTICHANNEL ADAPTIVE ALGORITHM FOR THE ACTIVE CONTROL OF START-UP TRANSIENTS"; Colloque Euromech 213, September 1986, Marseille (4 pages) *
S.J. Elliot and P.A. Nelson: "ALGORITHM FOR MULTICHANNEL LMS ADAPTIVE FILTERING", Electronics Letters, 10 October 1985, Vol. 21 (pages 979-981) *
Schwartz, Bennet u. Stein, "Communication Systems and Techniques", McGraw-Hill, 1966, pages 82-85 *
Taub und Schilling, "Principles of Communication Systems", McGraw-Hill, 1971, pages 157-195 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8160264B2 (en) 2008-07-30 2012-04-17 Fujitsu Limited Transfer function estimating device, noise suppressing apparatus and transfer function estimating method

Also Published As

Publication number Publication date
DE3786165D1 (en) 1993-07-15
JP2890196B2 (en) 1999-05-10
EP0285632A1 (en) 1988-10-12
JPH01501344A (en) 1989-05-11
WO1988002912A1 (en) 1988-04-21
DE3786165T2 (en) 1993-10-14

Similar Documents

Publication Publication Date Title
EP0285632B1 (en) Active vibration control
US5170433A (en) Active vibration control
US5222148A (en) Active noise control system for attenuating engine generated noise
EP0479367B1 (en) Method and apparatus for attenuating engine generated noise
US5359662A (en) Active noise control system
US5321759A (en) Active noise control system for attenuating engine generated noise
US5629986A (en) Method and apparatus for intelligent active and semi-active vibration control
GB2201858A (en) Active noise control
US5691893A (en) Adaptive control system
EP0724762B1 (en) Active control system for noise shaping
EP0712115A2 (en) Active noise and vibration control system accounting for time varying plant, using residual signal to create probe signal
Sommerfeldt et al. Adaptive control of a two‐stage vibration isolation mount
Kuo et al. Frequency-domain periodic active noise control and equalization
Tokhi et al. Design and implementation of self-tuning active noise control systems
JP2827374B2 (en) Active noise control device
JPH07133842A (en) Active vibration control device
Johansson et al. Convergence analysis of a twin-reference complex least-mean-squares algorithm
US5434925A (en) Active noise reduction
JP3579898B2 (en) Vehicle vibration control device and vibration control method
US5953428A (en) Feedback method of noise control having multiple inputs and outputs
Clark et al. A novel approach to feedforward higher‐harmonic control
EP0507829B1 (en) Active vibration reducing system
Kuo Adaptive active noise control systems: algorithms and digital signal processing (DSP) implementations
JPH06266373A (en) Noise cancellation system
Wright Innovative techniques in active noise control using active material-based actuators and signal processing

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: A1

Designated state(s): BE DE FR GB IT NL SE

17P Request for examination filed

Effective date: 19881013

17Q First examination report despatched

Effective date: 19910416

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): BE DE FR GB IT NL SE

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Effective date: 19930609

REF Corresponds to:

Ref document number: 3786165

Country of ref document: DE

Date of ref document: 19930715

ET Fr: translation filed
ITF It: translation for a ep patent filed
PLBI Opposition filed

Free format text: ORIGINAL CODE: 0009260

26 Opposition filed

Opponent name: NOISE CANCELLATION TECHNOLOGIES INC.

Effective date: 19940307

NLR1 Nl: opposition has been filed with the epo

Opponent name: NOISE CANCELLATION TECHNOLOGIES INC.

PLBL Opposition procedure terminated

Free format text: ORIGINAL CODE: EPIDOS OPPC

PLBM Termination of opposition procedure: date of legal effect published

Free format text: ORIGINAL CODE: 0009276

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: OPPOSITION PROCEDURE CLOSED

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: BE

Payment date: 19961004

Year of fee payment: 10

27C Opposition proceedings terminated

Effective date: 19960531

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: NL

Payment date: 19961031

Year of fee payment: 10

NLR2 Nl: decision of opposition
PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 19971031

BERE Be: lapsed

Owner name: ADAPTIVE CONTROL LTD

Effective date: 19971031

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: 19980501

NLV4 Nl: lapsed or anulled due to non-payment of the annual fee

Effective date: 19980501

REG Reference to a national code

Ref country code: GB

Ref legal event code: IF02

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20061002

Year of fee payment: 20

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: IT

Payment date: 20061031

Year of fee payment: 20

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DE

Payment date: 20061220

Year of fee payment: 20

REG Reference to a national code

Ref country code: GB

Ref legal event code: 732E

REG Reference to a national code

Ref country code: GB

Ref legal event code: PE20

REG Reference to a national code

Ref country code: FR

Ref legal event code: TP

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 EXPIRATION OF PROTECTION

Effective date: 20071005

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 20061016

Year of fee payment: 20