数值分析对分法和一般迭代法.ppt

举例 用一般迭代法求x3-x-1=0的正实根x* 容易得到：g′(x)在包含x*的某邻域U(x*) 内 连续，且|g′(x*)|1 例题 用一般迭代法求方程x-lnx＝2在区间（2，?）内的根，要求|xk-xk-1|/|xk|=10-8 解：令f(x)=x-lnx-2 f(2)0,f(4)0,故方程在（2,4）内至少有一个根 将方程化为等价方程：x＝2＋lnx 因此，? x0?（2，?），xk+1＝2＋lnxk产生的序列? xk ?收敛于X* 取初值x0＝3.0,计算结果如下： 7 3.146143611 8 3.146177452 9 3.146188209 10 3.146191628 11 3.146192714 12 3.146193060 13 3.146193169 14 3.146193204 k xi 0 3.000000000 1 3.098612289 2 3.130954362 3 3.141337866 4 3.144648781 5 3.145702209 6 3.146037143 另一种迭代格式： 0 3.000000000 1 3.147918433 2 3.146193441 3 3.146193221 程序演示 由此可见，对同一个非线性方程的迭代格式，在收敛的情形下，有的收敛快，有的收敛慢。 定义1. :设序列{xk}收敛于x*，若存在p≥1和正数c， 使得成立 则称{xk}为 p 阶收敛的 特别, p = 1，要求c1, 称线性收敛; 1p2,称超线性收敛 p=2,称平方收敛。 迭代法的收敛阶(收敛速度) 定理4.2.3 设x*为g的不动点，p≥2为正整数，g在x*的某邻域Ｕ(x*)内p阶连续可微，且 g′(x*)=g″(x*)=…=g (p-1)(x*)=0，而g(p)(x*)≠0， 则存在?0，当x0∈ ［x* - ? ，x*+ ?］ (x0≠x*)时，由迭代法(3)产生的序列｛xk｝以p阶收敛速度收敛于x*. Proof: (1)由g′(x*)=0?必存在?0，当x0∈ ［x* - ? ，x*+ ?］ ? U(x)时，由迭代格式(3)产生的序列{xk}收敛于x*,并有xk∈ ［x* - ? ，x*+ ?］ (2)由泰勒公式有xk+1=g(xk)=g(x* ）＋g′(x*)(xk- x*)+…+g (p-1) (x*)(xk-x*) p-1/(p-1)! + g (p)(x*+ ?(xk-x*))(xk-x*) p /p! ，0?1. 利用g在x*的各阶导数条件及g(x*)=x*,上式可改写成 (11) (3)由于g在x*处p阶连续可微且g(p)(x*)≠0，知必存在x*的某邻域Ｕ(x*)，当x∈U(x*)时，有g (p) (x)≠0. 由于x*+ ?(xk-x*) ∈ ［x* - ? ，x*+ ?］ ?U(x*)，故 g (p)(x*+ ?(xk-x*)) ≠0，k=0，1,2,…. 可见，当初值x0≠x*时，由(11)式可推出xk≠x* 于是由(11)式有 上式令k→∞取极限. 即{xk}有p阶收敛速度. * 方程是在科学研究中不可缺少的工具.方程求解是科学计算中一个重要的研究对象. 几百年前就已经找到了代数方程中二次至五次方程的求解公式.但是,对于更高次数的代数方程目前仍无有效的精确解法.对于无规律的非代数方程的求解也无精确解法.因此,研究非线性方程的数值解法成为必然. 第四章 非线性方程的数值解法  非线性方程求根 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 打 印 结