《第二章非线性方程求根省公共课一等奖全国赛课获奖课件.pptx》由会员分享，可在线阅读，更多相关《第二章非线性方程求根省公共课一等奖全国赛课获奖课件.pptx（53页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第二章第二章 非线性方程求根方法非线性方程求根方法第1页第二章第二章 非线性方程求根方法非线性方程求根方法n n引言引言n n方程求根二分法方程求根二分法n n迭代法及其收敛性迭代法及其收敛性n nNewton迭代法迭代法第2页n n方程是在科学研究中不可缺乏工具；方程是在科学研究中不可缺乏工具；方程是在科学研究中不可缺乏工具；方程是在科学研究中不可缺乏工具；n n方程求解是科学计算中一个主要研究对象；方程求解是科学计算中一个主要研究对象；方程求解是科学计算中一个主要研究对象；方程求解是科学计算中一个主要研究对象；n n几百年前就已经找到了代数方程中二次至五次方程几百年前就已经找到了代数方程中
2、二次至五次方程几百年前就已经找到了代数方程中二次至五次方程几百年前就已经找到了代数方程中二次至五次方程求解公式；求解公式；求解公式；求解公式；n n不过，对于更高次数代数方程当前仍无有效准确解不过，对于更高次数代数方程当前仍无有效准确解不过，对于更高次数代数方程当前仍无有效准确解不过，对于更高次数代数方程当前仍无有效准确解法；法；法；法；n n对于无规律非代数方程求解也无准确解法；对于无规律非代数方程求解也无准确解法；对于无规律非代数方程求解也无准确解法；对于无规律非代数方程求解也无准确解法；n n所以，研究非线性方程数值解法成为必定。所以，研究非线性方程数值解法成为必定。所以，研究非线性方程
3、数值解法成为必定。所以，研究非线性方程数值解法成为必定。2.1 引言引言第3页2.1 引言引言非线性方程普通形式：非线性方程普通形式：非线性方程普通形式：非线性方程普通形式：f f(x x)=0=0代数方程：代数方程：代数方程：代数方程：f f(x x)=a=a0 0+a+a1 1x+ax+an nx xn n （a an n 0 0）超越方程超越方程：f f(x x)中含三角函数、指数函数、或其中含三角函数、指数函数、或其它超越函数。它超越函数。用数值方法求解非线性方程步骤：用数值方法求解非线性方程步骤：（1）找出隔根区间；（只含一个实根区间称隔根区）找出隔根区间；（只含一个实根区间称隔根区
4、间）间）（2）近似根准确化。从隔根区间内一个或多个点出）近似根准确化。从隔根区间内一个或多个点出发，逐次迫近，寻求满足精度根近似值。发，逐次迫近，寻求满足精度根近似值。第4页2.2 方程求根二分法方程求根二分法n n二分法基本思想：二分法基本思想：假定假定f(x)=0在在a,b内有唯一单实根内有唯一单实根x*，考查考查有根区间有根区间a,b，取中点取中点x0=(a+b)/2，若若f(x0)=0，则则x*=x0，不然，不然，uu（1）若）若f(x0)f(a)0，则则x*在在x0右侧，令右侧，令a1=x0，b1=b;uu（2）若）若f(x0)f(a)0，则则x*在在x0左侧，令左侧，令a1=a，b
5、1=x0。n n定理定理1（介值定理）设函数（介值定理）设函数f(x)在区间在区间a,b连续，连续，且且f(a)f(b)0，则方程则方程f(x)=0在区间在区间a,b内最少有内最少有一个根。一个根。第5页以以a1,b1为新隔根区间，且仅为为新隔根区间，且仅为a,b二分之二分之一，对一，对a1,b1重复前过程，得新隔根区间重复前过程，得新隔根区间a2,b2，如此二分下去，得一系列隔根区间：，如此二分下去，得一系列隔根区间：a,b a1,b1 a2,b2 ak,bk 其中每个区间都是前一区间二分之一，故其中每个区间都是前一区间二分之一，故ak,bk 长度：长度：当当k趋于无穷时趋于趋于无穷时趋于0
6、。即若二分过程无限继续。即若二分过程无限继续下去，这些区间最终必收敛于一点下去，这些区间最终必收敛于一点x*，即方即方程根。程根。第6页n n性质：uuf(an)f(bn)0;uubn an=(b a)/2n2.2 方程求根二分法方程求根二分法每次二分后，取有根区间中点每次二分后，取有根区间中点每次二分后，取有根区间中点每次二分后，取有根区间中点x xk k=(a ak k+b bk k)/2)/2作作作作为根近似值，则可得一近似根序列：为根近似值，则可得一近似根序列：为根近似值，则可得一近似根序列：为根近似值，则可得一近似根序列：x x0 0,x x1 1,x x2 2,该序列必以根该序列必
7、以根该序列必以根该序列必以根x x*为极限。为极限。为极限。为极限。实际计算中，若给定充分小正数实际计算中，若给定充分小正数实际计算中，若给定充分小正数实际计算中，若给定充分小正数 0 0和允许误差和允许误差和允许误差和允许误差限限限限 1 1，当，当，当，当|f f(x xn n)|)|0 0或或或或b bn n-a an n 1 1时，均可取时，均可取时，均可取时，均可取x x*x xn n。第7页n n定理定理2 设设x*为方程为方程f(x)=0在在a,b内唯一根，内唯一根，且且f(x)满足满足f(a)f(b)0，则由二分法产生则由二分法产生第第n个个区间区间an,bn 中点中点xn满足
8、不等式满足不等式2.2 方程求根二分法方程求根二分法证实：证实：第8页n n计算过程简单，收敛性可确保；计算过程简单，收敛性可确保；n n对函数性质要求低，只要连续即可。对函数性质要求低，只要连续即可。n n收敛速度慢；收敛速度慢；n n不能求复根和重根；不能求复根和重根；n n调用一次求解一个调用一次求解一个a,b间多个根无法求间多个根无法求得。得。二分法求解非线性方程优缺点：二分法求解非线性方程优缺点：第9页n n不动点迭代法不动点迭代法n n不动点存在性与迭代法收敛性不动点存在性与迭代法收敛性n n迭代收敛加速方法迭代收敛加速方法2.3 迭代法及其收敛性迭代法及其收敛性第10页迭代法基本
9、思想：迭代法基本思想：迭代法是一个逐次迫近方法，用某个固定公式重复校正迭代法是一个逐次迫近方法，用某个固定公式重复校正迭代法是一个逐次迫近方法，用某个固定公式重复校正迭代法是一个逐次迫近方法，用某个固定公式重复校正根近似值，使之逐步准确化，最终得到满足精度要求结根近似值，使之逐步准确化，最终得到满足精度要求结根近似值，使之逐步准确化，最终得到满足精度要求结根近似值，使之逐步准确化，最终得到满足精度要求结果。果。果。果。例：求方程例：求方程例：求方程例：求方程 x x3 3-x x-1=0-1=0 在在在在 x x=1.5=1.5 附近一个根。附近一个根。附近一个根。附近一个根。解：将所给方程改
10、写成解：将所给方程改写成解：将所给方程改写成解：将所给方程改写成假设初值假设初值假设初值假设初值x x0 0=1.5=1.5是其根，代入得是其根，代入得是其根，代入得是其根，代入得x x1 1 x x0 0，再将，再将，再将，再将x x1 1代入得代入得代入得代入得x x2 2 x x1 1，再将，再将，再将，再将x x2 2代入得代入得代入得代入得第11页如此下去，这种逐步校正过程称为迭代过程。这里用如此下去，这种逐步校正过程称为迭代过程。这里用如此下去，这种逐步校正过程称为迭代过程。这里用如此下去，这种逐步校正过程称为迭代过程。这里用公式称为迭代公式，即公式称为迭代公式，即公式称为迭代公式
11、，即公式称为迭代公式，即k k=0,1,2,=0,1,2,迭代结果见下表：迭代结果见下表：迭代结果见下表：迭代结果见下表：k k x xk k k k x xk k0 01 12 23 34 41.51.51.357211.357211.330861.330861.325881.325881.324941.324945 56 67 78 81.324761.324761.324731.324731.324721.324721.324721.32472仅取六位数字，仅取六位数字，仅取六位数字，仅取六位数字，x x7 7与与与与x x8 8相同，即认为相同，即认为相同，即认为相同，即认为x x8
12、8是方程根。是方程根。是方程根。是方程根。x x*x x8 8=1.32472=1.32472第12页2.3.1 不动点迭代法不动点迭代法n n将连续函数方程将连续函数方程将连续函数方程将连续函数方程f f(x x)0 0改写为等价形式：改写为等价形式：改写为等价形式：改写为等价形式：x=x=(x x)其中其中其中其中 (x x)也是连续函数，称为迭代函数。也是连续函数，称为迭代函数。也是连续函数，称为迭代函数。也是连续函数，称为迭代函数。不动点不动点不动点不动点：若：若：若：若x x*满足满足满足满足f f(x*x*)=0 0，则则则则x*=x*=(x*x*)；反之，若反之，若反之，若反之，
13、若x*=x*=(x*x*)，则，则，则，则f f(x*x*)=0 0 ，称，称，称，称x x*为为为为 (x x)一个不动点。一个不动点。一个不动点。一个不动点。不动点迭代不动点迭代不动点迭代不动点迭代：(k k=0,1,)=0,1,)若对任意若对任意若对任意若对任意 x x0 0 a a,b b，由上述迭代得序列由上述迭代得序列由上述迭代得序列由上述迭代得序列 x xk k，有极限，有极限，有极限，有极限则称迭代过程收敛，且则称迭代过程收敛，且则称迭代过程收敛，且则称迭代过程收敛，且x x*=*=(x x*)*)为为为为 (x x)不动点。不动点。不动点。不动点。第13页n n但迭代法并不总
14、令人满意，如将前述方程但迭代法并不总令人满意，如将前述方程但迭代法并不总令人满意，如将前述方程但迭代法并不总令人满意，如将前述方程x x3 3-x x-1=0-1=0改改改改写为另一等价形式：写为另一等价形式：写为另一等价形式：写为另一等价形式：建迭代公式：建迭代公式：建迭代公式：建迭代公式：仍取初值仍取初值仍取初值仍取初值x x0 0=1.5=1.5，则有，则有，则有，则有x x1 1=2.375=2.375，x x2 2=12.39=12.396 6，x x3 3=1904,=1904,结果越来越大。此时称迭代过程发散。结果越来越大。此时称迭代过程发散。结果越来越大。此时称迭代过程发散。结
15、果越来越大。此时称迭代过程发散。x2 x1 x0y=x几何意义：几何意义：第14页收敛收敛 发散发散 第15页n n定理定理3（存在性）（存在性）设设(x)Ca,b且满足以下且满足以下两个条件：两个条件：uu（1）对于任意）对于任意x a,b，有有a(x)b；uu（2 2）若若(x)在在a,b一阶连续，且一阶连续，且存在常存在常数数0 0L1，使得对任意使得对任意x a,b，成立成立|(x)|L则则(x)在在a,b上存在唯一不动点上存在唯一不动点x*。2.3.2 不动点存在性与迭代法收敛性不动点存在性与迭代法收敛性第16页不动点存在性证实：不动点存在性证实：证：证：若若或或显然显然 (x)有不
16、动点；有不动点；不然，设不然，设则有则有（因（因a(x)b）记记则有则有故存在故存在x*使得使得即即x*即为不动点。即为不动点。第17页不动点存在唯一性证实：不动点存在唯一性证实：设有设有 x1*x2*,使得使得则则其中，其中，介于介于 x1*和和 x2*之间。之间。由定理条件由定理条件可得可得矛盾！矛盾！故故 x1*=x2*，不动点唯一存在。，不动点唯一存在。第18页定理定理4（全局收敛性）（全局收敛性）设设(x x)C a a,b b 且满足以下两个条件：且满足以下两个条件：则对任意则对任意x0 a a,b b，由由xn+1=(xn)得到迭代序列得到迭代序列xn 收敛到收敛到(x x)不动
17、点不动点x x*，并有误差预计：并有误差预计：（2）若若(x)在在a,b一阶连续，且存在常数一阶连续，且存在常数0L1，使得对任意使得对任意x a,b，成立，成立|(x)|L（1）对于任意）对于任意x a,b，有有a(x)b；第19页(0L1)故迭代格式收敛。故迭代格式收敛。所以所以证实：证实：第20页重复递推，可得误差预计式重复递推，可得误差预计式第21页n n定义定义 设设(x)有不动点有不动点x*，若存在若存在x*某邻域某邻域R：|x-x*|，对任意对任意x0 R，迭代过程迭代过程xk+1=(xk)产生序列产生序列xk R且收敛到且收敛到x*，则称不动则称不动点迭代法点迭代法xk+1=(
18、xk)局部收敛局部收敛。定理定理4给出收敛性称全局收敛性，实际应给出收敛性称全局收敛性，实际应用时通常只在不动点用时通常只在不动点x*邻近考查其收敛性，称邻近考查其收敛性，称局部收敛性。局部收敛性。第22页n n证实证实：依据连续函数性质，因：依据连续函数性质，因(x)连续，存在连续，存在x*某邻域某邻域R：|x-x*|，对任意对任意x R，|(x)|L1，且且|(x)-x*|=|(x)-(x*)|=|()|x-x*|L|x-x*|x-x*|即对任意即对任意x R，总有总有(x)R。由全局收敛性定义知，迭代过程由全局收敛性定义知，迭代过程 xk+1=(xk)对于对于任意初值任意初值x0 R均收
19、敛。均收敛。n n定理定理5（局部收敛性）（局部收敛性）设设x*为为(x)不动点，不动点，(x)在在x*某邻域连续，且某邻域连续，且|(x*)|11时称时称时称时称超线性收敛超线性收敛超线性收敛超线性收敛。且且且且r r 越大，收敛越快。越大，收敛越快。越大，收敛越快。越大，收敛越快。第26页n n定理定理定理定理：设：设：设：设x x*为为为为x x=(x x )不动点，若不动点，若不动点，若不动点，若 (x x )满足：满足：满足：满足：uu（1 1）(x x )在在在在x x*附近是附近是附近是附近是p p次连续可微次连续可微次连续可微次连续可微(p p1)1)；uu（2 2）则迭代过程
20、则迭代过程则迭代过程则迭代过程x xn n+1+1=(x xn n)在点在点在点在点x x*邻近是邻近是邻近是邻近是p p阶收敛。阶收敛。阶收敛。阶收敛。得得得得所以所以所以所以故故故故迭代过程迭代过程迭代过程迭代过程x xn n+1+1=(x xn n )p p阶收敛。阶收敛。阶收敛。阶收敛。n n证：证：证：证：由由Taylor公式公式第27页uu假定假定(x)改变不大，近似取某个近似值改变不大，近似取某个近似值L，则有则有2.3.3 迭代收敛加速方法迭代收敛加速方法一、一、Aitken加速收敛方法：加速收敛方法：uu由微分中值定理，有由微分中值定理，有同理同理两式相比，得两式相比，得第2
21、8页uu类推可得类推可得uu故故上式即为上式即为Aitken加速收敛方法迭代格式。加速收敛方法迭代格式。第29页将将Aitken加速技巧与不动点结合可得加速技巧与不动点结合可得或将其写为或将其写为二、二、Steffensen迭代法：迭代法：第30页解解：由前知，迭代格式：由前知，迭代格式 xk+1=xk3-1 是发散。现用是发散。现用steffensen迭代法计算。取迭代法计算。取(x)=x3-1，结果以，结果以下：下：k k x xk k y yk k z zk k0 01 12 23 34 45 51.51.51.416291.416291.355651.355651.328951.328
22、951.324801.324801.324721.324722.375002.375001.840921.840921.491401.491401.347101.347101.325181.3251812.396512.39655.238885.238882.317282.317281.444351.444351.327141.32714表明即使不动点迭代法不收敛，用表明即使不动点迭代法不收敛，用steffensen迭迭代法仍可能收敛。代法仍可能收敛。例例：用：用steffensen迭代法求解方程迭代法求解方程 x3-x-1=0。第31页求方程求方程在在3,4中解。中解。例例解解由由取对数取对
23、数结构迭代格式结构迭代格式故故当当x 3,4时，时，(x)3,4，且，且故迭代格式收敛。故迭代格式收敛。第32页取取x0=3.5，经计算可得迭代，经计算可得迭代16次后次后x16=3.73307，有有6位有效数字。位有效数字。若用若用steffensen迭代法加速，结果以下：迭代法加速，结果以下：k k xk k yk k zk k0 03.53.53.604143.604143.662023.662021 13.734443.734443.733813.733813.733473.733472 23.733073.73307说明说明steffensen迭代法收敛速度比不动点迭代迭代法收敛速度
24、比不动点迭代快得多。快得多。第33页n nNewton迭代法及其收敛性迭代法及其收敛性n n简化简化Newton迭代法（平行弦法）迭代法（平行弦法）n n弦截法弦截法n nNewton下山法下山法n n重根情形重根情形2.4 Newton迭代法迭代法第34页uu设方程设方程f(x)=0有近似根有近似根xk（f(xk)0），），将将f(x)在在xk展开：展开：（在在x和和xk之间）之间）2.4.1 Newton迭代法及其收敛性迭代法及其收敛性n n基本思想：基本思想：将非线性方程逐步归结为某种线将非线性方程逐步归结为某种线性方程求解。性方程求解。可设可设第35页记该线性方程根为记该线性方程根为x
25、k+1，则则故故f(x)=0可近似表示为可近似表示为即为即为Newton法迭代格式法迭代格式。（k=0,1,）第36页n nNewton迭代法几何意义迭代法几何意义：（亦称（亦称切线法切线法）xyx*xk+1xkPky=f(x)切线方程切线方程故故第37页n nNewton迭代法收敛性：迭代法收敛性：uu迭代函数：迭代函数：定理定理：假设：假设f(x)在在x*某邻域内含有连续二阶导数，某邻域内含有连续二阶导数，且设且设f(x*)=0，f(x*)0，则对充分靠近则对充分靠近x*初始初始值值x0，Newton迭代法产生序列迭代法产生序列xn最少平方收敛最少平方收敛于于x*。设设f(x*)=0，f(
26、x*)0，则则 (x*)=0，故故Newton迭代法在迭代法在x*附近最少平方收敛。附近最少平方收敛。第38页解解：f(x)=ex+xex，故，故newton迭代公式为迭代公式为k xk 00.51 0.5710220.5671630.56714迭代迭代3次即可得到精度为次即可得到精度为10-5近似解近似解0.56714。若用。若用不动点迭代，到达同一不动点迭代，到达同一精度需精度需17次。次。例例：用：用newton迭代法解方程迭代法解方程 xex-1=0。即即取迭代初值取迭代初值x0=0.5，结果以下：，结果以下：第39页Newton迭代法缺点：迭代法缺点：1.被零除错误被零除错误2.程序
27、死循环程序死循环y=arctan x方程方程:f(x)=x3 3x+2=0在重根在重根x*=1附近，附近，f(x)近似为零。近似为零。对对 f(x)=arctan x存在存在 x0，Newton迭迭代法陷入死循环。代法陷入死循环。第40页uu若若|(x)|=|1-cf(x)|1，即取即取0cf(x)2在在x*附近成立，则收敛。附近成立，则收敛。uu若取若取c=1/f(x0)，则称则称简化简化Newton法法。2.4.2 简化简化Newton法（平行弦法）法（平行弦法）迭代公式：迭代公式：（c 0,k=0,1,）迭代函数：迭代函数：第41页在在Newton迭代格式中，用差商近似导数，迭代格式中，
28、用差商近似导数，2.4.3 弦截法（割线法）弦截法（割线法）称称弦截法弦截法。得得第42页弦截法几何意义：弦截法几何意义：xyx*xk+1xk-1Pk-1y=f(x)xkPk弦线弦线PkPk-1方程：方程：当当y0时，时，第43页n n例例 用简化用简化Newton迭代法和弦截法计算方程迭代法和弦截法计算方程x3-3x+1=0根。根。n n解：解：设设f(x)=x3-3x+1，则，则f(x)=3x2-3由简化由简化Newton法，得法，得由弦截法，得由弦截法，得第44页x x0 0=0.5=0.5x x1 1=0.3333333333=0.3333333333x x2 2=0.34979423
29、87=0.3497942387x x3 3=0.3468683325=0.3468683325x x4 4=0.3473702799=0.3473702799x x5 5=0.3472836048=0.3472836048x x6 6=0.3472985550=0.3472985550 x x7 7=0.3472959759=0.3472959759x x8 8=0.3472964208=0.3472964208x x9 9=0.3472963440=0.3472963440 x x1010=0.3472963572=0.3472963572x x1111=0.3472963553=0.347
30、2963553x x0 0=0.5;=0.5;x x1 1=0.4;=0.4;x x2 2=0.3430962343=0.3430962343x x3 3=0.3473897274=0.3473897274x x4 4=0.3472965093=0.3472965093x x5 5=0.3472963553=0.3472963553x x6 6=0.3472963553=0.3472963553简化简化Newton法法弦截法弦截法要到达精度要到达精度10-8，简化，简化Newton法迭代法迭代11次，弦次，弦截法迭代截法迭代5次，次，Newton迭代法迭代迭代法迭代4次次。第45页不论前面哪种
31、迭代法：不论前面哪种迭代法：（NewtonNewton迭代法、迭代法、迭代法、迭代法、简化简化简化简化NewtonNewton法、法、法、法、弦截法）弦截法）弦截法）弦截法）Newton迭代法迭代法x x0 0=2=2x x1 1=-3.54=-3.54x x2 2=13.95=13.95x x3 3=-279.34=-279.34x x4 4=12=12是否收敛均与初值位置相关。是否收敛均与初值位置相关。如如x x0 0=1=1x x1 1=-0.5708=-0.5708x x2 2=0.1169=0.1169x x3 3=-0.0011=-0.0011x x4 4=7.9631e-010=
32、7.9631e-010 x x5 5=0=0收敛发散第46页n n为预防为预防Newton法发散，可增加一个条件法发散，可增加一个条件：|f(xk+1)|f(xk)|，满足该条件算法称下山法。满足该条件算法称下山法。n n可用下山法确保收敛，可用下山法确保收敛，Newton法加紧速度。法加紧速度。2.4.4 Newton下山法下山法称称Newton下山法下山法。（0 1,下山下山因子）因子）记记即即第47页从从=1开始，逐次减半计算。开始，逐次减半计算。次序，直到使下降条件次序，直到使下降条件|f(xk+1)|f(xk)|成立为成立为止。止。选取：选取：即按即按第48页n n例：求解方程例：求
33、解方程要求到达精度要求到达精度|xn-xn-1|10-5，取，取x0=-0.99。解：先用解：先用Newton迭代法：迭代法：f(x)=x2-1x x2 2=21.69118=21.69118 x x3 3=15.15689 =15.15689 x x4 4=9.70724 =9.70724 x x5 5=6.54091 =6.54091 x x6 6=4.46497 =4.46497 x x7 7=3.13384=3.13384x x8 8=2.32607 =2.32607 x x9 9=1.90230 =1.90230 x x1010=1.75248 =1.75248 x x1111=1.
34、73240 =1.73240 x x1212=1.73205 =1.73205 x x1313=1.73205=1.73205需迭代需迭代需迭代需迭代1313次才次才次才次才到达精度要求到达精度要求到达精度要求到达精度要求第49页n n用用Newton下山法，结果以下：下山法，结果以下：k k=0 =0 x x0 0=-0.99 =-0.99 f f(x x0 0)=0.666567)=0.666567k k=1 =1 x x1 1=32.505829 =32.505829 f f(x x)=11416.4)=11416.4 =0.5 =0.5 x x1 1=15.757915 =15.757
35、915 f f(x x)=1288.5)=1288.5 =0.25 =0.25 x x1 1=7.383958 =7.383958 f f(x x)=126.8)=126.8 =0.125 =0.125 x x1 1=3.196979 =3.196979 f f(x x)=7.69)=7.69 =0.0625 =0.0625 x x1 1=1.103489 =1.103489 f f(x x)=-0.655)=-0.655k k=2 =2 x x2 2=4.115071 =4.115071 f f(x x)=19.1)=19.1 =0.5 =0.5 x x2 2=2.60928 =2.6092
36、8 f f(x x)=3.31)=3.31 =0.25 =0.25 x x2 2=1.85638 =1.85638 f f(x x)=0.27)=0.27k k=3 =3 x x3 3=1.74352 =1.74352 f f(x x)=0.023)=0.023k k=4 =4 x x4 4=1.73216 =1.73216 f f(x x)=0.00024)=0.00024k k=5 =5 x x5 5=1.73205 =1.73205 f f(x x)=0.00000)=0.00000k k=6 =6 x x6 6=1.73205 =1.73205 f f(x x)=0.000000)=0
37、.000000k k 下山因子下山因子下山因子下山因子x xk k f f(x xk k)第50页n n设设f(x)=(x-x*)m g(x),m 2，m为整数，为整数，g(x*)0,则则x*为方程为方程f(x)=0m重根。此时有重根。此时有 f(x*)=f(x*)=f(m-1)(x*)=0,f(m)(x*)02.4.5 重根情形重根情形方法一方法一只要只要f(xk)0，仍可用仍可用Newton法计法计算，此时算，此时为线性收敛。为线性收敛。方法二方法二若取若取则则(x*)=0，用迭代法用迭代法求求m重根，则具重根，则具2阶阶收敛，但要知道收敛，但要知道m。第51页n n方法三方法三还可令还可令则则故故x*是是(x)=0单根，对单根，对(x)用用Newton法，可得法，可得它是二阶收敛。它是二阶收敛。2.4.5 重根情形重根情形第52页n nP43-44习题二：习题二：2，4，6，8，9，10。本章作业第53页