Cht7非线性方程求根选编.ppt

§1 方程求根与二分法 迭代法的几何意义 交点的横坐标 迭代法需解决的三个问题 迭代函数的构造 由迭代函数产生的解序列的收敛性 序列的收敛速度和误差估计 §3 迭代收敛的加速方法 Steffensen迭代格式几何解释 §4 牛顿法 （2） 牛顿下山法. §5 弦截法与抛物线法 7.5.2 抛物线法 7.6 解非线性方程组的牛顿迭代法 ，它不满足条件（4.10）. 通过 逐次取半进行试算，当 时可求得 . 此时有 ，而 显然 . 由 计算 时 , 均能使条件（4.10） 成立. 计算结果如下 ： 即为 的近似. 一般情况只要能使条件（4.10）成立， 则可得到 ，从而使 收敛. 四、重根情形 1.5 1.416666667 1.414215686 1.414213562 (2) 1.5 1.458333333 1.436607143 1.425497619 (1) x0 x1 x2 x3 xk 1.5 1.411764706 1.414211438 1.414213562 0 1 2 3 (3) k 图7-5 弦截法与切线法（牛顿法）都是线性化方法，但两者 有本质的区别. 切线法在计算 时只用到前一步的值 ，而弦截法 （5.2），在求 时要用到前面两步的结果 ，因 此使用这种方法必须先给出两个开始值 . 例10 用弦截法解方程 解 设取 作为 开始值，用弦截法求得的结果见表7-8， 比较例7牛顿法的计算结果可以看出， 弦截法的收敛速度也是相当快的. 实际上，弦截法具有超线性的收 敛性. 设已知方程 的三个近似根 ，以 这三点为节点构造二次插值多项式 ，并适当选取 的一个零点 作为新的近似根，这样确定的迭代过程称 抛物线法，亦称密勒（Müller）法. 在几何上，这种方法的基本思想是用抛物线 与 轴的交点 作为所求根 的近似位置（图7-6). 图7-6 插值多项式 有两个零点： （5.3） 式中 问题是该如何确定 ，假定在 三个近 似根中， 更接近所求的根 ，为了保证精度，选 （5.3） 中较接近 的一个值作为新的近似根 . 为此，只要取 根式前的符号与 的符号相同. 例11 用抛物线法求解方程 解 设用表7-8的前三个值 作为开始值，计算得 故 代入（5.3）式求得 以上计算表明，抛物线法比弦截法收敛得更快. 在一定条件下可以证明，对于抛物线法，迭代误差有 下列渐近关系式 可见抛物线法也是超线性收敛的，其收敛的阶 , 收敛速度比弦截法更接近于牛顿法. 从（5.3）看到，即使 均为实数， 也 可以是复数，所以抛物线法适用于求多项式的实根和复根. 考虑方程组 （6.1） 其中 均为 的多元函数. 用向量记号记 , (6.1）就可写成 （6.2） 当 ，且 中至少有一个是自变量 的非线性函数时，称方程组（6.1）为非线 性方程组. 非线性方程组求根问题是前面介绍的方程（即 ） 求根的直接推广，只要把前面介绍的单变量函数 看 成向量函数 则可将单变量方程求根方法推广到方程 组（6.2）. 若已给出方程（6.2）的一个近似根 ， 将函数 的分量 在 用多元函数泰 勒展开，并取其线性部分，则可表示为 令上式右端为零，得到线性方程组 （6.3） 其中 （6.4） 称为 的雅可比(Jacobi)矩阵. 求解线性方程组（6.3），并记解为 ，则得 （6.5） 这就是解非线性方程组（6.2）的牛顿迭代法. 例12 求解方程组 给定初值