第六章非线性方程求根试卷.pptVIP

m重根 若 2. 根的搜索 (1) 图解法（利用作图软件如 Matlab) §3 牛顿法 一、牛顿法的迭代公式 考虑非线性方程 3. 迭代法收敛阶的判别 定理6.3 如果 在 附近的某个领域内有p( )阶 连续导数，且 四、求m重根的牛顿法-修正牛顿法 修正格式一：构造迭代格式 (6.4.2) 注：重数m的确定 则迭代格式（6.4.2) 至少2阶收敛。 设 令 则 修正格式二： 令 则 是 的单根， 即 构造迭代格式 (6.4.3) 则迭代格式（6.4.3) 至少2阶收敛。 由于Newton迭代法的收敛性依赖于初值 的选取，如果 离方程的根 较远，则Newton迭代法可能发散。为了防止迭 代发散，可以将Newton迭代法与下山法结合起来使用，放宽 初值 的选取范围，即将（6.4.1）式修改为： 其中， 称为下山因子，选择下山因子时，希望 满足下 山法具有的单调性，即 这种算法称为Newton下山法。 在实际应用中，可选择 。 五、牛顿法的变形 1、牛顿下山法 2.牛顿下山法的计算步骤： （1）选取初始近似值x0； （2）取下山因子? = 1； （3）计算 （4）计算f (xk+1)，并比较 与 的大小，分以下二种情况 1）若 ，则当 时，取x*? xk+1，计算过程结束；当 时，则把 xk+1 作为新的 近似值，并返回到（3）。 2）若 ，则当?≤??且|f(xk+1)| ，取x*? xk，计算过程结束； 否则若?≤??，而 时，则把xk+1加上一个适当选定的小正数， 即取xk+1+?作为新的xk值，并转向（3）重复计算；当?＞??且 时，则将下山因子缩小一半，取?/2代入，并转向（3）重复计算。 【例9】求方程 f(x) = x3 – x – 1 = 0 的根。 x0 =0.6 牛顿下山法的计算结果： k ? xk 0 1 0.6 1 1/25 1.14063 2 1 1.36681 3 1 1.32628 4 1 1.32472 牛顿迭代法每迭代一次都需计算函数值 和导数值 计算量比较大；且迭代过程中计算 时，仅利用了 点的信息， 而没有充分利用已经求出的 ；在导数计算比较麻烦 或难以求出时， 迭代格式构造 (2) 构造方法：将Newton迭代格式中的导数用差商代替。 1、割线法： (1) 构造思想：用割线的斜率代替牛顿迭代法中切线的斜率； 设法避开导数值的计算，因此可以采用离散牛顿法（割线法）。 一个自然的想法就是在充分利用“旧信息”的同时， §4 弦截法与抛物线法 （1）割线法的几何意义 x0 x1 切线 割线 切线斜率　?　割线斜率 x2 割线法迭代格式： （2）割线法的局部收敛性与收敛速度 设 ， 在 ,且 的某一邻域 内二阶连续可微，当 时，由割线法产生的序列 收敛于 ，且收敛阶至少为1.312。 割线法是超线性收敛的！ 2、 双点弦截法 ： 切线斜率　?　割线斜率 初值 x0 和 x1。 x0 x1 x2 双点弦截法的局部收敛性与收敛速度 设 ， 在 ,且 的某一邻域 内二阶连续可微，且对任意的 则 当 时，由弦截法产生的序列 收敛于 且收敛阶至少为1.618。 定理6.5 ， 双点弦截法也是超线性收敛的！ 计算结果 k xk xk-xk-1 0 0.5 1 0.6 0.1 2 0.56532 -0.03468 3 0.56709 0.00177 4 0.56714 0.00005 从表中可以看出弦截法的收敛速度也是相当快的，迭代到第4步就得到精度 的结果. 例10 用弦截法解方程 由迭代格式（3） 取初值 得 结果精确到四位有效数