数值分析第二章new.pptVIP

简介（Introduction） 我们知道在实际应用中有许多非线性方程的例子，例如 （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的根 第二章 非线性方程求根 1．根的存在性。方程有没有根？如果有根，有几个根？ 定理1：设函数 f (x) 在区间[a, b]上连续，如果f (a) ? f (b) 0， 则方程 f (x) = 0 在[a, b]内至少有一实根x*。 2．这些根大致在哪里？如何把根隔离开来？ 3．根的精确化 a b x* f(x) 1．画出 f(x) 的略图，从而看出曲线与x 轴交点的位置。 2．从左端点x = a出发，按某个预先选定的步长h 一步一步地向右跨，每跨一步都检验每步起点x0 和终点x0 + h的函数值，若 那么所求的根x*必在x0与x0+h之间，这里可取x0或x0+h 作为根的初始近似。 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 开 始 读入a, h a ? x0 f (x0) ? y0 x0 + h ?x0 f (x0) ? y00 打 印 结 束 否 是 继续扫描 ? 例1：考察方程 x 0 0.5 1.0 1.5 f (x) 的符号 － － － ＋ a b x1 x2 a b 或 不能保证 x 的精度 x* ?2 x x* §1 二 分 法 执行步骤 1．计算f (x)在有解区间[a, b]端点处的值，f (a)，f (b)。 2．计算f (x)在区间中点处的值f (x1)。 3．判断若f (x1) = 0，则x1即是根，否则检验： （1）若f (x1)与f (a)异号，则知解位于区间[a, x1]， b1=x1, a1=a; （2）若f (x1)与f (a)同号，则知解位于区间[x1, b]， a1=x1, b1=b。 反复执行步骤2、3，便可得到一系列有根区间： (a, b), (a1, b1), …, (ak, bk), … 4、当 时 5、则 即为根的近似 ①简单; ② 对f (x) 要求不高(只要连续即可) . ①无法求复根及偶重根 ② 收敛慢 定义f (x) f (a)? f (b)0 f (a)? f (b)=0 f (a) =0 打印b, k 打印a, k 结束 是 是 是 否 否 否 m=(a+b)/2 |a-b|? f(a)?f(b)0 打印m, k a=m b=m 结束 k=K+1 是 是 否 否 输入 k = 0 误差 分析： 第1步产生的 有误差 第 k 步产生的 xk 有误差 对于给定的精度 ? ,可估计二分法所需的步数 k ： 例2： 求方程 k ak bk xk f (xk)的符号 0 1 1.5 1.25 - 1 1.25 1.5 1.375 + 2 1.25 1.375 1.3125 - 3 1.3125 1.375 1.3438 + 4 1.3125 1.3438 1.3281 + 5 1.3125 1.3281 1.3203 - 6 1.3203 1.3281 1.3242 - f (x) = 0 x = g (x) 等价变换 f (x) 的根 g (x) 的不动点 f(x)=0化为等价方程x=g(x)的方式是不唯一的,有的收敛,有的发散 For example：2x3-x-1=0 §2 迭 代 法 1．简单迭代法 (1) 如果将原方程化为等价方程 由此可见，这种迭代格式是发散的 取初值 (2) 如果将原方程化为等价方程 仍取初值 依此类推,得 x3 = 0.9940 x4 = 0.9990 x5 = 0.9998 x6 = 1.0000 x7 = 1.0000 已经收敛,故原方程的解为 x = 1.0000 同样的方程 ? 不同的迭代格式 有不同的结果 什么形式的迭代法能够收敛呢? x1 = 0.4771 x2 = 0.3939 … x6 = 0.3758 x7 =0.3758 例3：求方程 的一个根 迭代格式