数学建模第二章 非线性方程求根.ppt

第二章 非线性方程求根 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(m)0 打印m, k a=m b=m 结束 k=K+1 是 是 否 否 输入 k = 0 例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 - §2 迭 代 法 1．简单迭代法 x1 = 0.4771 x2 = 0.3939 … x6 = 0.3758 x7 =0.3758 2．迭代过程的收敛性 f (x) = 0 x = φ (x) 等价变换 例3：求方程 的一个根 迭代格式 x y y = x x y y = x x y y = x x y y = x x* x* x* x* y=φ(x) y=φ(x) y=φ(x) y=φ(x) x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? 定理2：如果? (x)满足下列条件 （1）当x?[a, b]时，?(x)?[a, b] （2）当任意x?[a, b]，存在0 L 1，使 则方程x = ? (x)在[a, b]上有唯一的根x*,且对任意初值 x0?[a, b]时，迭代序列xk+1=? (xk) (k = 0, 1, …)收敛于x*。 (2.1) 3．迭代法的结束条件 例4：求方程 在[0, 0.5]内的根，精确到10-5。 x1 = ? (0.25) = 0.3385416 (2.2) x2 = ? (x1) = 0.3462668 x3 = ? (x2) =0.3471725 x4 = ? (x3) =0.3472814 x5 = ? (x4) =0.3472945 x6 = ? (x5) =0.3472961 x7 = ? (x6) =0.3472963 4．迭代过程的收敛速度 设由某方法确定的序列{