WO2015125275A1 - 加速度検出器および能動騒音制御装置 - Google Patents

加速度検出器および能動騒音制御装置 Download PDF

Info

Publication number
WO2015125275A1
WO2015125275A1 PCT/JP2014/054174 JP2014054174W WO2015125275A1 WO 2015125275 A1 WO2015125275 A1 WO 2015125275A1 JP 2014054174 W JP2014054174 W JP 2014054174W WO 2015125275 A1 WO2015125275 A1 WO 2015125275A1
Authority
WO
WIPO (PCT)
Prior art keywords
acceleration
detected
detection
signal
acceleration signal
Prior art date
Application number
PCT/JP2014/054174
Other languages
English (en)
French (fr)
Inventor
敦仁 矢野
Original Assignee
三菱電機株式会社
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 三菱電機株式会社 filed Critical 三菱電機株式会社
Priority to US15/113,688 priority Critical patent/US9858914B2/en
Priority to CN201480075355.5A priority patent/CN106258000B/zh
Priority to JP2016503855A priority patent/JP5976251B2/ja
Priority to DE112014006393.2T priority patent/DE112014006393T5/de
Priority to PCT/JP2014/054174 priority patent/WO2015125275A1/ja
Publication of WO2015125275A1 publication Critical patent/WO2015125275A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
    • G10K11/1785—Methods, e.g. algorithms; Devices
    • G10K11/17853—Methods, e.g. algorithms; Devices of the filter
    • G10K11/17854—Methods, e.g. algorithms; Devices of the filter the filter being an adaptive filter
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P15/00—Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration
    • G01P15/18—Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration in two or more dimensions
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
    • G10K11/1787—General system configurations
    • G10K11/17879—General system configurations using both a reference signal and an error signal
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
    • G10K11/178—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound by electro-acoustically regenerating the original acoustic waves in anti-phase
    • G10K11/1787—General system configurations
    • G10K11/17879—General system configurations using both a reference signal and an error signal
    • G10K11/17883—General system configurations using both a reference signal and an error signal the reference signal being derived from a machine operating condition, e.g. engine RPM or vehicle speed
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/10—Applications
    • G10K2210/128—Vehicles
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/10—Applications
    • G10K2210/128—Vehicles
    • G10K2210/1282—Automobiles
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K2210/00—Details of active noise control [ANC] covered by G10K11/178 but not provided for in any of its subgroups
    • G10K2210/10—Applications
    • G10K2210/129—Vibration, e.g. instead of, or in addition to, acoustic noise

