(57)【要約】【課題】騒音抑制や速応性などに関する制御性能およびロバスト安定性の両面を考慮して複数の周波数帯域の振動の低減をそれぞれ対応して制御する複数のコントローラの設計方法等を提供することにある。【解決手段】複数のコントローラのうち１つのコントローラについて、１つのコントローラを除く他の１以上のコントローラでフィードバックした状態で１つのコントローラの制御対象からモデル誤差を求め、当該モデル誤差を用いてロバスト安定性を満たす重み関数を求めて１つのコントローラを設計する。これにより、周波数帯域を分割して各帯域のモデルを低次元で同定することができ、モデル化のための計算時間の短縮、モデル化精度の向上を達成することができる。各帯域用コントローラは、それぞれの帯域に対応した低次元モデルから設計されるため、制御系設計の高効率化を図ることができると同時に設計されるコントローラも低次元で構成することができる。【選択図】図１１
2008257348.tif 000002
【特許請求の範囲】
【請求項１】
振動源により発生する振動を低減する装置に備えられ、複数の周波数帯域の振動の低減をそれぞれ対応して制御する、複数のコントローラの設計方法であって、
前記複数のコントローラのうち１つのコントローラについて、前記１つのコントローラを除く他の１以上のコントローラでフィードバックした状態で前記１つのコントローラの制御対象からモデル誤差を求め、当該モデル誤差を用いてロバスト安定性を満たす重み関数を求めて前記１つのコントローラを設計する
ことを特徴とする振動低減装置のコントローラの設計方法。
【請求項２】
振動源により発生する振動を低減する装置に備えられ、複数の周波数帯域の振動の低減をそれぞれ対応して制御する、複数のコントローラの設計方法であって、
前記複数のコントローラのうち１つのコントローラについて、前記１つのコントローラを除く他の全てのコントローラでフィードバックした状態で前記１つのコントローラの制御対象からモデル誤差を求め、当該モデル誤差を用いてロバスト安定性を満たす重み関数を求めて前記１つのコントローラを設計する
ことを特徴とする振動低減装置のコントローラの設計方法。
【請求項３】
振動源により発生する振動を低減する装置に備えられ、第１の周波数帯域の振動の低減を制御する第１の周波数帯域コントローラと第２の周波数帯域の振動の低減を制御する第２の周波数帯域コントローラの設計方法であって、
前記第２の周波数帯域コントローラが無いとした状態で、前記第１の周波数帯域コントローラの制御対象から第１の周波数帯域モデル誤差を求め、当該第１の周波数帯域モデル誤差を用いてロバスト安定性を満たす第１の周波数重み関数を求めて前記第１の周波数帯域コントローラを設計する、第１のステップと、
当該第１の周波数帯域コントローラでフィードバックした状態の周波数応答を実制御対象の応答とみなして、前記第２の周波数帯域コントローラの制御対象から第２の周波数帯域モデル誤差を求め、当該第２の周波数帯域モデル誤差を用いてロバスト安定性を満たす第２の周波数重み関数を求めて前記第２の周波数帯域コントローラを設計する、第２のステップと、
当該第２の周波数帯域コントローラでフィードバックした状態の周波数応答を実制御対象の応答とみなして、前記第１の周波数帯域コントローラの制御対象から新たな第１の周波数帯域モデル誤差を求め、前記第１の周波数重み関数が当該第１の周波数帯域モデル誤差に対しロバスト安定性を満たしていれば前記各コントローラの設計を終了し、
一方、前記第１の周波数重み関数が当該第１の周波数帯域モデル誤差に対しロバスト安定性を満たしていないときは、当該第１の周波数帯域モデル誤差を用いてロバスト安定性を満たすように前記第１の周波数重み関数を修正して前記第１の周波数帯域コントローラを再設計し次の第４のステップに行く、第３のステップと、
更に、当該第１の周波数帯域コントローラでフィードバックした状態の周波数応答を実制御対象の応答とみなして、前記第２の周波数帯域コントローラの制御対象から新たな第２の周波数帯域モデル誤差を求め、前記第２の周波数重み関数が当該第２の周波数帯域モデル誤差に対しロバスト安定性を満たしていれば前記各コントローラの設計を終了し、
一方、前記第２の周波数重み関数が当該第２の周波数帯域モデル誤差に対しロバスト安定性を満たしていないときは、当該第２の周波数帯域モデル誤差を用いてロバスト安定性を満たすように前記第２の周波数重み関数を修正して前記第２の周波数帯域コントローラを再設計し上記第３のステップに戻る、第４のステップとを含み、
前記各コントローラが相互にロバスト安定性を満たすまで、前記各コントローラを設計する
ことを特徴とする振動低減装置のコントローラの設計方法。
【請求項４】
請求項１〜３の何れか一項に記載の振動低減装置のコントローラの設計方法により設計されたことを特徴とする振動低減装置のコントローラ。
【請求項５】
請求項４に記載のコントローラを備えたことを特徴とする振動低減装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、振動源により発生する振動を低減する装置に備えられ、複数の周波数帯域の振動の低減をそれぞれ対応して制御する、複数のコントローラの設計方法、振動低減装置のコントローラ及び振動低減装置に関する。
【背景技術】
【０００２】
優れた制御性能およびロバスト安定性能を達成させる制御法として、制御性能をＨ２ノルム、ロバスト安定性をＨ∞ノルムで評価する混合Ｈ２／Ｈ∞制御法がある。
【０００３】
図３５は、制御問題を表すブロック線図である。なお、図３５において、Ｐ（s）は制御対象、ｙは観測出力、ｕは制御入力である。外乱ｗから周波数重み付き制御応答Ｚ２までの伝達関数行列をＴＺ２ｗ（s）、外乱ｗから周波数重み付き制御入力Ｚ１までの伝達関数行列をＴＺ１ｗ（s）とすると、混合Ｈ２／Ｈ∞制御問題は次式（１）、（２）のようになる。
【０００４】
【数１】
2008257348.tif 000003
ここで、式（１）はＴＺ２ｗ（s）のＨ２ノルム、式（２）はＴＺ１ｗ（s）のＨ∞ノルムを表す。モデル誤差Δ∞（s）に対して適切に周波数重み関数Ｗ１（s）を与え、制振目標に対応した周波数重み関数Ｗ２（s）を設定した上で、制御問題である式（１）、（２）を解くことにより、モデル誤差に対しロバスト安定性を有するコントローラＫ（s）が設計される。本コントローラは、ロバスト安定性を満たしつつ、制御応答（振動応答）のエネルギーを最小化できるため、振動制御における有効性が高い。
【０００５】
騒音制御を行う方法として、フィードフォワード制御を基本とした能動的騒音制御手法がある。この方法は、騒音に関する参照信号とスピーカーで騒音を打ち消すべき計測点における誤差信号から騒音特性を学習により獲得し、スピーカーから騒音と逆位相の音を出力することにより騒音を相殺する制御法である。このような制御手法により、低周波から高周波にわたる広い周波数帯域で音響制御を行う場合、技術的な困難さが存在する。それは、広い周波数帯域でのモデル化は精度が低下することに加え次数が大きくなるため、制御回路における演算負荷の面から高周波数領域における制御は困難になる。これを解決する方法として、上記の能動的騒音制御手法において、低周波数帯域と高周波数帯域に分けて騒音特性の学習および制御回路を構成する方法がある。フィルターにより両制御回路の干渉を防ぐことにより、低周波数帯域と高周波数帯域の騒音制御を可能にしている（特許文献１参照）。
【０００６】
【特許文献１】特許第３４７１３７０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
上述した低周波数帯域と高周波数帯域に分けて騒音特性の学習および制御回路を構成する方法を用いる場合、フィードフォワード制御を基本としているため、騒音特性を学習により獲得した際のモデル誤差に対するロバスト安定性を考慮した設計はできない。さらに、騒音制御以外の制御、例えば一般的な振動制御に関してはロバスト安定性を考慮したフィードバック制御が主流であり、上述のフィードフォワード制御に基づく能動的騒音制御手法では騒音抑制や速応性などに関する制御性能およびロバスト安定性の両面で仕様を満足できないことが起こり得る。
【０００８】
本発明は、上記のような課題に鑑みなされたものであり、その目的は、騒音抑制や速応性などに関する制御性能およびロバスト安定性の両面を考慮して複数の周波数帯域の振動の低減をそれぞれ対応して制御する、複数のコントローラの設計方法、振動低減装置のコントローラ及び振動低減装置を提供することにある。
【課題を解決するための手段】
【０００９】
上記目的達成のため、本発明の振動低減装置のコントローラの設計方法では、振動源により発生する振動を低減する装置に備えられ、複数の周波数帯域の振動の低減をそれぞれ対応して制御する、複数のコントローラの設計方法であって、前記複数のコントローラのうち１つのコントローラについて、前記１つのコントローラを除く他の１以上のコントローラでフィードバックした状態で前記１つのコントローラの制御対象からモデル誤差を求め、当該モデル誤差を用いてロバスト安定性を満たす重み関数を求めて前記１つのコントローラを設計することを特徴としている。これにより、周波数帯域を分割して各帯域のモデルを低次元で同定することができ、モデル化のための計算時間の短縮およびモデル化精度の向上を達成することができる。そして、各帯域用のコントローラは、それぞれの帯域に対応した低次元モデルから設計されるため、制御系設計の高効率化を図ることができると同時に設計されるコントローラも低次元で構成することができる。
【００１０】
また、振動源により発生する振動を低減する装置に備えられ、複数の周波数帯域の振動の低減をそれぞれ対応して制御する、複数のコントローラの設計方法であって、前記複数のコントローラのうち１つのコントローラについて、前記１つのコントローラを除く他の全てのコントローラでフィードバックした状態で前記１つのコントローラの制御対象からモデル誤差を求め、当該モデル誤差を用いてロバスト安定性を満たす重み関数を求めて前記１つのコントローラを設計することを特徴としている。これにより、２つ以上のコントローラでも、他の全てのコントローラの影響を誤差として考慮した協調設計をすることで、大きな制御効果を得ることができ、さらに、確実に制御系を安定させることができる。
【００１１】
より具体的には、振動源により発生する振動を低減する装置に備えられ、第１の周波数帯域の振動の低減を制御する第１の周波数帯域コントローラと第２の周波数帯域の振動の低減を制御する第２の周波数帯域コントローラの設計方法であって、前記第２の周波数帯域コントローラが無いとした状態で、前記第１の周波数帯域コントローラの制御対象から第１の周波数帯域モデル誤差を求め、当該第１の周波数帯域モデル誤差を用いてロバスト安定性を満たす第１の周波数重み関数を求めて前記第１の周波数帯域コントローラを設計する、第１のステップと、当該第１の周波数帯域コントローラでフィードバックした状態の周波数応答を実制御対象の応答とみなして、前記第２の周波数帯域コントローラの制御対象から第２の周波数帯域モデル誤差を求め、当該第２の周波数帯域モデル誤差を用いてロバスト安定性を満たす第２の周波数重み関数を求めて前記第２の周波数帯域コントローラを設計する、第２のステップと、当該第２の周波数帯域コントローラでフィードバックした状態の周波数応答を実制御対象の応答とみなして、前記第１の周波数帯域コントローラの制御対象から新たな第１の周波数帯域モデル誤差を求め、前記第１の周波数重み関数が当該第１の周波数帯域モデル誤差に対しロバスト安定性を満たしていれば前記各コントローラの設計を終了し、一方、前記第１の周波数重み関数が当該第１の周波数帯域モデル誤差に対しロバスト安定性を満たしていないときは、当該第１の周波数帯域モデル誤差を用いてロバスト安定性を満たすように前記第１の周波数重み関数を修正して前記第１の周波数帯域コントローラを再設計し次の第４のステップに行く、第３のステップと、更に、当該第１の周波数帯域コントローラでフィードバックした状態の周波数応答を実制御対象の応答とみなして、前記第２の周波数帯域コントローラの制御対象から新たな第２の周波数帯域モデル誤差を求め、前記第２の周波数重み関数が当該第２の周波数帯域モデル誤差に対しロバスト安定性を満たしていれば前記各コントローラの設計を終了し、一方、前記第２の周波数重み関数が当該第２の周波数帯域モデル誤差に対しロバスト安定性を満たしていないときは、当該第２の周波数帯域モデル誤差を用いてロバスト安定性を満たすように前記第２の周波数重み関数を修正して前記第２の周波数帯域コントローラを再設計し上記第３のステップに戻る、第４のステップとを含み、前記各コントローラが相互にロバスト安定性を満たすまで、前記各コントローラを設計することを特徴としている。
【００１２】
上記目的達成のため、本発明の振動低減装置のコントローラでは、上記各振動低減装置のコントローラの設計方法により設計されたことを特徴としている。また、上記目的達成のため、本発明の振動低減装置では、上記コントローラを備えたことを特徴としている。これにより、上記各作用効果を奏する振動低減装置のコントローラ及び振動低減装置を提供することができる。
【発明を実施するための最良の形態】
【００１３】
本発明の実施形態について、図面を参照して説明する。尚、以下に説明する実施形態は特許請求の範囲に係る発明を限定するものではなく、また実施形態の中で説明されている特徴の組み合わせの全てが発明の解決手段に必須であるとは限らない。
【００１４】
本発明の実施形態に係る設計方法は、振動低減装置が備える複数のコントローラに対して有効に適用されるものであるが、説明を分かり易くするために、先ず単一のコントローラに対して適用する場合について説明する。
【００１５】
図１は、単一のコントローラを備えた振動低減装置を示す構成図である。この振動低減装置１０は、システム（制御対象）１１の振動を計測する加速度センサ１２、計測信号のノイズなどを除去するローパスフィルター１３、計測信号から制御則にしたがって制御入力を演算するためのデジタル制御システム１４、制御対象１１に制御入力を作用させるためのアクチュエータ１５により構成される。一方、データ処理により振動特性を評価するための測定系は、入出力信号とＦＦＴアナライザ１６から成る。この測定系は、制御対象１１の制御には何ら関与せず、データ処理・振動特性評価以外には用いる必要はない。制御対象１１には外力などの外乱入力が作用し、この外乱入力により制御対象１１が励振されることになるが、振動低減装置１０における加速度センサ１２で検出された信号がセンサーアンプ１７を介してアナログのローパスフィルター１３を通り、ＡＤ変換器１４ａを経て制御則を演算するコントローラ１４ｂに入力される。そして、制御入力信号がＤＡ変換器１４ｃで出力され、増幅器１８を経て制御用のアクチュエータ１５に入力される。アクチュエータ１５により制御対象１１に制御力が作用し、制御対象１１の振動が抑制される。
【００１６】
図２および図３は、振動低減装置１０の制御系の概略のブロック線図である。図２に示すように、制御対象１１には、外力に関する外乱ｗｄと制御入力に加算される外乱ｗｕが作用する。図３５および図２から、ｗｕからｕまでの伝達関数がロバスト安定性の評価指標となる。一方、外力外乱ｗｄに関しては、それを計測することや正確な特性を把握することが困難な場合がある。この場合は、図３に示すように、外乱をｗｕのみとし、ｗｄが作用したときと等価な外乱をｗｕとして作用させるという考えに基づいて対応することができる。このようにすれば、正確な外力外乱の作用点および特性がわからなくても、外力による大まかな出力特性を把握しておけば、それをｗｕの特性に反映させて、制御系を設計することができる。また、図３は図２の特殊なケースと見ることができるので、図２の制御系についてコントローラの設計法を構築しておけば、図３についても同様に対応することができる。
【００１７】
次に、周波数応答の測定方法について説明する。励振入力信号（外力外乱）と加速度センサ１２で検出された加速度信号をＦＦＴアナライザ１６に入力することにより、ＦＦＴアナライザ１６内で高速フーリエ変換に基づくデータ処理が行われ、入出力間の周波数応答がアクセレランス（｜加速度／力｜）やコンプライアンス（｜変位／力｜）として演算され結果が出力される。同様に、制御入力外乱として励振入力信号を与え、制御用のアクチュエータ１５から制御対象１１に励振力を作用させ、制御対象１１の振動応答を加速度センサ１２で計測し、制御入力外乱信号と加速度センサ１２の出力信号をＦＦＴアナライザ１６でデータ処理することにより周波数応答を計測する。一方、有限要素法などによる解析モデルが存在する場合、コンピュータを用いたシミュレーションにより得られた周波数応答を用いることができる。
【００１８】
次に、周波数応答のカーブフィッティングの方法について説明する。計測された周波数応答から制御対象１１のモデルを記述するためには、制御対象１１のモード特性（固有振動数、固有モードおよびモード減衰比など）を同定する必要がある。周波数応答（アクセレランスやコンプライアンスなど）からカーブフィッティングによりモード特性を同定する方法としては、偏分反復法などに代表される実験モード解析の手法が主流である。一方、有限要素法などによる解析モデルが存在する場合、コンピュータを用いたシミュレーションにより上記モード特性が計算されるので、同定手法によらず、これらの解析結果を用いることも可能である。
【００１９】
次に、モデル誤差と周波数重み関数の作成方法について説明する。上述のシステム同定手法によりモデル化を行った場合、モデルによる周波数特性と実際に計測で得られた（実システムの）周波数特性との間には誤差が存在するのが一般的である。制御系を構成する際のコントローラ１４ｂは、この誤差に対してロバスト安定性を有するように設計されなければならない。
【００２０】
図４は、モデル誤差を乗法的誤差として与え、外乱の作用する位置とあわせて表記したブロック線図である。図４において、Ｐ（s）は制御対象、Ｋ（s）はコントローラ、Δ∞（s）は制御対象の乗法的誤差、Ｐ~（s）は実際の制御対象（実システム）である。一方、図３５のブロック線図において、モデル誤差Δ∞（s）に対して適切に周波数重み関数Ｗ１（s）を与えることにより、ロバスト安定性を保証する制御系の設計が可能になる。すなわち、Δ∞（s）の周波数応答（σ−プロット）に対して次式（３）を満たすＷ１（s）によりスモールゲイン定理に基づくロバスト安定条件である次式（４）を満足するコントローラを設計する。この混合Ｈ２／Ｈ∞制御問題に基づくコントローラは、線形行列不等式（ＬＭＩ）の形で記述された制御問題に対し数値探索法により設計される。
【００２１】
【数２】
2008257348.tif 000004
以上のようにして設計されたコントローラは、Δ∞（s）という誤差の範囲内で、ロバスト安定性を満たすことができる。一方、制御応答の性能（振動の低減）に関しては、Δ∞（s）により規定されるロバスト安定条件の制約に依存する。すなわち、誤差Δ∞（s）が対象とする周波数領域において小さければ制御応答に関して高い制御性能を達成することは可能であるが、逆に誤差が大きい場合にはその性能は劣化する。広い周波数帯域を制御する場合、単一のコントローラで制御しようとすると、多くの振動モードを含む広帯域の特性を一つのモデルで表現する必要があり、必然的にシステム同定の精度低下（モデル誤差の増加）とモデル次数の増加をもたらす。このことは、結果として制御性能の低下につながる。
【００２２】
そこで、上述の単一コントローラによる広帯域制御の問題点を解決するために、各帯域用の複数コントローラをロバスト協調設計することにより構成される振動低減装置を説明する。
【００２３】
図５は、本発明の一実施の形態に係る複数のコントローラを備えた振動低減装置を示す構成図である。この振動低減装置２０は、低周波数帯域用のコントローラと高周波数帯域用のコントローラが並列結合されて広帯域用の制御器を構成しており、低周波数帯域コントローラ２１および高周波数帯域コントローラ２２以外の基本構成は前述の図１と同一であるため同一番号を付して説明を省略する。ここで、各周波数帯域用のコントローラ２１、２２は、それぞれ他の周波数帯域用のコントローラ２２、２１との干渉による不安定現象を起こさぬよう、ロバスト安定仕様を満足するように設計される。
【００２４】
次に、低周波数帯域コントローラ２１および高周波数帯域コントローラ２２のロバスト協調設計方法について説明する。
【００２５】
図６は、低周波数帯域用と高周波数帯域用の二つのコントローラを並列に結合した制御系のブロック線図である。ここで、ＫＬ（s）は低周波数帯域用コントローラ、ＫＨ（s）は高周波数帯域用コントローラである。フィードバック制御系において、両コントローラは互いに干渉するため、ロバスト安定性を考慮した両コントローラの協調設計を行う必要がある。そこで、一方のコントローラによるフィードバック制御による特性変化を他方のコントローラ設計におけるモデル誤差と考え、ロバスト制御系の設計を行う。ここで、各コントローラによるフィードバック制御系の周波数特性は、設計したコントローラと計測で得られた周波数応答から予測することができる。
【００２６】
以下に、複数コントローラの並列結合制御系において、ロバスト安定性を保証すべき協調設計法のプロセスを示す。
【００２７】
図７は、低周波数帯域コントローラのモデル誤差を乗法的誤差として与え、外乱の作用する位置とあわせて表記したブロック線図、図８は、低周波数帯域コントローラでフィードバックしたブロック線図、図９は、高周波数帯域コントローラのモデル誤差を乗法的誤差として与え、外乱の作用する位置とあわせて表記したブロック線図、図１０は、高周波数帯域コントローラでフィードバックしたブロック線図である。
【００２８】
<ステップ１>：図７において高周波数帯域コントローラが無い場合（ＫＨ（s）＝０）に、低周波数帯域制御対象ＰＬ（s）とモデル誤差Δ∞（s）を用い、次式（５）となるようにＷ∞Ｌ（s）を与え、ロバスト安定性を満たすように低周波数帯域コントローラＫＬ（s）を設計する。
【００２９】
【数３】
2008257348.tif 000005
<ステップ２>：低周波数帯域コントローラＫＬ（s）でフィードバックしたシステムＰ~（s）をＫＬ（s）と（計測した）実システムＰ~（s）から次式（６）により表す（図８参照）。
【００３０】
【数４】
2008257348.tif 000006
ここで、Ｐ~（ｊω）が計測された周波数応答に対応するので、式（２）の周波数応答を実システムの応答とみなした上で、高周波数帯域制御対象ＰＨ（s）とのモデル誤差Δ∞Ｈ（s）を次式（７）から求める。
【００３１】
【数５】
2008257348.tif 000007
そして、高周波数帯域制御対象ＰＨ（s）とモデル誤差Δ∞Ｈ（s）を用い、次式（８）となるようにＷ∞Ｈ（s）を与え、ロバスト安定性を満たすように高周波数帯域コントローラＫＨ（s）を設計する（図９参照）。
【００３２】
【数６】
2008257348.tif 000008
<ステップ３>：高周波数帯域コントローラＫＨ（s）でフィードバックしたシステムＰ~Ｌ（s）をＫＨ（s）と（計測した）実システムＰ~（s）から次式（９）により表す（図１０参照）。
【００３３】
【数７】
2008257348.tif 000009
この周波数応答を実システムの応答とみなした上で低周波数帯域制御対象ＰＬ（s）とのモデル誤差Δ∞Ｌ（s）を次式（１０）から求める。
【００３４】
【数８】
2008257348.tif 000010
もし、低周波数帯域周波数重み関数Ｗ∞Ｌ（s）がΔ∞Ｌ（s）に対しロバスト安定性（次式（１１））を満たしていれば終了する。そうでなければ、低周波数帯域制御対象ＰＬ（s）とモデル誤差Δ∞Ｌ（s）を用い、次式（１２）を満たすようにＷ∞Ｌ（s）を修正した上で低周波数帯域コントローラＫＬ（s）を再設計し、ステップ４に行く（図７参照）。
【００３５】
【数９】
2008257348.tif 000011
<ステップ４>：ステップ２の手順と同様に、低周波数帯域コントローラＫＬ（s）でフィードバックしたシステムＰ~Ｈ（s）の周波数応答をＫＬ（s）と計測した（実システムＰ~（s））周波数応答から計算し、それを実システムの応答とみなした上で高周波数帯域制御対象ＰＨ（s）とのモデル誤差Δ∞Ｈ（s）を求める。もし、高周波数帯域周波数重み関数Ｗ∞Ｈ（s）がΔ∞Ｈ（s）に対しロバスト安定性（次式（１３））を満たしていれば終了する。そうでなければ、高周波数帯域制御対象ＰＨ（s）とモデル誤差Δ∞Ｈ（s）を用い、次式（１４）を満たすようにＷ∞Ｈ（s）を修正した上で高周波数帯域コントローラＫＨ（s）を再設計し、ステップ３に戻る（図９参照）。
【００３６】
【数１０】
2008257348.tif 000012
ここで、本設計プロセスにおける設計終了条件は次式（１５）のとおりである。すなわち、低周波数帯域用コントローラと高周波数帯域用コントローラがそれぞれ互いにロバスト安定性を満たすときに設計終了となる。ここで、式（１５）の最左項は行列Δ∞Ｌ（ｊω）の最大特異値を示し、次式（１６）（尚、＊は複素転置行列）の最大固有値を表す。Ｗ∞Ｌ（ｊω）、Δ∞Ｈ（ｊω）およびＷ∞Ｌ（ｊω）の最大特異値についても同様である。
【００３７】
【数１１】
2008257348.tif 000013
図１１は、ステップ１からステップ４をフロー図にしたものである。このフロー図におけるプロセスは、上記終了条件を満たすまでコントローラの設計を繰り返すことを示している。上記プロセスにより設計されたそれぞれのコントローラは、お互いの干渉に対しロバスト安定性を満たすことになる。なお、低周波数帯域コントローラが無い場合（ＫＬ（s）＝０）に、高周波数帯域制御対象ＰＨ（s）とモデル誤差Δ∞（s）を用い、Ｗ∞Ｈ（s）を与え、ロバスト安定性を満たすように高周波数帯域コントローラＫＨ（s）を設計し、以降も設計する方法であっても良い。また、ステップ２〜ステップ４において、ＫＬ（s）（もしくはＫＨ（s））でフィードバック制御したときの実システムＰ~Ｈ（s）（もしくはＰ~Ｌ（s））の周波数応答として、制御実験による実測データを用いても構わない。
【００３８】
さらに、コントローラの数が３以上の場合でも、同様のステップでコントローラを設計できる。すなわち、初期状態のコントローラを全て０とし、各帯域のコントローラを他の帯域の全てのコントローラでフィードバックした状態においてロバスト安定性を満たすように上記ステップにしたがって設計する。例えば、コントローラの数が３の場合には、次のようになる。コントローラ１は、他のすべてのコントローラを０として設計する。コントローラ２は、コントローラ１をフィードバックし、コントローラ３を０として設計する。コントローラ３は、コントローラ１および２をフィードバックした状態で設計する。以降、再設計が不要になるまで、本設計ステップにしたがってコントローラ１、２、３の順に設計する。
【００３９】
以上のように、各帯域用の複数コントローラをロバスト協調設計する方法によれば、周波数帯域を分割することにより、各帯域のモデルを低次元で同定することができ、モデル化のための計算時間の短縮およびモデル化精度の向上が達成される。そして、各帯域用の制御器は、それぞれの帯域に対応した低次元モデルから設計されるため、制御系設計の高効率化が図られると同時に設計される制御器も低次元で構成される。設計の高効率化は、設計期間の短縮とそれによるコスト低減等をもたらすと期待できる。さらに、制御器が低次元で構成できるため、制御器設計における数値計算の誤差が小さく、設計される制御器の信頼性は高いと言える。実システムで用いる広帯域用の制御器は、各帯域用の低次元制御器を並列に結合したものであるため、広帯域用の大次元モデルから直接設計した単一の制御器よりも計算効率を著しく向上できる。この根拠は、デジタル制御システムにおける１サンプルあたりの乗算回数は、コントローラの次数の２乗にほぼ等しく、これに対し制御帯域を分割した場合は各コントローラの乗算回数の和に等しくなることから理解できる。これにより、デジタル制御システムのハードウェア上のサンプリング周波数を向上させることができ、高周波数領域を含む広帯域の実時間制御が実現可能となる。
【００４０】
図１２は、単一のコントローラを備えた振動低減装置の実施例を図１に対応させて示す構成図である。この振動低減装置１０の基本構成は、前述の図１と同一であるが、図１との違いは、制御対象１１を具体的なエンジン・エンジェクタシステムとして与えていることと、外力に関する外乱をインパルスハンマーにより作用させていることである。
【００４１】
先ず、周波数応答の測定方法と結果について説明する。外乱入力をインパルスハンマーにより与え、加速度センサ１２により制御対象１１の加速度を計測し、インパルスハンマーによる加振力の入力信号と加速度センサ１２からの出力信号をＦＦＴアナライザ１６に入力することにより、力入力に対する加速度出力の周波数応答を計測する。
【００４２】
図１３は、計測された外乱加振入力から加速度センサ１２の出力信号までの周波数応答（アクセレランス）を示す図である。図１３から、４００Ｈｚ、１６００Ｈｚ、２３００Ｈｚに大きなピークが存在することが分かる。同様に、制御入力外乱として励振入力信号を与え、制御用のアクチュエータ１５から制御対象１１に励振力を作用させ、制御対象１１の振動応答を加速度センサ１２で計測し、制御入力外乱信号と加速度センサ１２の出力信号をＦＦＴアナライザ１６でデータ処理することにより周波数応答を計測する。
【００４３】
図１４は、計測された制御入力外乱信号から加速度センサ１２の出力信号までの周波数応答（アクセレランス）を示す図である。図１４の結果から、図１３における結果と同一の周波数に共振ピークが存在していることが分かる。
【００４４】
次に、周波数応答のカーブフィッティングの方法と結果について説明する。計測された周波数応答からカーブフィッティングよりモード特性を同定する。対象とする周波数応答としてアクセレランスやコンプライアンスなどがあるが、ここでは後者に基づいて同定を行う。実験ではアクセレランス（｜加速度／力｜）を計測しているので、これに−１／ω２（ωは角振動数）を乗じることによりコンプライアンス（｜変位／力｜）に変換する。モード特性を同定する方法として、偏分反復法を用いる。比例粘性減衰系におけるシステムの周波数応答（コンプライアンス）Ｇ（s）は、システムの運動方程式をラプラス変換したｓ領域において次式（１７）のように記述される。
【００４５】
【数１２】
2008257348.tif 000014
ここで、Ωｒおよびζｒはそれぞれｒ次の固有振動数および減衰比であり、φｒｉおよびφｒｊはそれぞれｒ次モードにおける入力点ｉおよび出力点ｊの固有モード成分である。これらをモード特性と呼び、実験モード解析におけるカーブフィッティングにより同定される。ｒ次モードに関する入出力関係を状態方程式で表すと次式（１８）、（１９）のようになる。
【００４６】
【数１３】
2008257348.tif 000015
ここで、ｑｅ，ｒ、ｙｒはそれぞれｒ次モードに関する状態量および出力であり、Ａｅ，ｒ、Ｂｌｅ，ｒ、Ｂ２ｅ，ｒおよびＣ２ｅ，ｒはｒ次モードに関する状態方程式のシステム行列を表す。ｗおよびｕは、それぞれ外乱入力および制御入力である。φｒａ、φｒｂおよびφｒｃは、それぞれｒ次モードにおける応答出力点ａ、外乱入力点ｂおよび制御入力点ｃの固有モード成分である。上式を全モードについて重ね合わせると、次式（２０）、（２１）のシステム方程式が得られる。
【００４７】
【数１４】
2008257348.tif 000016
ここで、ｑｅ，ｒ、ｙはそれぞれ全モード成分を含む状態量および出力であり、Ａｅ、Ｂｌｅ、Ｂ２ｅおよびＣ２ｅは全モード成分を含む状態方程式のシステム行列を表す。また、ＤｂおよびＤｃはそれぞれ外乱入力および制御入力に対する応答をカーブフィッティングした際に、高次モードを省略したことによる補正項として同定される剰余剛性である。一方、低次モードを省略することによる補正項は剰余質量として同定されるが、一般に剰余質量を用いなくても対象とする周波数帯域の低域側における冗長なモードを採用することにより、その影響を補正することが可能である。
【００４８】
図１５〜図１８は、振動実験により計測したデータを用いて、実験モード解析によりカーブフィットした結果（アクセレランス）を示す図である。図１５および図１６は、低周波数帯域において採用モード数を４としたときの、それぞれ外乱入力および制御入力に対する応答の同定結果である。そして、図１７および図１８は、高周波数帯域において採用モード数を６としたときの、それぞれ外乱入力および制御入力に対する応答の同定結果である。両周波数帯域において、実験結果を正確にフィットする良好なシステム同定が行われていることがわかる。さらに、対象とする周波数帯域を分割してモデル化しているため、全周波数帯域を一括してモデル化する場合よりも、低次元でモデルが実現できることになる。
【００４９】
次に、モデル誤差と周波数重み関数の作成方法と結果について説明する。
【００５０】
図１９は、ロバスト安定性を考慮した周波数重みの設定について示すブロック線図である。図１６および図１８に示した制御入力に対する応答の同定結果から、次式（２２）の乗法的誤差を計算する。ここで、Ｐ（s）はモデル化された制御対象、Ｐ~（s）は実際の制御対象（実システム）、Δ∞（s）はモデルの乗法的誤差である。
【００５１】
【数１５】
2008257348.tif 000017
図２０および図２１は、低周波数帯域および高周波数帯域において、両帯域の乗法的誤差Δ∞Ｌ（s）およびΔ∞Ｈ（s）を計算した結果をそれぞれ破線で示し、両帯域において与えた周波数重み関数Ｗ∞Ｌ（s）およびＷ∞Ｈ（s）の周波数特性をそれぞれ実線で示す図である。前述のとおり、次式（２３）となるように、両帯域の周波数重み関数Ｗ∞Ｌ（s）およびＷ∞Ｈ（s）を与えた上で、ロバスト安定条件（次式（２４））を満たすように制御系を設計することになる。ここで、Ｔｕｗ（s）は図２０においてｗｕからｕまでの伝達関数行列である。そして、周波数重み関数Ｗ∞Ｌ（s）およびＷ∞Ｈ（s）をプロパーな伝達関数として与え、それらが次式（２５）を満たすように、すなわち１入出力系ではＷ∞Ｌ（s）およびＷ∞Ｈ（s）の周波数応答がそれぞれΔ∞Ｌ（s）およびΔ∞Ｈ（s）の周波数応答よりも全周波数帯域にわたって大きくなるように設定する。
【００５２】
【数１６】
2008257348.tif 000018
次に、コントローラの設計方法について説明する。上述のように両帯域の周波数重み関数Ｗ∞Ｌ（s）およびＷ∞Ｈ（s）を与えた上で、以下の混合Ｈ２／Ｈ∞制御問題（次式（２６）、（２７））を解くことによりコントローラを設計する。
【００５３】
【数１７】
2008257348.tif 000019
ここで、ＴＺ２ｗ（s）は外乱ｗから周波数重み付き制御応答ｚ２までの伝達関数行列、ＴＺ１ｗ（s）は外乱ｗから周波数重み付き制御入力ｚ１までの伝達関数行列であり、式（２６）はＴＺ２ｗ（s）のＨ２ノルム、式（２７）の左辺はＴＺ１ｗ（s）のＨ∞ノルムを表す。制御応答にかかる周波数重み関数Ｗ２（s）は、モード座標から成る状態ベクトルｑｅの成分に直接かかる重み係数として与える。こうすることにより、コントローラの次数を増加させることなく着目するモード（共振ピーク）を効果的に低減することができる。コントローラは、前述のとおり上記制御問題を線形行列不等式（ＬＭＩ）の形で記述した上で、数値探索法により設計される。設計されたコントローラは、次式（２８）の状態方程式として記述される。
【００５４】
【数１８】
2008257348.tif 000020
ここで、ｑｃｐはコントローラの状態ベクトル、Ａｃｐ、Ｂｃｐ、ＣｃｐおよびＤｃｐはコントローラのシステム行列であり、式（１２）は観測出力ｙから制御入力ｕの特性を表現したものである。行列Ａｃｐ、Ｂｃｐ、ＣｃｐおよびＤｃｐが制御系設計により決定されることになる。
【００５５】
図２２〜図２５は、シミュレーションおよび実験結果を示す図である。図２２および図２３は低周波数帯域におけるシミュレーションおよび実験結果、図２４および図２５は高周波数帯域におけるシミュレーションおよび実験結果を示す。シミュレーションに関しては、実験モード解析により同定されたモデルに対し、設計されたコントローラの出力を制御入力とする閉ループ系の拡大状態方程式を記述し、それに対して周波数応答を計算している。ここで、シミュレーションにおける実システムとして、低周波数帯域と高周波数帯域の両帯域に含まれるモードを採用した大次元の制御対象モデルを作成しておき、それに対して閉ループ系を構成することにより、対象としている全周波数帯域における閉ループ系の周波数応答解析が行われる。図２２〜図２５から、シミュレーション結果と実験結果はよく対応していることが分かる。
【００５６】
以上のようにして設計されたコントローラによれば、低周波数帯域での実験結果において、目標としていた４００Ｈｚ付近のピークが約１０ｄＢ低減している。また、図２０のようにロバスト安定性を考慮しているため、モデル化していない高周波数帯域でも発散せずに非制御時の応答にほぼ重なっている。高周波数帯域での実験結果に関しては、１６００Ｈｚ付近のピークが約１０ｄＢ、２３００Ｈｚ付近のピークが約５ｄＢ低減されている。先程と同様、図２１のようにロバスト安定性を考慮しているため、モデル化していない低周波数帯域でも発散せずに非制御時と同様の応答を示している。両帯域ともシミュレーション結果と実験結果が良く一致しており、各周波数帯域におけるモデル化および制御系設計が適切になされていることが分かる。
【００５７】
次に、各帯域用の複数コントローラをロバスト協調設計することにより構成される振動低減装置を説明する。
【００５８】
図２６は、複数のコントローラを備えた振動低減装置の実施例を図５に対応させて示す構成図である。この振動低減装置２０の基本構成は、前述の図５と同一であるが、図５との違いは、制御対象１１を具体的なエンジン・エンジェクタシステムとして与えていることと、外力に関する外乱をインパルスハンマーにより作用させていることである。
【００５９】
先ず、コントローラのロバスト協調設計方法と結果について説明する。図１１に示すステップにしたがって、コントローラのロバスト協調設計を行う。まず、ステップ１より、低周波数帯域モデル誤差Δ∞Ｌ（s）を求め、低周波数帯域周波数重み関数Ｗ∞Ｌ（s）を作成し、低周波数帯域コントローラＫＬ（s）を設計する。
【００６０】
図２７は、Δ∞Ｌ（s）とＷ∞Ｌ（s）の周波数特性を示す図である。そして、ステップ２において低周波数帯域コントローラＫＬ（s）でフィードバックした状態の周波数応答を設計したコントローラと計測で得られた周波数応答から計算し、それを実システムの応答とみなした上で高周波数帯域制御対象ＰＨ（s）とのモデル誤差Δ∞Ｈ（s）を求め、高周波数帯域周波数重み関数Ｗ∞Ｈ（s）を設計する。図９から、低周波数帯域コントローラＫＬ（s）でフィードバックした状態の高周波数帯域の実制御対象Ｐ~Ｈ（s）と高周波数帯域モデル誤差Δ∞Ｈ（s）は次式（２９）、（３０）で表せる。
【００６１】
【数１９】
2008257348.tif 000021
ここで、Ｐ~（s）はノミナルモデルであり、その周波数特性は実験で計測された周波数応答に対応する。
【００６２】
図２８は、ＰＨ（s）とＰ~Ｈ（s）を用いて作成した高周波数帯域モデル誤差Δ∞Ｈ（s）の周波数特性を低周波数帯域コントローラによる制御時および非制御時で比較したものを示す図である。図２８からわかるように低周波数帯域制御により高周波数帯域で誤差が増加することはなく、低周波数帯域で誤差が減少していることが分かる。高周波数帯域で誤差が増加しないことは、ロバスト安定性を満たすように低周波数帯域周波数重み関数Ｗ∞Ｌ（s）を与えて制御系を設計したことから理解できる。
【００６３】
図２９は、フィードバックした状態の実制御対象Ｐ~Ｈ（s）を用いて作成した高周波数帯域モデル誤差Δ∞Ｈ（s）に対してロバスト安定性を満たすように作成した高周波数帯域周波数重み関数Ｗ∞Ｈ（s）を示す図である。そして、高周波数帯域コントローラＫＨ（s）を設計する。次に、ステップ３において、高周波数帯域コントローラＫＨ（s）でフィードバックした状態の低周波数帯域モデル誤差Δ∞Ｌ（s）を求める。図７から、高周波数帯域コントローラＫＨ（s）でフィードバックした状態の低周波数帯域の実制御対象Ｐ~Ｈ（s）と低周波数帯域モデル誤差Δ∞Ｌ（s）は次式（３１）、（３２）で表せる。
【００６４】
【数２０】
2008257348.tif 000022
ステップ３で求めた低周波数帯域モデル誤差Δ∞Ｌ（s）に対し、ステップ１で作成した低周波数帯域周波数重み関数Ｗ∞Ｌ（s）がロバスト安定を満たしていれば終了である。
【００６５】
図３０は、ステップ１で求めた低周波数帯域モデル誤差Δ∞Ｌ（s）、ステップ３で再度求めた低周波数帯域モデル誤差Δ∞Ｌ（s）、およびステップ１で作成した低周波数帯域周波数重み関数Ｗ∞Ｌ（s）を比較した図である。モデル誤差は２３００Ｈｚ付近で悪化しているが、制御帯域である０〜８００Ｈｚ付近ではほぼ重なっている。これは、図２９において高周波数帯域周波数重み関数Ｗ∞Ｈ（s）が制御帯域外で十分にロバスト性を確保しているからである。ステップ１で作成した低周波数帯域周波数重み関数Ｗ∞Ｌ（s）はステップ３で再度求めた低周波数帯域モデル誤差Δ∞Ｌ（s）に対しロバスト安定を満たしている。したがって、ステップ１で作成した低周波数帯域コントローラＫＬ（s）は閉ループでもロバスト性をもつので設計終了である。
【００６６】
図３１および図３２は、上述のとおり設計されたコントローラを用いて、閉ループ系のシミュレーションおよび実験を行った結果を示す図である。なお、前述の単一コントローラの場合と同様、シミュレーションでは低周波数帯域と高周波数帯域の両帯域に含まれるモードを採用した大次元の制御対象モデルに対し閉ループ系を構成している。図３１および図３２から、シミュレーション結果と実験結果は良い対応を示していることが分かる。
【００６７】
以上のように、各帯域用の複数コントローラを協調設計する方法によれば、図３１および図３２の結果から、協調設計されたコントローラにより４００Ｈｚおよび１６００Ｈｚ付近のピークにおいて約１０ｄＢ、２３００Ｈｚ付近のピークにおいて約２ｄＢ低減させることができた。単一コントローラの場合と同様、制御帯域において効果的に制御しつつ、制御帯域外で応答が悪化していないことが分かる。特に、本制御手法により低周波数帯域と高周波数帯域の振動を同時に低減していることが最大の特徴である。ここで用いたコントローラは、各帯域用の低次元コントローラを並列に結合したものであり、広帯域用の大次元モデルから直接設計した単一のコントローラよりも計算効率を著しく向上させているため、本結果の達成につながっている。この根拠は、前述のとおり単一コントローラと並列コントローラの演算回数の比較から理解できる。本実施例のように約３ｋＨｚまでの振動の制御が要求される場合、一般にデジタル制御システムにおいて約３０ｋＨｚのサンプリング周波数が必要とされるが、広帯域用の大次元モデルから直接設計した単一の制御器で実現できるサンプリング周波数は１０ｋＨｚ程度であったのに対し、本手法を用いた場合は３０ｋＨｚが達成された。以上のことは、本ロバスト協調設計法の優れた有効性を示していることに他ならない。
【００６８】
次に、比較のため、各帯域用の複数コントローラを非ロバスト協調設計することにより構成される振動低減装置を説明する。
【００６９】
図３３は、複数のコントローラを備えた振動低減装置の実施例を図５に対応させて示す構成図である。この振動低減装置２０の基本構成は、前述の図５と同一であるが、図５との違いは、制御対象１１を具体的なエンジン・エンジェクタシステムとして与えていることと、外力に関する外乱をインパルスハンマーにより作用させていることである。
【００７０】
先ず、コントローラのロバスト協調設計方法と結果について説明する。図１９にて設計した低周波数帯域および高周波数帯域ごとの単一コントローラＫＬ（s）およびＫＨ（s）を並列に結合し、閉ループ系を構成する。したがって、両コントローラの干渉を考慮した協調設計はされておらず、ロバスト安定条件は保証されていない。
【００７１】
図３４は、非ロバスト協調設計による複数コントローラを用いた場合の実験結果を示す図である。図３４の結果から１６００Ｈｚ付近のピークで約３ｄＢ、２３００Ｈｚ付近のピークで約４ｄＢ低減できているが、４００Ｈｚ付近のピークでは約２ｄＢ増大している。また、１８００Ｈｚ付近に新たな大きなピークが存在しているのが分かる。
【００７２】
次に、ロバスト協調設計による複数コントローラの実験結果との比較について説明する。図３２および図３４より、ロバスト協調設計による複数コントローラを用いた場合は低周波数帯域、高周波数帯域ともにロバスト安定を保ちながら制振効果を得ているのに対し、非ロバスト協調設計による複数コントローラを用いた場合は制振効果にばらつきがあり、ロバスト安定が保たれていないのが分かる。これは、コントローラを並列に接続する際に互いの出力による干渉を考慮していないため、ロバスト安定性を満たしていないからであると考察できる。この比較から、コントローラを複数用いる場合、上記のロバスト協調設計を用いることが有効であると言える。
【００７３】
なお、上述した実施形態では、コントローラの数を２つとして説明したが、コントローラの数は基本的には２以上の任意の数で効果が期待できる。その根拠は、以下のとおりである。いま、単一のコントローラの次数をｎとし、それをｍ個のサブシステムに等分割できたとする。コントローラ演算における乗算回数は前者で約ｎ２、後者で約ｎ２／ｍとなり、ｍ≧２で後者が前者よりも必ず小さくなる。しかし、コントローラの数が多くなるにつれて制御系の協調設計の手間は大きくなる。すなわち、周波数重み関数を修正する回数は必然的に大きくなる。しかし、コントローラの分割数が多いほど、各コントローラの設計効率は高いという性質は不変である。以上のことから、コントローラの数は２個〜５個程度が適切であると考えられる。
【００７４】
また、本実施形態では、１軸方向の制御のみについて示したが、３次元方向の制御に拡張することは可能である。この場合、センサおよびアクチュエータの設置位置および方向に関して、３軸全ての方向について応答を検知できかつ駆動できること、もしくはこれを満たすように複数のセンサおよびアクチュエータを使用することが条件となる。これは、可制御性および可観測性という制御系に要求される基本的条件を満たすことに他ならない。これらの条件を満たした上で、システムのモデル化および制御系設計を行えばよいことになる。ここで、センサおよびアクチュエータの数が複数になった場合でも、各アクチュエータからセンサ出力までの周波数応答から、実験モード解析によりシステムがモデル化され、前述と同じプロセスで制御系が設計される。
【００７５】
また、本実施形態では、インパルスハンマーにより作用させた外乱（外部加振力）の数は１であったが、２以上の外乱が作用する場合でも制御系の設計は可能である。この場合、各外乱からセンサ出力までの周波数応答を計測するか、もしくは無相関な外乱を同時に作用させたときのセンサ出力までの周波数応答を計測すれば、実験モード解析によりモード特性が同定されることになる。そして、前述と同じプロセスにしたがって、同定されたモード特性から状態方程式を記述し制御系を設計すればよいことになる。
【００７６】
また、本実施形態の振動低減装置の適用対象としては、ディーゼルおよびガソリンエンジンの振動、自動車の車体など自動車全般に関わる振動、鉄道車両の振動、航空機の振動、磁気ディスク装置や光ディスク装置、ならびにプリンタ機器をはじめとする情報機器の振動などがある。
【００７７】
また、本実施形態は、混合Ｈ２／Ｈ∞制御方法に適用できるほか、Ｈ∞制御方法にも適用できる。すなわち、Ｈ∞制御の場合、制御応答の制御量とロバスト安定性の制御量を併せてＨ∞ノルムで評価するが、ロバスト安定性の評価指標は本発明と共通であるため、このＨ∞ノルムを評価関数としてとして与えるだけで、本発明と同じ設計方法によりＨ∞制御に対するコントローラのロバスト協調設計が可能である。
【００７８】
また、本実施形態は、音の特性をモデル化することにより、騒音低減制御に適用することができる。この場合、センサとしてマイクロホンを使用し、振動実験で計測される音圧の周波数応答を用いてカーブフィッティングすれば、入力から音圧までのモデルが同定される。そして、この音圧を制御応答として制御問題を定義すれば、本発明と同一の方法により、コントローラの協調設計が行われる。単一のコントローラを用いる場合において、システム同定による音のモデル化および騒音低減制御が可能なことは以下の論文で実証済みである。
（ａ）梶原逸朗、福田将直、下嶋浩、モード解析と線形行列不等式に基づくフィードバック音響制御、日本機械学会論文集（Ｃ編）、６４巻、６２１号、１９９８、ｐｐ．１６６８−１６７５
（ｂ）下嶋浩、松永義弘、小池関也、梶原逸朗、モード解析に基づく複数の制御音源を用いたフィードバック音響制御、日本機械学会論文集（Ｃ編）、６５巻、６３３号、１９９９、ｐｐ．１８４９−１８５６
論文（ａ）では，システムを比例粘性減衰系としてカーブフィッティングによるモデル化を行い、Ｈ∞制御方法および混合Ｈ２／Ｈ∞制御方法によりコントローラを設計し、マイクロホンとスピーカーを用いたフィードバック制御により騒音低減を実現している。論文（ｂ）では、複数のスピーカー（制御音源）をもつシステムを一般粘性減衰系としてカーブフィッティングによるモデル化を行い、論文（ａ）と同様の制御系設計およびフィードバック制御実験により優れた騒音低減効果を達成している。これらの論文では、単一の周波数帯域における単一のモデル化およびコントローラ設計法を示しているが、複数の周波数帯域に分割して考えることにより、本発明の方法をそのまま適用することができる。
【００７９】
また、本実施形態におけるモデル化では、式（１７）〜式（２１）に示すように比例粘性減衰系を扱ったが、系を一般粘性減衰系としたときのモデル化も同様に可能であり、本発明をそのまま適用できる。一般粘性減衰系におけるモード特性からシステムを状態方程式で記述する方法は上記（６）の論文（ｂ）に示されている。
【００８０】
また、本実施形態におけるモデル誤差として、乗法的誤差と同様に加法的誤差を扱うことも可能である。モデル誤差を加法的誤差として扱う場合には、システムの観測出力に作用する外乱から制御入力までの伝達関数に対し、加法的誤差特性を指標とする周波数重み関数を与えることにより、本実施形態をそのまま適用することができる。
【図面の簡単な説明】
【００８１】
【図１】単一のコントローラを備えた振動低減装置を示す構成図である。
【図２】図１の振動低減装置の制御系の概略の第１のブロック線図である。
【図３】図１の振動低減装置の制御系の概略の第２のブロック線図である。
【図４】図１の振動低減装置のモデル誤差を乗法的誤差として与え、外乱の作用する位置とあわせて表記したブロック線図である。
【図５】本発明の一実施の形態に係る複数のコントローラを備えた振動低減装置を示す構成図である。
【図６】図５の振動低減装置の低周波数帯域用と高周波数帯域用の二つのコントローラを並列に結合した制御系のブロック線図である。
【図７】図５の振動低減装置の低周波数帯域コントローラのモデル誤差を乗法的誤差として与え、外乱の作用する位置とあわせて表記したブロック線図である。
【図８】図５の振動低減装置の低周波数帯域コントローラでフィードバックしたブロック線図である。
【図９】図５の振動低減装置の高周波数帯域コントローラのモデル誤差を乗法的誤差として与え、外乱の作用する位置とあわせて表記したブロック線図である。
【図１０】図５の振動低減装置の高周波数帯域コントローラでフィードバックしたブロック線図である。
【図１１】本発明の一実施の形態に係る設計方法のステップ１からステップ４をフロー図にしたものである。
【図１２】単一のコントローラを備えた振動低減装置の実施例を図１に対応させて示す構成図である。
【図１３】図１２の振動低減装置の計測された外乱加振入力から加速度センサの出力信号までの周波数応答（アクセレランス）を示す図である。
【図１４】図１２の振動低減装置の計測された制御入力外乱信号から加速度センサの出力信号までの周波数応答（アクセレランス）を示す図である。
【図１５】図１２の振動低減装置の振動実験により計測したデータを用いて、実験モード解析によりカーブフィットした結果（アクセレランス）を示す第１の図である。
【図１６】図１２の振動低減装置の振動実験により計測したデータを用いて、実験モード解析によりカーブフィットした結果（アクセレランス）を示す第２の図である。
【図１７】図１２の振動低減装置の振動実験により計測したデータを用いて、実験モード解析によりカーブフィットした結果（アクセレランス）を示す第３の図である。
【図１８】図１２の振動低減装置の振動実験により計測したデータを用いて、実験モード解析によりカーブフィットした結果（アクセレランス）を示す第４の図である。
【図１９】図１２の振動低減装置のロバスト安定性を考慮した周波数重みの設定について示すブロック線図である。
【図２０】図１２の振動低減装置の低周波数帯域および高周波数帯域において、両帯域の乗法的誤差Δ∞Ｌ（s）およびΔ∞Ｈ（s）を計算した結果をそれぞれ破線で示し、両帯域において与えた周波数重み関数Ｗ∞Ｌ（s）およびＷ∞Ｈ（s）の周波数特性をそれぞれ実線で示す第１の図である。
【図２１】図１２の振動低減装置の低周波数帯域および高周波数帯域において、両帯域の乗法的誤差Δ∞Ｌ（s）およびΔ∞Ｈ（s）を計算した結果をそれぞれ破線で示し、両帯域において与えた周波数重み関数Ｗ∞Ｌ（s）およびＷ∞Ｈ（s）の周波数特性をそれぞれ実線で示す第１の図である。
【図２２】図１２の振動低減装置のシミュレーションおよび実験結果を示す第１の図である。
【図２３】図１２の振動低減装置のシミュレーションおよび実験結果を示す第２の図である。
【図２４】図１２の振動低減装置のシミュレーションおよび実験結果を示す第３の図である。
【図２５】図１２の振動低減装置のシミュレーションおよび実験結果を示す第４の図である。
【図２６】複数のコントローラを備えた振動低減装置の実施例を図５に対応させて示す構成図である。
【図２７】図２６の振動低減装置のΔ∞Ｌ（s）とＷ∞Ｌ（s）の周波数特性を示す図である。
【図２８】図２６の振動低減装置のＰＨ（s）とＰ~Ｈ（s）を用いて作成した高周波数帯域モデル誤差Δ∞Ｈ（s）の周波数特性を低周波数帯域コントローラによる制御時および非制御時で比較したものを示す図である。
【図２９】図２６の振動低減装置のフィードバックした状態の実制御対象Ｐ~Ｈ（s）を用いて作成した高周波数帯域モデル誤差Δ∞Ｈ（s）に対してロバスト安定性を満たすように作成した高周波数帯域周波数重み関数Ｗ∞Ｈ（s）を示す図である。
【図３０】図１１のフロー図のステップ１で求めた低周波数帯域モデル誤差Δ∞Ｌ（s）、ステップ３で再度求めた低周波数帯域モデル誤差Δ∞Ｌ（s）、およびステップ１で作成した低周波数帯域周波数重み関数Ｗ∞Ｌ（s）を比較した図である。
【図３１】図１１のフロー図により設計されたコントローラを用いて、閉ループ系のシミュレーションおよび実験を行った結果を示す第１の図である。
【図３２】図１１のフロー図により設計されたコントローラを用いて、閉ループ系のシミュレーションおよび実験を行った結果を示す第２の図である。
【図３３】複数のコントローラを備えた振動低減装置の実施例を図５に対応させて示す構成図である。
【図３４】図３３の振動低減装置の非ロバスト協調設計による複数コントローラを用いた場合の実験結果を示す図である。
【図３５】制御問題を表すブロック線図である。
【符号の説明】
【００８２】
１０、２０  振動低減装置、１１  制御対象、１２  加速度センサ、１３  ローパスフィルター、１４  デジタル制御システム、１４ａ  ＡＤ変換器、１４ｂ  コントローラ、１４ｃ  ＤＡ変換器、１５  アクチュエータ、１６  ＦＦＴアナライザ、１７  センサーアンプ、１８  増幅器、２１  低周波数帯域コントローラ、２２  高周波数帯域コントローラ
【図１】
2008257348.tif 000023
【図２】
2008257348.tif 000024
【図３】
2008257348.tif 000025
【図４】
2008257348.tif 000026
【図５】
2008257348.tif 000027
【図６】
2008257348.tif 000028
【図７】
2008257348.tif 000029
【図８】
2008257348.tif 000030
【図９】
2008257348.tif 000031
【図１０】
2008257348.tif 000032
【図１１】
2008257348.tif 000033
【図１２】
2008257348.tif 000034
【図１３】
2008257348.tif 000035
【図１４】
2008257348.tif 000036
【図１５】
2008257348.tif 000037
【図１６】
2008257348.tif 000038
【図１７】
2008257348.tif 000039
【図１８】
2008257348.tif 000040
【図１９】
2008257348.tif 000041
【図２０】
2008257348.tif 000042
【図２１】
2008257348.tif 000043
【図２２】
2008257348.tif 000044
【図２３】
2008257348.tif 000045
【図２４】
2008257348.tif 000046
【図２５】
2008257348.tif 000047
【図２６】
2008257348.tif 000048
【図２７】
2008257348.tif 000049
【図２８】
2008257348.tif 000050
【図２９】
2008257348.tif 000051
【図３０】
2008257348.tif 000052
【図３１】
2008257348.tif 000053
【図３２】
2008257348.tif 000054
【図３３】
2008257348.tif 000055
【図３４】
2008257348.tif 000056
【図３５】
2008257348.tif 000057