[論文レビュー] Computational Complexity Analysis of Interval Methods in Solving Uncertain Nonlinear Systems

Rudra Prakash, S. Janardhanan+1 more|arXiv (Cornell University)|Mar 20, 2026

Numerical Methods and Algorithms被引用数 0

ひとこと要約

本論文は、未確定非線形系を解く検証済み区間法のアルゴリズムレベルの最悪計算量フレームワークを構築し、区間分割ベースと微分ベースのアプローチを比較し、支配的なコスト要因を特定する。

ABSTRACT

This paper analyses the computational complexity of validated interval methods for uncertain nonlinear systems. Interval analysis produces guaranteed enclosures that account for uncertainty and round-off, but its adoption is often limited by computational cost in high dimensions. We develop an algorithm-level worst-case framework that makes the dependence on the initial search volume $\mathrm{Vol}(X_0)$, the target tolerance $\varepsilon$, and the costs of validated primitives explicit (inclusion-function evaluation, Jacobian evaluation, and interval linear algebra). Within this framework, we derive worst-case time and space bounds for interval bisection, subdivision$+$filter, interval constraint propagation, interval Newton, and interval Krawczyk. The bounds quantify the scaling with $\mathrm{Vol}(X_0)$ and $\varepsilon$ for validated steady-state enclosure and highlight dominant cost drivers. We also show that determinant and inverse computation for interval matrices via naive Laplace expansion is factorial in the matrix dimension, motivating specialised interval linear algebra. Finally, interval Newton and interval Krawczyk have comparable leading-order costs; Krawczyk is typically cheaper in practice because it inverts a real midpoint matrix rather than an interval matrix. These results support the practical design of solvers for validated steady-state analysis in applications such as biochemical reaction network modelling, robust parameter estimation, and other uncertainty-aware computations in systems and synthetic biology.

研究の動機と目的

• 探索初期ボリューム、目標公差、原始コストが検証済み区間法の複雑さに与える影響を定量化する。
• 不確実性の下で定常エンクロージャを包囲する代表的な区間法の最悪ケースの時間と空間の境界を提供する。
• 高次元で不確定な非線形系における主要な計算ボトルネックとコストドライバを特定する。
• 生化学的ネットワークやロバストパラメータ推定のような実用的なソルバ設計への指針を提供する。

提案手法

• 区間演算と検証済みプリミティブ（包含関数、ヤコビ行列、区間線形代数）のコストモデルを定義する。
• 区間二等分、分割+フィルタ、区間制約伝搬、区間ニュートン、区間クラウ百科の最悪ケースの時間/空間境界を導出する。
• Vol(X0)、ε、C_F、C_J、C_J^{-1}、N_it の観点で複雑さを分析する。
• 行列式と逆行列計算を含む区間線形代数が全体コストに与える影響を評価する。
• 区間ニュートンと区間クラウ百科の一階層上のコストを比較する。
• 観測された箱数と実行時間の数値 study を通じて実用的な洞察を提供する。

実験結果

リサーチクエスチョン

• RQ1 不確実性の下で検証済み定常エンクロージャを計算する一般的な区間法の最悪ケースの時間・空間複雑性はどれくらいか。
• RQ2 Vol(X0)、公差ε、検証済みプリミティブのコストが、異なる区間法全体の計算ボトルネックをどのように生み出すか。
• RQ3 大規模な場合における区間ニュートンと区間クラウ百科の一階層上の計算コストはどのように比較されるか。
• RQ4 行列式や逆行列計算のような区間線形代数のコストが、区間ソルバの実用性にどのような影響を与えるか。

主な発見

• 区間二等分の最悪ケース時間計算量は O((C_F + n) Vol(X0)/ε^n) で、空間計算量は O(n Vol(X0)/ε^n)。
• 分割+フィルタは時間で O(m^n C_F)、空間で O(m^n)、各次元あたりの分割数を m とする。
• 区間制約伝搬は、C_con = O(C_F) の仮定の下で、時間 O(m^n N_it (C_con + n))、空間 O(m^n) を与える。
• 区間行列式/逆行列をラプラス展開で計算すると n に対して階乗/指数的となり、代替的な区間線形代数のアプローチの必要性を示唆する。
• 区間ニュートンは複数解に対して最悪ケースで時間 O(N_it n^3 Vol(X0)/ε^n) を持ち、サブボックスごとの収束利点を強調する。
• 区間クラウ百科は、区間ニュートンと同様の一階層上のコストとされ、矩阵の逆算コスト次第で実用的な利点が変わることについて言及されている。