LU分解法的python代码

LU分解是一种线性代数中的矩阵分解技术，它将一个矩阵A分解为两个下三角矩阵L和上三角矩阵U的乘积，即A=LU。在Python中，我们可以使用NumPy库来进行LU分解。以下是一个简单的例子：

import numpy as np

def lu_decomposition(matrix):
    # 使用numpy的linalg.lu函数进行LU分解
    L, U = np.linalg.lu(matrix)
    return L, U

# 示例：对一个二维数组进行LU分解
example_matrix = np.array([[4, 1], [2, 5]])
L, U = lu_decomposition(example_matrix)

print("L Matrix:")
print(L)
print("\nU Matrix:")
print(U)

在这个代码里，`np.linalg.lu`函数返回的是两个矩阵，第一个是L（下三角），第二个是U（上三角）。注意，LU分解并不保证U矩阵是严格的上三角矩阵，因为可能存在非零的对角线下方元素。