CA2001199C - Signal processing means for identifying systems subject to periodic disturbances - Google Patents

Signal processing means for identifying systems subject to periodic disturbances

Info

Publication number
CA2001199C
CA2001199C CA 2001199 CA2001199A CA2001199C CA 2001199 C CA2001199 C CA 2001199C CA 2001199 CA2001199 CA 2001199 CA 2001199 A CA2001199 A CA 2001199A CA 2001199 C CA2001199 C CA 2001199C
Authority
CA
Canada
Prior art keywords
plant
equation
disturbances
given
accordance
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 - Fee Related
Application number
CA 2001199
Other languages
French (fr)
Other versions
CA2001199A1 (en
Inventor
Graham Paul Eatwell
Colin Fraser Ross
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.)
Noise Cancellation Technologies Inc
Original Assignee
Noise Cancellation Technologies Inc
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
Application filed by Noise Cancellation Technologies Inc filed Critical Noise Cancellation Technologies Inc
Publication of CA2001199A1 publication Critical patent/CA2001199A1/en
Application granted granted Critical
Publication of CA2001199C publication Critical patent/CA2001199C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • G05B13/045Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance using a perturbation signal
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B21/00Systems involving sampling of the variable controlled
    • G05B21/02Systems involving sampling of the variable controlled electric

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Feedback Control In General (AREA)

Abstract

In a signal processing system for a plant such as a motor or machine, subject to disturbances which are approximately periodic and wherein a signal (T/t) time related to the disturbances is available, the plant response is sampled in an ADC triggered by the signal (T(t)) in order to synchronise the sampling process with the cyclic disturbances, and a processing means derives from the sampled signals an estimate of the disturbances and thereby of that part of the plant output due to unknown inputs (Figure 2), the estimate of disturbances possibly being used to drive an active vibration control system employed in combination with the sampling system.

Description

2~Q~gg Title: Signal Processing Means for Identifying Systems Subject to Periodic Disturbances.

Field of invention This invention relates to a signal processing system for characterising a plant subject to periodic or almost periodic forcing.

The expression plant is employed herein to describe any closed system such as an industrial plant, a motor vehicle, or aeroplane, or machine etc.

Background to the invention In many situations a plant is required to produce some desired output, d(t). The plant is controlled by inputs u(t) but the measured plant output ~(t) is subject to disturbance y(t) and noise n(t), (see Figure 1). To achieve good control it is necessary to estimate a model or inverse-model of the plant as described in B. Widrow'&
S.D. Stearns, "Adaptive Signal Processing", Prentice Hall Inc., 1985 for example. If the plant characteristics are time varying the plant model must be updated. There exist several methods for doing this as described in G.
C. Goodwin'& R.L. Payne, "Dynamic System Identification:
Experiment Design and Data Analysis", Academic Press, 1977 and A.P. Sage'& J.L. Melsa, "System Identification", Academic Press, 1971 for example, but they assume that the disturbances are stochastic ie random, and are not suitable for application where the disturbance is substantially periodic.

Summary of the invention In a plant in which the disturbances (y(t)) are substantially periodic and a signal (T(t)), which is time-related to the disturbance (y(t)), is available, a signal processing system is provided to estimate the disturbances and to characterize the plant, wherein a sampling of the plant response is performed in an analogue to digital converter (ADC) and the signal T(t) is used as a trigger to synchronize the sampling process of the Analogue to Digital Converter.

In particular the signal processing system estimates that part of the measured plant output which is due to unknown inputs and estimates the relationship between the plant output and the known plant inputs.

Brief description of the drawinqs Figure 1 illustrates a block diagram of a typical environment for practicing the present invention.
Figure 2 illustrates an embodiment of a signal processing system according to the present invention.
Figure 3 illustrates a block diagram of the signal processing system of the present invention incorporated into an active control system.

The sampled output Yn due to an input (u(t)) is approximated by equation (3.1).

In (3.1) (al, a2, a3, .... ) is the sampled plant response.

The first simplification occurs when the input to the plant u(t) is periodic.

- 2a -If u is sampled M times per cycle then un~M = un so that Yn can be written as in equation (3.3).

Hence, instead of needing to estimate the whole impulse response, an, we only need the M terms of a in equation (3.4). Also because the disturbance y(t) is periodic, it has only M independent terms after sampling. We introduce the vectors given by equations (3.5) and (3.6) .- .

2~C~

