JPH06324686A - Method and device for active control - Google Patents

Method and device for active control

Info

Publication number
JPH06324686A
JPH06324686A JP5110170A JP11017093A JPH06324686A JP H06324686 A JPH06324686 A JP H06324686A JP 5110170 A JP5110170 A JP 5110170A JP 11017093 A JP11017093 A JP 11017093A JP H06324686 A JPH06324686 A JP H06324686A
Authority
JP
Japan
Prior art keywords
signal
filter
waveform
error
coefficient
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.)
Withdrawn
Application number
JP5110170A
Other languages
Japanese (ja)
Inventor
Takehiko Hiei
武彦 樋江井
Tetsuya Matsuura
哲哉 松浦
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.)
Daikin Industries Ltd
Original Assignee
Daikin Industries Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Daikin Industries Ltd filed Critical Daikin Industries Ltd
Priority to JP5110170A priority Critical patent/JPH06324686A/en
Publication of JPH06324686A publication Critical patent/JPH06324686A/en
Withdrawn legal-status Critical Current

Links

Abstract

PURPOSE:To provide algorithm which is equal in calculation quantity to LMS algorithm and can reduce a convergence error. CONSTITUTION:An inverted sound signal which is opposite in phase from the noise in a duct 2 and equal in amplitude is generated by an adaptive FIR filter 41 and an inverted sound corresponding to the inverted sound signal is radiated by a speaker 3b for sound elimination; and a monitor microphone 3c detects a reduced sound level at an observation point and the filter coefficient of the adaptive FIR filter 41 is updated so as to reduce the reduced sound level. A coefficient vector for updating the filter coefficient of the adaptive FIR filter 41 is calculated and added at intervals of detection sampling by the microphone and when the frequency of the addition matches the number of filter coefficients of the adaptive FIR filter 41, the filter coefficient of the adaptive FIR filter 41 is updated on the basis of the added coefficient vector.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、広帯域の変動波形信号
に対して該信号に基いて生成された付加信号を作用させ
るアクティブ制御方法及びその装置の改良に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an active control method and an apparatus therefor which act on a wide-band fluctuating waveform signal with an additional signal generated based on the signal.

【0002】[0002]

【従来の技術】一般に、この種のアクティブ制御装置に
は、例えば、ダクト内空間等の音波の伝播通路の非定常
的な広帯域の騒音に対して、これとは逆位相で且つ同振
幅の反転音を発生させて消音を行うアクティブ消音装置
がある。このアクティブ消音装置としては、従来より、
適応型FIR(Finite Inpulse Response)フィルタ
を使用し、検出マイクロフォンから出力される空間の騒
音信号に対して前記適応型FIRフィルタにより逆位相
で同振幅の反転音信号を生成した後、該反転音信号をス
ピーカに出力して反転音を空間に放射すると共に、この
空間の所定観測点に配置したモニタマイクロフォンによ
り該観測点にて騒音と前記スピーカから放射された反転
音との合成音を検出して、該合成音を低減音レベルとし
て入力し、該低減音レベルを小さくするように適応型F
IRフィルタのフィルタ係数を逐次更新することによ
り、前記観測点周辺の音圧レベルを低減するようにした
ものが知られている(例えば、電子情報通信学会の技術
研究報告の1988年EA−88−29参照)。
2. Description of the Related Art Generally, an active control device of this type has an inversion of the opposite phase and the same amplitude with respect to an unsteady broadband noise in a sound wave propagation path such as a space in a duct. There is an active muffling device that generates a sound to muffle the sound. As this active silencer,
An adaptive FIR (Finite Impulse Response) filter is used, and after generating a reversal sound signal of the same amplitude in antiphase by the adaptive FIR filter with respect to the noise signal in the space output from the detection microphone, the reversal sound signal Is output to the speaker to emit a reverse sound to the space, and a monitor microphone arranged at a predetermined observation point in this space detects a synthesized sound of the noise and the reverse sound emitted from the speaker at the observation point. , The synthetic sound is input as a reduced sound level, and the adaptive F is used to reduce the reduced sound level.
It is known that the sound pressure level around the observation point is reduced by successively updating the filter coefficient of the IR filter (for example, 1988 EA-88- of the Technical Research Report of the Institute of Electronics, Information and Communication Engineers). 29).

【0003】[0003]

【発明が解決しようとする課題】ところで、上述したア
クティブ消音装置のようなアクティブ制御装置にあって
は、適応型FIRフィルタのフィルタ係数を逐次更新す
る適応アルゴリズムとしてLMSアルゴリズムが広く用
いられている。しかしながら、このLMSアルゴリズム
は、瞬時予測値を用いて検出マイクロフォン及びモニタ
マイクロフォンでの騒音検出のサンプリング毎に全ての
フィルタ係数の更新を行うようにしているので、各フィ
ルタ係数の更新の度に収束誤差が生じており、このた
め、経時的に制御が収束された状態で、収束点に達した
状況にあっても、この収束誤差のために、適応制御の安
定且つ適切な収束状態を得ることができなかった。
By the way, in the active control device such as the above-mentioned active silencer, the LMS algorithm is widely used as an adaptive algorithm for successively updating the filter coefficient of the adaptive FIR filter. However, this LMS algorithm updates all the filter coefficients at every sampling of noise detection by the detection microphone and the monitor microphone using the instantaneous predicted value, so that the convergence error is updated every time each filter coefficient is updated. Therefore, even if the control is converged over time and the convergence point is reached, due to this convergence error, stable and appropriate convergence of adaptive control can be obtained. could not.

【0004】また、LMSアルゴリズムよりも収束誤差
を低減することができるアルゴリズムも種々存在する
が、これらのアルゴリズムはLMSアルゴリズムに比べ
て計算量が多いものである。例えば、フィルタ係数の数
Mに対して、LMSアルゴリズムは(2M+1)回の乗
算で済むのに対し、収束誤差を低減することができるア
ルゴリズムの代表的なものとしてのRLSアルゴリズム
にあっては{3M(3+M)/2}回の乗算が必要であ
る。このように、これまでの収束誤差を低減することが
できるアルゴリズムでは計算量が多いことから上述した
アクティブ消音のような実時間制御アプリケーションに
おいては実用的ではない。
There are various algorithms that can reduce the convergence error as compared with the LMS algorithm, but these algorithms have a larger amount of calculation than the LMS algorithm. For example, with respect to the number M of filter coefficients, the LMS algorithm only requires (2M + 1) times of multiplication, whereas the RLS algorithm as a typical algorithm that can reduce the convergence error is {3M (3 + M) / 2} multiplications are required. As described above, since the algorithm capable of reducing the convergence error so far has a large amount of calculation, it is not practical in the real-time control application such as the active muffling described above.

【0005】本発明は、斯かる点に鑑みてなされたもの
であり、計算量がLMSアルゴリズムと同等であって且
つ収束誤差を低減することができるアルゴリズムを提供
することを目的とする。
The present invention has been made in view of the above problems, and an object of the present invention is to provide an algorithm having a calculation amount equivalent to that of the LMS algorithm and capable of reducing a convergence error.

【0006】[0006]

【課題を解決するための手段】上記の目的を達成するた
め、本発明では、LMSアルゴリズムにおいて、複数回
サンプリング毎に適応型FIRフィルタのフィルタ係数
を更新するようにし、その更新されるフィルタ係数は複
数回サンプリング毎のフィルタ係数ベクトルを加算する
ことによって得るようにした。具体的に請求項1記載の
発明は、図1に示すように、波形検出手段(3a)により変
動波形信号を検出し、該波形検出手段(3a)から出力され
る波形信号に対して、該波形信号に基いて生成された制
御信号をデジタルフィルタ(41)から出力し、該制御信号
に対応した作動信号を信号出力手段(3b)から出力した
後、該信号出力手段(3b)が出力した作動信号と前記変動
波形信号との誤差を誤差検出手段(3c)により検出して、
該誤差検出手段(3c)から出力される誤差信号と前記波形
検出手段(3a)から出力される波形信号とを適応アルゴリ
ズム実行手段(44)が受け、該適応アルゴリズム実行手段
(44)により前記誤差が小さくなるように前記デジタルフ
ィルタ(41)のフィルタ係数を更新するようにしたアクテ
ィブ制御方法を前提としている。そして、前記適応アル
ゴリズム実行手段(44)に備えられたサンプリング回数カ
ウント手段(5) により、前記誤差検出手段(3c)からの誤
差信号と波形検出手段(3a)からの波形信号との検出サン
プリング回数をカウントし、このサンプリング回数が所
定値に達したときにサンプリング回数カウント手段(5)
からフィルタ更新信号を出力すると共に、前記適応アル
ゴリズム実行手段(44)に備えられた前記サンプリング回
数カウント手段(5) からのフィルタ更新信号が受信可能
な係数ベクトル加算手段(6) により、前記誤差検出手段
(3c)からの誤差信号と波形検出手段(3a)からの波形信号
との検出サンプリング毎に前記デジタルフィルタ(41)の
フィルタ係数を更新するための係数ベクトルを計算して
加算し、係数ベクトル加算手段(6) が、前記サンプリン
グ回数カウント手段(5) からのフィルタ更新信号を受け
たときに、この加算された係数ベクトルに基いて前記デ
ジタルフィルタ(41)のフィルタ係数を更新させるように
したものである。
In order to achieve the above object, the present invention updates the filter coefficient of the adaptive FIR filter in the LMS algorithm every sampling a plurality of times, and the updated filter coefficient is It is obtained by adding the filter coefficient vectors for each sampling. Specifically, as shown in FIG. 1, the invention according to claim 1 detects a fluctuation waveform signal by the waveform detection means (3a), and detects the fluctuation waveform signal with respect to the waveform signal output from the waveform detection means (3a). The control signal generated based on the waveform signal is output from the digital filter (41), the operation signal corresponding to the control signal is output from the signal output means (3b), and then the signal output means (3b) is output. Detecting the error between the operating signal and the fluctuation waveform signal by the error detecting means (3c),
The adaptive algorithm execution means (44) receives the error signal output from the error detection means (3c) and the waveform signal output from the waveform detection means (3a), and the adaptive algorithm execution means
It is premised on the active control method in which the filter coefficient of the digital filter (41) is updated so that the error is reduced by (44). Then, by the sampling number counting means (5) provided in the adaptive algorithm executing means (44), the number of detection samplings of the error signal from the error detecting means (3c) and the waveform signal from the waveform detecting means (3a) The sampling frequency counting means (5) when the sampling frequency reaches a predetermined value.
Output the filter update signal from the adaptive algorithm execution means (44) the sampling frequency counting means (5) filter update signal can be received by the coefficient vector addition means (6), the error detection means
(3c) error signal and waveform signal from the waveform detection means (3a) detection for each sampling sampling coefficient vector for updating the filter coefficient of the digital filter (41) and added, coefficient vector addition Means for updating the filter coefficient of the digital filter (41) based on the added coefficient vector when the means (6) receives the filter update signal from the sampling number counting means (5) Is.

