2025-11-27T07:13:19.186636

Local Dissipativity Analysis of Nonlinear Systems

Strong, Bridgeman

Dissipativity is an input-output (IO) characterization of nonlinear systems that enables compositional robust control through Vidyasagar's Network Dissipativity Theorem. However, determining the dissipativity of a system is an involved and, often, model-specific process. We present a general method to determine the local dissipativity properties of nonlinear, control affine systems. We simultaneously search for the optimal IO characterization of a system and synthesize a continuous piecewise affine (CPA) storage function via a convex optimization problem. To do so, we reformulate the relationship between the Hamilton-Jacobi inequality and the dissipation inequality as an linear matrix inequality (LMI) and develop novel LMI bounds for a triangulation. Further, we develop a method to synthesize a combined quadratic and CPA storage function to expand the systems the optimization problem is applicable to. Finally, we demonstrate that our method will always find a feasible IO characterization and a CPA or quadratic storage function given that the system is strictly locally dissipative.

academic

非線形システムの局所散逸性解析

基本情報

論文ID: 2511.20838
タイトル: Local Dissipativity Analysis of Nonlinear Systems
著者: Amy K. Strong, Leila Bridgeman (Duke University)
分類: eess.SY (Systems and Control), cs.SY
提出日: 2025年11月25日
論文リンク: https://arxiv.org/abs/2511.20838

要約

散逸性(Dissipativity)は非線形システムの入力-出力(IO)特性を記述し、Vidyasagar ネットワーク散逸性定理を通じて組合せロバスト制御を実現する。しかし、システムの散逸性を決定することは複雑で、通常はモデル依存のプロセスである。本論文は、非線形制御アフィンシステムの局所散逸性特性を決定するための汎用的な方法を提案する。本方法は凸最適化問題を通じてシステムの最適なIO特性を探索し、同時に連続区分アフィン(CPA)貯蔵関数を合成する。このため、著者らはHamilton-Jacobi不等式と散逸不等式の関係を線形行列不等式(LMI)として再定式化し、三角分割に対する新規なLMI誤差界を開発した。さらに、二次関数とCPA組合せ貯蔵関数を合成する方法を開発し、最適化問題の適用範囲を拡張した。最後に、システムが厳密に局所散逸である場合、本方法は常に実行可能なIO特性とCPAまたは二次貯蔵関数を見つけることができることを証明した。

研究背景と動機

1. 中核的な問題

入力-出力(IO)安定性理論はロバスト制御において重要であり、組合せ的アプローチを採用する。すなわち、ネットワークの閉ループIO安定性は各サブシステムの粗粒度の開ループIO特性に依存する。しかし、非線形システムの場合、IO特性を決定することは複雑で、通常はモデル依存のプロセスである。

2. 問題の重要性

組合せ的制御: IO理論により、異なるタイプのサブシステム(線形、非線形、時変など)を柔軟に扱うことが可能
ロバスト性: IO特性は様々な不確実性源に対してロバストであるように選択可能
広範な応用: 局所QSR散逸性は受動性、L2ゲイン、円錐セクタなど多くの有名なIO特性を包含

3. 既存手法の限界

線形システム手法の限界:

線形システムの場合、周波数領域解析、HJIから生じるLMIの求解、またはデータ駆動解析によってIO特性を確立可能
KYP補題を通じて、散逸性は正半定二次貯蔵関数の存在と等価であり、LMI最適化で解くことができる
しかし、非線形システムの貯蔵関数は固定された関数形式を欠き、体系的な方法がない

既存の非線形手法の問題:

動的計画法6,14: 反復法は良好な初期化が必要で、収束保証がない
変換法26,33: 変換が有効な限定された状態空間領域でのみ有効で、特定の技法が必要
SOS法25,27,1: 多項式システムにのみ適用可能で、非凸最適化に依存し過度に保守的である可能性
先行CPA法17: HJIとその誤差項は設計変数の多項式であり、保守的な非凸最適化につながり、線形および制御アフィン項を除外

4. 研究の動機

非線形制御アフィンシステムの局所散逸性特性を決定するための汎用的で体系的な方法が必要であり、以下を満たすべき:

広範な非線形システムクラスに適用可能
凸最適化を通じて実現
実行可能性保証を提供
貯蔵関数とIO特性を同時に合成

中核的な貢献

