养老保险问题建模分析ppt课件

1、科学计算与数学建模1；.第四章第四章 养老保险问题养老保险问题 非线性方程求根的数值解法非线性方程求根的数值解法养老保险问题养老保险问题4.1非线性方程求根的数值方法非线性方程求根的数值方法4.2养老保险模型的求解养老保险模型的求解4.32；.4.1.1 4.1.1 问题的引入问题的引入 养老保险是保险中的一种重要险种，保险公司将提供不同的保险方案供以选择，养老保险是保险中的一种重要险种，保险公司将提供不同的保险方案供以选择，分析保险品种的实际投资价值。也就是说，如果已知所交保费和保险收入，则按年分析保险品种的实际投资价值。也就是说，如果已知所交保费和保险收入，则按年或按月计算实际的利率是多少

2、？或者说，保险公司需要用你的保费实际至少获得多或按月计算实际的利率是多少？或者说，保险公司需要用你的保费实际至少获得多少利润才能保证兑现你的保险收益？少利润才能保证兑现你的保险收益？4.1 4.1 养老保险问题养老保险问题3；.4.1.2 4.1.2 模型分析模型分析 假设每月交费假设每月交费200200元至元至6060岁开始领取养老金，男子若岁开始领取养老金，男子若2525岁起投保，届时养老金每岁起投保，届时养老金每月月22822282元；如元；如3535岁起保，届时月养老金岁起保，届时月养老金10561056元；试求出保险公司为了兑现保险责任，元；试求出保险公司为了兑现保险责任，每月至少应

3、有多少投资收益率？这也就是投保人的实际收益率。每月至少应有多少投资收益率？这也就是投保人的实际收益率。4；.4.1.3 4.1.3 模型假设模型假设 这应当是一个过程分析模型问题。过程的结果在条件一定时是确定的。整个过程这应当是一个过程分析模型问题。过程的结果在条件一定时是确定的。整个过程可以按月进行划分，因为交费是按月进行的。假设投保人到第月止所交保费及收益可以按月进行划分，因为交费是按月进行的。假设投保人到第月止所交保费及收益的累计总额为，每月收益率为，用分别表示的累计总额为，每月收益率为，用分别表示6060岁之前和之后每月交费数和领取数，岁之前和之后每月交费数和领取数，N N表示停交保险