【0007】請求項2記載の発明は、前記請求項1記載
の方法を用いたアクティブ制御装置であって、具体的に
は、変動波形信号を検出して波形信号を出力する波形検
出手段(3a)と、該波形検出手段(3a)からの波形信号に対
し、該波形信号に基いて生成された制御信号を出力する
デジタルフィルタ(41)と、該デジタルフィルタ(41)より
出力された制御信号を受けて該制御信号に対応した作動
信号を出力する信号出力手段(3b)と、該信号出力手段(3
b)が出力した作動信号と前記変動波形信号との誤差を検
出して誤差信号を出力する誤差検出手段(3c)と、該誤差
検出手段(3c)から出力された誤差信号と前記波形検出手
段(3a)からの波形信号とを受けて前記誤差が小さくなる
ように前記デジタルフィルタ(41)のフィルタ係数を更新
する適応アルゴリズム実行手段(44)とを備えたアクティ
ブ制御装置を前提としている。そして、前記適応アルゴ
リズム実行手段(44)に、前記誤差検出手段(3c)からの誤
差信号と波形検出手段(3a)からの波形信号との検出サン
プリング回数をカウントし、このサンプリング回数が所
定値に達したときにフィルタ更新信号を出力するサンプ
リング回数カウント手段(5) と、該サンプリング回数カ
ウント手段(5) からのフィルタ更新信号が受信可能とさ
れ、前記誤差検出手段(3c)からの誤差信号と波形検出手
段(3a)からの波形信号との検出サンプリング毎に前記デ
ジタルフィルタ(41)のフィルタ係数を更新するための係
数ベクトルを計算して加算し、前記サンプリング回数カ
ウント手段(5) からのフィルタ更新信号を受けたとき
に、この加算された係数ベクトルに基いて前記デジタル
フィルタ(41)のフィルタ係数を更新させる係数ベクトル
加算手段(6) とを備えさせるような構成としている。
The invention according to claim 2 is an active control apparatus using the method according to claim 1, and more specifically, a waveform detecting means (3a) for detecting a fluctuation waveform signal and outputting the waveform signal. ), And a waveform signal from the waveform detecting means (3a), a digital filter (41) for outputting a control signal generated based on the waveform signal, and a control signal output from the digital filter (41) Signal output means (3b) for receiving an operation signal corresponding to the control signal, and the signal output means (3
error detecting means (3c) for detecting an error between the operation signal output by b) and the fluctuation waveform signal and outputting an error signal, and the error signal output from the error detecting means (3c) and the waveform detecting means It is premised on an active control device provided with an adaptive algorithm executing means (44) for receiving the waveform signal from (3a) and updating the filter coefficient of the digital filter (41) so as to reduce the error. Then, the adaptive algorithm executing means (44) counts the number of detection samplings of the error signal from the error detecting means (3c) and the waveform signal from the waveform detecting means (3a), and the sampling number becomes a predetermined value. Sampling number counting means (5) for outputting a filter updating signal when reaching, and the filter updating signal from the sampling number counting means (5) can be received, and the error signal from the error detecting means (3c) The coefficient vector for updating the filter coefficient of the digital filter (41) is calculated and added for each detection sampling with the waveform signal from the waveform detection means (3a), and the filter from the sampling frequency counting means (5) Coefficient vector adding means (6) for updating the filter coefficient of the digital filter (41) based on the added coefficient vector when receiving the update signal It is configured to be equipped with and.

【0008】請求項3記載の発明は、前記請求項2記載
のアクティブ制御装置において、サンプリング回数カウ
ント手段(5) を、サンプリング回数がデジタルフィルタ
(41)のフィルタ係数の数に一致したときにフィルタ更新
信号を出力するような構成としている。
According to a third aspect of the present invention, in the active control device according to the second aspect, the sampling number counting means (5) is provided with a sampling frequency digital filter.
The filter update signal is output when the number of filter coefficients in (41) matches.

【0009】請求項4記載の発明は、前記請求項2記載
の装置を消音装置に利用したものであって、具体的に
は、図2に示すように、空間(21)の騒音を検出して騒音
信号を出力する検出マイクロフォン(3a)と、該検出マイ
クロフォン(3a)からの騒音信号に対して逆位相で同振幅
の反転音信号を生成する適応型FIRフィルタ(41)と、
該適応型FIRフィルタ(41)より出力された反転音信号
を受けて該反転音信号に基く反転音を前記空間に放射す
るスピーカ(3b)と、前記空間(21)の所定観測点に配置さ
れ、該観測点の低減音レベルを検出して低減音信号をフ
ィードバック出力するモニタマイクロフォン(3c)と、前
記検出マイクロフォン(3a)からの騒音信号及びモニタマ
イクロフォン(3c)からフィードバックされた低減音信号
を入力し、騒音信号及び低減音信号により前記観測点周
辺の音圧レベルを低減するように前記適応型FIRフィ
ルタ(41)のフィルタ係数を更新する適応アルゴリズム実
行回路(44)とを備えたアクティブ制御装置を前提として
いる。そして、前記適応アルゴリズム実行回路(44)に、
前記モニタマイクロフォン(3c)からの低減音信号と検出
マイクロフォン(3a)からの騒音信号との検出サンプリン
グ回数をカウントし、このサンプリング回数が所定値に
達したときにフィルタ更新信号を出力するサンプリング
回数カウント手段(5) と、該サンプリング回数カウント
手段(5) からのフィルタ更新信号が受信可能とされ、前
記モニタマイクロフォン(3c)からの低減音信号と検出マ
イクロフォン(3a)からの波形信号との検出サンプリング
毎に前記適応型FIRフィルタ(41)のフィルタ係数を更
新するための係数ベクトルを計算して加算し、前記サン
プリング回数カウント手段(5) からのフィルタ更新信号
を受けたときに、この加算された係数ベクトルに基いて
前記適応型FIRフィルタ(41)のフィルタ係数を更新さ
せる係数ベクトル加算手段(6) とを備えさせるような構
成としている。
A fourth aspect of the present invention uses the device according to the second aspect as a silencer. Specifically, as shown in FIG. 2, the noise in the space (21) is detected. A detection microphone (3a) that outputs a noise signal, and an adaptive FIR filter (41) that generates an inverted sound signal of the same amplitude in antiphase with respect to the noise signal from the detection microphone (3a),
A speaker (3b) that receives a reversal sound signal output from the adaptive FIR filter (41) and emits a reversal sound based on the reversal sound signal to the space, and is arranged at a predetermined observation point in the space (21). , A monitor microphone (3c) that detects the reduced sound level at the observation point and outputs a reduced sound signal as feedback, and a noise signal from the detection microphone (3a) and a reduced sound signal fed back from the monitor microphone (3c). Active control including an adaptive algorithm execution circuit (44) for inputting and updating the filter coefficient of the adaptive FIR filter (41) so as to reduce the sound pressure level around the observation point by a noise signal and a reduced sound signal. It is based on the device. Then, in the adaptive algorithm execution circuit (44),
Counts the number of detection samplings of the reduced sound signal from the monitor microphone (3c) and the noise signal from the detection microphone (3a), and outputs a filter update signal when the number of samplings reaches a predetermined value. Means (5) and the filter update signal from the sampling number counting means (5) can be received, and the detection sampling of the reduced sound signal from the monitor microphone (3c) and the waveform signal from the detection microphone (3a) A coefficient vector for updating the filter coefficient of the adaptive FIR filter (41) is calculated and added for each time, and when the filter update signal from the sampling number counting means (5) is received, this addition is made. Coefficient vector adding means (6) for updating the filter coefficient of the adaptive FIR filter (41) based on the coefficient vector, It is configured to be equipped with.

【0010】[0010]

