基于Python实现Jacobi迭代法、Gauss-Seidel迭代法和SOR法方程组(一)

标题基于Python实现Jacobi迭代法方程组计算(一)

在这里插入图片描述
一、首先根据习题的方程求出迭代矩阵如下:

x11=round(-0.4x2-0.2x3-2.4,7)
x22=round(0.25x1-0.5x3+5,7)
x33=round(-0.2x1+0.3x2+0.3,7)

二、整个的代码如下:

x1=0
x2=0
x3=0
print('Jacobi')
for i in range(20):
    x11=round(-0.4*x2-0.2*x3-2.4,7)
    x22=round(0.25*x1-0.5*x3+5,7)
    x33=round(-0.2*x1+0.3*x2+0.3,7)
    print('|',i+1,'|',x11,'|',x22,'|',x33,'|')
    x1=x11
    x2=x22
    x3=x33

三、运行结果如下:

D:\Users\anaconda3\envs\tools1\python.exe D:/pythonProject/Jacobi迭代计算公式.py
Jacobi
| 1 | -2.4 | 5.0 | 0.3 |
| 2 | -4.46 | 4.25 | 2.28 |
| 3 | -4.556 | 2.745 | 2.467 |
| 4 | -3.9914 | 2.6275 | 2.0347 |
| 5 | -3.85794 | 2.9848 | 1.88653 |
| 6 | -3.971226 | 3.09225 | 1.967028 |
| 7 | -4.0303056 | 3.0236795 | 2.0219202 |
| 8 | -4.0138558 | 2.9814635 | 2.013165 |
| 9 | -3.9952184 | 2.9899536 | 1.9972102 |
| 10 | -3.9954235 | 3.0025903 | 1.9960298 |
| 11 | -4.0002421 | 3.0031292 | 1.9998618 |
| 12 | -4.001224 | 3.0000086 | 2.0009872 |
| 13 | -4.0002009 | 2.9992004 | 2.0002474 |
| 14 | -3.9997296 | 2.9998261 | 1.9998003 |
| 15 | -3.9998905 | 3.0001675 | 1.9998938 |
| 16 | -4.0000458 | 3.0000805 | 2.0000284 |
| 17 | -4.0000379 | 2.9999744 | 2.0000333 |
| 18 | -3.9999964 | 2.9999739 | 1.9999999 |
| 19 | -3.9999895 | 3.0000009 | 1.9999915 |
| 20 | -3.9999987 | 3.0000069 | 1.9999982 |
Process finished with exit code 0

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值