现代科学与工程计算-孟大志-第三篇、第9章.pptVIP

第三篇 方程求解 这种迭代格式称为割线法。割线法不但避免了计算导数，还充 分利用了“旧信息”，是对Newton迭代法的一个重要改进。 Newton迭代法和牛顿下山法每次迭代只需用到前一步信息 ，它们属于线性单步法。而割线法每计算一次 需要用到前 面两步的信息 格式是建立在线性插值多项式基础上的，所以自然会想到在二 次插值多项式基础上建立迭代格式，用构造的抛物线的零点作 为 ，因此它属于线性多步法。割线法的迭代 或许会得到更好的迭代方法，称这样得到的迭代法为抛 物线法，也称之为Muller法。一般一条抛物线有两个实零点时， 取距 较近的那个零点作为 。值得注意的是 使用Newton迭代法需给出一个初值 ，而使用 割线法必须给出两个初值 。 2．局部收敛性与收敛速度 设 ， ,且 在 的某一邻域 内二阶连续可微，当 时，由割线法产生的序列 收敛于 ，且收敛阶至少为1.618。 9.5 非线性方程组的数值解法 (1) 构造思想：用线性方程组近似非线性方程组，由线性 方程组解得的向量序列，逐步逼近非线性方程组的解向量。 (2) 构造方法：若记 9.5.1 非线性方程组的牛顿迭代法 则非线性方程组 其中 ，且 中至少有一个是 的非线 性函数。类似于 的情况，可将单变量方程求根的方法推广 到非线性方程组。若已给出方程组的一个近似根 。将函数 的分量 在 用多元函数泰勒展开， 并取其线性部分可表示为 （9.5.1） 令上式右端为零，得到线性方程组 （9.5.2） 其中 称为 的Jacobi矩阵，求解线性方程组（9.5.2），并记解为 ，则得 这就是解非线性方程组（9.5.1）的Newton迭代法。Newton迭 代法具有二阶的收敛速度，但对初值的要求很高，即充分靠近 解 。 图9.5.1 9.5.2 全球定位系统的求解: 解：卫星 的位置如图9.5.1所示，假设 表示R的当前位置，则它满足方程组 其中光速 ，上述方程组无疑是非线性的，但 很容易将所有二次项都消去，从而得到 * * * 第9章 非线性代数方程（组）求解 9.1 问题驱动：全球定位系统（GPS） 人类对导航和定位的需求是伴随着人类整个文明历史的进步而发展的，中国古代“四大发明”之一的指南针是最早的定位仪器和系统，其后还有经纬仪以及近代的雷达。如图9.1.1所示全球定位系统（GPS）是基于卫星的导航系统，最早由美国和前苏联分别在80年代研制，并于1993年正式投入使用。现代社会中全球定位系统越来越深入到人们生活的方方面面。例如市场上出售的手持型GPS，定位的精度可以达到10米以内，这无疑给旅行者提供了方便；安装有GPS的儿童手表，家长在家里的计算机上可以追踪到孩子的位置，防止儿童走 失；安装有GPS系统的汽车可以帮助新司机辨 识道路等等。 图9.1.1 卫星定位示意图 美国和前苏联的GPS都包括有24颗卫星，它们不断地向地 球发射信号报告当前位置和发出信号的时间，卫星分布如图 9.1.2所示。它的基本原理是：在地球的任何一个 位置，至少同时收到4颗以上卫星发射的信号。 设地球上一个点R，同时收到卫星 发射的信号， 假设接收的信息如表9.1.1所示。请设法确定R点的位置。 图9.1.2 卫星分布图 表9.1.1 GPS导航问题可归结为求解非线性代数数方程组 ,当 时就是单个方程. （9.1.1） 。 . 其中， 可以是代数方程，也可以是超越方程。使 成立的x 值称为方程的根，或称为 计算中，如电路和电力系统计算、非线性力学、非线性微（ 积分）方程、非线性规划（优化）等众多领域中，问题的求 解和模拟最终往往都要解决求根或优化问题。前一种情形要 求出方程（组）的根；后一种情形则要求找出函数取最大或 最小的点。即使是对实验数据进行拟合或数值求解微分方程， 也总是将问题简化成上述两类问题。上述除少数特殊方程外， 大多数非线性代数方程（组）很难使用解析法求解精确解， 一般需要通过一些数值方法逼近方程的解。这里主要介绍单 个方程的数值解法，方程组也可以采用类似的方法，将放在 后面讨论。 的零点。科学与工程 9.2 二 分 法 使用数值方法求解非线性代数方程（组）需要解决以下 几个问题： （1）根的存在性：方程是否有根？如果有根，有几个根？ （2）根的分布范围：如果有根，这些根大致分布在什么 位置？如何将有根区间隔离出来？ （3）根的精确化：如何选择求根法，使之简单、经济、 有效？ 搜索算法是方程求根的一种古老而有效的方法，二分法是 搜索