第四章 非线性方程.docVIP

第四章 非线性方程（组）的数值求解 求解非线性方程f(x)=0是工程技术上经常遇到的重要数学问题。在绝大多数情况下，由于函数f(x)的复杂性，往往得不到解x*的解析表达式，只能求得满足精度要求的近似数值解。 求解方程f(x)=0的数值法的基本思想是从某个初始近似值x0出发，按照某种数值过程模式进行重复，从而逐次改进前次的结果，直到达到规定的精度，便停止这种重复，并将最后所得结果xn作为方程f(x)=0的数值解。显然数值解是满足了精度要求的近似解。根据这种基本思想，数值法解方程必须做三件事：首先在函数f(x)的定义域内寻找初值x0；其次建立逐次逼近的数值过程模式；第三则是规定最终近似解的精度。关于解的精度取决于工程计算的实际要求而与算法无关。各种算法的实质是由数值过程的模式确定的，不同的模式具有不同的计算效果和速度。初值x0则影响计算过程的速度甚至成败。 4.1 解非线性方程Newton迭代法 4.1.1 Newton法 Newton法的基本思想：将非线性方程线性化，以线性方程的解逐步逼近非线性方程的解。 设f(x)在其零点x＊邻近一阶连续可微，且，当x0充分接近x*时。由Taylor公式有：，以方程作为方程的近似。其解可作为的近似解。重复以上过程得迭代公式：，如图所示，的曲线为过点（x0，f(x0)）的切线，x1为切线与x轴的交点，如此继续下去，xn+1为曲线上点（xn，f(xn)）处的切线与x轴的交点。因此Newton法是以曲线的切线与x轴的交点作为曲线与x轴的交点的近似。故Newton法又称为切线法。 4.1.2 弦截法 在Newton迭代公式中，为了避免计算导数，可用差商近似导数，得到迭代公式：，称为弦截法。由图示，过曲线上两点和的直线为，它与x轴的交点，于是，从几何上看。弦截法是以曲线上两点的割线与x轴的交点作为曲线与x轴的交点的近似，故弦截法又称为割线法。另一方面，割线的函数表达式恰为函数f(x)的以，为节点的线性插值多项式。用弦截法等价于用f(x)的线性插值多项式的零点近似函数f(x)的零点，因此又称作线性插值法。 用弦截法需要给出两个初始值，，通常取根所在区间的端点即可。 4.1.3 抛物线法（Muller法） 如果考虑用f(x)的二次插值多项式的零点来近似f(x)的零点，就导出抛物线法。设已知方程f(x)=0)的根的三个近似值，以这三点为节点的f(x)的二次插值多项式为： （4－1） 为方便起见，令： （4－2） 则（4-1）式可改写为： （4－3） 其零点为： （4－4） 按式(4-4)，f(x)的二次插值多项式P2(x)有两个零点，取哪个作为新的近似根？考虑到已是方程(4-1)的近似根，所的近似根自然应在的邻近，故选取新近似根的原则是使较小，于是有： （4－5） 按式(4-5)计算方程(4-1)的近似根称为抛物线法，也称Muller方法或二次插值法。 如图所示，y=P2(x)是过曲线上三点的抛物线．故抛物线法的几何意义是以过曲线上三点的抛物线与x轴的交点作为曲线与x轴交点的近似。 实际计算表明，抛物线法对初值要求并不苛刻。在初值不太好的情形下常常也能收敛。它的缺点是程序较复杂，并且在计算实根的过程中，也常常需要采用复数运算，增加／工作量。因此，抛物线法适用于当初值不太好时求方程的复根的情况。 4.2 非线性方程组的解法 设有非线性方程组： （4-6） 若记，则方程组（4-6）可简记为向量形式： 解非线性方程组的方法有许多种，其中一类是线性化的方法，它是将方程组中的每个方程线性化得到一个线性方程组，由此构造迭代格式，求得非线性方程组（4-6）的近似解。这类方法中常用的一种是Newton法；另一类方法是将解非线性方程组问题化成优化问题，然后以最优化方法求解。最速下降法是其中最基本的方法。 4.2.1 Newton-Raphson法 考虑一个由两个方程构成的非线性方程组： 考虑方程组在估计值（x1k，x2k）处的一阶Taylor展开式： 也可以将此方程组写成如下的矩阵-向量形式。 由此可以得到新的估计值： 若令，则方程求解过程可以写成： 这也可以看作4.1.1节中Newton法的直接推广。 考虑n维非线性方程组，设给定非线性方程组（4-6）的第k次近似解，以函数在处的一阶Tayler多项式近似函数，得非线性方程组： 即： （4-7） 其中：，称为向量函数的Jacobian矩阵，。线性方程组（4-7）称为Newton方程组。如果的Jacobi矩阵在处非奇异，方程组（4-7）有唯一解，则为非线性方程组的第k+1次近似解。此方法称为Newton法，也称作牛顿－拉夫森（Newton-Raphson）法。 计算过程的停止准则采用以下两种： ①