计算方法通选课课件第八章.pptVIP

第八章 非线性方程组解法 一、牛顿（Newton）迭代法 为方便起见，我们以二阶方程组为例作介绍，对三阶和三阶以上的方程组，其方法可类推。现在考虑非线性方程组 （8-1） 牛顿迭代法的基本思想是把非线性问题线性化，此思想可用于非线性方程，同样也可用于非线性方程组。 设已知非线性方程组（8-1）的一组初始近似值为(x0,y0)，把f1(x,y)和f2(x,y)都在(x0,y0)附近用二元泰勒公式展开，得 （8-2） 取其线性部分，并令其等于0，得到关于(x-x0)和(y-y0) 的线性方程组 （8-3） 只要系数矩阵的行列式 （8-4） 则方程组（8-3）的解为 （8-5） 其中 （8-6） 这样得到的(x,y)可作为比(x0,y0)更好的近似值(x1,y1) ，即 （8-7） 由此可得解非线性方程组的牛顿迭代法迭代公式 （8-8） 例1 设 用牛顿迭代法求其在(x0,y0)=(1,1)附近的解。 解：先写出和函数及其偏导数有关的三个行列式 由此可得牛顿迭代公式 取 ，计算结果如下表 n xn yn Jn Jxn Jyn 0 1 1 33 113 159 1 4.4242 5.818 6581.2 -7977.2 -11956.4 2 3.2121 4.001 1795.3 -1430.8 -1519.76 3 2.4152 3.155 705.4 -249.79 -112.676 4 2.0611 2.995 491.1 -29.235 2.23585 5 2.0015 2.9997 468.53 -0.7144 0.13 n xn yn Jn Jxn Jyn 6 2.00 2.99 468.00 -0.0005 7.361*10-5 7 2 3 468 -1.97*10-10 3.219*10-11 8 2 3 468 0 0 二、最速下降法 最速下降法是利用方程组（8-1）构造一个非负函数 （8-9） 用求F的零极小值点的方法得到方程组（8-1）的解。 假定F在解的附近具有二阶连续偏导数，并且它的Hessian矩阵行列式 （8-10） 函数F(x,y)的几何意义是一个空间曲面，它与xoy平面相切的点即是它的零极小值点，如下图8-1、图8-2所示。 图8-1 图8-2 对于空间曲面F=F(x,y)，如果用一系列平行于xoy平面的平面 F=常数 相截之，可以得到一族平面曲线，将它们投影于xoy平面上，得到如图8-2的曲线族，称为曲面的等高线族，处在同一条等高线上的F值是相同的，每条等高线对应一个F值，在极小值点(x*,y*)附近，其等高线形成以(x*,y*)为中心的封闭曲线族，且相应的F值由外往里不断下降，当达到(x*,y*)时，F值为零。 对于处在函数F(x,y)定义域D内的任何点，总存在族中的一条等高线通过它。如果从D内的某一点(x0,y0)出发，沿着使F值下降的方向逐步下降F值，一直降到它的零极小值，那么就可以得到所求 问题的解。 我们知道，