09-10方程求根.docVIP

第十章 非线性方程求根 知识点：根的隔离,根的精化,二分法，一般迭代法，牛顿迭代法，收敛性。 1.根的性态 （1）根的概念 方程f(x)=0的解叫做方程的根或f(x)的零点，当f(x)较复杂时，很难求解析解，但可以求近似根。求实根问题包括:根的存在性;根的分布;根的精确化. （2）根的存在性 理论:如果f(x)∈C[a,b]、f(a)f(b)0，则f(x)=0在（a,b)至少有一个根。 若f(x)∈C[a,b]、单调、f(a)f(b)0 ，则f(x)=0在(a,b)有惟一根。 （3）根的分布 ①(多项式)方程实根的上下界 ②根的隔离 在有根区间确定若干个子区间，使每个子区间有且只有f(x)=0的一个根的过程称根的隔离。一般情形下，隔根区间内任一点都可作为相应根的近似值，隔根区间越小，近似程度越好。 例 找 的根的位置。 解 因为 y=f(x)是单调递增函数。从f(0)=-7，f(1)=-5，f(2)=1知f(x)=0在（1，2）有且仅有一根（作图从略）。 （4）根的精确化 求出f(x)=0根的一个近似值，然后将根的近似值逐步精确化，直到满足精度要求为止。 2.二分法 （1）基本思想 对分有根区间，判断f(x)的符号，逐步将有根区间缩小，使得在足够小的区间内（方程有且仅有一根）取得满足预先给定精度度要求的近似值。 （2）具体做法 取[a,b]为隔根区间，二等分区间[a,b]。令x0=(a+b)/2，如果f(x0)=0,则求得实根α=(a+b)/2。否则，若f(a)f(x0)0,知α∈[a,x0],取a1=a,b1=x。；若f(b)f(x0)0,知α∈[x0,b],取a1=x0,b1=b。 在新的有根区间上重复二等分过程，得有根区间序列 [a,b] [a1,b1] [a2,b2] … [ak,bk] … 其中每个区间为前一个区间的一半，经过k次二分后，有根区间的长度为 无限二等细分区间[a,b](k→∞), 满足ak≤α≤bk，bk-ak→0(k→∞)，limak=limbk=α(k→∞)。 取有根[ak,bk]的中点xk=(ak+bk)/2为α的近似值，则以α为极限的近似值序列为 x0,x1,x2,…,xk,… 由于 对于预先给定精度ε0，只要 便有 这时的xk就是满足精度要求的近似值。 例 用二分法求 在(1,2)内的根，要求绝对误差不超过 解： 函数单调，在（1，2）内有唯一根，取为隔根区间。 f(1)=-50 有根区间 中点 xn f(2)=140 (1,2) x1=1.5 f(1.5)0 (1,1.5) x2=1.25 f(1.25)0 (1.25,1.5) x3=1.375 f(1.375)0 (1.25,1.375) x4=1.313 f(1.313)0 (1.313,1.375) x5=1.344 f(1.344)0 (1.344,1.375) x6=1.360 f(1.360)0 (1.360,1.375) x7=1.368 f(1.368)0 (1.360,1.368) x8=1.364 若取近似根α≈x8=1.364，则 （3）二分法求根算法 算法：设f(x)∈C[a,b]、f(a)f(b)0，计算f(x)=0在（a,b)内的根。 1:输入有根区间[a,b]的a、b值及精度控制参数ε1、ε2（两参数可以统一唯一个）。 2:if f(a)f(b)0 then 退出选用其他求根方法 else 转步骤 3 3：while |a-b|ε2时 令x=(a+b)/2，计算f(x)； if |f(x)| ε1 then 停止计算 x*=x，转步骤 5 else if f(a)f(x)0 then 修正区间[a,x]→[a,b] else 修正区间[x,b]→[a,b] endwhile 4:x*=(a+b)/2 5:输出近似值x*（ε1、ε2两参数统一为一个算法更简单） 优点：条件简单，编程容易。缺点：收敛慢，只能求实根，不能求偶数重根和复根。 3.迭代法 代法的基本思想 设法将方程f(x)=0改写为等价形式：x=??? x)（f和??连续），xn+1=??? x n)，n=0,12,…. 在根邻近任取初始值x0作为近似根，由迭代公式得一序列 x0，x1，x2，…