JP3517883B2 - Adaptive control device and active noise control device - Google Patents

Adaptive control device and active noise control device

Info

Publication number
JP3517883B2
JP3517883B2 JP06401292A JP6401292A JP3517883B2 JP 3517883 B2 JP3517883 B2 JP 3517883B2 JP 06401292 A JP06401292 A JP 06401292A JP 6401292 A JP6401292 A JP 6401292A JP 3517883 B2 JP3517883 B2 JP 3517883B2
Authority
JP
Japan
Prior art keywords
optimum value
filter
digital filter
filter coefficient
convergence
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP06401292A
Other languages
Japanese (ja)
Other versions
JPH05265467A (en
Inventor
健一郎 村岡
勉 浜辺
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.)
Nissan Motor Co Ltd
Original Assignee
Nissan Motor Co 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 Nissan Motor Co Ltd filed Critical Nissan Motor Co Ltd
Priority to JP06401292A priority Critical patent/JP3517883B2/en
Publication of JPH05265467A publication Critical patent/JPH05265467A/en
Application granted granted Critical
Publication of JP3517883B2 publication Critical patent/JP3517883B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Description

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

【0001】[0001]

【産業上の利用分野】この発明は、制御対象の振る舞い
が未知である場合或いは変動する場合であっても制御が
可能な適応制御装置、並びに騒音源から伝達される騒音
に制御音源から発せられる制御音を干渉させることによ
り騒音の低減を図る能動型騒音制御装置に関し、特に、
制御アクチュエータや制御音源を駆動する信号を生成す
るフィルタ係数可変のディジタルフィルタと、このディ
ジタルフィルタのフィルタ係数を所定の適応アルゴリズ
ムに応じて更新する適応処理部とを備えた装置におい
て、制御対象の状況変化や装置の構成部品の劣化等に基
づいて適応処理部の処理内容が実際の制御対象や装置の
特性に合致しない状態が生じても、ディジタルフィルタ
のフィルタ係数が最適値に素早く収束できるようにした
ものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an adaptive control device capable of controlling even if the behavior of a controlled object is unknown or fluctuates, and the noise transmitted from a noise source is emitted from a control sound source. Regarding an active noise control device for reducing noise by interfering with control sound,
In a device including a digital filter with a variable filter coefficient that generates a signal for driving a control actuator or a control sound source, and an adaptive processing unit that updates the filter coefficient of this digital filter according to a predetermined adaptive algorithm, the state of the controlled object Even if the processing contents of the adaptive processing unit do not match the actual control target or the characteristics of the device due to changes or deterioration of the component parts of the device, the filter coefficient of the digital filter can quickly converge to the optimum value. It was done.

【0002】[0002]

【従来の技術】この種の従来の装置として、英国特許第
2149614号や特公表1−501344号公報等に
記載のものがある。この従来の装置は、航空機の客室や
これに類する閉空間に適用される騒音低減装置であっ
て、閉空間の外部に位置するエンジン等の単一の騒音源
は、基本周波数f0 及びその高調波f1 〜fn を含む騒
音を発生するという条件の下において作動するものであ
る。
2. Description of the Related Art As a conventional device of this type, there are those described in British Patent No. 2149614 and Japanese Patent Publication No. 1-501344. This conventional device is a noise reduction device applied to a cabin of an aircraft or a similar closed space, and a single noise source such as an engine located outside the closed space has a fundamental frequency f 0 and its harmonics. It operates under the condition that noise including the waves f 1 to f n is generated.

【0003】具体的には、閉空間内の複数の位置に設置
され音圧を検出するマイクロフォンと、その閉空間に制
御音を発生する複数のラウドスピーカとを備え、騒音源
の周波数f0 〜fn 成分に基づき、それら周波数f0
n 成分と逆位相の信号でラウドスピーカを駆動させ、
もって閉空間に伝達される騒音と逆位相の制御音をラウ
ドスピーカから発生させて騒音を打ち消している。
[0003] More specifically, includes a microphone for detecting a plurality of the installed sound pressure to a location within the closed space, and a plurality of loudspeakers for generating a control sound to the closed space, the frequency f 0 of the noise source - Based on the f n component, those frequencies f 0 ~
The loudspeaker is driven by a signal having a phase opposite to that of the f n component,
Therefore, a control sound having a phase opposite to that of the noise transmitted to the closed space is generated from the loudspeaker to cancel the noise.

【0004】そして、ラウドスピーカから発せられる制
御音の生成方法として、PROCEEDINGS OF THE IEEE,VOL.
63 PAGE 1692−1975, “ADAPTIVE NOISE CANSELLATION
:PRINCIPLES AND APPLICATIONS ”で述べられている
‘WIDROW LMS’アルゴリズムを、多チャンネルに展開し
たアルゴリズムを適用している。その内容は、上記特許
の発明者による論文、“A MULTIPLE ERROR LMS ALGORIT
HM AND ITS APPLICATIONTO THE ACTIVE CONTROL OF SOU
ND AND VIBRATION ”,IEEE TRANS.ACOUST.,SPEECH,SIGN
AL PROCESSING,VOL.ASSP −35,PP.1423−1434,1987 に
も述べられている。
Then, as a method of generating the control sound emitted from the loudspeaker, PROCEEDINGS OF THE IEEE, VOL.
63 PAGE 1692-1975, “ADAPTIVE NOISE CANSELLATION
: The "WIDROW LMS" algorithm described in "PRINCIPLES AND APPLICATIONS" is applied to multiple channels. The contents are based on the paper "A MULTIPLE ERROR LMS ALGORIT" by the inventor of the above patent.
HM AND ITS APPLICATIONTO THE ACTIVE CONTROL OF SOU
ND AND VIBRATION ”, IEEE TRANS.ACOUST., SPEECH, SIGN
AL PROCESSING, VOL.ASSP −35, PP.1432-1434, 1987.

【0005】即ち、LMSアルゴリズムは、適応型ディ
ジタルフィルタのフィルタ係数を更新するのに好適なア
ルゴリズムの一つであって、例えばいわゆるFilte
red−X LMSアルゴリズムにあっては、ラウドス
ピーカからマイクロフォンまでの音響伝達特性を表すフ
ィルタ(伝達関数モデル)を全てのラウドスピーカとマ
イクロフォンとの組み合わせについて設定し、騒音源の
騒音発生状態を表す基準信号をそのフィルタで処理した
値と、各マイクロフォンが検出した残留騒音とに基づい
た所定の評価関数の値が低減するように、各ラウドスピ
ーカ毎に設けられたフィルタ係数可変のディジタルフィ
ルタのフィルタ係数を更新している。
That is, the LMS algorithm is one of the algorithms suitable for updating the filter coefficient of the adaptive digital filter, and is, for example, the so-called Filter.
In the red-X LMS algorithm, a filter (transfer function model) that represents the acoustic transfer characteristic from the loudspeaker to the microphone is set for all combinations of the loudspeaker and the microphone, and is a standard that represents the noise generation state of the noise source. A filter coefficient of a variable digital filter provided for each loudspeaker so that the value of a predetermined evaluation function based on the value obtained by processing the signal with the filter and the residual noise detected by each microphone is reduced. Have been updated.

【0006】[0006]

【発明が解決しようとする課題】確かに、上記のような
従来の装置であれば、ラウドスピーカを駆動する信号を
生成するディジタルフィルタのフィルタ係数が制御の進
行に伴って最適値に収束していくから、最終的には閉空
間内の騒音は低減するが、閉空間内の温度や湿度等の環
境や、ラウドスピーカやマイクロフォン等の劣化に伴
い、実空間の音響伝達特性が変化すると、ディジタルフ
ィルタのフィルタ係数を更新するアルゴリズムに必要な
伝達関数モデルの精度が相対的に低下し、ディジタルフ
ィルタのフィルタ係数の収束に時間がかかり適応動作が
遅くなるという不具合がある。
Certainly, in the conventional device as described above, the filter coefficient of the digital filter for generating the signal for driving the loudspeaker converges to the optimum value as the control progresses. Eventually, the noise in the closed space will be reduced, but if the acoustic transfer characteristics in the real space change due to changes in the environment such as temperature and humidity in the closed space, or deterioration of the loudspeaker or microphone, digital noise will change. There is a problem that the accuracy of the transfer function model required for the algorithm for updating the filter coefficient of the filter is relatively lowered, and it takes time to converge the filter coefficient of the digital filter, and the adaptive operation is delayed.

【0007】このような不具合を解決するのに、例え
ば、伝達関数モデルについての多次元のマップを、閉空
間内の温度や湿度、ラウドスピーカやマイクロフォンの
劣化状況等をパラメータとして構築しておき、閉空間内
の状況の変化等に応じて伝達関数モデルを適宜更新する
という方法が考えられるが、これでは、ラウドスピーカ
やマイクロフォンの数の増加に伴ってマップとして記憶
すべき伝達関数モデルの数が膨大になり、必要な記憶容
量が極端に大きくなって実用的でなく、しかも、予め設
定した伝達関数モデルを実際の音響空間の伝達関数に完
全に一致させることは困難である。
In order to solve such a problem, for example, a multidimensional map of a transfer function model is constructed with parameters such as temperature and humidity in a closed space and deterioration conditions of loudspeakers and microphones. A method of updating the transfer function model appropriately according to changes in the situation in the closed space can be considered, but with this, as the number of loudspeakers and microphones increases, the number of transfer function models to be stored as a map is increased. It becomes enormous and the storage capacity required becomes extremely large, which is not practical, and it is difficult to perfectly match a preset transfer function model with the actual transfer function of the acoustic space.

【0008】また、ラウドスピーカ及びマイクロフォン
間の音響伝達特性に影響を与える因子(例えば、閉空間
内の温度、窓の開閉状態、存在する人間の数等)が変化
した場合に、同定音を発生して音響伝達特性を測定し、
これに基づいて伝達関数モデルを更新するという方法も
考えられるが、これでは、同定を行う際に空間内に存在
する人間に同定音が聞かれてしまい不快感を与えるとい
う欠点がある。
An identification sound is generated when a factor that affects the acoustic transfer characteristics between the loudspeaker and the microphone (for example, the temperature in the closed space, the open / closed state of windows, the number of people present, etc.) changes. To measure the acoustic transfer characteristics,
A method of updating the transfer function model on the basis of this is also conceivable, but this has a drawback in that the identification sound is heard by a person existing in the space when performing identification, which causes discomfort.

【0009】本発明は、このような従来の技術が有する
未解決の課題に着目してなされたものであって、記憶容
量の極端な増加や人間に不快感を与えるという不具合を
招くことなく、制御対象の状況変化や装置の構成部品の
劣化等に基づいて適応処理部の処理内容が実際の制御対
象や装置の特性に合致しない状態が生じても、ディジタ
ルフィルタのフィルタ係数を最適値に素早く収束させる
ことができる適応処理装置及び能動型騒音制御装置を提
供することを目的としている。
The present invention has been made by paying attention to the unsolved problems of the prior art as described above, and does not cause a problem such as an extreme increase in storage capacity and an unpleasant feeling to humans. Even if the processing content of the adaptive processing unit does not match the actual control target or device characteristics due to changes in the status of the control target or deterioration of the components of the device, the filter coefficient of the digital filter can be quickly set to the optimum value. An object of the present invention is to provide an adaptive processing device and an active noise control device that can be converged.

【0010】[0010]

【課題を解決するための手段】上記目的を達成するため
に、請求項1記載の発明は、所定の基準信号に応じて制
御アクチュエータを駆動する信号を生成するフィルタ係
数可変のディジタルフィルタと、所定の評価関数の値が
目標値に一致するように前記ディジタルフィルタのフィ
ルタ係数を更新する適応処理手段と、を備えた適応制御
装置において、前記ディジタルフィルタのフィルタ係数
が最適値に収束したか否かを判定する最適値収束判定手
段と、この最適値収束判定手段が最適値に収束したと判
定した際の前記ディジタルフィルタのフィルタ係数から
判る最適値への最短の収束経路、及び最適値に収束し
ていない前記ディジタルフィルタのフィルタ係数から判
る実際の収束経路に基づいて、前記実際の収束経路が前
記最短の収束経路に近づくように、前記適応処理手段の
処理内容を補正する補正手段と、を設けた。
In order to achieve the above object, the invention according to claim 1 is a digital filter having a variable filter coefficient for generating a signal for driving a control actuator according to a predetermined reference signal, and a predetermined filter. In the adaptive control device, comprising: an adaptive processing unit that updates the filter coefficient of the digital filter so that the value of the evaluation function of ## EQU1 ## matches the target value. and the optimum value convergence determining means for determining, from the filter coefficient of the digital filter when the optimum value convergence determining means is determined to have converged to the optimum value
Shortest convergence path to the optimum value can be seen, and, determine the filter coefficients of the digital filter does not converge to the optimum value
Based on the actual convergence path that the actual convergence path before
The correction means for correcting the processing content of the adaptive processing means is provided so as to approach the shortest convergence path .

