JPH06337685A

JPH06337685A - 車両の振動低減装置

Info

Publication number: JPH06337685A
Application number: JP5147033A
Authority: JP
Inventors: Hiroshi Uchida; 博志内田; Tetsurou Butsuen; 哲朗仏圓; Norihiko Nakao; 憲彦中尾
Original assignee: Mazda Motor Corp
Current assignee: Mazda Motor Corp
Priority date: 1993-05-26
Filing date: 1993-05-26
Publication date: 1994-12-06
Anticipated expiration: 2017-05-13
Also published as: JP3281121B2

Abstract

(57)【要約】【目的】低減用振動生成に用いられるマイク／スピ−カ
間のインパルス応答デ−タ長を短くする。【構成】車室２に、車室内振動を低減するための低減用
振動を出力するスピ−カ１１と、車室内振動を検出する
マイク１２とが設けられる。マイク１２で検出される振
動が低減されるように、スピ−カ１１とマイク１２との
間のインパルス応答を考慮しつつ、低減用振動が最適化
される。インパルス応答デ−タから、スピ−カ１１の機
械的共振周波数以下の低周波数成分が、ハイパスフィル
タ等を利用してカットされる。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、車両の振動つまり騒音
を、低減用振動を利用した干渉作用によって低減するよ
うにした車両の振動低減装置に関するものである。

【０００２】

【従来の技術】車両、特にエンジンによる騒音振動つま
り第１振動が問題になる自動車等においては、スピ−カ
等から低減用振動つまり第２振動を発生させて、この第
１振動と第２振動との干渉により第１振動を低減するこ
とが提案されている。

【０００３】この種の振動低減装置にあっては、特表平
１−５０１３４４号公報に示すように、振動源からの振
動つまり第１振動に相当する信号をリファレンス信号と
して取り出すリファレンス信号発生器と、第１振動によ
る騒音が問題となる所定空間での振動をピックアップす
るマイクと、所定空間に向けて第２振動を発生させるス
ピ−カと、スピ−カから出力させる第２振動を生成する
ための適応型デジタルフィルタと、上記フィルタのフィ
ルタ係数を逐次的に最適化するためのアルゴリズム演算
装置と、を有する。すなわち、リファレンス信号に応じ
て適応型デジタルフィルタがリファレンス信号のゲイン
や位相等を調整して第２振動を生成する一方、マイクで
検出される振動が小さくなるように、適応型デジタルフ
ィルタのフィルタ係数がアルゴリズム演算装置によって
逐次的に最適化される。そして、低減用振動の最適化の
ために、マイクらとスピ−カとの間のインパルス応答が
加味されたものとされる。

【０００４】

【発明が解決しようとする課題】ところで、低減用振動
の最適化のために用いられるインパルス応答のデ−タ
は、そのデ−タ長が短いほど、振動の収束性向上、低減
用振動生成のための計算時間短縮化、メモリの負担軽減
となって好ましいものである。この一方、スピ−カは、
所定の振動に対応した信号を入力させた際、この信号の
入力が停止された後もしばらくの間は、スピ−カから低
周波数成分つまり残留低周波数成分が出力され続けるこ
とになる。この残留低周波数成分は、スピ−カの機械的
共振周に強く依存するものであり、これを無くすこと自
体は事実上不可能である。そして、従来は、この残留低
周波数成分を加味してインパルス応答デ−タを生成して
いたので、そのデ−タ長を短くするには限度があった。

【０００５】本発明は以上のような事情を勘案してなさ
れたもので、低減用振動を最適化するためのインパルス
応答のデ−タ長を十分短くできるようにした車両の振動
低減装置を提供することにある。

【０００６】

【課題を解決するための手段】上記目的を達成するた
め、本発明にあっては次のような構成としてある。すな
わち、車室内に低減用振動を出力させるためのスピ−カ
と、車室内の振動を検出するマイクと、該マイクで検出
される振動が低減されるように前記スピ−カとマイクと
の間のインパルス応答を考慮しつつ前記低減用振動を最
適化する最適化手段と、を備えた車両用振動制御装置に
おいて、前記インパルス応答のデ−タから、前記スピ−
カの機械的共振周波数以下の低周波数成分をカットする
修正手段を備えている、ような構成としてある。

