计算物理基础课件解读.ppt

绪论 计算物理基础 作为一种近似，可取 （a）两个相近的近似数相减时，有效数字会严重损失，实际计算时要尽量避免。[1.5] 例如，当x充分大时，计算下面表达式的值。 定积分 有两个递推公式 对公式（2）进行同样的分析有 4.1.引言 4.2.欧拉法 4.3.龙格-库塔法 4.4.线性多步法 4.5.高阶微分方程 方程中包含未知函数的导数或者微分的方程叫做微分方程,求解微分方程必须附加某种定解方程.微分方程和定解条件一起组成定解问题,定解条件分为初始条件 初值问题 和边界条件 边值问题 两种.未知函数为一元的微分方程叫做常微分方程,未知函数为多元函数,叫做偏微分方程.微分方程中导数的最高阶叫做微分方程的阶.本章主要讨论一阶常微分方程. 常微分方程在物理学中具有重要意义,牛顿第二定理就表示为一组二阶微分方程.前面提到的著名的三体问题就是由九个二阶常微分方程组成的方程组;在量子力学和电动力学中,由薛定谔及麦克斯韦方程分离变量得到本征值问题也是二阶常微分方程的求解问题. 现考虑上述一阶常微分方程初值问题.在常微分方程中，我们已经掌握了一些典型方程的解析解法。但解析方法只能解一些特殊类型的方程,求解从实际当中得出来的微分方程只能用数值方法求近似解.上述初值问题也就是由已知点 x0,y0 求解在某些点x上满足一定精度的y的近似解。 即使对有些函数可以求出解析解.但实际应用时,仍需求它在某一点的值,如果解析解的形式比较复杂,还是需要数值方法来求得.例 y’ 1-2xy y 0 0 其解析解如下,要计算具体的函数值y,还需要积分的方法,不如直接用数值法求解. 定理:设f x，y 在带形区域 a≤x≤b,-∞＜y＜+∞ 上为x，y的连续函数，且对任意的y满足李普希茨 Libusize 条件:存在常数L,使 ｜f x，y1 -f x，y2 ｜≤L｜y1-y2｜ 对所有x∈[a,b]及任意实数y1,y2均成立,则初值问题在[a,b]上有唯一解 在求初值问题的数值解时，我们通常采用离散化方法 数值微分、数值积分、泰勒展式等 ，求在自变量x的离散点a x0＜x1＜x2＜…＜xn b上的准确解y x 的近似值 y0，y1，y2，…，yn.常取离散点x0，x1，x2，…，xn为等距，即 x i+1-xi h，i 0，1，2，…，n-1. h称为步长。下图表示为初值问题在n+1个离散点上的准确解y x 的近似值。 4.2.1 欧拉法 折线法 若将函数y x 在点xi处的导数y’ xi 用 向前 差商来表示，即 例4.1:用欧拉法求初值问题的数值解 取h 0.1 。 隐式欧拉公式 若将函数y x 在点xi处的导数y’ xi 用 向后 差商来表示，即 欧拉公式变为: 两步欧拉公式 为了改善精度,将函数y x 在点xi处的导数y’ xi 用中心差商来表示，即 欧拉公式变为:误差正比与h3,是二阶精度 4.2.2 改进的欧拉法 欧拉法虽然形式简单，计算方便，但比较粗糙，精度也低。特别当y y x 的曲线曲率较大时，欧拉法的效果更差。为了达到较高精度的计算公式，对欧拉法进行改进，将在一点 xi，yi 的切线斜率f xi，yi 用两点的平均斜率来代替，即 代入欧拉公式得 不难发现，欧拉公式是关于yi+1的显式，只要已知yi，经一次计算可立即得到yi+1的值;而改进的欧拉公式中的yi+1以隐式给出，且yi+1含在函数f xi+1 ,yi+1 中，因此，通常用迭代法求解。具体做法是：先用欧拉公式求出一个y 0 i+1作为初始近似，然后再用改进的欧拉公式进行迭代，即 直到满足 则 4.2.3 预估校正法 改进的欧拉公式在实际计算时要进行多次迭代，因而计算量较大。所谓预估校正法，就是先用欧拉公式算出yi+1的预估值y p i+1，然后再用改进的欧拉公式式进行一次迭代便得到校正值y c i+1，即 虽然上式仅迭代一次，但因进行了预先估计，故精度却有较大的提高。 例4.2:求解初值问题 我们还是假设yi y xi ，利用泰勒级数展开求y xi+1 。若取右端前有限项作为y xi+1 的近似值，就可得到计算y xi+1 的各种不同截断误差的数值公式。 例如取前两项可得到 y xi+1 ≈y xi +hy’ xi y xi +hf xi，y xi yi+hf xi，yi ----欧拉法 若取前三项，可得到截断误差为O h3 的公式 类似地，若取前k项作为y xi+1 的近似值，便得到截断误差为O hk 的数值计算公式。这些公式的计算必须依赖于求y xi 的k阶导数，除非f x，y 足够简单，否则直接用泰勒展开法求解较为复杂。但是泰勒级数展开法的基本思想是许多数值方法的基础。 已经知道，初值问题等价于 二阶龙格库塔公式 将预估校正公式改写成更一般的形式