someng的博客

考虑n阶线性方程组：高斯消去法pythonimport numpy as npdef gaussin(a, b): m, n = a.shape c = np.zeros(n) for i in range(n): #限制条件 if (a[i][i] == 0): print("no answer") #k...

LU分解：将矩阵A分解成一个下三角矩阵L和一个上三角矩阵U的乘积 A=LU，称为对矩阵A的LU分解。LU分解代码：import numpy as npdef my_LU(b): a = np.array(b) n = len(a) L = np.zeros([n,n]) U = np.zeros([n,n]) for k in range(n-1)...