python实现迭代法求方程组的根过程解析

这篇文章主要介绍了python实现迭代法求方程组的根过程解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下

站在用户的角度思考问题,与客户深入沟通,找到合江网站设计与合江网站推广的解决方案,凭借多年的经验,让设计与互联网技术结合,创造个性化、用户体验好的作品,建站类型包括:做网站、成都网站制作、企业官网、英文网站、手机端网站、网站推广、主机域名、网页空间、企业邮箱。业务覆盖合江地区。

有方程组如下:

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)

结果如下:

python实现迭代法求方程组的根过程解析

python实现迭代法求方程组的根过程解析

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持创新互联。

网页标题:python实现迭代法求方程组的根过程解析
文章出自:https://www.cdcxhl.com/article44/ghoohe.html

成都网站建设公司_创新互联,为您提供关键词优化虚拟主机网站改版全网营销推广搜索引擎优化品牌网站制作

广告

声明:本网站发布的内容(图片、视频和文字)以用户投稿、用户转载内容为主,如果涉及侵权请尽快告知,我们将会在第一时间删除。文章观点不代表本网站立场,如需处理请联系客服。电话:028-86922220;邮箱:631063699@qq.com。内容未经允许不得转载,或转载时需注明来源: 创新互联

微信小程序开发