这篇文章主要介绍了python实现迭代法求方程组的根过程解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
有方程组如下:
迭代法求解x,python代码如下:
import numpy as np import matplotlib.pyplot as plt A = np.array([[8, -3, 2], [4, 11, -1], [6, 3, 12]]) b = np.array([[20, 33, 36]]) # 方法一:消元法求解方程组的解 result = np.linalg.solve(A, b.T) print('Result:\n', result) # 方法二:迭代法求解方程组的解 B = np.array([[0, 3/8, -2/8], [-4/11, 0, 1/11], [-6/12, -3/12, 0]]) f = np.array([[20/8, 33/11, 36/12]]) error = 1.0e-6 steps = 100 xk = np.zeros((3, 1)) # initialize parameter setting errorlist = [] for k in range(steps): xk_1 = xk xk = np.matmul(B, xk) + f.T print('xk:\n', xk) errorlist.append(np.linalg.norm(xk-xk_1)) if errorlist[-1] < error: print('iteration: ', k+1) break # 把误差画出来 x_axis = [i for i in range(len(errorlist))] plt.figure() plt.plot(x_axis, errorlist)
结果如下:
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。
免责声明:本站文章均来自网站采集或用户投稿,网站不提供任何软件下载或自行开发的软件!
如有用户或公司发现本站内容信息存在侵权行为,请邮件告知! 858582#qq.com
暂无“python实现迭代法求方程组的根过程解析”评论...
更新动态
2024年05月15日
2024年05月15日
- 松田聖子SeikoMatsuda《SEIKOJAZZ1》2017[WAV]
- 玩家发现《剑星》独特动画:出现概率极低!
- 消息称百度董事长李彦宏亲自开除璩静:对其非常生气
- 《纸片马里奥RPG重置版》新预告公开 5月23日发售
- 翁立友.2010-男性的坚持经典集【豪记】【WAV+CUE】
- 常宽.1994-宽【华星】【WAV+CUE】
- 草蜢.1998-快慢精选30首2CD【宝丽金】【WAV+CUE】
- 韩红/陈萨《沉默的光影》[FLAC/分轨][105.37MB]
- 群星《华纳至尊金曲》2CD[FLAC][865M]
- 周杰伦《JAY》[WAV分轨][624M]
- dnf110装备搭配
- 崩坏星穹铁道2.2美梦往事活动攻略大全 钟表小子美梦往事全流程选择
- 崩坏星穹铁道2.2新增钟表把戏任务怎么过 崩铁2.2梦境迷钟解密方法
- 刘思涵.2013-拥抱你【种子音乐】【WAV+CUE】
- 刘容嘉.2014-我的自卑感【亚神音乐】【WAV+CUE】