数值分析第七章 非线性方程与方程组的数值解法0607.ppt

三、迭代收敛的加速方法 埃特金加速收敛方法 加速方法 三、迭代收敛的加速方法 斯特芬森迭代法 加速方法 将埃特金加速技巧与不动点迭代结合 四、牛顿法 牛顿法的基本算法 设 x*是方程 f (x)=0 的根, xk是x*的近似值.将 f (x)在 xk 附近展开, 有 四、牛顿法 牛顿法的基本算法 牛顿迭代法的几何解释 x y x* x0 x 1 x 2 二、牛顿法应用举例 k xk 0 1 2 3 0.5 0.57102 0.56716 0.56714 k xk 0 1 2 3 4 10 10.750000 10.723837 10.723805 10.723805 四、牛顿法 收敛性 牛顿法就是对 的不动点迭代法 单根情形( f (x*)≠0 ) 因此牛顿迭代法局部平方收敛， 四、牛顿法 收敛性 牛顿法就是对 的不动点迭代法 单根情形( f (x*)≠0 ) 例 用牛顿法求f(x)=x2-a (a0)的根的收敛性及收敛速度. 解 f (√a)=2√a≠0，平方收敛 收敛速度 四、牛顿法 收敛性 牛顿法就是对 的不动点迭代法 重根情形( f (x*)=0 ,...,f (m-1)(x*)=0, f (m)(x*)≠0) 因此牛顿迭代法局部线性收敛， 四、牛顿法 收敛性 牛顿法就是对 的不动点迭代法 重根情形( f (x*)=0 ) 例 用牛顿法求f(x)=xm的根. 解 收敛速度？ 线性收敛. 四、牛顿法 收敛性 牛顿法就是对 的不动点迭代法 重根情形( f (x*)=0 ,...,f (m-1)(x*)=0, f (m)(x*)≠0) 四、牛顿法 习题： 1）求f(x)=sin x+x2cos x-x2-x的根x*=0的重数，并估计牛顿法收敛到6位正确小数所需步数(设x0=1). 四、牛顿法 简化牛顿法和牛顿下山法 简化牛顿法 牛顿法每步需计算f (xk),计算量大，且有时得不到，简化牛顿法可避免计算导数. 四、牛顿法 简化牛顿法和牛顿下山法 牛顿法下山法 牛顿法只有局部收敛性，若初始点偏离x*较远，牛顿法可能发散，牛顿法下山法可防止牛顿法发散. 四、牛顿法 简化牛顿法和牛顿下山法 牛顿法下山法 牛顿法只有局部收敛性，若初始点偏离x*较远，牛顿法可能发散，牛顿法下山法可防止牛顿法发散. 五、弦截法与抛物线法 避免牛顿法中对导数的计算，同时具有超线性收敛速度的两种算法 五、弦截法与抛物线法 弦截法 算法 五、弦截法与抛物线法 弦截法 收敛性估计 可以证明，在弦截法收 敛到x*，且f (x*)≠0的假设 下，近似误差关系 成立，也就意味着 弦截法超线性收敛. 五、弦截法与抛物线法 抛物线法 算法 其中，为以xk.xk-1,xk-2为插值节点的插值多项式P2(x)的零点. 抛物线法也适用于求解多项式的实根和复根. 五、弦截法与抛物线法 抛物线法 收敛性估计 迭代估计有下列渐近表达式 可见抛物线法也是超线性收敛的。 六、求根问题敏感性与多项式零点 求根问题敏感性与病态代数方程 求根问题敏感性 求根问题：f (x)=y，已知y，求x. 条件数： 若|f (x*)|非常小，求根问题病态. 六、求根问题敏感性与多项式零点 求根问题敏感性与病态代数方程 病态代数方程 对多项式 若系数有微小扰动，对根影响很大，则p(x)为病态代数方程. 条件数：令 例12 多项式 解 取 的根 由（6.4）可得 实际上，方程 的根 分别为 这说明方程是严重病态的. 六、求根问题敏感性与多项式零点 多项式的零点 使用牛顿法求解 对p(x)=0使用牛顿法计算出第一个根，计算过程使用秦九韶算法计算多项式及其导数值； 对 使用牛顿法计算第二个根； ... 为二次多项式可直接求解. 六、求根问题敏感性与多项式零点 多项式的零点 使用抛物线法求解 对p(x)=0使用抛物线法计算出第一个根x1； 若x1为实数， ，若x1为复数， 对q(x)再使用抛物线法重复以上过程，直到解 出所有根. 六、求根问题敏感性与多项式零点 多项式的零点 使用矩阵特征值算法 p(x)=0的根为以下矩阵的特征多项式的根 第七章