【0011】また、請求項2記載の発明は、上記請求項
1記載の発明において、補正手段は、前記実際の収束経
路が前記最短の収束経路に近づくように、ディジタルフ
ィルタのフィルタ係数の更新ゲインを変更するゲイン調
整手段を有するものである。そして、請求項3記載の発
明は、上記請求項1又は請求項2記載の発明において、
補正手段は、最適値収束判定手段が最適値に収束したと
判定した際のディジタルフィルタのフィルタ係数が変動
している場合には、適応処理手段の処理内容を補正しな
いこととした。
According to a second aspect of the present invention, in the above-mentioned first aspect of the invention, the correction means is the actual convergence process.
It has a gain adjusting means for changing the update gain of the filter coefficient of the digital filter so that the path approaches the shortest convergence path . And the invention of claim 3 is the same as the invention of claim 1 or claim 2,
The correction means does not correct the processing content of the adaptive processing means when the filter coefficient of the digital filter changes when the optimum value convergence determination means determines that the optimum value has converged to the optimum value.

【0012】さらに、請求項4記載の発明は、上記請求
項1乃至請求項3記載の発明において、補正手段は、最
適値収束判定手段が最適値に収束したと判定した際のデ
ィジタルフィルタのフィルタ係数及び制御開始時の前記
ディジタルフィルタのフィルタ係数から決まる方向と、
最適値に収束していない前記ディジタルフィルタのフィ
ルタ係数及び制御開始時の前記ディジタルフィルタのフ
ィルタ係数から決まる方向との間のずれが無くなるよう
適応処理手段の処理内容を補正する。
Further, in the invention described in claim 4, in the invention described in any one of claims 1 to 3, the correction means is a filter of the digital filter when the optimum value convergence judgment means judges that the optimum value has converged to the optimum value. A direction determined by the coefficient and the filter coefficient of the digital filter at the start of control,
To eliminate the deviation between the filter coefficient of the digital filter that has not converged to the optimum value and the direction determined by the filter coefficient of the digital filter at the start of control.
Then, the processing content of the adaptive processing means is corrected.

【0013】一方、上記目的を達成するために、請求項
5記載の発明は、騒音源から騒音が伝達される空間に制
御音を発生可能な制御音源と、前記騒音源の騒音発生状
態を検出し基準信号として出力する騒音発生状態検出手
段と、前記基準信号に応じて前記制御音源を駆動する信
号を生成するフィルタ係数可変のディジタルフィルタ
と、前記空間内の所定位置における残留騒音を検出する
残留騒音検出手段と、前記残留騒音に基づいた所定の評
価関数の値が低減するように前記ディジタルフィルタの
フィルタ係数を更新する適応処理手段と、を備えた能動
型騒音制御装置において、前記ディジタルフィルタのフ
ィルタ係数が最適値に収束したか否かを判定する最適値
収束判定手段と、この最適値収束判定手段が最適値に収
束したと判定した際の前記ディジタルフィルタのフィル
タ係数から判る最適値への最短の収束経路、及び最適
値に収束していない前記ディジタルフィルタのフィルタ
係数から判る実際の収束経路に基づいて、前記実際の収
束経路が前記最短の収束経路に近づくように、前記適応
処理手段の処理内容を補正する補正手段と、を設けた。
On the other hand, in order to achieve the above object, the invention according to claim 5 detects a control sound source capable of generating a control sound in a space where noise is transmitted from the noise source, and a noise generation state of the noise source. Noise generating state detecting means for outputting as a reference signal, a digital filter having a variable filter coefficient for generating a signal for driving the control sound source according to the reference signal, and a residual noise detecting residual noise at a predetermined position in the space. In an active noise control device comprising noise detection means and adaptive processing means for updating the filter coefficient of the digital filter so that the value of a predetermined evaluation function based on the residual noise is reduced, Optimal value convergence determining means for determining whether or not the filter coefficient has converged to an optimal value, and when it is determined that the optimal value convergence determining means has converged to an optimal value. Shortest convergence path to the optimum value can be seen from the filter coefficients of the digital filters, and, based on the actual convergence path seen from the filter coefficient of the digital filter does not converge to the optimum value, the actual yield
Correction means for correcting the processing content of the adaptive processing means so that the bundle path approaches the shortest convergence path .

【0014】また、請求項6記載の発明は、上記請求項
5記載の発明において、適応処理手段は、制御音源及び
残留騒音検出手段間の音響伝達特性に相当する伝達関数
モデルを有し、補正手段は、前記実際の収束経路が前記
最短の収束経路に近づくように、前記伝達関数モデルを
補正するものである。そして、請求項7記載の発明は、
上記請求項6記載の発明において、補正手段は、適応処
理手段に含まれる伝達関数モデルと直列に配設された伝
達特性可変の伝達関数補正フィルタと、最適値収束判定
手段が最適値に収束したと判定した際の前記ディジタル
フィルタのフィルタ係数及び最適値に収束していない前
記ディジタルフィルタのフィルタ係数に基づいて、前記
実際の収束経路が前記最短の収束経路に近づくように、
前記伝達関数補正フィルタの伝達特性を変更する伝達特
性制御手段と、を有するものである。
According to the invention of claim 6, in the invention of claim 5, the adaptive processing means has a transfer function model corresponding to an acoustic transfer characteristic between the control sound source and the residual noise detecting means, The means is that the actual convergence path is
The transfer function model is corrected so as to approach the shortest convergence path . The invention according to claim 7 is
In the invention described in claim 6, the correction means is a transfer function correction filter arranged in series with the transfer function model included in the adaptive processing means, the transfer characteristic variable filter, and the optimum value convergence determination means converges to an optimum value. Based on the filter coefficient of the digital filter and the filter coefficient of the digital filter that has not converged to the optimum value when it is determined that
As the actual convergence path approaches the shortest convergence path,
Transfer characteristic control means for changing the transfer characteristic of the transfer function correction filter.

【0015】さらに、請求項8記載の発明は、上記請求
項6記載の発明において、補正手段は、ディジタルフィ
ルタと直列に配設された伝達特性可変の伝達関数補正フ
ィルタと、最適値収束判定手段が最適値に収束したと判
定した際の前記ディジタルフィルタのフィルタ係数及び
最適値に収束していない前記ディジタルフィルタのフィ
ルタ係数に基づいて、前記実際の収束経路が前記最短の
収束経路に近づくように、前記伝達関数補正フィルタの
伝達特性を変更する伝達特性制御手段と、を有するもの
である。
Further, in the invention described in claim 8, in the invention described in claim 6, the correction means is a transfer function correction filter having a variable transfer characteristic arranged in series with a digital filter, and an optimum value convergence determination means. Based on the filter coefficient of the digital filter when it is determined that has converged to the optimum value and the filter coefficient of the digital filter that has not converged to the optimum value, the actual convergence path is the shortest.
Transfer characteristic control means for changing the transfer characteristic of the transfer function correction filter so as to approach the convergent path .

【0016】一方、請求項9記載の発明は、上記請求項
5記載の発明において、補正手段は、前記実際の収束経
路が前記最短の収束経路に近づくように、ディジタルフ
ィルタのフィルタ係数の更新ゲインを変更するゲイン調
整手段を有するものである。また、請求項10記載の発
明は、上記請求項5乃至請求項9記載の発明において、
補正手段は、最適値収束判定手段が最適値に収束したと
判定した際のディジタルフィルタのフィルタ係数が変動
している場合には、適応処理手段の処理内容を補正しな
いこととした。
On the other hand, in the invention described in claim 9, in the invention described in claim 5, the correction means is configured to perform the actual convergence process.
It has a gain adjusting means for changing the update gain of the filter coefficient of the digital filter so that the path approaches the shortest convergence path . The invention according to claim 10 is the invention according to any one of claims 5 to 9,
The correction means does not correct the processing content of the adaptive processing means when the filter coefficient of the digital filter changes when the optimum value convergence determination means determines that the optimum value has converged to the optimum value.

【0017】そして、請求項11記載の発明は、上記請
求項5乃至請求項10記載の発明において、補正手段
は、最適値収束判定手段が最適値に収束したと判定した
際のディジタルフィルタのフィルタ係数及び制御開始時
の前記ディジタルフィルタのフィルタ係数から決まる方
向と、制御開始直後の前記ディジタルフィルタのフィル
タ係数及び制御開始時の前記ディジタルフィルタのフィ
ルタ係数から決まる方向との間のずれが無くなるように
適応処理手段の処理内容を補正する。
According to an eleventh aspect of the present invention, in the invention according to the fifth to tenth aspects, the correction means is a filter of the digital filter when the optimum value convergence determination means determines that the optimum value has converged to the optimum value. The deviation between the direction determined by the coefficient and the filter coefficient of the digital filter at the start of control and the direction determined by the filter coefficient of the digital filter immediately after the start of control and the filter coefficient of the digital filter at the start of control is eliminated. The processing contents of the adaptive processing means are corrected.

【0018】[0018]

【作用】請求項1記載の発明にあっては、ディジタルフ
ィルタが所定の基準信号に応じて制御アクチュエータを
駆動する信号を生成すると、これによって制御アクチュ
エータが駆動され、所定の制御が実行されるが、制御開
始直後は、そのディジタルフィルタのフィルタ係数が最
適値に収束しているとは限らないので、必ずしも制御ア
クチュエータによって良好な制御が実行されるとはいえ
ない。
According to the first aspect of the present invention, when the digital filter generates a signal for driving the control actuator in response to a predetermined reference signal, the control actuator is driven by this signal and predetermined control is executed. Immediately after the start of control, the filter coefficient of the digital filter does not always converge to the optimum value, and therefore it cannot be said that good control is always executed by the control actuator.

【0019】しかし、適応処理手段が、所定の評価関数
が目標値に一致するようにディジタルフィルタのフィル
タ係数を更新するので、そのディジタルフィルタのフィ
ルタ係数が最適値に収束していき、制御アクチュエータ
によって良好な制御が実行されるようになる。つまり、
適応処理手段がディジタルフィルタのフィルタ係数を適
宜更新する結果、制御実行前は最適なディジタルフィル
タのフィルタ係数が未知であっても、最終的には良好な
制御が行えるフィルタ係数を得ることができる。
However, since the adaptive processing means updates the filter coefficient of the digital filter so that the predetermined evaluation function coincides with the target value, the filter coefficient of the digital filter gradually converges to the optimum value, and the control actuator operates. Good control will be performed. That is,
As a result of the adaptive processing means appropriately updating the filter coefficient of the digital filter, even if the optimum filter coefficient of the digital filter is unknown before the control is executed, it is possible to finally obtain the filter coefficient that allows good control.

【0020】そして、最適値に収束したディジタルフィ
ルタのフィルタ係数から最適値への最短の収束経路が判
り、また、最適値に収束していないディジタルフィルタ
のフィルタ係数から実際の収束経路が判ることから、最
適値収束判定手段が最適値に収束したと判定した際のデ
ィジタルフィルタのフィルタ係数と、最適値に収束して
いないディジタルフィルタのフィルタ係数とに基づい
て、最短の収束経路に対する実際の収束経路のずれが判
る。
Since the shortest convergence path from the filter coefficient of the digital filter that converges to the optimum value to the optimum value can be known, and the actual convergence path can be known from the filter coefficient of the digital filter that has not converged to the optimum value. Based on the filter coefficient of the digital filter when the optimum value convergence determination means determines that the optimum value has converged to the optimum value and the filter coefficient of the digital filter that has not converged to the optimum value, the actual convergence path for the shortest convergence path. You can see the deviation.

