非线性方程的求根方法1.pptVIP

Solutions of Nonlinear Equations 非线性方程的求根方法 Bisection Method方程求根的二分法 Fixed-Point Iteration迭代法及其收敛性 Newton Method of Nonlinear Equations Newton迭代法 在实际应用中有许多非线性方程的例子，例如 （1）在光的衍射理论（the theory of diffraction of light)中,需要求x-tanx=0的根 （2）在行星轨道（ planetary orbits）的计算中,对任意的a和b，需要求x-asinx=b的根 （3） 在数学中，需要求n次多项式xn + a1 xn-1+...+an-1 x + an ＝0的根 方程求解是科学计算中一个重要的研究对象; 几百年前就已经找到了代数方程中二次至五次方程的求解公式; 但是,对于更高次数的代数方程目前仍无有效的精确解法; 对于无规律的非代数方程的求解也无精确解法; 因此,研究非线性方程的数值解法成为必然。 非线性方程的一般形式： f(x)=0 代数方程： f(x)=a0+a1x+……+anxn （an?0） 超越方程 ：f(x)中含三角函数、指数函数、或其他超越函数。 用数值方法求解非线性方程的步骤： （1）找出有根区间；（只含一个实根的区间称隔根区间） （2）近似根的精确化。从隔根区间内的一个或多个点出发，逐次逼近，寻求满足精度的根的近似值。 Bisection Method 二分法 二分法的基本思想： 假定f(x)=0在[a,b]内有唯一单实根x*，考察有根区间[a,b]，取中点x0=(a+b)/2，若f(x0)=0，则x*= x0 ，否则， （1）若f(x0)f(a)0，则x*在x0右侧，令a1=x0， b1=b; （2）若f(x0)f(a)0，则x*在x0左侧，令a1=a， b1= x0。 介值定理 设函数f(x)在区间[a,b]连续，且f(a)f(b)0，则方程f(x)=0在区间(a,b)内至少有一个根。 以[a1, b1]为新的隔根区间，且仅为[a,b]的一半，对[a1, b1]重复前过程，得新的隔根区间[a2, b2]，如此二分下去，得一系列隔根区间： [a,b] ?[a1, b1] ? [a2, b2] ?…… ? [ak, bk] ? …… 其中每个区间都是前一区间的一半，故[ak, bk] 的长度： 当k趋于无穷时趋于0。即若二分过程无限继续下去，这些区间最后必收敛于一点x*，即方程的根。 性质： f(an)·f(bn)0; bn – an = (b – a)/ 2n Bisection Method 每次二分后，取有根区间的中点xk= (ak+bk) /2作为根的近似值，则可得一近似根序列： x0, x1, x2, …该序列必以根x*为极限。 实际计算中，若给定充分小的正数?0和允许误差限?1，当|f(xn)| ?0或bn- an ?1时，均可取x*? xn。 定理 设x*为方程f(x)=0在[a,b]内唯一根，且f(x)满足f(a)f(b)0，则由二分法产生的第n个区间[an, bn] 的中点xn满足不等式 证明： 1.先验误差估计： 利用误差估计定理，令 得 从而得到对分次数k，取xk作为根得近似值x*。 2.后验误差估计： 给定ε，每步检查 ，若成立，则取 ，否则继续对分。 例 用二分法求 在(1,2)内的根，要求绝对误差不超过 解： f(1)=-50 有根区间 中点 f(2)=140 -(1,2)+ f(1.25)0 (1.25,1.5) f(1.375)0 (1.25,1.375) f(1.313)0 (1.313,1.375) f(1.344)0 (1.344,1.375) f(1.360)0 (1.360,1.375) f(1.368)0 (1.360,1.368) f(1.5)0 (1,1.5) 计算过程简单，收敛性可保证； 对函数的性质要求低，只要连续即可。 收敛速度慢； 不能求复根和重根； 调用一次求解一个[a,b