EP0470153B1

EP0470153B1 - Active sound and/or vibration control

Info

Publication number: EP0470153B1
Application number: EP90907216A
Authority: EP
Inventors: Graham Paul Eatwell; Christopher Mark Dorling; William Richard Hodson
Original assignee: Noise Cancellation Technologies Inc
Current assignee: Noise Cancellation Technologies Inc
Priority date: 1989-04-25
Filing date: 1990-04-20
Publication date: 1996-02-28
Anticipated expiration: 2010-04-20
Also published as: DE69025604D1; AU5545690A; JPH04505221A; ES2084028T3; GB2230920B; DE69025604T2; CA2049332C; CA2049332A1; GB2230920A; GB8909433D0; WO1990013108A1; EP0470153A1; AU635266B2

Abstract

An active sound or vibration control system for compensating noise or vibration arising from a periodic source of changing periodicity, wherein sensors (6) sense the existing noise and vibration, the sensor signals are sampled several times per cycle by an ADC (8) triggered by a signal output from a sensor (3) detecting the position of the source in its cycle, the sampled signals are transformed into compensation actuator drive signal values by a series of transform and adaptation modules (11 to 13, 15) based on algorithms dependent on the frequency of the source, and the drive signal values are stored in a memory means (1) continually updated by the values derived from the sampled sensor signals.

Description

The present invention relates to an active sound or vibration cancelling apparatus for cancelling sound or vibration from a source thereof.
In contrast to previous systems aimed at controlling periodic sound or vibration, the system of the invention can be used even when the fundamental period of vibration is changing rapidly. For example, it can be used to control the engine noise in the interior of a vehicle.
The improved apparatus in accordance with the invention uses orthogonal transformations to reduce a multichannel control system to a series of single channel systems and provides a method by which the output of each such system can be adapted to maintain good performance of the control system even when the fundamental frequency of the vibration or sound source is changing.
The principles of active sound and vibration control have been known for many years and there is a wealth of published literature. Most patent specifications in this field relate to methods applicable to particular situations. The apparatus described herein relates to the control of periodic or almost periodic sound and vibration. The two main approaches to this problem are:

(i) Wave shaping or filtering, eg US-A-4506380 and GB-A-2201858, where a reference signal containing one or more frequencies of the unwanted sound and vibration is filtered to produce the signals to send to actuators which in turn produce the desired sound or vibration.
(ii) Waveform synthesis, where a waveform generator is triggered by a signal derived from the source, eg GB-A-1577322.

