数值分析—第2章非线性方程的数值解法.pptVIP

非线性方程（组）的数值解 单变量非线性方程求解：[x,fval]=fzero(f,x0,tol) 该函数采用迭代法计算函数f(x)的一个零点，迭代初值为x0，当两次迭代结果小于tol时停止迭代过程。tol的缺省值是eps。 注意: 在调用函数fzero 之前，要使用m文件建立自己要计算的函数f(x)，只有定义了函数f(x)的m文件后，才能在fzero函数的参数中使用自定义函数名。 例： 求f(x)=x-1/x+5 在x0=-5作为迭代初值时的零点。 先编制一个函数文件fz.m： function f=fz(x) f=x-1/x+5; 然后，在MATLAB命令窗口，输入命令： fzero(fz,-5) 2.3 加速收敛迭代法 2.3.1 Aitken （艾特肯）加速迭代法 设 为 的第k次迭代值， 迭代一次： 再迭代： 前述： 同理： 两者消去 ： 则序列{xk}满足 它表明序列{?xk}的收敛速度比{xk}的收敛速度快. 定理2-5 设序列{xk}线性收敛于x*，对一切k≥0 ， ek = x* -xk ≠0, 且 2.3.2 Steffensen(斯蒂芬森)迭代法 艾特肯方法不管原序列{xk}是怎样产生的，对{xk}进行加速计算，得到序列{?xk}. 如果把艾特肯加速技巧与不定点迭代结合，则可得到如下的迭代法： 实际上上述迭代是将不定点迭代法计算两步合并成一步得到的，可将它写成另一种不动点迭代 其中 ? Aitken 加速： x y y = x y = g(x) x* x0 P(x0, x1) x1 x2 P(x1, x2) 一般地，有： 比 收敛得略快。 ? Steffensen 加速： 定理2-6 （1）若x*为? 的不动点， ?? 在x*为领域连续且满足?? ( x*)≠1，则x*是Ψ(x)的不动点；反之，若x*为Ψ(x)的不动点，则x* 也是?(x) 的不动点。 （2）若x*是? 的不动点， ???? (x)在x*领域连续，且??( x*)≠1，则 Steffensen迭代法至少2阶收敛. 例2-6 求解方程 在x0 =1.5附近的根。 解 用迭代公式 求解方程是发散的。 采用Steffebsen迭代法： 12.3965 5.23888 2.31728 1.44435 1.32714 zk 2.37500 1.84092 1.49140 1.34710 1.32518 yk 1.5 1.41629 1.35565 1.32895 1.32480 1.32472 xk 0 1 2 3 4 5 k 2.4 Newton迭代法 2.4.1 牛顿迭代法 原理：将非线性方程线性化 —— Taylor 展开 /* Taylor’s expansion */ 取 x0 作为初始近似值，将 f (x)在 x0 做一阶Taylor展开: ，? 在 x0 和 x 之间。 线性 作为第一次近似值 重复上述过程 ? Newton 迭代公式 只要 f ?C1，每一步迭代都有f ’( xk ) ? 0， 而且 ，则 x*就是 f 的根。 高阶小量 　　设xk是根x*的某个近似值，过曲线y=f(x)上横坐标为xk的点Pk引切线，并将该切线与x轴交点的横坐标xk+1作为x*的新的近似值. 牛顿法的几何意义 x y x* x0 x 1 x 2 牛顿法也称为切线法 以直代曲 解 令f(x)=xex-1，则f?(x)=ex(1+x) ， 牛顿迭代公式为 取 x0=0.5，迭代得 x1=0.57102, x2=0.56716, x3=0.56714, … 例2-7 用牛顿迭代法求方程xex-1=0的近似根. 故方程的近似根为x*=0.5671 普通迭代法18次才能得到的计算结果。 定理2-7 设f?C2[a, b], 若x*为f(x)=0在[a, b]上的根，且f?(x*)?0，则牛顿迭代法是二阶收敛的，且 设x*是f(x)的一个单根，即f(x*)=0，f?(x*)≠0, 有 证明 牛顿迭代法实际上是一种特殊的不动点迭代，迭代函数为 平方收敛 由 Taylor 展开： 只要 f ’(x*) ? 0，则令 可得结论。 在单根 附近收敛快 由定理2-3 ? 问题 Q1: 若　　　 ，Newton’s Method 是否仍收敛？ 设 x* 是 f 的 m 重根，则： 且