[英]LU Decomp without pivoting in Python
我需要创建函数[L,U] = lr(A)
,该函数将无需旋转或使用inv
, lu
等来求解线性方程组即可计算矩阵A的LU反压缩。
现在就尝试找出伪代码并理解它。 关于如何入门的任何想法?
Kiusalaas的“ Python 3工程中的数值方法”是一个很好的资源。 以下是本书中Doolittle分解方法的代码。
# [L][U] = LUdecomp([A])
def LUdecomp(a):
n = len(a)
for k in range(0,n-1):
for i in range(k+1,n):
if a[i,k] != 0.0:
lam = a [i,k]/a[k,k]
a[i,k+1:n] = a[i,k+1:n] - lam*a[k,k+1:n]
a[i,k] = lam
return a
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.