Aitken和Steffensen方法的原理及代码实现

最新推荐文章于 2022-12-20 21:47:42 发布

c4fx

最新推荐文章于 2022-12-20 21:47:42 发布

阅读量2.3k

点赞数

分类专栏： Python数值分析文章标签： python

本文链接：https://blog.csdn.net/weixin_52116519/article/details/124109365

版权

Python数值分析专栏收录该内容

4 篇文章

订阅专栏

原理

在这里插入图片描述

题目

在这里插入图片描述

代码实现

def f1(x): 
    return 7*x**5-13*x**4-21*x**3-12*x**2+58*x+3
def f2(x): 
    return ((13*x**4+21*x**3+12*x**2-58*x-3)/7)**(1/5)
def f3(x): 
    return (13+21/x+12/(x**2)-58/(x**3)-3/(x**4))/7
def f4(x): 
    return ((12*x**2-58*x-3)/(7*x**2-13*x-21))**(1/3)
def f5(x): 
    return ((-58*x-3)/(7*x**3-13*x**2-21*x-12))**(1/2)
def Aitken(x0,f):
    k=0 
    print('第x_{}次迭代的结果为{}'.format(k, x0))
    k += 1
    print('第x_{}次迭代的结果为{}'.format(k, f(x0)))
    while True:
        x = x0-((f(x0)-x0)*(f(x0)-x0)/(f(f(x0))-2*f(x0)+x0))
        k += 1
        print('第x_{}次迭代的结果为{}'.format(k, x))
        if (abs(x-x0)<(1e-5)):  # abs()函数返回绝对值，1e-5为10的负5次方
            break
        x0=f(x0)
    print("方程的解为：")
    print(x)

def Steffensen(x0,f):
    k=0 
    print('第x_{}次迭代的结果为{}'.format(k,x0))
    while True:
        k+=1
        x=x0-((f(x0)-x0)*(f(x0)-x0)/(f(f(x0))-2*f(x0)+x0))
        print('第x_{}次迭代的结果为{}'.format(k, x))
        if (abs(x-x0)<(1e-5)):  # abs()函数返回绝对值，1e-5为10的负5次方
            break
        x0=x
    print("方程的解为：")
    print(x)

f=f1 #修改函数名
#Aitken(1.5,f)
Steffensen(1.5,f)