凸最適化フレームワーク: QSRパラメータとCPA貯蔵関数を同時に探索し、HJIが状態空間領域全体で成立することを保証する凸最適化問題を提案
新規なLMI誤差界(第3節):
- 三角分割上のCPA関数に対するLMI制約を課すための新型LMI誤差界を開発(定理10)
- 原点の閉球上で不等式を強制する方法を開発し、二次貯蔵関数の考慮を可能にする(定理11)
HJIのLMI再定式化(第4節):
- HJIをLMIとして再定式化し、先行の最適化問題を凸化
- CPA貯蔵関数がHJIに従う条件を提案(定理12)
- 二次関数とCPA貯蔵関数を組み合わせた合成方法を開発(定理13)
理論的保証(第5節):
- 厳密に局所散逸である非線形制御アフィンシステムに対して、最適化問題は常にCPAまたは二次貯蔵関数の散逸性特性を検証できることを証明(定理15, 16)
数値検証(第6節):
- 3つのシステムで本方法の有効性を実証し、解析的界および既存手法との比較を含む

方法の詳細

タスク定義

システムモデル: マッピング $G(x_0): \mathcal{L}^m_{2e} \to \mathcal{L}^p_{2e}$ を考える。ここで $y = G(x_0)u$ は以下で定義される:

\begin{cases} \dot{x} = f(x) + \bar{G}(x)u \\ y = h(x) + \bar{J}(x)u \end{cases}

ここで $\bar{G}(x) = B + G(x)$ , $\bar{J}(x) = D + J(x)$ , $f(0)=0, G(0)=0, h(0)=0, J(0)=0$ 。

目標: ロバスト正不変集 $\Omega \subset \mathbb{R}^n$ 上でシステムの局所QSR散逸性を決定する。すなわち、貯蔵関数 $V: \Omega \to \mathbb{R}_+$ が存在して:

V(x_0) + \int_{t_0}^{t_1} w(t)dt \geq V(x_1)

ここで供給率は:

w(u,y) = y^\top Q y + 2y^\top S u + u^\top R u

モデルアーキテクチャ

1. 連続区分アフィン(CPA)関数

三角分割: 状態空間領域 $\Omega$ を有限個のn-単体集合 $\mathcal{T} = \{\sigma_i\}_{i=1}^{m_T}$ に分割

CPA関数の定義: CPA関数は $\Omega$ 上で三角分割の頂点 $\mathcal{E}_\Omega$ での値によって一意に決定される。単体 $\sigma_i = \text{co}\{x_{i,j}\}_{j=0}^n$ に対して、関数値は:

W(x) = x^\top X_i^{-1} \bar{W}_i

ここで $X_i$ の第j行は $x_{i,j} - x_{i,0}$ 、 $\bar{W}_i$ の第j要素は $W_{x_{i,j}} - W_{x_{i,0}}$ 。

2. LMI誤差界(重要な革新)

定理10(単体上のLMI誤差界): LMI形式に対して:

M(x) = \begin{bmatrix} \phi(x) & 0 \\ 0 & -I \end{bmatrix} + \text{He}\left(\begin{bmatrix} 0 & \zeta(x)^\top \\ 0 & 0 \end{bmatrix}\begin{bmatrix} 0 & 0 \\ 0 & I \end{bmatrix}\right)

単体 $\sigma$ 上で、 $x = \sum_{j=0}^n \lambda_j x_j$ の場合:

M(x) - \sum_{j=0}^n \lambda_j M(x_j) \preceq E(x)

ここで誤差界行列:

E(x_{i,j}) = \begin{bmatrix} \frac{1}{2}\hat{\phi}_j & * & * \\ 0 & \frac{1}{2}\Pi^{-1} & * \\ \hat{\zeta} & 0 & -2\Pi^{-1} \end{bmatrix}

ここで $\hat{\phi}_j \geq \max_{j} |H_\phi(x_j, x)|$ , $\hat{\zeta}^{(k)} \geq \max_j |H_{\zeta^{(k)}}(x_j, x)|$ , $\Pi$ は対角正定行列。

定理11(球領域上の二次誤差界): 不等式に対して:

\zeta(x)^\top \zeta(x) + \frac{1}{2}(x^\top \theta(x) + \theta(x)^\top x) \leq 0

$x=0$ でのJacobian が以下を満たす場合:

M_\epsilon = \begin{bmatrix} \frac{1}{2}(J_\theta^\top + J_\theta + \hat{\theta}\epsilon I) & J_\zeta^\top & \hat{\zeta}\epsilon \mathbf{1}_{1\times m} \\ * & -I + \frac{1}{2}\Pi^{-1} & 0 \\ * & * & -2\pi I \end{bmatrix} \preceq 0

元の不等式は $B_\epsilon(0)$ 上で成立する。

3. 貯蔵関数合成条件

定理12(純CPA貯蔵関数, B=0, D=0): CPA関数 $V = \{V_x\}_{x \in \mathcal{E}_T}$ が以下を満たす場合:

$V_x \geq 0, \forall x \in \mathcal{E}_x$ (非負性)
$\|\nabla V_i\|_1 \leq l_i$ (勾配界)
$M(x_{i,j}) + E(x_{i,j}) \preceq 0$ (誤差界付きQSR条件)

$V$ は $\mathcal{X}^\circ$ 上でHJI散逸性を満たす。

定理13(組合せ二次-CPA貯蔵関数): 非ゼロの制御アフィン項を持つシステムに対して、組合せ貯蔵関数を定義:

V(x) = \begin{cases} V_p(x) & x \in \Psi \\ \min[V_P(x), V_C(x)] & x \in B_\epsilon(0) \setminus \Psi \\ V_C(x) & x \in \mathcal{X} \setminus B_\epsilon(0) \end{cases}

ここで $V_P(x) = x^\top P x$ は二次関数、 $V_C$ はCPA関数。追加制約により両者の境界での連続性を確保。

技術的革新点

LMI誤差界の構造化設計:
- Young関係とSchur補引理を利用して非対角項を処理可能な界に変換
- 設計変数 $\Pi$ を導入し、最適化で誤差界の緊密性を調整可能
- 原点単体と非原点単体の処理方法を区別
凸化戦略:
- 先行研究17ではHJIと誤差項は設計変数の多項式(非凸)
- 本論文はLMI再定式化を通じて問題を凸最適化に変換
- 鍵は $\nabla V$ を線形項として扱うこと
組合せ貯蔵関数:
- 原点近傍では二次関数を使用して制御アフィン項を処理
- 原点から遠い領域ではCPA関数で柔軟性を提供
- 制約(21a-b)を通じて関数の連続性を確保
理論的完全性:
- 数値方法を提供するだけでなく、実行可能性保証も証明
- 定理15と16は厳密に局所散逸なシステムに対して、本方法が常に解を見つけることを示す

実験設定

数値例

1. 1次元非線形円錐形システム

システム:

\begin{cases} \dot{x} = k_1 x^3 - (k_1 + k_2)x + Bu \\ y = Cx + Du \end{cases}

パラメータ: $k_1=1, k_2=2, B=1, C=1, D=0$

領域: $\mathcal{X} = [-1,1]$ , RPI集 $\Omega = [-0.8092, 0.8092]$

解析的界: $a = -5.8 \times 10^{-9}$ , $b = 0.50$

2. 2次元振り子システム

システム:

\begin{cases} \dot{x}_1 = x_2 \\ \dot{x}_2 = -\sin x_1 - x_2 + u \\ y = x_2 \end{cases}

入力制約: $|u| < 0.1942$

解析的Lyapunov関数: $V(x) = (1-\cos x_1) + \frac{1}{2}x_2^2$ , 解析的ゲイン界 $\gamma \geq 1$

3. 3次元多項式システム

システム:

\begin{cases} \dot{x}_1 = -x_1 - x_3 + x_2 - x_3 x_2^2 \\ \dot{x}_2 = -x_2 x_3^2 - x_2 + u \\ \dot{x}_3 = \frac{1}{2}(x_1 - x_3) \\ y = x_2 \end{cases}

領域: $\mathcal{X} = [-0.5, 0.5]^3$

評価指標

円錐形界: パラメータ $a, b$ で $G_\Omega \in \text{cone}(a,b)$
L2ゲイン: $\gamma$ で $\|y\|_2 \leq \gamma \|u\|_2$
解析的界との誤差: 数値解と理論界の接近度
収束性: 単体数増加に伴う性能改善

比較手法

解析的手法: 手工によるLyapunov関数構築(例1, 2)
SOS手法25: SOSAnalysisツールボックス使用(例3)
異なる単体密度: $\Delta = 0.1, 0.01, 0.001, 0.0001$

実装詳細

凸最適化ソルバーを使用して問題1と問題2を求解
RPI集は既存ツールで決定18,21
三角分割は均等間隔戦略を採用
組合せ貯蔵関数に対して $\epsilon = 2\Delta$

実験結果

主要な結果

1. 1次元円錐形システム(図1)

単体数	上界 $b$	下界 $a$	解析的界との誤差
10	0.535	0.035	大
100	0.510	0.005	中
1,000	0.502	0.001	小
10,000	0.501	$<10^{-3}$	非常に小
100,000	0.500	$\approx 0$	極小

主要な発見:

単体数増加に伴い、数値界は解析的界に急速に収束
比較的少ない単体数(100個)でも合理的な近似を得られる
上界 $b$ は下界 $a$ より速く収束

