5-非线性方程组的数值解法及最优化方法.pptVIP

非线性方程组的数值解法 若对任意 都有一个实数 与之对应，且满足： （1）非负性：当 时， ；当 时， ； （2）齐次性：对任何 ， ； （3）三角不等式：对任意 ，都有 （4）相容性：对任意 ，都有 ， 则称 为 上矩阵 的范数，简称矩阵范数。 非线性方程组的数值解法 考虑如下方程组 式中 均为 的多元函数，向量形式为 其中 非线性方程组的数值解法 当 ，并且 中至少有一个是自变量 的非线性实函数时，称方程组 为 非线性方程组。其求根问题就是确定方程组在指定范围内的 一组解，可以通过对单个非线性方程求根问题的直接推广得 到非线性方程组的求解算法。 非线性方程组的数值解法 常用解法分为两类：一类是线性化方法，将非线性方程组用 一个线性方程组来近似，由此构造一种迭代公式，逐次逼近 所求的解；另一类是属于求函数极小值的方法，即由非线性 函数 构造一个模函数，例如构造函数 然后通过各种下降法或优化算法求出模函数的极小值点，此 极小值点即为非线性方程组的一组解。 非线性方程组的数值解法 不动点迭代法：根据非线性方程求根的迭代法，将方程组改 写为如下等价方程组 构造迭代公式 选取初始向量 则由迭代公式可以得到一个向量序列 。如果 方程组有唯一解向量 ，并且 ，则 可 作为逐次逼近 的近似解。 非线性方程组的数值解法 如果把迭代公式写为向量形式 并记矩阵 为 则可以证明当 时，迭代公式是收敛的。 非线性方程组的数值解法 例题1：用迭代法解如下非线性方程组 取初值 。 解：构造迭代公式 非线性方程组的数值解法 式中 所以有 取初值 ，在 附近 ，所以迭 代公式收敛。 非线性方程组的数值解法 非线性方程组的数值解法 非线性方程组的数值解法 上述迭代公式与求解线性方程组的雅克比迭代公式形式相 同，可以对其进行改进，构造求解非线性方程组的高斯-赛 德尔迭代公式，即 对上例采用高斯-赛德尔迭代公式计算 迭代计算过程如下表所示 非线性方程组的数值解法 非线性方程组的数值解法 牛顿迭代法：根据求解非线性方程的牛顿迭代法，如果已经 给出方程组 的一个近似根 ，则 可把函数 的分量 在 处按多元函数泰 勒公式展开，取其线性部分做近似，得 则得到线性方程组 非线性方程组的数值解法 方程组的解为 上式即为求解非线性方程组的牛顿迭代公式。式中 称为 的雅克比矩阵。 非线性方程组的数值解法 例题2：用牛顿迭代法求解下面非线性方程组 计算时取初始值 。 解：先求雅克比矩阵 非线性方程组的数值解法 由牛顿迭代公式得到 迭代计算过程如下表所示。 非线性方程组的数值解法 非线性方程组的数值解法 非线性方程组的数值解法 练习题：用牛顿迭代法求解方程组 取 结果： 非线性方程组的数值解法 应用经过海底一次反射到达水听器阵的特征声线传播时间， 来反演海底参数。假设水中和沉积层声速都是恒定的，海底 沉积层上界面水平，下界面倾斜。特征声线由水中声源出发 折射进入沉积层，经过沉积层的下界面反射后，再折射进入 水中，由水中水听器阵接收。特征声线的传播时间为声线在 水中和沉积层中的传播时间之和。 三维坐标关系如图所示： 非线性方程组的数值解法 非线性方程组的数值解法 由声线传播的几何关系和折射定律，一个接收点可以导出如 下6个方程： 非线性方程组的数值解法 方程中 为声线的位置参数， 为海底倾角 为声源到沉积层的距离， 为声源处的沉积层厚