【作用】以上の構成により、本発明では以下に述べるよ
うな作用が得られる。請求項1または2記載の発明で
は、波形検出手段(3a)により変動波形信号が検出され
て、該波形検出手段(3a)から出力される波形信号に対し
て、該波形信号に基いて生成された制御信号がデジタル
フィルタ(41)から出力される。そして、この制御信号に
対応した作動信号を信号出力手段(3b)が出力することに
より、変動波形信号に作動信号の影響を与える。その
後、信号出力手段(3b)が出力した作動信号と前記変動波
形信号との誤差を誤差検出手段(3c)により検出して、該
誤差検出手段(3c)から出力される誤差信号と前記波形検
出手段(3a)から出力される波形信号とを適応アルゴリズ
ム実行手段(44)が受ける。そして、該適応アルゴリズム
実行手段(44)により前記誤差が小さくなるように前記デ
ジタルフィルタ(41)のフィルタ係数が更新される。そし
て、このデジタルフィルタ(41)のフィルタ係数の更新に
際し、サンプリング回数カウント手段(5) により、前記
誤差検出手段(3c)からの誤差信号と波形検出手段(3a)か
らの波形信号との検出サンプリング回数がカウントさ
れ、このサンプリング回数が所定値に達したときにサン
プリング回数カウント手段(5) からフィルタ更新信号が
係数ベクトル加算手段(6) に出力される。一方、係数ベ
クトル加算手段(6) にあっては、前記誤差検出手段(3c)
からの誤差信号と波形検出手段(3a)からの波形信号との
検出サンプリング毎に前記デジタルフィルタ(41)のフィ
ルタ係数を更新するための係数ベクトルが計算されて加
算されており、前記サンプリング回数カウント手段(5)
からのフィルタ更新信号が係数ベクトル加算手段(6) に
出力されたときに、この加算された係数ベクトルに基い
て前記デジタルフィルタ(41)のフィルタ係数が更新され
ることになる。そして、この際、係数ベクトル加算手段
(6) で加算されている係数ベクトルは収束誤差が相殺さ
れており、更新されたフィルタ係数は収束性の高い値と
なっている。また、この際の計算量は従来のLMSアル
ゴリズムと同等である。
With the above construction, the following effects can be obtained in the present invention. In the invention according to claim 1 or 2, the fluctuation waveform signal is detected by the waveform detection means (3a), and is generated based on the waveform signal with respect to the waveform signal output from the waveform detection means (3a). The control signal is output from the digital filter (41). The signal output means (3b) outputs an actuation signal corresponding to this control signal, so that the fluctuation waveform signal is affected by the actuation signal. Thereafter, the error between the actuation signal output by the signal output means (3b) and the fluctuation waveform signal is detected by the error detection means (3c), and the error signal output from the error detection means (3c) and the waveform detection The adaptive algorithm executing means (44) receives the waveform signal output from the means (3a). Then, the adaptive algorithm executing means (44) updates the filter coefficient of the digital filter (41) so as to reduce the error. When updating the filter coefficient of the digital filter (41), the sampling number counting means (5) detects the sampling of the error signal from the error detecting means (3c) and the waveform signal from the waveform detecting means (3a). The number of times is counted, and when the number of times of sampling reaches a predetermined value, the sampling number counting means (5) outputs a filter update signal to the coefficient vector addition means (6). On the other hand, in the coefficient vector addition means (6), the error detection means (3c)
The coefficient vector for updating the filter coefficient of the digital filter (41) is calculated and added for each detection sampling of the error signal from the waveform signal from the waveform detection means (3a), and the sampling frequency is counted. Means (5)
When the filter update signal from is output to the coefficient vector adding means (6), the filter coefficient of the digital filter (41) is updated based on the added coefficient vector. At this time, the coefficient vector adding means
The convergence error is canceled out in the coefficient vector added in (6), and the updated filter coefficient has a high convergence value. Further, the calculation amount at this time is equivalent to that of the conventional LMS algorithm.

【0011】請求項3記載の発明では、サンプリング回
数カウント手段(5) は、サンプリング回数がデジタルフ
ィルタ(41)のフィルタ係数の数に一致したときにフィル
タ更新信号を係数ベクトル加算手段(6) に出力し、これ
によってデジタルフィルタ(41)のフィルタ係数が更新さ
れる。
In the invention according to claim 3, the sampling number counting means (5) sends the filter update signal to the coefficient vector addition means (6) when the sampling number matches the number of filter coefficients of the digital filter (41). Then, the filter coefficient of the digital filter (41) is updated.

【0012】請求項4記載の発明では、適応型FIRフ
ィルタ(41)から出力された反転音信号がスピーカ(3b)に
入力され、該スピーカ(3b)から空間(21)に反転音が放射
される。また、モニタマイクロホン(3c)が観測点におい
て低減音レベルを検出して適応アルゴリズム実行回路(4
4)に低減音信号をフィードバックする。そして、この適
応アルゴリズム実行回路(44)により、観測点での音圧レ
ベルが低減されるように適応型FIRフィルタ(41)のフ
ィルタ係数が逐次更新される。そして、この適応型FI
Rフィルタ(41)のフィルタ係数の更新に際しては、上述
した請求項1または2記載の発明に係る作用と同様にし
て係数ベクトル加算手段(6) での係数ベクトルの加算に
よりその収束誤差が相殺されていることになり、更新さ
れたフィルタ係数は収束性の高い値となって、消音性能
が向上される。また、この消音動作における計算量は従
来のLMSアルゴリズムと同等となっている。
In the invention according to claim 4, the inverted sound signal output from the adaptive FIR filter (41) is input to the speaker (3b), and the inverted sound is radiated from the speaker (3b) to the space (21). It Also, the monitor microphone (3c) detects the reduced sound level at the observation point, and the adaptive algorithm execution circuit (4
The reduced sound signal is fed back to 4). Then, the adaptive algorithm execution circuit (44) sequentially updates the filter coefficient of the adaptive FIR filter (41) so that the sound pressure level at the observation point is reduced. And this adaptive FI
When the filter coefficient of the R filter (41) is updated, the convergence error is canceled by the addition of the coefficient vector by the coefficient vector addition means (6) in the same manner as the operation according to the invention described in claim 1 or 2. Therefore, the updated filter coefficient becomes a value with high convergence and the silencing performance is improved. Further, the amount of calculation in this muffling operation is equivalent to that of the conventional LMS algorithm.

【0013】[0013]

【実施例】以下、本発明の実施例を図面に基いて説明す
る。本例では、ダクト内の消音を行うアクティブ消音装
置に、本発明を採用した場合について説明する。図2
は、本例に係るアクティブ消音装置(1) の全体構成を示
し、ダクト(2) の内部には、図2の左端に位置する騒音
源(22)からの騒音音波の伝播通路としての空間(21)が形
成されている。また、前記ダクト(2) には、図2に矢印
で示す騒音の伝播方向の上流側から下流側に向って検出
マイクロフォン(3a)とスピーカ(3b)とモニタマイクロフ
ォン(3c)とが順に配置されている。そして、検出マイク
ロフォン(3a)は、前記空間(21)における変動波形信号で
ある騒音を検出し、波形信号である騒音信号をコントロ
ーラ(4) に出力する波形検出手段を構成している。ま
た、前記スピーカ(3b)は、前記コントローラ(4) からの
制御信号である反転音信号を受け、騒音とは逆位相で且
つ同振幅の反転音、つまり、作動信号を前記ダクト(2)
内の空間(21)に放射するための付加音源であって、信号
出力手段を構成している。また、前記モニタマイクロフ
ォン(3c)は、所定の観測点に位置し、前記スピーカ(3b)
から放射された反転音の作用により低減された低減音レ
ベルをその観測点で検出し、その低減音信号を出力して
いる。つまり、スピーカ(3b)からの作動信号である反転
音と前記変動波形信号である騒音との誤差である低減音
レベルを検出して誤差信号としての低減音信号を出力す
る誤差検出手段を構成している。
Embodiments of the present invention will be described below with reference to the drawings. In this example, a case will be described where the present invention is applied to an active muffling device that muffles the inside of a duct. Figure 2
Shows the overall structure of the active noise suppressor (1) according to the present example. Inside the duct (2), there is a space (as a propagation path for noise sound waves from the noise source (22) located at the left end of FIG. 21) has been formed. Further, in the duct (2), a detection microphone (3a), a speaker (3b), and a monitor microphone (3c) are sequentially arranged from the upstream side to the downstream side in the noise propagation direction shown in FIG. ing. The detection microphone (3a) constitutes a waveform detection means for detecting noise that is a fluctuating waveform signal in the space (21) and outputting the noise signal that is a waveform signal to the controller (4). Further, the speaker (3b) receives a reversal sound signal which is a control signal from the controller (4), and outputs a reversal sound having a phase opposite to that of noise and having the same amplitude, that is, an operation signal to the duct (2).
It is an additional sound source for radiating to the internal space (21) and constitutes a signal output means. Further, the monitor microphone (3c) is located at a predetermined observation point, and the speaker (3b)
The reduced sound level reduced by the action of the inversion sound radiated from is detected at the observation point, and the reduced sound signal is output. That is, an error detecting means is configured to detect a reduced sound level which is an error between the inversion sound which is the operation signal from the speaker (3b) and the noise which is the fluctuation waveform signal, and outputs the reduced sound signal as the error signal. ing.