【０００７】上記構成を前提とした本発明の好ましい態
様は、特許請求の範囲における請求項２以下に記載した
通りである。

【０００８】

【発明の効果】請求項１に記載された本発明によれば、
振動低減のためには事実上問題とならないスピ−カの残
留低周波数成分について、インパルス応答デ−タからカ
ットするようにしてあるので、このインパルス応答デ−
タ長を十分短くして、振動の収束性向上、低減用振動最
適化のための計算時間短縮、メモリの負担軽減等を得る
ことができる。

【０００９】請求項２に記載したような構成とすること
により、振動低減の効果をほとんど損なうことなく、イ
ンパルス応答デ−タ長をより短くして、請求項１に記載
された本発明の効果をより一層向上させることができ
る。

【００１０】請求項３に記載したような構成とすること
により、ハイパスフィルタを利用して、容易に修正手段
を構成することができる。

【００１１】請求項４に記載したような構成とすること
により、低減用振動生成のための計算量の低減ひいては
制御系の負担を軽減する上で好ましいものとなる。

【００１２】

【実施例】以下本発明の実施例を添付した図面に基づい
て説明する。全体の概要図１において、自動車１は、車室２内に運転席３と助手
席４と左右の後席５、６とを有する４人乗りの乗用車と
されている。車体前部に構成されたエンジンル−ム７に
は、直列４気筒のガソリンエンジン８が塔載され、その
イグニッションコイルが符号９で示される。

【００１３】エンジン８が、エンジン回転数に応じた周
期的な振動を発生する騒音発生源つまり第１振動源とさ
れている。そして、車室２が、エンジン８の振動を低減
すべき所定空間とされている。このため、所定空間とし
ての車室２には、５個のスピ−カ１１と、８個のマイク
１２とが設置されている。スピ−カ１１が、車室へエン
ジン騒音を低減するための第２振動を発生する第２振動
源とされ、マイク１２が、車室の実際の振動を検出する
振動検出手段とされる。そして、インストルメントパネ
ルには、既知のように、例えばチュ−ナとカセットデッ
キとＣＤデッキとの組み合わせからなるオ−ディオソ−
ス３１が配設されている。なお、スピ−カ１１は、オ−
ディオ用と兼用することなく、低減用振動発生用として
専用に設けてもよい。

【００１４】自動車１には、マイクロコンピュ−タを利
用して構成された制御ユニットＵが塔載されている。制
御ユニットＵに対する入出力関係を図２に示してあり、
制御ユニットＵは、ＣＰＵからなる制御部２０を有す
る。制御部２０には、イグニッションコイル９の一次コ
イルからの信号つまりエンジン回転数に応じた点火パル
ス信号が、波形整形回路２１、周期計算回路２２を経て
入力されると共に、各マイク１２からの信号が、アンプ
２３、ロ−パスフィルタ２４、Ａ／Ｄ変換器２５を介し
て入力される。また、制御部２０からの出力信号は、Ｄ
／Ａ変換器２６、ロ−パスフィルタ２７、アンプ２８を
介してスピ−カ１１へ出力される。

【００１５】制御部２０は、マイク１２で検出される振
動が低減されるように、スピ−カ１１から出力すべき低
減用振動としての第２振動を最適化する。以下、制御部
２０による第２振動の生成について説明するが、先ず、
第２振動の生成の基本部分ついて説明し、その後インパ
ルス応答のデ−タ長修正の点について説明する。

【００１６】なお、実施例では、第２振動は、エンジン
８の周期的な回転振動例えば回転２次成分を低減するも
ので、エンジン８の周期的な振動の１周期分まとめて生
成するようにして、第２振動最適化のための計算量の低
減ひいては制御系の負担が極力小さくなるようにしてあ
る。勿論、第２振動の生成は、既知の適宜の最適化手法
によりなし得るものであり、本発明は特定の最適化手法
に限定されないものではない。

【００１７】第２振動の生成（基本）図３は、制御部２０をブロック図的に示すものであり、
説明の簡単化のためにスピ−カ１１およびマイク１２を
それぞれ１個とした場合を示している。

