条件数

条件数(じょうけんすう、: condition number)は、問題のコンピュータでの数値解析しやすさの尺度であり、その問題がどれだけ数値解析に適しているかを表す。条件数が小さい問題は「良条件 (well-conditioned)」であり、条件数が大きい問題は「悪条件 (ill-conditioned)」である。

行列の条件数

たとえば A x = b {\displaystyle Ax=b} という方程式の条件数は、 x {\displaystyle x} を近似的に求める際の不正確さの上限を与える。なお、これには丸め誤差の影響は考慮しない。条件数は行列の属性であって、計算に使うシステムの浮動小数点数の精度やアルゴリズムとは無関係である。この場合(非常に大まかに言って)、 b {\displaystyle b} の変化によって解である x {\displaystyle x} が変化する率が条件数である。従って、条件数が大きければ b {\displaystyle b} の小さな誤差も x {\displaystyle x} の大きな誤差となって現れる。一方、条件数が小さければ、 x {\displaystyle x} における誤差は b {\displaystyle b} における誤差より大きくなることはない。

より正確に条件数を定義すると、 x {\displaystyle x} の相対誤差を b {\displaystyle b} の相対誤差で割った最大比率である。

b {\displaystyle b} の誤差を e {\displaystyle e} とする。すると解 A 1 b {\displaystyle A^{-1}b} の誤差は A 1 e {\displaystyle A^{-1}e} となる。解の相対誤差と b {\displaystyle b} の相対誤差の比率は、次のようになる。

A 1 e / A 1 b e / b {\displaystyle {\frac {\Vert A^{-1}e\Vert /\Vert A^{-1}b\Vert }{\Vert e\Vert /\Vert b\Vert }}}

これは容易に次のように書き換えられる。

( A 1 e / e ) ( b / A 1 b ) {\displaystyle (\Vert A^{-1}e\Vert /\Vert e\Vert )\cdot (\Vert b\Vert /\Vert A^{-1}b\Vert )}

b {\displaystyle b} e {\displaystyle e} がゼロでないとき)その最大値は明らかに2つの作用素ノルムの積となる。

κ ( A ) = A 1 A {\displaystyle \kappa (A)=\Vert A^{-1}\Vert \cdot \Vert A\Vert }

同様の定義は、任意の行列ノルムに当てはまる。この数は数値線型代数学にはよく使われるので、行列の条件数 (condition number of a matrix) と名づけられている。

もちろん、この定義はノルムの選択に依存している。

  • {\displaystyle \|\cdot \|} l 2 {\displaystyle l_{2}} ノルムなら、
    κ ( A ) = σ m a x ( A ) σ m i n ( A ) {\displaystyle \kappa (A)={\frac {\sigma _{\mathrm {max} }(A)}{\sigma _{\mathrm {min} }(A)}}} であり、ここで σ m a x ( A ) {\displaystyle \sigma _{\mathrm {max} }(A)} A {\displaystyle A} の最大特異値 σ m i n ( A ) {\displaystyle \sigma _{\mathrm {min} }(A)} は最小特異値である。したがって、
    • A {\displaystyle A} 正規なら
      κ ( A ) = | λ m a x ( A ) λ m i n ( A ) | {\displaystyle \kappa (A)=\left|{\frac {\lambda _{\mathrm {max} }(A)}{\lambda _{\mathrm {min} }(A)}}\right|} λ m a x ( A ) ,   λ m i n ( A ) {\displaystyle \lambda _{\mathrm {max} }(A),\ \lambda _{\mathrm {min} }(A)} はそれぞれ A {\displaystyle A} の最大および最小固有値
    • A {\displaystyle A} ユニタリなら
      κ ( A ) = 1 {\displaystyle \kappa (A)=1\,}
  • {\displaystyle \|\cdot \|} l {\displaystyle l_{\infty }} ノルムで、 A {\displaystyle A} 三角行列で特異値を持たない(すなわち、 a i i 0 i {\displaystyle a_{ii}\neq 0\;\forall i} )なら
    κ ( A ) max i ( | a i i | ) min i ( | a i i | ) {\displaystyle \kappa (A)\geq {\frac {\max _{i}(|a_{ii}|)}{\min _{i}(|a_{ii}|)}}}

それ以外の条件数

特異値分解の条件数、多項式の解を求める際の条件数、固有値の条件数など、様々な問題について条件数を定義できる。

一般に数値問題が良設定なら、それを関数 f {\displaystyle f} で表すことができ、 m {\displaystyle m} -タプルの実数 x {\displaystyle x} から n {\displaystyle n} -タプルの実数 f ( x ) {\displaystyle f(x)} への写像となる。

するとその条件数は、問題の定義域における解の相対誤差とデータの相対誤差の比の最大値と定義される。

max { | f ( x ) f ( x ) f ( x ) | / | x x x | : | x x | < ϵ } {\displaystyle \max \left\{\left|{\frac {f(x)-f(x^{*})}{f(x)}}\right|\left/\left|{\frac {x-x^{*}}{x}}\right|\right.:|x-x^{*}|<\epsilon \right\}}

ここで ϵ {\displaystyle \epsilon } は問題のデータの変化における何らかの適度に小さい値である。

f {\displaystyle f} が可微分であれば、次のように近似できる。

| f ( x ) f ( x ) | | x | {\displaystyle \left|{\frac {f'(x)}{f(x)}}\right|\cdot \left|x\right|}

関連項目

外部リンク

  • Condition Number of a Matrix at Holistic Numerical Methods Institute
  • Matrix condition number on PlanetMath
連立一次方程式
ベクトル
ベクトル空間
計量ベクトル空間
行列線型写像
演算・操作
不変量
クラス
行列式
多重線型代数
数値線形代数
基本的な概念
ソフトウェア
ライブラリ
反復法・技法
人物
行列値関数
その他
カテゴリ カテゴリ