第五讲方程求根的数值解法.pptVIP

第五讲方程求根的数值解法 第五讲主要知识点 二分法； 二分法的收敛速度； 误差公式； 迭代法的思想； 迭代法的收敛性； 迭代过程的局部收敛性； 迭代法的收敛速度. 求实根近似值的常用方法 1、二分法 2、迭代法 3、牛顿法 4、弦截法 根的存在性 二分法 二分法步骤 二分法步骤（续） 二分法误差估计 误差估计、对分次数 二分法程序框图 例题分析 例题分析（续1） 例题分析（续2） 二分法总结 迭代法 迭代法（续1） 迭代法（续2） 迭代法（续3） 迭代法的结束条件 迭代法例题 迭代过程的收敛性 定理3.2 设迭代函数 在 上具有连续的一阶导数，且 （1）当 时，有 ； （2）存在正数 ，使对任意的 　，有 成立。则 在 上存在唯一的解 ，并且对任意选取的初值 ，迭代过程 所产生的迭代数列收敛于 。 迭代过程的收敛性证明 迭代过程的收敛性证明（续1） 迭代过程的收敛性证明（续2） 迭代过程的收敛性 定理：设函数 例题分析 迭代法总结 迭代过程的局部收敛性 收敛速度的概念 迭代过程的收敛速度 迭代过程的收敛速度（续1） 迭代过程的收敛速度（续2） 迭代过程的收敛速度（续3） 本讲结束！ 谢谢大家！ 再见！ 定理 若 在 的某个邻近内有 阶连续导数，且 则对一个任意靠近 的初始值，迭代公式 是 阶收敛的，且有 如果 ，要求 证明：由定理3.5知，迭代格式（3.19）是局部收敛 的．应用泰勒展开式有 其中， 位于 与 之间，再由条件（3.18）得 即 两边取极限得 如果记 ，则上式还可以记为 * * a b x1 x2 a b 或 x* ?2 x x* 1．计算 在有解区间[a, b]端点处的值， 。 2．计算 在区间中点处的值 。 3．判断若 ，则 即是根，否则检验： 反复执行步骤2、3，便可得到一系列有根区间： 4、当 时 5、则 即为根的近似 ①简单; ② 对f (x) 要求不高(只要连续即可) . ①无法求复根及偶重根 ② 收敛慢 优点 缺点 b x y o a 对二分法，因为数值分析的结果允许含有一定的误差。所以对给定的误差精度 ，只要选择足够大的k，就可使得 这时即可用 作为 的近似值。上式也称为二分法 的误差估计式。 对于事先给定的精度 ，从式(3.4)易求得二分次数 定义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 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 - 例1：求下列方程位于【1，1.5】内的一个根。 例3 利用二分法求方程 ，在 内的根 解：令 取 有 利用上述二分法，可以得到如下的计算结果 （表3.1）。由表3.1可以看出，当 时， ，所以可以将 的中点作为所求根 的近似值。 0.85547 0.85937 0.85156 7 ＋ 0.85156 0.85937 0.84375 6 － 0.85937 0.875 0.84375 5 ＋ 0.84375 0.875 0.8125 4 ＋ 0.8125 0.875 0.75 3 － 0.875 1 0.75 2 ＋ 0.75 1 0.5 1 ＋ 0.5 1 0 0 的符号 取 ，即满足精度要求。 　　二分法是求实根的近似计算中行之有效的最