【００１８】制御部２０は、周期計測回路２２から入力
された結果によってスピ−カ１１に出力するスピ−カ入
力信号ｙのベクトルｙの周期を調整する（ステップ１、
以下ステップをＳと略す）と共に、内蔵しているプロセ
ッサで、マイク１２・スピ−カ２間の伝達特性であるイ
ンパルス応答ｈの行列ｈを、時系列ｈ変換する（Ｓ
２）。

【００１９】次に、制御部２０はプロセッサで、インパ
ルス応答ｈの時系列ｈとマイク１２から入力されるマイ
ク出力信号ｅとでベクトルｙを逐次的に最適化し（Ｓ
３）、その後、このベクトルｙを時系列ｙに変換してス
ピ−カ入力信号ｙとし（Ｓ４）、スピ−カ１１に出力す
る。

【００２０】スピ−カ１１は、このスピ−カ入力信号ｙ
をアンチ騒音Ｚとして再生する。一方、マイク１２は、
騒音ｄとアンチ騒音Ｚが打ち消し合って振動エネルギが
低減した騒音を検出して、この結果をディジタルのマイ
ク出力信号ｅとして制御部２０に内蔵されたプロセッサ
に出力する。以下、再びプロセッサは、上記ステップ３
およびステップ４を繰り返し行い、スピ−カ入力信号ｙ
のベクトルｙを逐次的に最適化して、最終的にマイク出
力信号ｅの値が０となるようにスピ−カ入力信号ｙのベ
クトルｙを設定する。

【００２１】次に、制御部２０で行われる上記ステップ
のアルゴリズムの演算について、以下に説明する。

【００２２】先ず、制御部２０によるマイク１２のマイ
ク出力信号ｅのサンプリング周期を△ｔとする。マイク
１２・スピ−カ１１間の伝達特性であるインパルス応答
ｈが有限時間Ｊ△ｔ以内で０に収束すると仮定し、イン
パルス入力が与えられてからｊ△ｔ時間経過後のインパ
ルス応答ｈの値をｈ_j とすると、エンジン８から発生し
た第１振動である騒音ｄ、スピ−カ入力信号ｙが与えら
れたときのスピ−カ１１から発生する第２振動であるア
ンチ騒音Ｚおよびそのときの時刻ｋにおけるマイク出力
信号ｅの第ｋサンプル値ｅ（k)の関係は、次式（１）で
表わすことができる。

【００２３】ｅ(k) ＝ｄ（k)＋Ｚ(k) ＝ｄ（k)＋行列ｈ^T ・行列ｙ（k) ・・・・（１）但し、行列ｈ＝［ｈ₀ ｈ₁ ｈ₂ ・・・・・ｈ_J-1 ］^T 行列ｙ（k)＝［ｙ（k) ｙ（k-1) ｙ（k-2)・・・・ｙ
（k-J+1)］^T ｄ(k):ｅ（k)に含まれている騒音ｄの成分Ｚ(k):ｅ（k)に含まれているアンチ騒音Ｚの成分ｙ(k):スピ−カ入力信号ｙの第ｋサンプル値従って、式（１）中のＺ（k)は、次の式（２）で示され
る。

【００２４】

【数１】

【００２５】ところで、騒音ｄは、ある周期Ｎ△ｔを持
っている周期性騒音であるので、この騒音ｄの振動エネ
ルギを低減させるアンチ騒音Ｚおよびスピ−カ入力信号
ｙ、騒音ｄと同じ周期Ｎ△ｔを持っている周期性振動お
よび周期性信号でなければならない。

【００２６】従って、スピ−カ入力信号ｙに関して次式
（３）が成立する。ｙ(k) ＝ｙ(K-qN)＝ｙ(k) ｙ (k-1)＝ｙ(k-qN-1)＝ｙ(k+N-1) ｙ(k-2) ＝ｙ(k-qN-2)＝ｙ(k+N-2) ・・・・（３）・・・・・・・・・ｙ(k-N+1) ＝ｙ(k-(q+1)N+1)＝ｙ(k+1) 但し、ｑ＝０，１，２，・・・・ゆえに、式（１）は、ｅ(k) ＝ｄ(k)+ベクトルｈ^T ・時系列ｙ(k) ・・・・（４）但し、時系列ｙ(k) ＝［ｙ (K) ｙ(K+N-1) ｙ(K+N-2) ・・
・・ｙ(K+1)]^T

