非线性方程组迭代解法DFP方法与BFS方法.doc

《非线性方程组迭代解法DFP方法与BFS方法.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、非线性方程组的拟牛顿法数理学院汪玉霞李锋摘要：本文主要介绍基于非线性方程组的牛顿法的秩1拟牛顿法和秩2拟牛顿法，并对每个方法举出相应的例子来说明在迭代过程中根的变化情况以及根随迭代序列的变化图，并对这两种迭代法进行总结。关键词：牛顿迭代法非线性方程组Broyden拟牛顿法Abstact:ThispapermainlyintroducestwokindofquasiNewtonmethodbasedonNewtoniterativemethodforsolvingnonlinearequationgroup,andforeachm

2、ethodwegivecorrespondingexampletoillustratethevariationofthechangesofrootwithiterativesequence.Keywords:NewtoniterativemethodnonlinearequationgroupBroydenmethod1.非线性方程组的牛顿迭代法设非线性方程组（1.1）其中是实变量，分别是个变量的元函数，且至少有一个是自变量的非线性函数，若记则（1.1）可以等价地记为（1.2）一般地，，当Jacobi矩阵非奇异时，可得（1.3）

3、这就是非线性方程组的（1.3）的Newton迭代法。利用迭代法（1.3）求解非线性方程组（1.1）的近似解，每一步迭代需要解一个线性方程组，且这一步与下一步的Jacobi矩阵的逆矩阵也不相同，因此计算量较大，可以证明Newton迭代法具有二阶收敛速度，但对迭代初始值的要求很高，即充分靠近精确解。一般地，（1.3）可改写为：2.Broyden秩1拟牛顿法一般地，（1.3）可改写为：考虑用比较简单的矩阵来逼近,将迭代公式改写为（2.1）这里依赖于及，为了避免都要重新计算，我们考虑下一步只是对进行修正，它要求满足方程（2.2）成为拟牛

4、顿方程。它表明矩阵关于点及具有“差商”性质，即当时，即为关于点及的差商，但当时并不确定，为此我们限制是由的一个低秩修正矩阵得到的，即（2.3）其中是秩为的修正矩阵，有（2.3.1），（2.3.2），（2.3.3）组成的迭代法就称为拟牛顿法。它通过给定初始近似和矩阵，根据这三式可以逐次计算得到和，从而可以避免每一步都要计算的雅可比矩阵，这样使得计算量减少，这种算法由于（2.3.3）中的有不同的算法，因此可以得到许多不同的拟牛顿法，常有的拟牛顿法是为秩1矩阵和为秩2矩阵的方法。秩1算法是指（2.4）中修正矩阵的秩时的拟牛顿法，因为对

5、任何的阶的秩1矩阵都可以表示为，其中，于是有（2.5）下面主要是选择合适的使得它们满足拟牛顿方程（2.6）若记，有（2.7）则可以得到即若，则得代入（2.2）中得到，（2.8）若取，只要，就有，则（2.9）由此可以得到一种秩1拟牛顿法（2.10）此算法称为布罗伊登（Broyden）秩1方法。3.秩2拟牛顿法现在考虑（3.1）中校正矩阵为秩2矩阵的情形，即的情形，此时，可以表示为其中，都是阶矩阵，将的第一、二列向量分别记作，，的第一、二列向量分别记作，，则（3.2）代入拟牛顿法方程记，即代入：中得或写成（3.3）现取，（3.4）则

6、可得（3.5）显然，如果取，，使得，（3.6）从而拟Newton方程也满足。令（3.7）（3.8）其中是一个实参数，经过整理后可以得到（3.9）选取不同的参数就可以得到不同的公式，从而得到解方程组的不同的秩2拟牛顿迭代法。4.拟牛顿法的实现例1：对非线性方程组，其中，应用Broyden秩1拟牛顿方法，取迭代格式为：其中：计算结果为：迭代次数为详细结果如下表所示：迭代次序x1x2x31.e+000.e+000.e+000.e+002.e+009.e-011.e+006.e-013.e+009.e-011.e+006.e-014.e

7、+009.e-011.e+006.e-015.e+009.e-011.e+006.e-016.e+009.e-011.e+006.e-017.e+009.e-011.e+006.e-018.e+009.e-011.e+006.e-01可以看出，整个迭代过程还是比较平稳的，在迭代进行到第二次的时候，数值解已经比较稳定，说明了这一方法的有效性。例2：对非线性方程组，初值取应用以下秩2拟牛顿方法（BFS方法），取迭代格式为：计算结果为：迭代次数为详细结果如下表所示：迭代次序xyz1.e+0000.e+0000.e+0000.e+000

8、2.e+0004.e-0018.e-002-5.e-0013.e+0005.e-0017.e-003-5.e-0014.e+0005.e-0012.e-004-5.e-0015.e+0005.e-0015.e-007-5.e-0016.e+0005.e-0011