[論文レビュー] Structured Singular Value of a Repeated Complex Full-Block Uncertainty

Talha Mushtaq, Diganta Bhattacharjee+2 more|arXiv (Cornell University)|Nov 11, 2022

Probabilistic and Robust Engineering Design被引用数 4

ひとこと要約

本稿では、流体力学に一般的に見られる繰り返し複素フルブロック不確実性を有するシステムにおける構造的固有値（μ）のより保守的でない上界および下界を計算するための効率的アルゴリズムを提案する。上界には勾配に基づく中心法を、下界には一般化されたパワー反復法を用いることで、Hessian行列の計算を回避し、大規模なシステムへのスケーラビリティを向上させるとともに、繰り返しブロックを非繰り返しとして近似する従来の手法に比べ、より正確な安定性マージンを明らかにする。

ABSTRACT

The structured singular value (SSV), or mu, is used to assess the robust stability and performance of an uncertain linear time-invariant system. Existing algorithms compute upper and lower bounds on the SSV for structured uncertainties that contain repeated (real or complex) scalars and/or non-repeated complex full blocks. This paper presents algorithms to compute bounds on the SSV for the case of repeated complex full blocks. This specific class of uncertainty is relevant for the input output analysis of many convective systems, such as fluid flows. Specifically, we present a power iteration to compute a lower bound on SSV for the case of repeated complex full blocks. This generalizes existing power iterations for repeated complex scalar and non-repeated complex full blocks. The upper bound can be formulated as a semi-definite program (SDP), which we solve using a standard interior-point method to compute optimal scaling matrices associated with the repeated full blocks. Our implementation of the method only requires gradient information, which improves the computational efficiency of the method. Finally, we test our proposed algorithms on an example model of incompressible fluid flow. The proposed methods provide less conservative bounds as compared to prior results, which ignore the repeated full block structure.

研究の動機と目的

• 繰り返し複素フルブロック不確実性を有するシステムにおける構造的固有値（μ）の境界を効率的に計算する手法の不足に取り組むこと。これは、流体力学および対流系において一般的に見られる。
• MATLABのRobust Control Toolboxなどで用いられるように、繰り返し構造を非繰り返しとして近似することによって生じるμ境界の保守的過剰を低減すること。
• Hessian行列の計算を回避し、勾配情報のみを活用することで、大スケールのシステムに適した計算効率の高いアルゴリズムを開発すること。
• 平面コアンエット流れのような系における不安定性メカニズムの物理的解釈を正確に行うには、繰り返しブロック構造を正しくモデル化することが不可欠であることを示すこと。

提案手法

• 上界を半正定値計画（SDP）として定式化し、勾配情報のみを用いて計算効率を高めた中心法により解く。
• パッカードらのパワー反復法を一般化し、繰り返し複素フルブロック不確実性に対するμの下界を計算する。これは、スカラーおよび非繰り返しブロックに対する既存の手法を拡張したものである。
• 上界アルゴリズムを内点法により実装し、Hessian行列の計算を回避する勾配のみの更新を用いる。
• 非圧縮性流体の4×4状態空間モデルおよび合成的な学術的例題に、提案手法を適用して性能と正確性を検証する。
• 繰り返しフルブロックを非繰り返しとして置き換える従来の手法と比較し、保守的過剰が低減することを示す。
• 構造的固有値を用いて時間周波数にわたる入出力ゲインを分析し、物理的不安定性特性を明らかにする。

実験結果

リサーチクエスチョン

• RQ1 繰り返し複素フルブロック不確実性を有するシステムにおける構造的固有値（μ）の上界および下界を、どのように効率的に計算できるか？
• RQ2 繰り返し複素フルブロックを非繰り返しブロックとして近似することで、μ境界および流体流れ系における物理的解釈にどのような影響が生じるか？
• RQ3 勾配に基づく最適化手法は、μ計算における計算効率を向上させつつ、正確性を損なわずに実現可能か？
• RQ4 提案手法は、繰り返しブロック構造を無視する従来の手法と比較して、より正確な時間的挙動および不安定性メカニズムをどのように明らかにするか？
• RQ5 単一の繰り返し複素フルブロックの場合、真のμと凸上界（Dスケーリング）との間の差はどの程度か？

主な発見

• 勾配に基づく中心法を用いた上界計算は、Hessian行列に基づく手法に比べて計算効率が高く、大規模システムへのスケーラビリティを実現する。
• 一般化されたパワー反復法により、繰り返し複素フルブロックに対するよりきつい下界が正確に計算可能であり、非繰り返しとして扱う手法を上回る性能を示す。
• ω = 1.896で、非繰り返し近似によるμ下界は、繰り返し構造を正しくモデル化した真の境界のおよそ1.7倍にのぼり、顕著な保守的過剰が確認された。
• 非繰り返しケースではμ上界のグローバルピークがω > 0に現れたが、繰り返しケースではω < 0に現れた。これは、繰り返し構造を無視すると、時間的不安定性挙動に関する誤った結論を導く可能性があることを示している。
• 提案手法は、既存の手法よりも保守的でない境界を提供し、安定性マージンの推定を改善するとともに、対流系における不安定性メカニズムのより正確な物理的解釈を可能にする。
• 平面コアンエット流れモデルの結果から、流体の入出力ゲインおよび不安定性特性を正しく捉えるには、繰り返しブロック構造をモデル化することが不可欠であることが確認された。