2017年11月22日水曜日

学習環境

ということで、行列(Matrix)の固有値(eigenvalue)をSymPy(Python)簡単に求めてみる。

コード(Emacs)

Python 3

#!/usr/bin/env python3
from sympy import pprint, Matrix

M = Matrix([[-2, -7, 2, -5],
            [1, 2, 0, 1],
            [3, 7, -1, 5],
            [1, 3, -1, 3]])

for t in [M, M.eigenvals()]:
    pprint(t)
    print()

入出力結果(Terminal, Jupyter(IPython))

$ ./sample.py
⎡-2  -7  2   -5⎤
⎢              ⎥
⎢1   2   0   1 ⎥
⎢              ⎥
⎢3   7   -1  5 ⎥
⎢              ⎥
⎣1   3   -1  3 ⎦

{0: 2, 1: 2}

$

ただ、これを簡単に求める方法としていいのかどうかは…

手計算した後の検算に役立つということで良しとすることに。

0 コメント:

コメントを投稿