4、费的月份，表示停交保险费的月份，M M表示停领养老金的月份。表示停领养老金的月份。5；.4.1.4 4.1.4 模型建立模型建立v在整个过程中，离散变量的变化规律满足：在整个过程中，离散变量的变化规律满足：v在这里实际上表示从保险人开始交纳保险费以后，保险人账户上的资金数值。在这里实际上表示从保险人开始交纳保险费以后，保险人账户上的资金数值。 MNkqrFFNkprFFkkkk,.,)1 (1,.,1 ,0,)1 (116；.4.1.4 4.1.4 模型建立模型建立v 我们关心的是，在第我们关心的是，在第M M月时，月时， FK FK 能否为非负能否为非负 数？如果为正，则表明保险公司数？如果

5、为正，则表明保险公司获得收益；如为负，则表明保险公司出现亏损。当为零时，表明保险公司最后一无获得收益；如为负，则表明保险公司出现亏损。当为零时，表明保险公司最后一无所有，所有的收益全归保险人，把它作为保险人的实际收益。所有，所有的收益全归保险人，把它作为保险人的实际收益。 v 从这个分析来看，引入变量从这个分析来看，引入变量FKFK，很好地刻画了整个过程中资金的变化关系；特，很好地刻画了整个过程中资金的变化关系；特别是引入收益率别是引入收益率 r r，虽然它不是我们所求的保险人的收益率，但从问题系统环境中，虽然它不是我们所求的保险人的收益率，但从问题系统环境中来看，必然要考虑引入另一对象来看，

6、必然要考虑引入另一对象保险公司的经营效益，以此作为整个过程中各保险公司的经营效益，以此作为整个过程中各量变化的表现基础。量变化的表现基础。7；.4.1.54.1.5 模型求解模型求解在在(4.1.4)(4.1.4)中两式，取初始值，我们可以得到：中两式，取初始值，我们可以得到：MNkrrqrFFNkrrprFFNkNkNkkkk,.,1,1)1()1(,.,2 , 1 , 0,1)1()1 (0再分别取，再分别取，k=Nk=N和和k=Mk=M，并利用，并利用FM=0FM=0可以求出：可以求出：它是一个非线性方程。它是一个非线性方程。0)1)(1 ()1 (pqrpqrNMM8；. 代数方程求根

7、问题是一个古老的数学问题。早在代数方程求根问题是一个古老的数学问题。早在1616世纪就找到了三次、四次方程世纪就找到了三次、四次方程的求根公式。但直到的求根公式。但直到1919世纪才证明了世纪才证明了 次的一般代数方程式是不能用代数公式求次的一般代数方程式是不能用代数公式求解的，因此需要研究用数值方法求得满足一定精度的代数方程式的近似解。解的，因此需要研究用数值方法求得满足一定精度的代数方程式的近似解。 在工程和科学技术中许多问题常归结为求解非线性方程式问题。正因为非线性方程在工程和科学技术中许多问题常归结为求解非线性方程式问题。正因为非线性方程求根问题是如此重要和基础，因此它的求根问题很早就

8、引起了人们的兴趣，并得到求根问题是如此重要和基础，因此它的求根问题很早就引起了人们的兴趣，并得到了许多成熟的求解方法。下面就让我们首先了解一下非线性方程的基本概念。了许多成熟的求解方法。下面就让我们首先了解一下非线性方程的基本概念。5n 9；.4.2.1 根的搜索相关定义根的搜索相关定义定义定义4.2.1 4.2.1 设有一个非线性方程设有一个非线性方程 ，其中，其中 为实变量为实变量 的非线性函数。的非线性函数。（1 1）如果）如果 有使有使 ，则称，则称 为方程的根，或为为方程的根，或为 的零点。的零点。（2 2）当）当 为多项式，即为多项式，即 则称则称 为为 次代数方程，次代数方程，

9、包含指数函数或者三角函数等特殊函数时，则称包含指数函数或者三角函数等特殊函数时，则称 为特殊为特殊方程。方程。（3 3）如果）如果 ，其中，其中 。 为正整数，则称为正整数，则称 为为 的的 重重根。当根。当 时称时称 为为 的单根。的单根。 0f x ( )f xx f xx( ) 0f xx f x110( )00nnnnnf x ax a xax aa ( ) 0f x n( )f x( ) 0f x ( ) ()( )mf xx xg x( )0g xmx( ) 0f x m1mx( ) 0f x 4.2 4.2 非线性方程求根的数值方法非线性方程求根的数值方法10；.定理定理4.2.

10、14.2.1 设设 为具有复系数的为具有复系数的 次代数方程，则次代数方程，则 在复数域上恰有在复数域上恰有 个根个根（ 重根计算重根计算 个）。如果个）。如果 为实系数方程，则复数根成对出现，即当为实系数方程，则复数根成对出现，即当: : 为为 的复根，则的复根，则 亦是亦是 的根。的根。定理定理4.2.24.2.2 设设 在在 连续连续, ,且且 ，则存在，则存在 ，使得，使得 ，即即 在在 内存在实零点。内存在实零点。( )f x,a b( )( ) 0f a f b,xa b( )0f x f x, )a b（( ) 0f x n( ) 0f x nrr( )0f x 0i ( ) 0

11、f x ( ) 0f x i11；.4.2.2 4.2.2 逐步搜索法逐步搜索法v 对于方程对于方程 ， ，为明确起见，设，为明确起见，设 , , ，从区间左端，从区间左端点点 ，出发按某个预定步长，出发按某个预定步长 （如取（如取 ， 为正整数），一步一步地向为正整数），一步一步地向右跨，每跨一步进行一次根的收索。即检查节点右跨，每跨一步进行一次根的收索。即检查节点 上的函数值上的函数值 的符号，若的符号，若 ，则，则 即为方程解。若即为方程解。若 ，则方，则方程根在区间程根在区间 中，其宽度为中，其宽度为 。 bahN 0f x ,xa b 0f a ( )0f b 0 xahNkxa k

12、h kf x 0kf x kx 0kf x 1,kkxxh12；.4.2.2 4.2.2 逐步搜索法逐步搜索法例例4.2.14.2.1 考察方程考察方程 由于由于 则则 在在 内至少有一个根，设从内至少有一个根，设从 出发，以出发，以 为步长向右进行根的搜索。列表记录各节点函数值的符号。可见在为步长向右进行根的搜索。列表记录各节点函数值的符号。可见在 内必有一内必有一根。根。 表表4.2.14.2.1 的符号的符号x00.51.01.5 的符号-+ 31 0f xx x 01 0,25 0ff f x0,20 x 0.5h1.0,1.5 f x( )f x13；.4.2.2 4.2.2 逐步搜

13、索法逐步搜索法 易见此方法应用关键在步长易见此方法应用关键在步长 的选择上。很明显，只要步长的选择上。很明显，只要步长 取得足够小，利取得足够小，利用此法就可以得到任意精度的根，但用此法就可以得到任意精度的根，但 缩小，搜索步数增多，从而使计算量增大，缩小，搜索步数增多，从而使计算量增大，用此方法对高精度要求不合适。用此方法对高精度要求不合适。hhh14；.4.2.3 4.2.3 二分法二分法对非线性方程：对非线性方程： 其中其中 在在 连续且连续且 无妨设无妨设 在在 内仅有一个零点。内仅有一个零点。 求方程（求方程（ ）的实根）的实根 的二分法过程，就是将的二分法过程，就是将 逐步分半，检

14、查函数值符逐步分半，检查函数值符号的变化，以便确定包含根的充分小区间。号的变化，以便确定包含根的充分小区间。 0 f x 4.2.1 f x,a b 0f a f b f x,a b4.2.1x,a b15；.二分法的步骤如下：记二分法的步骤如下：记 ， 第第1 1步：分半计算步：分半计算 ，将，将 分半。计算中点分半。计算中点 及及 。若。若 ，则根必在，则根必在 内，否则必在内，否则必在 内，（若内，（若 ，则，则 ），于是得到长），于是得到长度一半的区间度一半的区间 含根含根, ,即即 , ,且且 。 第第 步（步（* *分半计算）重复上述过程。分半计算）重复上述过程。1aa1bb2x1

15、1ba1122, ,a xa b1122, ,xbab22111()2baba1k11a,b 1f x11( )( ) 0f af x1( ) 0f x 1xx22,a b22() ( )0f af bk16；.设已完成第设已完成第1 1步步第第 步，分半计算得到含根区步，分半计算得到含根区间间 ，且满足，且满足 , ,即即 ， 即即 ， 则第则第k k步的分半计算：步的分半计算： ，且有：，且有：122,kkababab , kkxa b11()2kkkbaba2kkkabx1 22kkkkbaxxba 4 .2 .21k( ) ( ) 0kkf a f b 17；.确定新的含根区间确定新的

16、含根区间 ，即如果，即如果 ，则根必在，则根必在 内，否则必在内，否则必在 内，且有：内，且有： 。总之，由上述二分法得到序列。总之，由上述二分法得到序列 ，由，由(4.2.2)(4.2.2)有：有： 。可用二分法求方程可用二分法求方程 的实根的实根 的近似值到任意指定的精度，这是因为：设的近似值到任意指定的精度，这是因为：设 为给定精度要求，则由为给定精度要求，则由 ，可得分半计算次数，可得分半计算次数k k应满足：应满足：11,kkab( ) ( ) 0kkf a f x 11,kkkkabax11,kkkkabxb111()2kkkbaba kxlimkkxx0( ) 0f x x2kk

17、baxx18；. 二分法的优点是方法简单，且只要求二分法的优点是方法简单，且只要求 连续即可，可用二分法求出连续即可，可用二分法求出 在在 内全部实根，但二分法不能求复根及偶数重根，且收敛较慢，函数值内全部实根，但二分法不能求复根及偶数重根，且收敛较慢，函数值计算次数较多。计算次数较多。 lnlnln 2bak 4.2.3( )f x()0fx,ab19；.例例4.2.2 4.2.2 用二分法求用二分法求 在在1,21,2内一个实根，且要求精确到小数内一个实根，且要求精确到小数点后第三位。（即点后第三位。（即 ） 解解 由由 代入公式代入公式(4.2.3) (4.2.3) ，可确定所需分半次数

18、为，可确定所需分半次数为 ，计算结果部分如下表：（显然计算结果部分如下表：（显然 ）。）。6( )1f xxx*3121 0kxx30.5101,2)ab（11k(1)10,(2)0ff 20；.K81.132813 1.140625 1.136719 0.020619 91.132813 1.136719 1.134766 0.4268415 101.132813 1.134766 1.133789 111.133789 1.134766 1.134277 表表4.2.2 4.2.2 部分计算结果部分计算结果kakbkx()kf x00959799. 00045915. 021；.4.2.4