The two methods are equivalent only if the vibration source is exactly periodic. If the source characteristics are changing in time it is usual to use an adaptive control system in which sensors in the region to be controlled sense the residual sound or vibration and pass the information to a processor which alters the filter coefficients or the synthesized waveform so as to provide better control. GB-A-2201858 describes methods for adapting filter coefficients. GB-A-1577322 recognises the need for adaption and a later patent specification, GB-B-2107960, describes a simple technique for such a system using a single actuator and sensor. This latter patent specification does not explain how to control vibration where the period is changing, except to suggest that in this case the transform technique should produce frequency components from the lowest expected frequency to the highest, rather than just at frequencies corresponding to the harmonics of the source.
A further Patent Specification GB-B-2122052, uses a waveform synthesis technique for vibration control. In this method a sensor and actuator are placed at each of a number of locations. This results in a system with equal numbers of sensors and actuators and a method for adapting the waveform is presented for this special case. In most applications, however, the sources and sensors are not colocated and usually more sensors than sources are used in an effort to obtain a better measure of the resulting sound or vibration.
EP-A-0043565 discloses an active sound or vibration cancelling apparatus for cancelling sound or vibration from a source thereof, comprising: a source sensing means for producing a source signal related to the phase of unwanted sound or vibration at the source sensing means; sound or vibration sensing means; sound or vibration producing means; analogue-to-digital converter means for sampling the output signals of the sound or vibration sensing means in dependence on the source signal; and processing means responsive to the output of the analogue-to-digital converter means to produce drive signals for the sound or vibration producing means to effect cancellation of noise or vibration from the source thereof.
According to the present invention, there is provided an apparatus having features as defined with reference to EP-A-0043565 wherein the sound or vibration sensing means comprises a distributed plurality of sound or vibration sensors, the sound or vibration producing means comprises a distributed plurality of sound or vibration producing actuators and the processing means uses data in the form of singular values representing a singular value decomposition of a matrix representing transfer functions between said plurality of actuators and said plurality of sensors to calculate said drive signals.
The theoretical background to the present invention will now be described. The numbered mathematical equations referred to are set out in accompanying drawings.
The signal from each of a plurality of sensors is sampled using an analogue to digital converter (ADC) triggered by a signal related to the position of the source in its cycle. The data may be averaged over several cycles to improve accuracy. This gives an almost periodic sequence to which an orthogonal transform, such as the discrete Fourier transform, can be applied. This process is well known for the analysis of periodic signals, and is referred to as "order ratio analysis" or "order locked analysis".
The sample signal from the i-th sensor is given by equation (3.1), where I_ij(nT) is the response at sensor i, due to an impulse at the j-th controller output, Y_i(n) is the n-th value of sensor signal in the absence of any control and T is the sampling interval. J is the number of controller outputs. A slightly more complicated expression must be used if the length of the impulse response is comparable with the time over which the sampling period changes significantly. If r_i is sampled N times per cycle, then since x_j is periodic, equation (3.2) is applicable, where NT is the fundamental period. Equation (3.1) can then be written as equation (3.3), where equation (3.4) defines the cyclic impulse response.
An orthogonal transform can be used to simplify equation (3.3).
An example of this is a discrete Fourier transform defined by equation (3.5), where f=1/NT is the fundamental frequency.
Equation (3.3) then becomes equation (3.6).
It is to be noted that, since R_i, Y_i and X_j are sampled an exact number of times per cycle, they do not depend on the frequency, f. Equation (3.6) shows that each harmonic, k, of the system can be considered separately.
The control problem is to find the components X_j(k) which produce the desired values of R_i(k). This problem is complicated because all of the control outputs, X_j(k) interact to produce each sensor signal. It is possible, however, to use a technique which transforms the set of coupled equations (3.6) into a set of independent equations. The technique employs a singular value decomposition of the transfer function matrix A_ij(kf) for each kf. This gives equation (3.7), where the asterisk denotes complex conjugation. The matrices with complex components U_im and V_mj represent orthogonal transformations and so have the properties given by equations (3.8) and (3.9), where M is the number of sensors and δ_ℓm is the Kronecker delta. The term D _m (kf) is the m-th singular value at frequency kf. It is a real quantity. The method of decomposition is described in "Numerical recipes - the art of scientific computing" by W H Press and others, Cambridge University Press, 1986, pages 52 to 64. Equation (3.6) can be multiplied by U^* _ℓi and summed over i to give equation (3.10), to which equations (3.10.1) and (3.10.2) and (3.10.3) are applicable.
These quantities are called the principal components of the corresponding signals. Equation (3.10) is a single equation for the component
(kf) of the desired controller output, which can be solved directly if
and
are known or, since
may be changing, can be solved iteratively using standard adaption algorithms. If the explicit dependence on _ℓ and kf is dropped, equation (3.10) reduces to equation (3.11).
If the aim is to make as
small as possible, one algorithm, at the n-th step, results in equation (3.12), where µ is a real convergence factor.
Using equation (3.11) and (3.12) gives equation (3.12.1), and from equation (3.11), equation (3.12.2) results.
These can be combined to give equation (3.13) and this shows that the algorithm is stable provided equation (3.14) is applicable, whereby optimal convergence is obtained when µD = 1.
Hence it is desirable that equation (3.15) applies, that is, a different convergence factor is used for each frequency and each principal component.
In order to implement this algorithm it is necessary to measure the transfer functions A_ik(kf) at a number of different frequencies, kf. This can be done during an initial start-up or calibration phase and if necessary can be adapted using a parameter estimator as described in UK Patent Application 8825074.1. The transformation matrices U(kf) and V(kf) and the singular values D_ℓ(kf) are calculated from the measured transfer functions and stored for each frequency. During operation the frequency f (or, equivalently, the period T) is measured so that the appropriate transformation matrices and singular values can be used. Since kf is unlikely to correspond exactly to a value for which the transfer function was measured, the nearest value is used. Alternatively interpolation between nearby values could be used to obtain more accuracy. In order to maintain a given accuracy the former method uses more memory and the latter uses more computation time.
Once
(kf) has been found, equations (3.9) and (3.10.3) can be used to give equation (3.16).
It is then possible to apply an inverse discrete Fourier transform to obtain x_j(n). These control signals are sent to digital to analogue converters (DACs), then filtered and amplified to provide the drive signals for the actuators.
In some applications it is desirable that the actuators are not driven too hard, and it is important that the signals to the DAC's are within the correct range. One particular method of limiting the drive amplitudes is to use a minimisation constraint, λ in the algorithm given by equation (3.17). The constraint λ can be adapted after each iteration, that is λ is increased if any of the drive signals x_j is too large or reduced if they are all in the desired range.
The invention is exemplified with reference to the accompanying drawings, in which the single figure following the invention shows one embodiment of apparatus according to the invention.
Digital values are stored in a memory device (1), which may for example be a FIFO device. These values are sent to a set of digital to analogue converters (DACs) (2) which are triggered N times per cycle by a train of electrical pulses from a sensor (3). These pulses relate to the position of the source in its cycle. The analogue signals from the DACs are passed through signal conditioners (4) to provide the drive signals for a number of actuators (5). The resulting sound or vibration field is measured by sensors (6). The signals from these sensors are used to adapt the values stored in the memory device (1) so that the sensor signals approach the desired values. The sensor signals are passed through signal conditioners (7) and then sampled in synchrony with the source using analogue to digital converters (8) which are triggered by signals from the position sensor (3). These sampled values are placed in memory device (9) and may be averaged over a number of complete cycles to reduce the effects of signals unrelated to the source. A transform module (10), which may use a discrete Fourier transform, produces components related to the harmonic frequencies of the source for each sensor. The components from the different sensors are then combined in the transform module (11) so as to produce the principal components of sensor signals. Each of these independent components is modified in the adaption module (12) to produce the principal components of the new drive signals. These are combined with transform module (13) to produce the frequency components of each drive signal which are then converted to time values via an inverse transform module (15). The new time values then replace those in the memory device (1). The transform modules (11) and (13) and the adaption modules (12) require knowledge of the period or frequency of the source. This may be obtained from the position signal via a frequency counter (14) which contains a real time clock. This method can be used in aircraft cabins where the source of the noise is the propellers or propfans.
An important application of the method of active control described above is in the control of engine related noise in vehicles. A control system for controlling the "boom" in automobile interiors is described in published UK Patent Application 2,201,858. It uses the wave shaping or filtering technique described above. The system is designed to adapt on a time scale comparable with delays associated with the propagation time of sound from the actuators to the sensors. In an automobile interior, however there is sound from many sources which are not related to the engine: for example, road noise, wind noise, sound from the in-car entertainment system. This noise contaminates the sensor signals and degrades the performance of the system.
The method of this invention uses averaging of the synchronously sampled signals over several cycles. This reduces the level of contamination and improves the performance of the system. However, the time taken for averaging reduces the ability of the system to track changes in the sound field due to changes in engine speed and load. Therefore, for a given level of contaminating noise, there will be an optimum number of cycles for averaging which will depend upon the rate of change of engine speed and load. The rate of change of engine speed may be obtained from the position signal and engine load may be obtained from additional sensors, such as a pressure sensor in the inlet manifold or throttle position sensor. This information can be used to control the rate of adaption so that optimal performance of the system can be obtained. This enables good performance to be obtained over a whole range of conditions rather than just at "boom" where the unwanted sound is much louder than the contaminating noise. Most modern automobile engines use computer controlled engine management systems. Some of the sensors could be used both for the active control system and the engine management system. Additionally, the same microprocessor could be used to control both systems.

