(57)【特許請求の範囲】
【請求項１】
制御対象の全入出力パワーの収支の式が次式で近似される制御システムにおいて，
0004609767.tif 000078
（ここで，ｄ，ｅ，ｑ，ｕ，ｖ，ｚ∈Ｒｎ，Ｍ∈Ｒｎ×ｎは正定対称な慣性マトリクス，ｎは制御対象の自由度である．ｑは一般化座標，ｕは制御則，ｖは力入力の外乱，ｚは変位入力の外乱である．ｚはばねを介して慣性に作用するものとする．ｄはコリオリ力や遠心力やダンピング力など，ｅはポテンシャル力である．また制御対象は受動要素だけとは限らないため，内部にエネルギ源がありこれが運動に影響を与えていれば外乱入力として扱う．制御アクチュエータの数が制御対象の自由度ｎより少ない場合は、欠落したアクチュエータに対応するベクトルｕの要素をゼロとする．また一つのアクチュエータが複数の慣性に作用する場合には、対応するベクトルｕの要素にその条件を含めておく．変数ｑやｚの上に付された記号“・”とその数は時間微分とその回数を表す．）
次式であらわされる評価関数を最小にする制御方法．
0004609767.tif 000079
ここで、ｇはｑとｑの１階、２階の導関数に関する連続関数であって制御対象の制御性能を表現した関数であり、ｒは重み係数である．
【請求項２】
請求項１の制御システムにおいて、次のスカラー関数Ｌを定義し、
0004609767.tif 000080
（ここで、κは定数である．）
制御則ｕは次式から導かれた関係式を含むことを特徴とする制御方法．
0004609767.tif 000081
【請求項３】
請求項１において、制御則ｕが次式の各項の一部あるいは全部を含む制御方法．
0004609767.tif 000082
（ここで、κは定数である．）
【発明の詳細な説明】
【技術分野】
【０００１】  本発明は、主として機械力学系システムの実時間最適制御方法に関する．線形、非線形、ホロノミック、非ホロノミック系を問わず広範囲のシステムに適用できる．
【背景技術】
【０００２】  ポントリャーギンが随伴変数を導入して最適制御問題をハミルトン形式に書き換えて最大原理を見出したことは，その後の最適制御理論の発展に大きな影響を与えた．しかし，現在でも随伴変数と状態変数についての連立微分方程式を解かねばならない状況は変わらず，特に非線形システムの最適制御則を解析的に導くことは容易ではない．そのため様々な数値探索的方法や推論的手法が工夫されてはいる．しかし制御則が数値解の場合は、制御系設計者が得られた制御則を改善しようとしても、制御対象システムのパラメータとの対応関係がわからず試行錯誤的にならざるを得ない．このような状況から脱却するためには見通しの良い解析的手法が必要である．
【発明が解決しようとする課題】
【０００３】  このような背景から，本考案では，従来とは異なる手法で主として機械力学系システムの非線形系を含む比較的広範囲の最適制御問題に対して制御則を解析的に導く手法を提示する．制御対象のエネルギ収支に着目した受動性の概念によりシステムを漸近安定させる制御則が簡単に導かれるように（例えば非特許文献１参照），本考案でも制御対象のエネルギ収支に着目する．入力に関する評価関数を制御装置から制御対象に伝達されるエネルギとし，これに制御性能を表現する関数を加えて全体の評価関数とする．さらに付帯条件として制御対象の完全なエネルギ収支の式を加えて汎関数とする．こうすると汎関数は制御変数に関して１次式となるので，制御対象が線形，非線形を問わず汎関数を最小化する条件から微分方程式を解かずに最適制御則を解析的に導くことができる．
この汎関数を最小化する条件は，最適制御であるための必要条件であり，しかも未定定数を含む弱い条件であるが，制御対象の構造的性質をうまく活用することにより有益な制御則が得られる．後の事例によりその有用性を示す．
【０００４】
最適制御問題を定式化する．まず制御対象の特性表現について述べる．制御理論では制御対象の特性を運動方程式で記述するが，ここでは制御対象の全入出力パワーの収支の式を用いる．具体的には，システムの各自由度毎の運動方程式をベクトル表示しこれに速度ベクトルを乗じたものである．入力パワーは制御による入力だけでなく外乱入力も含まれ，また制御対象は受動要素だけとは限らないため，内部にエネルギ源がありこれが運動に影響を与えていれば外乱入力として扱う．
【式１】
0004609767.tif 000002
ここで，ｄ，ｅ，ｑ，ｕ，ｖ，ｚ∈Ｒｎ，Ｍ∈Ｒｎ×ｎ  は正定対称な慣性マトリクス，ｎは制御対象の自由度である．ｑは一般化座標，ｕは制御則，ｖは力入力の外乱，ｚは変位入力の外乱である．ｄはコリオリ力や遠心力やダンピング力など，ｅはポテンシャル力である．このような一般化座標ｑの説明やこれを用いた機械系の運動方程式については、例えば非特許文献２のＰ１１−１２や非特許文献３のｐ５４−５５に解説されている．
制御対象が非線形であっても制御装置を合理的に設計すれば，式（１）のようにｕを直接Ｍに作用させることができる．本論文では，このような合理的に設計された機械力学系システムを想定して最適制御則を導く．
上記のシステムに対し，次の評価関数を考える．
【式２】
0004609767.tif 000003
ここでｇは制御性能の評価を与えるスカラー関数であり，コスト関数の被積分項
0004609767.tif 000004
式（２）を最小化する制御ｕ（ｔ）を求めることが最適制御の課題である．
【非特許文献１】有本卓，新版ロボットの力学と制御，（２００２），ＰＰ．６８−６９．朝倉書店．
【非特許文献２】大貫義郎，解析力学，ＰＰ．１１−１２．朝倉書店．
【非特許文献３】有本卓，新版ロボットの力学と制御，（２００２），ＰＰ．５４−５５．朝倉書店．
【式１】
0004609767.tif 000005
ここで，ｄ，ｅ，ｑ，ｕ，ｖ，ｚ  ∈Ｒｎ  ，Ｍ  ∈Ｒｎ×ｎ  は正定対称な慣性マトリクス，ｎは制御対象の自由度である．ｑは一般化座標，ｕは制御入力，ｖは力入力の外乱，ｚは変位入力の外乱である．ｄはコリオリ力や遠心力やダンピング力など，ｅはポテンシャル力である．
制御対象が非線形であっても制御装置を合理的に設計すれば，式（１）のようにｕを直接Ｍに作用させることができる．本論文では，このような合理的に設計された機械力学系システムを想定して最適制御則を導く．
上記のシステムに対し，次の評価関数を考える．
【式２】
0004609767.tif 000006
0004609767.tif 000007
置が御対象に加えるパワーである．ｒは重み係数である．
式（２）を最小化する制御ｕ（ｔ）を求めることが最適制御の課題である．
【非特許文献１】有本卓，新版ロボットの力学と制御，（２００２），ＰＰ．６８−６９．朝倉書店．
【課題を解決するための手段】
【０００５】
従来の最適制御理論では，随伴変数を導入して微分方程式を導きこれに元の状態方程式を連立させて解くことになるが，本論文ではこのような方法とは全く異なる手法を示す．最適制御の必要条件を求めるため次のスカラー関数Ｌを定義する．
【式３】
0004609767.tif 000008
ここで，κは未定定数である．右辺の｛｝内は，式（１）の左辺と同じで制御対象の全パワー収支であるからエネルギ保存則を満たし常にゼロである．従って式（３）で表されるＬの積分を最小化する条件は，式（２）も最小化する．
従って関数Ｌにｑを変数とする変分原理を適用した次式はｕが最適制御であるための必要条件を与える．Ｌはｕに関して１次式であるから，∂Ｌ／∂ｕは意味がなく次式に制御に関するすべての情報が含まれる．
【式４】
0004609767.tif 000009
Ｌにはｑの２階の導関数が含まれるため一般的なオイラーの方程式に第３項が追加されている（例えば非特許文献４参照）．
式（４）を積分し積分定数をゼロとすると次式になる．
【式５】
0004609767.tif 000010
式（３）を式（５）に代入し左辺第１〜３項の順に記すと次のようになる．
【式６】
0004609767.tif 000011
（式５）はＬをｑおよびその導関数で偏微分したものであるが、注目されるのはＬの中
0004609767.tif 000012
ることである．このｕを左辺にもってきて他を右辺に移項すれば（式７）が得られる．
【式７】
0004609767.tif 000013
このように，本手法はＪを最小化するｑとｕの関係を直接導くことができるため，従来のように２点境界値問題を最適性の原理を用いて解くプロセスが不要になる．κは未定定数であるが，κ＝０の時はｕが制御の評価関数のパラメータのみで定まり，κ＝∞の時はｕが制御対象のパラメータのみで定まることからｕが最適であるためにはκはゼロでない有限値でなければならないことは明らかである．
上式の第１行は外力ｖと慣性のｑ依存性に対する制御，第２行はコリオリ力や遠心力やダンピング力に対する制御，第３行はポテンシャル力とそのｑ依存性および外力ｚに対する制御，第４行は評価関数を低減させる制御でありそれぞれ意味が明確である．
式（７）には未実行の微積分項が含まれているが，全ての外力と状態量の検出あるいは推定が可能とすればこれらの実時間での実行は可能である．式（５）において積分定数をゼロとしたが，上記の結果より積分定数は制御則に一定のバイアスを与えることになるためゼロとすることが妥当であることがわかる．式（７）中の積分も同様である．
式（７）ではアクチュエータの数と系の自由度が同じであることを想定しているが，アクチュエータの数が少ない場合には次のような処置が必要になる．例えばアクチュエータが二つの独立な慣性Ｍｉ，Ｍｉ＋１の間に置かれるような場合は，制御ベクトルｕにその拘束条件を含めておき最適制御則は二つの制御則を重みを付けて加算したものにする．
【式８】
ｕｏｐｔ，ｉ，ｉ＋１＝ρｉｕｉ＋ρｉ＋１ｕｉ＋１
ここでｕＩ，ｕｉ＋１はそれぞれが独立な制御として導いた場合の最適制御則であり，ρｉ，ρｉ＋１は重み係数である．これらの重み係数や未定定数κの決定法については，次の具体例で示す．
【非特許文献４】Ｌ．Ｅ．エルスゴルツ，瀬川富士訳，科学者・技術者のための変分法，（１９７２），ｐｐ．３４−３７．ブレイン図書出版．
【式３】
0004609767.tif 000014
ここで，κは未定定数である．右辺の｛｝内は，式（１）の左辺と同じで制御対象の全パワー収支であるからエネルギ保存則を満たし常にゼロである．従って式（３）で表されるＬの積分を最小化する条件は，式（２）も最小化する．
従って関数Ｌにｑを変数とする変分原理を適用した次式はｕが最適制御であるための必要条件を与える．Ｌはｕに関して１次式であるから，∂Ｌ／∂ｕは意味がなく次式に制御に関するすべての情報が含まれる．
【式４】
0004609767.tif 000015
Ｌにはｑの２階の導関数が含まれるため一般的なオイラーの方程式に第３項が追加されている（例えば非特許文献２参照）．
式（４）を積分し積分定数をゼロとすると次式になる．
【式５】
0004609767.tif 000016
式（３）を式（５）に代入し左辺第１〜３項の順に記すと次のようになる．
【式６】
0004609767.tif 000017
上式から制御則が次のように求まる．
【式７】
0004609767.tif 000018
このように，本手法はＪを最小化するｑとｕの関係を直接導くことができるため，従来のように２点境界値問題を最適性の原理を用いて解くプロセスが不要になる．κは未定定数であるが，κ＝０の時はｕが制御の評価関数のパラメータのみで定まり，κ＝∞の時はｕが制御対象のパラメータのみで定まることからｕが最適であるためにはκはゼロでない有限値でなければならないことは明らかである．
上式の第１行は外力ｖと慣性のｑ依存性に対する制御，第２行はコリオリ力や遠心力やダンピング力に対する制御，第３行はポテンシャル力とそのｑ依存性および外力ｚに対する制御，第４行は評価関数を低減させる制御でありそれぞれ意味が明確である．式（７）には未実行の微積分項が含まれているが，全ての外力と状態量の検出あるいは推定が可能とすればこれらの実時間での実行は可能である．式（５）において積分定数をゼロとしたが，上記の結果より積分定数は制御則に一定のバイアスを与えることになるためゼロとすることが妥当であることがわかる．式（７）中の積分も同様である．
式（７）ではアクチュエータの数と系の自由度が同じであることを想定しているが，アクチュエータの数が少ない場合には次のような処置が必要になる．例えばアクチュエータが二つの独立な慣性Ｍｉ，Ｍｉ＋１の間に置かれるような場合は，制御ベクトルｕにその拘束条件を含めておき最適制御則は二つの制御則を重みを付けて加算したものにする．
【式８】
0004609767.tif 000019
ここでｕｉ，ｕｉ＋１はそれぞれが独立な制御として導いた場合の最適制御則であり，ρｉ，ρｉ＋１は重み係数である．これらの重み係数や未定定数κの決定法については，次の具体例で示す．
【非特許文献２】Ｌ．Ｅ．エルスゴルツ，瀬川富士訳，科学者・技術者のための変分法，（１９７２），ｐｐ．３４−３７．ブレイン図書出版．
【発明の効果】
【０００６】  ここでは，従来の最適制御理論では解析解を得ることができない問題に本手法を適用しその有効性を明らかにする．最初は評価関数が２次形式にならないアクティブサスペンションのエネルギ回生制御を扱い，次に未だに最適制御則の導出が困難とされている非ホロノミック系の制御への適用を試みる．
【０００７】  まず、アクティブサスペンションのエネルギ回生制御への適用について具体的な定式化を行う．特許文献１では、図１に示すようなアクティブサスペンションのモデルを用いてエネルギ回生制御を考案し，近似解法により準最適制御則ｕ＊ｏｐｔを導いている．この事例では結果的に入力に関する評価関数が本論文の手法と同じく制御装置から制御対象に伝達されるエネルギ収支式となっているが，当時はまだ制御対象のエネルギ収支に着目した最適制御という認識はなく，従来の制御理論では解析的に解けないケースとして近似解法を試みたものであり，近似解法ながらほぼ最適制御であることを立証している．本考案では解析的に導いた最適制御則ｕｏｐｔが簡潔な表現ながらこの近似解ｕ＊ｏｐｔと同等の性能であることを確認する．
まず特許文献１の評価関数を変数ｑを用いて記述する．
【式９】
0004609767.tif 000020
ここで，ｒ１，ｒは重み係数である．重みｒ１はばね下質量の制振が従来車両と同程度になるようにｒ１＝５×１０−３  とし，ｒについては振動の評価とエネルギの評価が同程度の値になるようにｒ＝１０−５とした．
【０００８】  続いて制御則を導出する．まず図１のモデルから式（１）の各要素を求めると次のようになる．
【式１０】
0004609767.tif 000021
【式１１】
0004609767.tif 000022
【式１２】
0004609767.tif 000023
【式１３】
0004609767.tif 000024
【式１４】
0004609767.tif 000025
【式１５】
0004609767.tif 000026
式（１０）〜（１５）を式（３）に代入することでＬが求まる．
【式１６】
0004609767.tif 000027
式（１６）を式（５）に代入すると，
【式１７】
0004609767.tif 000028
より次式が得られる．
【式１８】
0004609767.tif 000029
これより次の制御則が得られる．
【式１９】
0004609767.tif 000030
同様にして，
【式２０】
0004609767.tif 000031
より次式が得られる．
【式２１】
0004609767.tif 000032
これより次の制御則が得られる．
【式２２】
0004609767.tif 000033
従って最終的な制御則は次の形になる．
【式２３】
0004609767.tif 000034
ρ１，ρ２はそれぞれの重み係数で，その和は次式とする．
【式２４】
0004609767.tif 000035
評価関数の重みを考慮したＭ１，Ｍ２へのｕ１，ｕ２の感度を等しいと考えると次のようになる．
【式２５】
0004609767.tif 000036
以上より最適制御則は次のように求まる．
【式２６】
0004609767.tif 000037
ここで，ρ１＝２Ｍ１／（Ｍ１＋ｒ１Ｍ２），ρ２＝２ｒ１Ｍ２／（Ｍ１＋ｒ１Ｍ２）である．本手法から求まる必要条件はここまでであって，κの最適化は探索的に求める必要がある．
【０００９】  ＭＡＴＬＡＢ／ＳＩＭＵＬＩＮＫによるランダム路面入力シミュレーションにより制御性能を評価する．以
0004609767.tif 000038
一定でレベルは極良路（ＩＳＯ／ＴＣ１０８／クラスＡ）を６０Ｋｍ／ｈ程度での走行に相当するように設定した．シミュレーションは，１０−４ｓｅｃの固定ステップで２ｓｅｃ間行った．
【００１０】上記シミュレーション条件において，Ｊを最小にするκ探索的に求めた結果を図２に示す．これより最適値は，κ＝−３．６×１０−７であることがわかる．
【００１１】  シミュレーション結果を図３，４に示す．制御性能比較のため準最適制御以外にパッシブダン
0004609767.tif 000039
した．パッシブダンパ制御は図１のパッシブダンパＣ２と同等の機能を制御で出すものである．スカイフックダンパ制御はばね下振動を押さえる機能がないためダンパが必要でありＣ２＝１２００Ｎｓｅｃ／ｍとし，他はすべて  Ｃ２＝０とした．評価関数Ｊを振動に関する部分Ｊ１と制御入力エネルギに関する部分Ｊ２に分離すると次のようになる．
【式２７】
0004609767.tif 000040
【式２８】
0004609767.tif 000041
図３は各制御についてＪ１，Ｊ２の比較を行ったものである．Ｊ２についてはｕｏｐｔとｕ＊ｏｐｔは同等であるが，Ｊ１ではｕｏｐｔの方が若干優れることがわかる．
図４は，本制御と準最適制御について制御対象に伝達されるパワーのシミュレーション結果を示す．ともに大半の時間帯で伝達パワーは負の値を示しており，エネルギ回生を効果的に行っていることがわかる．
以上により，本制御則の式（２６）は近似解法による準最適制御則（特許文献１の式（２０））と比べ簡潔でありながら同等以上の効果があることが確認できた．
【００１２】  次に非ホロノミック系システムの最適制御に適用した例を示す．近年，非ホロノミックな拘束をもつ非線形系のシステムへの関心が高まり各種の制御が研究されている．ここでは非特許文献３において解説されている転がるコインの制御を例題にとり，本手法により導いた最適制御則の検討を行う．
【００１３】  転がるコインのモデルを図５に示す．平面上に直行座標ｘ−ｙを取り，コインの転がっていく方向の角度をｑ１，平面上の位置を（ｑ２，ｑ３），コインの転がり角度をｑ４とする．ｑ１，ｑ４を増やす方向のトルクをｕ１，ｕ２とする．さらにコインの半径や慣性モーメントを１とし質量はゼロとする．
まず図５のモデルから式（１）の各要素を求めると次のようになる．
【式２９】
0004609767.tif 000042
【式３０】
0004609767.tif 000043
【式３１】
0004609767.tif 000044
【式３２】
0004609767.tif 000045
さらに非ホロノミックな拘束条件として次式が加わる．
【式３３】
0004609767.tif 000046
【式３４】
0004609767.tif 000047
このようにｕがＭを直接制御する系ならば，非ホロノミック系であっても本手法により最適制御則を解析的に導くことが出来ることを次に示す．
【００１４】  転がるコインの制御について具体的な定式化を行う．制御の目的は任意の初期状態からｑ１，ｑ２，ｑ３をゼロにすることであるが，モデルから明らかなように，ｑ１はｕ１を２回積分するだけで
0004609767.tif 000048
でよい．
【式３５】
0004609767.tif 000049
ｕ１は振動的にする必要があることが既に知られているため，ｋ１＝０．３，ｋ２＝１とした．この制御則ではｑ１の初期値が小さいと振動が小さすぎてｑ２，ｑ３の制御に支障をきたすが，この場合は例えば事前にｑ１をある値に設定してから全体の制御を開始するといった二段搆えの制御を行えば問題は回避できる．ｑ２，ｑ３に影響を与えずにｑ１を事前にある値に設定することは容易であるためこの方式は許容できよう．あるいはｑ１の初期値にかかわらずｑ１について式（３５）の制御則と同様な次のような位置制御をプログラムしておくといった方法でもよい．
【式３６】
0004609767.tif 000050
以上の事前検討を基にｕ２についての最適制御を考え，評価関数を次式とする．
【式３７】
0004609767.tif 000051
ここで，ｒ２，ｒ３，ｒ４，ｒ５  は重み係数である．上式右辺第１，２項はｑ２，ｑ３をゼロに収束させるという制御目的を表す項で，その意味をわかりやすくするため積分した結果を表示している．第３項がコインの回転運動に対する減衰を評価する項であり，第４項が制御対象への入力エネルギの項である．
【００１５】  続いて制御則を導出する．式（２９）〜（３４），（３７）を用いてＬを求める．ここでｇは，式（３７）右辺第１，２項を積分の形に戻した後に第１〜３項の被積分関数を加算したものに式（３３）（３４）の拘束条件を加味している．
0004609767.tif 000052
上式右辺第１行がパワー収支式から得られる項であり，第２，３行が評価関数から得られる項である．このケースではパワー収支式から得られる情報が少ないので評価関数に必要十分な情報を盛り込みこれを補っている．ここで，留意しなければならないのは，モデルから運動方程
0004609767.tif 000053
小化する条件の一部として既に式（３９）に含まれているためさらに代入することはＬに余計な拘束を強いることになる．
ｕ２を最適にするための必要条件は次の方程式から求めることができる．
【式３９】
0004609767.tif 000054
式（３８）を代入すると，
【式４０】
0004609767.tif 000055
となりこれより次の制御則が得られる．
【式４１】
0004609767.tif 000056
【００１６】  式（３５），（４１）の制御則を用いてシミュレーションによる評価を行う．
まず，重み係数ｒ２，ｒ３，ｒ４の選定を行う．図５のモデルからわかるようにｑ１を早めにゼロに収束させるとｑ３の調整が出来なくなる．従ってｑ３を収束させた後にｑ２を収束させることが望ましくこのためｒ２よりｒ３を大きく設定する．またコインの回転のオーバーシュートを抑えるためｒ４は大き目の値に設定する．以上よりｒ２＝０．２５，ｒ３＝２．５，ｒ４＝１．０とした．このケースでは入力エネルギは最終的にゼロになるため，ｒ５の値はさほど重要ではない．他の重みとレベルを合わせて，ｒ５＝０．５とした．
シミュレーション結果を図６〜９に示す．まず，κの最適化を行う．κを変化させた場合の式（３７）の値の変化を図６に示す．これよりκ＝−０．１が最適値であることがわかる．以後この値を用いてシミュレーションを行う．
図７は［ｑ１，ｑ２，ｑ３］の初期値を［１，１，−１］としてシミュレーションを行った結果である．ｑ１の振動を利用してまずｑ３が収束し続いてｑ２がゆっくり収束している様子がわかる．
図８は，このときのｑ２，ｑ３の軌跡を示す．ｑ２は一度の速度切り替えで目標に収束しており良好な結果が得られた．
図５のモデルを漸近安定させる一例として非特許文献３において紹介されている制御則を用いて，図７と同一条件でシミュレーションした結果を図９に示す．特に収束性を良くする工夫はされてなく，漸近安定させる例として示されたものなので，本提案の結果との直接性能比較の対象ではないが，参考として掲示する．
【００１７】  続いて、本制御はロバスト性能が優れていることを示す．再び図１のモデルを用いて従来の代表的な最適制御理論であるＬＱ制御と比較する．
【００１８】  ＬＱ制御則を求めるためには，図１のモデルを状態方程式で記述し直し，評価関数を２次形式にする必要がある．図１において，ｑ０をｘ０とし，ｑ１−ｑ０をｘ１とし，ｑ２−ｑ０をｘ２とし，ばね下速度をｘ３とし，ばね上速度をｘ４とすると状態方程式は以下のようになる．
【式４２】
0004609767.tif 000057
評価関数を次のような２次形式とする．
【式４３】
0004609767.tif 000058
式（４３）の被積分関数の第１項は式（２７）の性能評価項と同一であるが，第２項は異なりｕの２乗になっている．無限評価区間であるから，ＬＱ制御則は次のリカッチ代数方程式の解から求まる．
【式４４】
0004609767.tif 000059
本手法とＬＱ制御とは制御入力の評価関数が異なるため，性能比較を行うには性能評価値Ｊ１の値を基準値にとりこの値が同じになるようにして比較する必要がある．Ｊ１の値が図３に示す本手法の値と同じになるように，ｒ１ｑ＝１．５×１０−９とし，Ａ，Ｂ，Ｑに数値を入れると制御則は次のようになる．
【式４５】
0004609767.tif 000060
0004609767.tif 000061
である．これを式（２６）と同じ形式で記述し両者を数値で比較すると次のようになる．
【式４６】
0004609767.tif 000062
【式４７】
0004609767.tif 000063
【００１９】  本制御とＬＱ制御についてロバスト性能を比較する．式（４６），（４７）を比較すると，ＬＱ制御はｑ０とｑ１のゲインが大きいため信号ノイズに対するロバスト性の低下が懸念される．現実の制御では，路面入力ｑ０は推定値を用いることになるため，ある程度の推定誤差は前提とせざるをえない．ｑ０の推定誤差の影響についてのシミュレーション評価を行った．路面入力ｑ０とは無相関で統計的特性が同じ信号をｙ０とし，推定誤差相当のランダムノイズをαｙ０とした．αをパラメータとしたシミュレーション結果を図１０に示す．色塗りで示した本手法のプロットはｑ０の推定誤差による影響をほとんど受けずロバスト性が高いことがわかる．評価値Ｊは本手法が，評価値ＪＬＱではＬＱ制御が優れるのは当然であるが，αが０．２以上では評価値ＪＬＱにおいても本制御がＬＱ制御より優れることがわかる．従って従って，本手法は外乱入力の推定誤差に対するロバスト性が優れるといえる．
【００２０】  以上が今回の提案内容であるが、一般論としては式（２）が必ず最小値を持つことを証明しておく必要があるため、以下で確認する．
式（１）はパワー収支式であるが，変位入力の外乱によるパワーが陽に表されていない．これを明示し入力側を左辺に蓄積・散逸側を右辺に書き直せば，次式のようになる．
【式４８】
0004609767.tif 000064
ここで，ｆは変位入力外乱ｚに抗する制御対象側の反力，Ｕはポテンシャルエネルギである．式（１）との対応より次式の成立は明らかである．
【式４９】
0004609767.tif 000065
0004609767.tif 000066
動方程式の記述を用いて書き直し，残りのダンピング力を改めてｄ’とすれば，式（４８）の右辺は次のようになる．
【式５０】
0004609767.tif 000067
ここで，上式の第２項はコリオリ力によるもの，第３項は遠心力によるもの，第４項はダンピング力のみによるものである．式（５０）を積分すると，
【式５１】
0004609767.tif 000068
となる．ここで上式の第２項と３項は相殺され，結局式（５１）は次式になる．
【式５２】
0004609767.tif 000069
上式第１項は運動エネルギ，第２項は散逸エネルギ，第３項はポテンシャルエネルギであり，これは制御対象に蓄えられるエネルギと散逸するエネルギの総和である．式（４８）の左辺の積分は，制御対象に加えられるエネルギの総和であるから，結局式（４８）を積分したものがエネルギ保存則になっていることがわかる．
ここで式（５２）の第１，３項の和をＥ１とする．外乱入力パワーの和の積分をＥ２とすると，式（４８）において外乱入力項を右辺に移項して，左右辺を積分すれば，
【式５３】
0004609767.tif 000070
となる．ここで，対象が工学的に意味のあるシステムであれば外乱入力パワーは有界であるか
0004609767.tif 000071
後の項は散逸エネルギで非負であり時間とともに単調増加する．よって，次式が導かれる．
【式５４】
0004609767.tif 000072
また工学的に意味のある評価関数なら，有界な状態変数に対してｇは下限値を持つように選ばれているから，その有限評価区間の積分値も下限値ｃを持つ．
よって評価関数Ｊは次式を満たし最小値を持つ．
【式５５】
0004609767.tif 000073
【特許文献１】特願２００５−１７５６５６システムの最適制御方法
【非特許文献３】藤本健治，ハミルトニアンシステムの制御，計測と制御Ｖｏｌ．３９，Ｎｏ．２（２０００），ｐｐ．９９−１０４．
【発明を実施するための形態】
【００２１】  以下に、本発明を実施するための最良の形態を図面に基づいて説明する．
【００２２】  図１１にアクティブサスペンションの構造図を示す．アクチュエータは車輪１と車体２との間に設置された油圧シリンダ３とピストンロッド４とアキュームレータ５で構成されている．センサは，上下Ｇセンサ６，７とストロークセンサ８とサスペンション伝達力センサ９で構成され，これらの信号により制御回路１０により制御信号を発生させている．
双方向ポンプ１１は誘導モータ１２で回される．誘導モータ１２が右回転する場合はポンプ１１からシリンダ３の上室に作動油が送られ、左回転する場合は下室に作動油が送られる．誘導モータはインバータ１３の３相出力で駆動され，この３相出力は制御指令値によって制御される．バッテリーもしくはキャパシタ１５とインバータ１３の間には可逆チョッパ回路１４が設けられモータ駆動時にはバッテリー電圧を昇圧し，充電時にはモータからの電圧を調整しモータ電圧がバッテリー電圧を下回っても回生可能にする．
他の車輪も同様の構成である．
【００２３】  次に、制御回路１０について説明する．図１２はセンサと制御回路１０のブロック線図を示す．ばね上上下Ｇセンサ６、ばね下上下Ｇセンサ７、サスペンション伝達力センサ９、サスペンション変位センサ８からの信号が、制御回路１０に伝送される．ばね上速度
0004609767.tif 000074
【式５６】
0004609767.tif 000075
【００２４】  制御則演算部２０は、モータ制御部と回生制御部からなる．モータ制御部は上記ばね上速度
0004609767.tif 000076
じてモータの電流を変えれば、制御指令値に応じて力を発生させる制御系が実現できる．回生
0004609767.tif 000077
ｕｏｐｔとからモータ駆動時か回生時かを判断し、モータ駆動時にはバッテリー電圧を昇圧し，回生時にはモータ電圧を昇圧する指令値ｕｅを出力する．
【産業上の利用可能性】
【００２５】  本手法は、非ホロノミック拘束を含む非線形制御系の適用が可能であり、しかもシステムのエネルギ収支式から、微分方程式を解かずに最適制御則を得ることができるため有用性は高いものと考えられる．
【図面の簡単な説明】
【００２６】
【図１】  車両のアクティブサスペンションの１輪モデル図である．
【図２】  パラメータκの最適化の説明図である．
【図３】  アクティブサスペンションについて本発明の制御と従来の制御方法の性能比較図である．
【図４】  本発明と特許文献１との制御入力パワーの比較図である．
【図５】  転がるコインのモデル図である．
【図６】  パラメータκの最適化の説明図である．
【図７】  本発明の制御効果を説明する図である．
【図８】  本発明の制御効果を説明する図である．
【図９】  従来技術の制御効果を説明する図である．
【図１０】  本発明と従来制御のロバスト性の比較図である。
【図１１】  本発明のアクテイブサスペンションの構成図である．
【図１２】  本発明のアクティブサスペンションのセンサと制御回路の構成図である．
【符号の説明】
【００２７】１  車輪（ばね下質量）
２  車体（ばね上質量）
３  油圧シリンダ
４  ピストンロッド
５  アキュームレータ
６  上下Ｇセンサ（ばね上用）
７  上下Ｇセンサ（ばね下用）
８  サスペンション変位センサ
９  サスペンション伝達力センサ
１０  制御回路
１１  双方向ポンプ
１２  誘導モータ
１３  インバータ
１４  可変チョッパ回路
１５  バッテリもしくはキャパシタ
２０  制御則演算部
２１  ばね上速度演算部
２２  ばね下速度演算部
２３  タイヤ変位演算部
【図１】
0004609767.tif 000083
【図２】
0004609767.tif 000084
【図３】
0004609767.tif 000085
【図４】
0004609767.tif 000086
【図５】
0004609767.tif 000087
【図６】
0004609767.tif 000088
【図７】
0004609767.tif 000089
【図８】
0004609767.tif 000090
【図９】
0004609767.tif 000091
【図１０】
0004609767.tif 000092
【図１１】
0004609767.tif 000093
【図１２】
0004609767.tif 000094