【0021】そこで、最適値収束判定手段が最適値に収
束したと判定した際のディジタルフィルタのフィルタ係
から判る最適値への最短の収束経路と、最適値に収束
していないディジタルフィルタのフィルタ係数から判る
実際の収束経路とに基づいて、補正手段が適応処理手段
の処理内容を実際の収束経路が最短の収束経路に近づく
ように補正する結果、その後の適応処理手段の処理によ
るディジタルフィルタのフィルタ係数の収束経路は、補
正前の収束経路に比べて最短の収束経路に近くなる。
Therefore, the shortest convergence path to the optimum value which can be known from the filter coefficient of the digital filter when the optimum value convergence judgment means judges that the optimum value has converged to the optimum value, and the filter coefficient of the digital filter which has not converged to the optimum value. Understand from
Based on the actual convergence path , the correction means causes the processing content of the adaptive processing means to approach the convergence path where the actual convergence path is the shortest.
As a result of such correction, the convergence path of the filter coefficient of the digital filter by the subsequent processing of the adaptive processing means becomes closer to the shortest convergence path as compared with the convergence path before correction.

【0022】また、最短の収束経路に対する実際の収束
経路のずれが大きい場合に、ディジタルフィルタのフィ
ルタ係数の更新ゲインが大きいと、フィルタ係数の更新
に伴ってそれらの間のずれがさらに大きくなるおそれが
あるから、請求項2記載の発明のように、補正手段とし
て、ディジタルフィルタのフィルタ係数の更新ゲインを
実際の収束経路が最短の収束経路に近づくように変更す
るゲイン調整手段を有し、例えば、最短の収束経路に対
する実際の収束経路のずれが大きい場合に更新ゲインを
小さくすれば、収束経路のずれの拡大が防がれる。
Further, if the deviation of the actual convergence path from the shortest convergence path is large and the update gain of the filter coefficient of the digital filter is large, the deviation between them may become larger as the filter coefficient is updated. Therefore, as in the invention according to claim 2, as the correction means, the update gain of the filter coefficient of the digital filter is set.
There is a gain adjusting means for changing the actual convergence path so that it approaches the shortest convergence path . For example, if the deviation of the actual convergence path from the shortest convergence path is large, the deviation of the convergence path will be reduced if the update gain is reduced. Can be prevented from expanding.

【0023】なお、最適値収束判定手段が最適値に収束
したと判定した際のディジタルフィルタのフィルタ係数
が変動している場合は、最短の収束経路も正確に判明で
きないので、請求項3記載の発明のように、適応処理手
段の処理内容を補正しないこととすれば、適応処理手段
の処理内容が不明確な情報に基づいて変動してしまうよ
うなことが避けられる。
When the optimum value convergence determination means determines that the filter coefficient of the digital filter has converged to the optimum value, the shortest convergence path cannot be accurately determined. If the processing content of the adaptive processing means is not corrected as in the invention, it is possible to avoid that the processing content of the adaptive processing means fluctuates based on unclear information.

【0024】そして、最適値収束判定手段が最適値に収
束したと判定した際のディジタルフィルタのフィルタ係
数と、制御開始時のディジタルフィルタのフィルタ係数
(つまり、フィルタ係数の初期値)とから、最短の収束
経路の方向が決まり、また、最適値に収束していないデ
ィジタルフィルタのフィルタ係数(例えば、制御開始直
後のフィルタ係数)と、制御開始時のディジタルフィル
タのフィルタ係数とから、実際の収束経路の方向が判る
から、請求項4記載の発明のように、それら方向の間の
ずれが無くなるように補正手段が適応処理手段の処理内
容を補正すれば、最短の収束経路と実際の収束経路との
間のずれが無くなるような補正が行なわれる。
Then, from the filter coefficient of the digital filter when the optimum value convergence determining means determines that it has converged to the optimum value and the filter coefficient of the digital filter at the start of control (that is, the initial value of the filter coefficient), the shortest The direction of the convergence path is determined, and the actual convergence path is determined from the filter coefficient of the digital filter that has not converged to the optimum value (for example, the filter coefficient immediately after the start of control) and the filter coefficient of the digital filter at the start of control. Therefore, if the correction means corrects the processing content of the adaptive processing means so that the deviation between the directions is eliminated , the shortest convergence path and the actual convergence path can be obtained. The correction is performed so that the gap between the two is eliminated .

【0025】一方、請求項5記載の発明にあっては、騒
音発生状態検出手段から出力された基準信号に応じて、
フィルタ係数可変のディジタルフィルタが制御音源を駆
動する信号を生成するから、制御音源からは、騒音源に
発生する騒音に相関のある制御音が発生するが、制御開
始直後は、フィルタ係数が最適な値に収束しているとは
限らないので、必ずしも騒音が低減されるとはいえな
い。
On the other hand, according to the invention of claim 5, in accordance with the reference signal output from the noise generation state detecting means,
Since a digital filter with a variable filter coefficient generates a signal for driving the control sound source, the control sound source generates a control sound that is correlated with the noise generated in the noise source. The noise is not necessarily reduced because it does not necessarily converge to the value.

【0026】しかし、残留騒音検出手段が検出した残留
騒音に基づいた所定の評価関数の値が低減するように、
適応処理手段がディジタルフィルタのフィルタ係数を更
新するので、そのフィルタ係数が最適値に収束してい
き、制御音源から発せられる制御音によって騒音が打ち
消され、空間内の騒音が低減する。そして、最適値に収
束したディジタルフィルタのフィルタ係数から最適値へ
の最短の収束経路が判り、また、最適値に収束していな
いディジタルフィルタのフィルタ係数から実際の収束経
路が判ることから、補正手段が、最適値収束判定手段が
最適値に収束したと判定した際のディジタルフィルタの
フィルタ係数から判る最適値への最短の収束経路及び最
適値に収束していないディジタルフィルタのフィルタ係
から判る実際の収束経路に基づいて、適応処理手段の
処理内容を実際の収束経路が最短の収束経路に近づくよ
うに補正する結果、その後に例えば騒音の周波数の変動
に伴って最適値が変化して再び適応処理手段による処理
が行われると、ディジタルフィルタのフィルタ係数の収
束経路は、補正前の収束経路に比べて最短の収束経路に
近くなり、フィルタ係数は比較的早く最適値に収束する
ようになる。
However, in order to reduce the value of the predetermined evaluation function based on the residual noise detected by the residual noise detecting means,
Since the adaptive processing means updates the filter coefficient of the digital filter, the filter coefficient gradually converges to the optimum value, the control sound emitted from the control sound source cancels the noise, and the noise in the space is reduced. The shortest convergence path from the filter coefficient of the digital filter that converges to the optimum value to the optimum value is known, and the actual convergence path is known from the filter coefficient of the digital filter that has not converged to the optimum value. However, when the optimum value convergence determination means determines that the optimum value has converged to the optimum value, the shortest convergence path to the optimum value can be found from the filter coefficient of the digital filter, and the actual filter coefficient of the digital filter that has not converged to the optimum value can be found. Based on the convergence path , the processing contents of the adaptive processing means are adjusted so that the actual convergence path approaches the shortest convergence path.
Results urchin correction, the subsequent processing by the adaptive processing unit once again the optimum value is changed in accordance with the example the variation of the frequency of the noise is carried out, the convergence path of the filter coefficient of the digital filter is compared with a convergence path before correction Becomes closer to the shortest convergence path, and the filter coefficient converges to the optimum value relatively quickly.

【0027】ここで、適応処理手段が、制御音源及び残
留騒音検出手段間の音響伝達特性に相当する伝達関数モ
デルを有する場合、ディジタルフィルタのフィルタ係数
の適応動作が遅くなる主な原因は、実際の音響伝達特性
が変化して伝達関数モデルの精度が相対的に低下してし
まうことである。従って、請求項6記載の発明のよう
に、補正手段が適応処理手段に含まれる伝達関数モデル
実際の収束経路が最短の収束経路に近づくように補正
することとすれば、最短の収束経路に対する実際の収束
経路のずれが、的確に適応処理手段の処理内容の補正に
反映する。
Here, when the adaptive processing means has a transfer function model corresponding to the acoustic transfer characteristic between the control sound source and the residual noise detecting means, the main cause of the delay in the adaptive operation of the filter coefficient of the digital filter is actually That is, the acoustic transfer characteristics of the above change, and the accuracy of the transfer function model relatively decreases. Therefore, if the correction means corrects the transfer function model included in the adaptive processing means so that the actual convergence path approaches the shortest convergence path, the correction means for the shortest convergence path The deviation of the actual convergence path is accurately reflected in the correction of the processing content of the adaptive processing means.

【0028】そこで、請求項7記載の発明のように、伝
達特性制御手段が、伝達関数モデルと直列に配設された
伝達関数補正フィルタの周波数特性を、最適値収束判定
が最適値に収束したと判定した際のディジタルフィルタ
のフィルタ係数及び最適値に収束していないディジタル
フィルタのフィルタ係数に基づいて実際の収束経路が最
短の収束経路に近づくように変更する結果、最短の収束
経路と実際の収束経路との間のずれが、的確に適応処理
手段の処理内容の補正に反映する。
Therefore, the transfer characteristic control means converges the frequency characteristic of the transfer function correction filter arranged in series with the transfer function model to the optimum value by the optimum value convergence determination. The actual convergence path is determined based on the filter coefficient of the digital filter when it is determined that
As a result of the change so as to approach the short convergence path, the deviation between the shortest convergence path and the actual convergence path is accurately reflected in the correction of the processing content of the adaptive processing means.

【0029】同様に、請求項8記載の発明のように、伝
達特性制御手段が、基準信号に応じて制御音源を駆動す
る信号を生成するディジタルフィルタと直列に配設され
た伝達関数補正フィルタの周波数特性を実際の収束経路
が最短の収束経路に近づくように適宜変更する結果、最
短の収束経路と実際の収束経路との間のずれが、的確に
適応処理手段の処理内容の補正に反映する。また、請求
項9記載の発明にあっては、補正手段としてのゲイン調
整手段が、ディジタルフィルタのフィルタ係数の更新ゲ
インを変更する結果、上記請求項2記載の発明と同様
に、例えば、最短の収束経路に対する実際の収束経路の
ずれが大きい場合に更新ゲインを小さくすれば、収束経
路のずれの拡大が防がれる。
Similarly, the transfer characteristic control means of the transfer function correction filter arranged in series with the digital filter for generating the signal for driving the control sound source according to the reference signal. Frequency characteristics are actually converged
Is appropriately changed so as to approach the shortest convergent path, the deviation between the shortest convergent path and the actual convergent path is accurately reflected in the correction of the processing content of the adaptive processing means. In the invention according to claim 9, as a result of the gain adjusting means as the correcting means changing the update gain of the filter coefficient of the digital filter, as in the invention according to claim 2, for example, the shortest When the deviation of the actual convergence path with respect to the convergence path is large, the update gain can be reduced to prevent the deviation of the convergence path from increasing.

【0030】さらに、請求項10記載の発明にあって
は、上記請求項3記載の発明と同様に、適応処理手段の
処理内容が不明確な情報に基づいて変動してしまうよう
なことが避けられる。そして、請求項11記載の発明に
あっては、上記請求項4記載の発明と同様に、最適値収
束判定手段が最適値に収束したと判定した際のディジタ
ルフィルタのフィルタ係数及び制御開始時のディジタル
フィルタのフィルタ係数(つまり、フィルタ係数の初期
値)から決まる最短の収束経路の方向と、最適値に収束
していないディジタルフィルタのフィルタ係数(例え
ば、制御開始直後のフィルタ係数)及び制御開始時のデ
ィジタルフィルタのフィルタ係数から決まる実際の収束
経路の方向との間のずれが無くなるように、補正手段が
適応処理手段の処理内容を補正するから、最短の収束経
路と実際の収束経路との間のずれが無くなるような補正
が行われる。
Further, in the invention described in claim 10, as in the invention described in claim 3, it is avoided that the processing content of the adaptive processing means changes based on unclear information. To be In the invention described in claim 11, as in the invention described in claim 4, the filter coefficient of the digital filter when the optimum value convergence determining means determines that the optimum value has converged to the optimum value and the control coefficient at the start of control The direction of the shortest convergence path that is determined from the filter coefficient of the digital filter (that is, the initial value of the filter coefficient), the filter coefficient of the digital filter that has not converged to the optimum value (for example, the filter coefficient immediately after the start of control), and at the start of control. Since the correction means corrects the processing contents of the adaptive processing means so as to eliminate the deviation from the direction of the actual convergence path determined by the filter coefficient of the digital filter of, the shortest convergence path and the actual convergence path. The correction is performed so that the deviation of

【0031】[0031]

