第2章非线形方程及其非线性方程组解法解析.ppt

* * * §4 联立方程组的迭代解法 充分条件2: 若 ,则迭代过程收敛,且 列求和的最大值 * §4 联立方程组的迭代解法 若 ,则迭代过程收敛,且 充分条件3: 各元素平方和的开方。 * 通过画图法可以知道其中一个根的初值可取(x0,y0)=(3.4,2.2). 方程组可以改写为: §4 联立方程组的迭代解法 例2.11： 求下列方程组的根 解: * §4 联立方程组的迭代解法 则 * v1=a11+a21?2.3 v2=a12+a22?4.4 v=max(v1, v2) ?4.41 因为v 1，可能不收敛。方程组重新改写为: §4 联立方程组的迭代解法 按照充分条件2判断： * §4 联立方程组的迭代解法 则 * 计算结果为x*=3.487, y*=2.262。 则 v1=a11+a21?0.834 v2=a12+a22?0.248 v=max(v1, v2) ?0.8341 符合充分条件2，按照下式一定收敛 §4 联立方程组的迭代解法 按照充分条件2判断： * §5 联立方程组的延拓解法 5.1 同伦方程组及其建立方法 （扩大收敛域的求根方法） 同伦方程组 H(X,t)=0 （1）H(X(0),0)=0， （2） H(X,1)= F(X)=0, 满足上述条件的方程组为同伦方程组. X(0)=(x1(0), x2(0), …,xn(0))是任意取定的一组初值,即当t=0时,H(X,0)是具有已知解X(0)的联立方程组。 当t=1时,H(X,1)=0与原联立方程组F(X)相同。 * 5.1 同伦方程组及其建立方法 同伦方程组的构造 或 * §5.2 求解方法 (1)将t的值域[0,1]等距或不等距的划分为(为将已知解X(0)逐步引渡到F(X)=0的未知解) 0t0 t1…tN-1 tN=1 代入同伦方程组得 Hi=H(X, ti)=0 (i=1,2,…,N) (2)以下取X(0)作为H1=0的初值，用某种迭代法求解出H1=0的解X(1)。 * (4)当上述同伦方程组的解为t的连续时，则X(N)就是原方程F(X)=0的解 (3)再以X(1)作为H2=0的初值，求解出H2=0的解X(2)，继续下去，直到求出HN=0的解X(N)为止。 §5.2 求解方法 * §5.2 求解方法 例2.11： 求下列方程组的根 解: * §5.2 求解方法 例2.11： 求下列方程组的根 解: * §5.2 求解方法 * §5.2 求解方法 * §5.2 求解方法 * §5.2 求解方法 * §6 联立方程组的牛顿解法 * §6 联立方程组的牛顿解法 例子：试用牛顿迭代法求解下列方程组的根。 * §6 联立方程组的牛顿解法 * §6 联立方程组的牛顿解法 * * * * * * * * * * * * §3.3.2 牛顿迭代法的收敛性 x y 0 x0 x y 0 x0 不满足迭代条件时，可能导致迭代值远离根的情况而找不到根或死循环的情况。 （1）（2）（3）不满足。 （3）（4）不满足。 * §3.3.3 切线（牛顿）法的变形使用 取f ’(xn)为固定值，例如取f ’(x0) ，这时迭代公式成为xn+1=xn- f(xn)/ f ’(x0),称为简化切线法，或固定斜率的切线法 1. 简化切线法 y x 0 x0 y =f(x) x1 x2 ? 几何意义 * §3.3.3 切线法的变形使用 推广的简化切线法 当取f ’(xn)=C(任意不为0的常数),迭代公式成为xn+1=xn- f(xn)/ C,称为推广的简化切线法. C值应满足|φ’(x)|=|1- f ’(x)/ C |1,解不等式有: 1- f’(x)/ C 0 1- f’(x)/ C 1 -(1- f’(x)/ C) 0 -(1- f’(x)/ C) 1 解得0f ’(x)/ C2 可见,当C与f ’(x)同号,且满足上式时,推广的简化切线法是收敛的。 * §3.3.3 切线法的变形使用 引进参数λ，xn+1=xn- λ f(xn)/ f ’(xn)，并用尝试法修改λ值大小(即改变原切线的斜率)，使达到|f(x0)| |f(x1)|…单调下降的目的。称|f(xn+1)| |f(xn)|为下山条件,这种算法为下山法. 2.牛顿下山法 下山条件: * §3.3.3 切线法的变形使用 具体步骤 λ=1，即在xn基础上用切线法迭代一次得xn+1，然后检查下山条件是否满足； 若下山条件满足,则继续用切线法迭代,当满足精度要求时,停止迭代,获得最终