【0014】前記コントローラ(4) は、騒音とは逆位相
で且つ同振幅の反転音信号を生成するためのもので、前
記スピーカ(3b)に対して前記反転音信号がD/A変換器
(図示せず)等を介して出力するようになっている。ま
た、該コントローラ(4) には、前記検出マイクロフォン
(3a)の騒音信号とモニタマイクロフォン(3c)の低減音信
号がA/D変換器(図示せず)等を介して入力されてい
る。そして、該コントローラ(4) は、適応型FIRフィ
ルタ(41)と、第1フィルタ(42)及び第2フィルタ(43)と
を備えると共に、適応アルゴリズム実行回路(44)を備え
ている。
The controller (4) is for generating an inverted sound signal having a phase opposite to that of the noise and having the same amplitude, and the inverted sound signal for the speaker (3b) is a D / A converter ( (Not shown) or the like. Further, the controller (4) includes the detection microphone.
The noise signal (3a) and the reduced sound signal of the monitor microphone (3c) are input via an A / D converter (not shown) or the like. The controller (4) includes an adaptive FIR filter (41), a first filter (42) and a second filter (43), and an adaptive algorithm execution circuit (44).

【0015】該適応型FIRフィルタ(41)は、前記検出
マイクロフォン(3a)が出力する騒音信号を加算回路(46)
を介して受け、該騒音信号とは基本的に逆位相で同振幅
の反転音信号を生成するデジタルフィルタを構成してい
る。また、前記第1フィルタ(42)は、検出マイクロフォ
ン(3a)が出力する騒音信号を加算回路(46)を介して受
け、前記適応型FIRフィルタ(41)の反転音信号に基い
てスピーカ(3b)から放射された反転音がモニタマイクロ
フォン(3c)に入力されるのに要する伝播時間だけ検出マ
イクロフォン(3a)の騒音信号を予め所定時間遅延させる
と共に、波形減衰等を考慮したフィルタであって、スピ
ーカ(3b)とモニタマイクロフォン(3c)との間の伝達関数
を有するFIRフィルタで構成されている。また、前記
第2フィルタ(43)は、スピーカ(3b)から放射された反転
音が検出マイクロフォン(3a)に伝播してハウリングが発
生するのを防ぐために前記適応型FIRフィルタ(41)か
らの反転音信号を加算回路(46)にフィードバックするフ
ィルタであって、所定の伝達関数を有するFIRフィル
タで構成されている。更に、前記適応アルゴリズム実行
回路(44)は、最小二乗平均法(LMS;Least Mean
Square)アルゴリズムによる適応制御を行う適応アルゴ
リズム実行手段であって、前記第1フィルタ(42)を通し
て受ける騒音信号を遅延した信号と、モニタマイクロフ
ォン(3c)からフィードバックされる低減音信号とに基い
て、ダクト(2) 内の観測点周辺の音圧レベルを低減する
ようにLMSの制御パラメータとしての適応型FIRフ
ィルタ(41)のフィルタ係数を更新してその反転音信号を
適応制御し、補正するように構成されている。
The adaptive FIR filter (41) adds a noise signal output from the detection microphone (3a) to an adder circuit (46).
And a digital filter that generates a reversal sound signal having the same amplitude as that of the noise signal but having a phase opposite to that of the noise signal. The first filter (42) receives the noise signal output from the detection microphone (3a) via the adder circuit (46), and the speaker (3b) based on the inverted sound signal of the adaptive FIR filter (41). In addition to delaying the noise signal of the detection microphone (3a) by a predetermined time by the propagation time required for the inverted sound radiated from () to be input to the monitor microphone (3c), a filter considering waveform attenuation and the like, It is composed of an FIR filter having a transfer function between the speaker (3b) and the monitor microphone (3c). Further, the second filter (43) is an inversion from the adaptive FIR filter (41) in order to prevent the inverted sound radiated from the speaker (3b) from propagating to the detection microphone (3a) and causing howling. The filter is a filter for feeding back the sound signal to the adder circuit (46) and is composed of an FIR filter having a predetermined transfer function. Further, the adaptive algorithm execution circuit (44) uses the least mean squares method (LMS; Least Mean).
Square) adaptive algorithm execution means for performing adaptive control by an algorithm, which is based on a signal obtained by delaying a noise signal received through the first filter (42) and a reduced sound signal fed back from a monitor microphone (3c), Update the filter coefficient of the adaptive FIR filter (41) as the LMS control parameter so as to reduce the sound pressure level around the observation point in the duct (2), and adaptively control and correct the inverted sound signal. Is configured.

【0016】そして、本例の特徴としては、前記適応型
FIRフィルタ(41)のフィルタ係数を更新する適応アル
ゴリズム実行回路(44)における信号処理の手順にある。
具体的には、本適応アルゴリズム実行回路(44)では、各
マイクロフォン(3a),(3c) での騒音検出のサンプリング
毎に全てのフィルタ係数の更新を行うのではなく、複数
回のサンプリングが行われ、この各サンプリング夫々で
計算されたフィルタ係数ベクトルを合算した後、この合
算して得られたフィルタ係数ベクトルに基いて全てのフ
ィルタ係数を更新するようにしている。つまり、複数回
サンプリング毎に1回のフィルタ係数の更新が実行され
るようになっている。以下、この適応アルゴリズム実行
回路(44)における信号処理動作を図3のフローチャート
に沿って説明する。消音制御がスタートした後、先ず、
ステップST0においてフィルタ係数ベクトル(W) を初期
値0に設定し、その後、ステップST1において、フィル
タ係数の更新が行われるまでのサンプリング数をカウン
トするためのループ変数(LOOP)を初期値0に設定すると
共に、フィルタ係数の更新が行われるまでの各サンプリ
ング夫々で計算されたフィルタ係数ベクトルが加算され
てなるフィルタ係数ベクトルバッファ(WBUF)を初期値0
に設定する。その後、ステップST2において1回(今
回)のサンプリングにおける検出信号に応じたフィルタ
出力(y(n))を算出する。このフィルタ出力(y(n))の算出
はフィルタ係数ベクトル(W(n))と入力ベクトル(X(n-k))
との積の畳込み演算によって行われる。また、ステップ
ST3においてはステップST2において求められたフィル
タ出力(y(n))における誤差(e(n))が算出される。この誤
差(e(n))の算出は制御目標値(d(n))つまりアクティブ消
音にあっては消音しようとする騒音の出力から前記フィ
ルタ出力(y(n))を減算することによって行われる。そし
て、その後、ステップST4において、前記フィルタ係数
ベクトルバッファ(WBUF)に今回のサンプリングにおける
フィルタ係数ベクトルが加算される。具体的には、前回
までのフィルタ係数ベクトルバッファ(WBUF(n-1)) に、
誤差(e(n))の2倍とステップサイズパラメータ(myu) と
入力ベクトル(X(n))との積が加算されて今回のサンプリ
ングまでのフィルタ係数ベクトルバッファ(WBUF)が求め
られることになる。その後、ステップST5において前記
ループ変数(LOOP)に1が加算され、ステップST6におい
てループ変数(LOOP)が所定の値(NUMBER)に達したか否か
が判定される。尚、この所定の値(NUMBER)は、フィルタ
係数の更新を行うまでのサンプリング数を設定したもの
であって、例えばフィルタ係数の数と一致する値に設定
される。そして、このステップST6においてループ変数
(LOOP)が所定の値(NUMBER)に未だ達していない場合に
は、再び前記ステップST2に戻り、次のサンプリングに
おけるフィルタ出力(y(n))、誤差(e(n))及びフィルタ係
数ベクトルバッファ(WBUF)の算出が行われることにな
る。このような動作が複数回行われて、ループ変数(LOO
P)が所定の値(NUMBER)に達するとステップST7 に移り、
フィルタ係数ベクトル(W) として、それまでのフィルタ
係数ベクトル(W) にフィルタ係数ベクトルバッファ(WBU
F)が加算されて、それに基いて前記適応型FIRフィル
タ(41)の全てのフィルタ係数が更新されるようになって
いる。また、このように適応型FIRフィルタ(41)のフ
ィルタ係数が更新された後は、再びステップST1に戻っ
て上述した動作が繰り返し行われるようになっている。
A feature of this example is the signal processing procedure in the adaptive algorithm execution circuit (44) for updating the filter coefficient of the adaptive FIR filter (41).
Specifically, the adaptive algorithm execution circuit (44) does not update all the filter coefficients for each noise detection sampling in each microphone (3a), (3c), but performs multiple samplings. That is, after summing up the filter coefficient vectors calculated in each of the samplings, all the filter coefficients are updated based on the filter coefficient vector obtained by the summing. That is, the filter coefficient is updated once for each sampling. The signal processing operation of the adaptive algorithm execution circuit (44) will be described below with reference to the flowchart of FIG. After the silence control starts, first,
In step ST0, the filter coefficient vector (W) is set to an initial value 0, and then in step ST1, a loop variable (LOOP) for counting the number of samplings until the filter coefficient is updated is set to an initial value 0. At the same time, the filter coefficient vector buffer (WBUF) formed by adding the filter coefficient vectors calculated in each sampling until the filter coefficient is updated is initialized to 0.
Set to. After that, in step ST2, the filter output (y (n)) according to the detection signal in one (current) sampling is calculated. This filter output (y (n)) is calculated by the filter coefficient vector (W (n)) and the input vector (X (nk)).
It is performed by the convolution operation of the product of and. Also step
In ST3, the error (e (n)) in the filter output (y (n)) obtained in step ST2 is calculated. This error (e (n)) is calculated by subtracting the filter output (y (n)) from the control target value (d (n)), that is, the output of the noise to be silenced in active muffling. Be seen. Then, after that, in step ST4, the filter coefficient vector in the current sampling is added to the filter coefficient vector buffer (WBUF). Specifically, in the filter coefficient vector buffer (WBUF (n-1)) up to the last time,
Double the error (e (n)), the product of the step size parameter (myu) and the input vector (X (n)) are added to obtain the filter coefficient vector buffer (WBUF) up to this sampling. Become. After that, 1 is added to the loop variable (LOOP) in step ST5, and it is determined in step ST6 whether or not the loop variable (LOOP) has reached a predetermined value (NUMBER). The predetermined value (NUMBER) sets the number of samplings until the filter coefficient is updated, and is set to a value that matches the number of filter coefficients, for example. Then, in this step ST6, the loop variable
If (LOOP) has not reached the predetermined value (NUMBER) yet, the process returns to step ST2 again, and the filter output (y (n)), error (e (n)) and filter coefficient vector in the next sampling The buffer (WBUF) will be calculated. This kind of operation is performed multiple times, and the loop variable (LOO
When (P) reaches a predetermined value (NUMBER), move to step ST7,
As the filter coefficient vector (W), the filter coefficient vector buffer (WBU) is added to the previous filter coefficient vector (W).
F) is added, and all the filter coefficients of the adaptive FIR filter (41) are updated based on the addition. After the filter coefficient of the adaptive FIR filter (41) is updated in this way, the process returns to step ST1 and the above-described operation is repeated.

