(57)【要約】【課題】既存の制御システムの制御特性は変更せずに実機において操縦者の所望の操縦特性を容易に実現する操縦システムを提供する。【解決手段】制御対象である実機システム１１にフィードバック制御器１２を組み込んだ拡大実機システム１０の入力として、拡大実機システム１０のＨ∞ノルム定義型右逆システムであるフィードフォワード制御器２０を結合し、そのフィードフォワード制御器２０の入力に対し操縦者の所望の操縦特性を模擬したオーダモデル３０を結合したフライト操縦システム１００を構築する。なお、フィードバック制御器１２およびフィードフォワード制御器２０はデジタル変換され計算機内部に先に実装され、オーダモデル３０はデジタル変換され計算機内部または外部に後から実装される。【選択図】図１
2006327299.tif 000002
【特許請求の範囲】
【請求項１】
実機システムに対し該実機システムの出力に対応する入力を出力する右逆システムがフィードフォワード制御器として結合された操縦システムであって、前記右逆システムに対し所望の操縦特性が模擬されたオーダモデルが結合されて前記実機システムは前記オーダモデルと同様な操縦特性を有することを特徴とする操縦システム。
【請求項２】
前記右逆システムは、Ｇ(ｓ）を前記実機システムの入出力特性を示す伝達関数とし且つ、Ｗ(s)を低周波数帯域において大きなゲインを有する安定かつプロパーな対角重み伝達関数とする時、‖Ｗ(s)（Ｇ(s)Ｇ＊(s)−Ｉ）‖∞＜１（Ｉ：単位行列）を満たす安定かつプロパーなＨ∞ノルム定義型右逆システムＧ＊(s)である請求項１に記載の操縦システム。
【請求項３】
前記実機システムは、不安定である場合、フィードバック制御システムが結合され安定化される請求項１に記載の操縦システム。
【請求項４】
前記右逆システムおよび前記フィードバック制御システムは計算機内部に実装され、対する前記オーダモデルは計算機内部または外部に実装され応答特性を任意に変更することができ且つ交換可能である請求項１から３の何れかに記載の操縦システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、操縦システム、特に既存の制御システムの制御特性は変更せずに対象とするシステムにおいて各操縦者が希望する任意の操縦特性を実現することができる操縦システムに関するものである。
【背景技術】
【０００２】
あるシステムの入出力特性を示す伝達行列関数をＧ(s)と表記すると、理想的な右逆システムは、Ｇ(s)Ｇ＊(s)＝Ｉ（ここで、Ｉは単位行列）を満たすＧ＊(s)である。従来から研究されている逆システムは、Silvermanらの理想的な逆システム、Sainらの積分値を再現する逆システム、吉川らの低周波数成分のみを再現する逆システム（例えば、非特許文献１を参照。）の３つに分けられる。これらのうち、吉川らの「低域通過逆システム」は微分器を必要としないため、理想的な逆システムと考えられ、伝達関数上での設計法（例えば、非特許文献１を参照。）および状態空間上での設計法（例えば、非特許文献２を参照。）が研究されている。また、従来では設計することができなかった、開右半平面に不変零点を有するシステムを対象とする低域通過疑似逆システム（例えば、非特許文献３を参照。）についても研究されているが、位相を回復させることは出来ないという問題ならびに虚軸上に不変零点を有するシステムについては適用することは出来ないという問題があった。
ところで、システムには必ず遅れが存在し、その遅れを記述する方法として、Pade近似モデルは広く一般に使用されている。しかし、Pade近似モデルの導入は必ず不安定な零点を導入するという問題点がある。なお、不安定な零点は制御性を低下させることが一般に知られている。
他方、人間が操縦するシステムの場合、個々の人間に応じて所望の操縦特性が異なる。例えば、自動車の操縦システムを例にとると、ドライバーのステアリング操作力を軽減する油圧システムにおいては、女性のドライバーの場合は、油圧が大きい方が操作しやすい反面、男性のドライバーの場合は、油圧が小さい方が操作しやすいということも起こり得る。しかし、実際は時間およびコストの観点から、油圧システム等の操縦に係るシステムの設計が個々の人間に応じて行われることはなく、平均化した人間を対象とし最大公約数的に行われている。このように、人間が操縦するシステムは、コストと時間の制約から、個々のユーザの操縦特性が反映されたものとなっていないのが現状である。
一方で、人間が操縦する以上、その人間に適した所望の操縦特性が存在するため、その人間に適した制御系の設計を容易に行うことができれば、より快適な操縦が得られ、満足度が高まると考えられる。
【０００３】
【非特許文献１】吉川、杉江：高周波遮断システムとそのサーボ系への一応用、計測自動制御学会論文集、１８−８、７９２／７９９（１９８２）
【非特許文献２】山田、渡部：低域通過逆システムの状態空間法での構成、計測自動制御学会論文集、２８−８、９２３／９３０（１９９２）
【非特許文献３】山田、渡部：安定な低域通過逆システムの状態空間法による構成法法、計測自動制御学会論文集、３２−６、８６２／８７０（１９９６）
【発明の開示】
【発明が解決しようとする課題】
【０００４】
フィードバック制御器による制御性能（フィードバック制御器が存在しない場合は、システム固有の特性）を変えることなく、しかも、制御系を再設計することなく、後付けする任意のモデルの操縦特性を対象とするシステムにおいて実現するためには、対象とするシステムの右逆システムをフィードフォワード制御器として用いることが必要かつ十分な条件である。
しかし、上述した通り，不安定な不変零点を有するシステムに対して有効な右逆システムを設計する手法がなかったので，フィードバック制御器による制御性能，もしくはシステム固有の特性を変えることなく、対象とするシステムにおいて所望の操縦特性を実現する制御系は存在しなかった。
そこで、本発明は上記従来技術の問題点に鑑みなされたものであって、その解決しようとする課題は、既存の制御システムの制御特性は変更せずに対象とするシステムにおいて各操縦者が希望する任意の操縦特性を実現することができる操縦システムを提供することにある。
【課題を解決するための手段】
【０００５】
上記目的を達成するための第１の発明は、実機システムに対し該実機システムの出力に対応する入力を出力する右逆システムがフィードフォワード制御器として結合された操縦システムであって、前記右逆システムに対し所望の操縦特性が模擬されたオーダモデルが結合されて前記実機システムは前記オーダモデルと同様な操縦特性を有することを特徴とする。
上記第１の発明の操縦システムでは、実機システムの入力に対し実機システムの右逆システムであるフィードフォワード制御器が結合されているので、フィードフォワード制御器の入力として結合されているオーダモデルの出力が実機システムの出力としてほぼ同様に再現されることとなる。従って、オーダモデルの応答特性は操縦者の希望する操縦特性を模擬しているので、実機システムにおいて操縦者が希望する操縦特性が実現される。
【０００６】
第２の発明は、前記右逆システムは、Ｇ(ｓ）を前記実機システムの入出力特性を示す伝達関数とし且つ、Ｗ(s)を低周波数帯域において大きなゲインを有する安定かつプロパーな対角重み伝達関数とする時、‖Ｗ(s)（Ｇ(s)Ｇ＊(s)−Ｉ）‖∞＜１（Ｉ：単位行列）を満たす安定かつプロパーなＨ∞ノルム定義型右逆システムＧ＊(s)であることとした。
上記第２の発明の操縦システムでは、対角重み伝達関数Ｗ(s)は低周波数帯域においてはゲインが十分に大きいので、低周波帯域においてはＧ(s)Ｇ＊(s)≒Ｉが成り立つことがわかる。従って、Ｇ＊(s)は、低周波数帯域においてゲインが十分に大きい対角重み伝達関数Ｗ(s)対するＧ(s)の右逆システムとなり、そのＧ＊(s)を実機システムに対するフィードフォワード制御器として用い且つその入力に対し操縦者の所望の操縦特性を模擬したオーダモデルを結合することにより、実機システムにおいてオーダモデルの出力と同じ出力を得ることができ、実機システムにおいて所望の操縦特性を実現することができる。
【０００７】
第３の発明では、前記実機システムは、不安定である場合、フィードバック制御システムが結合され安定化されることとした。
上記第３の発明の操縦システムでは、制御対象である実機システムが安定になる。
【０００８】
第４の発明では、前記右逆システムおよび前記フィードバック制御システムは計算機内部に実装され、対する前記オーダモデルは計算機内部または外部に実装され応答特性を任意に変更することができ且つ交換可能であることとした。
上記第４の発明の操縦システムでは、右逆システムおよびフィードバック制御システムは計算機内部に先に実装され、対するオーダモデルは後から実装することによって、既存の制御システムの制御特性は変更せずに実機システムにおいて操縦者の所望の操縦特性を実現することが出来るようになる。
【発明の効果】
【０００９】
本発明の操縦システムによれば、実機システムの入力となるフィードフォワード制御器にＨ∞ノルム定義型右逆システムを用いることで、実機システムが不安定な不変零点を有する場合であっても実機システムの右逆システムが実現され、その結果、その右逆システムを実機システムに対するフィードフォワード制御器として用いることにより、実機システムの位相遅れの回復が好適に達成され、更にそのフィードフォワード制御器の入力に対し操縦者の希望する応答特性を模擬したオーダモデルを結合することにより、既存の制御システムを新たに再設計することなく、実機システムにおいて操縦者の希望する任意の操縦特性を実現することができる。また、実機システムが不安定の場合は、フィードバック制御システムを結合することにより、実機システムが安定化される。また、右逆システムおよびフィードバック制御システムは計算機内部に先に実装され、対するオーダモデルは後から実装することによって、既存の制御システムの制御特性は変更せずに実機システムにおいて操縦者の所望の操縦特性を容易に実現することが出来るようになる。
【発明を実施するための最良の形態】
【００１０】
以下、図に示す実施の形態により本発明をさらに詳細に説明する。
【００１１】
図１は、本発明の操縦システムを航空機に適用したフライト操縦システム１００を示す構成説明図である。
このフライト操縦システム１００は、対象プラントである実機システム１１とその実機システム１１の応答を安定化するために設計されたフィードバック制御器１２とから成る拡大実機システム１０と、拡大実機システム１０のＨ∞ノルム定義型右逆システムであるフィードフォワード制御器２０と、所望の操縦特性が実現されているオーダモデル３０とを具備して構成されている。なお、一点鎖線内部のフィードバック制御器１２およびフィードフォワード制御器２０は、例えばオンボードコンピュータ等の計算機内部に後述するデジタル変換によって先に実装され、対するオーダモデル３０は計算機内部または外部に同デジタル変換によって後から実装される。
【００１２】
実機システム１１は、航空機のシステムであり、本実施形態ではエルロン舵角δａおよびラダー舵角δｒを入力とし横方向速度ｖおよびロール角φを出力とする二入力二出力のシステムである。また、実機システム１１は不安定な不変零点を有するシステムであっても良く、不安定な場合は、上述したフィードバック制御器１２が結合され安定化される。例えば、図２にて示されるように、実機システム１１は、横／方向運動モデル１とアクチュエータモデル２と遅れモデル３が直列に結合したシステムとして構成される。
【００１３】
横／方向運動モデル１は、例えば入力変数がエルロン舵角δａ[rad]およびラダー舵角δｒ[rad]であり、状態変数が横方向速度ｖ[m/s]、ロールレートｐ[rad/s]、ロール角φ[rad]およびヨーレートｒ[rad/s]であり、制御したい出力（または出力変数）が横方向速度ｖ[m/s]およびロール角φ[deg]である二入力二出力のシステムである。なお、フィードバック制御器１２に帰還する変数は全ての状態変数とした。また、横／方向運動モデル１の状態空間表現は、以下の通り記述される。
【００１４】
【数１】
2006327299.tif 000003
【００１５】
アクチュエータモデル２は、アクチュエータとして電動アクチュエータを用いた。また、エルロンおよびラダーシステムは、ともに１次システムを用いて、例えば以下の通り記述される。
【００１６】
【数２】
2006327299.tif 000004
なお、δaｃはエルロンの舵角コマンドであり、δrｃラダーの舵角コマンドを表す。
【００１７】
遅れモデル３は、エルロンおよびラダーシステムとともに以下に示すPade一次近似モデルを用いた。
【００１８】
【数３】
2006327299.tif 000005
なお、遅れモデル３に対しPade一次近似モデルを適用したため、実機システム１１は不安定な不変零点を有している。
【００１９】
フィードバック制御器１２は、実機システム１１が不安定である場合、実機システム１１を安定化させる制御器である。例えば、横方向からの突風等の外乱に対する横方向速度ｖおよびロール角φの影響を最小化する。なお、フィードバック制御器１２の状態空間表現（Ｋ）は、例えば以下の通り記述される。
【００２０】
【数４】
2006327299.tif 000006
【００２１】
また、フィードバック制御器１２を加えた拡大実機システム１０の状態空間表現は、以下の通り記述される。
【００２２】
【数５】
2006327299.tif 000007
なお、ＢＰにおいて、０４×２は４行２列のゼロ行列を示し、同様にＣＰにおいて、０２×４は２行４列のゼロ行列を示している。また、拡大実機システム１０は、不安定な不変零点を有している。後述する通り、拡大実機システム１０が不安定な不変零点を有する場合であっても、拡大実機システム１０に対応する右逆システムを設計することができる。
【００２３】
フィードフォワード制御器２０は、拡大実機システム１０の入出力特性を示す伝達関数をＧ(s)とし且つ、Ｗ(s)を低周波数帯域において大きなゲインを有する安定かつプロパーな対角重み伝達関数とする時、‖Ｗ(s)（Ｇ(s)Ｇ＊(s)−Ｉ）‖∞＜１（Ｉ：単位行列）を満たす安定かつプロパーなＨ∞ノルム定義型右逆システムＧ＊(s)に相当し、例えば以下の手順により設計することができる。
［設定手順１］
【００２４】
先ず、対角重み関数Ｗ(s)を定める。すなわち、上記右逆システムＧ＊(s)が必要な周波数帯域において、ゲインが十分に大きく且つ上記右逆システムＧ＊(s)が不必要な周波数帯域において、ゲインが十分小さくなるように対角重み関数Ｗ(s)を定める。なお、本実施形態では１[rad/s]以下の低周波数帯域においてゲインが十分に大きく且つ５[rad/s]以下においてゲインが１（＝０[dB]）以上とする対角重み関数Ｗ(s)を下記の通り定めた。
【００２５】
【数６】
2006327299.tif 000008
【００２６】
上記対角重み関数Ｗ(s)のゲイン特性は、図３のゲイン線図にて示されるように、１[Hz]以下の低周波数帯域においてはゲインが十分大きく、５[Hz]以下の低周波数帯域ではゲインが減衰していない。この対角重み関数Ｗ(s)の状態空間表現｛Ａw,Ｂw,Ｃw,Ｄw｝と拡大実機システム１０の状態空間表現を用いて、一般化プラントを次のように定める。
【００２７】
【数７】
2006327299.tif 000009
【００２８】
この一般化プラントに対するＨ∞問題を、例えば市販の「Matlab LMI Control Toolbox」を用いて解き、Ｈ∞ノルム定義型右逆システムＧ＊(s)を求める。得られた右逆システムＧ＊(s)をフィードフォワード制御器２０として拡大実機システム１０に結合したシステムのボード線図を図４に示す。なお、左上のゲインおよび位相のボード線図が、二入力（フィードフォワード制御器２０の入力）に対する二出力（横方向速度ｖ[m/s]およびロール角φ[deg]）の伝達関数行列の（１，１）要素の特性を示し、左下の同線図が（２，１）要素の特性を示し、右上の同線図が（１，２）要素の特性を示し、右下の同線図が（２，２）要素の特性を示している。また、参考として、機体の横／方向運動モデル１、実機システム１１、拡大実機システム１０のボード線図も併せて示した。これらにより、実機システム１１は機体の横／方向運動モデル１に比べて遅れモデル３が加えられた分の位相が遅れていることがわかる。
【００２９】
また、２[Hz]以下の低周波数帯域において、拡大実機システム１０にＨ∞ノルム定義型右逆システムであるフィードフォワード制御器２０を付け加えることによって以下の３点が達成されていることがわかる。
（１）対角項のゲインがほぼ１であること。
（２）拡大実機システム１０の対角項の位相遅れ、または位相進みをフィードフォワード制御器２０が回復もしくは遅れさせ、対角項の位相がほぼ零となっていること。
（３）非対角項のゲインが十分小さいこと。
【００３０】
［設定手順２］
搭載コンピュータがデジタルコンピュータであるので、得られたフィードフォワード制御器２０を適当な方法、例えば双一次変換等のデジタル変換を用いて例えば０.０２[sec]（５０[Hz]）で離散化し実装する。なお、搭載はフィードフォワード制御器２０およびフィードバック制御器１２に相当する部分を、実機システム１１を制御するコンピュータに搭載する。そして、その後に、オーダモデル３０に相当するモデル（デジタルコンピュータを用いるのであれば、離散化したモデル）を加える。
【実施例１】
【００３１】
なお、オーダモデル３０の実施例１に係るモデルとしては例えば、モデルへの操舵入力がモデルのエルロン舵角δａ[deg]およびラダー舵角δｒ[deg]とした時、モデルのエルロン舵角δａに対するオーダモデル３０のロール角φ[deg]の伝達関数が１／(ｓ２＋２ｓ＋１）であり且つモデルのラダー舵角δｒに対するオーダモデル３０の横方向速度ｖ[m/s]の伝達関数が１／(ｓ２＋２ｓ＋１）であるシステムモデルとした。
【００３２】
図５は、本発明のフライト操縦システム１００に実施例１に係るオーダモデル３０を適用したhardware-in-the-loopシミュレーションの結果を示す説明図である。なお、hardware-in-the-loopシミュレーションとは、実験用航空機「ＭｕＰＡＬ−α」のアクチュエータモデル２及び遅れモデル３に相当するハードウェア機器と、横／方向運動モデル１に相当する同「ＭｕＰＡＬ−α」の運動計算を行う計算機を組み合わせたシミュレーションのことである。また、破線はオーダモデル３０の時歴であり、他方、実線は、オーダモデル３０の出力を得てフィードフォワード制御器２０が拡大実機システム１０の入力を計算し、その入力を拡大実機システム１０に入力することにより得られた拡大実機システム１０の出力の時歴である。
この図から、実際の出力である拡大実機システム１０のロール角φおよび横方向速度ｖの時歴がオーダモデル３０の出力であるロール角および横方向速度の時歴にほぼ重なっていることがわかる。つまり、拡大実機システム１０の応答特性がオーダモデル３０の応答特性にほぼ等しくなり、実機において操縦者の所望の操縦特性が実現されている。
【実施例２】
【００３３】
オーダモデル３０の実施例２に係るモデルとしては例えば、モデルのエルロン舵角δａに対するオーダモデル３０のロール角φの伝達関数が０.６２／(ｓ２＋０.４８ｓ＋０.６２）であり且つモデルのラダー舵角δｒに対するオーダモデル３０の横方向速度ｖの伝達関数が０.８２／(ｓ２＋０.８ｓ＋０.８２）であるシステムモデルとした。
【００３４】
図６は、本発明のフライト操縦システム１００に実施例２に係るオーダモデル３０を適用したhardware-in-the-loopシミュレーションの結果を示す説明図である。なお、破線がオーダモデル３０の時歴であり、実線が拡大実機システム１０の時歴である。
この図から、実際の出力である拡大実機システム１０のロール角φおよび横方向速度ｖの時歴がオーダモデル３０の出力である横方向速度およびロール角の時歴にほぼ重なっていることがわかる。つまり、拡大実機システム１０の応答特性がオーダモデル３０の応答特性にほぼ等しくなり、実機において操縦者の所望の操縦特性が実現されている。
【実施例３】
【００３５】
オーダモデル３０の実施例３に係るモデルとしては例えば、モデルのエルロン舵角δａに対するオーダモデル３０のロール角φの伝達関数が１２／(ｓ２＋０.６ｓ＋１２）であり且つモデルのラダー舵角δｒに対するオーダモデル３０の横方向速度ｖの伝達関数が０.５２／(ｓ２＋１.０ｓ＋０.５２）であるシステムモデルとした。
【００３６】
図７は、本発明のフライト操縦システム１００に実施例３に係るオーダモデル３０を適用したhardware-in-the-loopシミュレーションの結果を示す説明図である。なお、破線がオーダモデル３０の時歴であり、実線が拡大実機システム１０の時歴である。
この図から、実際の出力である拡大実機システム１０のロール角φおよび横方向速度ｖの時歴がオーダモデル３０の出力である横方向速度およびロール角の時歴にほぼ重なっていることがわかる。つまり、拡大実機システム１０の応答特性がオーダモデル３０の応答特性にほぼ等しくなり、実機において操縦者の所望の操縦特性が実現されている。
【００３７】
また、上記実施例１から実施例３に係るオーダモデルを用いた実フライトの時歴をそれぞれ図８、図９および図１０に示す。破線が、モデルの時歴、実線が拡大実機システム１０の時歴であるが、hardware-in-the-loopシミュレーションと同様に、オーダモデル３０の出力と拡大実機システム１０の「制御したい出力」が区別できないほど重なっている。
【００３８】
上記実施例１から実施例３において、フィードバック制御器１２およびフィードフォワード制御器２０は同一であり、オーダモデル３０の追加に伴う再設計を行っていないが、拡大実機システム１０の出力がオーダモデル３０の出力とほぼ同じであることが図５から図１０によって確認できる。すなわち、既存の制御器を新たに再設計することなく、オーダモデル３０を入れ替えることによって拡大実機システム１０において操縦者の所望の操縦特性が実現できることがわかる。
【００３９】
以上、本発明のフライト操縦システム１００によれば、実機システム１１を含む拡大実機システム１０の入力となるフィードフォワード制御器２０にＨ∞ノルム定義型右逆システムを用いることで、拡大実機システム１０が不安定な不変零点を有する場合であっても拡大実機システム１０の右逆システムが実現され、その結果、拡大実機システム１０の位相遅れの回復が達成され、更にそのフィードフォワード制御器２０の入力に対し操縦者の希望する応答特性を模擬したオーダモデル３０を結合することにより、既存の制御システムを新たに再設計することなく、拡大実機システム１０において操縦者の希望する任意の操縦特性を実現することができる。また、実機システム１１が不安定の場合は、フィードバック制御器１２を結合することにより、実機システム１１が安定化される。また、フィードフォワード制御器２０およびフィードバック制御器１２は計算機内部に先に実装され、対するオーダモデル３０は後から実装することによって、既存の制御システムの制御特性は変更せずに実機システム１１において所望の操縦特性を容易に実現することが出来るようになる。
【産業上の利用可能性】
【００４０】
本発明の操縦システムは、航空機に限らず人間の操縦を必要とするシステム、例えば自動車等に対しても好適に適用される。
【図面の簡単な説明】
【００４１】
【図１】本発明の操縦システムを航空機に適用したフライト操縦システムを示す構成説明図である。
【図２】図１の実機システムを示す構成説明図である。
【図３】対角重み伝達関数Ｗ(s)のゲイン特性を示すゲイン線図である。
【図４】フィードフォワード制御器２０と拡大実機システム１０を結合したシステムのゲイン特性および位相特性を示すボード線図である。
【図５】本発明のフライト操縦システムに実施例１に係るオーダモデルを適用したhardware-in-the-loopシミュレーションの結果を示す説明図である。
【図６】本発明のフライト操縦システムに実施例２に係るオーダモデルを適用したhardware-in-the-loopシミュレーションの結果を示す説明図である。
【図７】本発明のフライト操縦システムに実施例３に係るオーダモデルを適用したhardware-in-the-loopシミュレーションの結果を示す説明図である。
【図８】本発明のフライト操縦システムに実施例１に係るオーダモデルを適用した実フライト試験の結果を示す説明図である。
【図９】本発明のフライト操縦システムに実施例２に係るオーダモデルを適用した実フライト試験の結果を示す説明図である。
【図１０】本発明のフライト操縦システムに実施例３に係るオーダモデルを適用した実フライト試験の結果を示す説明図である。
【符号の説明】
【００４２】
１      横／方向運動モデル
２      アクチュエータモデル
３      遅れモデル
１０    拡大実機システム
１１    実機システム
１２    フィードバック制御器
２０    フィードフォワード制御器
３０    オーダモデル
１００  フライト操縦システム
【図１】
2006327299.tif 000010
【図２】
2006327299.tif 000011
【図３】
2006327299.tif 000012
【図４】
2006327299.tif 000013
【図５】
2006327299.tif 000014
【図６】
2006327299.tif 000015
【図７】
2006327299.tif 000016
【図８】
2006327299.tif 000017
【図９】
2006327299.tif 000018
【図１０】
2006327299.tif 000019