19、 4.2.4 迭代法迭代法 迭代法是一种逐次逼近法。它是求解代数方程，超越方程及方程组的一种基本迭代法是一种逐次逼近法。它是求解代数方程，超越方程及方程组的一种基本方法，但存在收敛性及收敛快慢的问题。方法，但存在收敛性及收敛快慢的问题。 用迭代法求解用迭代法求解 的近似根，首先需将此方程化为等价的方程：的近似根，首先需将此方程化为等价的方程： 然而将然而将 化为等价方程化为等价方程 的方法是很多的。的方法是很多的。 ( ) 0f x (4.2.4)( )0f x ( ) x g x (4.2.4)22；. 例例4.2.3 4.2.3 对方程对方程 可用不同的方法将其化为等价方程：可用不同的方法

20、将其化为等价方程： （1 1） （2 2）()sin0.50fxxx1sin0.5( )xxg x12sin0.5( )xxg x23；.定义定义4.2.2 （迭代法）（迭代法）设方程为设方程为 取方程根的一个初始近似取方程根的一个初始近似 ，且按下述逐，且按下述逐次代入法，构造一个近似解序列：次代入法，构造一个近似解序列： 这种方法称为迭代法（或称为单点迭代法），这种方法称为迭代法（或称为单点迭代法）， 称为迭代函数。称为迭代函数。若由迭代法产生序列若由迭代法产生序列 有极限存在，即有极限存在，即 ，称，称 为收敛或迭代过程为收敛或迭代过程 收敛，否则称迭代法不收敛。若收敛，否则称迭代法不收

21、敛。若 连续，且连续，且 ，则则 ， 即即 为方程为方程 的解（称的解（称 为函数为函数 的不动点），显然在由方程的不动点），显然在由方程 转化为等价方程转化为等价方程 时，选择不同的迭代函数时，选择不同的迭代函数 ，就会产生不同的序列，就会产生不同的序列 （即使初值（即使初值 选择一样）且这些序列的收敛情况也不会相同。选择一样）且这些序列的收敛情况也不会相同。( )xg x0 x 10211, kkxg xxg xxg x (4.2.5)( )g x kx*limkkxxkx (4.2.5)( )g x*limkkxx1limlim ()lim()kkkkkkxxg xgxg x (4.2.

22、4)x g xx( )0f x ( )xg x( )g x kx0 x24；.例例4.2.4 4.2.4 对例对例4.2.14.2.1中方程考查用迭代法求根中方程考查用迭代法求根 由计算可以看出，我们选取的两个函数由计算可以看出，我们选取的两个函数 ，分别构造序列，分别构造序列 收敛情形不收敛情形不一样（初值都取为一样（初值都取为1 1），在），在 中中 收敛且收敛且 ，在，在 中计算出中计算出 无定义。无定义。 111sin0.5,0,1,2,sin0.5 ,0,1,2,kkkkaxxkbxxk 12,g x g x kx( )a kx1.497300 x( )b114sin0.5sin1.

