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