(9)--2.2 一般迭代法(1)数值计算方法.pptVIP

迭代法的理论/*TheoryofIterationMethod*/一、不动点迭代/*Fixed-PointIteration*/f(x)=0x=g(x)等价变换f(x)的根思路从一个初值x0出发，计算x1=g(x0),x2=g(x1),…,xk+1=g(xk),…g(x)的不动点（迭代函数）（迭代格式）若收敛，即存在x*使得.若g连续，则由可知x*=g(x*)，即x*是g的不动点，也就是f的根。

看起来很简单，令人有点不敢相信！！那么需要注意什么问题呢——收敛？？？如何保证这种方法是收敛的呢？

例2.2已知方程在上有一个根（正根）下面选取5种迭代格式：即即即即即12345

取计算结果如下：法1法4法3法2法5

可以更换为Lipschitz条件成立考虑方程x=g(x),若(I)当x?[a,b]时，g(x)?[a,b]；(II)?0?L1使得对?x?[a,b]成立。则任取x0?[a,b]，由xk+1=g(xk)得到的序列收敛于g(x)在[a,b]上的唯一不动点。并且有误差估计式：??(k=1,2,…)且存在极限连续时(k=1,2,…)

的简洁证明：(I)g(x)在[a,b]上有不动点；??(II)g(x)在[a,b]上的不动点是惟一的；(III)由xk+1=g(xk)得到的序列是收敛的；(IV)推导误差估计式。

当x很大时，此处可改为算法:不动点迭代给定初始近似值x0，求x=g(x)的解.输入:初始近似值x0;容许误差TOL;最大迭代次数Nmax输出:近似解x或失败信息。Step1k=1;Step2x=g(x0);Step3whileabs(x-x0)TOLk?NmaxStep4k=k+1;Step5x0=x;Step6x=g(x0);Step7outputthesolutionofequation:xanditerativetimesk

休息一下