2. 振り子システム(図2)

単体数	L2ゲイン $\gamma$
100	3.8
500	3.2
1,000	2.8
5,000	2.4
10,000	2.3
50,000	2.221

主要な発見:

解析的界は $\gamma = 1.0$ 、本方法は $\gamma = 2.221$ を発見
保守性が存在するが、単体増加に伴い継続的に改善
SOS手法は不適用(非多項式システム)
図4は状態空間上の合成貯蔵関数の分布を示す

3. 3次元多項式システム

手法	L2ゲイン $\gamma$	単体/多項式次数
本論文手法	1.65	24,576単体
SOS手法21	1.01	線探索による

主要な発見:

本論文手法は見つけたゲインがより保守的(1.65 vs 1.01)
ただし本論文手法は凸最適化保証を提供し、SOSは非凸最適化に依存
本論文手法の適用性がより広い(多項式システムに限定されない)

ケース分析

図3と図4(振り子システム):

図3はRPI集 $\Omega$ (赤線)と三角分割を示す
図4は状態空間上の合成CPA貯蔵関数が滑らかで連続であることを示す
貯蔵関数は原点で最小で、外側に向かって増加し、物理的直感と一致

実験的発見

収束性: 3つの例すべてで、三角分割の細分化に伴い結果がより緊密な界に収束
適用性:
- 非多項式システムに適用可能(振り子システムの正弦項など)
- 異なる次元を処理可能(1D～3D)
- 複数のIO特性を分析可能(円錐形界、L2ゲイン)
保守性のトレードオフ:
- 手工解析的手法と比較して一定の保守性が存在
- しかし自動化された体系的分析フレームワークを提供
- 保守性は単体数増加により削減可能
計算可行性:
- 数万個の単体に対しても最適化問題は求解可能
- 凸最適化は全体最適解を保証

手法タイプ	適用システム	最適化タイプ	貯蔵関数	本論文の改善
SOS25	多項式	非凸	多項式	より広い適用性+凸最適化
先行CPA17	限定非線形	非凸	CPA	凸化+システムクラス拡張
変換26	局所有効	-	変換領域	直接状態空間解析
本論文	制御アフィン	凸	CPA/二次	体系的+保証

結論と議論

主要な結論

方法論的貢献: CPA貯蔵関数合成と局所QSR散逸性決定を凸最適化で同時に行う初の体系的方法を提案
理論的保証: 厳密に局所散逸なシステムに対して、本方法が実行可能解を見つけることを証明(定理15-16)
技術的革新:
- 新規なLMI誤差界により、CPA関数制約を有限個の頂点条件で実現可能
- 組合せ二次-CPA貯蔵関数が適用範囲を拡張
実用的価値: 数値実験により異なる次元とシステムタイプでの方法の有効性を検証

限界

保守性:
- 誤差界の存在により結果は常に保守的
- 保守性の程度は三角分割密度と関数の二階導数界に依存
- 実験は解析的界との差を示す(例: 振り子システム: 2.221 vs 1.0)
計算複雑性:
- 単体数は次元に対して指数的に増加(次元の呪い)
- 3次元システムで既に24,576個の単体が必要
- 高次元システムは計算不可能な可能性
適用範囲:
- $f, G, h, J \in C^2$ (二階連続微分可能)を要求
- RPI集 $\Omega$ の事前決定が必要
- 全体的特性解析には不適用
パラメータ選択:
- 三角分割設計(密度、形状)が結果に影響
- 重み行列 $\Pi$ の選択が保守性に影響する可能性
- 組合せ貯蔵関数の $\epsilon$ 選択は調整が必要
SOS比較:
- 多項式システムで最適化されたSOS手法より保守的である可能性
- ただしSOS は全体最適を保証せず、本論文手法は凸性を保証

今後の方向

保守性削減:
- より緊密な誤差界の開発
- 適応的三角分割戦略
- 高階補間法
適用性拡張:
- 時変システム
- ハイブリッドシステム
- 不確実性システム
計算効率:
- 分散/並列最適化アルゴリズム
- システム構造の疎性を利用
- 次元削減技術
データ駆動版:
- データからシステム動力学の界を学習
- 学習と最適化の統合

深度評価

利点

1. 方法の革新性 (★★★★★)

革新的貢献: CPA関数合成問題の凸化は初めてで、この分野の重要な進展
LMI誤差界: 定理10の構造化設計は行列構造を巧妙に利用し、先行研究の非凸性を回避
組合せ貯蔵関数: 原点での制御アフィン項の問題を優雅に解決

2. 理論的厳密性 (★★★★★)

