非线性方程解法.docVIP

第七章 非线性方程解法 ⒉ 迭代法 ⑴ 不动点的迭代法 类似于线性方程组的情况,先将f(x)=0化成x=φ(x),然后取x0,计算 xn+1=φ(x) n=0,1,2,… 若要求x满足f(x)=0,则x=，称x为函数的一个不动点。 求 f(x)的零点就等价于求（x）的不动点，选择一个初始近似值x0， 将它代入x=（x）右 端，即可求的可以如此反复迭代计算xn+1=φ(x), （x）称为迭代函数，如此得到的序列{x}有极限 则称迭代方程收敛，且x=为的不动点，故称xn+1=φ(x)为不动点迭代法。 迭代法实验⑵ 上机题目：迭代法的实现 实验目的：熟悉迭代法并在计算机上实现 实验要求：①上机前充分准备，复习有关内容，写出计算步骤，查对程序； ②完成实验后写出完整的实验报告，内容应该包括：所用的算法语言，算法步骤陈述，变量说明，程序清单，输出计算结果，结果分析等等；③用编好的程序在Ｍatlab环境中执行。 计算例题: 求f(x)=x3-x-1在[1,1.5]的零点. 取x0=1.5, Matlab程序：(希望学生们以自己的方法去编一个程序，把下面程序看作是理解的摸板，不看，等自己完成了进行比较) 由程序可以这样说明迭代法： function y=diedaifa(f,x0,w,max) f=inline(f); disp(k y) for k=1:max y=feval(f,x0); if abs(y-x0)w break end x0=y; disp(sprintf(%d %10e,k,y)) end k 下面来看个例子: 例2.????????????? 求f(x)=x3-x-1在[1,1.5]的零点. 取x0=1.5,迭代公式为xn+1=(1+xn)1/3,计算结果如下表: N x 0 1.50000000000000 1 1.35720880829745 2 1.33086095880143 3 1.32588377423235 4 1.32493936340188 5 1.32476001129270 6 1.32472594522689 7 1.32471947453436 8 1.32471824544894 9 1.32471801198820 若仅取六位有效数字，那么结果完全相同，这时可以认为x 实际上已满足方程，既为所求的根。若取x0=1.5,迭代公式为xn+1=xn3-1,则算得2.375,12.39,…,等等，表明迭代过程发散所以我们将从其几何解释中观察到迭代收敛的几个条件： ①不动点的几何解释 x=φ(x)的不动点ξ可解释为y= x和y=φ(x)两条曲线的交点.从x0出发向上到达y=φ(x)上点(x0,φ(x0)),由此点再作平行线到达y= x上点(φ(x0),φ(x0)),其横坐标即x1.如此做下去得一条阶梯形或环形折线,或向交点接近(收敛),或远离交点而去(不收敛).不难看出不动点ξ在∣φ′(x)∣1附近收敛,否则不收敛.下段将得到理论的证实. 不动点的收敛性定理 定理２：若对任何x1,x2都有x[a ,b]有ab,∣φ(x1)-φ(x2)∣≤L∣x1-x2∣,L1则迭代xn+1=φ(xn)收敛,在［a ,b］上存在唯一的不动点，对任意　　　x0[a ,b]，由得到的迭代序列｛x｝收敛到φ(x)的不动点，并有误差估计 ∣x-∣≤L/(1-L)∣x1-x0∣. 证明：由∣xn+1-xn∣≤L∣xn-xn-1∣≤L∣x1-x0∣则对任何自然数p有 ∣x- x∣≤∣x- xn+p-1∣+∣xn+p-1- xn+p-2∣+…+∣xn+1- x∣ ≤(L n+p-1+Ln+p-2+…+ L)∣x1-x0∣≤L/(1-L)∣x1-x0∣. 一次性证毕。 注1.????????? 加于φ(x)的条件Lipshitz条件(L称Lipshitz常数)它强与连续性. 注2.????????? 这φ(x)定义在().若是φ：[a, b]→[a, b],定理亦成立. 注3.????????? 这儿得到误差估计照例偏于保守.可在计算时用估计式: ∣xn+1-∣≤L/(1-L)∣x n+1-x n∣ 实践中常用∣φ′(x)∣≤L1,从画的图看的出规律.更好用的是: 1.?????假设在不动点附近连续,1,则xn+1=φ (x n) 在附近收敛。 2.????? 假设φ(x)在不动点附近有直到p阶连续导数并且 φ′()=…= φ(p-1)( )=0, φ(p)( )≠0 则迭代xn+1=φ(x n)在附近p阶收敛。 　　　下面来看个例子： 例3.????????????? 用迭代法求f(x)=4-4x+si