能動型振動騒音制御装置
【課題】周波数特性が変化する振動騒音に対してその変化に追従した相殺制御を実行可能な能動型振動騒音制御装置を提供する。
【解決手段】基準信号が入力され、第２制御信号を出力する第２適応ノッチフィルタのフィルタ係数Ｗ２（Ｒｗ２、Ｉｗ２）の複素平面上での位相角度θと、前回の更新の際に算出した前回位相角度θｏｌｄとの間の位相角度変化量ｄθを算出し、位相角度変化量ｄθに応じて基準信号の対象周波数Ｆｃを切り替える。
【解決手段】基準信号が入力され、第２制御信号を出力する第２適応ノッチフィルタのフィルタ係数Ｗ２（Ｒｗ２、Ｉｗ２）の複素平面上での位相角度θと、前回の更新の際に算出した前回位相角度θｏｌｄとの間の位相角度変化量ｄθを算出し、位相角度変化量ｄθに応じて基準信号の対象周波数Ｆｃを切り替える。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、例えば、車両の走行の際に発生する車室内の振動騒音を、相殺振動騒音の出力によって相殺する能動型振動騒音制御装置に関する。
【背景技術】
【０００２】
車両の走行の際、例えば、車輪の振動がサスペンションを介して車体側に伝達することで、ロードノイズ（振動、騒音を含む。以下、総称して「振動騒音」ともいう。）が車室内に発生する。そこで、この振動騒音と逆位相である相殺振動騒音をスピーカから出力し、前記振動騒音を相殺する能動型振動騒音制御装置が種々提案されている。
【０００３】
例えば、特許文献１では、車室内に設けられたマイクロフォンから得た誤差信号から、適応ノッチフィルタを用いて所定周波数の信号成分を抽出し、前記信号成分に基づき生成された制御信号の振幅及び位相を調整する能動型振動騒音制御装置が提案されている。これにより、振動騒音を相殺するための演算処理量を大幅に低減可能であり、当該装置の製造コストを抑制できる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００９−４５９５４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明は上記した特許文献１に開示されている技術的思想に関連してなされたものであって、周波数特性が変化する振動騒音に対してその変化に追従した相殺制御を実行可能な能動型振動騒音制御装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明に係る能動型振動騒音制御装置は、振動騒音に対する相殺信号に基づく相殺振動騒音を出力する振動騒音相殺手段と、前記振動騒音と前記相殺振動騒音との干渉による残留振動騒音を誤差信号として検出する誤差信号検出手段と、前記誤差信号が入力され、前記相殺信号を生成する能動型振動騒音制御手段とを有する装置であって、前記能動型振動騒音制御手段は、所定周波数の基準信号を生成する基準信号生成手段と、前記基準信号が入力され、前記相殺信号の生成に供される制御信号を出力する適応ノッチフィルタと、前記基準信号の周波数に応じた振幅又は位相の調整値を格納し、前記制御信号の振幅又は位相を調整することで前記相殺信号を生成する振幅位相調整手段と、前記誤差信号から前記制御信号を減算して補正誤差信号を生成する補正誤差信号生成手段と、前記基準信号と前記補正誤差信号とに基づいて、前記補正誤差信号が最小となるように前記適応ノッチフィルタのフィルタ係数を逐次更新するフィルタ係数更新手段と、前記振動騒音相殺手段から前記誤差信号検出手段までの伝達特性に基づいて前記相殺信号を補正して補正相殺信号を生成する補正相殺信号生成手段と、前記誤差信号から前記補正相殺信号を減算することで、非制御時検出信号を生成する非制御時検出信号生成手段と、複素平面上で定義されたフィルタ係数を備えており、前記基準信号が入力され、第２制御信号を出力する第２適応ノッチフィルタと、前記非制御時検出信号から前記第２制御信号を減算して第２補正誤差信号を生成する第２補正誤差信号生成手段と、前記基準信号と前記第２補正誤差信号とに基づいて、前記第２補正誤差信号が最小となるように前記第２適応ノッチフィルタのフィルタ係数を逐次更新する第２フィルタ係数更新手段と、前記第２適応ノッチフィルタのフィルタ係数の前記複素平面上での位相角度と、前回の更新の際に算出した位相角度との間の位相角度変化量を算出し、前記位相角度変化量に応じて前記基準信号の周波数を切り替える周波数切替手段とを備えることを特徴とする。
【０００７】
このように、第２適応ノッチフィルタのフィルタ係数の複素平面上での位相角度と、前回の更新の際に算出した位相角度との間の位相角度変化量を算出し、前記位相角度変化量に応じて基準信号の周波数を切り替える周波数切替手段を設けたので、前記フィルタ係数の複素平面上での位相角度の変化量を遂次監視可能であり、この変化量から周波数特性の変化の動向を簡便に且つ精度良く把握できる。これにより、振動騒音の周波数特性に変化があっても、その変化に追従した前記振動騒音の相殺制御を実行できる。
【０００８】
また、前記周波数切替手段は、前記誤差信号のサンプリング周期と、前記位相角度変化量とに基づいて周波数変化量を算出し、前記周波数変化量が下限闘値を下回った場合に前記基準信号の周波数を維持することが好ましい。これにより、周波数変化量が下限閾値を下回った場合、周波数の切り替えに伴う別異のノイズの発生を抑制できる。
【０００９】
さらに、前記周波数切替手段は、前記周波数変化量が、前記下限閾値よりも大きな上限闘値を上回った場合に前記基準信号の周波数を維持することが好ましい。これにより、周波数変化量が上限閾値を上回った場合、過剰な制御による別異のノイズの発生を抑制できる。
【００１０】
さらに、前記周波数切替手段が前記基準信号の周波数を切り替えたことに応じて、前記振幅位相調整手段が格納する前記調整値を切り替える振幅位相切替手段をさらに備えることが好ましい。これにより、基準信号の周波数を切り替えた状態を相殺信号に即時に反映でき、制御の追従性がさらに向上する。
【発明の効果】
【００１１】
本発明に係る能動型振動騒音制御装置によれば、第２適応ノッチフィルタのフィルタ係数の複素平面上での位相角度と、前回の更新の際に算出した位相角度との間の位相角度変化量を算出し、前記位相角度変化量に応じて基準信号の周波数を切り替える周波数切替手段を設けたので、前記フィルタ係数の複素平面上での位相角度の変化量を遂次監視可能であり、この変化量から周波数特性の変化の動向を簡便に且つ精度良く把握できる。これにより、周波数特性が変化する振動騒音に対してその変化に追従した相殺制御を実行できる。
【図面の簡単な説明】
【００１２】
【図１】本実施の形態に係る能動型振動騒音制御装置の概略ブロック図である。
【図２】図１に示すＡＮＣ装置の動作説明に供されるフローチャートである。
【図３】図１に示す能動型振動騒音制御部の第１詳細ブロック図である。
【図４】図１に示す能動型振動騒音制御部の第２詳細ブロック図である。
【図５】図２のステップＳ１１における対象周波数の更新方法を説明するための詳細フローチャートである。
【図６】第２適応ノッチフィルタのフィルタ係数の複素空間上での位相角度変化量の算出方法を表す概略説明図である。
【図７】図７Ａは、ＡＮＣ制御の実行前における誤差信号のスペクトラム図である。図７Ｂは、図７Ａに示す誤差信号に適したＳＡＮ型バンドパスフィルタの周波数特性図である。図７Ｃは、図７Ｂに示すＳＡＮ型バンドパスフィルタの周波数特性に対応する、第１適応ノッチフィルタの周波数特性図である。
【図８】図８Ａは、図７Ａに示す周波数特性が変化した誤差信号のスペクトラム図である。図８Ｂは、図８Ａに示す誤差信号に適したＳＡＮ型バンドパスフィルタの周波数特性図である。
【図９】図９Ａは、周波数切替処理を施さない場合での第１適応ノッチフィルタの周波数特性図である。図９Ｂは、周波数切替処理を施さない場合でのＡＮＣ装置の感度関数を表す周波数特性図である。図９Ｃは、周波数切替処理を施さない場合でのＡＮＣ制御の実行後における誤差信号のスペクトラム図である。
【図１０】図１０Ａは、周波数切替処理を施した場合での第１適応ノッチフィルタの周波数特性図である。図１０Ｂは、周波数切替処理を施した場合でのＡＮＣ装置の感度関数を表す周波数特性図である。図１０Ｃは、周波数切替処理を施した場合でのＡＮＣ制御の実行後における誤差信号のスペクトラム図である。
【発明を実施するための形態】
【００１３】
以下、本発明に係る能動型振動騒音制御装置について好適な実施形態を挙げ、添付の図面を参照して説明する。
【００１４】
図１に示すように、能動型振動騒音制御装置｛以下、ＡＮＣ（Adaptive Noise Control）装置１０という。｝は、車両１１に搭載されている。ＡＮＣ装置１０は、能動型振動騒音制御部１４（能動型振動騒音制御手段）と、マイクロフォン１６（誤差信号検出手段）と、スピーカ１８（振動騒音相殺手段）とを備える。
【００１５】
マイクロフォン１６は、車両１１の内外で発生する各種の音を入力する。入力される音には、路面１２から受ける車輪の振動に起因する振動騒音ＮＳや、この振動騒音ＮＳを相殺するための相殺振動騒音ＣＳが含まれる。この場合、マイクロフォン１６は、振動騒音ＮＳと相殺振動騒音ＣＳとの干渉による残留振動騒音を、能動型振動騒音制御部１４への入力信号（以下、誤差信号Ａという。）として検出する。本図例では、マイクロフォン１６は、車両１１の車室１３上方（具体的には、図示しない乗員の受聴点の近傍）に設けられている。
【００１６】
なお、本明細書における「振動騒音」は、弾性体を伝播する弾性波全般をいう。すなわち、可聴音（可聴周波数を有し、空気中を伝播する弾性波）のように狭義の意味に限定されるものでない。例えば、振動を検出する場合、マイクロフォン１６に代替して振動センサ等を用いてもよい。
【００１７】
スピーカ１８は、能動型振動騒音制御部１４からの出力信号（以下、相殺信号Ｂという。）に基づいて相殺振動騒音ＣＳを出力する。具体的には、スピーカ１８は、所定の周波数を主成分とする振動騒音ＮＳに対して逆位相である相殺振動騒音ＣＳを出力することで、波の干渉効果により振動騒音ＮＳの発生程度を抑制させる。本図例では、スピーカ１８は、車室１３内の座席周辺のキックパネルの近傍に設けられている。
【００１８】
能動型振動騒音制御部１４は、入力された誤差信号Ａに対して所定の信号処理を施して相殺信号Ｂを得た後、スピーカ１８を介して相殺振動騒音ＣＳを出力することで、振動騒音ＮＳを能動的に相殺する制御（以下、「ＡＮＣ制御」という。）を行う。能動型振動騒音制御部１４は、マイクロコンピュータ、ＤＳＰ（Digital Signal Processor）等により構成される。ＣＰＵ（中央演算処理装置）が各種信号の入力に基づき、ＲＯＭ等のメモリに記憶されているプログラムを実行することで、各種処理を実現可能である。
【００１９】
能動型振動騒音制御部１４は、任意の周波数帯域の中からＡＮＣ制御の対象である所定周波数（以下、対象周波数Ｆｃという。）を設定する周波数設定部２０と、周波数設定部２０により設定された対象周波数Ｆｃを主成分とする基準信号Ｘを生成する基準信号生成部２２（基準信号生成手段）と、基準信号生成部２２により生成された基準信号Ｘに対し、ＳＡＮ（Single Adaptive Notch）フィルタを施して第１制御信号Ｏ１（制御信号）を得る第１適応ノッチフィルタ２４（適応ノッチフィルタ）とを備える。
【００２０】
能動型振動騒音制御部１４は、マイクロフォン１６を介して入力された誤差信号Ａから、第１適応ノッチフィルタ２４を介して出力された第１制御信号Ｏ１を減算し、第１補正誤差信号Ｅ１（補正誤差信号）を得る減算器２６（補正誤差信号生成手段）と、基準信号生成部２２により生成された基準信号Ｘ、及び減算器２６から出力された第１補正誤差信号Ｅ１に基づいて、第１補正誤差信号Ｅ１が最小になるように第１適応ノッチフィルタ２４のフィルタ係数Ｗ１を遂次更新する第１フィルタ係数更新部２８（フィルタ係数更新手段）と、をさらに備える。なお、第１適応ノッチフィルタ２４及び減算器２６を組み合わせることで、ＳＡＮ型バンドパスフィルタ３０を構成する。すなわち、第１補正誤差信号Ｅ１は、誤差信号Ａに含まれた比較的広い帯域にわたる各周波数成分のうち、対象周波数Ｆｃを中心とする所定幅の周波数成分を除去した信号に相当する。
【００２１】
能動型振動騒音制御部１４は、振幅又は位相の調整値を用いて第１制御信号Ｏ１の振幅又は位相を調整することで相殺信号Ｂを生成する振幅位相調整部３６（振幅位相調整手段）と、対象周波数Ｆｃに応じて前記調整値を切り替える振幅位相切替部３８（振幅位相切替手段）と、スピーカ１８からマイクロフォン１６までの伝達特性に基づいて相殺信号Ｂを補正して補正相殺信号Ａｂを生成する補正相殺信号生成部４０（補正相殺信号生成手段）と、誤差信号Ａから補正相殺信号Ａｂを減算することで非制御時検出信号Ａｗｏｃを生成する減算器４２（非制御時検出信号生成手段）と、をさらに備える。
【００２２】
能動型振動騒音制御部１４は、基準信号生成部２２により生成された基準信号Ｘに対し、ＳＡＮフィルタを施して第２制御信号Ｏ２を得る第２適応ノッチフィルタ４４と、減算器４２を介して入力された非制御時検出信号Ａｗｏｃから、第２適応ノッチフィルタ４４を介して入力された第２制御信号Ｏ２を減算し、第２補正誤差信号Ｅ２を得る減算器４６（第２補正誤差信号生成手段）と、基準信号生成部２２からの基準信号Ｘ及び減算器４６からの第２補正誤差信号Ｅ２に基づいて、第２補正誤差信号Ｅ２が最小になるように第２適応ノッチフィルタ４４のフィルタ係数Ｗ２を遂次更新する第２フィルタ係数更新部４８（第２フィルタ係数更新手段）と、をさらに備える。なお、第２適応ノッチフィルタ４４及び減算器４６を組み合わせることで、ＳＡＮ型バンドパスフィルタ５０を構成する。すなわち、第２補正誤差信号Ｅ２は、非制御時検出信号Ａｗｏｃに含まれた比較的広い帯域にわたる各周波数成分のうち、対象周波数Ｆｃを中心とする所定幅の周波数成分を除去した信号に相当する。
【００２３】
能動型振動騒音制御部１４は、第２フィルタ係数更新部４８により遂次更新されたフィルタ係数Ｗ２を保持するフィルタ係数保持部５２と、フィルタ係数保持部５２から供給されたフィルタ係数Ｗ２に基づいて、対象周波数Ｆｃの更新要否若しくはその更新量を決定する周波数切替部５４（周波数切替手段）と、をさらに備える。
【００２４】
本実施の形態に係るＡＮＣ装置１０は、基本的には以上のように構成される。以下、この装置の動作について、図２のフローチャートを参照しながら説明する。ところで、図１に示す基準信号Ｘ、フィルタ係数Ｗ１、及びフィルタ係数Ｗ２は、複素平面上で定義されており、実部及び虚部の成分をそれぞれ有する。図２のステップＳ１〜Ｓ８について、図３の第１詳細ブロック図を参照しながら、実部成分及び虚部成分の信号処理流れに留意しつつ詳細に説明する。なお、図示の便宜のため、図３は、図１に示す構成要素の一部を省略して表記している。
【００２５】
ステップＳ１において、マイクロフォン１６は、車室１３内における残留振動騒音を検出し、誤差信号Ａとして入力する。誤差信号Ａには、上記した振動騒音ＮＳのみならず、この振動騒音ＮＳを相殺するためスピーカ１８から出力された相殺振動騒音ＣＳが含まれる。
【００２６】
ステップＳ２において、基準信号生成部２２は、対象周波数Ｆｃを主成分とする基準信号Ｘを生成する。基準信号Ｘの生成に先立ち、周波数設定部２０は、ＡＮＣ制御の対象とする周波数（すなわち、対象周波数Ｆｃ）を設定する。周波数設定部２０は、例えば、制御対象範囲を５０Ｈｚ〜３００Ｈｚとし、１Ｈｚ間隔で設定可能に構成されてもよい。その後、周波数設定部２０は、設定された対象周波数Ｆｃに従って基準信号生成部２２を駆動制御する。
【００２７】
基準信号生成部２２は、基準信号Ｘの実部に対応する実部基準信号Ｒｘ｛＝ｃｏｓ（２πＦｃ・ｔ）｝を生成する実部基準信号生成部６０と、基準信号Ｘの虚部に対応する虚部基準信号Ｉｘ｛＝ｓｉｎ（２πＦｃ・ｔ）｝を生成する虚部基準信号生成部６２とを備える。この場合、基準信号Ｘは、時間（ｔ）に対する三角関数、すなわち、Ｘ（ｔ）＝ｅｘｐ（ｉ２πＦｃ・ｔ）で表現される。なお、ｉは虚数単位である。
【００２８】
ステップＳ３において、第１適応ノッチフィルタ２４は、基準信号生成部２２からの基準信号Ｘに基づいて、減算器２６側及び振幅位相調整部３６側に供給する制御信号Ｏを生成する。以下、第１適応ノッチフィルタ２４の具体的な構成及び動作について説明する。
【００２９】
第１適応ノッチフィルタ２４は、実部フィルタ係数Ｒｗ１が可変に設定される第１フィルタ６４と、虚部フィルタ係数Ｉｗ１が可変に設定される第２フィルタ６６と、第１フィルタ６４側の出力信号から第２フィルタ６６側の出力信号を減算する減算器６８とを備える。第１フィルタ６４は、実部基準信号生成部６０側から入力された実部基準信号Ｒｘ（余弦波信号）の振幅成分をＲｗ倍だけ変倍し、減算器６８側に出力する。第２フィルタ６６は、虚部基準信号生成部６２から入力された虚部基準信号Ｉｘ（正弦波信号）の振幅成分をＩｗ倍だけ変倍し、減算器６８側に出力する。その後、減算器６８は、第１フィルタ６４側の出力信号（＝Ｒｗ１・Ｒｘ）から、第２フィルタ６６側の出力信号（＝Ｉｗ１・Ｉｘ）を減算する。その結果、第１適応ノッチフィルタ２４は、第１制御信号Ｏ１（＝Ｒｗ１・Ｒｘ−Ｉｗ１・Ｉｘ）を出力する。
【００３０】
ステップＳ４において、減算器２６は、マイクロフォン１６を介して入力された誤差信号Ａ（ステップＳ１参照）から、第１適応ノッチフィルタ２４を介して入力された第１制御信号Ｏ１（ステップＳ３参照）を減算することで第１補正誤差信号Ｅ１を生成する。この場合、ＳＡＮ型バンドパスフィルタ３０の作用により、対象周波数Ｆｃを中心とする所定幅の周波数成分を除去した第１補正誤差信号Ｅ１が得られる。
【００３１】
ステップＳ５において、第１フィルタ係数更新部２８は、第１適応ノッチフィルタ２４のフィルタ係数Ｗ１を更新させる。以下、第１適応ノッチフィルタ２４の具体的な構成及び動作について説明する。
【００３２】
第１フィルタ係数更新部２８は、フィルタ係数Ｗ１の実部に対応する実部フィルタ係数Ｒｗ１の更新に供される実部乗算器７０及びゲイン調整器７２、並びに、フィルタ係数Ｗ１の虚部に対応する虚部フィルタ係数Ｉｗ１の更新に供される虚部乗算器７４及びゲイン調整器７６を備える。すなわち、本構成例では、第１フィルタ係数更新部２８は、ＬＭＳ（Least Mean Square）アルゴリズムに従ってフィルタ係数Ｗ１、すなわち実部フィルタ係数Ｒｗ１及び虚部フィルタ係数Ｉｗ１をそれぞれ更新させる。なお、更新アルゴリズムは本手法に限られることなく、種々の最適化手法を採用してもよい。
【００３３】
実部乗算器７０は、実部基準信号生成部６０側から入力された実部基準信号Ｒｘに、減算器２６側から入力された第１補正誤差信号Ｅ１を乗算し、ゲイン調整器７２側に出力する。ゲイン調整器７２は、この乗算信号の振幅成分をμ１倍し、第１フィルタ６４側に出力する。ここで、定数μ１は、ステップサイズパラメータに相当する。そして、第１フィルタ６４は、現時点での実部フィルタ係数Ｒｗ１に対し、第１フィルタ係数更新部２８から取得した更新量（＝＋μ１・Ｒｘ・Ｅ１）を加算することで、新たな実部フィルタ係数Ｒｗ１を得る。すなわち、実部フィルタ係数Ｒｗ１は、次の（１）式に従って更新される。
Ｒｗ１←Ｒｗ１＋μ１・Ｒｘ・Ｅ１ ‥（１）
【００３４】
一方、虚部乗算器７４は、虚部基準信号生成部６２側から入力された虚部基準信号Ｉｘに、減算器２６側から入力された第１補正誤差信号Ｅ１を乗算し、ゲイン調整器７６側に出力する。ゲイン調整器７６は、この乗算信号の振幅成分をμ１倍し、位相を反転（πだけ調整した）した後、第２フィルタ６６側に出力する。その後、第２フィルタ６６は、現時点での虚部フィルタ係数Ｉｗ１に対し、第１フィルタ係数更新部２８から取得した更新量（＝−μ１・Ｉｘ・Ｅ１）を加算することで、新たな虚部フィルタ係数Ｉｗ１を得る。すなわち、虚部フィルタ係数Ｉｗ１は、次の（２）式に従って更新される。
Ｉｗ１←Ｉｗ１−μ１・Ｉｘ・Ｅ１ ‥（２）
【００３５】
ステップＳ６において、振幅位相調整部３６は、第１適応ノッチフィルタ２４側から入力された第１制御信号Ｏ１の振幅及び／又は位相を調整することで、相殺信号Ｂを生成する。
【００３６】
振幅位相調整部３６は、振幅を調整するパラメータである第１調整値Ｇｆｂを用いて第１制御信号Ｏ１の振幅を調整する振幅調整器７８と、位相を調整するパラメータである第２調整値θｆｂを用いて第１制御信号Ｏ１の位相を調整する位相調整器８０と、振幅調整器７８に供給する第１調整値Ｇｆｂを格納する第１格納部８２と、位相調整器８０に供給する第２調整値θｆｂを格納する第２格納部８４とを備える。すなわち、第１制御信号Ｏ１は、振幅調整器７８により振幅を調整され、位相調整器８０により位相を調整された後、相殺信号Ｂとしてスピーカ１８に供給される。
【００３７】
なお、三角関数の加法性を考慮すると、第１制御信号Ｏ１の振幅又は位相を調整した演算結果は、実部基準信号Ｒｘ、虚部基準信号Ｉｘの振幅又は位相をそれぞれ調整した後に合成した演算結果に一致する。そこで、振幅位相調整部３６は、基準信号生成部２２から実部基準信号Ｒｘ、虚部基準信号Ｉｘをそれぞれ取得し、これらの信号の振幅又は位相をそれぞれ調整した後に、合成して相殺信号Ｂを生成するようにしてもよい。
【００３８】
また、振幅位相切替部３８は、周波数切替部５４が基準信号Ｘの対象周波数Ｆｃを切り替えた（詳細は後述する。）ことに応じて、振幅位相調整部３６（第１格納部８２、第２格納部８４）が格納する調整値（第１調整値Ｇｆｂ、第２調整値θｆｂ）を切り替えてもよい。これにより、対象周波数ＦｃがＦｃ’に切り替えられた状態を、相殺信号Ｂを介して即時に反映可能であり、制御の追従性がさらに向上する。
【００３９】
ステップＳ７において、補正相殺信号生成部４０は、振幅位相調整部３６側から入力された制御信号Ｂを補正することで補正相殺信号Ａｂを生成する。補正相殺信号生成部４０は、例えばＦＩＲ（Finite Impulse Response）フィルタで構成され、スピーカ１８からマイクロフォン１６までの伝達経路を表す音の伝達特性を模擬したフィルタ特性を有する。これにより、補正相殺信号Ａｂは、マイクロフォン１６の位置での相殺信号を模擬的に再現する。
【００４０】
ステップＳ８において、減算器４２は、マイクロフォン１６を介して入力された誤差信号Ａ（ステップＳ１参照）から、補正相殺信号生成部４０を介して入力された補正相殺信号Ａｂ（ステップＳ７参照）を減算することで、非制御時検出信号Ａｗｏｃを生成する。この非制御時検出信号Ａｗｏｃは、仮に相殺振動騒音ＣＳを出力しなかった（ＡＮＣ制御を実行しなかった）状態下、マイクロフォン１６により検出されるであろう誤差信号Ａ（すなわち、振動騒音ＮＳに由来する信号）を模擬的に再現する。
【００４１】
続いて、図２のステップＳ９〜Ｓ１１について、図４の第２詳細ブロック図を参照しながら、実部成分及び虚部成分の信号処理流れに留意しつつ詳細に説明する。なお、図示の便宜のため、図４は、図１に示す構成要素の一部を省略して表記している。
【００４２】
ステップＳ９において、第２適応ノッチフィルタ４４は、基準信号生成部２２からの基準信号Ｘに基づいて第２制御信号Ｏ２を生成する。以下、第２適応ノッチフィルタ４４の具体的な構成及び動作について説明する。
【００４３】
第２適応ノッチフィルタ４４は、実部フィルタ係数Ｒｗ２が可変に設定される第１フィルタ９０と、虚部フィルタ係数Ｉｗ２が可変に設定される第２フィルタ９２と、第１フィルタ９０側の出力信号から第２フィルタ９２側の出力信号を減算する減算器９４とを備える。第１フィルタ９０は、実部基準信号生成部６０側から入力された実部基準信号Ｒｘ（余弦波信号）の振幅成分をＲｗ倍だけ変倍し、減算器９４側に出力する。第２フィルタ９２は、虚部基準信号生成部６２から入力された虚部基準信号Ｉｘ（正弦波信号）の振幅成分をＩｗ倍だけ変倍し、減算器９４側に出力する。その後、減算器９４は、第１フィルタ９０側の出力信号（＝Ｒｗ２・Ｒｘ）から、第２フィルタ９２側の出力信号（＝Ｉｗ２・Ｉｘ）を減算する。その結果、第２適応ノッチフィルタ４４は、第２制御信号Ｏ２（＝Ｒｗ２・Ｒｘ−Ｉｗ２・Ｉｘ）を出力する。
【００４４】
ステップＳ１０において、減算器４６は、減算器４２を介して入力された非制御時検出信号Ａｗｏｃ（ステップＳ８参照）から、第２適応ノッチフィルタ４４を介して入力された第２制御信号Ｏ２（ステップＳ９参照）を減算することで第２補正誤差信号Ｅ２を生成する。この場合、ＳＡＮ型バンドパスフィルタ５０の作用により、対象周波数Ｆｃを中心とする所定幅の周波数成分を除去した第２補正誤差信号Ｅ２が得られる。
【００４５】
ステップＳ１０において、第２フィルタ係数更新部４８は、第２適応ノッチフィルタ４４のフィルタ係数Ｗ２を更新させる。以下、第２適応ノッチフィルタ４４の具体的な構成及び動作について説明する。
【００４６】
第２フィルタ係数更新部４８は、フィルタ係数Ｗ２の実部に対応する実部フィルタ係数Ｒｗ２の更新に供される実部乗算器９６及びゲイン調整器９８、並びに、フィルタ係数Ｗ２の虚部に対応する虚部フィルタ係数Ｉｗ２の更新に供される虚部乗算器１００及びゲイン調整器１０２を備える。すなわち、本構成例では、第２フィルタ係数更新部４８は、ＬＭＳ（Least Mean Square）アルゴリズムに従ってフィルタ係数Ｗ２、すなわち実部フィルタ係数Ｒｗ２及び虚部フィルタ係数Ｉｗ２をそれぞれ更新させる。なお、更新アルゴリズムは本手法に限られることなく、種々の最適化手法を採用してもよい。
【００４７】
実部乗算器９６は、実部基準信号生成部６０側から入力された実部基準信号Ｒｘに、減算器４６側から入力された第２補正誤差信号Ｅ２を乗算し、ゲイン調整器９８側に出力する。ゲイン調整器９８は、この乗算信号の振幅成分をμ２倍し、第１フィルタ９０側に出力する。ここで、定数μ２は、ステップサイズパラメータに相当する。そして、第１フィルタ９０は、現時点での実部フィルタ係数Ｒｗ２に対し、第２フィルタ係数更新部４８から取得した更新量（＝＋μ２・Ｒｘ・Ｅ２）を加算することで、新たな実部フィルタ係数Ｒｗ２を得る。すなわち、実部フィルタ係数Ｒｗ２は、次の（３）式に従って更新される。
Ｒｗ２←Ｒｗ２＋μ２・Ｒｘ・Ｅ２ ‥（３）
【００４８】
一方、虚部乗算器１００は、虚部基準信号生成部６２側から入力された虚部基準信号Ｉｘに、減算器４６側から入力された第２補正誤差信号Ｅ２を乗算し、ゲイン調整器１０２側に出力する。ゲイン調整器１０２は、この乗算信号の振幅成分をμ２倍し、位相を反転（πだけ調整した）した後、第２フィルタ９２側に出力する。その後、第２フィルタ９２は、現時点での虚部フィルタ係数Ｉｗ２に対し、第２フィルタ係数更新部４８から取得した更新量（＝−μ２・Ｉｘ・Ｅ２）を加算することで、新たな虚部フィルタ係数Ｉｗ２を得る。すなわち、虚部フィルタ係数Ｉｗ２は、次の（４）式に従って更新される。
Ｉｗ２←Ｉｗ２−μ２・Ｉｘ・Ｅ２ ‥（４）
【００４９】
その後、フィルタ係数保持部５２（第１保持部１０４）は、ステップＳ１０において更新された実部フィルタ係数Ｒｗ２を保持する。また、フィルタ係数保持部５２（第２保持部１０６）は、ステップＳ１０において更新された虚部フィルタ係数Ｉｗ２を保持する。
【００５０】
ステップＳ１１において、周波数切替部５４は、ステップＳ２で設定された対象周波数Ｆｃを基準とし、次の対象周波数Ｆｃ’を決定する。本ステップでは、対象周波数Ｆｃを更新する場合（Ｆｃ’≠Ｆｃ）もあるし、更新しない場合（Ｆｃ’＝Ｆｃ）もある。本ステップの演算処理のことを「周波数切替処理」と称する場合がある。
【００５１】
以下、対象周波数Ｆｃの更新方法、換言すれば周波数切替部５４の具体的動作について、図５のフローチャート及び図６の概略説明図を参照しながら説明する。
【００５２】
ステップＳ２１において、周波数切替部５４は、第２適応ノッチフィルタ４４の現時点ｔでのフィルタ係数Ｗ２（ｔ）における、複素空間上での位相角度θ（０≦θ＜２π）を算出する。位相角度θは、具体的には、対象周波数Ｆｃでのフィルタ係数（Ｒｗ２，Ｉｗ２）を用いて、θ＝ｔａｎ−1（Ｉｗ２／Ｒｗ２）で算出される。
【００５３】
ステップＳ２２において、周波数切替部５４は、ステップＳ２１で算出された位相角度θ、及び前回に算出された位相角度（以下、前回位相角度θｏｌｄという。）から位相角度変化量ｄθを算出する。具体的には、次の（５）式で算出される。
ｄθ＝（θ−θｏｌｄ） ｍｏｄ ２π ‥（５）
【００５４】
ここで、前回位相角度θｏｌｄは、直近のフィルタ係数Ｗ２（ｔ−Ｔｓ）における位相角度θに相当する。なお、位相角度変化量ｄθは、位相角度θの差に限られることなく、前回位相角度θｏｌｄとの間の変化の程度を表すパラメータであれば種類を問わない。また、位相角度変化量ｄθは、前回に算出された位相角度のみならず、直近の複数回に算出された位相角度を用いて算出してもよい。
【００５５】
ステップＳ２３において、周波数切替部５４は、前回位相角度θｏｌｄに、ステップＳ２１で算出された位相角度θの値を代入する。この前回位相角度θｏｌｄは、次回のステップＳ２２での演算に用いられる。
【００５６】
ステップＳ２４において、周波数切替部５４は、ステップＳ２２で算出された位相角度変化量ｄθから周波数変化量ｄＦを算出する。具体的には、ｄＦ＝ｄθ／（２πＴｓ）で算出される。なお、Ｔｓは、誤差信号Ａを入力するサンプリング周期（単位：ｓ）に相当する。
【００５７】
ステップＳ２５において、周波数切替部５４は、対象周波数Ｆｃの更新条件を満たすか否かを判別する。周波数切替部５４は、ステップＳ２４で算出された周波数変化量ｄＦが、所定の範囲内に収まっているか否かで判別する。例えば、下限閾値Ｔｈ１として０．０５〜０．２Ｈｚのいずれかの値を選択し、上限閾値Ｔｈ２として１〜３Ｈｚのいずれかの値を選択しておく。
【００５８】
周波数変化量ｄＦが、Ｔｈ１≦｜ｄＦ｜≦Ｔｈ２の関係を満たす場合、新たな対象周波数Ｆｃ’を、更新式Ｆｃ’＝Ｆｃ＋γｄＦに従って決定する（ステップＳ２６）。なお、γは正値（例えば、０＜γ＜１）であり、本制御の追従速度を調整するためのパラメータに相当する。
【００５９】
一方、０≦｜ｄＦ｜＜Ｔｈ１を満たす場合、Ｆｃ’＝Ｆｃとし、対象周波数Ｆｃを更新しないで維持する（ステップＳ２７）。０≦｜ｄＦ｜＜Ｔｈ１を満たす場合、振動騒音ＮＳの周波数特性が安定していると想定される。対象周波数Ｆｃの切り替えを行わないことで、過度の制御による別異のノイズ（例えば、オーバーシュート）の発生を抑制できる。
【００６０】
あるいは、｜ｄＦ｜＞Ｔｈ２を満たす場合、Ｆｃ’＝Ｆｃとし、対象周波数Ｆｃを更新しないで維持する（ステップＳ２７）。｜ｄＦ｜＞Ｔｈ２を満たす場合、騒音信号ＮＳの挙動の予測が困難である場合や、ＡＮＣ装置１０の起動後からの経過時間が十分でない場合等が想定される。対象周波数Ｆｃの切り替えを行わないことで、過度の制御による別異のノイズ（例えば、オーバーシュート）の発生を抑制できる。
【００６１】
このようにして、周波数切替部５４は、所定のサンプリング周期Ｔｓで対象周波数Ｆｃを遂次決定する（ステップＳ１１）。
【００６２】
ステップＳ１２において、スピーカ１８は、振幅位相調整部３６からの相殺信号Ｂに基づいて相殺振動騒音ＣＳを出力する。以下、所定のサンプリング周期ＴｓでステップＳ１〜Ｓ１２を順次繰り返すことで、振動騒音ＮＳの相殺制御が可能になる。
【００６３】
続いて、上記した周波数切替処理を施すことで得られる作用効果について、図７Ａ〜図１０Ｃを参照しながら説明する。図７Ａ〜図１０Ｃはいずれも、横軸が周波数［Ｈｚ］であり、縦軸がゲイン［ｄＢ］（振幅対数）であるグラフを表す。
【００６４】
図７Ａは、ＡＮＣ制御の実行前における誤差信号Ａ、換言すれば振動騒音ＳＮ単体のスペクトラム図である。第１スペクトラムＳＰＣ１は、周波数４５Ｈｚ近傍に１つのピークを有し、周波数７０Ｈｚ近傍に１つのピークを有する。ここでは、ＡＮＣ制御を用いて、スペクトラム強度が最大である周波数７０Ｈｚ近傍のピークを抑制する場合を想定する。
【００６５】
図７Ｂは、図７Ａに示す誤差信号Ａに適したＳＡＮ型バンドパスフィルタ３０の周波数特性図である。周波数設定部２０（図１、図３及び図４参照）は、対象周波数ＦｃをＦｃ＝７０Ｈｚに設定することで、本図例のように周波数７０Ｈｚでのゲインが最大（信号損失が最小レベル）となるフィルタ特性が得られる。これにより、マイクロフォン１６から入力された振動騒音ＮＳのうち、相殺しようとする周波数成分を選択的に抽出できる。
【００６６】
図７Ｃは、図７Ｂに示すＳＡＮ型バンドパスフィルタ３０の周波数特性に対応する、第１適応ノッチフィルタ２４の周波数特性図である。本特性は、図７Ａに示す第１スペクトラムＳＰＣ１に対し、図７Ｂに示すＳＡＮ型バンドパスフィルタ３０のゲイン（単位：ｄＢ）を周波数毎に加算した結果に略一致する。
【００６７】
ところで、車両１１のサスペンション等を構成する各部品間の相互作用によって、共振ノイズの傾向が異なる場合がある。例えば、車両１１の走行状態に依存して、共振周波数が動的に変化する場合もある。
【００６８】
図８Ａに示すように、車両１１の走行中、誤差信号Ａの周波数特性（破線で図示する第１スペクトラムＳＰＣ１）が変化し、共振周波数が７０Ｈｚから６７Ｈｚにシフトしたとする。以下、変化後の周波数特性を、第２スペクトラムＳＰＣ２（実線で図示する。）と称する。
【００６９】
図８Ｂは、図８Ａに示す誤差信号Ａに適したＳＡＮ型バンドパスフィルタ３０の周波数特性図である。図７Ａ及び図７Ｂの場合と同様に、第２スペクトラムＳＰＣ２のピーク周波数６７Ｈｚでのゲインが最大となるフィルタを用いることで、マイクロフォン１６から入力された振動騒音ＮＳのうち、相殺しようとする周波数成分を選択的に抽出できる。
【００７０】
ところが、上述した周波数切替処理を施さない場合、ＳＡＮ型バンドパスフィルタ３０は、図７Ｂに示した周波数特性のままである。この場合、図９Ａに示すように、第１適応ノッチフィルタ２４の周波数特性は、図７Ｃの特性と比べて、６７Ｈｚ近傍のゲインが相対的に小さくなる。その結果、図９Ｂに示すようなＡＮＣ装置１０の感度関数が得られ、図９Ｃに示すような誤差信号Ａのスペクトラム図が得られる。
【００７１】
図９Ｃにおいて、実線のグラフは、第２スペクトラムＳＰＣ２を備える誤差信号Ａに対し、ＡＮＣ制御を実行した後の周波数特性である。また、破線のグラフは、第１スペクトラムＳＰＣ１を備える誤差信号Ａに対し、ＡＮＣ制御を実行した後の周波数特性である。このように、振動騒音ＮＳの共振周波数がシフトして対象周波数Ｆｃから僅かに外れた場合、その共振周波数近傍における振動騒音ＮＳを十分に相殺できない。
【００７２】
これに対して、本実施の形態に係るＡＮＣ装置１０では、能動型振動騒音制御部１４は、共振周波数のシフトに追従して、ＳＡＮ型バンドパスフィルタ３０の通過帯域を動的に変更する。具体的には、周波数切替部５４は、周波数変化量ｄＦ（＝−３Ｈｚ）を算出した後、対象周波数Ｆｃを７０Ｈｚから６７Ｈｚに切り替える。これにより、ＳＡＮ型バンドパスフィルタ３０の周波数特性は、図８Ｂの破線で示す特性から、同図の実線で示す特性に変更される。
【００７３】
すなわち、上述した周波数切替処理を施した場合、図１０Ａに示すように、第１適応ノッチフィルタ２４の周波数特性は、図９Ａの特性と比べて、６７Ｈｚ近傍のゲインが相対的に大きくなる。その結果、図１０Ｂに示すようなＡＮＣ装置１０の感度関数が得られ、図１０Ｃに示すような誤差信号Ａのスペクトラム図が得られる。
【００７４】
図１０Ｃにおいて、実線のグラフは、第２スペクトラムＳＰＣ２を備える誤差信号Ａに対し、ＡＮＣ制御を実行した後の周波数特性である。また、破線のグラフは、第１スペクトラムＳＰＣ１を備える誤差信号Ａに対し、ＡＮＣ制御を実行した後の周波数特性である。このように、振動騒音ＮＳの共振周波数がシフトした場合であっても、そのシフトの前後にわたって略同程度の相殺効果が得られた。
【００７５】
以上のように、第２適応ノッチフィルタ４４のフィルタ係数Ｗ２（Ｒｗ２、Ｉｗ２）の複素平面上での位相角度θと、前回の更新の際に算出した前回位相角度θｏｌｄとの間の位相角度変化量ｄθを算出し、位相角度変化量ｄθに応じて基準信号Ｘの対象周波数Ｆｃを切り替える周波数切替部５４を設けたので、フィルタ係数Ｗ２の複素平面上での位相角度θの変化量を遂次監視可能であり、位相角度変化量ｄθから周波数特性の変化の動向を簡便に且つ精度良く把握できる。これにより、振動騒音ＮＳの周波数特性に変化があっても、その変化に追従した振動騒音ＮＳの相殺制御を実行できる。
【００７６】
また、能動型振動騒音制御部１４は、２つの適応ノッチフィルタ、すなわち、ＡＮＣ制御に供される第１適応ノッチフィルタ２４と、対象周波数Ｆｃの更新制御に供される第２適応ノッチフィルタ４４とを備えている。位相角度θを算出するための第２適応ノッチフィルタ４４を別途設けたので、この算出処理の際、ＡＮＣ制御の状態の影響を受けにくくなる。これにより、位相角度θの計算精度がさらに高くなり、周波数切替制御の効果及び信頼性が高くなる。
【００７７】
なお、この発明は、上述した実施形態に限定されるものではなく、この発明の主旨を逸脱しない範囲で自由に変更できることは勿論である。
【符号の説明】
【００７８】
１０…ＡＮＣ装置 １１…車両
１４…能動型振動騒音制御部 １６…マイクロフォン
１８…スピーカ ２０…周波数設定部
２２…基準信号生成部 ２４…第１適応ノッチフィルタ
２６、４２、４６…減算器 ２８…第１フィルタ係数更新部
３０、５０…ＳＡＮ型バンドパスフィルタ ３６…振幅位相調整部
３８…振幅位相切替部 ４０…補正相殺信号生成部
４４…第２適応ノッチフィルタ ４８…第２フィルタ係数更新部
５４…周波数切替部
【技術分野】
【０００１】
この発明は、例えば、車両の走行の際に発生する車室内の振動騒音を、相殺振動騒音の出力によって相殺する能動型振動騒音制御装置に関する。
【背景技術】
【０００２】
車両の走行の際、例えば、車輪の振動がサスペンションを介して車体側に伝達することで、ロードノイズ（振動、騒音を含む。以下、総称して「振動騒音」ともいう。）が車室内に発生する。そこで、この振動騒音と逆位相である相殺振動騒音をスピーカから出力し、前記振動騒音を相殺する能動型振動騒音制御装置が種々提案されている。
【０００３】
例えば、特許文献１では、車室内に設けられたマイクロフォンから得た誤差信号から、適応ノッチフィルタを用いて所定周波数の信号成分を抽出し、前記信号成分に基づき生成された制御信号の振幅及び位相を調整する能動型振動騒音制御装置が提案されている。これにより、振動騒音を相殺するための演算処理量を大幅に低減可能であり、当該装置の製造コストを抑制できる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００９−４５９５４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明は上記した特許文献１に開示されている技術的思想に関連してなされたものであって、周波数特性が変化する振動騒音に対してその変化に追従した相殺制御を実行可能な能動型振動騒音制御装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明に係る能動型振動騒音制御装置は、振動騒音に対する相殺信号に基づく相殺振動騒音を出力する振動騒音相殺手段と、前記振動騒音と前記相殺振動騒音との干渉による残留振動騒音を誤差信号として検出する誤差信号検出手段と、前記誤差信号が入力され、前記相殺信号を生成する能動型振動騒音制御手段とを有する装置であって、前記能動型振動騒音制御手段は、所定周波数の基準信号を生成する基準信号生成手段と、前記基準信号が入力され、前記相殺信号の生成に供される制御信号を出力する適応ノッチフィルタと、前記基準信号の周波数に応じた振幅又は位相の調整値を格納し、前記制御信号の振幅又は位相を調整することで前記相殺信号を生成する振幅位相調整手段と、前記誤差信号から前記制御信号を減算して補正誤差信号を生成する補正誤差信号生成手段と、前記基準信号と前記補正誤差信号とに基づいて、前記補正誤差信号が最小となるように前記適応ノッチフィルタのフィルタ係数を逐次更新するフィルタ係数更新手段と、前記振動騒音相殺手段から前記誤差信号検出手段までの伝達特性に基づいて前記相殺信号を補正して補正相殺信号を生成する補正相殺信号生成手段と、前記誤差信号から前記補正相殺信号を減算することで、非制御時検出信号を生成する非制御時検出信号生成手段と、複素平面上で定義されたフィルタ係数を備えており、前記基準信号が入力され、第２制御信号を出力する第２適応ノッチフィルタと、前記非制御時検出信号から前記第２制御信号を減算して第２補正誤差信号を生成する第２補正誤差信号生成手段と、前記基準信号と前記第２補正誤差信号とに基づいて、前記第２補正誤差信号が最小となるように前記第２適応ノッチフィルタのフィルタ係数を逐次更新する第２フィルタ係数更新手段と、前記第２適応ノッチフィルタのフィルタ係数の前記複素平面上での位相角度と、前回の更新の際に算出した位相角度との間の位相角度変化量を算出し、前記位相角度変化量に応じて前記基準信号の周波数を切り替える周波数切替手段とを備えることを特徴とする。
【０００７】
このように、第２適応ノッチフィルタのフィルタ係数の複素平面上での位相角度と、前回の更新の際に算出した位相角度との間の位相角度変化量を算出し、前記位相角度変化量に応じて基準信号の周波数を切り替える周波数切替手段を設けたので、前記フィルタ係数の複素平面上での位相角度の変化量を遂次監視可能であり、この変化量から周波数特性の変化の動向を簡便に且つ精度良く把握できる。これにより、振動騒音の周波数特性に変化があっても、その変化に追従した前記振動騒音の相殺制御を実行できる。
【０００８】
また、前記周波数切替手段は、前記誤差信号のサンプリング周期と、前記位相角度変化量とに基づいて周波数変化量を算出し、前記周波数変化量が下限闘値を下回った場合に前記基準信号の周波数を維持することが好ましい。これにより、周波数変化量が下限閾値を下回った場合、周波数の切り替えに伴う別異のノイズの発生を抑制できる。
【０００９】
さらに、前記周波数切替手段は、前記周波数変化量が、前記下限閾値よりも大きな上限闘値を上回った場合に前記基準信号の周波数を維持することが好ましい。これにより、周波数変化量が上限閾値を上回った場合、過剰な制御による別異のノイズの発生を抑制できる。
【００１０】
さらに、前記周波数切替手段が前記基準信号の周波数を切り替えたことに応じて、前記振幅位相調整手段が格納する前記調整値を切り替える振幅位相切替手段をさらに備えることが好ましい。これにより、基準信号の周波数を切り替えた状態を相殺信号に即時に反映でき、制御の追従性がさらに向上する。
【発明の効果】
【００１１】
本発明に係る能動型振動騒音制御装置によれば、第２適応ノッチフィルタのフィルタ係数の複素平面上での位相角度と、前回の更新の際に算出した位相角度との間の位相角度変化量を算出し、前記位相角度変化量に応じて基準信号の周波数を切り替える周波数切替手段を設けたので、前記フィルタ係数の複素平面上での位相角度の変化量を遂次監視可能であり、この変化量から周波数特性の変化の動向を簡便に且つ精度良く把握できる。これにより、周波数特性が変化する振動騒音に対してその変化に追従した相殺制御を実行できる。
【図面の簡単な説明】
【００１２】
【図１】本実施の形態に係る能動型振動騒音制御装置の概略ブロック図である。
【図２】図１に示すＡＮＣ装置の動作説明に供されるフローチャートである。
【図３】図１に示す能動型振動騒音制御部の第１詳細ブロック図である。
【図４】図１に示す能動型振動騒音制御部の第２詳細ブロック図である。
【図５】図２のステップＳ１１における対象周波数の更新方法を説明するための詳細フローチャートである。
【図６】第２適応ノッチフィルタのフィルタ係数の複素空間上での位相角度変化量の算出方法を表す概略説明図である。
【図７】図７Ａは、ＡＮＣ制御の実行前における誤差信号のスペクトラム図である。図７Ｂは、図７Ａに示す誤差信号に適したＳＡＮ型バンドパスフィルタの周波数特性図である。図７Ｃは、図７Ｂに示すＳＡＮ型バンドパスフィルタの周波数特性に対応する、第１適応ノッチフィルタの周波数特性図である。
【図８】図８Ａは、図７Ａに示す周波数特性が変化した誤差信号のスペクトラム図である。図８Ｂは、図８Ａに示す誤差信号に適したＳＡＮ型バンドパスフィルタの周波数特性図である。
【図９】図９Ａは、周波数切替処理を施さない場合での第１適応ノッチフィルタの周波数特性図である。図９Ｂは、周波数切替処理を施さない場合でのＡＮＣ装置の感度関数を表す周波数特性図である。図９Ｃは、周波数切替処理を施さない場合でのＡＮＣ制御の実行後における誤差信号のスペクトラム図である。
【図１０】図１０Ａは、周波数切替処理を施した場合での第１適応ノッチフィルタの周波数特性図である。図１０Ｂは、周波数切替処理を施した場合でのＡＮＣ装置の感度関数を表す周波数特性図である。図１０Ｃは、周波数切替処理を施した場合でのＡＮＣ制御の実行後における誤差信号のスペクトラム図である。
【発明を実施するための形態】
【００１３】
以下、本発明に係る能動型振動騒音制御装置について好適な実施形態を挙げ、添付の図面を参照して説明する。
【００１４】
図１に示すように、能動型振動騒音制御装置｛以下、ＡＮＣ（Adaptive Noise Control）装置１０という。｝は、車両１１に搭載されている。ＡＮＣ装置１０は、能動型振動騒音制御部１４（能動型振動騒音制御手段）と、マイクロフォン１６（誤差信号検出手段）と、スピーカ１８（振動騒音相殺手段）とを備える。
【００１５】
マイクロフォン１６は、車両１１の内外で発生する各種の音を入力する。入力される音には、路面１２から受ける車輪の振動に起因する振動騒音ＮＳや、この振動騒音ＮＳを相殺するための相殺振動騒音ＣＳが含まれる。この場合、マイクロフォン１６は、振動騒音ＮＳと相殺振動騒音ＣＳとの干渉による残留振動騒音を、能動型振動騒音制御部１４への入力信号（以下、誤差信号Ａという。）として検出する。本図例では、マイクロフォン１６は、車両１１の車室１３上方（具体的には、図示しない乗員の受聴点の近傍）に設けられている。
【００１６】
なお、本明細書における「振動騒音」は、弾性体を伝播する弾性波全般をいう。すなわち、可聴音（可聴周波数を有し、空気中を伝播する弾性波）のように狭義の意味に限定されるものでない。例えば、振動を検出する場合、マイクロフォン１６に代替して振動センサ等を用いてもよい。
【００１７】
スピーカ１８は、能動型振動騒音制御部１４からの出力信号（以下、相殺信号Ｂという。）に基づいて相殺振動騒音ＣＳを出力する。具体的には、スピーカ１８は、所定の周波数を主成分とする振動騒音ＮＳに対して逆位相である相殺振動騒音ＣＳを出力することで、波の干渉効果により振動騒音ＮＳの発生程度を抑制させる。本図例では、スピーカ１８は、車室１３内の座席周辺のキックパネルの近傍に設けられている。
【００１８】
能動型振動騒音制御部１４は、入力された誤差信号Ａに対して所定の信号処理を施して相殺信号Ｂを得た後、スピーカ１８を介して相殺振動騒音ＣＳを出力することで、振動騒音ＮＳを能動的に相殺する制御（以下、「ＡＮＣ制御」という。）を行う。能動型振動騒音制御部１４は、マイクロコンピュータ、ＤＳＰ（Digital Signal Processor）等により構成される。ＣＰＵ（中央演算処理装置）が各種信号の入力に基づき、ＲＯＭ等のメモリに記憶されているプログラムを実行することで、各種処理を実現可能である。
【００１９】
能動型振動騒音制御部１４は、任意の周波数帯域の中からＡＮＣ制御の対象である所定周波数（以下、対象周波数Ｆｃという。）を設定する周波数設定部２０と、周波数設定部２０により設定された対象周波数Ｆｃを主成分とする基準信号Ｘを生成する基準信号生成部２２（基準信号生成手段）と、基準信号生成部２２により生成された基準信号Ｘに対し、ＳＡＮ（Single Adaptive Notch）フィルタを施して第１制御信号Ｏ１（制御信号）を得る第１適応ノッチフィルタ２４（適応ノッチフィルタ）とを備える。
【００２０】
能動型振動騒音制御部１４は、マイクロフォン１６を介して入力された誤差信号Ａから、第１適応ノッチフィルタ２４を介して出力された第１制御信号Ｏ１を減算し、第１補正誤差信号Ｅ１（補正誤差信号）を得る減算器２６（補正誤差信号生成手段）と、基準信号生成部２２により生成された基準信号Ｘ、及び減算器２６から出力された第１補正誤差信号Ｅ１に基づいて、第１補正誤差信号Ｅ１が最小になるように第１適応ノッチフィルタ２４のフィルタ係数Ｗ１を遂次更新する第１フィルタ係数更新部２８（フィルタ係数更新手段）と、をさらに備える。なお、第１適応ノッチフィルタ２４及び減算器２６を組み合わせることで、ＳＡＮ型バンドパスフィルタ３０を構成する。すなわち、第１補正誤差信号Ｅ１は、誤差信号Ａに含まれた比較的広い帯域にわたる各周波数成分のうち、対象周波数Ｆｃを中心とする所定幅の周波数成分を除去した信号に相当する。
【００２１】
能動型振動騒音制御部１４は、振幅又は位相の調整値を用いて第１制御信号Ｏ１の振幅又は位相を調整することで相殺信号Ｂを生成する振幅位相調整部３６（振幅位相調整手段）と、対象周波数Ｆｃに応じて前記調整値を切り替える振幅位相切替部３８（振幅位相切替手段）と、スピーカ１８からマイクロフォン１６までの伝達特性に基づいて相殺信号Ｂを補正して補正相殺信号Ａｂを生成する補正相殺信号生成部４０（補正相殺信号生成手段）と、誤差信号Ａから補正相殺信号Ａｂを減算することで非制御時検出信号Ａｗｏｃを生成する減算器４２（非制御時検出信号生成手段）と、をさらに備える。
【００２２】
能動型振動騒音制御部１４は、基準信号生成部２２により生成された基準信号Ｘに対し、ＳＡＮフィルタを施して第２制御信号Ｏ２を得る第２適応ノッチフィルタ４４と、減算器４２を介して入力された非制御時検出信号Ａｗｏｃから、第２適応ノッチフィルタ４４を介して入力された第２制御信号Ｏ２を減算し、第２補正誤差信号Ｅ２を得る減算器４６（第２補正誤差信号生成手段）と、基準信号生成部２２からの基準信号Ｘ及び減算器４６からの第２補正誤差信号Ｅ２に基づいて、第２補正誤差信号Ｅ２が最小になるように第２適応ノッチフィルタ４４のフィルタ係数Ｗ２を遂次更新する第２フィルタ係数更新部４８（第２フィルタ係数更新手段）と、をさらに備える。なお、第２適応ノッチフィルタ４４及び減算器４６を組み合わせることで、ＳＡＮ型バンドパスフィルタ５０を構成する。すなわち、第２補正誤差信号Ｅ２は、非制御時検出信号Ａｗｏｃに含まれた比較的広い帯域にわたる各周波数成分のうち、対象周波数Ｆｃを中心とする所定幅の周波数成分を除去した信号に相当する。
【００２３】
能動型振動騒音制御部１４は、第２フィルタ係数更新部４８により遂次更新されたフィルタ係数Ｗ２を保持するフィルタ係数保持部５２と、フィルタ係数保持部５２から供給されたフィルタ係数Ｗ２に基づいて、対象周波数Ｆｃの更新要否若しくはその更新量を決定する周波数切替部５４（周波数切替手段）と、をさらに備える。
【００２４】
本実施の形態に係るＡＮＣ装置１０は、基本的には以上のように構成される。以下、この装置の動作について、図２のフローチャートを参照しながら説明する。ところで、図１に示す基準信号Ｘ、フィルタ係数Ｗ１、及びフィルタ係数Ｗ２は、複素平面上で定義されており、実部及び虚部の成分をそれぞれ有する。図２のステップＳ１〜Ｓ８について、図３の第１詳細ブロック図を参照しながら、実部成分及び虚部成分の信号処理流れに留意しつつ詳細に説明する。なお、図示の便宜のため、図３は、図１に示す構成要素の一部を省略して表記している。
【００２５】
ステップＳ１において、マイクロフォン１６は、車室１３内における残留振動騒音を検出し、誤差信号Ａとして入力する。誤差信号Ａには、上記した振動騒音ＮＳのみならず、この振動騒音ＮＳを相殺するためスピーカ１８から出力された相殺振動騒音ＣＳが含まれる。
【００２６】
ステップＳ２において、基準信号生成部２２は、対象周波数Ｆｃを主成分とする基準信号Ｘを生成する。基準信号Ｘの生成に先立ち、周波数設定部２０は、ＡＮＣ制御の対象とする周波数（すなわち、対象周波数Ｆｃ）を設定する。周波数設定部２０は、例えば、制御対象範囲を５０Ｈｚ〜３００Ｈｚとし、１Ｈｚ間隔で設定可能に構成されてもよい。その後、周波数設定部２０は、設定された対象周波数Ｆｃに従って基準信号生成部２２を駆動制御する。
【００２７】
基準信号生成部２２は、基準信号Ｘの実部に対応する実部基準信号Ｒｘ｛＝ｃｏｓ（２πＦｃ・ｔ）｝を生成する実部基準信号生成部６０と、基準信号Ｘの虚部に対応する虚部基準信号Ｉｘ｛＝ｓｉｎ（２πＦｃ・ｔ）｝を生成する虚部基準信号生成部６２とを備える。この場合、基準信号Ｘは、時間（ｔ）に対する三角関数、すなわち、Ｘ（ｔ）＝ｅｘｐ（ｉ２πＦｃ・ｔ）で表現される。なお、ｉは虚数単位である。
【００２８】
ステップＳ３において、第１適応ノッチフィルタ２４は、基準信号生成部２２からの基準信号Ｘに基づいて、減算器２６側及び振幅位相調整部３６側に供給する制御信号Ｏを生成する。以下、第１適応ノッチフィルタ２４の具体的な構成及び動作について説明する。
【００２９】
第１適応ノッチフィルタ２４は、実部フィルタ係数Ｒｗ１が可変に設定される第１フィルタ６４と、虚部フィルタ係数Ｉｗ１が可変に設定される第２フィルタ６６と、第１フィルタ６４側の出力信号から第２フィルタ６６側の出力信号を減算する減算器６８とを備える。第１フィルタ６４は、実部基準信号生成部６０側から入力された実部基準信号Ｒｘ（余弦波信号）の振幅成分をＲｗ倍だけ変倍し、減算器６８側に出力する。第２フィルタ６６は、虚部基準信号生成部６２から入力された虚部基準信号Ｉｘ（正弦波信号）の振幅成分をＩｗ倍だけ変倍し、減算器６８側に出力する。その後、減算器６８は、第１フィルタ６４側の出力信号（＝Ｒｗ１・Ｒｘ）から、第２フィルタ６６側の出力信号（＝Ｉｗ１・Ｉｘ）を減算する。その結果、第１適応ノッチフィルタ２４は、第１制御信号Ｏ１（＝Ｒｗ１・Ｒｘ−Ｉｗ１・Ｉｘ）を出力する。
【００３０】
ステップＳ４において、減算器２６は、マイクロフォン１６を介して入力された誤差信号Ａ（ステップＳ１参照）から、第１適応ノッチフィルタ２４を介して入力された第１制御信号Ｏ１（ステップＳ３参照）を減算することで第１補正誤差信号Ｅ１を生成する。この場合、ＳＡＮ型バンドパスフィルタ３０の作用により、対象周波数Ｆｃを中心とする所定幅の周波数成分を除去した第１補正誤差信号Ｅ１が得られる。
【００３１】
ステップＳ５において、第１フィルタ係数更新部２８は、第１適応ノッチフィルタ２４のフィルタ係数Ｗ１を更新させる。以下、第１適応ノッチフィルタ２４の具体的な構成及び動作について説明する。
【００３２】
第１フィルタ係数更新部２８は、フィルタ係数Ｗ１の実部に対応する実部フィルタ係数Ｒｗ１の更新に供される実部乗算器７０及びゲイン調整器７２、並びに、フィルタ係数Ｗ１の虚部に対応する虚部フィルタ係数Ｉｗ１の更新に供される虚部乗算器７４及びゲイン調整器７６を備える。すなわち、本構成例では、第１フィルタ係数更新部２８は、ＬＭＳ（Least Mean Square）アルゴリズムに従ってフィルタ係数Ｗ１、すなわち実部フィルタ係数Ｒｗ１及び虚部フィルタ係数Ｉｗ１をそれぞれ更新させる。なお、更新アルゴリズムは本手法に限られることなく、種々の最適化手法を採用してもよい。
【００３３】
実部乗算器７０は、実部基準信号生成部６０側から入力された実部基準信号Ｒｘに、減算器２６側から入力された第１補正誤差信号Ｅ１を乗算し、ゲイン調整器７２側に出力する。ゲイン調整器７２は、この乗算信号の振幅成分をμ１倍し、第１フィルタ６４側に出力する。ここで、定数μ１は、ステップサイズパラメータに相当する。そして、第１フィルタ６４は、現時点での実部フィルタ係数Ｒｗ１に対し、第１フィルタ係数更新部２８から取得した更新量（＝＋μ１・Ｒｘ・Ｅ１）を加算することで、新たな実部フィルタ係数Ｒｗ１を得る。すなわち、実部フィルタ係数Ｒｗ１は、次の（１）式に従って更新される。
Ｒｗ１←Ｒｗ１＋μ１・Ｒｘ・Ｅ１ ‥（１）
【００３４】
一方、虚部乗算器７４は、虚部基準信号生成部６２側から入力された虚部基準信号Ｉｘに、減算器２６側から入力された第１補正誤差信号Ｅ１を乗算し、ゲイン調整器７６側に出力する。ゲイン調整器７６は、この乗算信号の振幅成分をμ１倍し、位相を反転（πだけ調整した）した後、第２フィルタ６６側に出力する。その後、第２フィルタ６６は、現時点での虚部フィルタ係数Ｉｗ１に対し、第１フィルタ係数更新部２８から取得した更新量（＝−μ１・Ｉｘ・Ｅ１）を加算することで、新たな虚部フィルタ係数Ｉｗ１を得る。すなわち、虚部フィルタ係数Ｉｗ１は、次の（２）式に従って更新される。
Ｉｗ１←Ｉｗ１−μ１・Ｉｘ・Ｅ１ ‥（２）
【００３５】
ステップＳ６において、振幅位相調整部３６は、第１適応ノッチフィルタ２４側から入力された第１制御信号Ｏ１の振幅及び／又は位相を調整することで、相殺信号Ｂを生成する。
【００３６】
振幅位相調整部３６は、振幅を調整するパラメータである第１調整値Ｇｆｂを用いて第１制御信号Ｏ１の振幅を調整する振幅調整器７８と、位相を調整するパラメータである第２調整値θｆｂを用いて第１制御信号Ｏ１の位相を調整する位相調整器８０と、振幅調整器７８に供給する第１調整値Ｇｆｂを格納する第１格納部８２と、位相調整器８０に供給する第２調整値θｆｂを格納する第２格納部８４とを備える。すなわち、第１制御信号Ｏ１は、振幅調整器７８により振幅を調整され、位相調整器８０により位相を調整された後、相殺信号Ｂとしてスピーカ１８に供給される。
【００３７】
なお、三角関数の加法性を考慮すると、第１制御信号Ｏ１の振幅又は位相を調整した演算結果は、実部基準信号Ｒｘ、虚部基準信号Ｉｘの振幅又は位相をそれぞれ調整した後に合成した演算結果に一致する。そこで、振幅位相調整部３６は、基準信号生成部２２から実部基準信号Ｒｘ、虚部基準信号Ｉｘをそれぞれ取得し、これらの信号の振幅又は位相をそれぞれ調整した後に、合成して相殺信号Ｂを生成するようにしてもよい。
【００３８】
また、振幅位相切替部３８は、周波数切替部５４が基準信号Ｘの対象周波数Ｆｃを切り替えた（詳細は後述する。）ことに応じて、振幅位相調整部３６（第１格納部８２、第２格納部８４）が格納する調整値（第１調整値Ｇｆｂ、第２調整値θｆｂ）を切り替えてもよい。これにより、対象周波数ＦｃがＦｃ’に切り替えられた状態を、相殺信号Ｂを介して即時に反映可能であり、制御の追従性がさらに向上する。
【００３９】
ステップＳ７において、補正相殺信号生成部４０は、振幅位相調整部３６側から入力された制御信号Ｂを補正することで補正相殺信号Ａｂを生成する。補正相殺信号生成部４０は、例えばＦＩＲ（Finite Impulse Response）フィルタで構成され、スピーカ１８からマイクロフォン１６までの伝達経路を表す音の伝達特性を模擬したフィルタ特性を有する。これにより、補正相殺信号Ａｂは、マイクロフォン１６の位置での相殺信号を模擬的に再現する。
【００４０】
ステップＳ８において、減算器４２は、マイクロフォン１６を介して入力された誤差信号Ａ（ステップＳ１参照）から、補正相殺信号生成部４０を介して入力された補正相殺信号Ａｂ（ステップＳ７参照）を減算することで、非制御時検出信号Ａｗｏｃを生成する。この非制御時検出信号Ａｗｏｃは、仮に相殺振動騒音ＣＳを出力しなかった（ＡＮＣ制御を実行しなかった）状態下、マイクロフォン１６により検出されるであろう誤差信号Ａ（すなわち、振動騒音ＮＳに由来する信号）を模擬的に再現する。
【００４１】
続いて、図２のステップＳ９〜Ｓ１１について、図４の第２詳細ブロック図を参照しながら、実部成分及び虚部成分の信号処理流れに留意しつつ詳細に説明する。なお、図示の便宜のため、図４は、図１に示す構成要素の一部を省略して表記している。
【００４２】
ステップＳ９において、第２適応ノッチフィルタ４４は、基準信号生成部２２からの基準信号Ｘに基づいて第２制御信号Ｏ２を生成する。以下、第２適応ノッチフィルタ４４の具体的な構成及び動作について説明する。
【００４３】
第２適応ノッチフィルタ４４は、実部フィルタ係数Ｒｗ２が可変に設定される第１フィルタ９０と、虚部フィルタ係数Ｉｗ２が可変に設定される第２フィルタ９２と、第１フィルタ９０側の出力信号から第２フィルタ９２側の出力信号を減算する減算器９４とを備える。第１フィルタ９０は、実部基準信号生成部６０側から入力された実部基準信号Ｒｘ（余弦波信号）の振幅成分をＲｗ倍だけ変倍し、減算器９４側に出力する。第２フィルタ９２は、虚部基準信号生成部６２から入力された虚部基準信号Ｉｘ（正弦波信号）の振幅成分をＩｗ倍だけ変倍し、減算器９４側に出力する。その後、減算器９４は、第１フィルタ９０側の出力信号（＝Ｒｗ２・Ｒｘ）から、第２フィルタ９２側の出力信号（＝Ｉｗ２・Ｉｘ）を減算する。その結果、第２適応ノッチフィルタ４４は、第２制御信号Ｏ２（＝Ｒｗ２・Ｒｘ−Ｉｗ２・Ｉｘ）を出力する。
【００４４】
ステップＳ１０において、減算器４６は、減算器４２を介して入力された非制御時検出信号Ａｗｏｃ（ステップＳ８参照）から、第２適応ノッチフィルタ４４を介して入力された第２制御信号Ｏ２（ステップＳ９参照）を減算することで第２補正誤差信号Ｅ２を生成する。この場合、ＳＡＮ型バンドパスフィルタ５０の作用により、対象周波数Ｆｃを中心とする所定幅の周波数成分を除去した第２補正誤差信号Ｅ２が得られる。
【００４５】
ステップＳ１０において、第２フィルタ係数更新部４８は、第２適応ノッチフィルタ４４のフィルタ係数Ｗ２を更新させる。以下、第２適応ノッチフィルタ４４の具体的な構成及び動作について説明する。
【００４６】
第２フィルタ係数更新部４８は、フィルタ係数Ｗ２の実部に対応する実部フィルタ係数Ｒｗ２の更新に供される実部乗算器９６及びゲイン調整器９８、並びに、フィルタ係数Ｗ２の虚部に対応する虚部フィルタ係数Ｉｗ２の更新に供される虚部乗算器１００及びゲイン調整器１０２を備える。すなわち、本構成例では、第２フィルタ係数更新部４８は、ＬＭＳ（Least Mean Square）アルゴリズムに従ってフィルタ係数Ｗ２、すなわち実部フィルタ係数Ｒｗ２及び虚部フィルタ係数Ｉｗ２をそれぞれ更新させる。なお、更新アルゴリズムは本手法に限られることなく、種々の最適化手法を採用してもよい。
【００４７】
実部乗算器９６は、実部基準信号生成部６０側から入力された実部基準信号Ｒｘに、減算器４６側から入力された第２補正誤差信号Ｅ２を乗算し、ゲイン調整器９８側に出力する。ゲイン調整器９８は、この乗算信号の振幅成分をμ２倍し、第１フィルタ９０側に出力する。ここで、定数μ２は、ステップサイズパラメータに相当する。そして、第１フィルタ９０は、現時点での実部フィルタ係数Ｒｗ２に対し、第２フィルタ係数更新部４８から取得した更新量（＝＋μ２・Ｒｘ・Ｅ２）を加算することで、新たな実部フィルタ係数Ｒｗ２を得る。すなわち、実部フィルタ係数Ｒｗ２は、次の（３）式に従って更新される。
Ｒｗ２←Ｒｗ２＋μ２・Ｒｘ・Ｅ２ ‥（３）
【００４８】
一方、虚部乗算器１００は、虚部基準信号生成部６２側から入力された虚部基準信号Ｉｘに、減算器４６側から入力された第２補正誤差信号Ｅ２を乗算し、ゲイン調整器１０２側に出力する。ゲイン調整器１０２は、この乗算信号の振幅成分をμ２倍し、位相を反転（πだけ調整した）した後、第２フィルタ９２側に出力する。その後、第２フィルタ９２は、現時点での虚部フィルタ係数Ｉｗ２に対し、第２フィルタ係数更新部４８から取得した更新量（＝−μ２・Ｉｘ・Ｅ２）を加算することで、新たな虚部フィルタ係数Ｉｗ２を得る。すなわち、虚部フィルタ係数Ｉｗ２は、次の（４）式に従って更新される。
Ｉｗ２←Ｉｗ２−μ２・Ｉｘ・Ｅ２ ‥（４）
【００４９】
その後、フィルタ係数保持部５２（第１保持部１０４）は、ステップＳ１０において更新された実部フィルタ係数Ｒｗ２を保持する。また、フィルタ係数保持部５２（第２保持部１０６）は、ステップＳ１０において更新された虚部フィルタ係数Ｉｗ２を保持する。
【００５０】
ステップＳ１１において、周波数切替部５４は、ステップＳ２で設定された対象周波数Ｆｃを基準とし、次の対象周波数Ｆｃ’を決定する。本ステップでは、対象周波数Ｆｃを更新する場合（Ｆｃ’≠Ｆｃ）もあるし、更新しない場合（Ｆｃ’＝Ｆｃ）もある。本ステップの演算処理のことを「周波数切替処理」と称する場合がある。
【００５１】
以下、対象周波数Ｆｃの更新方法、換言すれば周波数切替部５４の具体的動作について、図５のフローチャート及び図６の概略説明図を参照しながら説明する。
【００５２】
ステップＳ２１において、周波数切替部５４は、第２適応ノッチフィルタ４４の現時点ｔでのフィルタ係数Ｗ２（ｔ）における、複素空間上での位相角度θ（０≦θ＜２π）を算出する。位相角度θは、具体的には、対象周波数Ｆｃでのフィルタ係数（Ｒｗ２，Ｉｗ２）を用いて、θ＝ｔａｎ−1（Ｉｗ２／Ｒｗ２）で算出される。
【００５３】
ステップＳ２２において、周波数切替部５４は、ステップＳ２１で算出された位相角度θ、及び前回に算出された位相角度（以下、前回位相角度θｏｌｄという。）から位相角度変化量ｄθを算出する。具体的には、次の（５）式で算出される。
ｄθ＝（θ−θｏｌｄ） ｍｏｄ ２π ‥（５）
【００５４】
ここで、前回位相角度θｏｌｄは、直近のフィルタ係数Ｗ２（ｔ−Ｔｓ）における位相角度θに相当する。なお、位相角度変化量ｄθは、位相角度θの差に限られることなく、前回位相角度θｏｌｄとの間の変化の程度を表すパラメータであれば種類を問わない。また、位相角度変化量ｄθは、前回に算出された位相角度のみならず、直近の複数回に算出された位相角度を用いて算出してもよい。
【００５５】
ステップＳ２３において、周波数切替部５４は、前回位相角度θｏｌｄに、ステップＳ２１で算出された位相角度θの値を代入する。この前回位相角度θｏｌｄは、次回のステップＳ２２での演算に用いられる。
【００５６】
ステップＳ２４において、周波数切替部５４は、ステップＳ２２で算出された位相角度変化量ｄθから周波数変化量ｄＦを算出する。具体的には、ｄＦ＝ｄθ／（２πＴｓ）で算出される。なお、Ｔｓは、誤差信号Ａを入力するサンプリング周期（単位：ｓ）に相当する。
【００５７】
ステップＳ２５において、周波数切替部５４は、対象周波数Ｆｃの更新条件を満たすか否かを判別する。周波数切替部５４は、ステップＳ２４で算出された周波数変化量ｄＦが、所定の範囲内に収まっているか否かで判別する。例えば、下限閾値Ｔｈ１として０．０５〜０．２Ｈｚのいずれかの値を選択し、上限閾値Ｔｈ２として１〜３Ｈｚのいずれかの値を選択しておく。
【００５８】
周波数変化量ｄＦが、Ｔｈ１≦｜ｄＦ｜≦Ｔｈ２の関係を満たす場合、新たな対象周波数Ｆｃ’を、更新式Ｆｃ’＝Ｆｃ＋γｄＦに従って決定する（ステップＳ２６）。なお、γは正値（例えば、０＜γ＜１）であり、本制御の追従速度を調整するためのパラメータに相当する。
【００５９】
一方、０≦｜ｄＦ｜＜Ｔｈ１を満たす場合、Ｆｃ’＝Ｆｃとし、対象周波数Ｆｃを更新しないで維持する（ステップＳ２７）。０≦｜ｄＦ｜＜Ｔｈ１を満たす場合、振動騒音ＮＳの周波数特性が安定していると想定される。対象周波数Ｆｃの切り替えを行わないことで、過度の制御による別異のノイズ（例えば、オーバーシュート）の発生を抑制できる。
【００６０】
あるいは、｜ｄＦ｜＞Ｔｈ２を満たす場合、Ｆｃ’＝Ｆｃとし、対象周波数Ｆｃを更新しないで維持する（ステップＳ２７）。｜ｄＦ｜＞Ｔｈ２を満たす場合、騒音信号ＮＳの挙動の予測が困難である場合や、ＡＮＣ装置１０の起動後からの経過時間が十分でない場合等が想定される。対象周波数Ｆｃの切り替えを行わないことで、過度の制御による別異のノイズ（例えば、オーバーシュート）の発生を抑制できる。
【００６１】
このようにして、周波数切替部５４は、所定のサンプリング周期Ｔｓで対象周波数Ｆｃを遂次決定する（ステップＳ１１）。
【００６２】
ステップＳ１２において、スピーカ１８は、振幅位相調整部３６からの相殺信号Ｂに基づいて相殺振動騒音ＣＳを出力する。以下、所定のサンプリング周期ＴｓでステップＳ１〜Ｓ１２を順次繰り返すことで、振動騒音ＮＳの相殺制御が可能になる。
【００６３】
続いて、上記した周波数切替処理を施すことで得られる作用効果について、図７Ａ〜図１０Ｃを参照しながら説明する。図７Ａ〜図１０Ｃはいずれも、横軸が周波数［Ｈｚ］であり、縦軸がゲイン［ｄＢ］（振幅対数）であるグラフを表す。
【００６４】
図７Ａは、ＡＮＣ制御の実行前における誤差信号Ａ、換言すれば振動騒音ＳＮ単体のスペクトラム図である。第１スペクトラムＳＰＣ１は、周波数４５Ｈｚ近傍に１つのピークを有し、周波数７０Ｈｚ近傍に１つのピークを有する。ここでは、ＡＮＣ制御を用いて、スペクトラム強度が最大である周波数７０Ｈｚ近傍のピークを抑制する場合を想定する。
【００６５】
図７Ｂは、図７Ａに示す誤差信号Ａに適したＳＡＮ型バンドパスフィルタ３０の周波数特性図である。周波数設定部２０（図１、図３及び図４参照）は、対象周波数ＦｃをＦｃ＝７０Ｈｚに設定することで、本図例のように周波数７０Ｈｚでのゲインが最大（信号損失が最小レベル）となるフィルタ特性が得られる。これにより、マイクロフォン１６から入力された振動騒音ＮＳのうち、相殺しようとする周波数成分を選択的に抽出できる。
【００６６】
図７Ｃは、図７Ｂに示すＳＡＮ型バンドパスフィルタ３０の周波数特性に対応する、第１適応ノッチフィルタ２４の周波数特性図である。本特性は、図７Ａに示す第１スペクトラムＳＰＣ１に対し、図７Ｂに示すＳＡＮ型バンドパスフィルタ３０のゲイン（単位：ｄＢ）を周波数毎に加算した結果に略一致する。
【００６７】
ところで、車両１１のサスペンション等を構成する各部品間の相互作用によって、共振ノイズの傾向が異なる場合がある。例えば、車両１１の走行状態に依存して、共振周波数が動的に変化する場合もある。
【００６８】
図８Ａに示すように、車両１１の走行中、誤差信号Ａの周波数特性（破線で図示する第１スペクトラムＳＰＣ１）が変化し、共振周波数が７０Ｈｚから６７Ｈｚにシフトしたとする。以下、変化後の周波数特性を、第２スペクトラムＳＰＣ２（実線で図示する。）と称する。
【００６９】
図８Ｂは、図８Ａに示す誤差信号Ａに適したＳＡＮ型バンドパスフィルタ３０の周波数特性図である。図７Ａ及び図７Ｂの場合と同様に、第２スペクトラムＳＰＣ２のピーク周波数６７Ｈｚでのゲインが最大となるフィルタを用いることで、マイクロフォン１６から入力された振動騒音ＮＳのうち、相殺しようとする周波数成分を選択的に抽出できる。
【００７０】
ところが、上述した周波数切替処理を施さない場合、ＳＡＮ型バンドパスフィルタ３０は、図７Ｂに示した周波数特性のままである。この場合、図９Ａに示すように、第１適応ノッチフィルタ２４の周波数特性は、図７Ｃの特性と比べて、６７Ｈｚ近傍のゲインが相対的に小さくなる。その結果、図９Ｂに示すようなＡＮＣ装置１０の感度関数が得られ、図９Ｃに示すような誤差信号Ａのスペクトラム図が得られる。
【００７１】
図９Ｃにおいて、実線のグラフは、第２スペクトラムＳＰＣ２を備える誤差信号Ａに対し、ＡＮＣ制御を実行した後の周波数特性である。また、破線のグラフは、第１スペクトラムＳＰＣ１を備える誤差信号Ａに対し、ＡＮＣ制御を実行した後の周波数特性である。このように、振動騒音ＮＳの共振周波数がシフトして対象周波数Ｆｃから僅かに外れた場合、その共振周波数近傍における振動騒音ＮＳを十分に相殺できない。
【００７２】
これに対して、本実施の形態に係るＡＮＣ装置１０では、能動型振動騒音制御部１４は、共振周波数のシフトに追従して、ＳＡＮ型バンドパスフィルタ３０の通過帯域を動的に変更する。具体的には、周波数切替部５４は、周波数変化量ｄＦ（＝−３Ｈｚ）を算出した後、対象周波数Ｆｃを７０Ｈｚから６７Ｈｚに切り替える。これにより、ＳＡＮ型バンドパスフィルタ３０の周波数特性は、図８Ｂの破線で示す特性から、同図の実線で示す特性に変更される。
【００７３】
すなわち、上述した周波数切替処理を施した場合、図１０Ａに示すように、第１適応ノッチフィルタ２４の周波数特性は、図９Ａの特性と比べて、６７Ｈｚ近傍のゲインが相対的に大きくなる。その結果、図１０Ｂに示すようなＡＮＣ装置１０の感度関数が得られ、図１０Ｃに示すような誤差信号Ａのスペクトラム図が得られる。
【００７４】
図１０Ｃにおいて、実線のグラフは、第２スペクトラムＳＰＣ２を備える誤差信号Ａに対し、ＡＮＣ制御を実行した後の周波数特性である。また、破線のグラフは、第１スペクトラムＳＰＣ１を備える誤差信号Ａに対し、ＡＮＣ制御を実行した後の周波数特性である。このように、振動騒音ＮＳの共振周波数がシフトした場合であっても、そのシフトの前後にわたって略同程度の相殺効果が得られた。
【００７５】
以上のように、第２適応ノッチフィルタ４４のフィルタ係数Ｗ２（Ｒｗ２、Ｉｗ２）の複素平面上での位相角度θと、前回の更新の際に算出した前回位相角度θｏｌｄとの間の位相角度変化量ｄθを算出し、位相角度変化量ｄθに応じて基準信号Ｘの対象周波数Ｆｃを切り替える周波数切替部５４を設けたので、フィルタ係数Ｗ２の複素平面上での位相角度θの変化量を遂次監視可能であり、位相角度変化量ｄθから周波数特性の変化の動向を簡便に且つ精度良く把握できる。これにより、振動騒音ＮＳの周波数特性に変化があっても、その変化に追従した振動騒音ＮＳの相殺制御を実行できる。
【００７６】
また、能動型振動騒音制御部１４は、２つの適応ノッチフィルタ、すなわち、ＡＮＣ制御に供される第１適応ノッチフィルタ２４と、対象周波数Ｆｃの更新制御に供される第２適応ノッチフィルタ４４とを備えている。位相角度θを算出するための第２適応ノッチフィルタ４４を別途設けたので、この算出処理の際、ＡＮＣ制御の状態の影響を受けにくくなる。これにより、位相角度θの計算精度がさらに高くなり、周波数切替制御の効果及び信頼性が高くなる。
【００７７】
なお、この発明は、上述した実施形態に限定されるものではなく、この発明の主旨を逸脱しない範囲で自由に変更できることは勿論である。
【符号の説明】
【００７８】
１０…ＡＮＣ装置 １１…車両
１４…能動型振動騒音制御部 １６…マイクロフォン
１８…スピーカ ２０…周波数設定部
２２…基準信号生成部 ２４…第１適応ノッチフィルタ
２６、４２、４６…減算器 ２８…第１フィルタ係数更新部
３０、５０…ＳＡＮ型バンドパスフィルタ ３６…振幅位相調整部
３８…振幅位相切替部 ４０…補正相殺信号生成部
４４…第２適応ノッチフィルタ ４８…第２フィルタ係数更新部
５４…周波数切替部
【特許請求の範囲】
【請求項１】
振動騒音に対する相殺信号に基づく相殺振動騒音を出力する振動騒音相殺手段と、
前記振動騒音と前記相殺振動騒音との干渉による残留振動騒音を誤差信号として検出する誤差信号検出手段と、
前記誤差信号が入力され、前記相殺信号を生成する能動型振動騒音制御手段と
を有する能動型振動騒音制御装置であって、
前記能動型振動騒音制御手段は、
所定周波数の基準信号を生成する基準信号生成手段と、
前記基準信号が入力され、前記相殺信号の生成に供される制御信号を出力する適応ノッチフィルタと、
前記基準信号の周波数に応じた振幅又は位相の調整値を格納し、前記制御信号の振幅又は位相を調整することで前記相殺信号を生成する振幅位相調整手段と、
前記誤差信号から前記制御信号を減算して補正誤差信号を生成する補正誤差信号生成手段と、
前記基準信号と前記補正誤差信号とに基づいて、前記補正誤差信号が最小となるように前記適応ノッチフィルタのフィルタ係数を逐次更新するフィルタ係数更新手段と、
前記振動騒音相殺手段から前記誤差信号検出手段までの伝達特性に基づいて前記相殺信号を補正して補正相殺信号を生成する補正相殺信号生成手段と、
前記誤差信号から前記補正相殺信号を減算することで、非制御時検出信号を生成する非制御時検出信号生成手段と、
複素平面上で定義されたフィルタ係数を備えており、前記基準信号が入力され、第２制御信号を出力する第２適応ノッチフィルタと、
前記非制御時検出信号から前記第２制御信号を減算して第２補正誤差信号を生成する第２補正誤差信号生成手段と、
前記基準信号と前記第２補正誤差信号とに基づいて、前記第２補正誤差信号が最小となるように前記第２適応ノッチフィルタのフィルタ係数を逐次更新する第２フィルタ係数更新手段と、
前記第２適応ノッチフィルタのフィルタ係数の前記複素平面上での位相角度と、前回の更新の際に算出した位相角度との間の位相角度変化量を算出し、前記位相角度変化量に応じて前記基準信号の周波数を切り替える周波数切替手段と
を備えることを特徴とする能動型振動騒音制御装置。
【請求項２】
請求項１記載の能動型振動騒音制御装置において、
前記周波数切替手段は、前記誤差信号のサンプリング周期と、前記位相角度変化量とに基づいて周波数変化量を算出し、前記周波数変化量が下限闘値を下回った場合に前記基準信号の周波数を維持することを特徴とする能動型振動騒音制御装置。
【請求項３】
請求項２記載の能動型振動騒音制御装置において、
前記周波数切替手段は、前記周波数変化量が、前記下限閾値よりも大きな上限闘値を上回った場合に前記基準信号の周波数を維持することを特徴とする能動型振動騒音制御装置。
【請求項４】
請求項１〜３のいずれか１項に記載の能動型振動騒音制御装置において、
前記周波数切替手段が前記基準信号の周波数を切り替えたことに応じて、前記振幅位相調整手段が格納する前記調整値を切り替える振幅位相切替手段をさらに備えることを特徴とする能動型振動騒音制御装置。
【請求項１】
振動騒音に対する相殺信号に基づく相殺振動騒音を出力する振動騒音相殺手段と、
前記振動騒音と前記相殺振動騒音との干渉による残留振動騒音を誤差信号として検出する誤差信号検出手段と、
前記誤差信号が入力され、前記相殺信号を生成する能動型振動騒音制御手段と
を有する能動型振動騒音制御装置であって、
前記能動型振動騒音制御手段は、
所定周波数の基準信号を生成する基準信号生成手段と、
前記基準信号が入力され、前記相殺信号の生成に供される制御信号を出力する適応ノッチフィルタと、
前記基準信号の周波数に応じた振幅又は位相の調整値を格納し、前記制御信号の振幅又は位相を調整することで前記相殺信号を生成する振幅位相調整手段と、
前記誤差信号から前記制御信号を減算して補正誤差信号を生成する補正誤差信号生成手段と、
前記基準信号と前記補正誤差信号とに基づいて、前記補正誤差信号が最小となるように前記適応ノッチフィルタのフィルタ係数を逐次更新するフィルタ係数更新手段と、
前記振動騒音相殺手段から前記誤差信号検出手段までの伝達特性に基づいて前記相殺信号を補正して補正相殺信号を生成する補正相殺信号生成手段と、
前記誤差信号から前記補正相殺信号を減算することで、非制御時検出信号を生成する非制御時検出信号生成手段と、
複素平面上で定義されたフィルタ係数を備えており、前記基準信号が入力され、第２制御信号を出力する第２適応ノッチフィルタと、
前記非制御時検出信号から前記第２制御信号を減算して第２補正誤差信号を生成する第２補正誤差信号生成手段と、
前記基準信号と前記第２補正誤差信号とに基づいて、前記第２補正誤差信号が最小となるように前記第２適応ノッチフィルタのフィルタ係数を逐次更新する第２フィルタ係数更新手段と、
前記第２適応ノッチフィルタのフィルタ係数の前記複素平面上での位相角度と、前回の更新の際に算出した位相角度との間の位相角度変化量を算出し、前記位相角度変化量に応じて前記基準信号の周波数を切り替える周波数切替手段と
を備えることを特徴とする能動型振動騒音制御装置。
【請求項２】
請求項１記載の能動型振動騒音制御装置において、
前記周波数切替手段は、前記誤差信号のサンプリング周期と、前記位相角度変化量とに基づいて周波数変化量を算出し、前記周波数変化量が下限闘値を下回った場合に前記基準信号の周波数を維持することを特徴とする能動型振動騒音制御装置。
【請求項３】
請求項２記載の能動型振動騒音制御装置において、
前記周波数切替手段は、前記周波数変化量が、前記下限閾値よりも大きな上限闘値を上回った場合に前記基準信号の周波数を維持することを特徴とする能動型振動騒音制御装置。
【請求項４】
請求項１〜３のいずれか１項に記載の能動型振動騒音制御装置において、
前記周波数切替手段が前記基準信号の周波数を切り替えたことに応じて、前記振幅位相調整手段が格納する前記調整値を切り替える振幅位相切替手段をさらに備えることを特徴とする能動型振動騒音制御装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【公開番号】特開２０１３−１１２１４０（Ｐ２０１３−１１２１４０Ａ）
【公開日】平成２５年６月１０日（２０１３．６．１０）
【国際特許分類】
【出願番号】特願２０１１−２５９６８８（Ｐ２０１１−２５９６８８）
【出願日】平成２３年１１月２９日（２０１１．１１．２９）
【出願人】（０００００５３２６）本田技研工業株式会社 (23,863)
【Ｆターム（参考）】
【公開日】平成２５年６月１０日（２０１３．６．１０）
【国際特許分類】
【出願日】平成２３年１１月２９日（２０１１．１１．２９）
【出願人】（０００００５３２６）本田技研工業株式会社 (23,863)
【Ｆターム（参考）】
[ Back to top ]