非线性方程迭代法.pptVIP

第1页，共22页，星期日，2025年，2月5日迭代法/Fixed-PointIteration/f(x)=0x=g(x)等价变换f(x)的根g(x)的不动点思路从一个初值x0出发，计算x1=g(x0),x2=g(x1),…,xk+1=g(xk),…若收敛，即存在x*使得，且g连续，则由可知x*=g(x*)，即x*是g的不动点，也就是f的根。迭代法几何意义如下第2页，共22页，星期日，2025年，2月5日xyy=xxyy=xxyy=xxyy=xx*x*x*x*y=g(x)y=g(x)y=g(x)y=g(x)x0p0x1p1?x0p0x1p1?x0p0x1p1?x0p0x1p1?第3页，共22页，星期日，2025年，2月5日定理（充分条件）考虑方程x=g(x),g(x)?C[a,b],若(I)当x?[a,b]时，g(x)?[a,b]；(II)?0?L1使得|g’(x)|?L1对?x?[a,b]成立。则任取x0?[a,b]，由xk+1=g(xk)得到的序列收敛于g(x)在[a,b]上的唯一不动点。并且有误差估计式：??(k=1,2,…)且存在极限??????第4页，共22页，星期日，2025年，2月5日证明：①g(x)在[a,b]上存在不动点？令有根②不动点唯一？反证：若不然，设还有，则在和之间。而③当k??时，xk收敛到x*？???第5页，共22页，星期日，2025年，2月5日④?⑤?⑥?可用来控制收敛精度L越收敛越快小注：定理条件非必要条件，可将[a,b]缩小，定义局部收敛性：若在x*的某?领域B?={x||x?x*|??}有g?C1[a,b]且|g’(x*)|1，则由?x0?B?开始的迭代收敛。即调整初值可得到收敛的结果。第6页，共22页，星期日，2025年，2月5日设在区间[a,b]上方程x=φ(x)有根x*,且对一切x∈[a,b]都有|φ’(x)|≥1，则对于该区间上任意x0(≠x*),迭代公式xk+1=φ(xk)一定发散。证明:不可能收敛于0。定理第7页，共22页，星期日，2025年，2月5日改进、加速收敛/acceleratingconvergence/?待定参数法：若|g’(x)|?1，则将x=g(x)等价地改造为求K，使得例：求在(1,2)的实根。如果用进行迭代，则在(1,2)中有现令希望，即在(1,2)上可取任意，例如K=?0.5，则对应即产生收敛序列。第8页，共22页，星期日，2025年，2月5日?Aitken加速：xyy=xy=g(x)x*x0P(x0,x1)x1x2P(x1,x2)一般地，有：比收敛得略快。?Steffensen加速：详见P273第9页，共22页，星期日，2025年，2月5日定理（牛顿法收敛的充分条件）设f?C2[a,b]，若(1)f(a)f(b)0；(2)在整个[a,b]上f”不变号且f’(x)?0；(3)选取x0?[a,b]使得f(x0)f”(x0)0；则Newton’sMethod产生的序列{xk}收敛到f(x)在[a,b]的唯一根。有根根唯一产生的序列单调有界，保证收敛。定理（局部收敛性）设f?C2[a,b]，若x*为f(x)在[a,b]上的根，且f’(x*)?0，则存在x*的邻域使得任取初值，Newton’sMethod产生的序列{xk}收敛到x