【0017】次に、上記のアクティブ消音装置を用いて
ダクト(2) 内の騒音を消音する消音動作について説明す
る。先ず、ダクト(2) 内の空間(21)における騒音は、検
出マイクロフォン(3a)により検出され、この検出マイク
ロフォン(3a)によって検出された騒音信号は、加算回路
(46)を経て適応型FIRフィルタ(41)に入力される。ま
た、この騒音信号は、第1フィルタ(42)で所定時間だけ
遅延処理等が行われて適応アルゴリズム実行回路(44)に
も入力されている。そして、前記適応型FIRフィルタ
(41)においては、前記騒音信号とは逆位相で且つ同振幅
の反転音信号が生成され、この反転音信号がスピーカ(3
b)に出力されて、該スピーカ(3b)から反転音が空間(21)
に放射される。この放射により、ダクト(2) 内の騒音
は、スピーカ(3b)から放射された反転音により良好に低
減されることになる。更に、観測点では低減された騒音
レベル、つまり低減音レベルがモニタマイクロフォン(3
c)により検出され、この検出信号である低減音信号が適
応アルゴリズム実行回路(44)にフィードバックされて、
この低減音信号と第1フィルタ(42)からの騒音信号とに
基いて適応型FIRフィルタ(41)のフィルタ係数が更新
される。このフィルタ係数の更新により、前記適応型F
IRフィルタ(41)による反転音信号の生成がダクト(2)
内の騒音に対して経時的に精度良く行われ、ダクト(2)
内の観測点周辺の音圧レベルが最も良好に低減される。
そして、この適応アルゴリズム実行回路(44)による適応
型FIRフィルタ(41)のフィルタ係数の更新に際して
は、上述したように、複数回のサンプリングが行われ、
この各サンプリング夫々で計算されたフィルタ係数ベク
トルを合算した後、この合算して得られたフィルタ係数
ベクトルに基いて全てのフィルタ係数が更新されること
になる。これにより、各サンプリング夫々で計算された
フィルタ係数ベクトル、つまりフィルタ出力は収束誤差
が含まれているものであるが、これらを複数回のサンプ
リングに亘って合算することにより、各収束誤差が相殺
されることになって、フィルタ係数の更新時には極めて
収束性の高いフィルタ係数が得られることになる。つま
り、複数回サンプリングに1回のフィルタ係数の更新を
実行することにより、適応制御の安定且つ適切な収束状
態を得ることができることになる。また、この制御にお
ける計算量は従来のLMSアルゴリズムと同等であって
実用性も高いものである。つまり、従来のLMSアルゴ
リズムでは1回サンプリング毎にフィルタ係数の更新を
実行するのに対し、本例では、複数回のサンプリングに
1回のフィルタ係数の更新を実行するようにしているこ
とから計算量としては同等である。
Next, the silencing operation for silencing the noise in the duct (2) using the above active silencing device will be described. First, the noise in the space (21) in the duct (2) is detected by the detection microphone (3a), and the noise signal detected by this detection microphone (3a) is added by an addition circuit.
It is input to the adaptive FIR filter (41) via (46). The noise signal is also subjected to delay processing and the like by the first filter (42) for a predetermined time and is also input to the adaptive algorithm execution circuit (44). And the adaptive FIR filter
In (41), an inverted sound signal having the opposite phase and the same amplitude as the noise signal is generated, and the inverted sound signal is generated by the speaker (3
b) and the inverted sound is output from the speaker (3b) in the space (21).
Is emitted to. Due to this radiation, the noise in the duct (2) is favorably reduced by the reversal sound radiated from the speaker (3b). Furthermore, the noise level reduced at the observation point, that is, the reduced sound level
c), the reduced sound signal that is this detection signal is fed back to the adaptive algorithm execution circuit (44),
The filter coefficient of the adaptive FIR filter (41) is updated based on this reduced sound signal and the noise signal from the first filter (42). By updating this filter coefficient, the adaptive F
The generation of the inverted sound signal by the IR filter (41) is the duct (2)
The duct (2) is performed with high accuracy with time against the noise inside.
The sound pressure level around the observation point inside is best reduced.
Then, when the filter coefficient of the adaptive FIR filter (41) is updated by the adaptive algorithm execution circuit (44), sampling is performed a plurality of times as described above,
After the filter coefficient vectors calculated by the respective samplings are added up, all the filter coefficients are updated based on the filter coefficient vector obtained by this addition. As a result, the filter coefficient vector calculated in each sampling, that is, the filter output contains the convergence error, but by adding these over multiple samplings, the convergence errors are canceled out. As a result, when updating the filter coefficient, a filter coefficient having extremely high convergence is obtained. That is, it is possible to obtain a stable and appropriate convergence state of the adaptive control by executing the update of the filter coefficient once every plural times of sampling. Further, the amount of calculation in this control is equivalent to that of the conventional LMS algorithm and is highly practical. That is, in the conventional LMS algorithm, the filter coefficient is updated for each sampling, whereas in this example, the filter coefficient is updated once for a plurality of samplings. Are equivalent.

【0018】次に、上記の効果を確認するために行った
解析及び実験について説明する。今、図4に示すような
消音システム構造を考え、未知システム(2) からの信号
に対して適応型FIRフィルタ(41)のフィルタ係数を更
新することによって、前記未知システム(2) からの信号
レベルを低減させるようなものにおいて、未知システム
(2) の音響的伝達関数を設定するシステム同定を行う際
について解析する。具体的には、適応型FIRフィルタ
(41)のフィルタ係数の数を128 個とし、未知システム
(2) のインパルス応答として図5に示すようなAタイプ
の未知システム(ERROR PATH)と図6に示すようなBタイ
プの未知システム(SIMPLE DELAY)とにおいて夫々システ
ム同定を行い、その解析結果を図7及び図8に示す。
尚、図7に示すものはAタイプの未知システムに対して
の解析結果であって、図8に示すものはBタイプの未知
システムに対しての解析結果であり、横軸がフィルタ係
数更新の繰返し回数であり、縦軸が騒音レベルである。
また、図7及び図8において、2点鎖線(各図のA)で
示したものは本装置を作動させていない状態でのシステ
ム同定音のレベル、破線(各図のB)で示したものは従
来のLMSアルゴリズム(1サンプリング毎に全フィル
タ係数を更新するようにしたもの)におけるシステム同
定を行った場合の騒音レベル、実線(各図のC)で示し
たものは本発明に係るLMSアルゴリズムにおけるシス
テム同定を行った場合の騒音レベルであって、1回のフ
ィルタ係数の更新が実行されるサンプリング数を適応フ
ィルタ(41)のフィルタ係数の数と同一の値の128 に設定
した場合、1点鎖線(各図のD)で示したものは本発明
に係るLMSアルゴリズムにおけるシステム同定を行っ
た場合の騒音レベルであって、1回のフィルタ係数の更
新が実行されるサンプリング数を適応フィルタ(41)のフ
ィルタ係数の数の2倍の値の256 に設定した場合を夫々
示している。これら各図から判るように、本発明に係る
LMSアルゴリズムを使用した場合には、従来のLMS
アルゴリズムを使用した場合よりも騒音レベルが低く、
収束性に優れていることが判る。また、特に、図8に示
すものおいて、本図には示していないが、1回のフィル
タ係数の更新が実行されるサンプリング数を1から次第
に大きくしていくと収束性が向上していき、図8に実線
Cで示す128 に設定した場合に最も騒音レベルが低減さ
れ、この128 から更にサンプリング数を多くしていくと
騒音レベルが上昇することが判った。つまり、このBタ
イプの未知システムにあっては1回のフィルタ係数の更
新が実行されるサンプリング数を適応フィルタ(41)のフ
ィルタ係数の数と同一の値に設定することにより、本発
明の効果が最も大きく発揮されることが判る。一方、図
7に示すものでは、サンプリング数が128 の場合と256
の場合とでは騒音レベルが略同一であるが、サンプリン
グ数を、これら以外の値に設定しても騒音レベルが更に
低減するようなことはなかった。つまり、このAタイプ
の未知システムにあっては1回のフィルタ係数の更新が
実行されるサンプリング数を適応フィルタ(41)のフィル
タ係数の数と同一の値からフィルタ係数の数の2倍程度
までの範囲に設定することにより、本発明の効果が最も
大きく発揮されることが判る。
Next, the analysis and experiments conducted to confirm the above effects will be described. Now, considering the silencing system structure as shown in FIG. 4, the signal from the unknown system (2) is updated by updating the filter coefficient of the adaptive FIR filter (41) with respect to the signal from the unknown system (2). Unknown systems, such as those that reduce levels
The system identification for setting the acoustic transfer function of (2) is analyzed. Specifically, an adaptive FIR filter
The number of filter coefficients in (41) is 128 and unknown system
As the impulse response of (2), system identification is performed for each of the A type unknown system (ERROR PATH) as shown in FIG. 5 and the B type unknown system (SIMPLE DELAY) as shown in FIG. 6, and the analysis results are shown. This is shown in FIGS. 7 and 8.
Incidentally, what is shown in FIG. 7 is the analysis result for the A type unknown system, and what is shown in FIG. 8 is the analysis result for the B type unknown system, where the horizontal axis indicates the update of the filter coefficient. It is the number of repetitions, and the vertical axis is the noise level.
Further, in FIGS. 7 and 8, what is indicated by a two-dot chain line (A in each figure) is the level of the system identification sound in a state where this device is not operated, and what is indicated by a broken line (B in each figure). Is the noise level when system identification is performed using the conventional LMS algorithm (all filter coefficients are updated every sampling), and the solid line (C in each figure) indicates the LMS algorithm according to the present invention. Noise level when system identification is performed, and when the number of samplings at which the filter coefficient is updated once is set to 128, which is the same value as the number of filter coefficients of the adaptive filter (41), What is indicated by a dotted chain line (D in each figure) is the noise level when the system identification is performed in the LMS algorithm according to the present invention, and the noise level is updated once. The case set to 256 twice the value of the number of filter coefficients of the adaptive filter (41) the number of ring represents respectively. As can be seen from these figures, when the LMS algorithm according to the present invention is used, the conventional LMS algorithm is used.
The noise level is lower than with the algorithm,
It can be seen that the convergence is excellent. Further, in particular, although not shown in this figure in FIG. 8, convergence is improved by gradually increasing the number of samplings at which one update of the filter coefficient is executed from one. It was found that the noise level was reduced most when it was set to 128 shown by the solid line C in FIG. 8, and the noise level increased as the sampling number was further increased from 128. That is, in this B type unknown system, the effect of the present invention can be obtained by setting the number of samplings for which the filter coefficient is updated once to the same value as the number of filter coefficients of the adaptive filter (41). It can be seen that is maximized. On the other hand, in the case shown in FIG. 7, the number of samplings is 128 and 256.
Although the noise level was almost the same as that in the above case, the noise level was not further reduced even if the sampling number was set to a value other than these. That is, in this type A unknown system, the number of samplings at which the filter coefficient is updated once is from the same value as the number of filter coefficients of the adaptive filter (41) to about twice the number of filter coefficients. It can be seen that the effect of the present invention is maximized by setting the range to.

