[Pythonによる科学・技術計算] 連立一次方程式の解法, 数値計算, numpy - Qiita

はじめに Aを行列としたN次元連立一次方程式 $A x = b$ を$x$について解く。 内容 (1) Aの逆行列$A^-1$をnumpyのlinalg.invメソッドにより計算し， 解ベクトル$x=A^{-1} b$を求める。 (2) np.linalg.solveメソッドから直接$x$を計算する。こちらの方が一般的に速い。 コード(1) %%time import numpy as np """ A x = bを満たすxベクトルを linalg.solveを使って求める """ # 行列Aの生成 a1_lis = [1, 0, 0] a2_lis = [0, 2, 0] a3_lis = [0, 0, 4] A_matrix=np.array([a1_lis, a2_lis, a3_lis]) # b vector b = np.array([4,5,6]) # 行ベクトルとしてｂを生