本人是湖南大学的一名力学学生,2020年上数值分析课时写下的作业,希望对大家有所帮助,一起进步。
1、下面是自己整理的LU解线性方程组的基本理论。
2、接下来是python实现代码
import numpy as np
def LU(A):
'''
生成值全位0的U矩阵,和单位矩阵L
'''
L = np.eye(len(A))
U = np.zeros(np.shape(A))
for r in range(1, len(A)): # 求U的第一行和L的第一列
U[0, r - 1] = A[0, r - 1]
L[r, 0] = A[r, 0] / A[0, 0]
U[0, -1] = A