数值分析3.1.二分法、迭代法及收敛性.pptVIP

显然f(x)∈C[a, b],且满足f(a)=?(a)-a0, f(b)=?(b)-b0, 由连续函数性质可知存在 x*∈(a, b) 使 f(x*)=0，即x*=?(x*)，x*即为?(x)的不动点. 再证不动点的唯一性. 设x1*, x2*∈[a, b]都是?(x)的不动点，则由(2.4)得 引出矛盾，故?(x)的不动点只能是唯一的.　证毕. 定理2 设?(x)∈C[a, b]满足定理1中的两个条件，则对任意x0∈[a, b]，由(2.2)得到的迭代序列{xk}收敛到不动点x*，并有误差估计式 证明 设x*∈[a, b]是?(x)在[a, b]上的唯一不动点,由条件1o，可知{xk}∈[a, b]，再由(2.4)得 因0L1，故当k→∞时序列{xk}收敛到x*. 　　下面证明估计式(2.5)，由(2.4)有 于是对任意正整数 p 有 上述令p→∞, 注意到limxk+p=x* (p→∞)即得(2.5)式. 　　又由于对任意正整数 p 有 上述令p→∞, 及limxk+p=x* (p→∞)即得(2.6)式. 证毕. 注：误差估计式(2.5)原则上确定迭代次数，但它由于含有信息 L 而不便于实际应用. 而误差估计式(2.6)是实用的，只要相邻两次计算结果的偏差足够小即可保证近似值 xk 具有足够精度. 注： 对定理1和定理2中的条件2o可以改为导数，即在使用时如果?(x)∈C[a, b]且对任意x∈[a, b]有 则由微分中值定理可知对任意x, y∈[a, b]有 故定理中的条件2o是成立的.　 　　例如，在前面例2中采用的三种迭代公式，在有根区间(1, 1.2)内，有 故前两个迭代公式收敛，第三个迭代公式不收敛. 3.2.3 局部收敛性与收敛阶 上面给出了迭代序列{xk}在区间[a, b]上的收敛性,通常称为全局收敛性. 有时不易检验定理的条件，实际应用时通常只在不动点 x* 的邻近考察其收敛性，即局部收敛性. 定义1 设?(x)有不动点x*，如果存在 x* 的某个邻域U: |x-x*|≤δ，对任意x0∈U，迭代公式(2.2)产生的序列{xk}∈U，且收敛到 x*，则称迭代法(2.2)局部收敛. 定理3 设x*为?(x)的不动点，　　 在x*的某个邻域连续，且 ，则迭代法(2.2)局部收敛. 证明 由连续函数的性质，存在不动点 x* 的某个邻域 U: |x-x*|≤δ，使对于任意 x∈U 成立 此外，对于任意 x∈U，总有?(x)∈U，这是因为 于是依据定理2可以断定迭代过程 xk+1=?(xk) 对于任意初值 x0∈U 均收敛. 证毕. 例3 用不同迭代法求方程x2-3=0的根 . 解 这里 f(x)= x2-3，可以改写为各种不同的等价形式 x=?(x)，其不动点为 ，由此构造不同的迭代法. 取x0=2, 对上式4种迭代法, 计算三步所得结果入下表. 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.732361 ┆ 2 1.75 1.732143 1.732051 ┆ 定义2 设迭代过程 xk+1=?(xk) 收敛于方程 x=?(x)的根 x*，如果迭代误差 ek=xk-x*当 k→∞时成立下列渐近关系式 则称该迭代法是 p阶收敛的. 特别地，p=1时称线性收敛，p1时称超线性收敛，p=2 时称平方收敛. 定理4 对于迭代过程 xk+1=?(xk)，如果?(p)(x)在所求根 x* 的邻近连续（ p 1），并且 则该迭代过程在 x* 的邻近是 p 阶收敛的. 证明 由于??(x*)=0，根据定理3立即可以断定迭代过程xk+1=?(xk)具有局部收敛性. 再将?(xk)在根x*处做泰勒展开, 利用条件(2.8), 则有 注意到?(xk)=xk+1，?(x*)= x*，由上式得 因此对迭代误差，令k→∞时有 这表明迭代过程xk+1=?(xk)确实为 p 阶收敛. 证毕. 注：上述定理告诉我们，迭代过程的收敛速度依赖于迭代函数?(x)的选取. 如果x∈[a, b]但 ??(x)≠0 时，则该迭代过程只可能是线性收敛. 的三阶方法. 假设 x0 充分靠近 x*, 求 提示： 设?(x) =x(x2+3a)/(3x2+a)，