23、987761x25；.0 1.0 1.0 11.341471 0.523599 21.473820 0.023601 31.049530 -0.496555 41.497152 -1.487761 51.497289 61.497300 71.497300 k( )kax( ) kbx ( )kaf x73.6 10表表4.2.3 4.2.3 部分计算结果部分计算结果26；. 因此对用迭代法求方程因此对用迭代法求方程 的近似根，需要研究下述问题：的近似根，需要研究下述问题： （1 1）如何选取迭代函数）如何选取迭代函数 使迭代过程使迭代过程 收敛。收敛。 （2 2）若）若 收敛较慢时，怎样加速

24、收敛较慢时，怎样加速 收敛。收敛。 ( )0f x ( )g x 1kkxg x kx kx27；.迭代法的几何意义：迭代法的几何意义： 从几何意义看，求方程从几何意义看，求方程 根的问题，是求曲线根的问题，是求曲线 与直线与直线 交点的横坐标交点的横坐标 ，当迭代函数，当迭代函数 的导数函数的导数函数 在根在根 处满足下述几种条件处满足下述几种条件时，从几何上来看迭代过程时，从几何上来看迭代过程 的收敛情况如的收敛情况如图图4.2.14.2.1。 从曲线从曲线 上一点上一点 出发，沿着平行于出发，沿着平行于x x轴方向前进交轴方向前进交 于于一点一点 再从点再从点 沿平行于沿平行于y y轴方

25、向前进交轴方向前进交 于于 点，显然点，显然 的横坐标就的横坐标就是是 ，继续这过程就得到序列，继续这过程就得到序列 ，且从几何上观察知道在（，且从几何上观察知道在（1 1），（），（2 2）情况下情况下 收敛于收敛于 ，在（，在（3 3），（），（4 4）情况）情况 不收敛于不收敛于 。 ( )xg x( )yg xyxx( )g x xgx1kkxg x( )yg x000,px g x0Q0Q( )y gx1p1p10 xg xkx kx*x kx*xyx28；.图图4.2.1 4.2.1 迭代法的几何意义图迭代法的几何意义图29；. 由迭代法的几何定义知，为了保证迭代过程收敛，应该要求

26、迭代函数的导数满足由迭代法的几何定义知，为了保证迭代过程收敛，应该要求迭代函数的导数满足条件条件 。当。当 时，原方程在时，原方程在 中可能有几个根或迭代法不收敛，为此中可能有几个根或迭代法不收敛，为此有关于迭代收敛性的定理有关于迭代收敛性的定理4.2.34.2.3。 定理定理4.2.3 4.2.3 设有方程设有方程 ， (1) (1) 设设 于于 一阶导数存在，一阶导数存在， (2) (2) 当当 时，有时，有 ， ( (3 3) ) 满足条件：满足条件： 则有：则有： 在在 上有唯一解上有唯一解 ， 对任意选取初始值对任意选取初始值 ，迭代过程，迭代过程 收敛即收敛即 ， 误差估计误差估计

27、()1gx , xab , a bg( )xxg( ) xa,ba,bxg( ) a,bx ( )1, , ,g xLxa b g( ) x1g( )xx , a b*x3*11 1kkkxxxxL4*10 ,(1,2,.)1kkLxxxxkL20 , xab*klim xx1(),k0,1,.kkxg x30；.证明证明 只证只证 ， ， 由定理条件由定理条件 ，当取，当取 时，则有时，则有 记误记误差差 ，由中值定理有：，由中值定理有： ，其中，其中 在在 与与 之间，之间，即即 ，又由条件有：，又由条件有： ，由此递推可，由此递推可得：得： ，由，由 故故 。 由迭代公式由迭代公式 有：

28、有： ，其中，其中c c在在 与与 之间，于是：之间，于是： 即即 。 234*111*() -L x(1-L) kkkkkkkkkxxxxxxxxxxxxxxx*11111kkkkkLxxxxxxLL*klimxx20 , xa b , ,(1,2,.)kxa bk*kkexx*1()()( )()kkkxxg xg xg c xxc*xkx , ca b*2*120.kkkkxxL xxL xxL xx01L*1( )kkkxxg c xxL xx321()kkxg x1111( )()( )()kkkkkkkkxxg xg xg c xxLxx1kxkx331；. 由上面由上面 反复利用

29、代入上式中有反复利用代入上式中有： 由定理由定理 结果可知，当计算得到的相邻两次迭代满足条件结果可知，当计算得到的相邻两次迭代满足条件 时，则误时，则误差差 。 因此在计算机上可利用因此在计算机上可利用 来控制算法终止，但要注意来控制算法终止，但要注意 时，时，即使即使 很小，误差很小，误差 仍然可能很大。仍然可能很大。 另外，当已知另外，当已知 及及 及给定精度要求及给定精度要求 时，利用定理时，利用定理 结果可确定使结果可确定使误差达到给定精度要求时所需要迭代次数误差达到给定精度要求时所需要迭代次数k k，事实上，由，事实上，由 解得：解得： 411kkkkxxLxx*1121210111

30、L .11kkkkkkkkLxxxxxxLLLxxxxLL31kkxx*11kxxL*kxx1kkxx1L 1kkxx01,x x(1)L L4*101kkLxxxxL10(lnln)/ln 1xxkLL32；. 定理条件定理条件 ，在一般情况下，可能对大范围的含根区间不满足，而在，在一般情况下，可能对大范围的含根区间不满足，而在根的邻近是成立的，为此有如下迭代过程的局部收敛性结果。根的邻近是成立的，为此有如下迭代过程的局部收敛性结果。 定理定理4.2.4 4.2.4 （迭代法的局部收敛性）设给定方程（迭代法的局部收敛性）设给定方程 （1 1）设）设 为方程的解，为方程的解， （2 2）设）设

