计算方法第6章常微分方程数值解法讲解.ppt

第6章 常微分方程数值解法 §1 引言 §2 欧拉法和改进的欧拉法 §3 龙格库塔法 §4 阿达姆斯方法 §5 二阶线性常微分方程边值问题的数值解 微分的方程：含有未知函数的导数或微分的方程 §1 引言 §2 欧拉法和改进的欧拉法 2.1 欧拉法(折线法) 若将函数y(x)在点xi处的导数y (xi)用差商来表示，即 用yi近似地代替y (xi)，则： 例 用欧拉法求初值问题 欧拉法虽然形式简单，计算方便，但比较粗糙，精度也低。特别当y=y(x)的曲线曲率较大时，欧拉法的效果更差。 改进的欧拉法： 将在一点(xi，yi)的切线斜率f (xi，yi)用两个点的平均斜率来代替，即 改进的欧拉公式在实际计算时需多次迭代，计算量较大。 预估校正法：就是先算出yi+1的预估值y(p)i+1，然 后再用式进行一次迭代得到校正值y(c)i+1，即： 每迭代一次，都因进行了预先估计，使得精度有较大的提高。 在实际计算时，对预估校正法时常作下列变换： §3 龙格库塔法 第2章 一元线性方程的解法 §1 二分法 §2 迭代法 §3 切线法(牛顿法) §4 弦截法 §5 加速迭代法 第3章 线性代数计算方法 第4章 插值法 §1 插值问题 §2 线性插值与二次插值 §3 代数多项式插值的存在唯一性 §4 代数多项式的余项 §5 拉格朗日插值多项式 §6 牛顿均差插值多项式 §10 曲线拟合法 第5章 数值积分 §2 牛顿 ― 柯特斯(Newton―Cotes) 公式 §3 复合求积公式 第6章 常微分方程数值解法 §2 欧拉法和改进的欧拉法 §3 龙格库塔法 第6章 常微分方程数值解法 《计算方法》 微分方程的解是xoy平面上的一簇积分曲线，这簇积分曲线上任意点(x，y)的斜率为 即 f (x，y)。 确定一条特定的积分曲线就是解微分方程的初值问题。 x y 能够确定通解中某一定值的条件叫定解条件，常见的定解条件为初值条件，满足初值条件的特解，成为积分的初值问题。表示为： ,...) 2 , 1 , 0 ( ) ( = ? n y x y n n 微分方程的解析解是指解出y (x)，即：要解出所有x点上的准确函数值y (x)。 方程的数值解是指通过某种方法去获得解析解y (x)在点xi上的近似值yi，即用yi近似代替y (xi) ，表示为： ′ 欧拉公式有很明显的几何意义。 就是欧拉公式 的数值解(取h=0.1)。 由欧拉计算公式有： 解 因为 2.2 改进的欧拉法 改进的欧拉公式： 隐式 用迭代法求解改进的欧拉公式。 具体做法是：先用欧拉公式求出一个 作为初始近似值，再用改进的欧拉公式进行迭代求解，即： 预估： 校正： 取 2.3 预估校正法 k1 k2 进一步对预估校正法进行扩展： 利用泰勒级数展开求 第1章 误差 §1 计算方法的概念、意义及特点 §2 误差 §3 算术运算结果的误差 §4 算法的数值稳定性 误差定义、种类、绝对误差、相对误差、误差限、有效数字、两者关系 避免方法和注意事项 停止条件、二分次数等 求根步骤 收敛条件：李普希茨条件、迭代函数导数小于1 §1 高斯消去法 §4 矩阵的三角分解 §7 迭代法的收敛性 §8 矩阵的特征值与特征向量的计算 §2 高斯―约当消去法 §6 迭代法 ? 消去过程、上三角、下三角、列主元、全主元 LU分解、杜利特尔分解、克劳特分解 雅可比（Jacobi）迭代、赛德尔（Seidel）迭代 范数定义、迭代收敛条件 乘幂法 * * 第6章 常微分方程数值解法 《计算方法》