时间:2021-05-26
这篇文章主要介绍了python实现迭代法求方程组的根过程解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下
有方程组如下:
迭代法求解x,python代码如下:
import numpy as npimport 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-6steps = 100xk = np.zeros((3, 1)) # initialize parameter settingerrorlist = []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)结果如下:
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。
声明:本页内容来源网络,仅供用户参考;我单位不保证亦不表示资料全面及准确无误,也不保证亦不表示这些资料为最新信息,如因任何原因,本网内容或者用户因倚赖本网内容造成任何损失或损害,我单位将不会负任何法律责任。如涉及版权问题,请提交至online#300.cn邮箱联系删除。
本文实例讲述了javascript基于牛顿迭代法实现求浮点数的平方根。分享给大家供大家参考,具体如下:今天在网上看到一则利用牛顿迭代法求浮点数的平方根的方法,发
本文实例讲述了JS求解三元一次方程组值的方法。分享给大家供大家参考,具体如下://求用js码一段代码求三元一次方程组的值!//a1=y*b1+z*c1-d1*x
应该怎么使用lingo求解线性方程组?相信大家都就都觉的线性方程组很难解,接下来小编就以实例向大家介绍用lingo11优化解线性方程组的具体方法步骤。lingo
本文实例讲述了Python基于高斯消元法计算线性方程组。分享给大家供大家参考,具体如下:#!/usr/bin/envpython#coding=utf-8#以上
本文实例讲述了Python线性方程组求解运算。分享给大家供大家参考,具体如下:求解线性方程组比较简单,只需要用到一个函数(scipy.linalg.solve)