4.2 行列の固有値、固有ベクトルを求める

MATLAB の eigs() を使えば良い。

行列 $a$ の固有値と固有ベクトルを小さい方 (0 に近い方) から $100$ 個求める

» [v,d]=eigs(a,100,0);

norm(v'*a*v-d) とすると小さな値を出力することが確かめられる ( $V^\ast A V=D$ ということである)。