数值分析3(不动点迭代).pptVIP

*/25 迭代法回顾 二分法: 设f(x) = 0的根为 x*,通过迭代计算产生序列: * 设f(x) = 0的根为 x*,通过迭代计算,产生序列: 迭代的思想 迭代法研究包括方面: 迭代初值 迭代格式 判别收敛及收敛速度 * 《数值分析》3 不动点迭代法 不动点迭代的收敛性 迭代序列的收敛速度 序列收敛加速方法 ? ? ? ? * 引子 选择任意的数字x。Matlab中输入命令x=your number。例如x = 3。 Matlab中输入命令x=sqrt(1+x)。命令计算(1+x)1/2并用最新的结果替代以前的结果。重复上述过程(up-arrow键), 得到如下结果: 3 2 1.7321 1.6529 1.6288 1.6213 1.6191 1.6184 1.6181 1.6181 1.6180 1.6180 1.6180 Matlab中等号是赋值算子。即计算等式右边的值并将值存储到左边的变量。命令计算(1+x)1/2并用最新的结果替代以前的结果。重复上述过程得到如下最后结果为1.6180(黄金比例 )。 而数学中等号的意思有所不同。方程 的根称为函数 f(x)=(1+x)1/2的不动点。函数f(x)有一个不动点(1+(5)1/2)/2。 启示: 两种不同的思路计算方程的根: 1) 寻找根的显式计算公式; 2)通过重复简单的不动点计算并赋值来逐步逼近方程的解。 等价变换 f (x) 的根 的不动点 不动点迭代(Fixed Point Iteration) 对于f(x)=0, 不动点可以构造为: x*称为不 动点 * 不动点迭代(Fixed Point Iteration) 选择适当的初始值x0,按照如下的迭代格式计算: 如果数列{xn}有极限 , 则称迭代是收敛的。 是非线性方程的根和 的不动点。 基本思想是将非线性方程求解归结为一系列显式的函数值计算。 例1 方程 x3 + 4x2 – 10 = 0 在 [1, 2] 上有一个根, 构造求根的不动点迭代格式。 (1) ( n = 0, 1, 2, ······) ( n = 0, 1, 2, ······) (2) * 1.5000 0.8165 2.9969 (-8.65)1/2 0.5487 1.6317 0 1 2 3 n xn | xn – x* | 1.5000 1.3484 1.3674 1.3650 1.3653 1.3652 1.3652 0.0168 0.0022 2.0e-04 1.0e-04 0 0 0 1 2 3 4 5 6 n xn | xn - x* | 什么样的迭代格式收敛？ * 中值定理: 若函数f(x) 满足: (i)在[a,b]连续; (ii)在(a,b)可导; 则在 (a, b) 内至少存在一点 ,使得 * 证 若 或 ,显然 有不动点 设 , 则有 , 记 则有 所以, 存在 x*使得 即 , 故 x* 是 的不动点。 * * 如果 有两个不同的不动点 则有 两式相减得 由拉格朗日中值定理知, 存在 介于 和 之间 ? ? 与L1 条件矛盾 故不动点唯一。 * 证 ? ? * 压缩映像 ( 0L1 ) 所以 故迭代格式收敛 ? ? * * 局部收敛定理 收敛速度 数列的p阶收敛概念 设 , 若存在 a0 , p0 使得 则称数列{xn} p 阶收敛。 特别: (1) 收敛阶p =1时,称为线性收敛; (2) 收敛阶p 1时,称为超线性收敛; (3) 收敛阶p =2 时,称为平方收敛。 序列的收敛阶数越高, 则收敛速度越快。 * 定理2.6 设x*是 的不动点,且 而 则