【００２７】

【数２】

【００２８】尚、Ｑは、Ｊ≦(q+1)Nを満たす整数ｑの最
小値である。

【００２９】次に、時刻ｋからさらにｉだけ時間が経過
した時刻k+i のマイク出力信号ｅの第K+i サンプル値ｅ
（K+i)（但し、ｉ＝１，２，・・・・）は、次式（５）
で表わすことができる。

【００３０】ｅ（k+i)＝ｄ(k+i) ＋ベクトルｈ^T ・時系列ｙ(k+i) ＝ｄ(k+i) ＋時系列ｈ(i)^T・時系列ｙ(k) ・・・・・（５）但し、時系列ｙ(k+i) ＝［ｙ(k+i)^'ｙ(k+i^'-1 ）ｙ(k+N-1) ｙ(k+N-2) ・・・・・ｙ(k+i^'+1)]^T 時系列ｈ(i) ＝［バ−ｈ_i ^'バ−ｈ_i+1 ^'・・・・・バ−ｈ
_N+1 バ−ｈ₀バ−ｈ₁ ・・・・バ−ｈ_i ^' _-1］^T 尚、ｉ’は、ｉをＮで割ったときの整数剰余である。

【００３１】ところで、式（５）において、ｋはマイク
入力信号ｅの任意の初期時点を表わしているに過ぎな
い。よって、ｋ＝０と置き、ｉを改めてｋに置き直す
と、次式（６）が得られる。

【００３２】ｅ(k) ＝ｄ(k) ＋時系列ｈ(k)^T・時系列ｙ(0) ＝ｄ(k) ＋時系列ｈ(k)^T・ベクトルｙ但し、ベクトルｙ＝［ｙ(0) ｙ(N-1) ｙ(N-2) ・・・ｙ(1) ］^T ＝［ｙ₀ ｙ_N-1 ｙ_N-2 ・・・・ｙ₁ ］^T ここで、次の評価関数を導入する。Ｆ＝Ｅ［ｅ(k)²］＝Ｅ［ｄ(k) ＋時系列ｈ(k)^T・ベクトルｙ］＝Ｅ［ｄ(k)²］＋２ベクトルｙ^T ・Ｅ［ｄ(k) ・時系列ｈ(k) ］＋ベクトルｙ^T ・Ｅ［時系列ｈ(k) ・時系列ｈ(k)^T］ベクトルｙ・・・・・・（７）但し、Ｅ［］は、期待値を表わすものとする（Ｅは期
待演算子）。式（７）より、この評価関数のベクトルｙ
に関する勾配は、次式（８）で与えられる。 ∂Ｆ／∂ベクトルｙ＝２Ｅ［ｄ(k) ・時系列ｈ(k)] ＋２Ｅ［時系列ｈ(k) ・時系列ｈ(k)^T］ベクトルｙ＝２Ｅ［時系列ｈ(k){d(k)＋時系列ｈ(k)^Tベクトルｙ｝］＝２Ｅ［時系列ｈ(k) ・ｅ(k) ］・・・・・（８）ここで、Ｅ［時系列ｈ(k) ・ｅ(K)]の瞬時推定値とし
て、時系列ｈ(k)・e (K)を用いることにすれば、Ｆの最
小値を与える周期Ｎ△ｔ（すなわち要素数Ｎ）を持つス
ピ−カ出力信号ベクトルであるベクトルｙの値は、最急
降下法に基づく次の漸化式（９）を反復計算することに
にょり最適化することができる。

【００３３】ベクトルｙ(K+1) ＝ベクトルｙ(k) −μ・ｅ(k) ・時系列ｈ(k) ・・・（９）但し、μ／２は収束係数である。

【００３４】このようにして求めた漸化式（９）は、制
御部２０に内蔵されたデ−タ処理装置であるプロセッサ
が騒音の振動エネルギを低減させるアンチ騒音の振動エ
ネルギの設定を補正する際には、以下に示すような、よ
り簡単なアルゴリズムに置き換えられる。

