方程求根-read.docVIP

第二章　　方程求根 2.1　方程求根与二分法2.1.1　引言单个变量的方程(2.1.1) 求根是数值计算经常遇到的问题.当f(x)为一般连续函数时，称式(2.1.1)为超越方程，如果f为多项式若,为n次多项式，此时方程(2.1.1)称为代数(或多项式)方程.如果x*(实数或复数)使，则称x*为方程(2.1.1)的根，若，m≥1，且，当m＞1时，称x*为方程(2.1.1)的m重根或称x*是f的m重零点.若x*是f的m重零点，且g充分光滑，则。当f为式(2.1.2)表示的代数多项式时，根据代数基本定理可知方程(2.1.1)有n个根(含复根，m重根为m个根)，对n=2的代数方程它的根可由公式表示为而当n=3,4时方程的根虽可用公式表示，但表达式太复杂，一般不用，当n≥5已没有直接用公式表达的求根算法.因此对n≥3的代数方程求根方法与一般超越方程(2.1.1)一样都采用迭代方法求根，设(表示f在区间上连续)，若有f()f(b)＜0,则f(x)=0在区间上至少有一个实根，称为有根区间，通常可用逐次搜索法求得方程(2.1.1)的有根区间.例2.1　求方程的有根区间.解　根据有根区间定义，对方程的根进行搜索计算，结果如下表： 方程的三个有根区间为［1,2］,［3,4］,［5,6］. 讲解： 　　非线性方程f(x)=0求根，包括求超越方程和代数方程的根x*,方程的根也是f(x)的零点，即f(x*)=0,x*可以是实根也可以是复根，本章以求实根为主。求实根首先要确定根x*所在区间，称为有根区间。根据连续函数性质，若f(x)在上连续，当f()f(b)0时，为有根区间,为找到方程f(x)=0的有根区间,可用逐次搜索法,也就是在x的不同点上计算f(x),观察f(x)的符号,如例2.1表中所示,只要在相邻两点f反号,则得到有根区间,本例得到三个有根区间,分别为[1,2][3,4][5,6]. 2.1.2 二分法 　　设,且为有根区间，取中点，将它分为两半，检查与是否同号，若是，说明根x*仍在右侧，取,否则取,得到新的有根区间长度仅为的一半(见图2-1).重复以上过程，即取，将再分半，确定根在 的哪一侧，得到新区间，其长度为的一半，从而可得一系列有根区间 其中每一个区间长度都是前一个区间长度的一半，因此，的长度为 且,即为方程(2.1.1)的根x*的一个足够精确的近似根，且误差为 　　(2.1.3) 以上过程称为解方程的二分法.它计算简单且收敛性有保证，但收敛较慢，通常可用于求迭代法的一个足够好的初始近似. 例2.2 求方程在区间［1.0,1.5］内的一个实根，要求准确到小数点后第二位. 　　解 这里=1.0,b=1.5，而f()＜0,f(b)＞0，故在［1.0,1.5］中有根，由式(2.1.3)知 　　,即,当n=6时得,已达到精度要求，各次计算结果见表2-1. 故为方程的近似根，误差不超过0.005. 讲解： 　　从有根区间出发用二分法将它逐次分半，得到的近似根 ,其误差为(2.1.3)式．它是求实根近似的有效方法，但由于收敛慢，只用作迭代法的初始近似． 2.2 迭代法及其收敛性 2.2.1 不动点迭代法 求方程(2.1.1)的根时将方程改写为等价形式 　　　　　　　(2.2.1) 若x*满足，则称x*为φ的一个不动点，x*也是方程(2.1.1)的一个根，如果φ连续，可构造迭代法 　　 　　　(2.2.2) 称为不动点迭代法，φ称为迭代函数.若给定初始近似，由式(2.2.2)逐次迭代得到序列，如果，则由式(2.2.2)两端取极限，得 即x*为φ的不动点.? 从几何图象(参见图2-2)可知，（参见图2-2）可知，φ的不动点就是直线与曲线的交点，横坐标x*即为不动点，从它的某个初始近似出发，在曲线确定一点，引平行x轴直线，与直线交于点，其横坐标即为,由式(2.2.2)逐次求得即为如图2-2所示点的横坐标. 例2.3 求方程在附近的根.? 解 若将方程改写为,构造迭代法 　　　　　　(2.2.3) 由可知，显然不收敛.? 若将方程改为，构造迭代法 　　　　　　　(2.2.4) 计算结果见表2-2. 表 2－2 从结果看它是收敛的，且在6位有效数字时即为根x*的近似值.? 　　例题表明构造迭代法(2.2.2)，必须使迭代序列｛｝收敛，才能求得方程(2.2.1)的解x*.下面给出解的存在唯一性和迭代收敛性定理. 定理2.1 设，如果对有≤φ(x)≤b,且存在常数L(0,1)使 　　　　　　　　(2.2.5) 则φ在区间上存在唯一不动点x*，且由式(2.2.2)生成的迭代序列｛｝对任何收敛于x*，并有误差估计 　　　　　　　　　(2.2.6) 证明 先证明不动点存在性，记,由定理条件有及,若有一等号成立，则或,即φ有不动点，否则必有，因故必有,使,x*即为