2.6非线性方程近似根.ppt

Altken 迭代法 Altken迭代法 用 差商 近似 微商 设 x* 是方程的根，则由微分中值定理可得 Altken 迭代法 Altken迭代公式 k = 0, 1, 2, ... ... Altken 法同样具有较好的加速效果 * * * * * * * 问题背景和研究目的 次数大于等于5的代数方程没有用系数表达的解析解。 非线性方程求根的一般步骤： 1）根的搜索（根的隔离） 2）根的精确化 2.6 非线性方程近似根 非线性方程没有一般的解析方法。 解方程 是最常见的数学问题之一， 也是数学模型分析求解中不可避免的问题之一。 则可用搜索法求有根区间. x ?1 0 1 2 f(x)的符号 ? ? + + 问题： 求连续的非线性方程 的根的近似值。 根的搜索 根的搜索常用方法： 图解法； 近似方程法； 解析法； 定步长搜索法。 基本思想 二分法 将有根区间进行对分，根据分点处f(x)的符号逐步 将有根区间缩小，使得在足够小的区间内，方程有 且仅有一根。 适用范围 求有根区间内的 单根 或 奇数重实根。 数学原理：连续函数介值定理 设 f(x) 在 [a, b] 上连续，且 f(a) f(b)0，则由介值定理可得，在 (a, b) 内至少存在一点 ? 使得 f(?)=0。 算法 二分法 设方程在区间 [a,b] 内连续，且 f(a)f(b)0，给定精度要求 ? ，若有 |f(x)|? ，则 x 就是f(x) 在区间 (a,b) 内的 近似根。 k ak bk xk f(xk)符号 0 1 2 3 4 5 6 1.0 1.25 1.3125 1.3203 1.5 1.375 1.3438 1.3281 1.25 1.375 1.3125 1.3438 1.3281 1.3203 1.3242 ? + ? + + ? ? 收敛性分析 二分法收敛性 设方程的根为 x* ? (an , bn ) ，又 ，所以 0(n ?) 二分法总是收敛的 二分法的收敛速度较慢 通常用来给出根的一个 较为粗糙的近似。 简单（一步）迭代法 基本思想 从f (x) =0的某个近似值 x0 出发，用1个递推公式，反复校正这个初值，直至满足于先给定的精度要求。 迭代格式的构造 构造 f (x) = 0 的一个等价方程： 从某个初值 x0 出发，构造迭代格式（公式） 得到一个迭代序列 k = 0, 1, 2, ... ... ? (x) 的不动点 f (x) = 0 x = ? (x) 等价变换 f (x) 的零点 ?(x) 称为迭代函数 若 收敛，即 ，假设 ?(x) 连续，则 收敛性分析 迭代法的收敛性 即 注：若迭代序列发散，则迭代法失效！ k xk 0 1 2 3 4 5 6 7 1.5 1.35721 1.33086 1.32588 1.32494 1.32476 1.32473 1.32472 迭代法的收敛性判据 定理2.1：全局收敛 定理2.2：全局发散 定理2.3：局部收敛与发散 定理2.4：收敛速度 定义： 迭代法收敛性判断 定理 2：如果定理 1 的条件成立，则有如下估计 如果存在 x* 的某个邻域 ? =(x*-? , x* +? ), 使得对 ? x0 ? ? 开始的迭代 xk+1 = ?(xk) 都收敛, 则称该迭代法在 x* 附近局部收敛。 定理 1： 设 ?(x) 在某个邻域 ? 内连续，且对 ?x?? 都有 |?’(x)|?L 1, 则迭代局部收敛。 迭代法收敛性判断 定理 3： 设 ，且 对 ? x?[a, b]，有 ?(x)?[a, b]； 对 ? x?[a, b]，有|?’(x)|?L 1； 则 ?x0?[a, b] ，由迭代格式 xk+1 = ?(xk) 得到 的点列都收敛（全局收敛），且 L 越小，迭代收敛越快 k xk 迭代法(1) 迭代法(2) 迭代法(3) 迭代法(4) 0 1 2 3 ? x0 x1 x2 x3 ? 2 3 9 87 ? 2 1.5 2 1.5 ? 2 1.75 1.73475 1.732631 ? 2 1.75 1.732143 1.732051 ? 收敛阶 为了进一步研究收敛速度问题，引入