标题基于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