Mathematica 9 is now available
THIS IS DOCUMENTATION FOR AN OBSOLETE PRODUCT.
SEE THE DOCUMENTATION CENTER FOR THE LATEST INFORMATION.
Mathematica >

行列の高度な操作

SingularValueList[m]m の非零の特異値のリスト
SingularValueList[m,k]m の大きい方からk 番目までの特異値
SingularValueList[{m,a}]a についてのm の一般化された特異値
Norm[m,p]mp-ノルム
Norm[m,"Frobenius"]m のフローベニウス(Frobenius)ノルム

行列の特異値とノルムを求める

行列m の「特異値」は固有値m.m*の平方根である.ここで*はエルミート転置を表す.このような特異値の数は行列の小さい方の次元である.SingularValueListは特異値を大きいものから順に並べる.非常に小さな特異値は通常数値的に無意味である.オプション設定Tolerance->t とするとSingularValueListは最大特異値の端数t よりも小さな特異値は除去する.近似数値行列の許容率はデフォルトでゼロよりも若干大きくなっている.
n-次元の空間における単位球内の各点のベクトルにm×n 行列m を掛けるとm-次元の楕円体が得られる.この楕円体の主軸の長さはm の特異値で与えられる.
行列の「2ノルム」Norm[m, 2]は楕円体の最大主軸であり,行列の最大特異値に等しい.これはまた任意の可能な単位ベクトルv についてのm.v の最大の2ノルムの長さでもある.
行列の「p-ノルム」Norm[m, p]は一般にm.vp-ノルムの達することができる最大の長さである.最もよく考慮されているのはp=1p=2p=の場合である.フロベーニウスのノルムNorm[m, "Frobenius"]が考慮されることもある.これはm.m*のトレースの平方根である.
LUDecomposition[m]LU分解
CholeskyDecomposition[m]コレスキー(Cholesky)分解

正方行列の三角形への分解

LinearSolve[m]を使ってLinearSolveFunctionを作ると,これは行列m を三角形に分解することで作用する.そのような形を明示的に得られると便利なことがある.
「LU分解」は任意の正方行列を下三角行列あるいは上三角行列に効果的に分解する.「コレスキー(Cholesky)分解」は任意のエルミートの正定値行列を下三角行列とエルミート共役の積に効果的に分解する.これは行列の平方根を求めることに例えることができる.
PseudoInverse[m]擬似逆行列
QRDecomposition[m]QR分解
SingularValueDecomposition[m]特異値分解
SingularValueDecomposition[{m,a}]一般化された特異値分解

行列の直交分解

行列の逆行列の標準的な定義は,行列が正方行列あるいは特異行列ではない場合には成り立たない.しかし,行列m の「擬似逆行列」m (-1)はそれでも定義することができる.m.m (-1)-I 中のすべての要素の平方の和を最小にすることで設定できる.ここでI は単位行列である.擬似逆行列は一般化された逆行列あるいはムーア・ペンローズ(Moore-Penrose)型の逆行列ともいわれる.これは最小2乗フィットに関連した問題で特によく使われる.
「QR分解」は任意の行列m を積q*.r として表す.ここでq は正規直交行列であり*はエルミートの転置を表す.またrは主対角線以下のすべての要素がゼロの三角行列である.
「特異値分解(SVD)」は多くの数値行列のアルゴリズムの根底にある要素のひとつである.この基本概念は行列mu.s.v*という形式で書き表すことで,ここでs は対角線上にm の特異値を持つ行列,uv は正規直交行列で,v*v のエルミート型の転置である.
JordanDecomposition[m]ジョルダン(Jordan)分解
SchurDecomposition[m]シューア(Schur)分解
SchurDecomposition[{m,a}]一般化されたシューア(Schur)分解
HessenbergDecomposition[m]ヘッセンベルク(Hessenberg)分解

固有値問題に関連する関数

ほとんどの行列は相似変換としてその固有ベクトルの行列を適用することで固有値の対角行列に約すことができる.しかし,こうするのに十分な固有ベクトルがない場合でも,行列を対角線上に固有値とジョルダン(Jordan)ブロックの両方がある「ジョルダン(Jordan)型」に約すことができる.「ジョルダン(Jordan)分解」は一般に任意の正方行列をs.j.s-1として書く.
これよりも数値的により安定しているのは「Schur分解」である.これは任意の正方行列mq.t.q*として書く.ここでq は正規直交行列でありt はブロック上三角形である.これに関連するものとして,正方行列mp.h.p* の形式で表記する「ヘッセンベルク(Hessenberg)分解」が挙げられる.ここで,p は正規直交行列であり,h は主対角の下の対角まで非零要素を持つことができる.
Ask a question about this page  |  Suggest an improvement  |  Leave a message for the team