Jacobi在python中求解方程组的方法

问题描述

作为我现在正在进行的项目的一部分，我需要实现上述功能。我的实现无法正常工作，我无法发现错误。这是代码

def jacobi(A,b,x0 = None,iterations = 100):
  if x0==None:
    x0=np.zeros(A.shape[0])
  n = A.shape[0]
  x = x0
  x_prev = x0

  for i in range(iterations):
    for i in range(n):
      s = 0
      for j in range(i):
        s += A[i,j] * x_prev[j]
      for j in range(i+1,n):
        s += A[i,j] * x_prev[j]
      x[i] = (b[i]-s) / A[i,i]
    x_prev = x
  return x

为

x=np.array([[1,1,1],[1,2,3],3,4]])
b=[3,-2]
jacobi(x,b)

它返回一些疯狂的价值观。（解决方案应为[5，-1，-1]）

将寻求任何帮助。谢谢

解决方法

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

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

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

numerical-methods python