《计算机在材料科学与工程中的应用ch04 非线性方程的解法.pptVIP

第四章非线性方程的解法 主要内容 1、引言 引言 哪些情况下，非线性方程通常用计算机来求解？ 求根步骤 1.确定根的初始近似值。（粗略判断有根范围） 2.根的精确化（几种方法：二分法、迭代法、牛顿迭代法、弦截法等等） 二分法 基本思想 迭代法 基本思想 1,把一般方程f(x)=0化成等价形式x=g(x)，构造迭代式xk+1=g(xk)，k=0,1,2,… 2,由近似初始值x0，得到一系列迭代值：x0,x1,x2,…,xk,…。——记为{xk} 牛顿迭代法和弦截法 牛顿迭代法 §基本思想 本章小结 1/4、引言 * * Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 2、二分法 3、迭代法 (1),问题的提出 (2),求解步骤 计算机处理非线性方程的基本方法 常用的计算方法 4、牛顿迭代法和弦截法 内容3的两个例子 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. ——高阶代数方程 ——超越方程 不易用解析的方法求解，通常可用图解法获得近似解，在精确的场合则需要用数值方法 一般情况下，上述方程可写成f(x)=0的形式 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 第一步骤——粗略判断有根范围 (1) 判断f(x)在[a，b]连续，且f(a)*f(b)0 (2) f(x)在[a,b]单调，则仅有一个实根。 (3) 把[a,b]分成长度为h的小区间，x0=a, xk=a+k*h,计算f(xk). (4) 如果f(xk)f(xk+1)0，则判断根 x*在[xk,xk+1]内。 或者，如果f(x0)f(xk)0，且f(x0)f(xk+1)0 ，则判断根x* 在[xk,xk+1]内。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例： 方程 f(x)=x3-x-1=0 ，找出有根范围，要求范围大小为h=0.5 解：(a)先判断 (b)从x=0出发，h=0.5，判断f(x)的符号 + - - - f(x) 1.5 1 0.5 0 x (c)则粗略判断根 x* 在 [1,1.5] 内 (d)再利用其他方法进行精确化 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. a,平分含根的区间，判断根的位置在哪个区间 b,舍去无根的区间，再进行a的判断 c,重复a，b过程，直到达到事先要求的精度为止 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 计算过程 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例： 方程 f(x)=x3-2x2-4x-7=0在[3,4]内的根，精确到10-4 Y N N Y x=(a+b)/2 f(x)0? a=x N a,b,eps 开始 f(a) 结束 b=x |b-a|eps? 输出x x=(a+b)/2;f=f(x) f0*f0? a=x a,b,eps 开始 f0=f(a) float a,b,eps,x,f0,f; a=3;b=4;eps=0.0001; f0=f(a); while(1) { x=(a+b)/2; f=f(x); if(f0*f0) b=x; else a=x; if(fabs(b-a)eps) { printf(x=%f,a);