31、 在在 的邻域内连续可微，且有的邻域内连续可微，且有 ，则对任意初值，则对任意初值 （在（在 的邻的邻域内），迭代过程域内），迭代过程 ， 收敛于收敛于 。( )1, , g xLxa bg( )xx*xg( ) x*x*()1gx0 x*x1( )kxg x0,1,.k *x33；.例例4.2.5 4.2.5 由迭代法解方程由迭代法解方程( )ln(2)0f xxx34；.解解 (1)(1)显然有显然有 即知方程于即知方程于0,20,2及及-1.9,-1-1.9,-1内有根记内有根记为为 。 (2)(2)考察取初值考察取初值 迭代过程迭代过程 的收敛性，其中迭代函数的收敛性，其中迭代函数为为

32、 ，显然，显然 ， ，及，及 为增函数，则当为增函数，则当 时，时， ，又由，又由 则有则有 。 于是由定理于是由定理4.2.44.2.4可知，当初值可知，当初值 时，迭代过程时，迭代过程 收敛，如果收敛，如果要求要求 的近似根准确到小数点后第的近似根准确到小数点后第6 6位（即要求位（即要求 ）由计算结果可）由计算结果可知知 。且。且 ，则，则 ， 。 (0) (2)0, ( 1.9) ( 1)0ffff 12,xx00,2x 1ln(2)kkxx1( ) ln(2)g xx1(0) ln(2) 0.693 0g1(2) ln(4) 1.368 2g1( )g x02x10( )2g x1(

33、)2gxx111( )(0)1(0,2)22g xgxx 00,2x 1ln(2)kkxx1x611102kx x 7151410 xx12L *1141.1461931xx714()0.8 10f x35；.表表4.2.4 4.2.4 部分计算结果表部分计算结果表k1ln(2)kkxx00.010.6931471820.99071046141.1461931151.146193236；. (3)(3)为了求为了求-1.9,-1-1.9,-1内方程的根。由迭代方程内方程的根。由迭代方程 ，显，显然然 ，所以迭代过程，所以迭代过程 （初（初值值 ）不能保证收敛于）不能保证收敛于 。 (4)(4)

34、若将方程转化为等价方程若将方程转化为等价方程 或或 则则 ，且且 （ 时）时） 所以当选取所以当选取 时迭代过程时迭代过程 收敛。如取收敛。如取 ，则迭代，则迭代1212次有次有 ，且且 。 由此例可见，对于方程由此例可见，对于方程 ，迭代函数，迭代函数 取不同形式，相应的迭代法产生取不同形式，相应的迭代法产生的的 收敛情况也不一样，因此，我们应该选择迭代函数时构造的迭代过程收敛情况也不一样，因此，我们应该选择迭代函数时构造的迭代过程 收敛，且收敛较快。收敛，且收敛较快。1ln(2)kkxx*121( )()12gxgxx1ln(2)kkxx*002 1.9, 1,xxx *2xe2xx 22

35、( )xxegx2g ( )xxe22g ( )g ( 1)0.386 1x 1.9, 1x2( ) 1.9, 1g x 0 1.9, 1x k 12kxxe01x *2121.841405660 xx 812()0.2 10f x ( ) 0f x ( )g x kx1()kkxg x37；.关于迭代公式的加工：关于迭代公式的加工： 对于收敛的迭代过程，只要迭代足够多次，总可以使结果达到任意的精度。但对于收敛的迭代过程，只要迭代足够多次，总可以使结果达到任意的精度。但有时迭代收敛缓慢，从而使计算量变得很大，因此迭代过程的加速是一个很重要的有时迭代收敛缓慢，从而使计算量变得很大，因此迭代过程的

36、加速是一个很重要的课题。课题。 设设 为根为根 的某个预测值，用迭代公式校正一次得：的某个预测值，用迭代公式校正一次得： 由中值定理：由中值定理： ， 介于介于 之间，若之间，若 改变不大。近似地改变不大。近似地取某常数，则由取某常数，则由 10()xg x*10( )()x xgxx *0 xx，( )g x*10101()11Lx xLxxxxxLL 0 x*x 可以期望按上式右端求得的可以期望按上式右端求得的 是比是比更好的近似值。更好的近似值。 若将每得到一次改进值算作一步，并用若将每得到一次改进值算作一步，并用 和和 分别表示第分别表示第 步的校正值和改步的校正值和改进值，则加速迭代

37、计算方案如下：进值，则加速迭代计算方案如下： 校正：校正： 改进：改进： 由于使用参数由于使用参数 ，这在实际应用中不方便，下面进行改进计算。，这在实际应用中不方便，下面进行改进计算。2101101()111LLxxxxxxLLLxkxkk1 ()kxg xk 111k ( )1kkLxxxxLL38；. 设设 的某近似值的某近似值 ，将校正值，将校正值 再校正一次得：再校正一次得： ，由，由 与与 得：得： 由此得：由此得： 。这样将上式右端作为改进公式就不再含有。这样将上式右端作为改进公式就不再含有导数信息了。但需要用到两次迭代的结果进行加工。如果仍将得到一次改进值作为导数信息了。但需要用

