矩阵的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的)
本文详细介绍了矩阵LU分解的概念、过程,并通过实例演示如何使用SciPy库进行分解及验证分解结果。重点强调了分解步骤及其实现,旨在帮助读者深入理解矩阵分解的基本原理及其在数值计算中的应用。
矩阵的LU分解
from scipy.linalg import lu
import numpy as np
a = np.array([[2,1], [8,7]])
p, l, u = lu(a)np.mat(l) * np.mat(u)
900
251
4050

被折叠的 条评论
为什么被折叠?