Aitken加速法对收敛性的影响

最新推荐文章于 2023-10-15 09:18:41 发布

月月鸟语花香

最新推荐文章于 2023-10-15 09:18:41 发布

阅读量2.4k

点赞数 1

分类专栏：计算方法文章标签：算法

本文链接：https://blog.csdn.net/weixin_40166018/article/details/106249590

版权

本文探讨了Aitken加速法在处理不动点迭代法求解方程f(x)=x^3-2*x-1=0时的作用。通过选取迭代函数x=(x^3-1)/2，原本不收敛的迭代过程在应用Aitken加速法后表现出良好的收敛性。代码展示了Aitken加速法的实现过程，并统计了迭代次数，证明了该方法的有效性。

摘要由CSDN通过智能技术生成

在前文——“不动点迭代法计算三元一次方程及不同迭代函数对收敛性的影响”中，对于函数f(x)=x^3-2*x-1，用不动点迭代法求根，当选取迭代函数x=(x^3-1)/2时，计算结果不收敛。下面用Aitken加速法尝试，看看会发生什么：

题：

用不动点迭代法求方程f(x)=x^3-2*x-1=0的根，用Aitken加速法按等价方程x=(x^3-1)/2实现。要求输出每次的迭代结果并统计所用的迭代次数,取精度c=0.5*10^(-5)，x0=2.

迭代公式如下(迭代公式参考计算方法——江西高校出版社)：

程序流程图见本人前文——“不动点迭代法计算三元一次方程及不同迭代函数对收敛性的影响”，不同之处是在x=g(x)阶段，替换成上图的迭代公式。

代码：

#include <iostream>
#define MAXREPT 100
#define f(x) (x*x*x-1)/2//表示y与x0及z与y之间的函数
#define g(x,y,z) z-(z-y)*(z-y)/(z+x-2*y)//表示x1与x0关系的函数
using namespace std;
#define c 0.5e-5
#include<iomanip>
int main()
{
double x0 = 0, x1 = 2,y=0,z=0;
int i =