第三章迭代法.pptVIP

第1页，共36页，星期日，2025年，2月5日§3.1二分法根的估计二分法第2页，共36页，星期日，2025年，2月5日根的估计引理3.1（连续函数的介值定理）设f(x)在[a,b]上连续，且f(a)f(b)0，则存在x*?(a,b)使f(x*)=0。例3.1证明x3?3x?1=0有且仅有3个实根，并确定根的大致位置使误差不超过?=0.5。解:单调性分析和解的位置选步长h=2?,扫描节点函数值异号区间内有根第3页，共36页，星期日，2025年，2月5日f(x)=x3?3x?1第4页，共36页，星期日，2025年，2月5日二分法(更快的扫描法)条件:设f(x)在[a,b]上连续，f(x)=0在[a,b]上存在唯一解，且f(a)f(b)0。记Step1:Iff(a0)f(x0)0,thenx*?(a0,x0)leta1=a0,b1=x0;Elsex*?(x0,b0)leta1=x0,b1=b0;Letx1=(a1+b1)/2.Stepk:Iff(ak-1)f(xk-1)0,thenx*?(ak-1,xk-1)letak=ak-1,bk=xk-1;Elsex*?(xk-1,bk-1)letak=xk-1,bk=bk-1;Letxk=(ak+bk)/2.收敛性及截断误差分析:第5页，共36页，星期日，2025年，2月5日例3.2x3?3x?1=0,[1,2],精度0.5e-1第6页，共36页，星期日，2025年，2月5日二分法优点算法简单收敛有保证只要f(x)连续缺点对区间两端点选取条件苛刻收敛速度慢难以推广至多维情形第7页，共36页，星期日，2025年，2月5日§3.2迭代法原理迭代法的思想不动点原理局部收敛性收敛性的阶第8页，共36页，星期日，2025年，2月5日迭代法的思想条件:f(x)=0在x0附近有且仅有一个根设计同解变形x=g(x)迭代式xk=g(xk-1),k=1,2,…如果收敛xk?x*,则x*是f(x)=0的根第9页，共36页，星期日，2025年，2月5日不动点原理(迭代过程收敛)定理3.1(不动点原理)设映射g(x)在[a,b]上有连续的一阶导数且满足1o封闭性：?x?[a,b],g(x)?[a,b],2o压缩性：?L?(0,1)使对?x?[a,b],|g(x)|?L,则在[a,b]上存在唯一的不动点x*，且对?x0?[a,b],xk=g(xk-1)收敛于x*。进一步，有误差估计式后验估计先验估计算法设计中迭代结束条件:近似使用|xk-xk-1|?第10页，共36页，星期日，2025年，2月5日不动点原理例3.3x3?x?1=0，[1,2],x0=1.5,第11页，共36页，星期日，2025年，2月5日不动点原理证明步骤解的存在性;解的唯一性;解的收敛性;误差估计式。第12页，共36页，星期日，2025年，2月5日局部收敛性(格式收敛)定理3.2（局部收敛性）设g’(x)连续,则存在充分靠近x*的初值，使迭代收敛于x*。证明：利用定理3.1,取L=具有局部收敛性的迭代计算上不一定收敛，它是否收敛还要看初值是否取的恰当；而不具有局部收敛性的迭代对任何初值都不可能收敛。应用中:近似使用|g(x0)|1判断第13页，共36页，星期日，2025年，2月5日收敛性的阶(局部收敛速度)定义3.1当xk?x*，记ek=x*-xk，若存在实数p，使ek+1/epk?c?0，则称{xk}有p阶收敛速度。线性收敛p=1平方收敛p=2第14页，共36页，星期日，2025年，2月5日收敛性的阶(局部收敛速度)定理3.3设xk=g(xk-1)?x*，则(1)当g(x*)?0时，{xk}线性收敛；(2)当g(x*)=0,而g(x*)?0时，{xk}平方收敛。证明(2)第15页，共36页，星期日，2025年，2月5日3.3Newton迭代法和迭代加速牛顿(Newton)迭代法“迭代－加速”技术（略）第16页，共36页，星期日，2025年，2月5日牛顿(Newton)迭代法原理(1次近似,直线代替曲线)牛顿格式第17页，共36页，星期日，2025年，2月5日Newton法几何意义：切线法切线代替曲线第18页，共36页，星期日，20