python实现迭代法求方程组的根过程解析
这篇文章主要介绍了python实现迭代法求方程组的根过程解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
有方程组如下:
迭代法求解x,python代码如下:
importnumpyasnp
importmatplotlib.pyplotasplt
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))#initializeparametersetting
errorlist=[]
forkinrange(steps):
xk_1=xk
xk=np.matmul(B,xk)+f.T
print('xk:\n',xk)
errorlist.append(np.linalg.norm(xk-xk_1))
iferrorlist[-1]
结果如下:
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持毛票票。
声明:本文内容来源于网络,版权归原作者所有,内容由互联网用户自发贡献自行上传,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任。如果您发现有涉嫌版权的内容,欢迎发送邮件至:czq8825#qq.com(发邮件时,请将#更换为@)进行举报,并提供相关证据,一经查实,本站将立刻删除涉嫌侵权内容。