【実施例】以下、この発明の実施例を図面に基づいて説
明する。図1は、本発明の第1実施例の全体構成を示す
図であり、この実施例は、騒音源としてのエンジン4か
ら空間としての車室6内に伝達されるこもり音の低減を
図る能動型騒音制御装置1に本発明を適用したものであ
る。
Embodiments of the present invention will be described below with reference to the drawings. FIG. 1 is a diagram showing an overall configuration of a first embodiment of the present invention, which is an active system for reducing muffled noise transmitted from an engine 4 as a noise source into a vehicle interior 6 as a space. The present invention is applied to the mold noise control device 1.

【0032】先ず、構成を説明すると、車体3は、前輪
2a,2b,後輪2c,2d及び各車輪2a〜2dと車
体3との間に介在するサスペンションによって支持され
ている。なお、図1に示す車両は、前輪2a及び2bが
車体3前部に配置されたエンジン4によって回転駆動さ
れるいわゆる前置きエンジン前輪駆動車である。エンジ
ン4には、騒音発生状態検出手段としてのクランク角セ
ンサ5が取り付けられていて、このクランク角センサ5
は、エンジン4のクランク角に対応する(例えば、レシ
プロ4気筒の場合は、180度回転する度に一つの)パ
ルス信号でなる基準信号xをコントローラ10に供給す
る。
First, the structure will be described. The vehicle body 3 is supported by front wheels 2a, 2b, rear wheels 2c, 2d, and suspensions interposed between the wheels 2a to 2d and the vehicle body 3. The vehicle shown in FIG. 1 is a so-called front-mounted engine front-wheel drive vehicle in which the front wheels 2a and 2b are rotationally driven by an engine 4 arranged in the front part of the vehicle body 3. The engine 4 is provided with a crank angle sensor 5 as a noise generation state detecting means.
Supplies a reference signal x, which is a pulse signal corresponding to the crank angle of the engine 4 (for example, in the case of a reciprocating 4-cylinder, one pulse every 180 degrees) to the controller 10.

【0033】また、車体3内の車室6内には、制御アク
チュエータ又は制御音源としてのラウドスピーカ7a,
7b,7c及び7dが、前部座席S1 ,S2 及び後部座
席S 3 ,S4 のそれぞれに対向するドア部に配置されて
いる。さらに、各座席S1 〜S4 のヘッドレスト位置に
は、残留騒音検出手段としてのマイクロフォン8a〜8
hが、それぞれ二つずつ配設されていて、これらマイク
ロフォン8a〜8hが音圧として測定した残留騒音信号
1 〜e8 が、コントローラ10に供給される。
In the vehicle interior 6 of the vehicle body 3, a control actuator is installed.
A loudspeaker 7a as a chute or a control sound source,
7b, 7c and 7d are front seats S1, S2And rear seat
Seat S 3, SFourPlaced on the doors facing each of the
There is. Furthermore, each seat S1~ SFourAt the headrest position
Are microphones 8a to 8 as residual noise detecting means.
There are two h each, and these microphones
Residual noise signal measured as sound pressure by Rophone 8a-8h
e1~ E8Are supplied to the controller 10.

【0034】そして、コントローラ10は、クランク角
センサ5から供給される基準信号xと、マイクロフォン
8a〜8hから供給される残留騒音信号e1 〜e8 とに
基づいて、後述する演算処理を実行し、車室6内に伝達
されるこもり音を打ち消すような制御音がラウドスピー
カ7a〜7dから発せられるように、それらラウドスピ
ーカ7a〜7dに駆動信号y1 〜y4 を出力する。
Then, the controller 10 executes the arithmetic processing described later based on the reference signal x supplied from the crank angle sensor 5 and the residual noise signals e 1 to e 8 supplied from the microphones 8a to 8h. The drive signals y 1 to y 4 are output to the loudspeakers 7a to 7d so that the loudspeakers 7a to 7d generate a control sound for canceling the muffled sound transmitted to the vehicle interior 6.

【0035】コントローラ10は、図2に示すように、
クランク角センサ5から供給される基準信号xをディジ
タル値に変換するA/D変換器11と、このA/D変換
器11でディジタル値に変換された基準信号xが入力さ
れるディジタルフィルタ12と、同じく基準信号xが入
力されるフィルタ係数可変のディジタルフィルタとして
の適応ディジタルフィルタ13と、マイクロフォン8a
〜8hから供給され且つアンプ14a〜14hによって
増幅された残留騒音信号e1 〜e8 をディジタル値に変
換するA/D変換器15a〜15hと、ディジタルフィ
ルタ12でフィルタ処理された処理信号rlm及びA/D
変換器15a〜15hによってディジタル値に変換され
た残留騒音信号e1 〜e8 に基づいて所定の演算処理を
実行して適応ディジタルフィルタ13のフィルタ係数を
更新するマイクロプロセッサ16と、適応ディジタルフ
ィルタ13から出力された駆動信号y1 〜y4 をアナロ
グ値に変換してラウドスピーカ7a〜7dを駆動するア
ンプ18a〜18dに供給するD/A変換器17a〜1
7dと、を備えている。
The controller 10 is, as shown in FIG.
An A / D converter 11 for converting the reference signal x supplied from the crank angle sensor 5 into a digital value, and a digital filter 12 to which the reference signal x converted into a digital value by the A / D converter 11 is input. , An adaptive digital filter 13 as a filter filter variable digital filter to which the reference signal x is input, and a microphone 8a.
An A / D converter 15a~15h for converting the residual noise signal e 1 to e 8 amplified by the supplied and amplifier 14a~14h to a digital value from ~8H, processed signal r lm which is filtered by the digital filter 12 And A / D
A microprocessor 16 for executing a predetermined arithmetic processing based on the residual noise signals e 1 to e 8 converted into digital values by the converters 15a to 15h to update the filter coefficient of the adaptive digital filter 13, and the adaptive digital filter 13 output from the drive signal y 1 ~y 4 is converted into an analog value and supplies the amplifier 18a~18d driving the loudspeaker 7a to 7d D / a converter 17a~1
7d.

【0036】ここで、ディジタルフィルタ12は、ラウ
ドスピーカ7a〜7dとマイクロフォン8a〜8hとの
間の伝達関数を有限インパルス応答関数の形でモデル化
した伝達関数モデルとしてのディジタルフィルタC
lm(l=1,2,…,L、m=1,2,…,M)を、M
個(本実施例では、M=4)のラウドスピーカ7a〜7
d及びL個(本実施例では、L=8)のマイクロフォン
8a〜8hの全ての組み合わせ(L×M)について有し
ていて、基準信号xをそれらディジタルフィルタC lm
フィルタ処理した値rlmを生成し出力する。
Here, the digital filter 12 is
Of the speaker 7a-7d and the microphone 8a-8h
The transfer function between two models in the form of a finite impulse response function
Digital filter C as a transfer function model
lm(L = 1, 2, ..., L, m = 1, 2, ..., M)
(In this embodiment, M = 4) loudspeakers 7a to 7
d and L (L = 8 in this embodiment) microphones
8a-8h for all combinations (L × M)
And the reference signal x is converted to those digital filters C lmso
The filtered value rlmIs generated and output.

【0037】一方、適応ディジタルフィルタ13は、ラ
ウドスピーカ7a〜7dの個数に対応してM個のフィル
タ係数可変の適応ディジタルフィルタWm を有してい
て、基準信号xをそれら適応ディジタルフィルタWm
フィルタ処理することにより、駆動信号y1 〜y4 を生
成し出力する。そして、マイクロプロセッサ16は、デ
ィジタルフィルタ12から供給される処理信号rlmと、
マイクロフォン8a〜8hから供給される残留騒音信号
1 〜e8 とに応じて、適応ディジタルフィルタ13の
各適応ディジタルフィルタWmのフィルタ係数を、LM
Sアルゴリズムに基づいて更新する。
On the other hand, the adaptive digital filter 13 has M adaptive filter filters W m with variable filter coefficients corresponding to the number of loudspeakers 7a to 7d, and the reference signal x is adapted to the adaptive digital filter W m. The drive signals y 1 to y 4 are generated and output by performing the filtering process with. Then, the microprocessor 16 receives the processed signal r lm supplied from the digital filter 12,
According to the residual noise signals e 1 to e 8 supplied from the microphones 8a to 8h, the filter coefficient of each adaptive digital filter W m of the adaptive digital filter 13 is set to LM.
Update based on S algorithm.

【0038】ここで、LMSアルゴリズムは、適応ディ
ジタルフィルタのフィルタ係数を更新するのに好適なア
ルゴリズムの一つであって、l番目(l=1,2,…,
L:L=8)のマイクロフォンが検出した残留騒音信号
をel (n)、ラウドスピーカから制御音が発生してい
ない時のl番目のマイクロフォンが検出した残留騒音信
号をepl (n)、m番目(m=1,2,…,M:M=
4)のラウドスピーカとl番目のマイクロフォンとの間
の伝達関数を有限インパルス応答関数の形でモデル化し
たディジタルフィルタClmのj番目(j=0,1,2,
…,IC −1:IC は定数)のフィルタ係数をClmj
基準信号をx(n)、基準信号x(n)が入力されたm
番目のラウドスピーカを駆動する適応ディジタルフィル
タWm のi番目(i=0,1,2,…,IK −1:IK
は定数)のフィルタ係数をWmiとすると、 が成立する。
Here, the LMS algorithm is one of the algorithms suitable for updating the filter coefficient of the adaptive digital filter, and is the l-th (l = 1, 2, ...,).
L: L = 8) the residual noise signal detected by the microphone is e l (n), the residual noise signal detected by the l-th microphone when no control sound is generated from the loudspeaker is e l (n), m-th (m = 1, 2, ..., M: M =
J-th digital filter C lm of the transfer function between the loudspeaker and the l-th microphone modeled in the form of a finite impulse response function 4) (j = 0,1,2,
, I C −1: I C is a constant), and the filter coefficient is C lmj ,
X (n) is the reference signal, and m is the reference signal x (n)
The i-th (i = 0, 1, 2, ..., I K −1: I K of the adaptive digital filter W m for driving the th loudspeaker
Is a constant) and the filter coefficient is W mi , Is established.

【0039】なお、(n)がつく項は、いずれもサンプ
リング時刻nにおけるサンプル値を表し、また、IC
フィルタClmのタップ数(フィルタ次数)、IK は適応
ディジタルフィルタWm のタップ数(フィルタ次数)で
ある。上記(1)式中、右辺の「ΣWmix(n−j−
i)」の項は適応ディジタルフィルタに基準信号x
(n)を入力した時の出力ym (n)を表し、「ΣC
lmj {ΣWmix(n−j−i)}」の項はm番目のラウ
ドスピーカに入力された信号y m (n)がそこから制御
音として空間に出力され伝達関数Clmを経てl番目のマ
イクロフォンに到達した時の信号を表し、さらに、「Σ
ΣClmj {ΣWmix(n−j−i)}」の項はl番目の
マイクロフォンへ到達した信号を足し合わせているか
ら、l番目のマイクロフォンに到達する制御音の総和を
表している。
The terms with (n) are all sump
Represents the sampled value at ring time n, and ICIs
Filter ClmNumber of taps (filter order), IKIs adapted
Digital filter WmBy the number of taps (filter order)
is there. In the above formula (1), “ΣW” on the right sidemix (n-j-
i) ”is the reference signal x to the adaptive digital filter.
Output y when (n) is inputm(N) is represented by “ΣC
lmj{ΣWmix (n-j-i)} ”is the m-th Lau
Signal y input to the speaker m(N) controls from there
Transfer function C is output to the space as soundlmThrough the l-th
It shows the signal when it reaches the icrophone, and
ΣClmj{ΣWmix (n−j−i)} ”is the l-th term
Is the signal reaching the microphone added up?
The sum of control sounds reaching the l-th microphone
It represents.

【0040】次いで、評価関数Jeを、 とする。Next, the evaluation function Je is And

【0041】そして、評価関数Jeを最小にするフィル
タ係数Wmiを求めるのがLMSアルゴリズムであり、具
体的には、評価関数Jeを各フィルタ係数Wmiについて
偏微分した値で、フィルタ係数Wmiを更新する。そこ
で、上記(2)式より、 となるが、上記(1)式より、 となるから、この(4)式の右辺をrlm(n−i)とお
けば、フィルタ係数の更新は、下記の(5)式のように
なる。
The LMS algorithm finds the filter coefficient W mi that minimizes the evaluation function Je. Specifically, the filter coefficient W mi is a value obtained by partially differentiating the evaluation function Je with respect to each filter coefficient W mi. To update. Therefore, from equation (2) above, From equation (1) above, Therefore, if the right side of the equation (4) is set to r lm (n−i), the filter coefficient is updated by the following equation (5).

