LINEスタンプ制作代行サービス・LINEスタンプの作り方!

お電話でのお問い合わせ:03-6869-8600

stampfactory大百科事典

コレスキー分解

コレスキー分解(コレスキーぶんかい)とは、正定値エルミート行列"A"を下三角行列"L"と"L"の共役転置"L"との積に分解すること。Aのエルミート性を利用したLU分解の特別な場合である。"L"の対角成分は実数にとることができて(符号・位相の自由度があるが)通常は、対角成分を正の実数に採り、その場合には、"L"が一意に定まる。"アンドレ=ルイ・コレスキー"にちなんで名づけられた。"A"が実対称行列の場合、上式の共役転置は転置に単純化される。エルミート対称行列Aが正定値であることと、Aのコレスキー分解が存在することは同値になる。コレスキー法はガウスの消去法の改良版である。ガウスの消去法は"A"の左方から順次"L"を作用させ前進消去(LU分解に対応。)するが、"A"="L"A"    ( または、"A"="L"A" )コレスキー法は"A"を順次"L"と"L"で挟んで前進消去していくと考えればよい。"A"="L"A" "L"   ( または、"A"="L"A" "L" )このとき"A"のエルミート性は保たれる。詳細は以下の解説を参照。"A"が実行列の場合は単純に、エルミート→対称、共役転置→転置と読み替えればよい。コレスキー分解の再帰的アルゴリズムでは、まず最初に"A"を下のように置く(定義する)。以下、i回目のステップ。エルミート性を保ちながら"A"のi行とi列を前進消去して"A"を生成することを考える。 "A"はi-1行・列まで前進消去されたエルミート行列であるので、下式のように書ける。ここで、Iはi-1次の単位行列、aはi番目の対角要素、bはi列目の下三角部、Bは、Aのi+1行・列以降の部分でやはりエルミートある。次に、Liをと定義するとAは、と書ける。( A = Li A Li より。)ここで、Aはである。( 注.ここでb b は、行列の積。 )以上が、i回目のステップ。AからAが計算出来たことになる。nをAの次数として、このステップをn回繰り返すと"A" = Iとなりコレスキー分解は終了する。であり、と置くと、(これが最終的に求める"L"である。)であることが確認できる。コレスキー・バナキエヴィッツ法 は直接下三角行列 "L" の各エントリを計算するための式を与える。行列 "L" の左上隅から始め行ごとに計算を進める。コレスキー・クラウト法 はコレスキー・バナキエヴィッツ法とは少し異なる方法で、下三角行列Lの各エントリを計算する。すなわち、行列Lの左上隅から始め列ごとに計算を進める。使用する計算式はコレスキー Banachiewicz 法と同一である。上述した分解法では、計算に平方根演算を用いるため,分解後の行列Lに無理数が現れることが普通であり、コレスキー分解の結果を利用した後の計算が面倒となる。そこで、この欠点を解消するために考え出された方法が修正コレスキー分解である(改訂コレスキー分解と呼ぶことがある)。修正コレスキー分解では、"A"が正定値行列である必要はなくA=LDLの形に平方根演算を使わずに四則演算だけを用いて分解の計算を行なうことができる。ここで、Dは対角行列で、行列Lの対角成分はすべて1とする。注意:対称行列は正則であってもその修正コレスキー分解が存在しない場合がある(たとえば対角要素が0で非対角要素が1である2次の対称行列は、正則でも修正コレスキー分解が存在しない例である)。正定値行列であれば,分解は必ず存在する。不完全コレスキー分解は、修正コレスキー分解により行列Aをと分解するところ、行列Lを後の計算が簡略化されるものに変更し、と分解する手法である(Nはある行列で,この不完全分解の残差の行列である)。共役勾配法(傾斜法)の前処理の1つとして採用されることがある。

出典:wikipedia

LINEスタンプ制作に興味がある場合は、
下記よりスタンプファクトリーのホームページをご覧ください。