第二章—非线性方程求根.ppt

§6 解非线性方程组的迭代法 * 2021/9/12 在迭代点x(k),将各函数作一阶Taylor展开 另展开式等于零,得到下一个迭代点x(k+1) * 2021/9/12 k x (k) 0 1 2 3 (1.5, 1.0)T (1.5, 0.75)T (1.488095, 0.755952)T (1.488034, 0.755983)T * 2021/9/12 §5.6 迭代法的收敛阶和Aitken加速方法 一. 迭代过程的收敛速度 二. 常见迭代过程的收敛阶 三. 埃特金(Aitken)加速法 四. 斯蒂芬森(Steffensen)法 * 2021/9/12 §6 迭代法的收敛阶和Aitken加速方法 一. 迭代过程的收敛速度 一种迭代法要具有实用价值,不但要肯定它是收敛的,还要求它收敛的比较快。所谓迭代过程的收敛速度，是指在接近收敛时迭代误差的下降速度，具体地说，如果迭代误差 当 时成立 则称迭代过程是 p 阶收敛的。 * 2021/9/12 特别地, p = 1 (且 0 C 1 )时称线性收敛, p 1 时称超线性收敛， 其中 p = 2 时称二次收敛 （或称为平方收敛）。 p 越大，{ xk } 收敛于 x* 的速度就越快。 p 值的大小是衡量一个迭代过程优劣的标志之一。 * 2021/9/12 二. 常见迭代过程的收敛阶 一般迭代法:p＝1,线性收敛 牛顿法:p＝2,二阶收敛（当x*为二重根时，线性收敛） 弦割法：p=1.618，超线性收敛 对收敛较慢数列{ xk }，一个补救的方法是采用加速公式。 * 2021/9/12 三. 埃特金(Aitken)加速法 * 2021/9/12 §5.6 迭代法的收敛阶和Aitken加速方法 x y y = x y = g(x) x* x0 P(x0, x1) x1 x2 P(x1, x2) * 2021/9/12 例 用埃特金方法求方程 在初值 附近的一个根, 精度要求 , 取迭代格式 解 埃特金方法迭代格式为 只迭代二次就得到满足精度要求的解。 §2.2 二分法 0 x y X* x0 a b y=f(x) a1 b1 0 x y X* x0 a b y=f(x) a1 b1 0 x y X* x0 a b y=f(x) a1 b1 0 x y X* x0 a b y=f(x) a1 b1 §5.2 二分法 * 2021/9/12 §5.2 二分法 * 2021/9/12 §5.2 二分法 * 2021/9/12 二分法优缺点 优点:计算简单,方法可靠,只要求f (x)连续，在两个点上异号。 缺点:不能求偶数重根, 也不能求复根, 收敛速度不算太快(与以1/2为比值的等比级数相同)。 因此，一般在求方程近似根时，不单独使用，常用来为其它方法提供好的初值。 对于方程求根，最常用方法是迭代法。 * 2021/9/12 f (x) = 0 x = g (x) 等价变换 f (x) 的根 g (x) 的不动点 §2.3 迭代法 迭代法 /* Fixed-Point Iteration */ * 2021/9/12 一. 迭代格式的构造 * 2021/9/12 §5.3 迭代法 二. 迭代过程的几何表示 O x* x2 x1 x0 x y y=x y=g(x) P0 P1 P2 P* Q1 Q2 * 2021/9/12 例1:用迭代法求方程f(x)=x2-2x-3=0的根（x1=3,x2= -1） 解:(1)方程改写成 x=(2x+3)1/2 建立迭代公式 xk+1=(2xk+3)1/2 (k=0,1,2…), 取x0=4，x1=3.316，x2= 3.104，x3=3.034， x4= 3.011，x5= 3.004 当k?∞， xk ?3，收敛; (2)方程改写成 x=1/2*(x2-3) 建立迭代公式 xk= 1/2 *(xk 2-3) (k=0,1,2…)， 取x0=4， x1=6