[論文レビュー] Exact Volume of Zonotopes Generated by a Matrix Pair

Mingwang Zhao|arXiv (Cornell University)|Apr 12, 2020

Point processes and geometric inequalities参考文献 13被引用数 4

ひとこと要約

本稿では、離散時間線形システムにおける行列対 {A, B} によって生成されるゾノトープの正確な体積を再帰的・解析的に計算する手法を提示する。この手法は、有限時刻および無限時刻における可制御領域および到達可能領域への応用を含む。本手法は固有値構造、特に実数で相異なる非負の固有値を活用し、計算複雑性を顕著に低減する。数値実験により、正確な体積一致と、最適化された再帰による桁違いの高速化が確認された。

ABSTRACT

In this article, we define a class of special zonotopes generated by a matrix pair with finite-interval parameters. We discuss the relationship between the volume of these zonotopes and the controllability of one aspect (the volume of the controllable region) of the dynamic systems. We present a corollary and develop an effective recursive method to compute the volume of the special zonotopes. Furthermore, we develop two recursive and analytical volume-computation methods for the finite- and infinite-time controllable regions with real eigenvalues. We conduct numerical experiments to demonstrate the effectiveness of these new volume-computation methods for zonotopes and regions.

研究の動機と目的

• 離散時間線形システムにおける有限区間パrameterを有する行列対 {A,B} によって生成されるゾノトープのクラスを定義・分析すること。
• これらのゾノトープの体積と動的システムの可制御性、特に到達可能領域および可制御領域のサイズとの関係を確立すること。
• 有限時刻および無限時刻における可制御領域および到達可能領域に関連するゾノトープの正確な体積を計算するための効率的で再帰的なアルゴリズムを開発すること。
• 実数で相異なる非負の固有値を有するシステムに対して、固有値構造を活用することで体積計算の計算複雑性を低減すること。
• 標準的な行列式に基づく計算との比較において、正確性と効率性の両面で妥当性を数値実験により検証すること。

提案手法

• 可制御行列 P_N = [B, AB, ..., A^{N-1}B] からのベクトルのミンコフスキー和として定義されるゾノトープモデル E_q(P_N) を提案。パラメータ c_i ∈ [0,1] を有する。
• P_N のランク q および行列 A の固有値分布に基づく再帰的アルゴリズムを導入。特に、n 個の相異なる実数固有値 λ_i ≥ 0 または λ_i ∈ [0,1) を有する場合を想定。
• 有限時刻における到達可能領域体積を再帰的行列式更新を用いて計算するための補題 1 および定理 2 を導出。これにより、n×n 行列式計算の回数を削減する。
• スペクトル半径 ρ(A) < 1 を仮定し、N→∞ の極限における到達可能領域の極限を分析することにより、無限時刻における可制御領域体積を計算するための定理 3 を開発。
• 式 (19) および (20) の再帰的式を用いて、各ステップで全行列式を再計算せずに体積を計算。これにより計算コストを顕著に低減。
• 可制御行列からの生成子を有するゾノトープに本手法を適用し、可制御性の幾何的解析に向けた正確な体積計算を可能にする。

実験結果

リサーチクエスチョン

• RQ1 有界入力を有する離散時間線形システムにおける行列対 {A,B} によって生成されるゾノトープの体積を、正確かつ効率的に計算する方法は何か？
• RQ2 線形時不変システムにおける到達可能領域／可制御領域の体積と可制御性の間にはどのような関係があるか？
• RQ3 可制御行列から生じるゾノトープの体積計算の計算複雑性を低減するための再帰的アルゴリズムを設計できるか？
• RQ4 特に実数で相異なる非負の固有値を有する行列 A の固有値分布は、体積計算アルゴリズムの効率性および構造にどのように影響を与えるか？
• RQ5 標準的な行列式に基づく体積計算と比較して、提案手法の再帰的アプローチの性能向上（時間および演算回数の観点）はどの程度か？

主な発見

• 提案された再帰的手法（補題 1 および定理 2）は、標準的な行列式ベースの手法（定理 1）と同一の正確な体積を有限時刻における到達可能領域に対して計算し、正しさを確認した。
• N=800 の場合、n×n 行列式計算回数は、定理 1 の 8.501×10^7 回から、補題 1 の 3.188×10^5 回にまで 99% を超える削減が達成された。
• N=800 の場合、計算時間は定理 1 の 3.585×10^2 秒から、補題 1 の 1.692 秒にまで短縮され、約 212 倍の高速化が達成された。
• 無限時刻における可制御領域体積は、わずか 1.328×10^{-3} 秒で計算され、たった 26 回の乗算で達成された（定理 3）。極限ケースにおいて極めて高い効率性を示した。
• スケーリング性能も良好である。N=500 の場合、再帰的手法（定理 2）では 1.243×10^5 回の行列式計算と 0.6586 秒を要したが、標準的手法では 2.071×10^7 回と 85.19 秒を要した。
• 結果は、再帰的アプローチが数値的正確性を維持しながら、特に大きな N 値および高次元システムにおいて計算複雑性を顕著に低減できることを確認した。