Claims

An active sound or vibration cancelling apparatus for cancelling sound or vibration from a source thereof, comprising:
a source sensing means (3) for producing a source signal related to the phase of unwanted sound or vibration at the source sensing means;

sound or vibration sensing means (6);

sound or vibration producing means (5);

analogue-to-digital converter means (8) for sampling the output signals of the sound or vibration sensing means in dependence on the source signal; and

processing means (1,9,10,11,12,13,15) responsive to the output of the analogue-to-digital converter means to produce drive signals for the sound or vibration producing means (5) to effect cancellation of noise or vibration from the source thereof,
characterized in that the sound or vibration sensing means (6) comprises a distributed plurality of sound or vibration sensors (6), the sound or vibration producing means (5) comprises a distributed plurality of sound or vibration producing actuators (5) and the processing means (1,9,10,11,12,13,15) uses data in the form of singular values representing a singular value decomposition of a matrix representing transfer functions between said plurality of actuators (5) and said plurality of sensors (6) to calculate said drive signals.
An apparatus according to claim 1, wherein the processing means (1,9,10,11,12,13,15) produces said drive signals with reference to a look up table containing data representing the relevant transfer functions between said actuators (5) and said sensors (6) for a plurality of discrete source frequencies.
An apparatus according to claim 2, wherein the processing means (1,9,10,11,12,13,15) interpolates between the lookup table data when the table does not contain data for the current source fundamental frequency.
An apparatus according to claim 2 or 3, wherein said data comprises transfer function data at harmonics of said discrete source fundamental frequencies.