常微分方程与偏微分方程数值方法比较.docxVIP

常微分方程与微分方程数值方法比较1. 微分方程数值方法的有关概念首先回顾微分方程的定义与分类。含有自变量、未知函数及其导数（微分或偏导数）的方程称为微分方程；如果未知函数只含有一个变量，则称为常微分方程；如果未知函数含有若干个变量，则称为偏微分方程。微分方程中未知函数的导数或偏导数的最高阶次称为微分方程的阶。例如：微分方程是一阶常微分方程，而是二阶偏微分方程。所有使微分方程成为等式的函数，都是微分方程的解；在 n 阶微分方程中，将微分方程的含有 n 个任意常数的解称为该微分方程的通解。为确定微分方程通解中的任意常数而需要的条件称为定解条件；定解条件可以分为初始条件和边界条件两类。由微分方程和定解条件一起构成的问题称为微分方程定解问题。根据定解条件的不同，常微分方程分为初值问题和边值问题；若定解条件是描述函数在一点（或初始点）处状态的，则称为初值问题，一阶常微分方程初值问题的一般形式为：若定解条件描述了函数在至少两点（或边界）处状态的称为边值问题，例如：2.常微分方程数值方法有限差分法是常微分方程中数值解法中通常有效的方法，建立差分算法的两个基本的步骤：1. 建立差分格式，包括：a. 对解的存在域剖分；b. 采用不同的算法可得到不同的逼近误差—截断误差（相容性）；c.数值解对真解的精度—整体截断误差（收敛性）；d.数值解收敛于真解的速度；e. 差分算法—舍人误差（稳定性）.2.差分格式求解，将积分方程通过差分方程转化为代数方程求解，一般常用递推算法。差分方法的基本思想“就是以差商代替微商”，差分形式如下：1.2.3.4.对初值问题中的导数进行不同的离散化处理。对于初值问题的数值解法，首先要解决的问题就是如何对微分方程进行离散化，建立求数值解的递推公式。递推公式通常有两类，一类是计算yi+1时只用到xi+1, xi 和yi，即前一步的值，因此有了初值以后就可以逐步往下计算，此类方法称为单步法；其代表是龙格—库塔法。另一类是计算yi+1时，除用到xi+1,xi和yi以外，还要用到，即前面k步的值，此类方法称为多步法；其代表是亚当斯法。欧拉（Euler）方法是解初值问题的最简单的数值方法。初值问题的解y=y(x)代表通过点的一条称之为微分方程的积分曲线。积分曲线上每一点的切线的斜率等于函数在这点的值。通常取 (常数),则Euler法的计算格式 i=0,1,…,n 带入不同的差分方法，构造出了隐式欧拉法、梯形公式、中点欧拉公式方法优点缺点显式欧拉简单精度低隐式欧拉稳定性最好精度低, 计算量大梯形公式精度提高计算量大中点公式精度提高, 显式多一个初值,可能影响精度为了提高数值解法的精度，又发展出改进欧拉法，然后，通过在这一步内多预报几个点的斜率值，然后将其加权平均作为平均斜率，则可构造出更高精度的计算格式，这就是龙格—库塔（Runge-Kutta）法的基本思想，亚当姆斯（Adams）方法是设用xi,xi+1两点的斜率值加权平均作为区间上的平均斜率，有计算格式将在xi处Taylor展开代入计算因此有这样导出的计算格式称之为二阶亚当姆斯格式。类似地可以导出三阶亚当姆斯格式。常微分方程初值问题的基本数值解法。包括单步法和多步法。单步法主要有欧拉法、改进欧拉法和龙格—库塔方法。多步法是亚当姆斯法。它们都是基于把一个连续的定解问题离散化为一个差分方程来求解，是一种步进式的方法。用多步法求常微分方程的数值解可获得较高的精度。3.偏微分方程数值方法差分方法又称为有限差分方法或网格法，是求偏微分方程定解问题的数值解中应用最广泛的方法之一。它的基本思想是：先对求解区域作网格剖分，将自变量的连续变化区域用有限离散点（网格点）集代替；将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替；通过用网格点上函数的差商代替导数，将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组（称为差分格式）。如果差分格式有解，且当网格无限变小时其解收敛于原微分方程定解问题的解，则差分格式的解就作为原问题的近似解（数值解）。因此，用差分方法求偏微分方程定解问题一般需要解决以下问题：（i）选取网格；（ii）对微分方程及定解条件选择差分近似，列出差分格式；（iii）求解差分格式；（iv）讨论差分格式解对于微分方程解的收敛性及误差估计。下面我们只对偏微分方程的差分解法作一简要的介绍。（二维）一、二阶偏导数的有限差分近似1.2.3.5.如上所示，其中有向前差分，向后差分，中心差分，在实际计算中，会通过对三种方式的组合差分来对不同维方向进行差分，如迎风格式、蛙跳格式、Lax格式。要得到偏微分方程的唯一解，需要定解条件，即问题的初始和边界条件。边界条件有三类：第一类是在边界上直接给出未知函数的数值,也称为 Dirichlet 条件；第二类是在边界上给