科学与计算工程办法第2章.pdfVIP

第二章 方程求根 §2.0 引 言 §2.1 二分法 §2.2 简单迭代法 §2.3 牛顿（Newton ）法 §2.4 其它求根方法（迭代过程的加速方法） §2.5 作业讲评 §2.0 引 言 非线性科学是当今科学发展的一个重要研究方向，非线性方程的求根也成为 其中一个重要内容。一般而言，非线性方程的求根非常复杂。 在实际应用中 有许多非线性方程的例子，例如 （1）在光的衍射理论（the theory of diffraction of light)中,需要求 x-tanx=0的根 （2）在行星轨道（ planetary orbits）的计算中,对任意的 a 和b，需要求 x-asinx=b 的根 （3）在数学中，需要求n 次多项式 a xn + a xn -1 +... +a x + a 0 的根。 n n -1 1 0 非线性方程的一般形式 f (x ) 0 这里 f (x ) 是单变量 x 的函数，它可以是代数多项式 f (x) a xn =+ a xn-1 +...+a x + a ( a ≠0 ) n n-1 1 0 n 也可以是超越函数，即不能表示为上述形式的函数。满足方程 f (x ) 0 的 x 值通常叫做方程的根或解，也叫函数 f (x ) 0 的零点。 §2.1 二分法(Bisection Method) 1 概念: 二分法也称对分区间法、对分法等，是最简单的求根方法，属于区 间法求根类型。 在用近似方法时，需要知道方程的根所在区间。 若区间[a,b]含有方程f (x)=0 的根，则称[a,b]为f (x)=0 的有根区间； 若区间[a,b]仅含方程f (x)= 0 的一个根，则称[a,b]为f (x)= 0 的一个 单根区间。 2 ．基本思想 根的存在定理(零点定理) ： f(x)为[a,b]上的连续函数，若f (a)·f (b)0 ，则[a,b] 中至少有一个实 根。如果f (x)在[a,b]上还是单调递增或递减的，则f (x)=0 仅有一个 实根。 x1 x* a x2 b 图示：二分法 3 ．构造原理 直接取区间[a,b]的中点x=(a +b)/2 作为问题的近似解，那么我 们可以估计出绝对误差限仅为区间长的一半，即 e=(b-a)/2 。如果 这个结果能满足精度要求，我们就停止进一步的计算；如果不能，就 求出 f(x),结果只能是下面三种情况之一： （1） f(a)·f(x)0 ，此时我们有x* ∈[a,x]； （2 ） f(x)·f(b)0 ，此时我们有x* ∈[x,b]； （3 ） f(x)=0 ，此时x 即为问题的精确解。 在前两种情况下，我们可以用 x 分别替换原问题中的 b 或 a,从而把 求解的区间减小了一半。这样我们又可以取新区间[a,b]的中点。经过N 次迭代后，剩下的区间长为(b- a