【0042】 ここで、更新ゲインとしての収束係数αは、フィルタが
最適に収束する速度や、その安定性に関与する係数であ
る。
[0042] Here, the convergence coefficient α as the update gain is a coefficient relating to the speed at which the filter converges optimally and its stability.

【0043】つまり、マイクロプロセッサ16は、上記
(5)式に基づいて、適応ディジタルフィルタWm のフ
ィルタ係数Wmiを更新する。さらに、本実施例では、デ
ィジタルフィルタ12を構成する各ディジタルフィルタ
lmを、LMSアルゴリズムの演算に必要な次数IC
りも多い次数とし、且つ、それらフィルタ係数Clmj
マイクロプロセッサ16によって変更可能とし、もっ
て、各ディジタルフィルタClmの進み・遅れを任意に調
整できるようにしている。
That is, the microprocessor 16 updates the filter coefficient W mi of the adaptive digital filter W m based on the above equation (5). Further, in this embodiment, each digital filter C lm constituting the digital filter 12 is set to an order larger than the order I C necessary for the operation of the LMS algorithm, and the filter coefficient C lmj can be changed by the microprocessor 16. Therefore, the lead / lag of each digital filter C lm can be arbitrarily adjusted.

【0044】ここで、エンジン4で発生する騒音が単一
周波数であって、適応ディジタルフィルタWm のフィル
タ次数が2次(IK =2)であると仮定し、この場合の
適応ディジタルフィルタWm のフィルタ係数Wm0,Wm1
の最適値への収束の様子を、図3を伴って説明する。な
お、図3は制御開始時のフィルタ係数を原点とした直交
座標系であって、最適値に対応する座標を(Wm00 ,W
m10 )としている。
Here, it is assumed that the noise generated in the engine 4 has a single frequency and the filter order of the adaptive digital filter W m is the second order (I K = 2), and the adaptive digital filter W in this case is filter coefficients of m W m0, W m1
The state of convergence to the optimum value of is described with reference to FIG. Note that FIG. 3 is an orthogonal coordinate system in which the filter coefficient at the start of control is the origin, and the coordinates corresponding to the optimum value are (W m00 , W
m10 ).

【0045】即ち、各ディジタルフィルタClmがラウド
スピーカ7a〜7d及びマイクロフォン8a〜8h間の
音響伝達特性を正確に表しているのならば、原点から始
まるフィルタ係数Wm0,Wm1の収束経路は、Aで示すよ
うに一直線に最適点(Wm00,Wm10 )へ向かう最短経
路を採るはずである。しかし、車室6内の温度や湿度等
の状況の変化や、乗員の増減、或いはラウドスピーカ7
a〜7d,マイクロフォン8a〜8hの劣化等に伴い、
ラウドスピーカ7a〜7d及びマイクロフォン8a〜8
h間の音響伝達特性が変動し、各ディジタルフィルタC
lmの精度が相対的に低下すると、そのディジタルフィル
タClmが表す伝達関数が、実際の音響空間の伝達関数に
対して位相誤差を持つようになり、最適点(Wm00 ,W
m10 )とは別の誤った最適点(Wm00 ’,Wm10 ’)の
影響を受けてしまい、Bで示すように最適な収束経路A
から大きく外れた曲線の軌跡を描きながら最適点(W
m00 ,Wm10 )に収束するようになり、最適値への収束
はなされても、それに至るまでに比較的長時間を要する
ことになる。
That is, if each digital filter C lm accurately represents the acoustic transfer characteristic between the loudspeakers 7a to 7d and the microphones 8a to 8h, the convergence path of the filter coefficients W m0 and W m1 starting from the origin is , A, the shortest route that goes straight to the optimum point (W m00 , W m10 ) should be taken. However, changes in conditions such as temperature and humidity in the passenger compartment 6, increase / decrease in occupants, or the loudspeaker 7
a to 7d and microphones 8a to 8h are deteriorated,
Loudspeakers 7a-7d and microphones 8a-8
The acoustic transfer characteristic between h varies and each digital filter C
When the accuracy of lm relatively decreases, the transfer function represented by the digital filter C lm has a phase error with respect to the transfer function of the actual acoustic space, and the optimum point (W m00 , W
It is affected by the wrong optimum points (W m00 ', W m10 ') different from m10 ), and as shown by B, the optimum convergence path A
While drawing a curved trajectory that deviates greatly from the optimum point (W
m00 , W m10 ), and even if the convergence is made to the optimum value, it takes a relatively long time to reach it.

【0046】そこで、本実施例では、マイクロプロセッ
サ16において後述する演算処理を実行し、ディジタル
フィルタClmが表す伝達関数と、実際の音響空間の伝達
関数との間の位相誤差を求め、その位相誤差がなくなる
ように、ディジタルフィルタClmのフィルタ係数Clmj
を変更するものとする。図4は、コントローラ10内で
実行される処理の概要を示すフローチャートであり、以
下、図4に従って本実施例の動作を説明する。
Therefore, in the present embodiment, the microprocessor 16 executes the arithmetic processing described later to obtain the phase error between the transfer function represented by the digital filter C lm and the actual transfer function in the acoustic space, and the phase error is obtained. The filter coefficient C lmj of the digital filter C lm is set so that the error is eliminated.
Shall be changed. FIG. 4 is a flow chart showing the outline of the processing executed in the controller 10. The operation of this embodiment will be described below with reference to FIG.

【0047】先ず、ステップ001において、フィルタ
係数Wmiが最適値に収束しているか否かを判定する。な
お、このステップ001では、例えば、上記(5)式に
基づくフィルタ係数Wmiの更新量が、所定回数以上、所
定値を超えない場合には、フィルタ係数Wmiが最適値に
収束していると判定する。ステップ001でフィルタ係
数Wmiが未だ最適値に収束していないと判定された場合
には、ステップ002に移行し、上記(5)式に基づい
て、フィルタ係数Wmiを更新し、次いで、ステップ00
3に移行し、駆動信号y1 〜y4 を生成しラウドスピー
カ7a〜7dに出力する。
First, in step 001, it is determined whether or not the filter coefficient W mi has converged to the optimum value. In this step 001, for example, when the update amount of the filter coefficient W mi based on the above equation (5) does not exceed the predetermined value a predetermined number of times or more, the filter coefficient W mi converges to the optimum value. To determine. When it is determined in step 001 that the filter coefficient W mi has not yet converged to the optimum value, the process proceeds to step 002, the filter coefficient W mi is updated based on the equation (5), and then step 00
3 proceeds to generate a drive signal y 1 ~y 4 outputs to the loudspeaker 7a to 7d.

【0048】すると、ラウドスピーカ7a〜7dから車
室6内に制御音が発生するが、制御開始直後は適応ディ
ジタルフィルタWm の各フィルタ係数Wmiが最適値に収
束しているとは限らないので、必ずしも車室6内に伝達
されたこもり音が低減されるとはいえない。しかし、上
記ステップ001〜003の処理を繰り返し実行する結
果、各フィルタ係数Wmiは最適値に向かって収束してい
き、車室6内に伝達されるこもり音がラウドスピーカ7
a〜7dから発せられる制御音によって打ち消され、車
室6内の騒音の低減が図られる。
[0048] Then, the control sound into the passenger compartment 6 from loudspeakers 7a~7d occurs immediately after start of control is not necessarily the filter coefficient W mi of the adaptive digital filter W m is converged to an optimum value Therefore, it cannot be said that the muffled sound transmitted to the vehicle interior 6 is necessarily reduced. However, as a result of repeatedly executing the processing of steps 001 to 003, each filter coefficient W mi converges toward the optimum value, and the muffled sound transmitted into the vehicle interior 6 is transmitted to the loudspeaker 7.
It is canceled by the control sounds emitted from a to 7d, and the noise in the vehicle interior 6 is reduced.

【0049】そして、フィルタ係数Wmiが最適値に収束
すると、上記ステップ001の判定が「YES」となる
から、ステップ004に移行し、その時のフィルタ係
数、即ち、最適値(Wm00 ,Wm10 )を記憶する。な
お、説明の便宜上、適応ディジタルフィルタWm のフィ
ルタ次数は2次とする。次いで、ステップ005に移行
し、再び初期値から適応ディジタルフィルタW m のフィ
ルタ係数Wmiの更新を所定の回数行い、制御開始直後に
おける適応ディジタルフィルタWm のフィルタ係数(W
m01 ,Wm11 )を求める。
Then, the filter coefficient WmiConverges to the optimal value
Then, the determination in step 001 is “YES”.
To step 004, the filter section at that time
Number, that is, the optimum value (Wm00, Wm10) Is remembered. Na
For convenience of explanation, the adaptive digital filter WmThe fi
Ruta order is quadratic. Then move to step 005
And again from the initial value, the adaptive digital filter W mThe fi
Coefficient WmiIs updated a predetermined number of times, and immediately after the start of control
Adaptive digital filter W inmFilter coefficient (W
m01, Wm11).

【0050】このフィルタ係数(Wm01 ,Wm11 )は、
図3に示す実際の収束経路Bの原点近傍の点であるか
ら、原点からこのフィルタ係数(Wm01 ,Wm11 )に向
かう方向は、原点から誤った最適点(Wm00 ’,W
m10 ’)に向かう方向と一致するとしても大きな差異は
ない。そして、原点から最適点(Wm00 ,Wm10 )に向
かう方向と、原点から誤った最適点(Wm00 ’,W
m10 ’)に向かう方向との間のずれは、結局はディジタ
ルフィルタClmが表す伝達関数と、実際の音響空間の伝
達関数との間の位相誤差に起因するものであるから、そ
れら方向のずれに基づいてディジタルフィルタClmを補
正すれば、伝達関数間の位相誤差を低減することができ
る。
This filter coefficient (W m01 , W m11 ) is
Since the point is near the origin of the actual convergence path B shown in FIG. 3, the direction from the origin to the filter coefficient (W m01 , W m11 ) is incorrect from the origin to the optimum point (W m00 ', W m00 ).
There is no big difference even if it agrees with the direction toward m10 '). The direction from the origin to the optimum point (W m00 , W m10 ) and the wrong optimum point (W m00 ', W m from the origin)
The deviation from the direction toward m10 ') is due to the phase error between the transfer function represented by the digital filter C lm and the actual transfer function in the acoustic space, and therefore the deviation from those directions. If the digital filter C lm is corrected based on, the phase error between transfer functions can be reduced.

