简单迭代法求解非线性方程组

最新推荐文章于 2021-01-15 12:06:50 发布

leetteel

最新推荐文章于 2021-01-15 12:06:50 发布

阅读量2.0k

点赞数 2

分类专栏：数值分析文章标签：非线性方程组的求解

本文链接：https://blog.csdn.net/weixin_40653652/article/details/112134766

版权

数值分析专栏收录该内容

46 篇文章 16 订阅

订阅专栏

高斯-塞德尔迭代

import numpy as np

x = [0, 0]
err = 1
k = 0
while err > 0.0000001:
    k = k+1
    x[0] = 0.25*(1+x[1]-0.1*np.exp(x[0]))
    x[1] = 0.25*(x[0]-(1/8)*x[0]*x[0])
    err = max(abs(4*x[0]-x[1]+0.1*np.exp(x[0])-1), abs(-x[0]+4*x[1]+x[0]*x[0]/8))
    print('第{}次迭代，误差为{}'.format(k, err))
    print('迭代结果为',  x)

在这里插入图片描述

简单迭代法

import numpy as np

x = [0, 0]
err = 1
k = 0
while err > 0.0000001:
    k = k+1
    x1 = [0, 0]
    x1[0] = 0.25*(1+x[1]-0.1*np.exp(x[0]))
    x1[1] = 0.25*(x[0]-(1/8)*x[0]*x[0])
    x = x1
    err = max(abs(4*x[0]-x[1]+0.1*np.exp(x[0])-1), abs(-x[0]+4*x[1]+x[0]*x[0]/8))
    print('第{}次迭代，误差为{}'.format(k, err))
    print('迭代结果为',  x)