数值分析 第三篇方程求解.pptVIP

付出多少得到多少！ 计算两次函数值还能提高误差阶数吗？多展开一项： 四阶Rung-Kutta算法 R-K方法的绝对稳定区域 §5. 常微分方程组与高阶方程的数值解 2) 方程组的 R-k 法 数值实验3： 蝴蝶效应与混沌解 1963年 Lorentz 给出简化大气环流模型： 偏微分方程 物理定律：物理量与其偏导数的关系。 例1，物体内部不具有热源时，它的温度分布与变化 用物理量 u=u（t，x，y，z）表示，则规律为： 溶解物在溶液中的密度 u = u（t，x，y，z）也满足(1). 又称扩散方程. 波动方程与位势方程 这里 a 是声音在空气中的传播速度，也可以表示电磁波与一般波动，称（2）为波动方程。 除表示热平衡外，还可表示静电场、静磁场、经典引力场或稳态流动等。称 Laplace 方程或调和方程。 例2, 声波在空气中传播时，如果 u 表示压强，则规律： 偏微分方程数值方法简介（双曲型方程的差分解法） 差分解法是偏微分方程的最基本解法。 即将偏导数用差商来代替。 差分格式建立 差分方程 差分格式 若对 t 用向前差商，对 x 用向后差商，得差分格式： 隐藏格式 对 t 用向后差商，对 x 采用中心差商，得差分格式： 差分格式收敛性 库郎条件：差分格式的依赖区域应包含微分方程的依赖区域. 差分格式的稳定性 一点引入的误差,在逐层计算时累计误差是否淹没真值——稳定性。 Newton法解方程组例 用 Newton 法求（1，1）附近的根。 当相邻两组的最大差小于给定精度而求出解。 §６ 最速下降法 最速下降法 几何解释 负梯度方法 等高线上的函数值相同而闭曲线内部值小于外部值。 负梯度分析的步长选取（1） 负梯度分析的步长选取（2） 负梯度分析的步长选取（3） 最速下降下是线性收敛的，对初始点依赖小；而Newton法收敛快，对初始点依赖大。常结合使用。 第十章 常微分方程初值问题的数值解法 最为普遍的模型，最为复杂的问题 建模框架： 明确地叙述问题→分析其中的物理、化学、 生物定律→寻找数量关系： 变化率=输入率 - 输出率 世界本质上是非线性的：混沌与现实 一阶常系数微分方程的初值问题——基本模型 微分方程的解形成曲线族，或解流形，初值条件定出 其中一条曲线。 数值解的基本方法 §1 简单方法——解法的基础 2，Euler法的分析与解释 C. 数值积分 3, Euler方法的误差估计 2), 总体方法误差(1) 2), 总体方法误差(2) 微分方程数值解的稳定性 Euler法的绝对稳定区域 （二）向后Euler 方法 （二）向后Euler 法的稳定性 （三）梯形公式 梯形公式的稳定性 §2. Runge-Kutta法 A. 一步法分析：Euler 法与梯形法 用待定系数法来求这些参数，使上述3，中相同项更多 B. 一般公式 理论的简单例——梯形法 第三篇．方程求数值解 方程是什么？ 含未知变量的等式； 函数是什么？ 变量之间的关系：世界上最普遍的结构； 方程是特定的函数表示（方程建模）： 定态关系－代数方程与超越方程， 动态关系－微分方程（OD,ED）; 方程求解：逻辑逆问题和隐含关系 第九章 非线性方程及方程组的解法 §1. 非线性方程实根的对分法(二分法) 二分法的收敛性 §2. 迭代法 迭代过程的几何表示 收敛充分性定理 收敛定理证明 §3. 迭代加速收敛方法 Aitken方法 Aitken方法的几何解释 §4. Newton 法 非线性问题的最简单解法是线性近似. Newton 法的几何解释 §５. 解非线性方程组的Newton迭代法 Newton迭代线性方程组。 行列式在（xn, yn)处取值。