scipy.linalg.detを使ってみる

とにかく使ってみる

行列の逆行列は以下で得られる。

今回は二次正方行列の逆行列を求める。

from scipy import linalg
import numpy as np

A = np.array([[1,2],[3,4]])
A_tilde = np.array([[4,-2],[-3,1]])

A_tilde/linalg.det(A)
#array([[-2. ,  1. ],
#       [ 1.5, -0.5]])

linalg.inv(A)
#array([[-2. ,  1. ],
#       [ 1.5, -0.5]])

転置の性質

行列式の性質として、がある。
折角なので、デカい行列で試してみる。

A = np.random.randint(1, 10, (10000, 10000))
linalg.det(A) == linalg.det(A.T)
#True

数秒で答え返ってきた。PCすごい！