[論文レビュー] A Pseudospectral Approach to High Index DAE Optimal Control Problems

Harleigh C. Marsh, Mark Karpenko+1 more|arXiv (Cornell University)|Nov 30, 2018

Numerical methods for differential equations参考文献 24被引用数 5

ひとこと要約

本稿は、特にDIDO MATLABツールボックスを通じた擬スペクトル（PS）法が、高インデックス微分代数方程式（DAE）最適制御問題——例えばインデックス3の振り子問題——を、インデックス還元を経ることなく直接解くことができることを示している。これは、かつて長らく信じられてきた「このような問題は本質的に難しい」という考えに挑戦するものである。主な貢献は、PS法が従来の計算的障壁を回避できることの証拠であり、DAEの難易度に関する新しい理論の必要性を示唆している。

ABSTRACT

Historically, solving optimal control problems with high index differential algebraic equations (DAEs) has been considered extremely hard. Computational experience with Runge-Kutta (RK) methods confirms the difficulties. High index DAE problems occur quite naturally in many practical engineering applications. Over the last two decades, a vast number of real-world problems have been solved routinely using pseudospectral (PS) optimal control techniques. In view of this, we solve a "provably hard," index-three problem using the PS method implemented in DIDO, a state-of-the-art MATLAB optimal control toolbox. In contrast to RK-type solution techniques, no laborious index-reduction process was used to generate the PS solution. The PS solution is independently verified and validated using standard industry practices. It turns out that proper PS methods can indeed be used to "directly" solve high index DAE optimal control problems. In view of this, it is proposed that a new theory of difficulty for DAEs be put forth.

研究の動機と目的

• 高インデックスDAE最適制御問題が本質的に難しいのか、それとも特定の計算手法に起因する困難なのかを調査すること。
• ルンゲ＝クッタ（RK）法がインデックス還元なしに解けなかったベンチマーク用のインデックス3 DAE問題に対して、擬スペクトル（PS）法の有効性を検証すること。
• 実世界の工学的問題に見られる特徴的制御（singular arcs）と代数的制約を伴う高インデックスDAEを、PS法が直接処理できることを検証すること。
• 高インデックスDAE最適制御問題を解く際にインデックス還元が必須であるという一般的な常識に挑戦すること。
• 問題の本質的複雑性ではなく、計算手法の依存性に基づいたDAEの難易度に関する新しい理論を提唱すること。

提案手法

• 本研究では、必要条件の導出に共変量マッピング原理（CMP）を用いる擬スペクトル最適制御を実装するDIDO MATLAB最適制御ツールボックスが使用された。
• インデックス3 DAE最適制御問題は、状態変数 $x_1, x_2, x_3, x_4, x_5$ を用いて標準的な状態空間形式で定式化され、$x_5$ は代数的制約 $x_1^2 + x_3^2 = L^2$ のラグランジュ乗数を表す。
• 2次擬スペクトル微分行列が状態および共役状態の軌道をスペクトル精度で近似するために用いられ、最適制御問題の高精度な解法が可能になった。
• 解の妥当性は、コストレートと制御条件（式9および10）の数値的満たしを確認することで検証された。その際、コストレートと制御式のプロットが用いられた。
• 独立した検証として、業界標準の実務的手法を用いて、インデックス還元や問題変換なしにPS解の信頼性を確認した。
• PS法の収束理論は、ストーン＝ヴェイエルシュトラスの定理およびアーヴェラ＝アスコリの定理に基づいており、インデックスや特異性に関する仮定に依存しない。

実験結果

リサーチクエスチョン

• RQ1 擬スペクトル法は、インデックス還元や問題変換なしに高インデックスDAE最適制御問題を直接解くことができるか？
• RQ2 なぜルンゲ＝クッタ法は、擬スペクトル法が成功する高インデックスDAE問題では失敗するのか？
• RQ3 高インデックスDAE最適制御問題の難易度は、問題自体に内在するものなのか、それとも数値的手法の選択によるものなのか？
• RQ4 特徴的制御と高インデックスDAEは、収束性や精度を損なわずに擬スペクトル法で効果的に取り扱えるか？
• RQ5 複雑な現実世界のDAE問題を擬スペクトル法が解く成功事例は、現在のDAE難易度理論の見直しを要請するものか？

主な発見

• DIDOに実装された擬スペクトル法は、インデックス還元や問題変換なしに、ベンチマーク用のインデックス3 DAE最適制御問題を正常に解いた。
• 解は、すべての必要条件——終端透過性および制御最適性——を満たしており、コストレートと制御式の数値的プロットによって確認された。
• DIDOが生成したコストレートは終端時刻でゼロであり、終端透過性条件の満たされていることを確認した。
• 必要条件から導かれた制御則（式9）は、両辺が図8で一致することにより、数値的に検証された。
• 特徴的制御の条件（式10）も、図9に示されるように数値的に満たされており、解が最適性条件と整合していることが確認された。
• 本研究の結論として、適切に実装された擬スペクトル法を用いれば、高インデックスDAE問題は本質的に難しいものではなく、計算手法の依存性に基づいたDAE難易度の新しい理論の必要性が示唆された。