【0019】更に、制御の安定性を確認するための解析
として、図9に示すように、104繰返し回数毎に上記
の各未知システムを交互に変更したような場合について
の騒音レベルを求めた場合でも、未知システムが変更さ
れた直後から騒音レベルは低減されて収束方向に向うこ
とが確認された。つまり、システムの状態が変化したよ
うな場合であっても発散するようなことなしにシステム
に追従した消音動作を行うことができることが確認され
たことになる。
Further, as an analysis for confirming the control stability, as shown in FIG. 9, the noise level was obtained in the case where the above unknown systems were alternately changed every 10 4 repetitions. Even in this case, it was confirmed that immediately after the unknown system was changed, the noise level was reduced and headed toward the convergence direction. In other words, it has been confirmed that even if the state of the system changes, the silencing operation that follows the system can be performed without diverging.

【0020】また、上述したような1回のフィルタ係数
の更新が実行されるサンプリング数と収束時の誤差との
関係を解析した結果を図10に示す。この図10におけ
る横軸は1回のフィルタ係数の更新が実行されるサンプ
リング数(タップ数)であり、縦軸は収束時における誤
差、つまり理想出力に対する該理想出力とフィルタ出力
との差である。また、この図10におけるEは図5に示
すAタイプの未知システムを対象としたもので、Fは図
6に示すBタイプの未知システムを対象としたものであ
る。この図から判るように、1回のフィルタ係数の更新
が実行されるサンプリング数を増大させていくとある値
までは収束時の誤差が低減されることになる。そして、
Aタイプの未知システムにあっては適応フィルタ(41)の
フィルタ係数の数と同一の値(128) に設定した場合に最
も誤差が小さくなり、256 以上に設定した場合には発散
してしまう。一方、Bタイプの未知システムにあっては
適応フィルタ(41)のフィルタ係数の数の2倍の値(256)
に設定した場合に最も誤差が小さくなり、それ以上に設
定した場合には発散してしまう。このことから、本例の
LMSアルゴリズムでは、1回のフィルタ係数の更新が
実行されるサンプリング数を適応フィルタ(41)のフィル
タ係数の数と同一の値に設定するようにすれば、未知シ
ステムに拘りなく且つ発散を招くこともなく高い収束性
を確保することができることが判る。
FIG. 10 shows the result of analysis of the relationship between the number of samplings at which the filter coefficient is updated once as described above and the error at the time of convergence. The horizontal axis in FIG. 10 represents the number of samplings (the number of taps) at which the filter coefficient is updated once, and the vertical axis represents the error at the time of convergence, that is, the difference between the ideal output and the filter output. . Further, E in FIG. 10 is intended for the A type unknown system shown in FIG. 5, and F is intended for the B type unknown system shown in FIG. As can be seen from this figure, if the number of samplings in which the filter coefficient is updated once is increased, the error at the time of convergence is reduced to a certain value. And
In the unknown system of type A, the error becomes the smallest when the value is set to the same value (128) as the number of filter coefficients of the adaptive filter (41), and the error occurs when the value is set to 256 or more. On the other hand, in an unknown system of type B, the value (256) that is twice the number of filter coefficients of the adaptive filter (41)
The error is the smallest when set to, and diverges when set to more. From this, in the LMS algorithm of the present example, if the number of samplings at which the filter coefficient is updated once is set to the same value as the number of filter coefficients of the adaptive filter (41), an unknown system can be obtained. It can be seen that high convergence can be secured without concern and without causing divergence.

【0021】次に、従来の構成のものと本発明に係る消
音装置との消音性能について行った実験結果を説明す
る。つまり、アクティブ消音を行わないもの、従来のL
MSアルゴリズムを使用したもの、本発明のLMSアル
ゴリズムを使用したものとの夫々について各周波数帯で
の音圧レベルを測定し、その結果を図11に示す。本図
の棒グラフにおいて、最も左側に位置する白抜きはアク
ティブ消音を行わないもの、中央の右斜下ハッチングの
ものは本発明のLMSアルゴリズムを使用したもの、最
も右側に位置する左斜下ハッチングのものは従来のLM
Sアルゴリズムを使用したものを夫々示している。この
ように、広い帯域おいて本発明のLMSアルゴリズムを
使用したものにおける消音効果が確認された。尚、本図
11におけるLは各周波数帯での音圧レベルを線形加算
したものを示し、AはこれにA補正(人の聴感に対応さ
せて各周波数帯での加算割合を補正)したものを示して
いる。
Next, the result of an experiment conducted on the silencing performance between the conventional construction and the silencing apparatus according to the present invention will be described. In other words, the one that does not perform active muting, the conventional L
The sound pressure level in each frequency band was measured for each of the one using the MS algorithm and the one using the LMS algorithm of the present invention, and the results are shown in FIG. In the bar graph of this figure, the white areas on the far left are those for which active muffling is not performed, those in the center for the right slanted bottom hatching are those using the LMS algorithm of the present invention, and those for the leftmost shaded left slanted bottom hatch are Things are conventional LM
Those using the S algorithm are shown respectively. As described above, the sound deadening effect was confirmed in the wide band using the LMS algorithm of the present invention. Note that L in FIG. 11 indicates that the sound pressure level in each frequency band is linearly added, and A is corrected by A (the addition ratio in each frequency band is corrected according to human hearing). Is shown.

【0022】また、上述した効果に加えて、本発明のよ
うなLMSアルゴリズムを採用したものでは、例えば適
応型FIRフィルタ(41)及び適応アルゴリズム実行回路
(44)を固定小数点演算型のプロセッサを使用したような
場合に、桁落ちによる打ち切り誤差が相殺されることに
なって、この打ち切り誤差の蓄積が抑制されることにな
るので、この種のプロセッサを用いた場合での消音性能
の向上を図ることができることになる。
In addition to the above-mentioned effects, in the case of adopting the LMS algorithm like the present invention, for example, an adaptive FIR filter (41) and an adaptive algorithm execution circuit are provided.
If (44) is used with a fixed-point arithmetic type processor, the truncation error due to cancellation will be canceled out, and the accumulation of this truncation error will be suppressed. It is possible to improve the silencing performance in the case of using.

【0023】尚、本実施例では、本発明をダクト内の消
音を行うアクティブ消音装置に採用した場合について説
明したが、請求項1乃至3記載の発明はこれに限らず、
その他のアクティブ制御装置(例えば振動低減用のアク
ティブ制御装置)やエコーキャンセラ、ノイズキャンセ
ラへの適応も可能である。
In the present embodiment, the case where the present invention is applied to the active muffling device for muffling the inside of the duct has been described. However, the invention according to claims 1 to 3 is not limited to this.
It is also possible to adapt to other active control devices (for example, active control devices for vibration reduction), echo cancellers, and noise cancellers.

【0024】[0024]

