WO1990013108A1

WO1990013108A1 - Active sound and/or vibration control

Info

Publication number: WO1990013108A1
Application number: PCT/GB1990/000617
Authority: WO
Inventors: Graham Paul Eatwell; Christopher Mark Dorling; William Richard Hodson
Original assignee: Active Noise And Vibration Technologies Inc.
Priority date: 1989-04-25
Filing date: 1990-04-20
Publication date: 1990-11-01
Also published as: DE69025604D1; AU5545690A; JPH04505221A; ES2084028T3; GB2230920B; DE69025604T2; CA2049332C; CA2049332A1; GB2230920A; GB8909433D0; EP0470153A1; EP0470153B1; 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

Title: Active Sound and/or Vibration Control

Field of invention

This invention relates generally to systems for controlling sound or vibration, and more especially to active control systems which use a plurality of actuators to produce a—controlling -sound or vibration field-and—a— plurality of sensors to measure the residual field.

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 method 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.

Background to the invention

The principles of active sound and vibration control have been known for many years and there is a wealth of published literature. Most patent sp-ecifications in this field relate to methods applicable to particular situations. The method and system described herein relate 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 Patent No. 4,506,380 and published UK Patent Application No. 2,201,858, 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) ci-veforn.^~syrrt_h.esis, -where a waveform generator is triggered by a signal derived from the source, eg UK Patent Specification No. 1,577,322.

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. Published UK Patent Application No. 2,201,858 describes methods for adapting filter coef icients. UK Patent Specification No. 1,577,322 recognises the need for adaption and a later patent specification, UK Patent No. 2,107,960, 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 No. 2,122,052, 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 co- located and usually more sensors than sources are used in an effort to obtain a better measure of the resulting sound or vibration.

The invention

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. -(nT) is the response at sensor i, due to an impulse at the j-th controller output, x.(m) is the m-th value of the j-th controller output, y. (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. is sampled N times per cycle, then since x. 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 or-thog nal tra-nsform «an- be used to simpli y equation (3.3).

An example of this is a discrete Fourier transform defined by equation (3.5), where f=l/NT is the fundamental frequency.

Equation (3.3) then becomes equation (3.6).

It is to be noted that, since R 1. , Y1. and Xj. are samp^led 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.(k) which produce the desired values of R.(k). This problem is complicated because all of the control outputs, X.(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..(kf) for each kf. This gives equation (3.7), where the asterisk denotes complex conjugation. The matrices with complex components U. and V . represent orthogonal transformations and so have the properties given by equations (3.8) and (3.9), where M is the number of sensors and δ„ is the ronecker delta. The term D ⁽kf) St, m m 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

H Press and others, Cambridge University Press, 1986,

* pages 52 to 64. Equation (3.6) can be multiplied by U . 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 X (kf) of the desired controller output, which can be solved directly if Y and

____* ""^** «

R are known or, since Y 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 R 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., (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 X (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.(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. is too large or reduced if they are all in the desired range.

Description of embodiment

The invention-is--exempli-fied— ith -reference to the - - — - accompanying drawings, in which the single figure following the invention shows one embodiment of apparatus for implementing the method.

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 ffequency 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

1. An active sound or vibration control system for controlling periodic sound or vibration produced by a source having a changing fundamental period of vibration, comprising a distributed plurality of noise or vibration sensors, a distributed plurality of noise or vibration compensation actuators, at least one analogue to digital converter (ADC) for sampling the sensor signals, and a source sensor for producing a signal related to the position of the source in its cycle and which is used to trigger the at least one ADC a plurality of times per cycle, the outputs of the at least one ADC being utilised to produce drive signals for the actuators.

2. A system according to claim 1, wherein a memory means is used to store digital values, and the values stored in the .memory means are continually adapted or modified by sampled signals derived through the at least one ADC from the outputs of the noise or vibration sensors.

3. A system according to claim 2, wherein the memory means is a FIFO device.

A. A system according to claim 2 or claim 3, wherein the sampled values of the noise and vibration sensors are placed in a memory device and are averaged over a plurality of cycles.

5. A system according to claim 2 or claim 3 or claim , wherein signals output from the memory device are utilised in a signal analysing and processing means to produce the values fed to the memory means.

6. A system according to claim 5, wherein the signal analysing and processing means comprises a first transform module producing components related to the harmonic frequencies of the source for each noise and vibration sensor, a second transform module producing the principal independent components of the sensor signals, an adaptation processor for producing from said sensor signal principal components the principal components of the drive signals, and a third transform module for producing the frequency components of each drive signal from said drive signal principal components. ^~~

7. A system according to claim 6, wherein the signal analysing and processing means also includes an inverse transform module for producing time values from the frequency components, and said time values are fed to the memory means.

8. A system according to claim 6 or claim 7, wherein the period or frequency of the source is derived from the output of the source sensor by a frequency counter incorporating a clock and the source period or frequency signal is fed to the second and third transform modules and to the adaptation processor to control generation of the frequency components.

9. A system according to any of claims 2 to 8, wherein output signals from the memory means are fed to at least one digital to analogue converter (DAC) for producing the actuator drive signals.

10. A system according to claim 5, wherein the signal analysing and producing means operates in accordance with equations ⁽3.1⁾ to (3.16) or (3.17) hereinbefore set forth.