【００３５】先ず、一対のスピ−カ１１およびマイク１
２を用いる場合には、漸化式（９）は次式（１０）に置
き換えられる。ｙ_(k-j+QN) ^'(k+1)＝ｙ_(k-j+QN) ^' ・(k) −μ・ｅ(k) ・ｈ_j ・・・（１０）このときプロセッサは、時刻ｋにおいては、例えば以下
に示す４つの動作手順を行っている。

【００３６】動作１：スピ−カ入力信号ｙ_k ^' (k)をスピ
−カ１１に対して出力する。動作２：マイク出力信号ｅ(K) をマイク１２から入力す
る。動作３：周期計測回路２２から入力されたエンジン２２
の回転周期にＯrd／△ｔまたは１／（Ｏrd・△ｔ）を乗
じた値に最も近い整数値をＮとする。動作４：ｊ＝０，１，２，・・・・，Ｊ−１について漸
化式（１０）の計算を行う。但し、ｋ’，（ｋ−ｊ＋ＱＮ）’は、それぞれｋ（ｋ−
ｊ＋ＱＮ）をＮで，割ったときの整数剰余であり、ま
た、Ｏrdは、低減させようとしている騒音のエンジン回
転数に対する最低次数を設定するための任意の一定の整
数である。

【００３７】次に、複数のスピ−カ１１・・・とマイク
１２・・・とを用いる場合には、例えば、最急降下法に
基づき、

【００３８】

【数３】

【００３９】の瞬時推定値として、

【００４０】

【数４】

【００４１】を用いると、評価関数

【００４２】

【数５】

【００４３】を最小化する第１スピ−カ出力信号ベクト
ルであるベクトルｙ₁ の最適値は、次の漸化式（１１）
を反復計算することにより求められる。

【００４４】

【数６】

【００４５】但し、ｙ_lk ^' ：時刻ｋにおける第１スピ−カ入力信号ｅ _m ：第ｍマイク出力信号ｈ_lmj ：第１スピ−カ・第ｍマイク間のインパルス応答
のｊ△ｔ時間後の値Ｌ：スピ−カの個数Ｍ：マイクの個数Ｊ：全てのスピ−カ・マイク間のインパルス応答が有限
時間△ｔ以内で０に収束することを示す整数値また、ベクトルｙ_l ＝［ｙ_l ₀ ｙ_l _N-1 ｙ_l _N-2 ・・・ｙ
_l ₁ ］^T 時系列ｈ_lm(k) ＝［バ−ｈ_lm _k' バ−ｈ_lm _k'+1 ・・・
バ−ｈ_lm _N+1バ−ｈ_lm ₀ バ−ｈ_lm ₁・・・バ−ｈ_lm
_k'-1 ］^T さらに、バ−ｈ_lm ₀＝ｈ_lm ₀ ＋ｈ_lm _N ＋・・・・ｈ_lm _QN バ−ｈ_lm ₁＝ｈ_lm ₁ ＋ｈ_lm _N+1＋・・・ｈ_lm _QN+1 ・・・・・・・・・・・・・・・バ−ｈ_lm _j-QN-1 ＝ｈ_lm _j-QN-1 ＋ｈ_lm _j-(Q-1)N-1 ＋・・・＋ｈ_lm _j-1 バ−ｈ_lm _j-QN ＝ｈ_lm _j-QN ＋ｈ_lm _j-(Q-1)N ＋・・・＋０・・・・・・・・・・・・・・・・・・バ−ｈ_lm _N-1 ＝ｈ_lm _N-1 ＋ｈ_lm _2N-1 ＋・・・＋０ｌ＝１，２，・・・・，Ｌｍ＝１，２，・・・・，Ｍ従って、漸化式（９）は次式（１２）に置き換えられ
る。

【００４６】

【数７】

【００４７】このときプロセッサは、時刻ｋにおいて
は、例えば以下に示す４つの動作手順を行っている。

