解非线性方程python实现弦截法

安心不心安

已于 2023-11-29 15:22:36 修改

阅读量281

点赞数

分类专栏：数值分析文章标签：算法

于 2023-10-15 09:06:49 首次发布

本文链接：https://blog.csdn.net/anxinbuxinan/article/details/133838067

版权

数值分析专栏收录该内容

14 篇文章 7 订阅

订阅专栏

1.基本概念

弦截法基本实现与牛顿迭代法相同，但是将牛顿迭代法中的导数，用差商来代替，即

与牛顿迭代法相比，弦截法更容易实现，但可能会出现收敛慢、振荡或发散的情况。通常可以通过调整初始点的选择、迭代控制和收敛准则等方式来提高算法的可靠性和收敛速度。

2.代码实现

求以下方程在0.5附近的根

$x e^x-1=0$

"""
@Time ： 2023/10/31 0012 13:18
@Auth ： yeqc
"""


import math

x0, x1 = 0.5, 0.6
N = 1000
for i in range(0, N):
    x2 = x1 - ((x1 - math.e ** (-x1)) * (x1 - x0)) / ((x1 - x0) - (math.e ** (-x1) - math.e ** (-x0)))
    print("结果：", x2)
    x0 = x1
    x1 = x2
    if x2 * (math.e ** x2) - 1 < 10 ** (-8):
        break

牛顿迭代法：python实现牛顿迭代法_安心不心安的博客-CSDN博客