使用numpy和Python使用后退/前进子程序从LU分解中求解系统

问题描述

试图解决使用此LU分解系统。这是到目前为止的代码，出现错误：

import numpy as np     
from copy import copy

L=np.array([(1,0),(-1,1,(2,-1,(-3,2,-2,1)])
U=np.array([(2,3),(0,3)])
b=np.array([(12,-8,21,-26)])

def solve_system_LU(L,U,b):
    x = np.zeros(b.shape)
    n=len(L)
    y=[0 for i in range(n)]
    for i in range(0,n,1):
        y[i]=b[i]/float(L[i][i])
        for k in range(0,i,1):
            y[i]-=y[k]*L[i][k]

    n=len(U)
    x=[0 in range(n)]
    for i in range(n-1,-1):
        x[i]=y[i]/float(U[i][i])
        for k in range (i-1,-1):
            U[i]-=x[i]*U[i][k]

    return(x)

ans=solve_system_LU(L,b)
print(ans)

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）