【0051】そこで、ステップ006に移行し、下記の
(6)式に基づいて、原点から最適点(Wm00
m10 )に向かう方向と、原点から制御開始直後のフィ
ルタ係数(Wm01 ,Wm11 )に向かう方向との間の角度
Δθを演算する。ただし、原点(制御開始時)のフィル
タ係数は0としている。 Δθ=cos-1{(Wm00 m01 +Wm10 m11 )/((Wm00 2+Wm10 21/ 2 (Wm01 2+Wm11 21/2 )} ……(6) 次いで、ステップ007に移行し、ディジタルフィルタ
lmの位相が上記ステップ006で求めた角度Δθに応
じた分だけ進む又は遅れるように、そのフィルタ係数C
lmj を適宜調整する。
Therefore, the process proceeds to step 006, and the optimum point (W m00 , W m00 ,
The angle Δθ between the direction toward W m10 ) and the direction from the origin toward the filter coefficient (W m01 , W m11 ) immediately after the start of control is calculated. However, the filter coefficient at the origin (at the start of control) is 0. Δθ = cos -1 {(W m00 W m01 + W m10 W m11) / ((W m00 2 + W m10 2) 1/2 (W m01 2 + W m11 2) 1/2)} ...... (6) Then, step 007, the filter coefficient C of the digital filter C lm is advanced or delayed by an amount corresponding to the angle Δθ obtained in step 006.
Adjust lmj accordingly.

【0052】即ち、図5に示すように、任意の周波数に
おける位相誤差(∝角度Δθ)が判れば、全周波数域の
位相誤差は推定できるから、その位相誤差がなくなるよ
うにディジタルフィルタClmの位相を調整すればよい。
このようにディジタルフィルタClmを補正した結果、デ
ィジタルフィルタClmが表す伝達関数と実際の音響空間
の伝達関数との間の位相誤差が小さくなるから、その後
に例えば騒音の周波数の変動に伴って最適値が変化して
再び適応ディジタルフィルタWm のフィルタ係数Wmi
更新する際には、そのフィルタ係数Wmiは、図3のAに
示すような最短の収束経路の近傍を通過して最適値に収
束するようになり、素早い適応動作が可能となる。
That is, as shown in FIG. 5, if the phase error (∝ angle Δθ) at an arbitrary frequency is known, the phase error in the entire frequency range can be estimated, so that the digital filter C lm can eliminate the phase error. Adjust the phase.
As a result of correcting the digital filter C lm in this way, the phase error between the transfer function represented by the digital filter C lm and the transfer function of the actual acoustic space becomes small, and thereafter, for example, with the fluctuation of the frequency of noise. When the optimum value changes and the filter coefficient W mi of the adaptive digital filter W m is updated again, the filter coefficient W mi passes through the vicinity of the shortest convergence path as shown in A of FIG. The value converges to the value, and quick adaptive operation becomes possible.

【0053】そして、適応ディジタルフィルタWm の適
応動作が遅くなる主な原因は、実空間におけるラウドス
ピーカ7a〜7dとマイクロフォン8a〜8hとの間の
音響伝達特性が変化した結果、ディジタルフィルタ12
の精度が相対的に低下してしまうことであるから、本実
施例のように、ディジタルフィルタClmのフィルタ係数
lmj を適宜変更する構成であれば、音響伝達特性の変
化を的確に制御内容に反映することができる。
The main cause of the delay in the adaptive operation of the adaptive digital filter W m is that the digital filter 12 is changed as a result of the change in the acoustic transfer characteristics between the loudspeakers 7a to 7d and the microphones 8a to 8h in the real space.
Therefore, if the configuration is such that the filter coefficient C lmj of the digital filter C lm is appropriately changed as in the present embodiment, the change of the acoustic transfer characteristics can be controlled accurately. Can be reflected in.

【0054】しかも、予め音響伝達特性の変化を考慮し
て複数の伝達関数を求めてこれを記憶する必要がないか
ら、特に記憶容量の増大を招くことがなく、実用的であ
るし、また、同定音を発生させる必要もないから、乗員
に不快感を与えるおそれもない。なお、ステップ005
の処理を実行する際にはステップ002,003の通常
の騒音低減制御を停止することになるが、ステップ00
4〜007の処理を、例えば、ステップ002,003
の処理を100ループ行うのに対して1〜2回程度実行
するという制御内容にすることにより、騒音低減制御の
劣化を極めて低く抑えることができるから、特に不具合
はない。
Moreover, since it is not necessary to obtain a plurality of transfer functions in advance in consideration of changes in the acoustic transfer characteristics and store the transfer functions, the storage capacity is not particularly increased, which is practical. Since it is not necessary to generate the identification sound, the occupant is not uncomfortable. Note that step 005
The normal noise reduction control of steps 002 and 003 is stopped when the processing of step 00 is executed.
The processing of 4 to 007 is performed, for example, in steps 002 and 003.
There is no particular problem because the deterioration of the noise reduction control can be suppressed to an extremely low level by setting the control content to be executed once or twice for 100 loops.

【0055】また、本実施例では、ディジタルフィルタ
lmのフィルタ次数をLMSアルゴリズムの演算に必要
な次数IC よりも多い次数とすることにより各ディジタ
ルフィルタClmの進み・遅れを任意に調整できるように
しているが、例えば、当初設定された伝達関数にexp
(Δθ(ω))(ωは位相誤差が生じた周波数)を畳み
込むことにより伝達関数を演算し直し、その演算結果で
ディジタルフィルタC lmを置き換えることにより進み・
遅れを任意に調整できるようにしてもよい。
In this embodiment, the digital filter is used.
ClmRequired filter order for LMS algorithm operation
Na degree ICBy setting each digit to a higher degree than
Le Filter ClmTo be able to arbitrarily adjust the lead / lag of
However, for example, the transfer function initially set is exp
Convolve (Δθ (ω)) (ω is the frequency where the phase error occurs)
The transfer function is recalculated by inserting the
Digital filter C lmProceed by replacing
The delay may be adjusted arbitrarily.

【0056】ここで、本実施例では、ディジタルフィル
タ12,マイクロプロセッサ16及びステップ002,
ステップ003の処理によって適応処理手段が構成さ
れ、ステップ001の処理が最適値収束判定手段に対応
し、ステップ004〜ステップ007の処理によって補
正手段が構成される。図6乃至図8は、本発明の第2実
施例を示す図であり、上記第1実施例と同様にエンジン
から車室内に伝達されるこもり音の低減を図る能動型騒
音制御装置に本発明を適用したものである。
Here, in this embodiment, the digital filter 12, the microprocessor 16 and the step 002 are used.
The processing of step 003 constitutes the adaptive processing means, the processing of step 001 corresponds to the optimum value convergence determination means, and the processing of steps 004 to 007 constitutes the correction means. FIGS. 6 to 8 are views showing a second embodiment of the present invention. The present invention relates to an active noise control device for reducing muffled noise transmitted from an engine into a vehicle compartment as in the first embodiment. Is applied.

【0057】即ち、本実施例では、図6に示すように、
ディジタルフィルタ12とマイクロプロセッサ16との
間に、図7に示すような周波数特性を有し且つその周波
数特性が可変の伝達関数補正フィルタ20をディジタル
フィルタ12と直列の関係になるように介在させ、その
伝達関数補正フィルタ20の周波数特性をマイクロプロ
セッサ16により制御可能としたものである。
That is, in this embodiment, as shown in FIG.
A transfer function correction filter 20 having a frequency characteristic as shown in FIG. 7 and a variable frequency characteristic is interposed between the digital filter 12 and the microprocessor 16 in a serial relationship with the digital filter 12. The frequency characteristic of the transfer function correction filter 20 can be controlled by the microprocessor 16.

【0058】この伝達関数補正フィルタ20は、図7
(a)に示すように、所定周波数fa〜fb 間の位相を
任意の値に制御することができ、且つ、図7(b)に示
すように振幅は変化させない特性を有するフィルタであ
る。つまり、ディジタルフィルタ12から出力された処
理信号rlmをこの伝達関数補正フィルタ20を通じてマ
イクロプロセッサ16に供給することにより、周波数f
a 〜fb 間の処理信号rlmの位相のみが変化することに
なるから、結果として、ディジタルフィルタ12の処理
内容が補正されたことになる。
This transfer function correction filter 20 is shown in FIG.
(A), the phase between the predetermined frequency f a ~f b can be controlled to an arbitrary value, and is a filter having the characteristic of not amplitude was varied as shown in FIG. 7 (b) . That is, by supplying the processed signal r lm output from the digital filter 12 to the microprocessor 16 through the transfer function correction filter 20, the frequency f
Only the phase of the processed signal r lm between a and f b changes, and as a result, the processing content of the digital filter 12 is corrected.

【0059】図8は、コントローラ10における処理の
概要を表すフローチャートであり、以下、図8に従って
本実施例の動作を説明する。なお、図8に示すステップ
001〜006の処理は、上記第1実施例と同様である
ため、その重複する説明は省略する。即ち、ステップ0
06で角度Δθが求められたら、ステップ010に移行
し、ディジタルフィルタClmから出力された処理信号r
lmの位相がその角度Δθに応じた分だけ進む又は遅れる
ように、伝達関数補正フィルタのフィルタ係数を変更す
る。なお、ステップ010では、騒音の周波数に応じて
図7(a)に示す周波数fa ,fb をも変更する。
FIG. 8 is a flow chart showing an outline of processing in the controller 10. The operation of this embodiment will be described below with reference to FIG. Since the processing of steps 001 to 006 shown in FIG. 8 is the same as that of the first embodiment, the duplicated description will be omitted. That is, step 0
When the angle Δθ is obtained in 06, the process proceeds to step 010, and the processed signal r output from the digital filter C lm is output.
The filter coefficient of the transfer function correction filter is changed so that the phase of lm leads or lags by an amount corresponding to the angle Δθ. In step 010, the frequencies f a and f b shown in FIG. 7A are also changed according to the noise frequency.

【0060】この結果、ディジタルフィルタClm及び伝
達関数補正フィルタ20が表す伝達関数と実際の音響空
間の伝達関数との間の位相誤差が小さくなるから、上記
第1実施例と同様の作用効果が得られる。ここで、本実
施例では、マイクロプロセッサ16及びステップ007
の処理によって伝達特性制御手段が構成される。
As a result, the phase error between the transfer function represented by the digital filter C lm and the transfer function correction filter 20 and the transfer function of the actual acoustic space becomes small, so that the same effect as the first embodiment can be obtained. can get. Here, in this embodiment, the microprocessor 16 and step 007 are used.
The transfer characteristic control means is constituted by the processing of.

【0061】図9は、本発明の第3実施例を示す図であ
り、本実施例も、上記第1実施例と同様にエンジンから
車室内に伝達されるこもり音の低減を図る能動型騒音制
御装置に本発明を適用したものである。即ち、本実施例
では、伝達関数補正フィルタ20を、ディジタルフィル
タ12の次段ではなく、適応ディジタルフィルタ13と
直列の関係に配設したことを除いては、上記第2実施例
と同様である。
FIG. 9 is a diagram showing a third embodiment of the present invention. In this embodiment as well, similar to the first embodiment, the active noise for reducing the muffled noise transmitted from the engine into the vehicle compartment is used. The present invention is applied to a control device. That is, this embodiment is the same as the above-mentioned second embodiment except that the transfer function correction filter 20 is arranged in series with the adaptive digital filter 13 instead of the next stage of the digital filter 12. .

【0062】このような構成でも、駆動信号y1 〜y4
がディジタルフィルタClmが表す伝達関数と、実際の音
響空間の伝達関数との間の位相誤差に応じて補正される
から、上記第1実施例と同様の作用効果が得られる。図
10は本発明の第4実施例を示す図であり、本実施例
も、上記第1実施例と同様にエンジンから車室内に伝達
されるこもり音の低減を図る能動型騒音制御装置に本発
明を適用したものである。
Even with such a configuration, the drive signals y 1 to y 4
Is corrected according to the phase error between the transfer function represented by the digital filter C lm and the transfer function of the actual acoustic space, so that the same effect as the first embodiment can be obtained. FIG. 10 is a diagram showing a fourth embodiment of the present invention. This embodiment is also applied to an active noise control device for reducing muffled noise transmitted from the engine into the vehicle compartment as in the first embodiment. The invention is applied.

【0063】図10は、コントローラ10における処理
の概要を表すフローチャートであり、ステップ001〜
006の処理は、上記第1実施例と同様であるため、そ
の重複する説明は省略する。即ち、本実施例では、ステ
ップ006で角度Δθが求められたら、ステップ011
に移行し、角度Δθの大きさに応じて収束係数αを補正
するものである。
FIG. 10 is a flow chart showing the outline of the processing in the controller 10.
Since the processing of 006 is the same as that of the first embodiment, the duplicated description will be omitted. That is, in this embodiment, when the angle Δθ is obtained in step 006, step 011
Then, the convergence coefficient α is corrected according to the magnitude of the angle Δθ.

【0064】例えば、角度Δθが大きいということは、
ディジタルフィルタClmが表す伝達関数と、実際の音響
空間の伝達関数との間の位相誤差が大きいということで
あるから、収束係数αが大きいと図3に示す経路Bの膨
らみがさらに大きくなって、最適値への収束にさらに時
間を要することになるから、収束係数αを小さくするこ
とが望ましい。
For example, a large angle Δθ means that
Since the phase error between the transfer function represented by the digital filter C lm and the transfer function of the actual acoustic space is large, if the convergence coefficient α is large, the bulge of the path B shown in FIG. 3 is further increased. Since it takes more time to converge to the optimum value, it is desirable to reduce the convergence coefficient α.

【0065】逆に、角度Δθが小さいということは、デ
ィジタルフィルタClmが表す伝達関数と、実際の音響空
間の伝達関数との間の位相誤差が小さいということであ
るから、収束係数αを大きくすれば、最適値への適応動
作がさらに短時間で行えるようになる。ここで、本実施
例ではステップ011の処理がゲイン調整手段に対応す
る。
On the contrary, the small angle Δθ means that the phase error between the transfer function represented by the digital filter C lm and the transfer function in the actual acoustic space is small, so that the convergence coefficient α is increased. Then, the adaptive operation to the optimum value can be performed in a shorter time. Here, in this embodiment, the processing of step 011 corresponds to the gain adjusting means.

【0066】なお、特に図示はしないが、収束係数αの
補正と、上記第1乃至第3実施例で説明したディジタル
フィルタClmの補正とを両方実行すれば、両方の効果を
得ることができる。図11は本発明の第5実施例を示す
図であり、本実施例の構成は、ステップ004の前にス
テップ012の処理を行うことを除いては、上記第1実
施例と同様である。
Although not shown in particular, both effects can be obtained by executing both the correction of the convergence coefficient α and the correction of the digital filter C lm described in the first to third embodiments. . FIG. 11 is a diagram showing a fifth embodiment of the present invention, and the configuration of the present embodiment is the same as that of the above-mentioned first embodiment except that the processing of step 012 is performed before step 004.

【0067】即ち、本実施例では、ステップ001で
「YES」と判定されたら、ステップ012に移行し、
適応ディジタルフィルタWm のフィルタ係数Wmiの最適
値が、変動しているか否かを判定し、その最適値が変動
していないと判定された場合にのみ、ステップ004以
降の処理を実行する。ステップ004以降の処理は、通
常の騒音低減制御に対して例えば100ループ毎の割り
込み処理として実行し、ディジタルフィルタClmを補正
するものであるが、上記(6)式で示したように、フィ
ルタ係数Wmiの最適値に基づいて位相誤差を求めている
ため、そのフィルタ係数Wmiの最適値が変動する状況、
例えば、エアコンディショナを作動させた直後のように
車室6内の温度が急峻に変化するような状況の下でディ
ジタルフィルタClmの補正を行っても、良好な効果が期
待できない。
That is, in this embodiment, if it is determined "YES" in step 001, the process proceeds to step 012,
It is determined whether or not the optimum value of the filter coefficient W mi of the adaptive digital filter W m is fluctuating, and only when it is determined that the optimum value is not fluctuating, the processing from step 004 is executed. The process after step 004 is executed as an interrupt process for every 100 loops with respect to the normal noise reduction control to correct the digital filter C lm , but as shown in the above formula (6), Since the phase error is obtained based on the optimum value of the coefficient W mi, the optimum value of the filter coefficient W mi changes,
For example, even if the digital filter C lm is corrected in a situation where the temperature inside the vehicle compartment 6 changes sharply immediately after the air conditioner is operated, a good effect cannot be expected.

【0068】そこで、ステップ012で最適値が変動し
ていると判定された場合には、ステップ002側に戻っ
て通常の騒音低減制御を実行することとしたのである。
その他の作用効果は、上記第1実施例と同様である。な
お、上記各実施例では、フィルタ係数Wmiが最適値に収
束したと判定された後に、ステップ005で初期値から
再びフィルタ係数Wmiの更新を数回行って制御開始直後
のフィルタ係数(Wm01 ,Wm11 )を求めているが、状
況変化がさほど激しくない場合には、ステップ004以
降の処理を実行する直前の騒音低減制御の制御開始直後
のフィルタ係数Wmiを記憶しておき、そのフィルタ係数
miに基づいて角度Δθを演算するようにしてもよく、
このようにすれば、ステップ005の処理が不要となっ
て、通常の騒音低減制御と並行してディジタルフィルタ
12の補正処理が行えるという利点がある。
Therefore, if it is determined in step 012 that the optimum value fluctuates, the process returns to step 002 and the normal noise reduction control is executed.
Other functions and effects are similar to those of the first embodiment. In each of the above embodiments, after it is determined that the filter coefficient W mi has converged to the optimum value, the filter coefficient W mi is updated again from the initial value in step 005 several times, and the filter coefficient (W m01 , W m11 ), but when the situation change is not so drastic, the filter coefficient W mi immediately after the start of the noise reduction control immediately before the processing of step 004 and thereafter is stored, and The angle Δθ may be calculated based on the filter coefficient W mi ,
By doing so, there is an advantage that the processing of step 005 becomes unnecessary and the correction processing of the digital filter 12 can be performed in parallel with the normal noise reduction control.

【0069】また、上記各実施例では、本発明にかかる
能動型騒音制御装置を車室6内に伝達される騒音の低減
を図る装置に適用した場合について説明したが、これ以
外の空間の騒音の低減を図る装置であってもよい。ま
た、上記各実施例では、上記(6)式で求められる角度
Δθに基づいて補正処理を行う構成としているが、これ
に限定されるものではなく、例えば、図3に示す経路B
の曲率や、経路A及びBで囲まれる面積等に基づいて補
正処理を行う構成としてもよい。
In each of the above embodiments, the case where the active noise control device according to the present invention is applied to the device for reducing the noise transmitted to the vehicle interior 6 has been described. May be a device for reducing the above. Further, in each of the above embodiments, the correction process is performed based on the angle Δθ obtained by the above equation (6), but the present invention is not limited to this, and for example, the route B shown in FIG.
The correction process may be performed on the basis of the curvature, the area surrounded by the paths A and B, and the like.

【0070】そして、上記各実施例では、エンジン4か
ら車室6内に伝達されるこもり音の低減を図る場合につ
いて説明しているが、本発明によって低減できる騒音は
これに限定されるものではなく、例えば、ロード・ノイ
ズに相関のあるサスペンション振動のピックアップ信
号、ドアミラー付近における風切り音のピックアップ信
号、ディファレンシャルギアやトランスミッションギア
のケース振動に対するピックアップ信号(駆動力伝達系
のケース振動に起因した騒音に相関のある信号)、車速
計測用のトランスミッションの出力軸の回転に応じたパ
ルス信号(ディファレンシャルギアやトランスミッショ
ンの噛み合いによる騒音に相関のある信号)等を検出す
る構成とすれば、それら各騒音の低減を図る装置とする
ことも可能である。
In each of the above-described embodiments, the case where the muffled noise transmitted from the engine 4 into the vehicle interior 6 is reduced has been described, but the noise that can be reduced by the present invention is not limited to this. Instead, for example, pickup signals for suspension vibration that correlates with road noise, pickup signals for wind noise near door mirrors, pickup signals for case vibration of differential gears and transmission gears (such as noise caused by case vibration of the driving force transmission system). Correlation signals) and pulse signals (signals that correlate with noise due to meshing of differential gears and transmission) according to rotation of the output shaft of the vehicle speed measurement transmission can be reduced It is also possible to use a device for

【0071】また、上記実施例では、本発明に係る適応
制御装置を車両用の能動型騒音制御装置に適用した場合
について説明しているが、本発明の適用対象は騒音の低
減を図る装置に限定されるものではなく、それ以外の適
応制御系であっても適用可能である。そして、上記実施
例では、制御音源としての四つのラウドスピーカ7a〜
7dと、残留騒音検出手段として八つのマイクロフォン
8a〜8hを設けているが、これらの数は任意である。
Further, in the above embodiment, the case where the adaptive control device according to the present invention is applied to the active noise control device for a vehicle is explained, but the application target of the present invention is to a device for reducing noise. The present invention is not limited to this, and other adaptive control systems can be applied. Then, in the above embodiment, four loudspeakers 7a to
7d and eight microphones 8a to 8h as residual noise detecting means are provided, but the number of these is arbitrary.

【0072】さらに、上記実施例では、適応ディジタル
フィルタWm のフィルタ係数Wmiを更新するアルゴリズ
ムとして時間領域のLMSアルゴリズムを適用している
が、これに限定されるものではなく、例えば周波数領域
のLMSアルゴリズム等を適用してもよい。
Further, in the above embodiment, the LMS algorithm in the time domain is applied as the algorithm for updating the filter coefficient W mi of the adaptive digital filter W m , but the present invention is not limited to this and, for example, in the frequency domain. LMS algorithm etc. may be applied.

【0073】[0073]

【発明の効果】以上説明したように、本発明によれば、
最適値に収束した際のディジタルフィルタのフィルタ係
から判る最適値への最短の収束経路と、最適値に収束
する前のフィルタ係数から判る実際の収束経路とに基づ
いて、適応処理手段の処理内容を実際の収束経路が最短
の収束経路に近づくように補正する構成としたため、必
要な記憶容量の増大や人間に不快感を与えるという不具
合を招くことなく、制御対象の状況変化や装置の構成部
品の劣化等に基づいて適応処理部の処理内容が実際の制
御対象や装置の特性に合致しない状態が生じても、ディ
ジタルフィルタのフィルタ係数を最適値に素早く収束さ
せることができるという効果がある。
As described above, according to the present invention,
The processing contents of the adaptive processing means based on the shortest convergence path to the optimum value found from the filter coefficient of the digital filter when converged to the optimum value and the actual convergence path found from the filter coefficient before convergence to the optimum value. The actual convergence path is the shortest
Since the configuration is such that the correction is performed so as to approach the convergence path of the above, it is adapted based on the change in the condition of the control target and the deterioration of the component parts of the device without causing the trouble of increasing the necessary storage capacity and causing discomfort to humans. Even if the processing content of the processing unit does not match the actual control target or the characteristics of the device, the filter coefficient of the digital filter can be quickly converged to the optimum value.

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

【図1】第1実施例の全体構成を示す図である。FIG. 1 is a diagram showing an overall configuration of a first embodiment.

【図2】第1実施例におけるコントローラの構成を示す
ブロック図である。
FIG. 2 is a block diagram showing a configuration of a controller in the first embodiment.

【図3】適応ディジタルフィルタのフィルタ係数の適応
の様子を示す図である。
FIG. 3 is a diagram showing how filter coefficients of an adaptive digital filter are adapted.

【図4】第1実施例における処理の概要を示すフローチ
ャートである。
FIG. 4 is a flowchart showing an outline of processing in the first embodiment.

【図5】位相誤差と角度Δθとの関係を説明する図であ
る。
FIG. 5 is a diagram illustrating a relationship between a phase error and an angle Δθ.

【図6】第2実施例の構成を示すブロック図である。FIG. 6 is a block diagram showing a configuration of a second embodiment.

【図7】伝達関数補正フィルタの周波数特性を示す図で
ある。
FIG. 7 is a diagram showing frequency characteristics of a transfer function correction filter.

【図8】第2実施例における処理の概要を示すフローチ
ャートである。
FIG. 8 is a flowchart showing an outline of processing in the second embodiment.

【図9】第3実施例の構成を示すブロック図である。FIG. 9 is a block diagram showing a configuration of a third embodiment.

【図10】第4実施例における処理の概要を示すフロー
チャートである。
FIG. 10 is a flowchart showing an outline of processing in a fourth embodiment.

【図11】第5実施例における処理の概要を示すフロー
チャートである。
FIG. 11 is a flowchart showing an outline of processing in the fifth embodiment.

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

1 能動型騒音制御装置(適応制御装置) 4 エンジン(騒音源) 5 クランク角センサ(騒音発生状態検出手
段) 6 車室(空間) 7a〜7d ラウドスピーカ(制御アクチュエータ、制
御音源) 8a〜8d マイクロフォン(残留騒音検出手段) 10 コントローラ 12 ディジタルフィルタ(適応処理手段、伝達
関数モデル) 13 適応ディジタルフィルタ(フィルタ係数可
変のディジタルフィルタ) 16 マイクロプロセッサ(適応処理手段) 20 伝達関数補正フィルタ
DESCRIPTION OF SYMBOLS 1 Active noise control device (adaptive control device) 4 Engine (noise source) 5 Crank angle sensor (noise generation state detection means) 6 Cabin (space) 7a-7d Loudspeaker (control actuator, control sound source) 8a-8d Microphone (Residual noise detection means) 10 Controller 12 Digital filter (adaptive processing means, transfer function model) 13 Adaptive digital filter (digital filter with variable filter coefficient) 16 Microprocessor (adaptive processing means) 20 Transfer function correction filter

Claims (11)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 所定の基準信号に応じて制御アクチュエ
ータを駆動する信号を生成するフィルタ係数可変のディ
ジタルフィルタと、所定の評価関数の値が目標値に一致
するように前記ディジタルフィルタのフィルタ係数を更
新する適応処理手段と、を備えた適応制御装置におい
て、 前記ディジタルフィルタのフィルタ係数が最適値に収束
したか否かを判定する最適値収束判定手段と、 この最適値収束判定手段が最適値に収束したと判定した
際の前記ディジタルフィルタのフィルタ係数から判る最
適値への最短の収束経路、及び最適値に収束していな
い前記ディジタルフィルタのフィルタ係数から判る実際
の収束経路に基づいて、前記実際の収束経路が前記最短
の収束経路に近づくように、前記適応処理手段の処理内
容を補正する補正手段と、を設けたことを特徴とする適
応制御装置。
1. A digital filter having a variable filter coefficient for generating a signal for driving a control actuator according to a predetermined reference signal, and a filter coefficient of the digital filter so that the value of a predetermined evaluation function matches a target value. In an adaptive control device including an adaptive processing unit for updating, an optimum value convergence determining unit that determines whether or not the filter coefficient of the digital filter has converged to an optimum value, and the optimum value convergence determining unit sets the optimum value to an optimum value. The maximum that can be known from the filter coefficient of the digital filter when it is determined that it has converged.
The shortest convergence path to an appropriate value, and the actual fact that can be understood from the filter coefficient of the digital filter that has not converged to the optimum value
Based on the convergence path, the actual convergence path the shortest
An adaptive control device for correcting the processing content of the adaptive processing means so as to approach the convergent path .
【請求項2】 補正手段は、前記実際の収束経路が前記
最短の収束経路に近づくように、ディジタルフィルタの
フィルタ係数の更新ゲインを変更するゲイン調整手段を
有する請求項1記載の適応制御装置。
2. The correction means is such that the actual convergence path is the
2. The adaptive control device according to claim 1, further comprising gain adjusting means for changing the update gain of the filter coefficient of the digital filter so as to approach the shortest convergence path .
【請求項3】 補正手段は、最適値収束判定手段が最適
値に収束したと判定した際のディジタルフィルタのフィ
ルタ係数が変動している場合には、適応処理手段の処理
内容を補正しない請求項1又は請求項2記載の適応制御
装置。
3. The correction means does not correct the processing content of the adaptive processing means when the filter coefficient of the digital filter changes when the optimum value convergence determination means determines that the optimum value has converged to the optimum value. 1 or the adaptive control device according to claim 2.
【請求項4】 補正手段は、最適値収束判定手段が最適
値に収束したと判定した際のディジタルフィルタのフィ
ルタ係数及び制御開始時の前記ディジタルフィルタのフ
ィルタ係数から決まる方向と、最適値に収束していない
前記ディジタルフィルタのフィルタ係数及び制御開始時
の前記ディジタルフィルタのフィルタ係数から決まる方
向との間のずれが無くなるように適応処理手段の処理内
容を補正する請求項1乃至請求項3のいずれかに記載の
適応制御装置。
4. The correction means converges to an optimum value and a direction determined by the filter coefficient of the digital filter when the optimum value convergence determination means determines that the optimum value has converged to the optimum value and the filter coefficient of the digital filter at the start of control. 4. The processing content of the adaptive processing means is corrected so that the deviation between the filter coefficient of the digital filter which is not performed and the direction determined by the filter coefficient of the digital filter at the start of control is eliminated. An adaptive control device according to claim 1.
【請求項5】 騒音源から騒音が伝達される空間に制御
音を発生可能な制御音源と、前記騒音源の騒音発生状態
を検出し基準信号として出力する騒音発生状態検出手段
と、前記基準信号に応じて前記制御音源を駆動する信号
を生成するフィルタ係数可変のディジタルフィルタと、
前記空間内の所定位置における残留騒音を検出する残留
騒音検出手段と、前記残留騒音に基づいた所定の評価関
数の値が低減するように前記ディジタルフィルタのフィ
ルタ係数を更新する適応処理手段と、を備えた能動型騒
音制御装置において、 前記ディジタルフィルタのフィルタ係数が最適値に収束
したか否かを判定する最適値収束判定手段と、 この最適値収束判定手段が最適値に収束したと判定した
際の前記ディジタルフィルタのフィルタ係数から判る最
適値への最短の収束経路、及び最適値に収束していな
い前記ディジタルフィルタのフィルタ係数から判る実際
の収束経路に基づいて、前記実際の収束経路が前記最短
の収束経路に近づくように、前記適応処理手段の処理内
容を補正する補正手段と、を設けたことを特徴とする能
動型騒音制御装置。
5. A control sound source capable of generating a control sound in a space in which noise is transmitted from a noise source, a noise generation state detecting means for detecting a noise generation state of the noise source and outputting it as a reference signal, and the reference signal. A digital filter having a variable filter coefficient for generating a signal for driving the control sound source according to
Residual noise detecting means for detecting residual noise at a predetermined position in the space, and adaptive processing means for updating the filter coefficient of the digital filter so that the value of a predetermined evaluation function based on the residual noise is reduced. In an active noise control device provided with, an optimum value convergence determination means for determining whether or not the filter coefficient of the digital filter has converged to an optimum value, and when it is determined that this optimum value convergence determination means has converged to an optimum value. most apparent from the filter coefficients of the digital filter
The shortest convergence path to an appropriate value, and the actual fact that can be understood from the filter coefficient of the digital filter that has not converged to the optimum value
Based on the convergence path, the actual convergence path the shortest
And a correction unit that corrects the processing content of the adaptive processing unit so as to approach the convergence path of the active noise control device.
【請求項6】 適応処理手段は、制御音源及び残留騒音
検出手段間の音響伝達特性に相当する伝達関数モデルを
有し、補正手段は、前記実際の収束経路が前記最短の収
束経路に近づくように、前記伝達関数モデルを補正する
請求項5記載の能動型騒音制御装置。
6. The adaptive processing means has a transfer function model corresponding to the acoustic transfer characteristic between the control sound source and the residual noise detection means, and the correction means has the actual convergence path with the shortest convergence.
The active noise control device according to claim 5 , wherein the transfer function model is corrected so as to approach the bundle path .
【請求項7】 補正手段は、適応処理手段に含まれる伝
達関数モデルと直列に配設された伝達特性可変の伝達関
数補正フィルタと、最適値収束判定手段が最適値に収束
したと判定した際の前記ディジタルフィルタのフィルタ
係数及び最適値に収束していない前記ディジタルフィル
タのフィルタ係数に基づいて、前記実際の収束経路が前
記最短の収束経路に近づくように、前記伝達関数補正フ
ィルタの伝達特性を変更する伝達特性制御手段と、を有
する請求項6記載の能動型騒音制御装置。
7. The correction means includes a transfer function correction filter arranged in series with a transfer function model included in the adaptive processing means and having a variable transfer characteristic, and an optimum value convergence judgment means when it judges that the optimum value converges to an optimum value. on the basis of the filter coefficient of the digital filter in the filter coefficient and the optimum value it does not converge to the digital filter, the actual convergence path before
7. The active noise control device according to claim 6 , further comprising a transfer characteristic control unit that changes a transfer characteristic of the transfer function correction filter so as to approach the shortest convergence path .
【請求項8】 補正手段は、ディジタルフィルタと直列
に配設された伝達特性可変の伝達関数補正フィルタと、
最適値収束判定手段が最適値に収束したと判定した際の
前記ディジタルフィルタのフィルタ係数及び最適値に収
束していない前記ディジタルフィルタのフィルタ係数に
基づいて、前記実際の収束経路が前記最短の収束経路に
近づくように、前記伝達関数補正フィルタの伝達特性を
変更する伝達特性制御手段と、を有する請求項6記載の
能動型騒音制御装置。
8. The correction means comprises a transfer function correction filter having a variable transfer characteristic arranged in series with a digital filter,
Based on the filter coefficient of the digital filter when the optimum value convergence determination means determines that it has converged to the optimum value and the filter coefficient of the digital filter that has not converged to the optimum value, the actual convergence path is the shortest convergence. On the route
The active noise control device according to claim 6 , further comprising a transfer characteristic control unit that changes a transfer characteristic of the transfer function correction filter so as to approach .
【請求項9】 補正手段は、前記実際の収束経路が前記
最短の収束経路に近 づくように、ディジタルフィルタの
フィルタ係数の更新ゲインを変更するゲイン調整手段を
有する請求項5記載の能動型騒音制御装置。
9. The correction means is such that the actual convergence path is the
As approached, the shortest convergence path, active noise control apparatus according to claim 5, further comprising a gain adjustment means for changing the update gain of the filter coefficient of the digital filter.
【請求項10】 補正手段は、最適値収束判定手段が最
適値に収束したと判定した際のディジタルフィルタのフ
ィルタ係数が変動している場合には、適応処理手段の処
理内容を補正しない請求項5乃至請求項9のいずれかに
記載の能動型騒音制御装置。
10. The correction means does not correct the processing content of the adaptive processing means when the filter coefficient of the digital filter changes when the optimum value convergence determination means determines that the optimum value has converged to the optimum value. The active noise control device according to any one of claims 5 to 9.
【請求項11】 補正手段は、最適値収束判定手段が最
適値に収束したと判定した際のディジタルフィルタのフ
ィルタ係数及び制御開始時の前記ディジタルフィルタの
フィルタ係数から決まる方向と、制御開始直後の前記デ
ィジタルフィルタのフィルタ係数及び制御開始時の前記
ディジタルフィルタのフィルタ係数から決まる方向との
間のずれが無くなるように適応処理手段の処理内容を補
正する請求項5乃至請求項10のいずれかに記載の能動
型騒音制御装置。
11. The correcting means determines the direction determined by the filter coefficient of the digital filter when the optimum value convergence determining means determines that the optimum value has converged to the optimum value and the filter coefficient of the digital filter at the start of control, and immediately after the start of control. 11. The processing content of the adaptive processing means is corrected so that there is no deviation between the filter coefficient of the digital filter and the direction determined by the filter coefficient of the digital filter at the start of control. Active noise control device.
JP06401292A 1992-03-19 1992-03-19 Adaptive control device and active noise control device Expired - Fee Related JP3517883B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP06401292A JP3517883B2 (en) 1992-03-19 1992-03-19 Adaptive control device and active noise control device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP06401292A JP3517883B2 (en) 1992-03-19 1992-03-19 Adaptive control device and active noise control device

Publications (2)

Publication Number Publication Date
JPH05265467A JPH05265467A (en) 1993-10-15
JP3517883B2 true JP3517883B2 (en) 2004-04-12

Family

ID=13245842

Family Applications (1)

Application Number Title Priority Date Filing Date
JP06401292A Expired - Fee Related JP3517883B2 (en) 1992-03-19 1992-03-19 Adaptive control device and active noise control device

Country Status (1)

Country Link
JP (1) JP3517883B2 (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3384493B2 (en) * 1992-04-03 2003-03-10 富士重工業株式会社 Interior noise reduction device
JPH0764567A (en) * 1993-08-23 1995-03-10 Mitsubishi Electric Corp Active noise controller
US5978473A (en) * 1995-12-27 1999-11-02 Ericsson Inc. Gauging convergence of adaptive filters
JPH1011075A (en) * 1996-06-19 1998-01-16 Toa Corp Active silencer
JPH1049175A (en) * 1996-07-29 1998-02-20 Fuji Heavy Ind Ltd Active noise control device
US9240819B1 (en) * 2014-10-02 2016-01-19 Bose Corporation Self-tuning transfer function for adaptive filtering
JP2022088859A (en) * 2020-12-03 2022-06-15 いすゞ自動車株式会社 Parameter adjustment device

Also Published As

Publication number Publication date
JPH05265467A (en) 1993-10-15

Similar Documents

Publication Publication Date Title
JP3094517B2 (en) Active noise control device
JP2882170B2 (en) Active noise control device
JP2748626B2 (en) Active noise control device
JPH0561483A (en) Active type noise controller
JP3517883B2 (en) Adaptive control device and active noise control device
JPH0573073A (en) Active type noise controller
JPH07248784A (en) Active noise controller
JP2529745B2 (en) Active noise control device
JP3328946B2 (en) Active uncomfortable wave control device
JPH07281676A (en) Active type vibration and noise controller
JP2751685B2 (en) Active noise control device
JP3355706B2 (en) Adaptive control device
JP3264951B2 (en) Active noise control device
JP3489137B2 (en) Active noise control device
JP2757514B2 (en) Active noise control device
JPH06250674A (en) Active noise controller
JP3500643B2 (en) Active noise control device
JP3674963B2 (en) Active noise control device and active vibration control device
JPH06130970A (en) Active noise controller
JP3517886B2 (en) Active noise control device
JPH0784585A (en) Active noise controller
JPH05289679A (en) Active type noise controller
JPH05232971A (en) Active type noise control device
JP3275449B2 (en) Active noise control device and active vibration control device
JP2572644Y2 (en) Vehicle interior noise reduction device

Legal Events

Date Code Title Description
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20040106

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20040119

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090206

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100206

Year of fee payment: 6

LAPS Cancellation because of no payment of annual fees