LU分解python

LU分解是一种将一个矩阵分解为一个下三角矩阵(L)和一个上三角矩阵(U)的方法。在Python中，可以使用SciPy库中的linalg.lu函数进行LU分解。以下是一个示例代码：

import numpy as np
from scipy.linalg import lu

A = np.random.rand(4, 4)
p, l, u = lu(A)

在这个例子中，我们生成了一个随机的4x4矩阵A，并使用lu函数进行LU分解。分解的结果是一个置换矩阵p、一个下三角矩阵l和一个上三角矩阵u。通过np.allclose函数可以验证分解的正确性。