(57)【特許請求の範囲】
【請求項１】
航空機の飛行中における航空機の飛行パラメータのリアルタイム決定を改善するための方法であって、
以下のステップ、即ち、
−  前記航空機の推定されるべき少なくとも２つの飛行パラメータが選択されるステップと、
−  選択された飛行パラメータにそれぞれ関係しているとともに、選択された前記飛行パラメータ間に依存性関係が存在する飛行力学方程式が特定されるステップと、
−  選択された前記飛行パラメータと関連した前記特定された飛行力学方程式に基づき、少なくとも選択された前記飛行パラメータを含んだ入力パラメータの値を受け取る拡張カルマンフィルタ(１０)が構成されるステップと、
−  前記航空機の飛行中に、前記拡張カルマンフィルタ(１０)が実行されて、該フィルタが、出力として、選択された前記飛行パラメータの合同推定値を供するステップと
が実行されることを特徴とする方法。
【請求項２】
−  前記カルマンフィルタは以下の行列、即ち、
・  測定ノイズに関係すると共に測定ノイズＶの対角行列と関連した共分散行列Ｒと、
・  展開ノイズに関係すると共に展開ノイズＷの対角行列と関連した共分散行列Ｑと
によって定義され、かつ
−  以下の付加的ステップ、即ち
・  前記入力パラメータの値が許容し得るものである旨が検証されるステップと、
・  入力パラメータの値に欠陥が検出された場合に、前記共分散行列ＲおよびＱの少なくともいずれか一方の少なくともいずれか１つの要素の現在値がリアルタイムで適合化されるステップと
が実行されることを特徴とする、請求項１に記載の方法。
【請求項３】
−  予備的ステップにおいて、それぞれ測定ノイズと展開ノイズとに関係する共分散行列ＲおよびＱの複数の設定値が定義され、こうして定義された前記設定値は前記入力パラメータのうちのいずれか１個のパラメータの欠陥がある値に関連し、かつ、
−  入力パラメータの値に欠陥が検出された場合、測定ノイズと展開ノイズとに関係する前記共分散行列ＱおよびＲの現在値を適合化すべく、検出された欠陥がある値に対応する所定の設定値が測定ノイズと展開ノイズとに関係した共分散行列ＲおよびＱに割り当てられる
ことを特徴とする、請求項２に記載の方法。
【請求項４】
前記航空機に搭載の１つまたは複数のセンサによって測定された前記入力パラメータのうちのいずれか１個のパラメータの値に欠陥が検出された場合、前記欠陥がある測定値は前記拡張カルマンフィルタ(１０)の出力として供された対応する推定値によって置換されることを特徴とする、請求項２または３に記載の方法。
【請求項５】
−  前記拡張カルマンフィルタ(１０)がそのために推定値を供する選択された前記飛行パラメータのうちの少なくともいずれか１個のパラメータが考慮されるステップと、
−  前記拡張カルマンフィルタ(１０)の入力パラメータの値のうち、前記航空機に搭載のセンサに由来する、前記当該飛行パラメータに対応するものが選択されるステップと、
−  選択された前記値のうち少なくともいずれか１つと結びついた不整合が検出されるステップと、
−  選択された前記飛行パラメータの現在値が、残りの選択された１つまたは複数の値と選択された前記飛行パラメータの推定値とに基づき決定され、他方、検出された不整合な１つまたは複数の値は除外されるステップと
が実行されることを特徴とする、請求項１〜４のいずれか１項に記載の方法。
【請求項６】
前記拡張カルマンフィルタ(１０)と関連した状態ベクトルは以下の１２個の状態、即ち、
−  迎角αと、
−  地球座標系で速度νと、
−  対地速度Ｖと、
−  縦揺れ割合ｑと、
−  姿勢θと、
−  高度ｈと、
−  有効原動推力ＴＢと、
−  縦加速度ｎｘに対応する方角に投影されたバイアスｂｎｘと、
−  横加速度ｎｙに対応する方角に投影されたバイアスｂｎｙと、
−  地球座標系(ｘ、ｙ、ｚ)のｘ軸に沿った風速Ｗｘと、
−  地球座標系(ｘ、ｙ、ｚ)のｙ軸に沿った風速Ｗｙと、
−  地球座標系(ｘ、ｙ、ｚ)のｚ軸に沿った風速Ｗｚと、
によって定義されることを特徴とする、請求項１〜５のいずれか１項に記載の方法。
【請求項７】
前記拡張カルマンフィルタ(１０)の入力パラメータが、慣性パラメータ(ｑ、ｐ、ｒ
、φ、θ、ｈ、Ｖｚ、Ｖ、ｎｘ、ｎｙ、Ｎｚ)と、風速のパラメータ(α、β、αａ、βａ、Ｖａ)と、前記航空機に固有のパラメータ(Ｍ、Ｉyy)と、搭載モデリングから発生する中間パラメータ(ＦＺａ、Ｍａ、ＴＢＳ、ＭＴＢ)を含むものである請求項１〜６のいずれか１項に記載の航空機の飛行パラメータの決定方法。
【請求項８】
航空機の飛行中における航空機の飛行パラメータをリアルタイムで決定するための、入力パラメータの値を受け取ることのできる拡張カルマンフィルタ(１０)を含んだシステムであって、
−  前記拡張カルマンフィルタ(１０)は、前記入力パラメータに属する前記航空機の推定されるべき少なくとも２個の事前選択された飛行パラメータ間の相互依存関係を定める飛行力学方程式に基づき構成され、かつ、
−  前記拡張カルマンフィルタ(１０)は、前記航空機の飛行中に、選択された前記飛行パラメータの合同推定値を供するように形成されている
ことを特徴とするシステム。
【請求項９】
−  前記カルマンフィルタが以下の行列、即ち、
・ 測定ノイズと関係するとともに、測定ノイズＶの対角行列に関連した共分散行列Ｒ
と、
・ 展開ノイズと関係するとともに、展開ノイズＷの対角行列に関連した共分散行列Ｑ
と、によって定義され、
−  かつ、そのシステム(３)は以下の手段、即ち、
・  前記拡張カルマンフィルタ(１０)の前記入力パラメータ値が許容しうるものであるか否かを検証するための手段(１４)と、
・  前記検証手段(１４)によって入力パラメータの値に欠陥が検出された場合、共分散行列ＲおよびＱの少なくともいずれか一方の現在値をリアルタイムで適合化するための手段(１５)と
を含む
ことを特徴とする、請求項８に記載のシステム。
【請求項１０】
さらに、欠陥がある１つまたは複数の値を、前記拡張カルマンフィルタ(１０)によって推定されたそれらの値が得られる場合に、該推定値によって置換するための手段(１５)を含むことを特徴とする、請求項９に記載のシステム。
【請求項１１】
請求項８〜１０のいずれか１項に記載された飛行パラメータの決定システム(３)を少なくとも１つの具備する航空機。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は飛行中における航空機の飛行パラメータを決定するための方法およびシステムに関する。
【背景技術】
【０００２】
電気飛行制御の導入ならびに最新の輸送機の自動化レベルの一層の進展は輸送機の安全性を著しく改善することを可能にした。こうした飛躍的な技術の進歩によって、特に、異常な状況時、飛行システムの障害発生時、危険な周囲環境条件等に際する航空機の安全確保に積極的に関与する、特に自動操縦装置によって実施される、より安全な航法を定めることが可能になった。
【０００３】
但し、多数の航法の使用は、航空機の状態パラメータの測定に対する飛行システムの依存度を高めることになった。とりわけ、今や、航空機の飛行の安全確保にとって特定の飛行パラメータが不可欠であることが判明していることから、これらのパラメータについて信頼できる値を得ることが必要となっている。これらの不可欠なパラメータとは、特に、風速および慣性パラメータであると共に、当該航空機の速度、高度および慣性データ(特にトリム)に関係する情報を供するＡＤＩＲＳ(「Air Data and Inertial Reference System」を意味する頭字語)によって実施されるそれらのパラメータの処理を含む。
【従来技術の説明】
【０００４】
現在、飛行パラメータの値は適合化されたセンサに由来する測定値に基づき定められていることから、これらの値の妥当性は、基本的にこれらのセンサによって実施された測定の妥当性に依存している。しかしながら、外部センサにあっては、外乱(例えば、センサの着氷あるいは汚損)によって、これらの外部センサによって行われる測定の精度ならびに正確性が損なわれることがある。
【０００５】
それゆえ、航空安全当局によって課される規制上ならびに安全上の義務を満たすため、飛行パラメータの値の信頼度、精度および可用性を確かなものとすることが必要である。
【０００６】
したがって、現在行われている公知の解決方法の１つはハードウェア重複に基づく。これは、所定の飛行パラメータの値を得ることを可能にするセンサおよびコンピュータの数の増強を伴うことになる。
【０００７】
但し、所定の飛行パラメータの不整合なおよび／または誤りのある値の検出は、ある種の極端な場合にあっては、センサによって得られた値がまったく顧慮されない結果になることがあり、そのため、こうした場合には上記の当該飛行パラメータと関連した値はまったく得られないことになる。こうした値の欠如は搭載飛行システムによって使用される航法の変化を招くことがある。
【０００８】
こうしたことから、１つまたは複数のパラメータの値が完全に欠如しているような極端なケースを回避すべく、カルマンフィルタを用いて不可欠な飛行パラメータのうちの少なくともいくつかを推定することが知られている。例えば、迎角は飛行中の迎角の発展を定義する飛行力学方程式に基づき構成されたカルマンフィルタによって推定可能であることは公知に属する。こうした推定を実施するため、カルマンフィルタは入力として、搭載センサによって得られた、対気速度の測定値(迎角は後者と密接に関係している)を受け取り、出力として、迎角の推定値を供する。
【０００９】
但し、対気速度の測定値に(例えば、センサの障害に起因して)誤りがあれば、迎角の推定値も欠陥があることになる。したがって、カルマンフィルタリングによって行われる固有の飛行パラメータの推定は、当該飛行パラメータの推定に必要な入力パラメータの測定または測定値の欠陥の回避を可能とするものではない。というのも、そうした場合には、得られた推定値も誤りがあり、したがって使用不能だからである。
【００１０】
本発明の目的は、これらの欠陥を解消し、とりわけ、カルマンフィルタの１つまたは複数の入力パラメータの測定欠陥を回避して、カルマンフィルタの出力による当該飛行パラメータの信頼できる値の可用性を保証することである。
【発明の概要】
【００１１】
上記目的のため、本発明によれば、航空機の飛行中における航空機の飛行パラメータのリアルタイム決定を改善するための方法は、以下のステップ、即ち、
−  上記航空機の推定されるべき少なくとも２個の飛行パラメータが選択されるステップと、
−  選択された上記飛行パラメータとそれぞれ関連していると共に、選択された上記飛行パラメータ間の相互依存関係が存在する当該飛行力学方程式が特定されるステップと、
−  選択された上記飛行パラメータと関連した上記特定された飛行力学方程式に基づき、少なくとも選択された上記飛行パラメータを含んだ入力パラメータの値を受け取る拡張カルマンフィルタが構成されるステップと、
−  上記航空機の飛行中に、上記拡張カルマンフィルタが実行されて、該フィルタが、出力として、選択された上記飛行パラメータの合同推定値を供するステップと
が遂行されることを特徴としている。
【００１２】
こうして、本発明により、特別に特定された飛行力学方程式によって相互に関係付けられた少なくとも２個の飛行パラメータは適切に構成された拡張カルマンフィルタによって同時に推定される。上記入力パラメータの１つまたは複数の値の誤りはこうして回避される。実際、選択された飛行パラメータの推定中に使用される、上記拡張カルマンフィルタの入力パラメータの誤った値は、本発明により、とりわけ、特定された上記飛行力学方程式の媒介によって結合された、選択された他方の飛行パラメータの推定に基づき、選択された上記飛行パラメータの推定を再構成することができるために、もはや必ずしも選択された飛行パラメータの誤った推定を招来するものではない。
【００１３】
換言すれば、本発明は、フィルタの１つまたは複数の入力値に誤りがあっても、飛行力学方程式によって相互に関係付けられた、前もって選択された少なくとも２個の飛行パラメータの共同推定を可能にするものである。これによって、航空機航法で使用される選択された上記飛行パラメータの値の可用性が保証される。したがって、本発明は、上記航空機のセンサがもはや通常の方法でそれを可能としない場合に、選択された上記飛行パラメータのバックアップ推定値の提供を可能とするものである。
【００１４】
好ましくは、上記カルマンフィルタは以下の行列、即ち、
−  測定ノイズに関係すると共に測定ノイズＶの対角行列と関連した共分散行列Ｒと、
−  展開ノイズに関係すると共に展開ノイズＷの対角行列と関連した共分散行列Ｑと
によって定義され、
以下の付加的ステップ、即ち、
−  上記入力パラメータの値が許容し得るものである旨が検証されるステップと、
−  入力パラメータの値に欠陥が検出された場合に、上記共分散行列ＲおよびＱの少なくともいずれか一方の少なくともいずれか１つの要素の現在値がリアルタイムで適合化されるステップと
が遂行される。
【００１５】
こうして、ゲインの適合化を通じて上記拡張カルマンフィルタの設定を修正して、上記入力パラメータの値に影響する−上記搭載センサによって得られた−一定の測定値の(万一、外部の現象が上記航空機に搭載のセンサの作動に影響を及ぼす場合の)誤りを斟酌することが可能である。上記共分散行列ＲおよびＱの値を適合化することにより、上記搭載センサによって行われる測定または上記拡張カルマンフィルタの推定に、より多くの信頼が置かれることになる。
【００１６】
好ましくは、
−  予備的ステップにおいて、それぞれ測定ノイズと展開ノイズとに関係する上記共分散行列ＲおよびＱの複数の設定値が定義され、こうして定義された上記設定値はそれぞれ上記入力パラメータのうちのいずれか１個のパラメータの欠陥がある値に関連し、かつ、
−  入力パラメータの値に欠陥が検出された場合、測定ノイズと展開ノイズとに関係する上記共分散行列ＱおよびＲの現在値を適合化すべく、検出された欠陥がある値に対応する所定の設定値が測定ノイズと展開ノイズとに関係した上記共分散行列ＲおよびＱに割り当てられる。
【００１７】
さらに、変形例としてまたは補足例として、上記航空機搭載の１つまたは複数のセンサによって測定された上記入力パラメータのうちのいずれか１個のパラメータの値に欠陥が検出された場合、欠陥がある当該測定値は上記拡張カルマンフィルタの出力として供された対応する推定値によって置換される。
【００１８】
こうして、欠陥がある当該測定値は、選択された上記飛行パラメータの推定にあたり上記拡張カルマンフィルタによってもはや考慮されることはない。
【００１９】
さらに、以下のステップ、即ち、
−  上記カルマンフィルタがそのために推定値を供する選択された上記飛行パラメータのうちの少なくともいずれか１個のパラメータが考慮されるステップと、
−  上記カルマンフィルタの入力パラメータの値のうち、上記航空機に搭載のセンサに由来する、上記当該飛行パラメータに対応するものが選択されるステップと、
−  選択された上記の値の少なくともいずれか１つと結びついた不整合が検出されるステップと、
−  選択された上記飛行パラメータの現在値が、残りの選択された１つまたは複数の値と選択された上記飛行パラメータの推定値とに基づき決定され、他方、検出された不整合な１つまたは複数の値は除外されるステップと
が遂行可能であるのが有利である。
【００２０】
こうして、選択された上記飛行パラメータの現在値の決定は、不整合な１つまたは複数の測定値を参照することなく行われることができる。情報の１項目つまり選択された当該飛行パラメータの推定が加えられ、こうして、１つまたは複数の対応する測定値が使用不可もしくは不整合な場合にあっても、上記パラメータの現在値の可用性が高められる。こうして、選択された飛行パラメータの決定を可能とする値のタイプが多様化され、こうして、上記パラメータと関連した現在値の供給が保証される。関連するすべてのセンサに障害が生じた場合にあっても、選択された上記飛行パラメータの現在値の供給を確かなものとすることにより(この場合、上記拡張カルマンフィルタによって行われる推定を使用することが可能である)、上記当該飛行パラメータを含む航法の変化のリスクが低減される。それゆえ、いわゆる標準航法の可用性が高められ、こうして、上記航空機の制御パフォーマンスの継続性が確保される。
【００２１】
さらに、本発明による方法の実施にあたり、上記拡張カルマンフィルタと関連した状態ベクトルは以下の１２個の状態、即ち、
−迎角αと、
−地球座標系における速度νと、
−対地速度Ｖと、
−縦揺れ割合ｑと、
−姿勢θと、
−高度ｈと、
−有効原動推力ＴＢと、
−縦加速度ｎｘに対応する方向に投射されたバイアスｂｎｘと、
−横加速度ｎｙに対応する方向に投射されたバイアスｂｎｙと、
−地球座標系(ｘ、ｙ、ｚ)のｘ軸に沿った風速Ｗｘと、
−地球座標系(ｘ、ｙ、ｚ)のｙ軸に沿った風速Ｗｙと、
−地球座標系(ｘ、ｙ、ｚ)のｚ軸に沿った風速Ｗｚと
によって定義される。
【００２２】
さらに、この実施態様によれば、上記拡張カルマンフィルタの上記入力パラメータは、慣性パラメータと、風速パラメータと、上記航空機に固有のパラメータと、搭載モデリングから生ずる中間パラメータとを含む。
【００２３】
さらに、本発明はまた、航空機の飛行中における航空機の飛行パラメータをリアルタイムで決定するための、入力パラメータの値を受け取ることのできる拡張カルマンフィルタを含んだ注目すべきシステムであって、
−  上記拡張カルマンフィルタは、上記入力パラメータに属する上記航空機の推定されるべき少なくとも２個の事前選択された飛行パラメータ間の相互依存関係を定める飛行力学方程式に基づき構成され、かつ
−  上記拡張カルマンフィルタは、上記航空機の飛行中に、選択された上記飛行パラメータの合同推定値を供するように形成されている
ことを特徴とするシステムに関する。
【００２４】
好ましくは、上記カルマンフィルタは以下の行列、即ち、
−  測定ノイズに関係すると共に測定ノイズＶの対角行列と関連した共分散行列Ｒと、
−  展開ノイズに関係すると共に展開ノイズＷの対角行列と関連した共分散行列Ｑと
によって定義され、
かつ、上記システムは以下の手段、即ち、
−  上記拡張カルマンフィルタの上記入力パラメータの値が許容し得るものであるか否かを検証するための手段と、
−  上記検証手段によって入力パラメータの値に欠陥が検出された場合、共分散行列ＲおよびＱの少なくともいずれか一方の現在値をリアルタイムで適合化するための手段と
を含む。
【００２５】
さらに、上記システムは、欠陥がある１つまたは複数の値を、上記拡張カルマンフィルタによって推定されたそれらの値が得られる場合に、該推定値によって置換するための手段を含む。
【００２６】
本発明は、さらに、上記に説明した少なくとも１つのシステムを含んだ航空機に関する。
【図面の簡単な説明】
【００２７】
添付図面は本発明が実施され得る態様を説明したものである。これらの図中、同一符号は同一要素を示す。
【図１】本発明に係る飛行パラメータを決定するためのシステムを示す、航空機の空力操縦表面のための制御系の概略構成図ある。
【図２】図１の飛行パラメータを決定するためのシステムの推定ユニットの概略構成図である。
【好ましい実施形態の説明】
【００２８】
図１には、概略的な形で、航空機(図示せず)の空力操縦翼面２(例えば補助翼)の制御チェーン１が表されており、該制御系において、操縦翼面２の作動にかかわる航空機の飛行パラメータをリアルタイムで決定するためのシステム３が実行される。
【００２９】
同図によって示すように、制御系１は以下を含む：
−  測定された航空機のパラメータ値(例えば、迎角、縦揺れ角、横揺れ角、対気速度、高度等)を標準的な方法で供することのできる、複数の測定センサおよびコンピュータによって形成される搭載情報ユニット４、
−  情報ユニット４によって測定された値を、リンクＬ０を経て、入力値として受け取る、飛行パラメータをリアルタイムで決定するためのシステム３。これについては以下に詳述する、
−  飛行制御系、即ち、深度、補助翼、方位、エンジンによって供されるパワー等を管理するための多数の航法を実行する航空機の自動操縦装置５。この自動操縦装置５の特定の自動操縦モードは、自動操縦モードの採用が対応する航法を起動するように、それぞれの航法と関連させられている。自動操縦装置５は、決定システム３によって決定されて、リンクＬ２を介して接続されたスイッチング手段Ｃを経て伝送される飛行パラメータ値を受け取ることかできる。自動操縦装置はさらに、出力として、操縦翼面２の制御システム７(アクチュエータを含む)向けの制御コマンドを供することができる、
−  飛行中に航空機を制御するための、特に操縦設備(例えば、操縦桿)を含んだ手動操縦手段６。これらの設備は、決定システム３によって決定されて、リンクＬ３を介して接続されたスイッチング手段Ｃを経て伝送される飛行パラメータ値を受け取ることができる。手動操縦手段６はさらに、出力として、操縦翼面２の制御システム７向けの制御コマンドを供することができる、
−  決定システム３、自動操縦装置５および手動操縦手段６にそれぞれリンクＬ１、Ｌ２およびＬ３を介して接続されたスイッチング手段Ｃ。このスイッチング手段Ｃによって、自動操縦モードと手動操縦モードとの切り換えが可能である。このスイッチング手段Ｃは、決定システム３によって決定された飛行パラメータ値を自動操縦装置５または手動操縦手段６に伝送することができる、および
−  操縦翼面２のポジション調節を行うための制御システム７。この制御システム７は、リンクＬ４を経て、自動操縦装置５または手動操縦手段６に由来する制御コマンドを受け取り、続いて、操縦翼面２の傾きを調節することができる。
【００３０】
本発明によれば、本発明の決定システム３は以下を含む：
−  あらかじめ選択された飛行パラメータと、測定ベクトルＺ(ｔ)を定義するこれらの選択された飛行パラメータのセットとを推定するためのユニット８、および
−  測定ベクトルＺ(ｔ)の推定値 
0005982213.tif 000002
(記号“＾”は推定を意味する)を受け取ることのできる、航空機の現在の飛行パラメータ値を決定するためのユニット９。
【００３１】
図２に示すように、推定ユニット８は、状態ベクトルＸ(ｔ)、測定ベクトルＺ(ｔ)および制御ベクトルＵ(ｔ)と関連した拡張カルマンフィルタ１０を含む。このカルマンフィルタ１０は以下の方程式によって定義される：
0005982213.tif 000003
式中、
−Ｆは状態行列である、
−Ｈは測定ノイズと関連した行列である、
−Ｗ(ｔ)は展開ノイズベクトルである、
−Ｖ(ｔ)は測定ノイズベクトルである、そして
−記号“・”は時間に関する導関数を意味している。
【００３２】
カルマンフィルタ１０の代表的な実施形態において、状態ベクトルＸ(ｔ)は以下の１２個の飛行パラメータによって定義される：
−  迎角αと、
−  地球座標系(ｘ、ｙ、ｚ)における航空機の速度νと、
−  航空機の対地速度Ｖと、
−  縦揺れ割合ｑと、
−  縦方向姿勢θと、
−  高度ｈと、
−  有効原動推力ＴＢと、
−  縦加速度ｎｘに対応する方向に投射されたバイアスｂｎｘと、
−  横加速度ｎｙに対応する方向に投射されたバイアスｂｎｙと、
−  地球座標系(ｘ、ｙ、ｚ)のｘ軸に沿った風速Ｗｘと、
−  地球座標系(ｘ、ｙ、ｚ)のｙ軸に沿った風速Ｗｙと、
−  地球座標系(ｘ、ｙ、ｚ)のｚ軸に沿った風速Ｗｚ。
【００３３】
測定ベクトルＺ(ｔ)は、その一部につき、以下の９個の飛行パラメータによって定義される：
−  空力迎角αａと、
−  空力横滑り角βａと、
−  空気の速度Ｖａと、
−  縦揺れ割合ｑと、
−  縦方向姿勢θと、
−  高度ｈと、
−  垂直速度Ｖｚと、
−  標準荷重倍数Ｎｚと、
−  航空機の対地速度Ｖ。
【００３４】
本発明によれば、測定ベクトルＺ(ｔ)を定義する飛行パラメータが選択されるが、それはこれらのパラメータが特別に特定された飛行力学方程式によって相互に関係付けられるからである。このようにして、この特定の飛行力学方程式によって、測定ベクトルＺ(ｔ)を形成する選択されたさまざまな飛行パラメータ間の相互依存関係を定めることが可能になる。
【００３５】
拡張カルマンフィルタ１０は、測定ベクトルＺ(ｔ)を形成する選択された飛行パラメータと関連した特定の飛行力学方程式に基づいて構成されている。フィルタ１０の初期パラメータ化も実施される。
【００３６】
さらに、制御ベクトルＵ(ｔ)は以下の１６個の飛行パラメータによって形成される：
−  横揺れ割合ｐと、
−  縦揺れ割合ｑと、
−  偏揺れ割合ｒと、
−  横揺れ角φと、
−  姿勢θと、
−  縦加速度ｎｘと、
−  横加速度ｎｙと、
−  迎角αと、
−  横滑り角βと、
−  空気の速度Ｖａと、
−  航空機の質量Ｍと、
−  慣性Ｉｙｙと、
−  揚力ＦＺａと、
−  横揺れモーメントＭａと、
−  実有効静原動推力ＴＢＳと、
−  有効原動推力によるモーメントＭＴＢ。
【００３７】
測定ベクトルＺ(ｔ)と制御ベクトルＵ(ｔ)とのパラメータ値のセットは拡張カルマンフィルタ１０の入力値を定義する。
【００３８】
決定システム３は、リンクＬ５を経て受け取った情報ユニット４に由来する測定値のうち、推定ユニット８特に拡張カルマンフィルタ１０の入力値を形成するものを選択するための第１のモジュール１１を含む。換言すれば、第１の選択モジュール１１は、出力として、選択された入力パラメータの測定値を供する。第１の選択モジュール１１は制御ベクトルＵ(ｔ)の測定値を、リンクＬ８を経て直接に、拡張カルマンフィルタ１０に供する。
【００３９】
図２に図示されているように、カルマンフィルタ１０の演算は２つの異なった相を呈する。それらは以下の通りである：
−  推定値
0005982213.tif 000004
および
0005982213.tif 000005
が得られる予測相(ブロック１２によってシンボリックな形で表示)、および
−  測定ベクトルＺ(ｔ)を形成する、リンクＬ９を経て受け取られ、測定されたパラメータ値を使用して、推定相(ブロック１２)中に供された推定値
0005982213.tif 000006
および
0005982213.tif 000007
が補正される更新相(ブロック１３によってシンボリックな形で表示)。
【００４０】
本発明によれば、拡張カルマンフィルタ１０は、上記航空機の飛行中に、測定ベクトルＺ(ｔ)を形成する選択された飛行パラメータの合同推定値(例えば、空力迎え角αａおよび対気速度Ｖａ)をリアルタイムで供するために形成されている。
【００４１】
代表的な当該実施形態において、拡張カルマンフィルタ１０は以下の方程式によって定義される：
0005982213.tif 000008
および
0005982213.tif 000009
式中、
−  インデックスｍは対応するパラメータの測定値を表している、
−  (ｕ、ｖ、ｗ)はそれぞれ地球座標系(ｘ、ｙ、ｚ)における速度、例えば
0005982213.tif 000010
である、
−  (ｕａ、ｖａ、ｗａ)はそれぞれ地球座標系(ｘ、ｙ、ｚ)における対気速度、例えば
0005982213.tif 000011
である。
【００４２】
決定システム３はまた、入力として、測定ベクトルＺ(ｔ)のパラメータの測定値を受け取って、受け取った測定値が所定の許容基準の関数として許容し得るものであるか、あるいは逆に、(例えば、誤りがあるかまたは不整合であるために)欠陥があるものと見なされるか否かを検証する第１の検証モジュール１４を含む。この第１の検証モジュール１４は、例えば、コンパレータとボーターで形成される(図２には図示せず)。これは、出力として、欠陥がある測定値がログされる信号を供することができる。
【００４３】
本実施例において、第１の選択モジュール１１ならびに検証モジュール１４は推定ユニット８の外部に配置されているが、これらは、変形実施例として、ユニット８に組み込むことも可能であろう。
【００４４】
さらに、推定ユニット８は、第１の選択モジュール１１によって選択された測定値と、第１の検証モジュール１４の出力信号と、フィルタ１０によって出力として供された測定ベクトルＺ(ｔ)の推定値 とを、それぞれリンクＬ１０、Ｌ１１およびＬ１２を経て、入力として受け取る、拡張カルマンフィルタ１０のパラメータ化を適合化するためのモジュール１５を含む。
【００４５】
第１の検証モジュール１４によって、測定ベクトルＺ(ｔ)の１つまたは複数の測定値に欠陥が検出された場合には、適合化モジュール１５は、出力として、リンクＬ１３を経て、欠陥がある測定値に対応する、カルマンフィルタ１０によって供された推定値だけでなく、測定ベクトルＺ(ｔ)の許容可能な測定値も供する。したがって、欠陥がある測定値は、拡張カルマンフィルタ１０への入力として、対応する推定値によって置換される。但し、欠陥がある測定値が再び許容可能なものになるや直ちに、適合化モジュール１５は対応する推定値に代えて、新たな測定値を供することができることは言うまでもない。
【００４６】
さらに、拡張カルマンフィルタ１０は以下の行列、即ち、
−  測定ノイズに関係すると共に測定ノイズＶの対角行列と関連した共分散行列Ｒと、
−  展開ノイズに関係すると共に展開ノイズＷの対角行列と関連した共分散行列Ｑによって定義される。
【００４７】
測定ノイズに関係した共分散行列Ｒは以下の関係式
0005982213.tif 000012
［式中、Ｅは数学的期待値を示す］によって定義される。同様に、展開ノイズに関係した共分散行列Ｑは以下の関係式
0005982213.tif 000013
によって定義される。行列ＲおよびＱはそれぞれ、情報ユニット４に由来する測定値および拡張カルマンフィルタ１０によって供された推定値に置かれた信頼度を担っている。
【００４８】
こうして、第１の検証モジュール１４によって、測定ベクトルＺ(ｔ)の１つまたは複数の測定値に欠陥が検出された場合には、適合化モジュール１５は、フィルタ１０と関連して、共分散行列ＲおよびＱの現在値をリアルタイムで適合化することができる。
【００４９】
測定値に誤りが出現すると、欠陥によって影響される行列Ｒの１行または複数の行に高い値が設定され、欠陥によって影響される行列Ｑの１行または複数の行に低い値が設定されることに注目すべきである。実際、行列Ｒの高い値は、カルマンフィルタ１０に達する測定値により多くの信頼が置かれ、そのため、推定値にさらなる信頼が与えられることを示している。
【００５０】
さらに、適合化モジュール１５は、共分散行列ＲおよびＱに関するプリセットペアが保存される記憶装置１６を含む。代表的な実施形態において、各々のプリセットペアは、測定ベクトルＺ(ｔ)の所定の欠陥がある測定値に関連させられている。各々のプリセットペアは、リンクＬ１４を経て、共分散行列ＲおよびＱの選択的適合化を、例えば、欠陥がある測定値が影響を及ぼすこれらの行列の係数に限って調整することによって可能とすることができる。
【００５１】
代表的な実施形態において、２つの共分散行列ＲおよびＱは同時に適合化される。但し、変法として、行列Ｒのみの適合化または行列Ｑのみの適合化を思料することも可能であることはいうまでもない。
【００５２】
行列ＱおよびＲは、拡張カルマンフィルタ１０の初期パラメータ化中に、所定の値によって初期化されることに注目すべきである。
【００５３】
さらに、決定システム３は以下を含む：
−  情報ユニット４(リンク６)から受け取った飛行パラメータの測定値のうち、そのために拡張カルマンフィルタ１０によって対応する飛行パラメータが推定されるものを選択するための第２のモジュール１７(従って、これは、測定ベクトルＺ(ｔ)を形成するパラメータを必要とする)と、
−  入力として、リンクＬ１５を経て、第２の選択モジュール１７によって選択された測定値を受け取る第２の検証モジュール１８で、例えば、コンパレータとボーターで形成されるこの第２の検証モジュール１８は、受け取った測定値が所定の許容基準の関数として許容し得るものであるか、あるいは逆に、(例えば、誤りがあるかまたは不整合であるために)欠陥があるものと見なされるか否かを検証することができる。選択された測定値の少なくともいずれか１つに欠陥が検出される場合には、第２の検出モジュール１８は、出力として、欠陥がある１つまたは複数の測定値がログされる信号を供することができる。
【００５４】
本実施例において、第２の選択モジュール１７および検証モジュール１８は決定ユニット９の外部に配置されているが、これらは、変形実施例として、ユニット９に組み込むことも可能であろう。
【００５５】
さらに、変形実施例として、第１の選択モジュール１１と第２の選択モジュール１７は単に単一の選択モジュールを形成することもできる。同様に、第１の検証モジュール１４と第２の検証モジュール１８は、変形実施例として、互いに統合されて単に単一の検証モジュールを形成することもできる。
【００５６】
さらに、決定ユニット９は、第２の選択モジュール１７によって選択された測定値と、第２の検証モジュール１８によって供された欠陥信号と、推定ユニット８によって供された測定ベクトルＺ(ｔ)の推定値 とを、それぞれＬ１６、Ｌ１７およびＬ７を経て受け取って、測定ベクトルＺ(ｔ)を形成する飛行パラメータの現在値を決定するためのモジュール１９を含む。こうして、測定ベクトルＺ(ｔ)の当該飛行パラメータにつき、上記当該飛行パラメータと関連した少なくともいずれか１つの測定値に欠陥がある場合、モジュール１９は、この飛行パラメータと関連した残りの選択された非欠陥値と前記当該飛行パラメータの推定値とに基づき、前記当該飛行パラメータの現在値を決定し、それを、出力として、Ｌ２およびＬ３を経て、自動操縦装置５および手動操縦手段６に供することができる。したがって、決定モジュール１９もコンパレータおよび／またはボーターを組み入れることが可能である。
【００５７】
このようにして、本発明により、測定ベクトルＺ(ｔ)を形成する上記飛行パラメータのうちのいずれか１個のパラメータに関連した測定値が全面的に利用不能な場合にあっても、推定ユニット８による上記パラメータの推定値に対応する現在値が自動操縦装置５または手動操縦手段６に供されることができるため、それらの装置または手段が入力値を欠くことはない。
【００５８】
さらに、本発明による飛行パラメータ決定システムは空力操縦翼面に関する制御系に関連して説明したが、こうしたシステムはその他の制御系、例えば、エンジンスピードの制御系にも等しく組み込むことが可能であろう。
【図１】
0005982213.tif 000014
【図２】
0005982213.tif 000015