and the matrix U with elements given by equation (3.7). U
is a cyclic matrix, and the plant output E can be written as given in equation (3.8).

The noise, n(t), can be reduced by averaging over many cycles so that it can be ignored in the following description.

If the plant, the disturbance and the inputs are slowly varying compared with the cyclic variations, the output value of equation (3.8) at time t can be well approximated by equation (3.9).

Alternatively, if the plant is constant but the period of the disturbance and the input is changing, the impulse response a(t) will be sampled at sliqhtly different times so that a will change as the period changes. Thus the plant can again be approximated by (3.9). Therefore a constant plant subject to inputs of slowly changing periodicity can be modelled as a time-varying plant subject to periodic inputs, provided that the sampling of the ADC is synchronised by the external trigger signal T(t)- This has considerable advantages over methods which take no account of periodicity. For example, i) Averaging can be used to reject the noise n(t) ii) The vector a describes the plant response and contains far fewer parameters than would be required to characterise the whole impulse response.

iii) The cyclic nature and low order of the equation (3.9) reduces the computationa~ effort in the estimation of a and y(t), particularly since transforms, such as the 2Q(~

finite Fourier transform for example, can be used to diagonalise the matrix U.

In order to estimate a and y we require that small changes are made to U on a time-scale shorter than the plant or periodicity variations, so that (3.9) becomes the equation given at (3.10) where X(t) is given by expression (3.11) and b by (3.12) and I is the MxM identity matrix.

The small changes to U can either be as the result of changing control signals or, to avoid ill-conditioning, the result of psuedo-random changes made by the processor.

An orthogonal transform of (3.10) gives the expression shown in (3.13).

In (3.13), n is the transform variable, which corresponds to the harmonic number if a finite Fourier transform is used, and xnT and bn are as given in (3.14) and (3.15).

In (3.14) and (3.15) the superposed bar denotes the transform, and the superposed T denotes the transpose.

The extension to a multichannel system with S inputs is given by expression (3.16) for each plant output i. The vectors on the right hand side are defined in (3.17) and (3.18).

~xpression (3.16) and the single channel version, (3.13), are in a form which allows the use of standard algorithms for parameter estimation. One solution is as given by ( 3 . 1 9 ? as amplified by ~3.20) and (3.21).

2S~0~1 ~9 The solutions to (3.20) and (3.21) can be obtained recursively as in G.C. Goodwin -~ R.L. Payne, "Dynamic System Identification: Experiment Design and Data Analysis", Academic Press, 1977 and A.P. Sage & J.L.
Melsa, "System Identification", Academic Press, 1971.
\
One embodiment of the invention is shown in Figure 2.
The system shown differs from previous systems in that the sampling in the ADCs is controlled by the external trigger signal T(t) which is time related to the disturbance _(t).

A preferred configuration employs transform and inverse-transform modules. The transform has the property that it diagonalises circular matrices such as that occuring in equation (3.8). The solution can be obtained without the use of orthogonal transforms but this would generally require more computation.

A signal processing system embodying the invention can be in equation (3.8). The solution can be obtained without the use of orthogonal transforms but this would generally require more computation.

A signal processing system embodying the invention can be used in conjuction with an active control system. The desired signal is then zero, the disturbance y(t) is the output from sensors due to the sound or vibration to be cancelled, the signals u(t) are the drive signals sent to the actuators and the plant represents the acoustic or vibration path from the actuators to the sensors characteristics. The required actuator drives are as ~iven by equation (3.22).

2n~ ss In (3.22) the matrix An has components given by (3.23) and the vector Yn has components given by (3.24).

An is a generalised inverse of An.

If a complex finite Fourier transform is used then A+ is as given by equation (3.25).

In (3.25) the star denotes conjugate transpose and it has been assumed that the number of sensors is not less than the number of actuators. The drives can also be obtained iteratively for example using the LMS algorithm. At the k-th step the value of un is given by equation ! 3.~6) where JU is a convergence factor.

