非线性方程的求根方法3.ppt

Other Kinds of Newton Method Newton 迭代法几个变形 简化Newton法（平行弦法） 弦截法，又称割线法，弦位法，线性插值法 Newton下山法 Newton法对重根的求解 抛物线法（muller法），又称二次插值法 代数多项式求根 简化Newton法（平行弦法） 若| ? `(x)|=|1-cf `(x)|1，即取0cf `(x)2在x*附近成立，则收敛。 若取c=1/f `(x0) or 1/f `(xk)，则称简化Newton法。 弦截法（割线法），弦位法 在Newton迭代格式中，用差商近似导数， 弦截法的几何意义： 例 用简化Newton法和弦截法计算方程 x3-3x+1=0的根。 Newton下山法 为防止Newton法发散，可增加一个条件： |f(xk+1)||f(xk)|，满足该条件的算法称下山法。 可用下山法保证收敛，Newton法加速收敛。 下山因子选取法 ?的选取：从?=1开始，逐次减半计算，即按 例 求解方程 用Newton下山法，结果如下： Newton法对重根的求解 设f(x)=(x-x*)m g(x) ,m ?2，m为整数，g(x*)?0,则x*为方程f(x)=0的m重根。此时有 f(x*)=f `(x*)=……= f(m-1) (x*)=0, f(m) (x*) ?0 方法三：还可令 则 故x*是μ(x)=0的单根，对μ(x)用Newton法，可得 它是平方阶收敛的。 代数多项式求解 Splitting Method 劈因子法 秦九韶算法 设n次代数方程 用Newton迭代法求有限区间的实根，则要计算 ，一般采用秦九韶算法。 由Taylor展式 同理 比较x的同次幂系数得： 故代数方程的Newton迭代公式 代数方程的Newton迭代法算法 迭代公式： （c?0,k=0,1,……） 迭代函数： 称弦截法. 得 x y x* xk+1 xk-1 Pk-1 y=f(x) xk Pk 弦线PkPk-1的方程： 当y＝0时， 证明略，因弦截法非单步法，不能用前述定理判别 证明参考（关治，陆金甫 《数值分析基础》）。 解 设f(x)=x3-3x+1，则f `(x)=3x2-3 由简化的Newton法，得 由弦截法，得 x0=0.5 x1= 0.3333333333 x2 = 0.3497942387 x3 = 0.3468683325 x4 = 0.3473702799 x5 = 0.3472836048 x6 = 0.3472985550 x7 = 0.3472959759 x8 = 0.3472964208 x9 = 0.3472963440 x10 = 0.3472963572 x11 = 0.3472963553 x0=0.5; x1=0.4; x2 = 0.3430962343 x3 = 0.3473897274 x4 = 0.3472965093 x5 = 0.3472963553 x6 = 0.3472963553 简化Newton法 弦截法 要达到精度10-8，简化Newton法迭代11次，弦截法迭代5次（ Newton迭代法迭代4次）。 无论前面哪种迭代法： （Newton迭代法、 简化Newton法、 弦截法） Newton迭代法 是否收敛均与初值的位置有关。 例 x0 = 2 x1 = -3.54 x2 = 13.95 x3 = -279.34 x4 = 122017 x0 =1 x1 = -0.5708 x2 = 0.1169 x3 = -0.0011 x4 = 7.9631e-010 x5 = 0 收敛 发散 （0??1,下山因子） 记 称Newton下山法。 即 的顺序，直到使下降条件|f(xk+1)||f(xk)|成立为止。 要求达到精度|xn-xn-1|≤10-5，取x0=-0.99. 解：先用Newton迭代法：f `(x)=x2-1 x2=21.69118 x3=15.15689 x4 = 9.70724 x5 = 6.54091 x6 = 4.46497 x7 = 3.13384 x8 = 2.32607 x9 = 1.90230 x10= 1.75248 x11= 1.73240 x12= 1.73205 x13= 1.73205 需迭代13次才达到精度要求 k=0 x0 =-0.99