详细实验指导见上一篇,此处只写内容啦

  实验内容: 求解如下4元线性方程组的近似解。

      

Jacobi迭代过程

复制代码
 1 import numpy as np
 2 
 3 A = np.array([[10,-1,2,0],[-1,11,-1,3],[2,-1,10,-1],[0,3,-1,8]])
 4 B = np.array([6, 25, -11, 15])
 5 x0 = np.array([0.0, 0, 0, 0])
 6 x = np.array([0.0, 0, 0, 0])
 7 
 8 times = 0
 9 
10 while True:
11     for i in range(4):
12         temp = 0
13         for j in range(4):
14             if i != j:
15                 temp += x0[j] * A[i][j]
16         x[i] = (B[i] - temp) / A[i][i]
17     calTemp = max(abs(x - x0))
18     times += 1
19     if calTemp < 1e-5:
20         break
21     else:
22         x0 = x.copy()
23 
24 print(times)
25 print(x)
复制代码

    运行结果:

   

 

Gauss-Seidel迭代

 

 

复制代码
 1 import numpy as np
 2 
 3 A = np.array([[10,-1,2,0],[-1,11,-1,3],[2,-1,10,-1],[0,3,-1,8]])
 4 B = np.array([6, 25, -11, 15])
 5 x0 = np.array([0.0, 0, 0, 0])
 6 x = np.array([1.0, 2, -1, 1])
 7 times = 0
 8 
 9 while True:
10     for i in range(4):
11         temp = 0
12         tempx = x0.copy()
13         for j in range(4):
14             if i != j:
15                 temp += x0[j] * A[i][j]
16         x[i] = (B[i] - temp) / A[i][i]
17         x0[i] = x[i].copy()
18     calTemp = max(abs(x - tempx))
19     times += 1
20     if calTemp < 1e-5:
21         break
22     else:
23         x0 = x.copy()
24 
25 print(times)
26 print(x)
复制代码

    运行结果:

   

 

SOR迭代法

复制代码
 1 import numpy as np
 2 
 3 A = np.array([[10,-1,2,0],[-1,11,-1,3],[2,-1,10,-1],[0,3,-1,8]])
 4 B = np.array([6, 25, -11, 15])
 5 x0 = np.array([0.0, 0, 0, 0])
 6 x = np.array([1.0, 2, -1, 1])
 7 w = 1.2
 8 times, MT = 0, 1000
 9 
10 while times < MT:
11     tempx = x0.copy()
12     for i in range(4):
13         temp = 0
14         for j in range(4):
15             if i != j:
16                 temp += x0[j] * A[i][j]
17         x[i] = (B[i] - temp) / A[i][i]
18         x0[i] = x[i]
19     x = w * x + (1-w) * tempx
20     calTemp = max(abs(x - tempx))
21     times += 1
22     if calTemp < 1e-5:
23         break
24 print(times)
25 print(x)
复制代码

    运行结果:

   

 

 posted on   ぺあ紫泪°冰封ヤ  阅读(1037)  评论(0编辑  收藏  举报
编辑推荐:
· 用 C# 插值字符串处理器写一个 sscanf
· Java 中堆内存和栈内存上的数据分布和特点
· 开发中对象命名的一点思考
· .NET Core内存结构体系(Windows环境)底层原理浅谈
· C# 深度学习:对抗生成网络(GAN)训练头像生成模型
阅读排行:
· 趁着过年的时候手搓了一个低代码框架
· 本地部署DeepSeek后,没有好看的交互界面怎么行!
· 为什么说在企业级应用开发中,后端往往是效率杀手?
· 用 C# 插值字符串处理器写一个 sscanf
· 乌龟冬眠箱湿度监控系统和AI辅助建议功能的实现
Live2D
欢迎阅读『数值分析实验之线性方程组的迭代求解(Python实现)』

喜欢请打赏

扫描二维码打赏

了解更多

点击右上角即可分享
微信分享提示