第3章 一元非线性方程的解法1.ppt

第3章 方程的近似解法 §3.1 二分法 §3.2 迭代法 §3.3 切线法(牛顿法) §3.4 弦截法 §3.5 加速迭代法  §3.1 根的搜索与二分法 我们已经熟悉求解一元一次方程、一元二次方程以及某些特殊类型的高次代数方程或非线性方程的方法。这些方法都是代数解法,求出的根是方程的准确根。但是在许多实际问题中遇到的方程,例如代数方程 x3-x-1=0 ? 或超越方程 等等,看上去形式简单,但却不易求其准确根。为此,只能求方程达到一定精度的近似根。 方程的形式很多,我们主要讨论一元线性方程,也即 f(x)=0 (3―1) 方程(3―1)可以有实根,也可以有复根或者重根等。本章主要讨论它的实根的数值计算问题。 方程根的数值计算大致可分三个步骤进行: (1) 判定根的存在性。 (2)确定根的分布范围,即将每一个根用区间隔离开来。 (3)根的精确化,即根据根的初始近似值按某种方法逐步精确化,直至满足预先要求的精度为止。 3.1.1 根的搜索 设f(x)为定义在某[a,b]上的连续函数,有f(a).f(b)0,则方程f(x)=0 在[a,b]内至少有一个实根.这时我们称[a,b]为方程f(x)=0的有根区间。 设f(x)为定义在某区间上的连续函数,方程(3―1)存在实根。虽然方程(3―1)的根的分布范围一般比较复杂,但我们不难将函数f(x)的定义域分成若干个只含一个实根的区间。 例如考虑方程 x2-2x-1=0 由图2.1所示,该方程的一个负实根在-1和0之间,另一个正实根在2和3之间。 这样,我们总可以假设方程(3―1)[a,b]内有且仅有一个单实根x*。由连续函数的介值定理知 f(a)·f(b)＜0 从x0=a出发，按步长h=(b-a)/N (N为正整数) ，逐点计算 Xk=a+kh处的函数值f(xk),当f(xk)与f(xk+1)的值异号时，[xk,xk+1] 就是 方程f(x)=0的有根区间。 若数值b-a较小,那么我们可在[a,b]上任取一点x0作为方程的初始近似根。 f(x)=xn+a1xn-1+a2xn-2+…+an-1x+an=0的根的绝对值上下界估计 (1)若aa=max{|a1|,|a2|,…,|an|} 则 |x|aa+1 (2)若bb=(1/|an|) .max{1,|a1|,|a2|,…,|an-1|}, 则 |x| 3.1.2 二分法 下面在有根区间[a,b]内介绍二分法的基本思想。 取x0=(a+b)/2 计算f(a)与f(x0),若 ? f(a)·f(x0)＜0 ? 则根x∈(a,x0),令 a1=a,b1=x0 否则x∈(x0,b),令 ? a1=x0,b1=b 如此逐次往复下去,便得到一系列有根区间 [a,b],[a1,b1],[a2,b2],…,[ak,bk],… 其中 我们把每次二分后的有根区间[ak,bk]的中点 ? 则结果xk就是方程(3―1)满足预给精度ε的近似根,也即 1.计算步骤 ①输入有根区间的端点a,b及预先给定的精度ε; ②(a+b)/2 x; ③若f((a+b)/2)=0,则(a+b)/2为f(x)=0的根。 ④若f(a)f(x)＜0,则x→b,转向④;否则x→a,转向⑤ 。 ⑤若|b-a|＜ε,则输出方程满足精度的根(a+b)/2,结束;否则转向②。 例1 求方程