微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

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

如何解决使用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 举报,一经查实,本站将立刻删除。

相关推荐


Selenium Web驱动程序和Java。元素在(x,y)点处不可单击。其他元素将获得点击?
Python-如何使用点“。” 访问字典成员?
Java 字符串是不可变的。到底是什么意思?
Java中的“ final”关键字如何工作?(我仍然可以修改对象。)
“loop:”在Java代码中。这是什么,为什么要编译?
java.lang.ClassNotFoundException:sun.jdbc.odbc.JdbcOdbcDriver发生异常。为什么?
这是用Java进行XML解析的最佳库。
Java的PriorityQueue的内置迭代器不会以任何特定顺序遍历数据结构。为什么?
如何在Java中聆听按键时移动图像。
Java“Program to an interface”。这是什么意思?