38、到两次迭代的结果进行加工。如果仍将得到一次改进值作为一步，则计算过程如下：一步，则计算过程如下： 上述处理过程称为（埃特金）方法。上述处理过程称为（埃特金）方法。*x0 x10()xg x21( )xg x*21()x x Lx x *10()()x xLx x2*02121201201222x xxxxxxxxxxxx2（）*01*21xxxxxxxxk11k1k1k1k12k1k1k1 () () 2kkkxgxxgxxxxxxxx校 正 ：再 校 正 ：（）改 进 ：39；.4.2.5 Newton4.2.5 Newton公式公式 对于方程对于方程 ，应用迭代法时先要改写成，应用迭代法时

39、先要改写成 ，即需要针对，即需要针对 构构造不同的合适的迭代函数造不同的合适的迭代函数 ，显然可以取迭代函数为，显然可以取迭代函数为 ，相应迭代，相应迭代公式为公式为 。 一般地，这种迭代公式不一定收敛，或者速度很慢。对此公式应用前面的加速一般地，这种迭代公式不一定收敛，或者速度很慢。对此公式应用前面的加速技术具体格式为：技术具体格式为： ( )0f x ( )xg x( )f x( )g x( )( )g xxf x 1()kkkxxf x1111()()1kkkkkkkxxf xLxxxxL40；. 记记 ，则上二式可合并写为：，则上二式可合并写为： 。此公式称为简单的。此公式称为简单的N

40、ewtonNewton公式，公式，其迭代函数为：其迭代函数为： 。又由于。又由于 为为 的近似值，而的近似值，而 ，因此，因此 实际上是实际上是 的近似值，故用的近似值，故用 代替上式中的代替上式中的 即得到下面的迭代函数：即得到下面的迭代函数： 。 相应的迭代公式为：相应的迭代公式为： ， 即为即为NewtonNewton公式。公式。1ML1()kkkf xxxM()()fxg xxML( )g x()()g xxfx1ML( )f x( )fxM( )( )( )f xg xxfx1()()kkkkf xxxfx41；.4.2.6 Newton4.2.6 Newton法的几何意义法的几何意

41、义 NewtonNewton法的基本思想就是将非线性方程法的基本思想就是将非线性方程 逐步线性化求解，设逐步线性化求解，设 有有近似的根近似的根 ，将，将 在在 处处 展开得：展开得： ，从而，从而 近似地表为：近似地表为： 。方程。方程 的根的根 即为曲线即为曲线 与与 轴焦点的横坐标。设轴焦点的横坐标。设 为为 的的 一个近似，过曲线一个近似，过曲线 上横坐标上横坐标 为为 的点的点 作曲线的切线，该切线作曲线的切线，该切线 与与 轴焦点的横坐标即为轴焦点的横坐标即为 的新近似的新近似 值值 ，它与，它与 轴交点的横坐标为：轴交点的横坐标为： ，因此，因此 NewtonNewton法亦称切

42、线法。法亦称切线法。( )0f x ( )0f x kx( )f xkxTaylor( )( )( )()kkkf xf xf x x x( )0f x ()()()0kkkf xfxxx( ) 0f x *x( )yf xxkx*x( )yf xkxkpx*x1kxx1( )/( )kkkkxxf xf x42；.4.2.7 Newton4.2.7 Newton法的局部收敛性法的局部收敛性定义定义4.2.34.2.3 设迭代过程设迭代过程 收敛于方程收敛于方程 的根的根 ，如果迭代误，如果迭代误差差 ，当，当 时有：时有： 则称该迭代过程为则称该迭代过程为 阶收敛的。阶收敛的。定理定理4.2

43、.54.2.5 对迭代过程对迭代过程 如果如果 在在 附近连续，且：附近连续，且： 且且 ，则该迭代过程在，则该迭代过程在 附近是附近是 阶收敛的。阶收敛的。1()kkxg x( )x gx*x*kkexxk 1 (0,)kpkecce为常数pk 1( ),kxg x( )( )pgx*x*(p-1)*g( ) g( ) . g( ) 0 xxx (p)*g ( ) 0 x xp43；. 证明证明 由于由于 ，则有前面关于迭代法的局部收敛性定理知：此迭代过程，则有前面关于迭代法的局部收敛性定理知：此迭代过程 具有局部收敛性。即具有局部收敛性。即 。将。将 在在 处处 展开，并注意到展开，并注意

44、到 有：有： 而，而， 从而上式化为：从而上式化为： *g( )0 1x 1( )kkxg x0kx ( )kg x*xTaylor*(p-1)*g( ). g( )0 xx( )*()()()(), , !pkkkkkgg xg xxxx xp*k 1( ), ( )kxg xxg x()*k1()()!ppkkgxxxxp44；. 即：即： 故知迭代过程具有故知迭代过程具有 阶收敛性。阶收敛性。 *()()*1k1p*k()e()e()!ppkkpkxxggxxxppp 定理定理4.2.54.2.5 表明迭代过程的收敛速度依赖于迭代函数表明迭代过程的收敛速度依赖于迭代函数 的选取，如果的选

45、取，如果 时时 。则迭代过程只可能是线性收敛的。则迭代过程只可能是线性收敛的。( )g x , xa b( ) 0g x 对于对于NewtonNewton法，由迭代函数为：法，由迭代函数为： 则则 ， 若若 为为 的一个单根。即的一个单根。即 ，则由上式知，则由上式知 。 由上面定理可知由上面定理可知NewtonNewton法在根法在根 的邻域内是平方收敛的（二阶收敛的邻域内是平方收敛的（二阶收敛 的）。的）。() ()-()fxgxxfx222( )( )( )( )( )g ( )1( )( )fxf x fxf x fxxfxfx*x( )f x *()00f xfx，*( ) 0g x

