第3章_非线性方程求根解析.ppt

　　下面证明估计式(3.6)，由(3.4)有 于是对任意正整数p有 上述令p→∞, 注意到limxk+p=x* (p→∞)即得(3.6)式. 　　又由于对任意正整数p有 上述令p→∞, 及limxk+p=x* (p→∞)即得(3.7)式. 证毕. 迭代过程是个极限过程. 在用迭代法进行时，必须按精度要求控制迭代次数. 误差估计式(3.6)原则上确定迭代次数，但它由于含有信息L而不便于实际应用. 而误差估计式(3.7)是实用的，只要相邻两次计算结果的偏差足够小即可保证近似值xk具有足够精度. * 注1：(3.6)式可用来估计迭代次数，但结果偏保守，次数偏大，一般用得不多。称为事前误差估计。 注2：(3.7)式可用在程序中置退出迭代的条件， 即当|xk-xk-1|?时，认为|x*-xk-1|? ，编程计算时常用.称为事后误差估计. 注3：若 ，则 收敛于x速度慢；若 ，则收 敛于x速度快。 注4：若 ，则 不收敛于[a,b]内的不动点x。 　　例如，在前面例3中采用的三种迭代公式，在隔根区间(1, 1.2)内，有 故前两个迭代公式收敛，第三个迭代公式不收敛. 3.4.3 局部收敛性与收敛阶 上面给出了迭代序列{xk}在区间[a, b]上的收敛性,通常称为全局收敛性. 有时不易检验定理的条件，实际应用时通常只在不动点x*的邻近考察其收敛性，即局部收敛性. 定义3-3 设?(x)有不动点x*，如果存在x*的某个邻域R: |x-x*|≤δ，对任意x0∈R，迭代公式(3.2)产生的序列{xk}∈R，且收敛到x*，则称迭代法(3.2)局部收敛. 定理3-4 设x*为?(x)的不动点，　　 在x*的某个邻域连续，且 ，则迭代法(3.2)局部收敛. 证明 由连续函数的性质，存在不动点x*的某个邻域R: |x-x*|≤δ，使对于任意x∈R成立 此外，对于任意x∈R，总有?(x)∈R，这时因为 于是依据定理3.2可以断定迭代过程xk+1=?(xk)对于任意初值x0∈R均收敛. 证毕. 例4 只用四则运算不用开放求方程x2-3=0的根 . 解 这里f(x)= x2-3，可以改写为各种不同的等价形式x=?(x)，其不动点为 ，由此构造不同的迭代法. 取x0=2, 对上式4种迭代法, 计算三步所得结果入下表. 2 1.75 1.732143 1.732051 ┆ 2 1.75 1.73475 1.732361 ┆ 2 1.5 2 1.5 ┆ 2 3 9 87 ┆ x0 x1 x2 x3 ┆ 0 1 2 3 ┆ 迭代法(4) 迭代法(3) 迭代法(2) 迭代法(1) xk k 注意 ，从计算结果看到迭代法(1)及(2)均不收敛，且它们均不满足定理3中的局部收敛条件，迭代法(3)和(4)均满足局部收敛条件，且迭代法(4)比(3)收敛快，因在迭代法(4)中??(x*)=0. 为了衡量迭代法收敛速度的快慢先讲3.5节. 定义3-4 设迭代过程xk+1=?(xk)收敛于方程x=?(x)的根x*，如果迭代误差ek=xk-x*当k→∞时成立下列渐近关系式 则称该迭代法是p阶收敛的. 特别地，p=1时称线性收敛，p1时称超线性收敛，p=2时称平方收敛. 3.5 迭代法的收敛阶 定理3-6 对于迭代过程xk+1=?(xk)，如果?(p)(x)在所求根x*的邻近连续，并且 则该迭代过程在x*的邻近是p阶收敛的. 证明 由于??(x*)=0，根据定理3立即可以断定迭代过程xk+1=?(xk)具有局部收敛性. 再将?(xk)在根x*处做泰勒展开, 利用条件(3.8), 则有 注意到?(xk)=xk+1，?(x*)= x*，由上式得 因此对迭代误差，令k→∞时有 这表明迭代过程xk+1=?(xk)确实为p阶收敛. 证毕. 上述定理告诉我们，迭代过程的收敛速度依赖于迭代函数?(x)的选取. 如果x∈[a, b]但??(x)≠0时，则该迭代过程只可能是线性收敛. 的三阶方法. 假设 x0 充分靠近 x*, 求 证明 首先由泰勒展式可得 例子 证明迭代公式 xk+1=xk(xk2+3a)/(3xk2+a)是求 而1/4a≠0，故此迭代