Definitions

  • the present invention relates to an acceleration detector that detects acceleration of a vibration noise source, and an active noise control device that silences vibration noise by active noise control with reference to the detected acceleration.
  • a technique is disclosed in which an acceleration sensor is installed in the vibration noise source and the detected acceleration is referred to when the vibration noise generated by the vibration of a vibration noise source such as a machine or a vehicle is silenced by active noise control (for example, Patent Document 1).
  • the vibration of the object does not necessarily have a single movement direction, and may have vibration elements having a plurality of different movement directions such as vertical and horizontal directions.
  • the motion direction of the vibration to be detected does not necessarily coincide with the axis of the acceleration sensor (that is, the acceleration detection direction) for design reasons.
  • a technique is disclosed in which the deviation of the axis of the acceleration sensor with respect to the direction of vibration to be detected is measured in advance, and the output value of the acceleration sensor is corrected according to the deviation (for example, Patent Document 2). reference).
  • Patent Document 3 discloses a technique for measuring the direction of motion using an arctangent function from an acceleration ratio of each axis using an acceleration sensor that detects accelerations of two orthogonal axes. Yes.
  • a complex structure such as a machine or a vehicle has an infinite number of vibration modes, and the motion direction and frequency in each vibration mode vary. Therefore, the acceleration signal of each axis includes vibration elements that mainly contribute to noise. There is a problem that various vibration components other than the above are mixed as disturbances.
  • the present invention has been made to solve the above-described problems. First, it detects the direction of vibration that contributes to noise or its acceleration, and secondly effectively silences the vibration noise. Objective.
  • the acceleration detector receives acceleration signals of respective axes observed in a three-axis or two-axis orthogonal coordinate system orthogonal to each other, and outputs a detected acceleration signal converted into a target detection direction. Based on the detected declination that defines the target detection direction in the Cartesian coordinate system, the direction vector setting unit that sets the direction vector, and the inner product of the acceleration signal of each axis and the direction vector are calculated to obtain the detected acceleration signal. And an inner product calculation unit to be obtained.
  • the active noise control device outputs a control signal for canceling the vibration noise from the vibration noise source from the speaker, an error signal between the vibration noise and the control signal detected by the microphone, and the vibration noise detected by the acceleration sensor.
  • a direction vector setting unit for setting a direction vector based on a detection declination defining a target detection direction in a three-axis or two-axis orthogonal coordinate system orthogonal to each other;
  • An inner product calculation unit that calculates an inner product of an acceleration signal of each axis of the orthogonal coordinate system input from the acceleration sensor and a direction vector, and obtains a detected acceleration signal obtained by converting the acceleration signal of each axis into a target detection direction.
  • the control signal is controlled based on the error signal and the detected acceleration signal.
  • the direction vector is set in accordance with the detected deflection angle, and the detected acceleration signal is obtained by calculating the inner product with the three-axis or two-axis acceleration signal orthogonal thereto. It is possible to reduce the sensitivity to acceleration other than the motion direction of the vibration to be detected and to accurately detect the acceleration of the target vibration. In addition, since the direction vector is set according to the detected deflection angle, even if the movement direction of the target vibration changes, the acceleration of the target vibration is detected promptly based on the new detection deflection angle. be able to.
  • vibration that does not contribute to noise is obtained by obtaining the acceleration of a vibration element having a high contribution to noise as a detected acceleration signal among the vibration elements of the vibration noise source, and performing active noise control based on this acceleration signal. It is possible to suppress the influence of elements and enhance the silencing effect.
  • 6 is a graph for explaining the operation of the acceleration detector according to the first embodiment. 6 is a three-dimensional graph showing an example of detection sensitivity with respect to acceleration in each direction of 0 ⁇ a ′ ⁇ 2 ⁇ and 0 ⁇ b ′ ⁇ of the acceleration detector when the detection declination is determined.
  • Embodiment 1 FIG.
  • an acceleration detector that accurately detects an acceleration in a motion direction when the motion direction of the vibration to be detected is previously known will be described.
  • FIG. 1 is a block diagram showing a configuration of an acceleration detector 10 according to Embodiment 1 of the present invention.
  • an acceleration detector that detects acceleration in a target direction from three-axis accelerations orthogonal to each other will be described.
  • those skilled in the art use the same means as in the first embodiment to perform orthogonality.
  • An acceleration detector that detects acceleration in a target direction from biaxial acceleration can be configured, and such a configuration is naturally included in the scope of the present invention.
  • an acceleration detector 10 includes a direction vector setting unit 11 that sets a direction vector u according to detection deviation angles a and b given in advance, a direction vector u and triaxial acceleration signals x (n), y.
  • the inner product calculation unit 12 calculates an inner product of (n) and z (n) and obtains an acceleration (hereinafter referred to as a detected acceleration signal) v (n) in the motion direction of the vibration to be detected.
  • the acceleration detector 10 is composed of, for example, a CPU (Central Processing Unit) (not shown), and when the CPU executes a program stored in an internal memory, the direction vector setting unit 11 and the inner product calculation unit 12 are executed. As a function.
  • the acceleration signals x (n), y (n), z (n) input to the acceleration detector 10 are obtained from the output of the acceleration sensor, for example.
  • FIG. 2 is a graph for explaining the operation of the acceleration detector 10, and the detected deflection angles a and b and the direction vector indicated by a three-dimensional coordinate system constituted by three axes of x, y and z orthogonal to each other.
  • u the observed values x (n), y (n), z (n) of the acceleration of each axis, the observed value vector p, and the detected acceleration v (n) are illustrated.
  • the detected deflection angles a and b input to the direction vector setting unit 11 are angles that determine the direction of vibration to be detected in the three-dimensional coordinate system.
  • the detected deflection angle a is an angle formed by an orthogonal projection vector of the direction vector u on the xy plane with the x axis
  • the detected deflection angle b is an angle formed by the direction vector u with the z axis. It is.
  • these detection declinations a and b values determined in accordance with the motion direction of the vibration to be detected are input. However, in the description here, 0 ⁇ a ⁇ 2 ⁇ and 0 ⁇ b ⁇ are assumed for simplicity.
  • the detected declination angles a and b are determined according to the direction of motion of the vibration element having the greatest contribution to noise.
  • the direction vector setting unit 11 sets the direction vector u as shown in Expression (1) using the input detected deflection angles a and b, and outputs the result to the inner product calculation unit 12.
  • the inner product calculation unit 12 includes an observation value vector p having elements of acceleration signals x (n), y (n), and z (n) corresponding to each of the three orthogonal xyz axes shown in Expression (2). Then, the inner product with the direction vector u output from the direction vector setting unit 11 is calculated according to the equation (3), and the detected acceleration signal v (n) that is the acceleration in the motion direction indicated by the direction vector u is calculated.
  • the acceleration signals x (n), y (n), and z (n) are obtained from the output of the acceleration sensor, for example. If the obtained acceleration signal has only two orthogonal xy axes, the detected deflection angle b is limited to ⁇ / 2, so that the inner product calculation unit 12 detects the detected acceleration signal v (n) according to the equation (4). Is calculated.
  • the sign of the detected acceleration signal v (n) calculated by the above equation (3) or (4) is positive when the displacement is in the same direction as the direction vector u, and negative when the displacement is in the reverse direction. Become. If the direction vector u ′ of the equation (5) obtained by turning the direction vector u in the opposite direction is used instead of the direction vector u, the detected acceleration signal v (n) whose sign is reversed is obtained. If the sign of the detected acceleration signal v (n) is not important, either of the direction vectors u and u ′ may be used.
  • the detection sensitivity is calculated by Equation (6).
  • the direction vector setting unit 11 sets the direction vector according to the input detected declination, and the inner product calculation unit 12 is orthogonal to each other.
  • the inner product calculation unit 12 is orthogonal to each other.
  • the direction vector setting unit 11 sets the direction vector according to the input detection declination, so that the motion direction of the vibration to be detected is set. Even when is changed, the acceleration of the target vibration can be detected promptly by inputting a new detected deflection angle.
  • Embodiment 2 the direction vector is set according to the input detected deflection angle. However, in the second embodiment, the direction vector is determined based on the triaxial or biaxial acceleration signals orthogonal to each other. The structure to perform is demonstrated.
  • FIG. 4 is a block diagram illustrating a configuration of the acceleration detector 20 according to the second embodiment.
  • the acceleration detector 20 includes a direction vector setting unit 11, an inner product calculation unit 12, and a detected declination a from three-axis acceleration signals x (n), y (n), z (n) orthogonal to each other. , B and the detection declination setting unit 21.
  • the direction vector setting unit 11 and the inner product calculation unit 12 denoted by the same reference numerals as those in FIG. 1 are the same components as the direction vector setting unit 11 and the inner product calculation unit 12 of the acceleration detector 10 according to the first embodiment. Therefore, explanation is omitted.
  • the detection declination setting unit 21 receives triaxial or biaxial acceleration signals orthogonal to each other, and determines the detection declinations a and b so that the average power of the detection acceleration signal v (n) is maximized.
  • the vibration element having the maximum power is considered to have the highest contribution to noise. Therefore, if the detection declinations a and b are determined as described above, the acceleration of the vibration element can be detected. In other applications, the acceleration detector 20 can be used when detection of the vibration element having the maximum power is important.
  • the detected deflection angles a such that the average power of the detected acceleration signal v (n) is maximized with respect to the triaxial acceleration signals x (n), y (n), z (n) orthogonal to each other.
  • a method for obtaining b will be described.
  • the gradient of E [v 2 (n)] with respect to a and b is obtained.
  • the detected deflection angles a and b may be sequentially updated by Expressions (7) and (8) using a value obtained by multiplying this by a predetermined constant as an update amount.
  • ⁇ is a predetermined constant that satisfies ⁇ > 0, ⁇ 2 x , ⁇ 2 y , and ⁇ 2 z are average powers of x (n), y (n), and z (n), respectively, and ⁇ xy is x (N) and y (n), ⁇ xz is the correlation coefficient of x (n) and z (n), ⁇ yz is the correlation coefficient of y (n) and z (n), and the following equations (9), ( 10).
  • Equation (11) represents an average operation.
  • Equation (11) a method of calculating these average values using a moving average as shown in Equation (11) is known.
  • P x (n) is a calculated value of ⁇ 2 x
  • ⁇ is a smoothing coefficient that satisfies 0 ⁇ ⁇ 1.
  • the above formulas (7) and (8) are obtained by calculating the gradient of the detected deflection angles a (n) and b (n) with respect to E [v 2 (n)] and multiplying this by an arbitrary constant ⁇ .
  • the update amount is added to each original value to update the value.
  • a (n) and b (n) are recursively updated so that E [v 2 (n)] becomes larger, and finally a [in which E [v 2 (n)] becomes maximum. It converges to (n) and b (n).
  • the detected deflection angle b is limited to ⁇ / 2, and therefore the detected deflection angle setting unit 21 is detected by the equation (12).
  • the declination a (n) may be updated.
  • Equation (13) which is an analytical solution of the detected deflection angle a that maximizes E [v 2 (n)] without sequentially updating, is obtained.
  • the detection deviation angle a may be determined by using this.
  • the direction vector u is determined by the direction vector setting unit 11 in accordance with the detected declination angles a (n) and b (n) determined by any one of the methods described above, whereby the detected acceleration signal output from the inner product calculation unit 12 v (n) is the acceleration of the vibration element having the maximum power.
  • the detection deflection angle setting unit 21 determines the detection deflection angle according to the triaxial or biaxial acceleration signals orthogonal to each other. Thus, even if the motion direction of the vibration to be detected is unknown, the acceleration can be detected.
  • the detection deviation angle setting unit 21 determines the detection deviation angle based on the average power and the correlation coefficient of the acceleration signal. It is possible to suppress the fluctuation of the detected declination due to instantaneous disturbance and to stably detect the motion direction and acceleration of the target vibration.
  • the detection deviation angle setting unit 21 determines a detection deviation angle that maximizes the average amplitude of the detection acceleration signal. The direction and acceleration of the vibration element having the maximum power can be obtained.
  • the detected deviation angle calculated based on the average power and the correlation coefficient of the acceleration signal is based on the gradient of the magnitude of the detected acceleration signal with respect to the detected deviation angle. Since the update amount is set in this way, and the detection declination is updated sequentially, the motion direction and acceleration of the target vibration can be accurately detected even if steady disturbance is mixed in the acceleration signal. Will be able to.
  • FIG. 5 is a block diagram illustrating a configuration of the active noise control device 30 according to the third embodiment.
  • the active noise control device 30 includes an acceleration detector 20, a FIR (Finite Impulse Response) filter 31, a secondary path filter 32, and a LMS (Least Mean Square) processing unit 33.
  • the active noise control device 30 is connected to an acceleration sensor 101, an error detection microphone 102, and a speaker 103 installed in a vibration noise source 100 that generates vibration noise.
  • the acceleration detector 20 is the acceleration detector described in the second embodiment, and the same reference numerals as those in FIG. 4 are given. If appropriate detection deviation angles a and b can be input from the outside, the acceleration detector 20 may be replaced with the acceleration detector 10 described in the first embodiment.
  • the active noise control device 30 includes three xyz acceleration signals x (n), y (n), and z (n) orthogonal to each other that are output from the acceleration sensor 101 installed in the vibration noise source 100 that generates vibration noise. And outputs a control signal d (n) for canceling noise (indicated by a broken line arrow in FIG. 5) emitted from the vibration noise source 100.
  • the control signal d (n) is output from the speaker 103 and interferes with noise in space. If the output sound of the speaker 103 is correctly in the opposite phase in the space, the noise is canceled out and a silencing effect is obtained.
  • the error detection microphone 102 detects this cancellation error and outputs an error signal e (n).
  • the error signal e (n) is input to the active noise control device 30, and the active noise control device 30 controls the control signal d (n) so that the cancellation error is reduced.
  • the acceleration signals x (n), y (n), and z (n) output from the acceleration sensor 101 are input to the acceleration detector 20 of the active noise control device 30.
  • the acceleration detector 20 detects the acceleration of the vibration element having the maximum power in the vibration of the vibration noise source 100 and outputs it as a detected acceleration signal v (n). .
  • the detected acceleration signal v (n) is input to the FIR filter 31.
  • the FIR filter 31 filters the detected acceleration signal by the filter coefficient that is held, and outputs a control signal d (n).
  • the secondary path filter 32 filters the detected acceleration signal v (n) with a filter coefficient simulating the acoustic transfer characteristic from the speaker 103 to the error detection microphone 102, and outputs v ′ (n) to the LMS processing unit 33. To do.
  • the LMS processing unit 33 Based on the error signal e (n) input from the error detection microphone 102 and the detected acceleration signal v ′ (n) filtered by the secondary path filter 32, the LMS processing unit 33 reduces the cancellation error.
  • the filter coefficient of the FIR filter 31 is updated.
  • a series of signal processing composed of the FIR filter 31, the secondary path filter 32, and the LMS processing unit 33 is known as a Filtered X-LMS system. For example, “Active Noise Control” written by Masaharu Nishimura et al. (Corona, 2006) Issued on July 7, pp. 74-76).
  • the present invention is not necessarily limited to the active noise control device using the Filtered X-LMS system, and the filter coefficient of the active noise control device using another adaptive algorithm or the FIR filter 31 is optimized in advance.
  • An active noise control device that does not perform updating may be used. It is considered that such an active noise control device may have a configuration different from that of the above-described active noise control device 30, but as long as it has the configuration of the acceleration detectors 10 and 20 of the present invention, they are also within the scope of the present invention. Contained within.
  • the active noise control device 30 of the third embodiment the acceleration of the vibration element having a high contribution to the noise among the vibration elements of the vibration noise source 100 is detected and active based on this.
  • noise control it is possible to suppress the influence of vibration elements that do not contribute to noise, and to enhance the silencing effect.
  • Embodiment 4 In the third embodiment, the detection deviation angle is determined so that the power of the detected acceleration signal is maximized. However, in the fourth embodiment, the detection deviation angle is set so that the power of the error signal detected by the microphone is minimized. The structure which determines is described.
  • FIG. 6 is a block diagram showing the configuration of the active noise control device 40 according to the fourth embodiment.
  • the active noise control device 40 includes an acceleration detector 41, an FIR filter 31, a secondary path filter 32, and an LMS processing unit 33.
  • the active noise control device 40 is connected to an acceleration sensor 101, an error detection microphone 102, and a speaker 103 installed in the vibration noise source 100.
  • the acceleration detector 41 includes a detection deviation angle setting unit 42, a direction vector setting unit 11, and an inner product calculation unit 12.
  • FIG. 6 the same components as those in FIGS. 1 to 5 in the first to third embodiments are denoted by the same reference numerals, and detailed description thereof is omitted.
  • a detection declination setting unit 42 which is a part of the acceleration detector 41 includes acceleration signals x (n), y (n), z (n) output from the acceleration sensor 101 and errors output from the error detection microphone 102.
  • the signal e (n) and the filter coefficient h N held by the FIR filter 31 are input.
  • the detection deviation angle setting unit 42 obtains the gradient of the power e 2 (n) of the error signal e (n) with respect to the detection deviation angles a and b, and based on this, the detection deviation angle a is reduced so that the power e 2 (n) decreases.
  • B are determined, both are updated and output to the direction vector setting unit 11.
  • c M is a transfer characteristic of the secondary path corresponding to the coefficient sequence of the secondary path filter 32.
  • X M (n) is an acceleration signal x (n)
  • y M (n) is an acceleration signal y (n)
  • z M (n) is an acceleration signal z (n). It is a vector as an element, and is expressed as in the following formulas (16), (17), (18), and (19).
  • T is a transposition symbol.
  • the detection declination setting unit 42 recursively detects the recursions of the detection declinations a and b by multiplying the purchase represented by the equations (20) and (21) by a predetermined constant ⁇ (0 ⁇ ). By updating the deflection angles a and b, the detected deviation angles a and b are converged so as to minimize the power e 2 (n) of the error signal e (n).
  • the detected deflection angle b is limited to ⁇ / 2, and therefore the equation (20) is replaced by the following equation (22). Is done.
  • ⁇ is a delay time and satisfies ⁇ ⁇ 0.
  • the detection deviation angle setting unit 42 is input from the acceleration signal of each axis input from the acceleration sensor 101 and the error detection microphone 102. Based on the error signal power gradient with respect to the detected deflection angle calculated based on the error signal and the filter coefficient of the FIR filter 31, an update amount of the detected deflection angle is determined so as to reduce the power of the error signal, As a result, the detection declination is sequentially updated, so that the detection declination that minimizes the power of the error signal can be automatically selected to enhance the silencing effect.
  • the detection deviation angle setting unit 42 sets the update amount of the detection deviation angle based on at least the transfer characteristic of the secondary path from the speaker 103 to the error detection microphone 102.
  • Embodiment 5 FIG.
  • the detection deflection angle is updated so that the power of the error signal is minimized.
  • the detection deflection angle update and the FIR filter coefficient update interfere with each other. A configuration for adjusting both update processes will be described.
  • the active noise control device of the fifth embodiment has the same configuration as the active noise control device 40 of FIG. 6 in the drawing, the fifth embodiment will be described with reference to FIG.
  • the LMS processing unit 33 and the detection declination setting unit 42 both detect the filter coefficient or detection of the FIR filter 31 so that the power e 2 (n) of the error signal e (n) decreases.
  • the declination angles a and b are performed, if these are updated at the same time, interference may occur with each other, and the power e 2 (n) of the error signal may not decrease smoothly, and the silencing effect may be impaired.
  • each update process is performed alternately, or a constant value that is multiplied by one update amount is set smaller than the other, and one of the updates is slowed down. Can be stabilized.
  • update of the detection declination and coefficient update of the FIR filter are alternately performed, or a constant multiplied by one update amount is set as the other update amount.
  • any combination of each embodiment, any component of each embodiment can be modified, or any component can be omitted in each embodiment. .
  • the acceleration detector according to the present invention can accurately detect the acceleration of a target vibration, the acceleration detector silences the vibration noise of a vibration noise source having an infinite number of vibration modes such as machines and vehicles. Suitable for use in noise control devices.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • General Physics & Mathematics (AREA)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)
  • Circuit For Audible Band Transducer (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

　加速度検出器１０は、互いに直交したｘｙｚの３軸の直交座標系において目的の検出方向を定義した検出偏角ａ，ｂを基に方向ベクトルｕを設定する方向ベクトル設定部１１と、その直交座標系において観測された各軸の加速度信号ｘ（ｎ），ｙ（ｎ），ｚ（ｎ）と方向ベクトルｕの内積を計算して検出加速度信号ｖ（ｎ）を求める内積計算部１２とを備える。

Description

加速度検出器および能動騒音制御装置
　この発明は、振動騒音源の加速度を検出する加速度検出器、および検出された加速度を参照して能動騒音制御によって振動騒音を消音する能動騒音制御装置に関するものである。
　機械、車両などの振動騒音源の振動によって生じる振動騒音を能動騒音制御によって消音する際、この振動騒音源に加速度センサを設置し、検出された加速度を参照する技術が開示されている（例えば、特許文献１参照）。
　しかし、物体の振動は、必ずしも単一の運動方向だけでなく、例えば縦、横など異なる複数の運動方向の振動要素を持つ場合がある。このとき、設計上の都合により、必ずしも検出したい振動の運動方向と加速度センサの軸線（即ち、加速度の検出方向）とが一致しない場合が考えられる。このような場合には、検出したい振動の運動方向に対する加速度センサの軸線のずれを事前に計測し、ずれに応じて加速度センサの出力値を補正する技術が開示されている（例えば、特許文献２参照）。
　上記特許文献２の技術は、騒音に寄与する振動の運動方向と加速度センサの軸線とのずれが一定でなく、条件によって変化する場合には対応できない。この問題に関して、例えば特許文献３では、直交する２軸の加速度を検出する加速度センサを用いて、各軸の加速度の比率から逆正接関数（アークタンジェント）により運動方向を計測する技術が開示されている。
特開平６－１１０４７４号公報 特開２０１０－１１２７９３号公報 特開２００９－１２８１６４号公報
　機械または車両のような複雑な構造物には無数の振動モードが存在し、各振動モードでの運動方向および周波数は様々であるため、各軸の加速度信号には、騒音に主として寄与する振動要素以外の様々な振動成分が外乱として混入するという課題があった。
　このような場合には、外乱の振動成分の影響を抑え、騒音に主として寄与する振動要素のみを精度よく検出する必要があるが、上記特許文献１～３のいずれにもこのような技術は開示されていない。例えば特許文献３の技術は、目的とする振動と外乱との合成結果に基づき運動方向と加速度を時々刻々と計測するのみであるので、この技術を能動騒音制御に用いた場合、外乱の影響によって消音効果が減じるという問題が生じる。
　この発明は、上記のような課題を解決するためになされたもので、第一に騒音に寄与する振動の方向またはその加速度を検出すること、第二に効果的に振動騒音を消音することを目的とする。
　この発明に係る加速度検出器は、互いに直交した３軸または２軸の直交座標系において観測された各軸の加速度信号が入力され、目的の検出方向に変換された検出加速度信号を出力するものであって、直交座標系において目的の検出方向を定義した検出偏角を基に、方向ベクトルを設定する方向ベクトル設定部と、各軸の加速度信号と方向ベクトルの内積を計算して検出加速度信号を求める内積計算部とを備えるものである。
　この発明に係る能動騒音制御装置は、振動騒音源からの振動騒音を消音する制御信号をスピーカから出力させ、マイクが検出する振動騒音と制御信号との誤差信号、および加速度センサが検出する振動騒音源の加速度信号が入力されるものであって、互いに直交した３軸または２軸の直交座標系において目的の検出方向を定義した検出偏角を基に、方向ベクトルを設定する方向ベクトル設定部と、加速度センサから入力された直交座標系の各軸の加速度信号と方向ベクトルの内積を計算して、各軸の加速度信号を目的の検出方向に変換した検出加速度信号を求める内積計算部とを備え、誤差信号および検出加速度信号に基づいて制御信号を制御するようにしたものである。
　この発明によれば、検出偏角に応じて方向ベクトルを設定し、これと直交する３軸または２軸の加速度信号との内積を計算して検出加速度信号を求めるようにしたので、検出目的とする振動の運動方向以外の加速度に対する感度を減少させ、目的とする振動の加速度を精度よく検出することができる。また、検出偏角に応じて方向ベクトルを設定するようにしたので、目的とする振動の運動方向が変化した場合でも、新たな検出偏角を基に、すみやかに目的の振動の加速度を検出することができる。
　この発明によれば、振動騒音源の振動要素のうち、騒音への寄与度の高い振動要素の加速度を検出加速度信号として求め、これに基づいて能動騒音制御を行うことにより、騒音に寄与しない振動要素の影響を抑え、消音効果を高めることができる。
この発明の実施の形態１に係る加速度検出器の構成を示すブロック図である。 実施の形態１に係る加速度検出器の動作を説明するためのグラフである。 検出偏角を定めたときの、加速度検出器の、０≦ａ’＜２π、０≦ｂ’＜πの各方向の加速度に対する検出感度の例を示した３次元グラフである。 この発明の実施の形態２に係る加速度検出器の構成を示すブロック図である。 この発明の実施の形態３に係る能動騒音制御装置の構成を示すブロック図である。 この発明の実施の形態４に係る能動騒音制御装置の構成を示すブロック図である。
　以下、この発明をより詳細に説明するために、この発明を実施するための形態について、添付の図面に従って説明する。
実施の形態１．
　実施の形態１では、検出したい振動の運動方向が予め明らかである場合において、当該運動方向の加速度を精度よく検出する加速度検出器について述べる。
　図１は、この発明の実施の形態１に係る加速度検出器１０の構成を示すブロック図である。実施の形態１では、例として、互いに直交する３軸の加速度から目的方向の加速度を検出する加速度検出器について述べるが、当業者は本実施の形態１と同様の手段を用いることにより、直交する２軸の加速度から目的方向の加速度を検出する加速度検出器を構成することが可能であり、当然このような構成も本発明の範囲内に含まれる。
　図１において、加速度検出器１０は、予め与えられている検出偏角ａ，ｂに従い方向ベクトルｕを設定する方向ベクトル設定部１１と、方向ベクトルｕと３軸の加速度信号ｘ（ｎ），ｙ（ｎ），ｚ（ｎ）の内積を計算し、検出目的とする振動の運動方向の加速度（以下、検出加速度信号）ｖ（ｎ）を求める内積計算部１２とによって構成されている。
　この加速度検出器１０は、例えば、不図示のＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）で構成されており、このＣＰＵが内部メモリに格納されたプログラムを実行することによって、方向ベクトル設定部１１および内積計算部１２としての機能を実現する。
　また、加速度検出器１０に入力される加速度信号ｘ（ｎ），ｙ（ｎ），ｚ（ｎ）は、例えば加速度センサの出力などから得られる。
　図２は、加速度検出器１０の動作を説明するためのグラフであり、互いに直交するｘ，ｙ，ｚの３軸によって構成される３次元座標系で示される検出偏角ａ，ｂ、方向ベクトルｕ、各軸の加速度の観測値ｘ（ｎ），ｙ（ｎ），ｚ（ｎ）とその観測値ベクトルｐ、および検出加速度ｖ（ｎ）を例示する。
　方向ベクトル設定部１１に入力される検出偏角ａ，ｂは、３次元座標系において検出目的とする振動の方向を定める角度である。即ち、検出偏角ａは、方向ベクトルｕのｘｙ平面上への直交射影ベクトルがｘ軸との間になす角度であり、検出偏角ｂは、方向ベクトルｕがｚ軸との間になす角度である。これら検出偏角ａ，ｂは、検出目的とする振動の運動方向に従って定められた値が入力されるものとする。ただし、ここでの説明では簡便のため０≦ａ＜２π、０≦ｂ＜πとする。
　例えば能動騒音制御においては、最も騒音への寄与度の大きい振動要素の運動方向に従って検出偏角ａ，ｂが定められる。
　方向ベクトル設定部１１は、入力された検出偏角ａ，ｂを用いて、方向ベクトルｕを式（１）のように設定し、内積計算部１２へ出力する。
Figure JPOXMLDOC01-appb-I000001
　内積計算部１２は、式（２）で示される、直交するｘｙｚの３軸のそれぞれに対応する加速度信号ｘ（ｎ），ｙ（ｎ），ｚ（ｎ）を要素とする観測値ベクトルｐと、方向ベクトル設定部１１から出力された方向ベクトルｕとの内積を式（３）に従って計算し、方向ベクトルｕが示す運動方向の加速度である検出加速度信号ｖ（ｎ）を算出する。
Figure JPOXMLDOC01-appb-I000002

Figure JPOXMLDOC01-appb-I000003
　ここで、加速度信号ｘ（ｎ），ｙ（ｎ），ｚ（ｎ）は、例えば加速度センサの出力などから得る。なお、得られる加速度信号が直交するｘｙの２軸のみである場合は、検出偏角ｂはπ／２に限定されるので、内積計算部１２は式（４）に従って検出加速度信号ｖ（ｎ）を算出する。
Figure JPOXMLDOC01-appb-I000004
　以上の式（３）または式（４）によって算出された検出加速度信号ｖ（ｎ）の符号は、方向ベクトルｕと同じ方向の変位である場合が正、逆方向の変位である場合が負となる。
　もし方向ベクトルｕの代わりに、これを逆方向に転じた式（５）の方向ベクトルｕ’を用いれば、単純に符号が逆転した検出加速度信号ｖ（ｎ）が得られる。
　検出加速度信号ｖ（ｎ）の符号が重要でない場合は、方向ベクトルｕ，ｕ’のいずれを用いてもよい。
Figure JPOXMLDOC01-appb-I000005
　図３は、検出偏角（ａ，ｂ）＝（１．４π，０．２π）［ｒａｄ］と定めたときの、加速度検出器１０の、０≦ａ’＜２π、０≦ｂ’＜πの各方向の加速度に対する検出感度の例を示した３次元グラフである。
　ここで、検出感度は式（６）によって計算している。
（検出感度）＝（検出加速度）÷（偏角ａ’，ｂ’方向の加速度）　（６）
　図３に示されるように、検出感度は（ａ’，ｂ’）＝（１．４π，０．２π）［ｒａｄ］で極大値１、（ａ’，ｂ’）＝（０．４π，０．８π）［ｒａｄ］で極小値－１ｔなり、これら以外の方向は検出感度の絶対値が低い。
　このことから、方向ベクトルｕおよびその逆の－ｕ方向の加速度の検出感度が高く、これら以外の方向の検出感度が低くなることが明らかである。
　以上説明したように、実施の形態１の加速度検出器１０によれば、方向ベクトル設定部１１が、入力された検出偏角に応じて方向ベクトルを設定し、内積計算部１２が、互いに直交する３軸または２軸の加速度信号と方向ベクトルとの内積を計算して検出加速度信号を求めるようにしたことで、検出目的とする振動の運動方向以外の加速度に対する感度を減少させ、目的とする振動の加速度を精度よく検出することができるようになる。
　また、実施の形態１の加速度検出器１０によれば、方向ベクトル設定部１１が、入力された検出偏角に応じて方向ベクトルを設定するようにしたことで、検出目的とする振動の運動方向が変化した場合でも、新たな検出偏角が入力されることで、すみやかに目的の振動の加速度を検出することができるようになる。
実施の形態２．
　上記実施の形態１では、入力された検出偏角に応じて方向ベクトルを設定するようにしていたが、実施の形態２では、互いに直交する３軸または２軸の加速度信号に基づき方向ベクトルを決定する構成を説明する。
　図４は、実施の形態２に係る加速度検出器２０の構成を示すブロック図である。図４において、加速度検出器２０は、方向ベクトル設定部１１と、内積計算部１２と、互いに直交する３軸の加速度信号ｘ（ｎ），ｙ（ｎ），ｚ（ｎ）から検出偏角ａ，ｂを定める検出偏角設定部２１とによって構成される。
　ここで、図１と共通の符号を付した方向ベクトル設定部１１および内積計算部１２は、上記実施の形態１の加速度検出器１０の方向ベクトル設定部１１および内積計算部１２と同一の構成要素であるので説明は省略する。
　次に、加速度検出器２０の動作を説明する。
　検出偏角設定部２１は、互いに直交する３軸または２軸の加速度信号を受け、検出加速度信号ｖ（ｎ）の平均パワーが最大となるように検出偏角ａ，ｂを定める。例えば能動騒音制御においては、振動騒音の音圧は振動する物体の加速度の大きさに比例するので、最大のパワーを持つ振動要素が最も騒音への寄与度が高いと考えられる。したがって上述のように検出偏角ａ，ｂを定めれば、当該振動要素の加速度を検出することができるようになる。またその他の用途においても、最大のパワーを持つ振動要素の検出が重要である場合に、この加速度検出器２０を用いることができる。
　ここで、互いに直交する３軸の加速度信号ｘ（ｎ），ｙ（ｎ），ｚ（ｎ）に対して、検出加速度信号ｖ（ｎ）の平均パワーが最大となるような検出偏角ａ，ｂを求める方法を説明する。
　検出偏角ａ，ｂに対して検出加速度信号ｖ（ｎ）の平均パワーＥ［ｖ2（ｎ）］を最大とするには、ａ，ｂに対するＥ［ｖ2（ｎ）］の勾配を求め、これに所定の定数を乗じた値を更新量として式（７）および式（８）により逐次的に検出偏角ａ，ｂを更新すればよい。
Figure JPOXMLDOC01-appb-I000006

Figure JPOXMLDOC01-appb-I000007
　ただし、μはμ＞０を満たす所定の定数であり、σ2 x，σ2 y，σ2 zはそれぞれｘ（ｎ），ｙ（ｎ），ｚ（ｎ）の平均パワー、σxyはｘ（ｎ）とｙ（ｎ）、σxzはｘ（ｎ）とｚ（ｎ）、σyzはｙ（ｎ）とｚ（ｎ）の相関係数であって、以下の式（９）、（１０）で表される。
Figure JPOXMLDOC01-appb-I000008

Figure JPOXMLDOC01-appb-I000009
　上式において、Ｅ［・］は平均操作を表す。これら平均値は、実際の装置では例えば式（１１）のような移動平均を用いて算出する方法が知られている。
Figure JPOXMLDOC01-appb-I000010
　ここで、Ｐx（ｎ）はσ2 xの算出値であり、λは０＜λ＜１を満たす平滑化係数である。
　上式（７）および上式（８）は、Ｅ［ｖ2（ｎ）］に対する検出偏角ａ（ｎ），ｂ（ｎ）の勾配を求め、これに任意の定数μを乗じた値を更新量としてもとの値それぞれに加算し、値を更新する式となっている。その結果、ａ（ｎ），ｂ（ｎ）はＥ［ｖ2（ｎ）］がより大きくなるよう再帰的に更新され、最終的にＥ［ｖ2（ｎ）］が最大となるようなａ（ｎ），ｂ（ｎ）に収束する。
　なお、加速度センサなどから得られる加速度信号が直交するｘｙの２軸のみである場合は、検出偏角ｂはπ／２に限定されるので、検出偏角設定部２１は式（１２）によって検出偏角ａ（ｎ）を更新すればよい。
Figure JPOXMLDOC01-appb-I000011
　また、得られる加速度信号が直交するｘｙの２軸のみの場合は、逐次更新を行わずＥ［ｖ2（ｎ）］を最大にする検出偏角ａの解析解である、式（１３）を用いて検出偏角ａを定めてもよい。
Figure JPOXMLDOC01-appb-I000012
　上述のいずれかの方法で定められた検出偏角ａ（ｎ），ｂ（ｎ）に従って、方向ベクトル設定部１１により方向ベクトルｕが定められることで、内積計算部１２から出力される検出加速度信号ｖ（ｎ）は最大のパワーを持つ振動要素の加速度となる。
　以上説明したように、実施の形態２の加速度検出器２０によれば、検出偏角設定部２１が、互いに直交する３軸または２軸の加速度信号に応じて検出偏角を定めるようにしたことで、検出目的とする振動の運動方向が不明であっても、その加速度を検出することができるようになる。
　また、実施の形態２の加速度検出器２０によれば、検出偏角設定部２１が、加速度信号の平均パワーおよび相関係数を基に検出偏角を定めるようにしたことで、加速度信号への瞬間的な外乱による検出偏角のぶれを抑制し、目的とする振動の運動方向と加速度を安定して検出することができるようになる。
　また、実施の形態２の加速度検出器２０によれば、検出偏角設定部２１が、検出加速度信号の平均的な振幅の大きさが最大となるような検出偏角を定めるようにしたことで、最大のパワーをもつ振動要素の方向と加速度が得られるようになる。
　また、実施の形態２の加速度検出器２０によれば、加速度信号の平均パワーおよび相関係数に基づいて算出される検出偏角について、当該検出偏角に対する検出加速度信号の大きさの勾配を基に更新量を定め、これによって検出偏角を逐次的に更新するようにしたことで、加速度信号に定常的な外乱が混入しても、目的とする振動の運動方向と加速度を精度よく検出することができるようになる。
実施の形態３．
　実施の形態３では、この発明に係る加速度検出器を用いた能動騒音制御装置の構成例を説明する。
　図５は、実施の形態３に係る能動騒音制御装置３０の構成を示すブロック図である。図５において、能動騒音制御装置３０は、加速度検出器２０と、ＦＩＲ（Ｆｉｎｉｔｅ　Ｉｍｐｕｌｓｅ　Ｒｅｓｐｏｎｓｅ）フィルタ３１と、二次経路フィルタ３２と、ＬＭＳ（Ｌｅａｓｔ　Ｍｅａｎ　Ｓｑｕａｒｅ）処理部３３とによって構成されている。また、能動騒音制御装置３０は、振動騒音を発生する振動騒音源１００に設置された加速度センサ１０１と、誤差検出マイク１０２と、スピーカ１０３とに接続されている。
　ここで、加速度検出器２０は、上記実施の形態２で説明した加速度検出器であり、図４と共通の符号を付した。もし適切な検出偏角ａ，ｂが外部より入力可能であれば、加速度検出器２０を上記実施の形態１で説明した加速度検出器１０に置き換えてもよい。
　次に、能動騒音制御装置３０の外部動作を説明する。
　能動騒音制御装置３０は、振動騒音を発する振動騒音源１００に設置された加速度センサ１０１から出力される互いに直交するｘｙｚの３軸の加速度信号ｘ（ｎ），ｙ（ｎ），ｚ（ｎ）を受け、振動騒音源１００から発せられる騒音（図５に破線矢印で示す）を打ち消すための制御信号ｄ（ｎ）を出力する。制御信号ｄ（ｎ）はスピーカ１０３から出力され、空間上で騒音に干渉する。スピーカ１０３の出力音が空間上で正しく逆相音になっていれば、騒音が打ち消され消音効果が得られる。誤差検出マイク１０２は、この打消し誤差を検出し、誤差信号ｅ（ｎ）を出力する。誤差信号ｅ（ｎ）は能動騒音制御装置３０に入力され、能動騒音制御装置３０は、打消し誤差が減少するように制御信号ｄ（ｎ）を制御する。
　次に、能動騒音制御装置３０の内部動作を説明する。
　加速度センサ１０１の出力する加速度信号ｘ（ｎ），ｙ（ｎ），ｚ（ｎ）は、能動騒音制御装置３０の加速度検出器２０に入力される。加速度検出器２０は、上記実施の形態２で説明したように、振動騒音源１００の振動の中で、最大のパワーをもつ振動要素の加速度を検出し、検出加速度信号ｖ（ｎ）として出力する。
　検出加速度信号ｖ（ｎ）はＦＩＲフィルタ３１に入力され、ＦＩＲフィルタ３１は、保持するフィルタ係数によってこれを濾波し、制御信号ｄ（ｎ）を出力する。また、二次経路フィルタ３２は、スピーカ１０３から誤差検出マイク１０２までの音響伝達特性を模擬したフィルタ係数で検出加速度信号ｖ（ｎ）を濾波し、ｖ’（ｎ）をＬＭＳ処理部３３へ出力する。
　ＬＭＳ処理部３３は、誤差検出マイク１０２から入力された誤差信号ｅ（ｎ）と、二次経路フィルタ３２によって濾波された検出加速度信号ｖ’（ｎ）に基づき、打消し誤差が減少するようにＦＩＲフィルタ３１のフィルタ係数を更新する。
　これらＦＩＲフィルタ３１、二次経路フィルタ３２、ＬＭＳ処理部３３によって構成される一連の信号処理は、Ｆｉｌｔｅｒｅｄ　Ｘ－ＬＭＳシステムとして知られ、例えば、西村正治他著「アクティブノイズコントロール」（コロナ社、２００６年７月７日発行、ｐ．７４～７６）に開示されている。
　ただし、本発明は必ずしもＦｉｌｔｅｒｅｄ　Ｘ－ＬＭＳシステムを用いる能動騒音制御装置に限定されるものではなく、その他の適応アルゴリズムを用いる能動騒音制御装置、あるいはＦＩＲフィルタ３１のフィルタ係数が予め最適化されていて更新を行わないような能動騒音制御装置などであってもよい。そのような能動騒音制御装置は、上述の能動騒音制御装置３０と異なる構成となる場合があると考えられるが、本発明の加速度検出器１０，２０の構成を有する限り、それらも本発明の範囲内に含まれる。
　以上説明したように、実施の形態３の能動騒音制御装置３０によれば、振動騒音源１００の振動要素のうち、騒音への寄与度の高い振動要素の加速度を検出し、これに基づいて能動騒音制御を行うことにより、騒音に寄与しない振動要素の影響を抑え、消音効果を高めることができる。
実施の形態４．
　上記実施の形態３では、検出加速度信号のパワーが最大となるよう検出偏角を定めるようにしていたが、実施の形態４では、マイクで検出する誤差信号のパワーが最小となるよう検出偏角を定める構成を説明する。
　図６は、実施の形態４に係る能動騒音制御装置４０の構成を示すブロック図である。図６において、能動騒音制御装置４０は、加速度検出器４１と、ＦＩＲフィルタ３１と、二次経路フィルタ３２と、ＬＭＳ処理部３３とによって構成されている。また、能動騒音制御装置４０は、振動騒音源１００に設置された加速度センサ１０１と、誤差検出マイク１０２と、スピーカ１０３とに接続されている。また、加速度検出器４１は、検出偏角設定部４２と、方向ベクトル設定部１１と、内積計算部１２とによって構成されている。図６において、上記実施の形態１～３の図１～図５と同一の構成要素については同一符号を付し、詳細な説明は省略する。
　次に、能動騒音制御装置４０の動作を説明する。
　加速度検出器４１の一部である検出偏角設定部４２には、加速度センサ１０１の出力する加速度信号ｘ（ｎ），ｙ（ｎ），ｚ（ｎ）と、誤差検出マイク１０２の出力する誤差信号ｅ（ｎ）と、ＦＩＲフィルタ３１の保持するフィルタ係数ｈNとが入力される。
　検出偏角設定部４２は、検出偏角ａ，ｂに対する誤差信号ｅ（ｎ）のパワーｅ2（ｎ）の勾配を求め、これに基づきパワーｅ2（ｎ）が減少するよう検出偏角ａ，ｂの更新量を定め、両者を更新して方向ベクトル設定部１１へ出力する。
　検出偏角ａ，ｂに対する誤差信号パワーｅ2（ｎ）の勾配は、式（１４）および式（１５）によって表される。
Figure JPOXMLDOC01-appb-I000013

Figure JPOXMLDOC01-appb-I000014
　ここで、ｈk（ｋ＝０，１，・・・，Ｎ－１）はＦＩＲフィルタ３１のフィルタ係数であり、ｃMは二次経路フィルタ３２の係数列に相当する二次経路の伝達特性を示すインパルス応答列であり、ｘM（ｎ）は加速度信号ｘ（ｎ）を、ｙM（ｎ）は加速度信号ｙ（ｎ）を、ｚM（ｎ）は加速度信号ｚ（ｎ）をそれぞれ要素とするベクトルであり、以下の式（１６）、（１７）、（１８）、（１９）のように表される。
Figure JPOXMLDOC01-appb-I000015

Figure JPOXMLDOC01-appb-I000016

Figure JPOXMLDOC01-appb-I000017

Figure JPOXMLDOC01-appb-I000018
　なお、Ｔは転置記号である。
　検出偏角設定部４２では、式（２０）および式（２１）によって表される購買に所定の定数μ（０＜μ）を乗じて検出偏角ａ，ｂの更新量とし、再帰的に検出偏角ａ，ｂを更新することで誤差信号ｅ（ｎ）のパワーｅ2（ｎ）を最小とするような検出偏角ａ，ｂに収束させる。
Figure JPOXMLDOC01-appb-I000019

Figure JPOXMLDOC01-appb-I000020
　また、加速度センサ１０１から得られる加速度信号が直交するｘｙの２軸のみである場合は、検出偏角ｂはπ／２に限定されるので、式（２０）は以下の式（２２）で代替される。
Figure JPOXMLDOC01-appb-I000021
　また、二次経路が単純な遅延で近似できる場合は、上式（２０）および上式（２１）を以下の式（２３）および式（２４）のように単純化することができ、この場合は二次経路の伝達特性は計算式に明示的に現われなくなる。
Figure JPOXMLDOC01-appb-I000022

Figure JPOXMLDOC01-appb-I000023
　なお、τは遅延時間であり、τ≧０を満たす。
　以上説明したように、実施の形態４の能動騒音制御装置４０によれば、検出偏角設定部４２が、加速度センサ１０１から入力される各軸の加速度信号と、誤差検出マイク１０２から入力される誤差信号と、ＦＩＲフィルタ３１のフィルタ係数とに基づいて算出される、検出偏角に対する誤差信号のパワーの勾配を基に、誤差信号のパワーを減少させるように検出偏角の更新量を定め、これによって検出偏角を逐次的に更新するようにしたことで、誤差信号のパワーが最小となるような検出偏角を自動的に選択し、消音効果を高めることができるようになる。
　また、実施の形態４の能動騒音制御装置４０によれば、検出偏角設定部４２が、少なくともスピーカ１０３から誤差検出マイク１０２までの二次経路の伝達特性に基づいて検出偏角の更新量を定めるようにしたことで、二次経路が単純な遅延で近似できないような複雑な伝達特性を持っていたとしても、誤差信号のパワーが最小となるような検出偏角を自動的に選択し、消音効果を高めることができるようになる。
実施の形態５．
　上記実施の形態４では、誤差信号のパワーが最小となるように検出偏角を更新するようにしていたが、実施の形態５では、検出偏角の更新とＦＩＲフィルタの係数更新とが互いに干渉しないよう、両者の更新処理を調整する構成を説明する。
　実施の形態５の能動騒音制御装置は、図面の上では図６の能動騒音制御装置４０と同一の構成であるため、図６を援用して実施の形態５を説明する。
　図６に示したように、ＬＭＳ処理部３３および検出偏角設定部４２は、共に誤差信号ｅ（ｎ）のパワーｅ2（ｎ）が減少するように、ＦＩＲフィルタ３１のフィルタ係数、もしくは検出偏角ａ，ｂを行うものであるが、これらを同時に更新すると、互いに干渉を生じて誤差信号のパワーｅ2（ｎ）が順調に減少せず、消音効果が損なわれる場合がある。
　このような問題に対しては、それぞれの更新処理を交互に行うか、または一方の更新量に乗じる定数の値を、他方より小さく設定して、一方の更新を緩慢にすることで、消音効果を安定させることができる。
　以上説明したように、実施の形態４の能動騒音制御装置４０によれば、検出偏角の更新とＦＩＲフィルタの係数更新とを交互に行うか、一方の更新量に乗じる定数を他方の更新量に乗じる定数より小さく設定して一方の更新を緩慢にすることで、それぞれの更新処理が互いに干渉して消音効果が損なわれることを防ぐことができる。
　なお、本発明はその発明の範囲内において、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。
　以上のように、この発明に係る加速度検出器は、目的とする振動の加速度を精度よく検出できるので、機械および車両のように無数の振動モードが存在する振動騒音源の振動騒音を消音する能動騒音制御装置に用いるのに適している。
　１０，２０，４１　加速度検出器、１１　方向ベクトル設定部、１２　内積計算部、２１，４２　検出偏角設定部、３０，４０　能動騒音制御装置、３１　ＦＩＲフィルタ、３２　二次経路フィルタ、３３　ＬＭＳ処理部、１００　振動騒音源、１０１　加速度センサ、１０２　誤差検出マイク、１０３　スピーカ。

Claims (10)

  1. 　互いに直交した３軸または２軸の直交座標系において観測された各軸の加速度信号が入力され、目的の検出方向に変換された検出加速度信号を出力する加速度検出器であって、
    　前記直交座標系において前記目的の検出方向を定義した検出偏角を基に、方向ベクトルを設定する方向ベクトル設定部と、
    　前記各軸の加速度信号と前記方向ベクトルの内積を計算して前記検出加速度信号を求める内積計算部とを備えることを特徴とする加速度検出器。
  2. 　前記各軸の加速度信号に基づいて前記検出偏角を定め、前記方向ベクトル設定部へ通知する検出偏角設定部を備えることを特徴とする請求項１記載の加速度検出器。
  3. 　前記検出偏角設定部は、少なくとも前記各軸の加速度信号の平均パワーおよび軸間の加速度信号の相関係数に基づいて前記検出偏角を定めることを特徴とする請求項２記載の加速度検出器。
  4. 　前記検出偏角設定部は、前記検出加速度信号の平均パワーが最大となるように前記検出偏角を定めることを特徴とする請求項３記載の加速度検出器。
  5. 　前記検出偏角設定部は、少なくとも前記各軸の加速度信号の平均パワーおよび軸間の加速度信号の相関係数に基づいて算出される前記検出偏角に対する前記検出加速度信号の平均パワーの勾配を基に、前記検出加速度信号の平均パワーを増大させるように前記検出偏角の更新量を定め、当該更新量により前記検出偏角を逐次的に更新していき、前記検出加速度信号の平均パワーが最大となるときの前記検出偏角を求めることを特徴とする請求項４記載の加速度検出器。
  6. 　振動騒音源からの振動騒音を消音する制御信号をスピーカから出力させ、マイクが検出する前記振動騒音と前記制御信号との誤差信号、および加速度センサが検出する前記振動騒音源の加速度信号が入力される能動騒音制御装置であって、
    　互いに直交した３軸または２軸の直交座標系において目的の検出方向を定義した検出偏角を基に、方向ベクトルを設定する方向ベクトル設定部と、
    　前記加速度センサから入力された前記直交座標系の各軸の加速度信号と前記方向ベクトルの内積を計算して、前記各軸の加速度信号を前記目的の検出方向に変換した検出加速度信号を求める内積計算部とを備え、
    　前記誤差信号および前記検出加速度信号に基づいて前記制御信号を制御することを特徴とする能動騒音制御装置。
  7. 　前記検出加速度信号を濾波して前記制御信号を制御するフィルタと、
    　少なくとも前記各軸の加速度信号、前記誤差信号、および前記フィルタのフィルタ係数に基づいて算出される、前記検出偏角に対する前記誤差信号のパワーの勾配を基に、前記誤差信号のパワーを減少させるように前記検出偏角の更新量を定め、当該更新量により前記検出偏角を逐次的に更新していき、前記誤差信号のパワーが最小となるときの前記検出偏角を求めて前記方向ベクトル設定部へ通知する検出偏角設定部とを備えることを特徴とする請求項６記載の能動騒音制御装置。
  8. 　前記検出偏角設定部は、前記各軸の加速度信号、前記誤差信号、および前記フィルタのフィルタ係数に加えて、前記スピーカから前記マイクへの伝達特性に基づいて、前記検出偏角の更新量を定めることを特徴とする請求項７記載の能動騒音制御装置。
  9. 　前記誤差信号および前記検出加速度信号に基づいて前記フィルタのフィルタ係数を更新するＬＭＳ処理部を備え、
    　前記検出偏角設定部と前記ＬＭＳ処理部は、前記検出偏角の更新処理と前記フィルタ係数の更新処理を交互に行うことを特徴とする請求項７記載の能動騒音制御装置。
  10. 　前記誤差信号および前記検出加速度信号に基づいて前記フィルタのフィルタ係数を更新するＬＭＳ処理部を備え、
    　前記検出偏角設定部および前記ＬＭＳ処理部のいずれか一方は、前記検出偏角の更新量および前記フィルタ係数の更新量のいずれか一方に乗じる定数を他方に乗じる定数より小さい値にして緩慢に更新処理を行うことを特徴とする請求項７記載の能動騒音制御装置。
PCT/JP2014/054174 2014-02-21 2014-02-21 加速度検出器および能動騒音制御装置 WO2015125275A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
US15/113,688 US9858914B2 (en) 2014-02-21 2014-02-21 Acceleration detector and active noise-control device
CN201480075355.5A CN106258000B (zh) 2014-02-21 2014-02-21 加速度检测器和有源噪声控制装置
JP2016503855A JP5976251B2 (ja) 2014-02-21 2014-02-21 加速度検出器および能動騒音制御装置
DE112014006393.2T DE112014006393T5 (de) 2014-02-21 2014-02-21 Beschleunigungssensor und aktive Geräuschsteuervorrichtung
PCT/JP2014/054174 WO2015125275A1 (ja) 2014-02-21 2014-02-21 加速度検出器および能動騒音制御装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2014/054174 WO2015125275A1 (ja) 2014-02-21 2014-02-21 加速度検出器および能動騒音制御装置

Publications (1)

Publication Number Publication Date
WO2015125275A1 true WO2015125275A1 (ja) 2015-08-27

Family

ID=53877809

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/054174 WO2015125275A1 (ja) 2014-02-21 2014-02-21 加速度検出器および能動騒音制御装置

Country Status (5)

Country Link
US (1) US9858914B2 (ja)
JP (1) JP5976251B2 (ja)
CN (1) CN106258000B (ja)
DE (1) DE112014006393T5 (ja)
WO (1) WO2015125275A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110830862A (zh) * 2019-10-10 2020-02-21 广东思派康电子科技有限公司 一种自适应降噪的降噪耳机
US11232778B1 (en) * 2020-07-22 2022-01-25 Bose Corporation Systems and methods for detecting divergence in an adaptive system
CN113436599A (zh) * 2021-05-31 2021-09-24 新东方教育科技集团有限公司 信号处理方法、装置、存储介质及电子设备
CN114684159A (zh) * 2022-03-21 2022-07-01 潍柴动力股份有限公司 一种车辆质量估算方法、装置、电子设备和存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08509823A (ja) * 1993-04-27 1996-10-15 ノイズキャンセレイション テクノロジイーズ インコーポレーテッド 能動的音響及び振動制御のための単一及び多重チャネルブロック適応方法と装置
JP2009508133A (ja) * 2005-09-16 2009-02-26 ヴェーテーイー テクノロジーズ オサケユキチュア 加速度の微小機械式計測方法および微小機械式加速度センサー
JP5222897B2 (ja) * 2009-05-14 2013-06-26 ハーマン インターナショナル インダストリーズ インコーポレイテッド 適応スピーカ選択を用いたアクティブノイズ制御のためのシステム

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06110474A (ja) * 1992-09-30 1994-04-22 Matsushita Electric Ind Co Ltd 消音装置
JP2009128164A (ja) 2007-11-22 2009-06-11 Alps Electric Co Ltd 加速度・角速度・磁気方位検出用複合センサ及びこれを用いた装置
JP2011121534A (ja) * 2009-12-14 2011-06-23 Honda Motor Co Ltd 能動型騒音制御装置
JP5388997B2 (ja) 2010-11-25 2014-01-15 三菱電機株式会社 速度計測装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08509823A (ja) * 1993-04-27 1996-10-15 ノイズキャンセレイション テクノロジイーズ インコーポレーテッド 能動的音響及び振動制御のための単一及び多重チャネルブロック適応方法と装置
JP2009508133A (ja) * 2005-09-16 2009-02-26 ヴェーテーイー テクノロジーズ オサケユキチュア 加速度の微小機械式計測方法および微小機械式加速度センサー
JP5222897B2 (ja) * 2009-05-14 2013-06-26 ハーマン インターナショナル インダストリーズ インコーポレイテッド 適応スピーカ選択を用いたアクティブノイズ制御のためのシステム

Also Published As

Publication number Publication date
US9858914B2 (en) 2018-01-02
JP5976251B2 (ja) 2016-08-23
JPWO2015125275A1 (ja) 2017-03-30
DE112014006393T5 (de) 2016-11-17
CN106258000A (zh) 2016-12-28
US20170004816A1 (en) 2017-01-05
CN106258000B (zh) 2019-05-21

Similar Documents

Publication Publication Date Title
JP6650570B2 (ja) 能動型騒音低減装置
JP5976251B2 (ja) 加速度検出器および能動騒音制御装置
JP6073454B2 (ja) 能動振動騒音制御装置
Beijen et al. Disturbance feedforward control for active vibration isolation systems with internal isolator dynamics
WO2014068624A1 (ja) 能動振動騒音制御装置
JP6180680B2 (ja) 能動振動騒音制御装置
JP2008250131A (ja) 能動型騒音制御装置
Van De Ridder et al. Active vibration isolation feedback control for coriolis mass-flow meters
JP5962842B2 (ja) 空間安定装置、空間安定方法、及び空間安定プログラム
Fleming et al. Spatial system identification of a simply supported beam and a trapezoidal cantilever plate
WO2019227353A1 (en) Method and device for estimating a direction of arrival
JP6433630B2 (ja) 雑音除去装置、エコーキャンセリング装置、異音検出装置および雑音除去方法
CN116438597A (zh) 用于适应所估计的次级路径的系统和方法
KR20130082327A (ko) 가속도 센서가 구비된 오브젝트의 진동 제어 방법 및 장치
JP5141351B2 (ja) 能動型騒音制御装置
JP6627563B2 (ja) 音響センサ及び音波検出方法
Krause et al. Active Noise Control Experiment Minimising Radiation of Active Energy
Abir et al. A novel accelerometer based feedback concept for improving machine dynamic performance
Yuan Self-learning active noise control
Kobayashi et al. A control source structure of single loudspeaker and rear sound interference for inexpensive active noise control
JP2012168282A (ja) 能動型振動騒音抑制装置
JP2018205059A (ja) 慣性計測装置およびその加速度データ処理構造
JP2008250130A (ja) 能動型騒音制御装置
JP2007140249A (ja) 能動騒音低減装置
JP2013011696A (ja) 能動型振動騒音抑制装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14882971

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2016503855

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 15113688

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 112014006393

Country of ref document: DE

122 Ep: pct application non-entry in european phase

Ref document number: 14882971

Country of ref document: EP

Kind code of ref document: A1