完全な証明: すべての定理に厳密な数学的証明がある
実行可能性保証: 定理15-16は理論的保証を提供し、単なるヒューリスティック手法ではない
連続性処理: 組合せ貯蔵関数の連続性証明(定理13)は細密で厳密

3. 実験の充分性 (★★★★☆)

多様性: 1D～3D システム、異なるIO特性を網羅
比較: 解析的界およびSOS手法との比較
収束性分析: 単体増加に伴う改善傾向を示す
不足: 計算時間分析と高次元システムテストの欠如

4. 執筆の明確性 (★★★★★)

論理的明確性: 背景→誤差界→貯蔵関数条件→最適化→保証と層状構造
記号規範: 詳細な記号説明部分
可読性: 複雑な数学導出に直感的説明を付加

不足

1. 実験設計の欠陥

実行時間分析の欠如: 異なる規模の問題の計算時間を報告していない
パラメータ感度: $\Pi, \epsilon, \Delta$ などのパラメータの影響を体系的に分析していない
高次元の欠如: 最高3次元システムのみテスト、方法の次元限界を未探索

2. 保守性分析の不足

定量分析の欠如: 保守性の理論的界または経験的公式を提供していない
改善戦略: 実践で保守性を減らす方法の十分な議論がない
トレードオフ分析: 計算コストと保守性のトレードオフが定量化されていない

3. 比較の不完全性

基線の限定: SOS手法との比較は1つの例のみ
最新手法の欠如: 2020年以降の関連研究との比較がない
公平性問題: 異なる手法の計算リソースが統一されていない

4. 実用性の議論

RPI集決定: 外部ツールに依存し、このステップの困難さを未議論
初期化: 最適化問題の初期化戦略が明記されていない
失敗事例: 方法が失敗または性能不良の場合が未議論

影響力評価

1. 学術的影響 (★★★★☆)

理論的貢献: CPA貯蔵関数合成の凸化は重要な理論的進展
引用可能性: 非線形制御およびロバスト制御分野での高い引用が予想される
開拓性: 後続研究に新しいツールと思想を提供

2. 実用的価値 (★★★☆☆)

適用シナリオ: 低次元非線形システムのIO解析
ツール化の可能性: 方法は体系的でソフトウェアツール化が容易
産業応用: 計算複雑性削減後に広範な応用が可能

3. 再現性 (★★★☆☆)

利点: 数学記述が完全で、アルゴリズムが明確
欠点:
- コードが提供されていない
- 実装細部(三角分割生成など)が未詳述
- 最適化ソルバーの選択と設定が未説明

適用シナリオ

最適な場合

低次元非線形システム(n ≤ 3)のIO特性解析
制御アフィンシステムの散逸性検証
理論的保証が必要な応用場面
非多項式システム(SOS不適用の場合)

不適切な場合

高次元システム(n > 5): 計算不可能
リアルタイム応用: 最適化求解時間が長い可能性
全体的特性解析: 方法は局所特性に対応
非滑らかシステム: $C^2$ 連続性を要求

潜在的応用分野

ロボット制御: マニピュレータ、移動ロボットのロバスト制御
電力システム: 局所安定性とダンピング解析
航空宇宙: 航空機姿勢制御
プロセス制御: 化学プロセスの安定性解析

総合評価

次元	評点	説明
革新性	9/10	凸化は重要な突破
理論的厳密性	10/10	証明は完全で厳密
実験の充分性	7/10	時間と高次元テストの欠如
実用的価値	7/10	次元制限
執筆品質	9/10	明確で専門的
総合	8.4/10	優秀な理論研究

参考文献(精選)

12 Hill & Moylan (1976): "The stability of nonlinear dissipative systems" - 散逸性理論の基礎的研究
8,9 Giesl & Hafstein (2012, 2014): CPA Lyapunov関数合成の開拓的研究、本論文の重要な基礎
25 Summers et al. (2013): "Quantitative local L2-gain and reachability analysis" - SOS手法の代表作
17 Lavaei & Bridgeman (2022): 先行CPA貯蔵関数研究、本論文の直接的改善
28 Van Der Schaft (1992): "L2-gain analysis of nonlinear systems" - 非線形システムL2ゲイン解析の古典文献

総括: これは非線形システムの散逸性解析において重要な貢献をした高品質の理論研究である。CPA貯蔵関数合成問題の凸化は重要な革新であり、理論的保証は方法の信頼性を強化している。主な限界は高次元システムの計算複雑性と手工手法との比較における保守性である。本研究は非線形ロバスト制御に価値のある体系的ツールを提供し、制御理論界で大きな影響を与えることが予想される。