《数值分析》10_YQJ_Ch4非线性方程求根Q.ppt

4.1.5 牛顿法;此式称为简化的Newton公式。 ;? 牛顿公式的另一种形式：线性化;x;例 用牛顿法求方程 x3 -x -1= 0在x0=1.5附近内的根 s.;注：牛顿法收敛性依赖于 x0 的选取。如果 x0 偏离所求的根 x*比较远，则牛顿法可能发散。;定理4.6 (局部收敛性) 设 f ? C 2[a, b]，若 s 为 f (x) 在[a, b]上的根，且 f (s) ? 0，则存在 s 的邻域 B? (s) , 使得任取初值 x0∈B? (s) ，牛顿法产生的序列{ xk } 收敛到 s，且满足;(2)至少平方收敛 由 Taylor 展开;定理4.7 （收敛的充分条件）设 f ?C 2[a, b]，若 (1) f (a) f (b) 0； (2)在整个[a, b]上 f 不变号且 f (x) ? 0； (3)选取 x0 ? [a, b] 使得 f (x0) f (x0) 0；;定理4.7证明 牛顿法的迭代函数为;例3 用牛顿法求方程 x -lnx = 2 在区间(2, ∞)内的根. 解: 在例1中,已经分析此方程在(2, ∞)内仅有一个根 s, 且 s 在 (2,4) 内. 因为, 在[2,4]内取初值 x0=4, 易知 f (x) 在[2,4]内满足定理4.7的条件.下面的牛顿迭代序列收敛于s.;例3 的 Matlab 计算程序与结果;x0 = 3 x1 = 3.147 918 433 002 16 x2 = 3.146 193 440 797 91 x3 = 3.146 193 220 620 59 x4 = 3.146 193 220 620 58 x5 = 3.146 193 220 620 58 x6 = 3.146 193 220 620 58; f (x)的图像; f (x)的图像;4.1.6 求方程 m 重根的牛顿法;?牛顿法的迭代函数 g(x) 为;由此知， f (x)=0 的 m 重根 s 是 x = g(x) 的根；又因 0 g (s) 1， 所以， 依定理4.3知，当 x0 充分靠近 s 时， 牛顿法产生的序列{ xk } 收敛于 s. 但仅有线性速度.;依定理4.3知, 下面变形的牛顿法产生的序列{ xk } 最少以二阶的速度收敛于 s. ;由于 f ( x ) = ( x–s )m g(x), g(s) ? 0, 知;知 s 为 u(x) 的根，且为单根。对 u(x)=0 用牛顿法。 ;例 已知方程 f (x) = 0 有一个三重根 s=0.8, 比较(4.7)和(4.9) 的收敛速度.;子函数 function y=f(x) y=x^4 -1.4*x^3 -0.48*x^2 +1.408*x -0.512; return function y=f1(x) y=4*x*x*x -4.2*x*x -0.96*x +1.408; return function y=f2(x) y=12*x*x -8.4*x -0.96; return;公式(4.7)迭代的结果: Lmnewton(1) x1 = 0.93548387096774 x2 = 0.89135231690659 x3 = 0.86138403395984 x4 = 0.84114516594029 x5 = 0.82753152136394 x6 = 0.81840020161855 x7 = 0.81228741928093 x8 = 0.80820084858679 x9 = 0.80547135880118 x10 = 0.80364941296261 x11 = 0.80243376187056 x12 = 0.80162287288720 x13 = 0.80108207764138 x14 = 0.80072145730991 x15 = 0.80048100366194;使用迭代公式(4.9) 求解, 也取初值 x0=1. ;公式(4.9)迭代的结果： Lmnewton(1) x1 = 0.79999333336173 x2 = 0.79999749269611 x3 = 0.79999749269611 ;4.1.7 割线法 (弦截法);x0;引理4.1 设 f (s)=0, 又 f (x)∈C 2[s- ?, s+ ? ], 且在s 的此邻域内?(x) ? 0, 及 xk-1 , xk?[s- ?, s+ ? ]且 s, xk-1 , xk 互异, 记 ek = s - xk ,