【発明の効果】以上説明してきたように、本発明のアク
ティブ制御装置によれば以下に述べるような効果が発揮
される。請求項1記載の発明によれば、適応アルゴリズ
ム実行手段(44)に備えられたサンプリング回数カウント
手段(5) により、誤差検出手段(3c)からの誤差信号と波
形検出手段(3a)からの波形信号との検出サンプリング回
数をカウントし、このサンプリング回数が所定値に達し
たときにサンプリング回数カウント手段(5) からフィル
タ更新信号を出力すると共に、前記適応アルゴリズム実
行手段(44)に備えられた前記サンプリング回数カウント
手段(5) からのフィルタ更新信号が受信可能な係数ベク
トル加算手段(6) により、前記誤差検出手段(3c)からの
誤差信号と波形検出手段(3a)からの波形信号との検出サ
ンプリング毎にデジタルフィルタ(41)のフィルタ係数を
更新するための係数ベクトルを計算して加算し、係数ベ
クトル加算手段(6) が、前記サンプリング回数カウント
手段(5) からのフィルタ更新信号を受けたときに、この
加算された係数ベクトルに基いて前記デジタルフィルタ
(41)のフィルタ係数を更新させるようにしたために、係
数ベクトル加算手段(6) における加算により係数ベクト
ルの収束誤差を相殺されて該誤差を低減させることがで
き、計算量を従来のLMSアルゴリズムと同等にしなが
らフィルタ係数の収束性の向上が図れる方法を得ること
ができる。
As described above, according to the active control device of the present invention, the following effects are exhibited. According to the invention of claim 1, the sampling frequency counting means (5) provided in the adaptive algorithm executing means (44) causes the error signal from the error detecting means (3c) and the waveform from the waveform detecting means (3a). Counts the number of detection samplings with a signal, outputs a filter update signal from the sampling number counting means (5) when the sampling number reaches a predetermined value, and the adaptive algorithm executing means (44) is provided with Detection of the error signal from the error detection means (3c) and the waveform signal from the waveform detection means (3a) by the coefficient vector addition means (6) capable of receiving the filter update signal from the sampling number counting means (5) The coefficient vector for updating the filter coefficient of the digital filter (41) for each sampling is calculated and added, and the coefficient vector adding means (6) When receiving the filter update signal from the count means (5), the digital filter based on the summed coefficient vector
Since the filter coefficient of (41) is updated, the convergence error of the coefficient vector can be canceled by the addition in the coefficient vector addition means (6) and the error can be reduced, and the calculation amount can be reduced by the conventional LMS algorithm. It is possible to obtain a method capable of improving the convergence of the filter coefficient while making the same.

【0025】請求項2記載の発明によれば、上記請求項
1記載の発明に係る方法を利用したアクティブ制御装置
により、収束誤差を低減することができるアルゴリズム
を利用した装置を得ることができる。
According to the second aspect of the present invention, the active control apparatus utilizing the method according to the first aspect of the present invention makes it possible to obtain an apparatus utilizing an algorithm capable of reducing the convergence error.

【0026】請求項3記載の発明によれば、サンプリン
グ回数カウント手段(5) を、サンプリング回数がデジタ
ルフィルタ(41)のフィルタ係数の数に一致したときにフ
ィルタ更新信号を出力するようにしたことにより、制御
の発散を招くことなしに高い収束性を確保することがで
きる。
According to the third aspect of the present invention, the sampling number counting means (5) outputs the filter update signal when the sampling number matches the number of filter coefficients of the digital filter (41). As a result, high convergence can be secured without causing divergence of control.

【0027】請求項4記載の発明によれば、上述した請
求項2記載の発明に係る装置を消音装置に利用したこと
により、計算量を従来のLMSアルゴリズムと同等にし
ながら収束性の高いフィルタ係数を得ることができ、装
置の消音性能の向上を図ることができる。
According to the invention described in claim 4, since the device according to the invention described in claim 2 is used as a silencer, a filter coefficient having a high convergence while making the amount of calculation equivalent to that of the conventional LMS algorithm. Therefore, the sound deadening performance of the device can be improved.

【図面の簡単な説明】[Brief description of drawings]

【図1】本発明の構成を示すブロック図である。FIG. 1 is a block diagram showing a configuration of the present invention.

【図2】アクティブ消音装置全体の概略構成を示す図で
ある。
FIG. 2 is a diagram showing a schematic configuration of the entire active silencer.

【図3】適応アルゴリズム実行回路における信号処理の
手順を示すフローチャート図である。
FIG. 3 is a flowchart showing a procedure of signal processing in an adaptive algorithm execution circuit.

【図4】本発明の効果を確認するために行った解析に使
用した消音システム構造の概略図である。
FIG. 4 is a schematic diagram of a silencing system structure used in the analysis performed to confirm the effect of the present invention.

【図5】Aタイプの未知システムのインパルス応答を示
す図である。
FIG. 5 is a diagram showing an impulse response of an unknown system of type A.

【図6】Bタイプの未知システムのインパルス応答を示
す図である。
FIG. 6 is a diagram showing an impulse response of a B type unknown system.

【図7】Aタイプの未知システムにおいてシステム同定
を行った場合の係数更新繰返し回数と騒音レベルとの関
係を示す図である。
FIG. 7 is a diagram showing a relationship between a coefficient update repetition number and a noise level when system identification is performed in an A type unknown system.

【図8】Bタイプの未知システムにおいてシステム同定
を行った場合の係数更新繰返し回数と騒音レベルとの関
係を示す図である。
FIG. 8 is a diagram showing the relationship between the number of times coefficient updating is repeated and the noise level when system identification is performed in a B type unknown system.

【図9】未知システムを経時的に変更した場合の係数更
新繰返し回数と騒音レベルとの関係を示す図である。
FIG. 9 is a diagram showing the relationship between the coefficient update repetition number and the noise level when the unknown system is changed over time.

【図10】1回のフィルタ係数の更新が実行されるサン
プリング数と収束時の誤差との関係を示す図である。
FIG. 10 is a diagram showing the relationship between the number of samplings at which one update of the filter coefficient is executed and the error at the time of convergence.

【図11】各周波数帯での音圧レベルを従来のものと比
較した図である。
FIG. 11 is a diagram comparing the sound pressure level in each frequency band with a conventional one.

【符号の説明】[Explanation of symbols]

(21) 空間 (3a) 検出マイクロフォン(波形検出手段) (3b) スピーカ(信号出力手段) (3c) モニタマイクロフォン(誤差検出手段) (41) 適応型FIRフィルタ(デジタルフィルタ) (44) 適応アルゴリズム実行回路(適応アルゴリズ
ム実行手段) (5) サンプリング回数カウント手段 (6) 係数ベクトル加算手段
(21) Space (3a) Detection microphone (waveform detection means) (3b) Speaker (signal output means) (3c) Monitor microphone (error detection means) (41) Adaptive FIR filter (digital filter) (44) Adaptive algorithm execution Circuit (adaptive algorithm execution means) (5) Sampling count counting means (6) Coefficient vector addition means

───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.5 識別記号 庁内整理番号 FI 技術表示箇所 H03H 21/00 7037−5J ─────────────────────────────────────────────────── ─── Continuation of the front page (51) Int.Cl. 5 Identification code Office reference number FI technical display location H03H 21/00 7037-5J

Claims (4)

