计算方法 | Jacobi迭代法求线性方程组

很简单，类似简单迭代法。

import numpy as np
# Jacobi迭代法
# 测试数据 增广矩阵
data = [[10,-1,-2, 7.2], [-1,10,-2,8.3], [-1,-1,5, 4.2]]
print(np.array(data))
# 取初始值为0,0,0
x = [[0] for i in range(len(data))]
_x = [[0] for i in range(len(data))]
print(np.array(x))
# 开始迭代
print("开始迭代：")
for n in range(10):
    print("round{}:".format(n))
    for i in range(len(data)):
        _sum = data[i][-1] / data[i][i]
        for j in range(len(data[0])-1):
            if i != j:
                _sum -= (data[i][j] / data[i][i]) * x[j][0] # 用-=来移项
        _x[i][0] = _sum
    for i in range(len(x)):
        x[i][0] = _x[i][0]
    print(x)

over。