JP3536313B2 - 車両の振動低減装置 - Google Patents

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 0} ｙ_{l N-1} ｙ_{l N-2} ・・・ｙ_{l 1} ］^T 時系列ｈ_lm(k) ＝［バ−ｈ_{lm k'} バ−ｈ_{lm k'+1} ・・・バ−ｈ_{lm N+1} バ−ｈ_{lm 0} バ−ｈ_{lm 1}・・・バ−ｈ_{lm k'-1} ］^T さらに、 バ−ｈ_{lm 0}＝ｈ_{lm 0} ＋ｈ_{lm N} ＋・・・・ｈ_{lm QN} バ−ｈ_{lm 1}＝ｈ_{lm 1} ＋ｈ_{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・△ｔ）を
乗じた値に最も近い整数値をＮとする。 動作２４：１＝１、２、・・・・、Ｌおよびｊ＝０、
１、２・・・・、Ｊ−１について漸化式（１４）の計算
を行う。 従って、上記アルゴリズムの演算は、漸化式（９）、
（１１）および（１３）、あるいはこれら漸化式を単純
化した漸化式（１０）、（１２）および（１４）を反復
計算するだけで良いので、スピ−カ入力制御の計算時間
を短縮することが可能となる。

【００６１】ニュ−トン法について これまでと同様に、Ｊ＝Ｅ［ｅ（ｋ）² ］を最小化する
ことを目的としつつ、ニュ−トン法によりベクトルｙを
最適化するための漸化式は次式（１５）となる。

【００６２】 ｙ（ｋ＋１）＝ｙ(k) −μ・Ｅ［ｈ(k）ｈ(k)^T］^-1・Ｅ［ｈ(k)・ｅ(k) ］ ・・・（１５） ただし、ｅ(k) ：マイク信号。 ｙ(k) ：出力波形データベクトル。 ｈ(k) ：スピ−カノマイク間インパルス応答系列ｈ
（ｋ）を出力波形データベクトルと同じデータ長の領域
内で順次上方向にロ−テ−トしながら得られるベクト
ル。 μ ：収束係数。

【００６３】上記Ｅ［ｈ(k)・ｅ(k) ］の瞬時推定値とし
てｈ(k)・ｅ(k) を用いるならば、（１５）式は次式（１
６）で示される。 ｙ（ｋ＋１）＝ｙ(k) −μ・Ｅ［ｈ(k）ｈ(k)^T］^-1・ｈ(k)・ｅ(k) ＝ｙ(k) −μ・ｇ(k)・ｅ(k) ・・・（１６） ただし、ｇ(k) ：２次勾配ベクトルの瞬時推定値。 ｇ(k) ＝Ｅ［ｈ(k）ｈ(k)^T］・ｈ(k) 以下これを単に２次勾配ベクトルと呼ぶ。

【００６４】ところで、ｈ(k) が一般的な信号系列であ
るならば、Ｅ［ｈ(k）ｈ(k)^T］は対称Toepritzと呼ばれ
る構造になるが、本方式の場合には、ｈ(k）は周期性を
有するために、Ｅ［ｈ(k）ｈ(k)^T］はさらに特殊な構造
を有する。すなわち、出力波形データ長ｎが偶数の場合
には、次式（１７）で示される。

【００６５】

【数１０】

【００６６】ただし、δ_T ＝Ｅ［ｈ(k）ｈ(k+ τ)^T］ ｍ＝ｎ／２ また、ｎが奇数の場合には、次式（１８）で示される。

【００６７】

【数１１】

【００６８】ただし、ｍ＝（ｎ−１）／２ つまりこれらの行列は、対称Toepritzであるという特徴
に加えて、第ｉ行は第ｉ−１行を一要素分だけ右ロ−テ
−トしたものになるという特徴を持っている。以下この
ようなデータ構造を『周期対称Toepritz』と呼ぶことに
する。なお、周期対称Toepritzの逆行列はやはり周期対
称Toepritzになるという重要な性質がある。（一般的な
対称Toepritzの逆行列は対称Toepritzとはならない）こ
のことから、Ｅ［ｈ(k）ｈ(k)^T］^-1を、次のような逆共
分散関数で次式（１９）として表すことにする。

【００６９】

【数１２】

【００７０】Ｅ［ｈ(k）ｈ(k)^T］においてδ₀,δ₁,δ₂,
…δ_m はｈ(k）の自己共分散関数（の一部）を表してい
るが、Ｅ［ｈ(k）ｈ(k)^T］^-1におけるγ₀,γ₁,γ₂,…γ
_m は、自己共分散関数との相互共分散関数が単位インパ
ルス系列を与えるものである。以下この系列（γ₀,γ₁,
γ₂,…γ_m ）をｈ(k）の逆共分散関数と呼ぶことにす
る。

【００７１】さて、（１６）式を実際の制御に用いるた
めには、以下の２点について検討しなければならない。 勾配ベクトルの各要素を順次ロ−テ−トさせつつ計算
するという、従来の計算アルゴリズムがそのまま適用可
能であるか否か（つまり、出力波の位相を変えてもｇ
(k) に含まれる各要素の値は変化せず、単にロ−テ−ト
するだけで済むかどうか）。 エンジン回転数が変化したとき、若干の修正を加える
だけで（例えば単に０をいくつか挿入するだけで）新し
いｇ(k) を得ることができるかどうか。以下、この２点
についての考察結果を述べる。

【００７２】の問題について 詳細については省略するが、ｇ(k) ＝Ｅ［ｈ(k）ｈ
(k)^T］^-1・ｈ(k) における各ベクトルと行列の中に含ま
れる要素を詳しく検討してみると、ｈ(k) がロ−テ−ト
したとき、ｇ(k) もやはり同じ方向にロ−テ−トするだ
けであることがわかる。これは、Ｅ［ｈ(k）ｈ(k)^T］^-1
が、周期対称Toepritzという特殊な構造を有しているた
めである。このことから、（１６）式のニュ−トン法の
採用にあたっては、従来の計算アルゴリズムにおいて、
単にｈ(k）をｇ(k) に置き換えるだけで、全く同じアル
ゴリズムを用いることができる。したがって、計算量を
全く増やすことなく、収束速度を大幅に改善することが
できる。

【００７３】の問題について 理論的にいえば、エンジン回転周期が変動することによ
って、逆共分散関数、ひいては逆共分散行列（Ｅ［ｈ
(k）ｈ(k)^T］^-1）、さらには勾配ベクトル（ｇ(k) ）の
いずれの要素も悉く変化する。変化前の値に若干の修正
を加えるだけで変化後の値を得るということはできな
い。従って厳密にいえば各エンジン回転数毎の勾配ベク
トルの値をあらかじめ全て求めておくか、そうでなけれ
ばその時々のエンジン回転数に応じてその都度新たに勾
配ベクトルを求めなければならない。

【００７４】しかしながら、エンジン回転数の変化に伴
って逆共分散関数や勾配ベクトルの各要素がどのように
変化するかを、インパルス応答系列の有効データ長（０
でない部分の長さ）が、出力波形データ長の１／２以下
の場合に限って詳しく見てみれば、ある種の傾向がある
ことがわかる。すなわち、逆共分散関数の傾向について
述べれば、 ｉ）エンジン回転数が低くなるにつれて、γ₀,γ₁,γ₂,
…は、この順序で（左端から順に）一定値に収束してゆ
く（なお、収束の模様は単調ではなく振動的である）。 ii）γ₀,γ₁,γ₂,…の収束値は、この順序で（右端に近
くなるほど）０に近付いてゆく。 iii) i) の収束は、インパルス応答の有効データ長と出
力波形データ長との比が小さいほど早く進む。

【００７５】つまり、インパルス応答持続時間がエンジ
ン回転周期に比べて十分短い場合には、逆共分散関数
は、γ₀,γ₁,γ₂,…γ_?,０，０，…０という形になり、
エンジン回転数がさらに低下したとき、後に０をいくつ
か追加するだけで済むようになる（理論的に言えばそう
ではないのだが、微小な変化分を無視するならば、こう
いう取扱いが可能である。コンピュ−タ上では離散値し
か扱わないので、この方法は厳密な処理を行なった場合
の結果と一致する）。この時、逆共分散行列Ｅ［ｈ(k）
ｈ(k)^T］^-1は図４に示すような０部分と非０部分を持つ
行列となる。

【００７６】この図４に示す行列において０部分の帯の
幅Ｗは、エンジン回転数が低下するに従って広がってゆ
くのであるが、さらにこのＷがインパルス応答の有効デ
ータ長よりも大きくなった時、勾配ベクトルｇ(k) はや
はり図５に示すような０部分と非０部分を持つベクトル
となる。ただし、図５において、どこに０部分ができる
かは、ｈの有効データがｈ中のどこに位置するかによ
る。このベクトルでは、エンジン回転数がさらに低下す
るに従って０部分の長さが増してゆくだけである。

【００７７】以上の事柄をまとめれば、インパルス応答
の有効データ長が出力波形データ長に比べて十分に短い
場合には、勾配ベクトルｇ(k) はエンジン回転数に依存
せず一定である（エンジン回転が低下した時、一定の個
所に０を挿入してゆくだけで済む）が、そうした好条件
が得られない場合には、勾配ベクトルｇ(k) はエンジン
回転数の変動に伴って要素毎に変化する。

【００７８】振動低減制御（ＡＮＣ制御）の実際 次に、図６以下を参照しつつ、前述した振動低減の制御
つまりＡＮＣ制御を具体的にどのように行なうかについ
て説明する。先ず、図６は、ニュ−トン法における２次
勾配ベクトルについて、エンジン回転数に応じたデ−タ
数（第１振動の周期をサンプリング周期で除した整数）
と、デ−タ値の全体傾向を示す。全てのエンジン回転域
について常にニュ−トン法を適用する場合（最急降下法
は用いない場合）は、全てのエンジン回転数域に渡って
２次勾配ベクトルを記憶してもよいが、図６に示すよう
なデ−タ値の傾向からして、図６ハッチングを付した領
域となる所定回転数以上の領域については２次勾配ベク
トルを全て記憶する一方、該所定回転数未満の領域では
２次勾配ベクトルそのものを記憶することなく、記憶さ
れているデ−タ値を利用して（オフセットしたり０を付
け加える）、記憶されていない回転域での２次勾配ベク
トルを決定するようにしてもよい。このようにすること
によって、記憶量の低減が図られる。

【００７９】エンジン回転数に応じて、最急降下法とニ
ュ−トン法とを使い分けることができる。すなわち、図
６に対応した図７に示すように、２次勾配ベクトルの記
憶が要求されるハッチングを付した領域のうち、もっと
も低回転域についての２次勾配ベクトルのみを記憶し、
この記憶された２次勾配ベクトルに対応したエンジン回
転数よりも高回転域では最急降下法を用い、記憶された
２次勾配ベクトルに対応したエンジン回転数以下の回転
域ではニュ−トン法を用いる。そして、ニュ−トン法を
用いる場合、２次勾配ベクトルが記憶されてない回転域
での２次勾配ベクトルは、記憶されている２次勾配ベク
トルのデ−タ値を利用して決定される。

【００８０】図８は、最急降下法とニュ−トン法とを上
述のようにエンジン回転数に応じて使い分ける場合の制
御ブロック図を示す。すなわち、適応アルゴリズムを選
択する選択部１１によって、エンジン回転数に応じてニ
ュ−トン法を用いるか最急降下法を用いるかの選択がな
される。また、図７に示すようにして所定の回転域につ
いてのみ２次勾配ベクトルを記憶したメモリＳ１３と、
２次勾配ベクトルを選択するための２次勾配ベクトル選
択部Ｓ１２とが設けられる。２次勾配ベクトル選択部Ｓ
１２は、ニュ−トン法を用いたＡＮＣ制御を行なう場合
に、メモリＳ１３から現在のエンジン回転数に対応し
て、２次勾配ベクトルの選択（記憶値そのものの選択、
あるいは記憶されている２次勾配ベクトルを利用した２
次勾配ベクトルの調整）を行なう。勿論、Ｓ１２で選択
された２次勾配ベクトルに関するデ−タは、ベクトルｙ
の最適化を行なう制御部Ｓ３へ転送される。なお、全て
の回転域について常にニュ−トン法を用いる場合は、Ｓ
１１は不用である。

【００８１】図９、図１０は、２次勾配ベクトルを逆共
分散関数の形式で記憶しておくときの好ましい記憶の仕
方を示し、ハッチングを付した領域のみ２次勾配ベクト
ル（のデ−タ値）が記憶されて、これ以外の領域に関す
る２次勾配ベクトルは、記憶されているデ−タ値に基づ
いて、デ−タ値のオフセットや０値追加によって決定さ
れる。なお、図９、図１０のいずれの場合も、常にニュ
−トン法を用いてもよく、あるいは最急降下法とニュ−
トン法とをエンジン回転数に応じて使い分けるようにし
てもよい。

【００８２】以上実施例について説明したが、低減用振
動発生源としては、スピ−カに限らず、エンジンと車体
との間に介在される容量可変式のアクチュエ−タとする
こともできる。また、周期的な第１振動の発生源として
は、エンジンに限らず、車両に塔載された周期的振動を
発生する適宜の機器を対象とすることができる。勿論、
振動低減されるべき所定空間としては、車室に限らず、
エンジンル−ムや排気経路等適宜選択できるものであ
る。

【図面の簡単な説明】

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

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

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

【図４】逆共分散行列の全体の傾向を示す図。

【図５】図４に示す傾向が強調される状態のときを示す
図。

【図６】エンジン回転数に応じた２次勾配ベクトルのデ
−タ値の傾向と記憶が必要な回転域とを示す図。

【図７】２次勾配ベクトルのデ−タ値の記憶領域の別の
設定変例、およびニュ−トン法と最急降下法との使い分
け領域の設定例を示すもので、図６に対応した図。

【図８】ニュ−トン法と最急降下法とを使い分ける場合
の制御系統をブロック図的に示す図。

【図９】逆共分散関数の形式での好ましい２次勾配ベク
トルの記憶例を示す図。

【図１０】逆共分散関数の形式での他の好ましい２次勾
配ベクトルの記憶例を示す図。

【符号の説明】

１：自動車 ２：車室（所定空間） ８：エンジン（第１振動源） １１：スピ−カ（第２振動源） １２：マイク（振動検出手段） ２０：低減用振動生成回路 Ｓ３：第２振動最適化のための制御部 Ｓ４：第１振動の周期（エンジン回転数）計測部 Ｓ１１：適応アルゴリズム選択部 Ｓ１２：２次勾配ベクトル選択部 Ｓ１３：２次勾配ベクトルを記憶したメモリ

フロントページの続き (56)参考文献 特開 平４−347559（ＪＰ，Ａ) 特開 平１−177604（ＪＰ，Ａ) 特開 平５−11783（ＪＰ，Ａ) 特開 平６−348287（ＪＰ，Ａ) 特表 平１−501344（ＪＰ，Ａ) 浜田晴夫，外２名，同期式適応フィル タとそのアクティブ・ノイズコントロー ルへの応用，第２回環境工学総合シンポ ジウム’92講演論文集，社団法人日本機 械学会，1992年 ７月15日，ｐ．46−50 (58)調査した分野(Int.Cl.⁷，ＤＢ名) G10K 11/178 G10K 11/16 B60R 11/02 F16F 15/02

Claims (5)

(57)【特許請求の範囲】
1. 【請求項１】第１振動源により発生される周期的な第１
   振動を車両の所定空間において低減させる車両の振動低
   減装置であって、 前記第１振動の周期を検出する周期検出手段と、 前記第１振動の振動エネルギを低減させる低減用振動と
   しての第２振動を出力するための第２振動源と、 前記所定空間の振動を検出する振動検出手段と、 前記第２振動源から出力させる第２振動の振動エネルギ
   を、前記周期検出手段で検出される１周期毎に設定する
   設定手段と、 前記設定手段の出力を、前記振動検出手段の出力、およ
   び該振動検出手段と前記第２振動源との間の伝達特性に
   基づいて補正する補正手段と、 ニュ−トン法における２次勾配ベクトルを記憶した記憶
   手段と、 を備え、前記補正手段における補正のための適応アルゴ
   リズムとして、前記第１振動の周期が大きい時に前記記
   憶手段に記憶されている２次勾配ベクトルを利用したニ
   ュ−トン法が設定されると共に、該第１振動の周期が小
   さい時には最急降下法が設定される、 ことを特徴とする車両の振動低減装置。
2. 【請求項２】請求項１において、 前記第１振動源がエンジンとされて、エンジンの運転に
   起因して生じる周期的な第１振動を第２振動で低減す
   る、ことを特徴とする車両の振動低減装置。
3. 【請求項３】請求項２において、 前記記憶手段に、エンジン高回転時用の２次勾配ベクト
   ルのみが記憶され、 エンジン低回転時用の２次勾配ベクトルが、前記記憶さ
   れている２次勾配ベクトルに基づいて設定される、こと
   を特徴とする車両の振動低減装置。
4. 【請求項４】請求項２において、 ２次勾配ベクトルが、逆共分散関数の形式で、かつエン
   ジン回転数に応じて変化される値についてのみ前記記憶
   手段に記憶され、 前記記憶手段に記憶されていない２次勾配ベクトルの値
   は、前記記憶されている２次勾配ベクトルに基づいて設
   定される、 ことを特徴とする車両の振動低減装置 。
5. 【請求項５】請求項２において、 ２次勾配ベクトルが、逆共分散関数の形式で、かつエン
   ジン高回転時用についてのみ前記記憶手段に記憶され、 エンジン低回転時用の２次勾配ベクトルは、前記記憶さ
   れている２次勾配ベクトルに基づいて設定される、 ことを特徴とする車両の振動低減装置 。