【００４８】動作１１：スピ−カ入力信号ｙ_1k ^' （k)，
ｙ_2k ^' （k)，・・・・，ｙ_lk ^'(k ）をそれぞれ第１スピ
−カ、第２スピ−カ、・・・、第Ｌスピ−カに対して出
力する。動作１２：マイク出力信号ｅ₁(k), ｅ₂(k),・・・, ｅ_M(k)
をそれぞれ第１マイク、第２マイク、・・・・、第Ｍマ
イクから入力する。動作１３：周期計測回路２２から入力されたエンジン２
２の回転周期にＯrd／△ｔまたは１／（Ｏrd・ △ｔ）を
乗じた値に最も近い整数値をＮとする。動作１４：１＝１、２、・・・・・Ｌおよびｊ＝０，
１，２，・・・・Ｊ−１について漸化式（１２）の計算
を行う。また、上記の複数のスピ−カ１１・・・とマイク１２・
・・とを用いる場合について、

【００４９】

【数８】

【００５０】の瞬時推定値として、α・時系列ｈ_1k ^'(k)
・ｅ_k ^' (k)を用いると、最急降下法に基づいて評価関数

【００５１】

【数９】

【００５２】を最小化する第１スピ−カ出力信号ベクト
ルであるベクトルｙ₁ の最適値は、次の漸化式（１３）
を反復計算することにより求められる。ベクトルｙ₁ (k+ 1)＝ベクトルｙ₁ (k) −μ・α・時系列ｈ_1k ^"(k)・ｅ_k ^"(k) ・・・・（１３）但し、ｋ”は、ｋをＭで割ったときの整数剰余に１を加
えた値であり、また、αは任意の定数である。この漸化
式（１３）は、漸化式（１１）よりも短時間で演算でき
る。

【００５３】従って、漸化式（９）は次式（１４）に置
き換えられる。ｙ_1(k-J+QN) ^'(k+1) ＝ｙ_1(K-j+QN) ^'（k)−μ・α・ｅ_k(k)・ｈ_1k ^" _j ・・・・・（１４）このときプロセッサは、時刻においては、例えば以下に
示す４つの動作手順を行っている。

【００５４】動作２１：スピ−カ入力信号ｙ_1k ^'(k), ｙ
_2k ^'(k), ・・・・、ｙ_Lk ^'(k ）をそれぞれ第１スピ−
カ、第２スピ−カ、・・・・・、第Ｌスピ−カに対して
出力する。動作２２：マイク出力信号ｅ_k ^"(k) を第ｋ”マイクから
入力する。動作２３：周期計測回路２２から入力されたエンジン２
２の回転周期にＯrd／△ｔまたは１／（Ｏrd・△ｔ）を
乗じた値に最も近い整数値をＮとする。動作２４：１＝１、２、・・・・、Ｌおよびｊ＝０、
１、２・・・・、Ｊ−１について漸化式（１４）の計算
を行う。従って、上記アルゴリズムの演算は、漸化式（９）、
（１１）および（１３）、あるいはこれら漸化式を単純
化した漸化式（１０）、（１２）および（１４）を反復
計算するだけで良いので、スピ−カ入力制御の計算時間
を短縮することが可能となる。

【００５５】インパルス応答デ−タ長の修正次に、図４以下を参照しつつ、インパルス応答デ−タ長
を修正する点について説明する。先ず、図４は、修正前
のインパルス応答の例を示してあり、残留低周波数成分
がかなり長い時間に渡って生じることが理解される。図
５は、周波数とスピ−カゲインとの関係を示してあり、
図中ｆ１で示す周波数が、スピ−カ１１の機械的共振周
波数よりも若干高い周波数を示しており、この周波数ｆ
１以下の低周波数成分が、インパルス応答デ−タからカ
ットされる対象となる。さらに、図６は、エンジン音ト
ラッキングデ−タで、この図６における周波数ｆ２以下
の低周波数成分は、エンジン騒音として車室内で問題と
ならないものとなる。そして、この周波数ｆ２は、上記
周波数ｆ１よりも大きい点を考慮して、インパルス応答
デ−タの修正に際しては、この周波数ｆ２以下の低周波
数成分をカットするようにしてある。

【００５６】図７は、前記周波数ｆ２以下の低周波数成
分を、ハイパスフィルタ４１、４２を利用してカットす
る場合の回路例を示してある。この図７において、Ｓ
２、Ｓ３は図２に示すものと対応しており、時系列に変
換された測定直後のインパルス応答デ−タｈ₀ が、先ず
ハイパスフィルタ４１を通過されて時系列ｈ^' とされ、
このｈ^' が逆の方向からハイパスフィルタ４２を通過さ
れて、Ｓ３で用いる最終的な（修正後の）時系列ｈとさ
れる。このように、ハイパスフィルタ４２で逆方向から
通過させるのは、ハイパスフィルタ４１を通過したとき
に生じる位相ずれを、後に通過されるハイパスフィルタ
４２で元の位相に戻すためである。

【００５７】上記時系列ｈ^' は次の（15) 式で示され
る。

【００５８】

【数１０】

【００５９】ただし、測定されたインパルス応答デ−タ
ｈ₀ は（16) 式に、ハイパスフィルタ４１、４２の特性
ｇは（17) 式に、修正されたインパルス応答デ−タｈ^'
は（18) 式に示す通りである。ｈ₀ ＝［ｈ₀₍₀₎，ｈ₀₍₁₎，ｈ₀₍₂₎・・・ｈ_0(2J-1) ］・・・（１６）