46、 *x45；. 特别地考察特别地考察NewtonNewton公式：设公式：设 二次连续可微，二次连续可微，则则 ， 在在 之间，特别地取之间，特别地取 ，注，注意意 ，则，则 设设 。两边同除以。两边同除以 ，得：，得： （注：（注： ），利用），利用NewtonNewton公式，即有：公式，即有： 当当 ，则，则 ， 或或( )f x2( )( )()()()()2!kkkkff xf xfxxxxx,kx x*xx*( ) 0f x *2( )0( )( )( )()()2!kkkkff xf xf xxxxx( )0kf x( )kf x*2()()()()2()kkkkkfxfxxxx

47、fxfx1()()kkkkf xxxfx*1*2()()2()kkkkxxfmxxfx k *( )()2()2()kffxfxfx *2211( )()2 ( )kkkk kkfxxexxmef x46；. 可见可见 （误差）与（误差）与 的误差的误差 的平方成比例。当初始误差的平方成比例。当初始误差 充分小时，以后迭代的误差将减少得非常快。反之充分小时，以后迭代的误差将减少得非常快。反之 ，则放大。，则放大。NewtonNewton法法每计算一步，需要计算一次函数值每计算一步，需要计算一次函数值 和一次导数值和一次导数值 。 例例4.2.6 4.2.6 用用NewtonNewton法求解法

48、求解 。 解解 显然显然 。则在。则在0,20,2内方程有一个根，求导内方程有一个根，求导 则则NewtonNewton公式为：公式为： 取取 ，迭代，迭代6 6次得近似根为次得近似根为 , , 。这表明，当初值。这表明，当初值 取值靠近取值靠近 时，时，NewtonNewton法收敛且收敛较快，否则法收敛且收敛较快，否则NewtonNewton法可能不收敛。法可能不收敛。1kekxke*00exx01e ()kf x( )kf x4( )(2)10 xf xex(0) (2)0ff4( )(6) / 4xfxex414()(2)1()(6) / 4kkxkkkkkxkkf xexxxxfxe

49、x01.0 x *0.783596x *6( )3.8 10f x0 x*x47；. 下面考虑下面考虑NewtonNewton法的误差估计，由中值定理有：法的误差估计，由中值定理有： ，当当 充分接近充分接近 时，有时，有 因此，用因此，用NewtonNewton法求方程单根法求方程单根 的近似根的近似根 的误差的误差 可用可用 来估计。来估计。 *( )( )( )( )()kkkkf xf xf xfxxkx*1( )( )( )( )Newtonkkkkkkkf xf xxxxxff x*x*xkx*kkexx 1kkxx48；.4.2.8 Newton4.2.8 Newton法应用举例

50、法应用举例1. 1. 对给定的正数对给定的正数 ，应用，应用NewtonNewton法解二次方程法解二次方程 可导出求开方值可导出求开方值 的计的计算格式：算格式： 可证明公式可证明公式 对任意函数初值对任意函数初值 都是收敛的。这是因为：都是收敛的。这是因为： C20 xCC11() 2kkkCxxx (4.2.7)(4.2.7)00 x 21211()21()2kkkkkkxCxCxxCxCx49；. 两式相除得：两式相除得： 利用此式递推可得：利用此式递推可得： （ 由由 可知：可知： ，则：，则： ）而）而 ， 故由公式知故由公式知 即迭代法恒收敛。）即迭代法恒收敛。）211()kkk

51、kxCxCxCxC222102102()1kkkkkkkxCxCqqxCCxCxCq21010()kkkxCxCxCxC2211kkkqxCq2222()1kkkkkqxCxC qCq00 x 001xCqxC()kxC k50；. 例例4.2.7 4.2.7 求求 的近似值，要求的近似值，要求 终止迭代。终止迭代。 解解 取取 经经6 6次迭代后：次迭代后： ， ， ，故，故 。 对给定正数对给定正数 ，应用，应用NewtonNewton法求解法求解 ，由此式可导出求，由此式可导出求 而不用除法的计而不用除法的计算程序：算程序： 。 这个算法对于没有设置除法操作的电子计算机是有用的。可以证明

52、，此算法初值满足这个算法对于没有设置除法操作的电子计算机是有用的。可以证明，此算法初值满足 时是收敛的，这是因为：时是收敛的，这是因为： 即：即： ，令，令 ，有递推公式：，有递推公式： ，反复递推，反复递推得：得： 。 当当 ，即，即 时，有时，有 即即 ，从而迭代法收敛。，从而迭代法收敛。 106110kkxx1110()2kkkxxx01.0 x 53.16227767x 63.16227766x 7650.1 10 xx103.16227766C10Cx1C1(2)kkkxxCx020 xC21111(2)()kkkkxxCxC xCCC211(1)kkCxCx1kkrC x21kkr

