方程求根q.pptVIP

function y=JDDDL31(x0,epsilon) ek=1;xk=x0; %n=1; for n=1:16 yk=g(xk); xk=yk; ek=abs(yk-xk)/abs(xk); fprintf(\n y%d = %15.14f \n,n,yk) end fprintf(\n ek=|yk-xk|/|xk|= %f \n,ek) y=yk; %%%%%%%%%%%%%%% function Y=g(x) Y=exp(-x); 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）估计收敛速度。 4.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越大，收敛越快，所以说收敛阶是对迭代法收 敛速度的一种度量。 不一定整数 迭代误差 定理4.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} 是线性收敛的. 证 用定理4.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 问：会出现 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