[論文レビュー] Convex Optimization of Linear Impulsive Rendezvous

Boris Benedikter, Alessandro Zavoli|arXiv (Cornell University)|Dec 17, 2019

Spacecraft Dynamics and Control参考文献 26被引用数 5

ひとこと要約

本稿では、Tschauner-Hempel力学を用いた時間固定線形インパulse的ランデブー問題に対して、凸最適化手法を提案する。空間的インパulse位置のグリッド上での離散化により、非凸な最適制御問題を凸問題に変換する。本手法は最小限の計算時間で最適解へのグローバル収束を保証し、インパulseの数とタイミングを正確に捉えることができ、既知のベンチマークとほぼ完全に一致することが検証された。

ABSTRACT

This note aims at providing a concise and self-contained document that describes a clear and easy-to-understand method, that could be useful for a reader that is approaching the linear-impulsive rendezvous topic for the first time, but that is also flexible enough to accommodate problem variations, such as additional constraints like bounded-magnitude DV, or direction constraints with minor modifications.

研究の動機と目的

• 楕円軌道における時間固定線形インパulse的ランデブー問題を解く信頼性が高く、計算効率の良い手法を開発すること。
• 事前にインパulse数と位置に関する仮定を必要とする従来手法の限界を克服すること。
• 初期値依存性のない推進剤消費（Δv）のグローバル最適化を可能とすること。
• 自律的ランデブーを含む時間的に制限のある宇宙ミッション応用における高い精度と高速収束を確保すること。
• 凸変換におけるグリッド解像度、解の精度、数値的安定性の間のトレードオフを調査すること。

提案手法

• 独立変数（真近点角）を有限個のグリッドノードに離散化することで、元の最適制御問題を凸最適化問題に変換する。
• インパULSE位置は、これらのグリッドノードでのみ発生すると制約され、状態遷移行列を用いて連続時間のダイナミクスを代数的状態遷移の系列に変換する。
• 数値的条件数と安定性を向上させるために、変換された変数（例：スケーリングされた位置・速度）を用いて問題を再定式化する。
• 得られた凸問題は内点法により解かれる。この手法は初期値に依存せず、グローバル最適解への収束を保証する。
• インパULSE数が最適化の過程で暗黙的に決定されることを活用し、繰り返し調整の必要がなくなる。
• 丸め誤差と切り捨て誤差を低減するため、問題のスケーリングと次元なし化が施される。特に細かいグリッドでは顕著である。

実験結果

リサーチクエスチョン

• RQ1 Tschauner-Hempel力学を用いた楕円軌道における時間固定線形インパULSE的ランデブー問題に対して、凸定式化を導出可能か？
• RQ2 グリッドベースのインパULSE位置局在化は、初期値に依存せず最適解へのグローバル収束を可能にするか？
• RQ3 グリッド解像度の選択がΔvコストとインパULSEタイミングの精度にどのように影響するか？
• RQ4 グリッド密度が増加するに従い、計算コストと解の品質の間のトレードオフはいかなるものか？
• RQ5 問題スケーリングと許容誤差制御により、インパULSEの拡散などの数値的問題を緩和可能か？

主な発見

• 凸最適化手法は、円軌道から円軌道へのランデブーにおいて、2インパULSE解を完全に再現し、Arzelierらの基準解と100%の一致を示した。
• 全体のΔvコストは粗いグリッドでも速やかに最適値に収束し、8～16個のグリッドノードでほぼ最適性能を達成した。
• 計算時間はグリッドノード数にほぼ線形に増加するが、凸ソルバの効率性のおかげで、大規模グリッドに対しても非常に低く保たれた。
• 細かいグリッドでは、有限精度とソルバの許容誤差の影響で、インパULSEが複数ノードにまたがって拡散する数値的不安定性が生じた。
• 中程度から粗いグリッドが、精度と計算効率のバランスを最良に保ち、誤差蓄積を回避しながら高い解の忠実性を維持した。
• 収束の保証と低い計算オーバーヘッドのおかげで、本手法はリアルタイム応用に対しても頑健かつ適している。