【特許請求の範囲】[Claims] 【請求項1】 波形検出手段(3a)により変動波形信号を
検出し、該波形検出手段(3a)から出力される波形信号に
対して、該波形信号に基いて生成された制御信号をデジ
タルフィルタ(41)から出力し、該制御信号に対応した作
動信号を信号出力手段(3b)から出力した後、該信号出力
手段(3b)が出力した作動信号と前記変動波形信号との誤
差を誤差検出手段(3c)により検出して、該誤差検出手段
(3c)から出力される誤差信号と前記波形検出手段(3a)か
ら出力される波形信号とを適応アルゴリズム実行手段(4
4)が受け、該適応アルゴリズム実行手段(44)により前記
誤差が小さくなるように前記デジタルフィルタ(41)のフ
ィルタ係数を更新するようにしたアクティブ制御方法に
おいて、 前記適応アルゴリズム実行手段(44)に備えられたサンプ
リング回数カウント手段(5) により、前記誤差検出手段
(3c)からの誤差信号と波形検出手段(3a)からの波形信号
との検出サンプリング回数をカウントし、このサンプリ
ング回数が所定値に達したときにサンプリング回数カウ
ント手段(5) からフィルタ更新信号を出力すると共に、
前記適応アルゴリズム実行手段(44)に備えられた前記サ
ンプリング回数カウント手段(5) からのフィルタ更新信
号が受信可能な係数ベクトル加算手段(6) により、前記
誤差検出手段(3c)からの誤差信号と波形検出手段(3a)か
らの波形信号との検出サンプリング毎に前記デジタルフ
ィルタ(41)のフィルタ係数を更新するための係数ベクト
ルを計算して加算し、係数ベクトル加算手段(6) が、前
記サンプリング回数カウント手段(5) からのフィルタ更
新信号を受けたときに、この加算された係数ベクトルに
基いて前記デジタルフィルタ(41)のフィルタ係数を更新
させるようにしたことを特徴とするアクティブ制御方
法。
1. A waveform detecting means (3a) detects a fluctuation waveform signal, and a control signal generated based on the waveform signal is digitally filtered with respect to the waveform signal output from the waveform detecting means (3a). (41), the operation signal corresponding to the control signal is output from the signal output means (3b), and then the error between the operation signal output by the signal output means (3b) and the fluctuation waveform signal is detected. The error detection means by detecting by means (3c)
The error signal output from (3c) and the waveform signal output from the waveform detection means (3a) are adapted to adaptive algorithm execution means (4
4) received, in the active control method adapted to update the filter coefficient of the digital filter (41) so that the error is reduced by the adaptive algorithm execution means (44), in the adaptive algorithm execution means (44) By the sampling number counting means (5) provided, the error detecting means
The number of detection samplings of the error signal from (3c) and the waveform signal from the waveform detection means (3a) is counted, and when the number of samplings reaches a predetermined value, the filter update signal is output from the sampling number counting means (5). While outputting
By the coefficient vector addition means (6) capable of receiving the filter update signal from the sampling number counting means (5) provided in the adaptive algorithm execution means (44), the error signal from the error detection means (3c) and The coefficient vector for updating the filter coefficient of the digital filter (41) is calculated and added for each detection sampling with the waveform signal from the waveform detection means (3a), and the coefficient vector addition means (6) An active control method characterized in that, when receiving a filter update signal from the number counting means (5), the filter coefficient of the digital filter (41) is updated based on the added coefficient vector.
【請求項2】 変動波形信号を検出して波形信号を出力
する波形検出手段(3a)と、 該波形検出手段(3a)からの波形信号に対し、該波形信号
に基いて生成された制御信号を出力するデジタルフィル
タ(41)と、 該デジタルフィルタ(41)より出力された制御信号を受け
て該制御信号に対応した作動信号を出力する信号出力手
段(3b)と、 該信号出力手段(3b)が出力した作動信号と前記変動波形
信号との誤差を検出して誤差信号を出力する誤差検出手
段(3c)と、 該誤差検出手段(3c)から出力された誤差信号と前記波形
検出手段(3a)からの波形信号とを受けて前記誤差が小さ
くなるように前記デジタルフィルタ(41)のフィルタ係数
を更新する適応アルゴリズム実行手段(44)とを備えたア
クティブ制御装置において、 前記適応アルゴリズム実行手段(44)には、前記誤差検出
手段(3c)からの誤差信号と波形検出手段(3a)からの波形
信号との検出サンプリング回数をカウントし、このサン
プリング回数が所定値に達したときにフィルタ更新信号
を出力するサンプリング回数カウント手段(5) と、該サ
ンプリング回数カウント手段(5) からのフィルタ更新信
号が受信可能とされ、前記誤差検出手段(3c)からの誤差
信号と波形検出手段(3a)からの波形信号との検出サンプ
リング毎に前記デジタルフィルタ(41)のフィルタ係数を
更新するための係数ベクトルを計算して加算し、前記サ
ンプリング回数カウント手段(5) からのフィルタ更新信
号を受けたときに、この加算された係数ベクトルに基い
て前記デジタルフィルタ(41)のフィルタ係数を更新させ
る係数ベクトル加算手段(6) とが備えられていることを
特徴とするアクティブ制御装置。
2. A waveform detecting means (3a) for detecting a varying waveform signal and outputting the waveform signal, and a control signal generated based on the waveform signal from the waveform detecting means (3a). A digital filter (41) for outputting the signal, a signal output means (3b) for receiving a control signal output from the digital filter (41) and outputting an operation signal corresponding to the control signal, and the signal output means (3b Error detection means (3c) for detecting an error between the operation signal output by the) and the fluctuation waveform signal and outputting an error signal, and the error signal output from the error detection means (3c) and the waveform detection means (3c). In the active control device comprising an adaptive algorithm executing means (44) for updating the filter coefficient of the digital filter (41) so as to reduce the error by receiving the waveform signal from 3a), the adaptive algorithm executing means In (44), the error detection The sampling frequency counting means for counting the number of detection samplings of the error signal from the stage (3c) and the waveform signal from the waveform detecting means (3a), and outputting the filter update signal when the sampling frequency reaches a predetermined value ( 5) and the filter update signal from the sampling number counting means (5) can be received, and each sampling of the error signal from the error detecting means (3c) and the waveform signal from the waveform detecting means (3a) The coefficient vector for updating the filter coefficient of the digital filter (41) is calculated and added to, and when the filter update signal from the sampling number counting means (5) is received, the added coefficient vector is An active control device comprising: a coefficient vector adding means (6) for updating the filter coefficient of the digital filter (41) based on the above.
【請求項3】 サンプリング回数カウント手段(5) は、
サンプリング回数がデジタルフィルタ(41)のフィルタ係
数の数に一致したときにフィルタ更新信号を出力するよ
うに構成されていることを特徴とする請求項2記載のア
クティブ制御装置。
3. The sampling number counting means (5) comprises:
3. The active control device according to claim 2, wherein the filter update signal is output when the number of samplings matches the number of filter coefficients of the digital filter (41).
【請求項4】 空間(21)の騒音を検出して騒音信号を出
力する検出マイクロフォン(3a)と、 該検出マイクロフォン(3a)からの騒音信号に対して逆位
相で同振幅の反転音信号を生成する適応型FIRフィル
タ(41)と、 該適応型FIRフィルタ(41)より出力された反転音信号
を受けて該反転音信号に基く反転音を前記空間に放射す
るスピーカ(3b)と、 前記空間(21)の所定観測点に配置され、該観測点の低減
音レベルを検出して低減音信号をフィードバック出力す
るモニタマイクロフォン(3c)と、 前記検出マイクロフォン(3a)からの騒音信号及びモニタ
マイクロフォン(3c)からフィードバックされた低減音信
号を入力し、騒音信号及び低減音信号により前記観測点
周辺の音圧レベルを低減するように前記適応型FIRフ
ィルタ(41)のフィルタ係数を更新する適応アルゴリズム
実行回路(44)とを備えたアクティブ制御装置において、 前記適応アルゴリズム実行回路(44)には、前記モニタマ
イクロフォン(3c)からの低減音信号と検出マイクロフォ
ン(3a)からの騒音信号との検出サンプリング回数をカウ
ントし、このサンプリング回数が所定値に達したときに
フィルタ更新信号を出力するサンプリング回数カウント
手段(5) と、該サンプリング回数カウント手段(5) から
のフィルタ更新信号が受信可能とされ、前記モニタマイ
クロフォン(3c)からの低減音信号と検出マイクロフォン
(3a)からの波形信号との検出サンプリング毎に前記適応
型FIRフィルタ(41)のフィルタ係数を更新するための
係数ベクトルを計算して加算し、前記サンプリング回数
カウント手段(5) からのフィルタ更新信号を受けたとき
に、この加算された係数ベクトルに基いて前記適応型F
IRフィルタ(41)のフィルタ係数を更新させる係数ベク
トル加算手段(6) とが備えられていることを特徴とする
アクティブ制御装置。
4. A detection microphone (3a) that detects noise in the space (21) and outputs a noise signal, and a reversal sound signal having the same amplitude and opposite phase to the noise signal from the detection microphone (3a). An adaptive FIR filter (41) for generating, a speaker (3b) for receiving an inverted sound signal output from the adaptive FIR filter (41) and radiating an inverted sound based on the inverted sound signal to the space, A monitor microphone (3c) arranged at a predetermined observation point in the space (21), which detects a reduced sound level at the observation point and outputs a reduced sound signal as feedback, and a noise signal from the detection microphone (3a) and a monitor microphone. The reduced sound signal fed back from (3c) is input, and the adaptive coefficient that updates the filter coefficient of the adaptive FIR filter (41) so as to reduce the sound pressure level around the observation point by the noise signal and the reduced sound signal. In the active control device including the algorithm execution circuit (44), the adaptive algorithm execution circuit (44) detects the reduced sound signal from the monitor microphone (3c) and the noise signal from the detection microphone (3a). A sampling number counting means (5) for counting the number of sampling times and outputting a filter updating signal when the sampling number reaches a predetermined value, and a filter updating signal from the sampling number counting means (5) can be received. , Reduced sound signal from the monitor microphone (3c) and detection microphone
The coefficient vector for updating the filter coefficient of the adaptive FIR filter (41) is calculated and added for each detection sampling with the waveform signal from (3a), and the filter is updated from the sampling number counting means (5). When a signal is received, the adaptive F based on the added coefficient vector
An active control device comprising: a coefficient vector adding means (6) for updating the filter coefficient of the IR filter (41).
JP5110170A 1993-05-12 1993-05-12 Method and device for active control Withdrawn JPH06324686A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP5110170A JPH06324686A (en) 1993-05-12 1993-05-12 Method and device for active control

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP5110170A JPH06324686A (en) 1993-05-12 1993-05-12 Method and device for active control

Publications (1)

Publication Number Publication Date
JPH06324686A true JPH06324686A (en) 1994-11-25

Family

ID=14528830

Family Applications (1)

Application Number Title Priority Date Filing Date
JP5110170A Withdrawn JPH06324686A (en) 1993-05-12 1993-05-12 Method and device for active control

Country Status (1)

Country Link
JP (1) JPH06324686A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0816215A1 (en) 1996-07-03 1998-01-07 Yamaha Hatsudoki Kabushiki Kaisha Electric motor-assisted bicycle

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0816215A1 (en) 1996-07-03 1998-01-07 Yamaha Hatsudoki Kabushiki Kaisha Electric motor-assisted bicycle

Similar Documents

Publication Publication Date Title
JP3089082B2 (en) Adaptive digital filter
JP4790843B2 (en) Active silencer and active silencer method
EP0578212B1 (en) Active control apparatus with an adaptive digital filter
US8189800B2 (en) Active silencer and method for controlling active silencer
JP5707663B2 (en) Active silencer
JP6343970B2 (en) Signal processing device, program, range hood device
WO2016178309A1 (en) Signal processing device, signal processing method, program, and rangehood apparatus
JPH03274897A (en) Electronic muffling method and its device
JPH06324686A (en) Method and device for active control
JP2989379B2 (en) Active control device using adaptive IIR digital filter
JPH0540486A (en) Active muffling device
JP2014235240A (en) Active noise controller
JPH0934470A (en) Adaptive filter
JPH06318083A (en) Active muffler device
JPH0774590A (en) Electronic muffler
JP3671232B2 (en) Echo canceller
JPH0635482A (en) Method and device for active noise elimination
JP6213868B2 (en) Signal processing apparatus, active noise control apparatus, and active noise control method
JP3441753B2 (en) Active vibration control device
JPH06332468A (en) Active silencer
JP2002333887A (en) Transfer function identifying device and active noise eliminating device
JPH06250676A (en) Active controller
JPH07219557A (en) Active noise elimination device
JPH05333870A (en) Active sound eliminating device
JPH04247497A (en) Active type silencing device

Legal Events

Date Code Title Description
A300 Withdrawal of application because of no request for examination

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20000801