第四章 非线性代数方程与方程组求解.pptVIP

第四章 非线性方程和非性方程组的解法; 　　１.非线性方程的解法 　　２.非线性方程组的线性化解法 　　　－－牛顿迭代法 　　３.非线性方程组的极值求解法 　　　－－最速下降法　|　单纯形法 　　　－－共轭梯度法　|　Powell 方法 　　　－－变尺度法 | 　　　（可变矩阵方法）|　直接法 　　　　DFP 方法 |; ４.1 引言 　　在科学研究中，常常会遇到非线性方程 或非线性方程组的问题。例如解方程 　或 　一般的，我们记非线性方程为;4.1 　　 非线性方程组的一般形式是：;4.1;4.1 ? 方程的解亦称方程的根或函数的零点。 ? 　根可能是实数或复数。 ? 　若　　　　　　　　　　则　称为单根； 　　若　 　　而　　　　　，则　称为 k 重根。 　 常见的求解问题有两种： 　(1) 要求定出在给定范围内的某个解。 　(2) 要???定出在给定范围内的全部解。 ? 非线性问题，除少数情况外，一般不能 　　不利用公式求解。而要采用某种迭代解法。 　　即构造出一近似值序列　　　　　　　　 　　逼近真解 。;4.1 ? 　 迭代过程的收敛性一般与初值的选取和方 　　程的性态有关，某些解法仅与初值有关。 ? 收敛速度一般由迭代方法所决定，方程的 　　性态也会起一些作用。 　　 本章主要介绍非线性方程组的解法， 而方程的解法用较少的篇幅在4.2中扼要介绍。 　　解非线性方程和方程组有很大区别。后者 　　要困难得多。主要的区别在于一维情形可 　　以找到一个根的范围，然后缩小，最终找 　　到根。而多维情况则很难确定根的存在。 　　直到你求得它的解。;4.2 非线性方程的解法;4.2.1 用图来表示这个过程：;a;a; 　　例如，在［a,b］之间寻找 f(x) 可能有 　　的根可以用等分试探法：; 用二分法求函数FUNC位于（x1,x2）之间的根，准确性为?XACC。 FUNCTION RTBIS (FUNC,X1,X2,XACC) PARAMETER (JMAX=40) FMID=FUNC(X2) F=FUNC(X1) IF (F*FMID.GE.0.) PAUSE 函数FUNC在x1,x2处不异号 IF (F.LT.0.) THEN RTIBIS=X1 DX=X2-X1 ELSE RTBIS=X2 DX=X1-X2 ENDIF DO 11 J=1, JMAX DX=DX*0.5 XMID=RTBIS+DX FMID=FUNC(XMID) IF (FMID.LE.0.) RTBIS=XMID IF (ABS(DX) .LT. XACC .OR. FMID .EQ. 0.) RETURN 11 CONTINUE PAUSE 迭代次数越界 END;FUNCTION FF(X) FF = X*X + 2.5*X + 0.5+SIN(X) END PROGRAM ROOTFIND EXTERNAL FF X1=-1.0 X2=0.0 ROOT=RTBIS(FF, X1, X2, 1.0E-5) PRINT *, 方程在（-1，0）区间内有一个根，X=, ROOT STOP END; 4.2.2 线性插值法（又称弦位法）;4.2.2;4.2.2;4.2.2;4.2.2; 4.2.3 Newton 法;4.2.3; 　　　　每次函数求值相当的收敛阶为： 　　　　　　 　　b. 求 fk 有时工作量大，甚至不可能。 　(4) 选用收敛域较大的方法（如二分法） 　　 先进行迭代，然后再用Newton法。 　　　　－－组合方法。 4.2.4 二次插值法 　　　　设 f(x)=0 的三个近似解及函数值 　　 构造二次函数g(y)使得：;4.2.4;4.2.4;4.2.4 (1) 要有三个初始值 (2) 当 。且收敛速度 是 1.84 阶。（单根） 二重根的收敛阶是1.23。 (3) (4) 发生超射、越界。 ; 4.2.5 组合方法（Brent Method） 　　　　能否有一种方法综合上述方法的优 　　点呢？Brent 做了一些工作。 　　　　Brent 把二分法和二次插值法结合 　　起来。 　　（１）一定收敛。 　　（２）收敛