逆行列 - Qiita

前回 前回では行列をLU分解した。 今回はLU分解した行列の逆行列を求める。 前回作成したlu_decomposition関数は引き続き使用する。 下三角行列Lの逆行列 ソースコード def inverse_l(l): length = len(l) il = np.eye(length) for i in range(1,length): for j in range(i): il[i,j] = -l[i,j] for k in range(j+1,i): il[i,j] -= l[i,k] * il[k,j] return il np.set_printoptions(linewidth=150) a = np.random.rand(5,5) l, u = lu_decomposition(a) il = inverse_l(l) print('L\n', l) print('L^{