MIT线性代数课程笔记对应代码-【lecture 4】

矩阵的LU分解

from scipy.linalg import lu
import numpy as np

a = np.array([[2,1], [8,7]])
p, l, u = lu(a)

因为 PA = LU， 所以LU分解返回了三个参数，为了验证是否正确

np.mat(l) * np.mat(u)

发现就是a。（要转换为矩阵再操作，因为原始是array，array相乘是element的）