The signal processing system can be rearranged to give An directly. Equation (3.16) can be rewritten using (3.173 and (3.18) as (3.27), which leads to (3.28).

This can be solved in a similar fashion to give estimates of An and -Anyn and this gives an inverse model of the plant.

The incorporation of the signal processing system in an active control system is shown in Figure 3.

The input to the control system is also sampled by an ~DC.
This is triggered by the signal T(t), as is the Digital to Analogue Converter (DAC) which converts the digital output from the control system into the analogue drive signals to the actuators. The ADC can be common to both the control system and the signal processor which produces the estimates of the plant and the disturbance. The transform and inverse transform modules could also be common to both z~ 9 systems.

Any sensor failure will result in one row of the matrix A
becoming close to zero, while an actuator failure will result in a column of A becoming small. The system can monitor changes in the matrix A and identify and indicate failures in sensors or actuators. This is an important aid to system maintenance and repair.

If more than one periodic source is present an extra trigger signal for each additional source can be incorporated to allow the system to discriminate against these sources.

Z~l.99 FORMULA SHE~r 1/6 Y ~m an mUm ~m--~ Un-m m (3.1) Yn k- - ~ m-l Un- mam+kM ( 3 . 3 ) - . u a m-~ n-m m am k--~~ am+kM (3.4) ~ ~Yl ~ Y2 ~ - - - Ym) (3 ~ 5) Zt~ 9 FORMU LA SHEET ~/6 a - (al, a2, ...... am) (3.6) Umn Un-m (3. 7) - ua + y + n(t) (3 . 8 ~(t) U(t)a(t) + y(t) (3.9 (t) ~ U(t) a + y (3.10) - X(t)b 2~ 9 X(t) [U(t), I], (3.11) b ~1 (3.12) en(t) ~ xn(t) bn ,n -- 1, ........... ,M (3.13) --n [U(n, t) 1] (3 .14) b [ a (n) ] ( 3 .15 ) 2(~ 9 ei(t) - xT(t) bi (3.16) xT _ [Ul(n,t) ,U2(n,t)- ~ - - - - - - - - ~Us( ~ ) ~ ] (3.17) bn [ail(n), ai2(n), ......... ,ai (n) ~ yi(n)]T (3.18) bn ~xx ~ex (3.19) ; -- k~l xn(tk) xn(tk) (3 . 20) 20011 ~9 rORMULA SHEET 5/6 ~exN k-l Xn(tk) en(tk) (3 . 21) ~n ~ -AnX (3 . 22) [A ] ~ ( - ( 3 . 2 3 ) [~n]i Yi(n) . (3.24) 2Q(~1 99 + * 1 *
A -- (A A ) A (3 . 25) n n n n k k-l *
u u -~ A ~ (3 . 26) -n ~ n n~n en(t) - An--n(t) + Yn (3 . 27) n en(t) ~ An Xn(t) -- u (t) [ A - A Xn ] [ 1 ] ~ Un ( t ) ( 3 . 2 8 )

Claims (10)

  1. Claims l. A signal processing system for a plant subject to disturbances (y(G)) which are substantially periodic and wherein a signal (T(t)) time related to the disturbances is available, characterised in that, in order to estimate the disturbances and thereby characterise the plant, the plant response is sampled in an analogue to digital converter (ADC) and the signal T(t) is employed as a trigger and to synchronise the sampling process of the ADC.
  2. 2. A system according to claim 1, wherein, by means of the sampling process, the system estimates that part of a measured plant output due to unknown inputs and thereafter estimates the relationship between the total plant output and known plant inputs.
  3. 3. A system according to claim 1 or claim 2, wherein the sampled plant output (Yn) due to an input (u(t)) is approximated by equation (3.1) hereinbefore given and defined, is simplified for a periodic input (u(t)) which is periodic, whereby (Yn) is estimated in accordance with equation (3.3) hereinbefore given, which reduces to equation (3.4) hereinbefore given and defined.
  4. 4. A system according to claim 3, wherein vectors given by equations (3.5) and (3.6) together with a cyclic matrix (U) with elements given by equation (3.7) are introduced into equation (3.4) so that the plant output is measured in accordance with equation (3.8), all said equations being hereinbefore given.
  5. 5. A system accoding to claim 4, wherein the plant output (a(t)) at a time (t) is estimated in accordance with equation (3.9) hereinbefore given, both in a plant where the disturbances and/or inputs vary slowly compared with the cyclic variations, the sampling periodicity then being constant, and in a constant plant where the disturbances and/or inputs are changing, the sampling periodicity then being slowly varied.
  6. 6. A system according to claim 5, wherein noise (n(t)) is reduced to a negligible factor in the estimated output by averaging over many cycles.
  7. 7. A system according to claim 5 or claim 6, wherein transforms, such as finite Fourier transforms, are employed to diagonalise the matrix U.
  8. 8. A system according to any of claims 1 to 7, employed in conjunction with an active vibration control system having sensors for measuring unwanted sounds and vibrations to be cancelled, and wherein the signals (U(t)) of equation (3.9) are used as drive signals defined in accordance with equation (3.22) hereinbefore given and defined for actuators for actively cancelling said sounds and vibrations.
  9. 9. A system according to claim 8, wherein an inverse transform of the matrix (An) of equation (3.22) is developed to provide an input to the control system, in accordance with equation (3.25) hereinbefore given.
  10. 10. A system in accordance with claim 8, wherein the input to the control system is also sampled by an ADC
    triggered by the signal T(t) and a digital to analogue converter (DAC) is provided to convert the digital output of the control system into analogue drive signals for the actuators.
CA 2001199 1988-10-26 1989-10-23 Signal processing means for identifying systems subject to periodic disturbances Expired - Fee Related CA2001199C (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
GB888825074A GB8825074D0 (en) 1988-10-26 1988-10-26 Signal processing means for identifying systems subject to periodic disturbances
GB8825074.1 1988-10-26

Publications (2)

Publication Number Publication Date
CA2001199A1 CA2001199A1 (en) 1990-04-26
CA2001199C true CA2001199C (en) 1998-08-11

Family

ID=10645832

Family Applications (1)

Application Number Title Priority Date Filing Date
CA 2001199 Expired - Fee Related CA2001199C (en) 1988-10-26 1989-10-23 Signal processing means for identifying systems subject to periodic disturbances

Country Status (6)

Country Link
EP (1) EP0440720A1 (en)
JP (1) JPH04501476A (en)
AU (1) AU4493389A (en)
CA (1) CA2001199C (en)
GB (1) GB8825074D0 (en)
WO (1) WO1990004820A1 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2289957A (en) * 1994-05-28 1995-12-06 British Aerospace Adaptive control systems
US5992383A (en) * 1996-05-28 1999-11-30 U.S. Philips Corporation Control unit having a disturbance predictor, a system controlled by such a control unit, an electrical actuator controlled by such a control unit, and throttle device provided with such an actuator

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61234409A (en) * 1985-04-10 1986-10-18 Hitachi Ltd Flow control system

Also Published As

Publication number Publication date
EP0440720A1 (en) 1991-08-14
CA2001199A1 (en) 1990-04-26
WO1990004820A1 (en) 1990-05-03
JPH04501476A (en) 1992-03-12
AU4493389A (en) 1990-05-14
GB8825074D0 (en) 1988-11-30

Similar Documents

Publication Publication Date Title
EP0515518B1 (en) Repetitive sound or vibration phenomena cancellation arrangement with multiple sensors and actuators
CA2159589C (en) Frequency domain adaptive control system
US5691893A (en) Adaptive control system
Roure et al. The remote microphone technique for active noise control
US7343016B2 (en) Linear independence method for noninvasive on-line system identification/secondary path modeling for filtered-X LMS-based active noise control systems
Bouchard New recursive-least-squares algorithms for nonlinear active control of sound and vibration using neural networks
Kuo et al. Review of DSP algorithms for active noise control
EP0660958B1 (en) Sampled-data filter with low delay
Das et al. Nonlinear active noise control for headrest using virtual microphone control
EP0789900B1 (en) Reduction of computational burden of adaptively updating control filter(s) in active systems
CA2001199C (en) Signal processing means for identifying systems subject to periodic disturbances
EP0373943A2 (en) Vibration control system
US5652770A (en) Sampled-data filter with low delay
JP3732227B2 (en) Adaptive control system for controlling repetitive events
Douglas Analysis of the multiple-error and block least-mean-square adaptive algorithms
Chen et al. Effects of multiple secondary paths on convergence properties in active noise control systems with LMS algorithm
Kewley et al. Feedforward control using the higher‐harmonic, time‐averaged gradient descent algorithm
EP0694234B1 (en) Control system for periodic disturbances
JP3220418B2 (en) Signal control system and signal control method
JPH04505221A (en) Active sound/vibration cancellation device
Scott et al. Robust multivariable control of a double beam cantilever smart structure
JP2827603B2 (en) Active noise control device
Gilson et al. IV methods for closed-loop system identification
JPH0816176A (en) Noise cancel system
JPH05232969A (en) Vehicle vibration control device

Legal Events

Date Code Title Description
EEER Examination request
MKLA Lapsed