53、r20kkrr0011rCx020 xC0kr 1kxC51；.4.2.9 Newton4.2.9 Newton下山法下山法 NewtonNewton法收敛性依赖于法收敛性依赖于 初值的选取，如果初值的选取，如果 偏离偏离 较远，则较远，则NewtonNewton法可能法可能发散。发散。 例如，对方程例如，对方程 。求在。求在 附近的一个根附近的一个根 。若取初值。若取初值 ，则由则由NewtonNewton法：法： 计算得计算得 ，仅迭代，仅迭代3 3次即得有次即得有6 6位有效数字的近位有效数字的近似值似值 。但若取初值。但若取初值 则由同一则由同一NewtonNewton公式计算得公式计

54、算得 ，这反而比，这反而比 更远离所求根更远离所求根 ，因此发散。为防止发散，对迭代过程加一下降要求：，因此发散。为防止发散，对迭代过程加一下降要求： 满足这项要求的算法称为下山法。满足这项要求的算法称为下山法。0 x0 x*x31 0 xx 1.5x *x01.5x 312131kkkkkxxxxx1231.34783,1.32520,1.32472xxx3x00.6x 117.9x 00.6x*1.32472x 1()()kkfxfx52；. 将将NewtonNewton法与下山法结合，即在下山法保证函数下降条件下，用法与下山法结合，即在下山法保证函数下降条件下，用NewtonNewton

55、法加速收敛。法加速收敛。为此，可将为此，可将NewtonNewton计计 算结果算结果 与每一步近似值与每一步近似值 作加权均：作加权均： ，其中，其中 （ ）称为下山）称为下山 因子。选择下山因子因子。选择下山因子 以保证下降性。以保证下降性。 的选择方法是：由的选择方法是：由 反复减半的试探法，若能找到反复减半的试探法，若能找到 使下降性成立，则下山使下降性成立，则下山成功，否则下山失败，改变初值成功，否则下山失败，改变初值 重新开始。重新开始。 1()()kkkkf xxxfxkx11(1)kkkxxx0110 x53；.4.2.10 4.2.10 弦截法与拋物法弦截法与拋物法 Newt

56、onNewton法法 每迭代一次计算函数值每迭代一次计算函数值 ，导数值，导数值 各一次，各一次，当当 函数本身比较复杂时，求导数值更加困难。函数本身比较复杂时，求导数值更加困难。 下面方法多利用以前各次计算的函数值下面方法多利用以前各次计算的函数值 来回避导数值来回避导数值 的计算，的计算，导出这种求根方法的基本原理是插值法。导出这种求根方法的基本原理是插值法。 设设 是是 的一组近似值，利用对应的函数的一组近似值，利用对应的函数值值 ，构造插值多项式，构造插值多项式 ，适当选取，适当选取 的一个根作为的一个根作为 的新的近似根的新的近似根 。这样就确定了一个迭代过程，记迭代函数。这样就确定

57、了一个迭代过程，记迭代函数为为 ，则，则 ，下面具体考察，下面具体考察 （弦截法），（弦截法）， （拋物（拋物法）两种情形。法）两种情形。 1()()kkkkf xxxfx( )kf x( )kf xf1( ), (),kkf xf x()kfx1,kkk rx xx( ) 0f x 1( ), (), ()kkk rf xf xf x( )rp x( ) 0rp x ( )0f x 1kxg11( ,)kkkk rxg x xx1r 2r 54；.4.2.11 4.2.11 弦截法弦截法 设设 为为 的近似根，过点的近似根，过点 ， 构造一次插值多项构造一次插值多项式式 ，并用，并用 的根作

58、为的根作为 的新的近似根的新的近似根 。由于。由于 则由则由 可得：可得： 另外，公式另外，公式(4.2.9)(4.2.9)也可以用导数也可以用导数 的差商的差商 近似取代近似取代NewtonNewton公公式中的式中的 ，同样得公式，同样得公式 。1,kkxx( ) 0f x ( , ( )kkx f x11(, ()kkxf x1( )p x1( ) 0p x ( )0f x 1kx111()()( )()() kkkkkkf xf xp xf xxxxx(4.2.8)1( )0p x 111()() ()()kkkkkkkf xxxxxf xf x (4.2.9)( )f x11()()

59、kkkkf xf xxx( )fx (4.2.9)55；.v弦截法的几何意义：弦截法的几何意义： 如图，曲线如图，曲线 上横坐标为上横坐标为 的点分别记为的点分别记为 ，则弦线，则弦线 的斜的斜率等于差商率等于差商 。 的方程为：的方程为： 则按则按 求得的近似根求得的近似根 实际上是弦实际上是弦线线 与与 轴交点的横坐标。因此这种轴交点的横坐标。因此这种算法称为弦截法，又称割线法。算法称为弦截法，又称割线法。 ( )yf x1,kkxx1,kkpp1kkp p11()()kkkkf xf xxx1kkp p11( )()( )()0kkkkkf xf xf xx xxx (4.2.9)1kx

60、1kkp px56；. 弦截法与切线法（弦截法与切线法（NewtonNewton法）都是线性化方程，但两者有本质区别。法）都是线性化方程，但两者有本质区别。NewtonNewton切线法切线法在计算在计算 时只用到前一步的时只用到前一步的 及及 ，但要计算，但要计算 ，而弦截法在计算，而弦截法在计算 时要用前面两步的结果时要用前面两步的结果 ，而不须计算导数。这种方法必须有两个启动值，而不须计算导数。这种方法必须有两个启动值 。 例例4.2.8 4.2.8 用割线法求解方程用割线法求解方程 在在 的根。的根。 解解 取初值取初值 ，则迭代，则迭代5 5次后有次后有 ， 。例子表明弦截法仍具有较