数值分析方程求根的迭代法.ppt

必要性的证明必要性.设迭代xk+1=?(xk)是p阶收敛。迭代两边取极限，由?(x)的连续性可知x*=?(x*)。设p0是满足的最小正整数。由充分性的证明过程可知迭代p0阶收敛。又若p0p，与迭代p阶收敛矛盾p0=p第29页,共43页，星期六，2024年，5月迭代过程的加速设有不动点迭代：设：缺点：每次迭代需计算第30页,共43页，星期六，2024年，5月埃特金算法设：Aitken加速当x收敛到x*时，修正项分子趋于零。第31页,共43页，星期六，2024年，5月一点注记第32页,共43页，星期六，2024年，5月Newton迭代基本思想：将非线性方程线性化设xk是f(x)=0的近似根，将f(x)在xk一阶Taylor展开:，?在xk和x之间。xyx*xkxk+1条件：f’(x)?0第33页,共43页，星期六，2024年，5月Newton迭代Newton法可以看作下面的不动点迭代：其中?’(x*)=0Newton法至少二阶局部收敛定理设f(x)在其零点x*的某个邻域内二阶连续可导且f’(x)?0，则存在x*的某个?邻域N(x*)=[x*-?,x*+?],使得对?x0?N(x*)，Newton法产生的序列以不低于二阶的收敛速度收敛到x*。第34页,共43页，星期六，2024年，5月Newton迭代Newton法也可以看作一类特殊的加速迭代取?(x)=x+f(x)第35页,共43页，星期六，2024年，5月收敛性定理定理设f?C2[a,b]，且f满足(1)f(a)f(b)0(2)对?x?[a,b]，f’(x)?0且f”(x)?0；(3)初始点x0?[a,b]满足f(x0)f”(x0)0；则Newton法产生的序列收敛到f在[a,b]的唯一零点x*。第36页,共43页，星期六，2024年，5月全局收敛性定理定理设f?C2[a,b]，且f满足(1)f(a)f(b)0(2)对?x?[a,b]，f’(x)?0且f”(x)?0；则对任意初始点x0?[a,b]，Newton法产生的序列收敛到f在[a,b]的唯一零点x*。(3)第37页,共43页，星期六，2024年，5月举例（一）例：设计一个二阶收敛算法计算(a0)。解：转化为求x2-a=0的正根Newton迭代：二阶收敛第38页,共43页，星期六，2024年，5月重根情形设x*是f(x)的m(m?2)重根，Newton法是否收敛？Taylor展式Newton迭代：线性收敛。且重数m越高，收敛越慢。第39页,共43页，星期六，2024年，5月提高收敛阶提高收敛速度但m通常无法预先知道！法一：取二阶收敛法二：将求f(x)的重根转化为求另一个函数的单根。构造针对?(x)的具有二阶收敛的Newton迭代：令，则x*是?(x)的单重根。第40页,共43页，星期六，2024年，5月降低初始点的要求例：求sin(x)-x/6=0的正根。Newton下山法：?k?k为数列中满足的最大数。算法7.2(Newton下山法)给定初始点x0，精度要求?1.如果|f(xk)|?,停机，输出xk2.计算，?=13.如果|f(xk+?dk)||f(xk)|，令xk+1=xk+?dk，返回第1步；

否则?折半，重新计算第3步Newton法的收敛依赖于初始点的选取。第41页,共43页，星期六，2024年，5月割线法Newton法的缺点：每步迭代都要计算导数值只需计算函数值，避免计算导数；切线斜率?割线斜率xk-1xkxk+1xk+1切线割线需要两个初始点；收敛比Newton法稍慢，但对初始点要求同样高。第42页,共43页，星期六，2024年，5月割线法公式两点割线法单点割线法第43