1方程求根分解.pptVIP

x0=0.5 x0 = 0.5000 epsilon=0epsilon = 1.0000e-008 y=JDDDL31(x0,epsilon) y1 = 0.60653065971263 y2 = 0.54523921189261 y3 = 0.57970309487807 y4 = 0.56006462793890 y5 = 0.57117214897722 y6 = 0.56486294698032 y7 = 0.56843804757007 y8 = 0.56640945274692 y9 = 0.56755963426224 y10 = 0.56690721293547 y11 = 0.56727719597078 y12 = 0.56706735185373 y13 = 0.56718636008764 y14 = 0.56711886425699 y15 = 0.56715714370764 y16 = 0.56713543365927 ek=|yk-xk|/|xk|= 0.000000 y = 0.5671 练习 为求方程 在区间[1.3，1.6]内的一个根，将方程改写成下列等价形式，并建立相应的迭代公式： （1） 迭代公式 （2） 迭代公式 求：（i）分析每种迭代公式的收敛性； （ii）估计收敛速度。(接下来讲到) ? 不能用判别法， 由计算知道是发散 2.1.3 简单迭代法的收敛速度 ?对于一种迭代过程，为了保证它是有效的，需要肯定它的收敛性，同时考察它的收敛速度。 ?收敛阶：在接近收敛的过程中，迭代误差的下降速度。 收敛阶是对迭代法收敛速度的一种描述。 定义 设迭代 xk+1 = g(xk) 收敛到g(x) 的不动点 s. 又 ek = xk ? s，若有常数 C 和 p≥1 , 使得 ?若 p = 1 , 称{xk}为线性收敛, 这时 0 C ≤1。 ? p1, 称{xk}为超线性收敛; p=2, 称其为平方收敛. 则称序列{xk}为 p 阶收敛。常数 C 叫收敛因子. 注：数p的大小反映了迭代法的收敛速度的快慢， P越大，收敛越快，所以说收敛阶是对迭代法收 敛速度的一种度量。 不一定整数 迭代误差 定理2.3 设函数 g(x)?C[a, b], g(x)?C(a, b), 且 (I) 当 x?[a, b]时, g(x)?[a, b]; 即, (II) 当 x?(a, b) 时, g(x) ≠0, | g(x) | ? L 1, 其中L 是一个常数. 则对任取的 x0?[a, b], 由 xk+1 = g(xk) 得到的序列 { xk} 收敛于 g(x) 在[a, b]上的唯一不动点 s .并且, 当 x0 ≠ s 时, { xk} 是线性收敛的. 证 用定理2.1, 知{ xk} 收敛于 g(x) 在[a, b]上的唯一不动点 s. 再证, 当 x0 ≠ s 时, 序列{ xk} 是线性收敛的. 收敛阶的判别定理1 由Taylor 公式, 有 xk+1=g(xk)= g(s) + g[s+? (xk - s )](xk-s) , 记 ek = xk ? s, 由上式可得, ek+1 = g[s+? ek]ek , 故有 因为| g(s) | 0, 依定义知, { xk} 是线性收敛的. ek≠ 0 ? s 1.问：会出现 ek = xk ? s= 0 吗? 由于 x?(a, b)， g (x) ≠0， 故当 x0 ≠ s 时， 必有 xk ≠ s ， k = 0,1, 2, ...， 这是因为， x1-s = g(x0) -g(s) = g(?1)(x0-s) ≠0,?1?(a, b), x1≠ s . x2 - s = g(x1) -g(s) = g(?2)(x1-s) ≠0, x2≠ s . ...... 2.从定理可知，迭代函数在s处的一阶导不等0，相应的简单迭代法只是线性收敛 定理2.4 设 s为g(x)的不动点， g(p)(x) 在含 s的某个开区间内连续 ( p≥2 )。 如果 g(k)(x) =0，( k =1, 2, ..., p-1)， g(p)(x) ≠ 0， 则存在 ? 0，对于 x0?B?=[s- ?, s+ ?]但 x0 ≠ s ，使得序列 xk+1= g(xk) 在 B? 内 p 阶收敛