［Octave］固有値と固有ベクトルを求める

eig関数を使う。固有値と固有ベクトルが行列でそれぞれ戻り値になるため、戻り値は2つ指定する。

>> mxaa = [3, 1, 1; 1, 2, 0; 1, 0, 2]
mxaa =
   3   1   1
   1   2   0
   1   0   2
>> [mxv, mxl] = eig(mxaa)
mxv =
  -5.7735e-01  -1.7174e-16   8.1650e-01
   5.7735e-01  -7.0711e-01   4.0825e-01
   5.7735e-01   7.0711e-01   4.0825e-01
mxl =
Diagonal Matrix
   1   0   0
   0   2   0
   0   0   4

この例では固有値と固有ベクトルは3つあり、例えば3番目の固有値と固有ベクトルを取り出すには、以下のようにする。

>> mxl(3, 3)
ans = 4.0000
>> mxv(:, 3)
ans =
   0.8165
   0.4082
   0.4082