《第八章非线性方程组求根》由会员分享,可在线阅读,更多相关《第八章非线性方程组求根(71页珍藏版)》请在金锄头文库上搜索。
1、第八章第八章非线性非线性方程(组)求根方程(组)求根问题驱动:全球定位系统(问题驱动:全球定位系统(GPS) 人类对导航和定位的需求是伴随着人类整个文明历史的进人类对导航和定位的需求是伴随着人类整个文明历史的进步而发展的,中国古代步而发展的,中国古代“四大发明四大发明”之一的指南针是最早的定之一的指南针是最早的定位仪器和系统,其后还有经纬仪以及近代的雷达。如图位仪器和系统,其后还有经纬仪以及近代的雷达。如图8.1.1所所示全球定位系统(示全球定位系统(GPS)是基于卫星的导航系统,最早由美国)是基于卫星的导航系统,最早由美国和前苏联分别在和前苏联分别在80年代研制,并于年代研制,并于1993年
2、正式投入使用。现年正式投入使用。现代社会中全球定位系统越来越深入到人们生活的方方面面。例代社会中全球定位系统越来越深入到人们生活的方方面面。例如市场上出售的手持型如市场上出售的手持型GPS,定位的精度可以达到,定位的精度可以达到10米以内,米以内,这无疑给旅行者提供了方便;安装有这无疑给旅行者提供了方便;安装有GPS的儿童手表,家长在的儿童手表,家长在家里的计算机上可以追踪到孩子的位置,防止儿童走家里的计算机上可以追踪到孩子的位置,防止儿童走失;安装有失;安装有GPS系统的汽车可以帮助新司机辨识道路等等。系统的汽车可以帮助新司机辨识道路等等。 图图8.1.1 卫星定位示意图卫星定位示意图 美国
3、和前苏联的美国和前苏联的GPS都包括有都包括有24颗卫星,它们不断地向地颗卫星,它们不断地向地球发射信号报告当前位置和发出信号的时间,卫星分布如图球发射信号报告当前位置和发出信号的时间,卫星分布如图8.1.2所示。它的基本原理是:在地球的任何一个位置,至少所示。它的基本原理是:在地球的任何一个位置,至少同时收到同时收到4颗以上卫星发射的信号。颗以上卫星发射的信号。发射的信号,发射的信号, 设地球上一个点设地球上一个点R R,同时收到卫星,同时收到卫星 假设接收的信息如表假设接收的信息如表8.1.1所示。请设法确定所示。请设法确定R点的位置。点的位置。 图图8.1.2 卫星分布图卫星分布图表表8
4、.1.1GPS导航问题可归结为求解非线性代数数方程组导航问题可归结为求解非线性代数数方程组 , 当当 时就是单个方程时就是单个方程. (8.1.1)。.其中,其中, 可以是代数方程,也可以是超越方程。使可以是代数方程,也可以是超越方程。使 成立的成立的x 值称为方程的根,或称为值称为方程的根,或称为 计算中,如电路和电力系统计算、非线性力学、非线性微(计算中,如电路和电力系统计算、非线性力学、非线性微(积分)方程、非线性规划(优化)等众多领域中,问题的求积分)方程、非线性规划(优化)等众多领域中,问题的求解和模拟最终往往都要解决求根或优化问题。前一种情形要解和模拟最终往往都要解决求根或优化问题
5、。前一种情形要求出方程(组)的根;后一种情形则要求找出函数取最大或求出方程(组)的根;后一种情形则要求找出函数取最大或最小的点。即使是对实验数据进行拟合或数值求解微分方程,最小的点。即使是对实验数据进行拟合或数值求解微分方程,也总是将问题简化成上述两类问题。上述除少数特殊方程外,也总是将问题简化成上述两类问题。上述除少数特殊方程外,大多数非线性代数方程(组)很难使用解析法求解精确解,大多数非线性代数方程(组)很难使用解析法求解精确解,一般需要通过一些数值方法逼近方程的解。这里主要介绍单一般需要通过一些数值方法逼近方程的解。这里主要介绍单个方程的数值解法,方程组也可以采用类似的方法,将放在个方程
6、的数值解法,方程组也可以采用类似的方法,将放在后面讨论。后面讨论。的零点。科学与工程的零点。科学与工程1根的存在性。方程有没有根?如果有,有几个根?根的存在性。方程有没有根?如果有,有几个根?2根的搜索。根的搜索。这些根大致在哪里?如何把根隔离开?这些根大致在哪里?如何把根隔离开?3根的精确化。根的精确化。 f (x) = 0 (8.2.1)1.1.根的存在性根的存在性定理定理1:设函数设函数 f (x) 在区间在区间a, b上连续上连续,如果如果f (a) f (b) 0打打 印印结结 束束否否是是继续扫描继续扫描例例1 1:考察方程:考察方程x00.51.01.5f (x) 的的符符号号a
7、bx1x2ab或或不能保证不能保证 x 的精的精度度x* 2xx*1 1 1 1 二二二二 分分分分 法法法法执行步骤执行步骤1计算计算f (x)在有解区间在有解区间a, b端点处的值端点处的值,f (a),f (b)。2计算计算f (x)在区间中点处的值在区间中点处的值f (x1)。3判断若判断若f (x1) = 0,则则x1即是根,否则检验即是根,否则检验:(1)若若f (x1)与与f (a)异号异号,则知解位于区间则知解位于区间a, x1, b1=x1, a1=a;(2)若若f (x1)与与f (a)同号同号,则知解位于区间则知解位于区间x1, b, a1=x1, b1=b。反复执行步骤
8、反复执行步骤2、3,便可得到一系列有根区间便可得到一系列有根区间: 4、当当时,停止;时,停止;即为即为根的近似。根的近似。当当 时,时, ,即这些区间必将收缩于一点,也就是,即这些区间必将收缩于一点,也就是方程的根。在实际计算中,只要方程的根。在实际计算中,只要 的区间长度小于预定容的区间长度小于预定容许误差许误差就可以停止搜索,即就可以停止搜索,即 然后取其中点然后取其中点 作为方程的一个根的近似值。作为方程的一个根的近似值。 注:注: 例例1 证明方程证明方程 存在唯一的实根存在唯一的实根 用二分法用二分法求出此根,要求误差不超过求出此根,要求误差不超过 。解:记解:记 ,则对任意,则对
9、任意 ,因而,因而, 是严格单调的,是严格单调的, 最多有一个根,最多有一个根,所以,所以, 有唯一实根有唯一实根 又因为又因为 用二分法求解,要使用二分法求解,要使 ,只要,只要 解得解得 ,取,取 。所以只要二等分。所以只要二等分7次,即可求得满次,即可求得满足精度要求的根。计算过程如表足精度要求的根。计算过程如表6.2.1所示所示 k f(ak)及符号及符号f(xk)及符号及符号f(bk)及符号及符号01234 5670()0()0()0()0.0625()0.0625()0.078125()0.0859375()0.5(+)0.25(+)0.125(+)0.0625()0.09375(
10、+)0.078125()0.0859375()1( + )0.5( + )0.25( + )0.125( + )0.125( + )0.09375( + )0.09375( + )0.09375( + )表表6.2.1所以,所以, 简单简单; 对对f (x) 要求不高要求不高(只要连续即可只要连续即可) .无法求复根及偶重根无法求复根及偶重根 收敛慢收敛慢 二分法的优缺点二分法的优缺点 问题问题 虽然二分法计算简单,能够保证收敛,但是它对于方程虽然二分法计算简单,能够保证收敛,但是它对于方程单根存在区域信息要求太高,一般情况下很难实现,并单根存在区域信息要求太高,一般情况下很难实现,并且不能求
11、偶重根、复根和虚根。在实际应用中,用来求且不能求偶重根、复根和虚根。在实际应用中,用来求解方程根的主要方法是迭代法。解方程根的主要方法是迭代法。使用迭代法求解非线性代数方程的步骤为:使用迭代法求解非线性代数方程的步骤为:(1) 迭代格式的构造;迭代格式的构造;(2) 迭代格式的收敛性分析;迭代格式的收敛性分析;(3) 迭代格式的收敛速度与误差分析。迭代格式的收敛速度与误差分析。 2 迭迭 代代 法法1 1简单迭代法简单迭代法f (x) = 0x = (x)等价变换等价变换其中其中 (x)是连续函数。方程(是连续函数。方程(8.3.1)称为不动点方程,满足)称为不动点方程,满足(8.3.1)式的
12、点称为)式的点称为不动点,不动点,这样就将求这样就将求 (8.3.1)的零点问题转的零点问题转化为求化为求的不动点问题。的不动点问题。称这种迭代格式为称这种迭代格式为不动点迭代。不动点迭代。以不动点方程为原型构造迭代格式以不动点方程为原型构造迭代格式例例3:求方程求方程的一个根的一个根.构造迭代格式构造迭代格式x1 = 0.4771x2 = 0.3939x6 = 0.3758x7 =0.3758解:解:给定初始点给定初始点xyy = xxyy = xxyy = xxyy = xx*x*x*x*y= (x)y= (x)y= (x)y= (x)x0p0x1p1 x0p0x1p1 x0p0x1p1x
13、0p0x1p1定理定理8.1:如果如果 (x)满足下列条件满足下列条件(1 1)当)当x a, b时,时, (x) a, b(2 2)当任意当任意x a, b,存在存在0 L 1,使使 则方程则方程x = (x)在在a, b上有唯一的根上有唯一的根x*, ,且对任意初值且对任意初值 x0 a, b,迭代序列,迭代序列xk+1= (xk) (k = 0, 1, )收敛于收敛于x*。(8.3.2) 注注 此处此处L L可以看成是可以看成是 在区间在区间 a,ba,b 内的上界内的上界。2迭代过程的收敛性迭代过程的收敛性3迭代法的误差估计迭代法的误差估计 在定理在定理8.18.1的条件下,简单迭代法
14、产生的迭代点列的条件下,简单迭代法产生的迭代点列有如下误差估计式:有如下误差估计式:停机准则。停机准则。(8.3.3) 求方程求方程在在内的根内的根例:例:。解:解:原方程可以等价变形为下列三个迭代格式原方程可以等价变形为下列三个迭代格式由迭代格式由迭代格式 (1) 取初值取初值得得 结果是发散结果是发散的?!的?!由迭代格式由迭代格式 (2) 取初值取初值得得 结果精确到四位有效数字,迭代到结果精确到四位有效数字,迭代到得到收敛结果。得到收敛结果。 十步才能得到十步才能得到收敛的结果!收敛的结果! 由迭代格式(由迭代格式(3) 取初值取初值得得 结果精确到四位有效数字,迭代到结果精确到四位有
15、效数字,迭代到得到收敛结果。得到收敛结果。四步就能得到四步就能得到收敛的结果了!收敛的结果了!迭代格式(迭代格式(1 1)的迭代函数为)的迭代函数为 求导得求导得 当当时时故迭代格式(故迭代格式(1 1)是发散的。)是发散的。分析分析: : 迭代格式(迭代格式(2 2)的迭代函数为)的迭代函数为 当当时时由由知知当当时,时, 所以迭代格式(所以迭代格式(2 2)是收敛的。)是收敛的。迭代格式(迭代格式(3 3)的迭代函数为)的迭代函数为当当时时 由由时,时, 知知当当所以迭代格式(所以迭代格式(3 3)也是收敛的。)也是收敛的。结论结论: : 通过以上算例可以看出对迭代函数通过以上算例可以看出
16、对迭代函数所得到的所得到的若小于若小于1 1,则收敛;且上界越小收敛速度越快。,则收敛;且上界越小收敛速度越快。求导,求导,的上界若是大于的上界若是大于1 1,则迭代格式发散;,则迭代格式发散;4迭代过程的局部收敛性迭代过程的局部收敛性定义定义8.18.1定理定理8.28.2若存在若存在 的某个邻域的某个邻域R R: , 使迭代过程使迭代过程在根在根 附近具有局部收敛性。附近具有局部收敛性。对于任意初值对于任意初值 均收敛,均收敛,则称迭代过程则称迭代过程设设 为方程为方程 的根,的根, 在在 的邻近连续且的邻近连续且 ,则迭代过程,则迭代过程 在在 邻近具有局部收敛性。邻近具有局部收敛性。5
17、迭代过程的收敛速度迭代过程的收敛速度 设由某方法确定的序列设由某方法确定的序列xk收敛于方程的根收敛于方程的根x*,如果存在正实数如果存在正实数p,使得使得(C C为非零常数)为非零常数)定义定义8.2:则称序列则称序列xk收敛于收敛于x*的收敛速度是的收敛速度是p阶的阶的,或称该方法,或称该方法具有具有p 阶阶收收敛速度敛速度。当。当p = 1时,称该方法为时,称该方法为线性(一次)收敛线性(一次)收敛;当当p = 2时,称方法为时,称方法为平方(二次)收敛平方(二次)收敛;当;当1 p 2或或C=0,p=1时,称方法为时,称方法为超线性收敛超线性收敛。 6. 6. 迭代法收敛阶的判别迭代法
18、收敛阶的判别定理定理8.38.3 如果如果 在在 附近的某个领域内有附近的某个领域内有p( )p( )阶阶 连续导数,且连续导数,且则迭代格式则迭代格式 在在 附近是附近是p p阶局部收阶局部收敛的,且有敛的,且有 7. 加速收敛技术加速收敛技术 L越小迭代法的收敛速度越快,因此,可以从越小迭代法的收敛速度越快,因此,可以从寻找较小的寻找较小的L来改进迭代格式以加快收敛速度。来改进迭代格式以加快收敛速度。思思路路(1). 松弛法松弛法引入待定参数引入待定参数 ,将将 作等价变形为作等价变形为 (8.3.4) 将方程右端记为将方程右端记为 ,则得到新的迭代格式,则得到新的迭代格式 由定理由定理2
19、知知 为了使新的迭代格式比原来迭为了使新的迭代格式比原来迭代格式收敛得更快,只要满足代格式收敛得更快,只要满足且且 越小,所获取的越小,所获取的L就越小,就越小,迭代法收敛的就越快,因此我们希望迭代法收敛的就越快,因此我们希望 。可取可取 ,若记,若记 则(则(8.3.4)式可改写为)式可改写为 称为称为松弛因子松弛因子,这种方法称为,这种方法称为松弛法松弛法。为使迭代速度加快,。为使迭代速度加快,需要边计算边调整松弛因子。由于计算松弛因子需要用到微商,需要边计算边调整松弛因子。由于计算松弛因子需要用到微商,在实际应用中不便使用,具有一定局限性。若迭代法是线性收在实际应用中不便使用,具有一定局
20、限性。若迭代法是线性收敛的,当计算敛的,当计算 不方便时,可以采用不方便时,可以采用埃特金加速公式埃特金加速公式。 (2). 埃特金加速公式埃特金加速公式设迭代法是线性收敛,由定义知设迭代法是线性收敛,由定义知成立,故当成立,故当 时有时有 由此可得由此可得 的近似值的近似值 (8.3.5) 由此获得比由此获得比 和和 更好的近似值更好的近似值 ,利用利用(8.3.5)序列序列 的方法称为的方法称为(3). Steffensen 加速法:加速法: 将将Aitken加速公式与不动点迭代相结合,可得加速公式与不动点迭代相结合,可得(8.3.6) 式构造式构造埃特金(埃特金(AitkenAitken
21、)加速方法)加速方法。利用(利用(8.3.6)式构造序列)式构造序列 的方法称为的方法称为Steffensen加速方法。加速方法。即每进行两次不动点迭代,就执行一次即每进行两次不动点迭代,就执行一次Aitken加速。加速。 例例2 试用简单迭代法和试用简单迭代法和Steffensen加速法求方程加速法求方程在在 附近的根,精确至四位有效数。附近的根,精确至四位有效数。 解:记解:记 ,简单迭代法公式为:,简单迭代法公式为: 计算得计算得kxkkxkkxk00.570.55844140.5671210.6065380.56641150.5671620.5452490.56756160.56714
22、30.57970100.56691170.5671540.56006110.56728180.5671450.57117120.5670760.56486130.56719Aitken加速公式加速公式计算得计算得所以所以, 。3 3 牛顿法牛顿法一、牛顿法的迭代公式一、牛顿法的迭代公式 考虑非线性方程考虑非线性方程 原理:原理:将非线性方程线性化将非线性方程线性化 Taylor 展开展开取取 x0 x*,将将 f (x)在在 x0 做一阶做一阶Taylor展开展开:, 在在 x0 和和 x 之间。之间。将将 (x* x0)2 看成高阶小量,则有:看成高阶小量,则有:只要只要 f C1,且每步迭
23、代都有,且每步迭代都有 , 而且而且则则 x*就是就是 f (x)的根。的根。公式(公式(9.4.19.4.1)称为)称为牛顿迭代公式。牛顿迭代公式。(8.4.1)构造迭代公式构造迭代公式x*x0x1x2xyf(x)二、牛顿法的几何意义二、牛顿法的几何意义三、牛顿法的收敛性三、牛顿法的收敛性定理定理8.4: 设设f (x)在在a, b上存在二阶连续导数且满足下列条上存在二阶连续导数且满足下列条件件:(1)f (a) f (b) 0则由则由(8.4.1)确定的牛顿迭代序列确定的牛顿迭代序列xk二阶收敛于二阶收敛于f (x)在在a, b上的唯一单根上的唯一单根x*。注:注:注:注:Newton法的
24、收敛性依赖于法的收敛性依赖于x0 的选取。的选取。x*x0 x0 x0四四. 牛顿迭代法的局部收敛性与收敛速度牛顿迭代法的局部收敛性与收敛速度 ,, ,且且 设设 在包含在包含 的一个区间二阶连续可的一个区间二阶连续可导,则导,则Newton迭代法至少二阶收敛,即迭代法至少二阶收敛,即 值得注意的是,当值得注意的是,当 充分光滑且充分光滑且 是是 的重根时,牛顿的重根时,牛顿法在法在的附近是线性收敛的。且的附近是线性收敛的。且Newton迭代法在迭代法在 上的上的收敛性依赖于初值收敛性依赖于初值的选取。即初值的选取。即初值 的选取充分靠近的选取充分靠近 时,时,一般可保证一般可保证Newton
25、迭代法收敛。迭代法收敛。 并得出了并得出了 是该方程的一个根,无人知道他用什么是该方程的一个根,无人知道他用什么方法得出的,在当时这是一个非常有名的结果,试用牛顿法求方法得出的,在当时这是一个非常有名的结果,试用牛顿法求出此结果。出此结果。 解:解: 记记则则当当 时,时, ,又,又所以所以 有唯一实根有唯一实根 ,并改写,并改写 例例3 Leonardo于于1225年研究了方程年研究了方程 用牛顿迭代格式用牛顿迭代格式所以,所以, 。五、求五、求m m重根的牛顿法重根的牛顿法1 1、迭代格式、迭代格式(8.4.2)2 2、重数、重数m m的确定的确定3 3、迭代格式(、迭代格式(9.4.2)
26、9.4.2)的收敛阶(至少的收敛阶(至少2 2阶收敛)阶收敛)由于由于Newton迭代法的收敛性依赖于初值迭代法的收敛性依赖于初值 的选取,如果的选取,如果 离方程的根离方程的根 较远,则较远,则Newton迭代法可能发散。为了防止迭迭代法可能发散。为了防止迭代发散,可以将代发散,可以将Newton迭代法与下山法结合起来使用,放宽迭代法与下山法结合起来使用,放宽初值初值的选取范围,即将(的选取范围,即将(8.4.1)式修改为:)式修改为: 其中,其中, 称为下山因子,选择下山因子时,希望称为下山因子,选择下山因子时,希望 满足下满足下山法具有的单调性,即山法具有的单调性,即这种算法称为这种算法
27、称为Newton下山法。下山法。在实际应用中,可选择在实际应用中,可选择 。六、六、牛顿法的变形牛顿法的变形1 1、牛顿下山法、牛顿下山法牛顿下山法的计算步骤:牛顿下山法的计算步骤:(1 1)选取初始近似值)选取初始近似值x x0 0;(2)取下山因子取下山因子 = 1= 1;(3)计算计算(4)计算)计算f (xk+1),并比较并比较 与与 的大小,分以下二种情况的大小,分以下二种情况1)若)若 ,则当,则当 时,取时,取x* xk+1,计算过程结束;计算过程结束;当当 时,则把时,则把 xk+1 作为新的作为新的 近似值,并返回到(近似值,并返回到(3)。)。 2)若)若 ,则当,则当 且
28、且|f(xk+1)| ,取,取x* xk,计算过程结束;计算过程结束;否则若否则若 ,而,而 时,则把时,则把xk+1加上一个适当选定的小正数,加上一个适当选定的小正数,即取即取xk+1+ 作为新的作为新的xk值,并转向(值,并转向(3)重复计算;当)重复计算;当 ;且;且 时时,则将下山因子缩小一半,取,则将下山因子缩小一半,取 /2代入,并转向(代入,并转向(3)重复计算。)重复计算。 例例5 5:求方程:求方程f f ( (x x) = ) = x x3 3 x x 1 = 0 1 = 0 的根。的根。k xk010.611/251.14063211.36681311.32628411.
29、32472牛顿下山法的计算结果:牛顿下山法的计算结果:牛顿迭代法每迭代一次都需计算函数值牛顿迭代法每迭代一次都需计算函数值 和导数值和导数值 计算量比较大;且迭代过程中计算计算量比较大;且迭代过程中计算 时,仅利用了时,仅利用了 点的信息,点的信息,而没有充分利用已经求出的而没有充分利用已经求出的 ;在导数计算比较麻烦;在导数计算比较麻烦或难以求出时,或难以求出时, 迭代格式构造迭代格式构造 (2) 构造方法:将构造方法:将Newton迭代格式中的导数用差商代替。迭代格式中的导数用差商代替。 2、割线法、割线法:(1) (1) 构造思想:用割线的斜率代替牛顿迭代法中切线的斜率;构造思想:用割线
30、的斜率代替牛顿迭代法中切线的斜率;设法避开导数值的计算,因此可以采用设法避开导数值的计算,因此可以采用离散牛顿法(割线法)离散牛顿法(割线法)。 一个自然的想法就是在充分利用一个自然的想法就是在充分利用“旧信息旧信息”的同时,的同时,割线法的几何意义割线法的几何意义x0x1切线切线割线割线切线斜率切线斜率 割线斜率割线斜率x2割线法迭代格式割线法迭代格式:割线法的局部收敛性与收敛速度割线法的局部收敛性与收敛速度 设设 ,在在 , ,且且 的某一邻域的某一邻域 内二阶连续可微,当内二阶连续可微,当 时,由割线法产生的序列时,由割线法产生的序列 收敛于收敛于 ,且收敛阶至少为,且收敛阶至少为1.6
31、18。 3 3、 双点弦截法双点弦截法 :切线斜率切线斜率 割线斜率割线斜率初值初值 x0 和和 x1。x0x1x24 非线性方程组的数值解法非线性方程组的数值解法(1) 构造思想:用线性方程组近似非线性方程组,由线性构造思想:用线性方程组近似非线性方程组,由线性方程组解得的向量序列,逐步逼近非线性方程组的解向量。方程组解得的向量序列,逐步逼近非线性方程组的解向量。 (2) 构造方法:若记构造方法:若记一、一、 非线性方程组的牛顿迭代法非线性方程组的牛顿迭代法则非线性方程组则非线性方程组其中其中 ,且,且 中至少有一个是中至少有一个是 的非线的非线性函数。类似于性函数。类似于的情况,可将单变量
32、方程求根的方法推广的情况,可将单变量方程求根的方法推广到非线性方程组。若已给出方程组的一个近似根到非线性方程组。若已给出方程组的一个近似根 。将函数。将函数 的分量的分量 在在 用多元函数泰勒展开,用多元函数泰勒展开,并取其线性部分可表示为并取其线性部分可表示为 (8.5.1)令上式右端为零,得到线性方程组令上式右端为零,得到线性方程组(8.5.2) 其中其中称为称为 的的Jacobi矩阵,求解线性方程组(矩阵,求解线性方程组(8.5.2),并记解为),并记解为 ,则得,则得 这就是解非线性方程组(这就是解非线性方程组(8.5.1)的)的Newton迭代法。迭代法。Newton迭迭代法具有二阶
33、的收敛速度,但对初值的要求很高,即充分靠近代法具有二阶的收敛速度,但对初值的要求很高,即充分靠近解解 。图8.5.1二、全球定位系统的求解二、全球定位系统的求解:解:卫星解:卫星 的位置如图的位置如图9.5.1所示,假设所示,假设 表示表示R的当前位置,则它满足方程组的当前位置,则它满足方程组 其中,光速其中,光速 ,上述方程组无疑是非线性的,但,上述方程组无疑是非线性的,但很容易将所有二次项都消去,从而得到很容易将所有二次项都消去,从而得到 由求解非线性方程组的由求解非线性方程组的Newton迭代法知迭代格式为迭代法知迭代格式为其中其中, 使用使用Matlab求解得迭代求解得迭代4次就可以得
34、到相当精确的结果。次就可以得到相当精确的结果。 nxyzt00.00.00.010010.00716.3687270.374601-2.4039719.98521824.9840633.0182411.04630310000.26635.0001602.9998080.9995859999.99845.0000003.0000001.00000010000.000它的解是:它的解是:历史与注记历史与注记 艾萨克艾萨克牛顿(牛顿(Isaac Newton 16421727) 是英国物理学家、数学家、天文学家和自是英国物理学家、数学家、天文学家和自然哲学家。然哲学家。1643年诞生于英格兰林肯郡乌
35、尔索普年诞生于英格兰林肯郡乌尔索普镇。镇。1727年卒于伦敦。年卒于伦敦。1665年他发现了二项式定年他发现了二项式定理,理,1669年担任卢卡斯讲座的教授,年担任卢卡斯讲座的教授,1696年牛顿年牛顿任造币厂监督,任造币厂监督,1699年升任厂长年升任厂长,1705年因改革币年因改革币制制 牛顿是有史以来最伟大的科学家之一,他在力学、数学、牛顿是有史以来最伟大的科学家之一,他在力学、数学、光学、热学、天文学和哲学方面都有突出的贡献。他在数学方光学、热学、天文学和哲学方面都有突出的贡献。他在数学方面的贡献为:牛顿将古希腊以来求解无穷小问题的种种特殊方面的贡献为:牛顿将古希腊以来求解无穷小问题的
36、种种特殊方法统一为两类算法:正流数术(微分)和反流数术(积分),法统一为两类算法:正流数术(微分)和反流数术(积分),与此同时,他还在与此同时,他还在1676年首次公布了他发明的二项式展开定理。年首次公布了他发明的二项式展开定理。并和并和G.W.莱布尼茨几乎同时创立了微积分学。牛顿在数值计算莱布尼茨几乎同时创立了微积分学。牛顿在数值计算上的主要贡献有:牛顿插值法、牛顿积分法、牛顿迭代法等。上的主要贡献有:牛顿插值法、牛顿积分法、牛顿迭代法等。 有功受封为爵士,有功受封为爵士,16721672年起他被接纳为皇家学会会员,年起他被接纳为皇家学会会员,17031703年被选为皇家学会主席直到逝世。年
37、被选为皇家学会主席直到逝世。 关于特殊的非线性方程求根问题的迭代法最早出现在古关于特殊的非线性方程求根问题的迭代法最早出现在古希腊、巴比伦和印第安人的著作中。牛顿法的只有一部分属希腊、巴比伦和印第安人的著作中。牛顿法的只有一部分属于牛顿本人,于牛顿本人,1669年牛顿第一次提出了与现在牛顿法基本等年牛顿第一次提出了与现在牛顿法基本等价的方法,但令人惊讶的是该方法并没有使用导数,而是基价的方法,但令人惊讶的是该方法并没有使用导数,而是基于二项展开式,因此只适用于多项式。于二项展开式,因此只适用于多项式。1690年,拉弗森对牛年,拉弗森对牛顿法作了简化和改进,称为牛顿顿法作了简化和改进,称为牛顿拉
38、弗森法。在牛顿法中使拉弗森法。在牛顿法中使用导数是由辛普森用导数是由辛普森1740年首次提出的,并将其从一维空间推年首次提出的,并将其从一维空间推广到多维空间,这才是现在所称的牛顿法。广到多维空间,这才是现在所称的牛顿法。1798年,拉格朗年,拉格朗日提出了牛顿法简单而精巧的表达式,并在日提出了牛顿法简单而精巧的表达式,并在1831年由傅立叶年由傅立叶作了推广。非线性方程组的大多数方法都采用了牛顿法的设作了推广。非线性方程组的大多数方法都采用了牛顿法的设计原理,并以牛顿法为度量标准。计原理,并以牛顿法为度量标准。“牛顿法牛顿法”已成为将不同已成为将不同类型非线性方程组局部线性化的一个经典范例。
39、关于一维非类型非线性方程组局部线性化的一个经典范例。关于一维非线性方程组的解法的经典方法的详细资料可参考文献线性方程组的解法的经典方法的详细资料可参考文献1,2,3,若想进一步了解近期的研究成果可参考文献,若想进一步了解近期的研究成果可参考文献4,5 参考文献参考文献1J.F.Traub.IterativeMethodsfortheSolutionofEquations.PrenticeHall,EnglewoodCliffs,NJ,1964.2A.M.Ostrowski.SolutionofEquationsandSystemsofEquations.Academic,NewYork,2dedition,1966.3A.S.Householder.TheNumericalTreatmentofaSingleNonlinearEquation.McGraw-Hill,NewYork,1970.4C.T.Kelley.IterativeMethodsforLinearandNonlinearEquations.SIAM,Philadelphia,PA,1995.5W.C.Rheinboldt.MethodsforSolvingSystemsofNonlinearEquations.SIAM,Philadelphia,PA,2dedition,1998.