104753150634吴丽焕用有限差分格式解抛物型方程.doc

用有限差分格式解抛物型方程 姓名：吴丽焕 学号：104753150634 专业：应用数学 用有限差分格式解抛物型方程 摘要 本文以偏微分方程在初值问题的数值解法为核心研究内容，主要研究了向前差分格式，向后差分格式，六点对称格式以及格式的构造，并用不同方法对这几种不同格式的稳定性进行了判断。 关键词：微分方程数值解 有限差分法 稳定性 1前言 微分方程是数学的一个重要分支，同样也是微积分、变分法、控制论、复变函数、组合数学、组合拓扑等学科的基础。如果一个微分方程中出现的未知数只含有一个自变量，这个方程叫做常微分方程；如果一个微分方程中出现多元函数的偏导数，或者说未知函数和几个变量有关，而且方程中出现未知数对几个变量的导数，那么这种方程就是偏微分方程。 求偏微分方程的定解问题可以先求出它的通解，然后再用定解条件确定出函数。但一般来说，在实际中通解是不易求解的，用定解条件确定函数更是比较困难的。应该指出偏微分方程的定解虽然各有解法，但是我们不能忽视由于某些原因有许多定解条件是不能严格解出的，只可以用近似方法求出满足实际需要的近似程度的近似解。常用的方法有变分法和有限差分法。有限差分法是把定解问题转化成代数方程，然后用计算机进行计算。 2有限差分法简介 有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线性组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是有所有单元上的近似解构成。 有限差分方法是计算机数值模拟最早采用的方法，至今仍被广泛使用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。 从差分的空间形式来考虑，可分为中心格式和迎风格式。 考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。 目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。 构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式有三种：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写为由各个变量或其导数的节值点与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解，采用不同的权函数和插值函数形式，便构成不同的有限元方法。 在采用数值计算方法求解偏微分方程时，若将每一处导数由有限差分近似公式代替，从而把求解偏微分方程的问题转换成求解代数方程的问题，即所谓的有限差分法。有限差分法求解偏微分方程的步骤如下： 区域离散化，即把所给偏微分方程的求解区域细分为由有限个点组成的网格； 近似代替，即采用有限差分法代替每一个格点的导数； 逼近求解。换而言之，这一过程可以看作是用一个插值多项式及其微分来替代偏微分方程的解的过程。 3差分逼近的基本概念 3.1 差分方程 以最简单的以为对流问题为例，引入用差分方法求偏微分方程数值解的一些概念。考虑对流方程的初值问题 （3.1） 区域的剖分： 网格剖分可以采用两组平行于x轴和t轴的直线形成的网覆盖区域，它们的交点称为网格节点节点记为。间距称为空间步长，间距称为时间步长。 Taylor公式：设在的某个邻域内具有直到阶的导数，则有 是余项，且 微商与差商之间的关系式 （向前差商） （3.2） （向前差商） （3.3） （向后差商） （3.4） （中心差商） （3.5） 由于u是方程（3.1）的解，所以 （3.6