第一章节常用数值分析方法§1非线性方程求根.pptVIP

第一章 材料科学研究中的常用数值分析方法 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 主 要 内 容 §1 非线性方程求解 §2 线性方程组的数值解法 §3 插值法与曲线拟合 §4 有限差分法与有限单元法 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. §1 非线性方程求解 1.1 概 述 1.2 对分法 1.3 迭代法? 1.4 Newton法 1.5 弦截法 其他方法： Aitken加速法 、Steffensen加速法、重根加速收敛法、 抛物线法、牛顿下山法、劈因子法等。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 1.1 非线性方程求解概述 很多科学计算问题常常归结为求解方程： Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 非线性方程求解概述(续） 例如，从曲线y = ?x和y = lg x的简单草图可看出方程 lg x + x = 0有唯一的正根x*，但是没有求x*的准确值的已知方法，即使是对代数方程，要求其精确解也是困难的。对于二次方程ax2 + bx+c = 0，我们可以用熟悉的求根公式： 对于三、四次代数方程，尽管存在求解公式，但并不实用。而对于大于等于五次的代数方程，它的根不能用方程系数的解析式表示，至于一般的超越方程，更没有求根公式。因此，为求解一个非线性方程，我们必须依靠某种数值方法来求其近似解。 对于方程（1-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)在区间[a,b]上连续，严格单调，且f (a ) f (b) 0，则在[a,b]内方程f (x) = 0有且仅有一个实根。 根据此结论，我们可以采用如下两种方法求出根的隔离区间。 1. 根的存在性：方程是否有根？如果有根，有几个根？ 2. 根的隔离：确定根所在的区间，使方程在这个小区间内有且仅有一个根，这一过程称为根的隔离，完成根的隔离，就可得到方程的各个根的近似值。 3. 根的精确化：已知一个根的粗略近似值后，建立计算方法将近似解逐步精确化，直到满足给定精度为止。 关于根的存在性是纯数学问题，不详细介绍，可查阅有关代数学内容。 根的隔离主要依据如下结论： Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 求根的隔离区间的两种方法 1. 描图法： 画出y = f (x)的草图，由f (x)与x轴交点的大概位置来确定有根区间。也可利用导函数f ?(x)的正、负与函数f (x)的单调性的关系来确定根的大概位置。 例1 求f (x) = 3x ? 1? cos x = 0的有根区间 解：将方程变形为3x ? 1= cos x 绘出曲线 y =3x?1及 y = cos x， 由图1-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)=4x3?12x2 = 0，可得驻点x1 = 0, x2 = 3,由此而得到三个区间 (??,0) (0,3),(3,??), f ?(x)在此三个区间上的正负号分别为“?”，“?”，“+”，由此可见，函数f