第章计算方法精要.pptVIP

* 第二章 一元非线性方程的解法 设非线性方程为 f (x)=0 (2-1) 方程（2-1）的解 x* 称为方程的根或函数f (x)的零点。若f (x)可表示为 f (x)=(x-x*)m g (x) 其中m为大于1的整数,且g(x) ≠0,称x*为方程(2-1)的m重根,或函数 f (x) 的m重零点. 若 f (x)为n次多项式,则称 f (x)=0为n次代数方程;若 f (x)为超越函数,则称f (x)=0为超越方程 。 若 f (x)在[a,b]内连续, 且 f(a)?f(b)0, 则 f(x)=0 在[a,b]内必有根; 若f(x)在[a,b]内还严格单调,则 f(x)=0在[a,b]内只有一根。据此可得求隔根区间的两种方法。 1. 做图法 画出 y= f(x) 的草图, 由 f (x) 与横轴交点的大概位置来确定隔根区间; 或者利用导函数 f (x)的正、负与函数 f (x) 的单调性的关系确定根的大概位置。 第一节 二分法 一、求隔根区间的一般方法 若 f (x) 比较复杂, 还可将方程 f(x)=0 化为一个等价方程 φ(x)=ψ(x), 则曲线 y=φ(x)与 y=ψ(x) 之交点 A(x*,y*)的横坐标 x* 即为原方程之根, 据此也可通过作图求得x*的隔根区间. 例1 判别下列方程有几个实根，并求隔根区间。(1) f(x)=x3-x-1=0 (2) f(x)=x4-4x3+1=0 解 将方程变形为x3=x+1, 绘曲线 y=x3 及 y=x+1，由图可知,方程只有一个实根x*?(1,1.5)，所以(1,1.5) 即为其隔根区间。 方程 f (x)=x4-4x3+1=0,由 f ?(x)= 4 x2(x-3)=0 得驻点x1=0, x2=3. 该二点将实轴分为三个区间: (-∞,0), (0,3)，(3, +∞) f ?(x) 在此三区间的符号分别为“-”、“- ”、“+”. 又知 f (-∞)0, f (0)=10, f (3)=-260, f (+∞)0 可见 f (x) 仅有两个实根, 分别位于(0, 3) , (3,+∞), 又 f (4)=10, 所以第二根的隔根区间可缩小为 (3,4)。 以上分析可用下表表示 (3,4) (0,3) 隔根区间 + ↗ ++ + ↗ 0 - - ↘ 0 + - ↘ f ?(x) f (x) (4,+∞) 4 (3,4) 3 (0,3) 0 (-∞,0) x 2. 逐步搜索法 从区间[a,b]的左端点 a 出发, 按选定的步长h 一步步向右搜索，若 f(a+jh) ·f(a+(j+1)h)0 (j=0,1,2, …) 则区间[ a+jh , a+(j+1)h ]内必有根。搜索过程也可从b开始，这时应取步长 h 0。 二、二分法 设f (x)在区间[a,b ]上连续, f (a) ·f (b)0, 则[a, b] 内有方程的根。取[ a , b ]的中点 将区间一分为二。若 f ( x0 ) = 0, 则 x0 就是方程的根,否则判别根 x*在 x0 的左侧还是右侧。 若f (a) ·f (x0) 0, 则x*∈( a , x0 ), 令 a1= a , b1=x0; 若f (x0) ·f (b) 0, 则x*∈(x0 , b ), 令 a1= x0 , b1=b. 不论出现哪种情况,( a1 , b1 )均为新的有根区间, 它的长度只有原有根区间长度的一半, 达到了压缩有根区间的目的。 对压缩了的有根区间, 又可实行同样的步骤, 再压缩。如此反复进行, 即可的一系列有根区间套 由于每一区间都是前一区间的一半，因此区间[ an , bn ]的长度为 若每次二分时所取区间中点都不是根，则上述过程将无限进行下去。当 n→∞ 时，区间必将最终收缩为一点x* ，显然x*就是所求的根。若 取区间[ an , bn ]的中点 作为x*的 近似值，则有下述误差估计式 只要 n 足够大, 即区间二分次数足够多 ），误差就可足够小。 由于在偶重根附近曲线 y=f(x) 为上凹或下凸，即 f(a) 与 f(b) 的符号相同，因此不能用二分法求偶重根. 例 2 用二分法求例1中（1）方程的实根,要求误差不超过0.005。 解 由例1可知x* ∈ (1,1