【００６０】ｇ＝［ｇ₍₀₎ ，ｇ₍₁₎ ，ｇ₍₂₎ ・・・・ｇ_(J-1) ］・・・（１７）

【００６１】ｈ＝［ｈ₍₀₎ ，ｈ₍₁₎ ，ｈ₍₂₎ ・・・・ｈ_(2J-1)］・・・（１８）

【００６２】前記ハイパスフィルタ４１、４２の処理内
容を、図８のフロ−チャ−トに示してある。この図８に
おいて、Ｑはステップを示し、Ｑ１０までの処理がハイ
パスフィルタ４１による処理内容に相当し、Ｑ１１〜Ｑ
１８の処理がハイパスフィルタ４２の処理内容に相当
し、Ｑ１９以下が通常のデ−タに戻すための処理とな
る。

【００６３】以上実施例について説明したが、低減すべ
き振動としては、エンジン振動に限らず、車室内で問題
となる適宜の信号とすることができる。

【図面の簡単な説明】

【図１】本発明が適用された車両を上方から見た図。

【図２】低減用振動生成のための全体制御系統図。

【図３】図２のうち低減用振動の最適化部分の構成をブ
ロック図的に示す図。

【図４】スピ−カの残留低周波数成分を示す図。

【図５】周波数とスピ−カゲインと機械的共振周波数と
の関係を示す図。

【図６】エンジン音トラッキングデ−タとスピ−カ出力
との関係を示す図。

【図７】インパルス応答デ−タ長を修正するための処理
をブロック図的に示す図。

【図８】インパルス応答デ−タ長を修正するためのフロ
−チャ−ト。

【符号の説明】

１：自動車２：車室８：エンジン（振動源）１１：スピ−カ１２：マイク２０：低減用振動生成回路４１：ハイパスフィルタ４２：ハイパスフィルタ

Claims

【特許請求の範囲】

【請求項１】車室内に低減用振動を出力させるためのス
ピ−カと、車室内の振動を検出するマイクと、該マイク
で検出される振動が低減されるように前記スピ−カとマ
イクとの間のインパルス応答を考慮しつつ前記低減用振
動を最適化する最適化手段と、を備えた車両の振動低減
装置において、前記インパルス応答のデ−タから、前記スピ−カの機械
的共振周波数以下の低周波数成分をカットする修正手段
を備えている、ことを特徴とする車両の振動低減装置
【請求項２】請求項１において、前記修正手段が、前記機械的共振周波数よりも高い周波
数域でかつエンジン騒音が問題とならない低周波数成分
をもカットするもの。
【請求項３】請求項１において、前記修正手段による低周波数成分のカットが、インパル
ス応答デ−タを、１回ハイパスフィルタを通過させた後
さらに逆方向からハイパスフィルタを通過させることに
より行なわれるもの。
【請求項４】請求項１ないし請求項３のいずれか１項に
おいて、車室内で低減すべき振動が周期的な振動とされ、前記最適化手段が、前記周期的な振動の１周期毎に前記
低減用振動を最適化するもの。