详细实验指导见上一篇,此处只写内容啦
实验内容: 求解如下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)
运行结果:
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 用 C# 插值字符串处理器写一个 sscanf
· Java 中堆内存和栈内存上的数据分布和特点
· 开发中对象命名的一点思考
· .NET Core内存结构体系(Windows环境)底层原理浅谈
· C# 深度学习:对抗生成网络(GAN)训练头像生成模型
· 为什么说在企业级应用开发中,后端往往是效率杀手?
· 本地部署DeepSeek后,没有好看的交互界面怎么行!
· 趁着过年的时候手搓了一个低代码框架
· 用 C# 插值字符串处理器写一个 sscanf
· 推荐一个DeepSeek 大模型的免费 API 项目!兼容OpenAI接口!