第二章有限差分法初步-1题材.ppt

一、差商与微商;;（ii）偏差---Taylor级数展开;;（iii）微商与差商的几何意义 ;图2.1表示了差商与微商之间的关系。应当指出，用不同方法得到的差商去代替微商，它们带来的偏差是不同的。;中心差商，取向右差商与向左差商的平均值：; 将Taylor级数写成：; 由式（2.7）可得;（2.9）+（2.10），得到;用向右差商与向左差商代替微商，其截断误差为与 ;上述一阶差商一般仍是x的函数，对它们还可以求差商。这种一阶差商的差商称为二阶差商，它是二阶微商的近似，常用向右差商的向左差商来近似二阶微商，即：;根据式（2.7）+（2.8），可得; 结论： 由于用差商代替微商必然带来截断误差，相应地用差分方程代替微分方程也必然带来截断误差。这是有限差分法固有的。因此，在应用有限差分法进行数值解时，必须对差分的构成及其对方程造成的误差引起注意。 ;§2 从微分形式出发的差??格式;问题是求图2.2所示的边值问题的解，其数学表达如下，方程：;边界条件：; 在问题的提法已经明白之后，差分格式的构成： ; 1．区域离散化 所谓离散化，就是把几何上连续的区域用一系列网格线把它划分开。一般说来，网格形式应视几何区域的不同而不同，对于矩形区域而言，用矩形的网格，如图2.2，用五条水平网线与五条垂直网线把矩形区域离散掉。网线与网线的交点称之为“节点”，节点与节点的距离称之为步长，x方向的步长表示为 ，y方向的步长表示为 。; 节点编号：为便于计算，需对节点逐个编号。常用（i，j） 表示节点位置，其中，i、j是与网线相对应的正整数。;在区域内的节点称“内节点”，在边界上的节点称“边界节点”。图2.2所示边界是规则的,则节点或在区域内，或正好落在边界上。 ;从物理方面对区域离散化可作这样的理解，即认为区域内离散的每个节点，都集中着它周围区域（尺度为步长）的热容，或者说，区域内连续分布的热容都被分别地集中到离散的节点上去了。这样，节点的温度代表着它周围区域的某种平均温度。一系列离散的节点温度值代表着连续区域内的温度分布。 ; 2．差分方程代替微分方程 在上节我们已对有限差分法的数学基础作了简要 的介绍，说明了如何用差商代替微商，以及由此带来的误差。这里介绍用差商代替微商的办法来处理导热方程（2.13），得到相应的差分方程。 ;在（i，j）处存在二阶偏微商 与 这些二阶偏微商所对应的差商可表示成：;;将式（2.18）与（2.19）代入方程（2.13），得;式（2.20）中去掉 项，得到;式（2.21）被称为对应于方程（2.13）的差分方程。方程（2.21）被改写成：;若 ， ，则式（2.22）又被 改写成： ;;3．边界条件的差分形式;; 用T 对 y 的向前差商代替式（2.15）中T 对 y的一阶偏微商，使式（2.15）变成为如下差分形式: ;绝热边界条件： ; 至此，我们对全部节点，包括内节点与边界节点，都用差分形式代替了原来的函数形式。对于内节点上差分形式，我们通称差分方程，因为内节点上温度都是未知的。对于边界节点的差分形式，在边界节点的温度为未知量时，它是差分方程。而对于边界节点为给定的温度时，得到的就不是差分方程了。但在实际应用中，人们往往习惯地把由内节点与边界节点建立起来的差分形式，都统称为差分方程。笼统地讲，一个节点对应一个差分方程。;在边界节点的处理方面还有几点需要强调： ;4．差分格式的构成;也就是说，在组成方程组时，不必把每个节点方程都写出来，而只要写出几个规格化了的方程就可以了。因此，人们把规格化了的，由内节点与边界节点全部差分方程所构成的线性代数方程组，称之为“差分格式”。 ;一般地说，差分格式被写成如下的形式： ;方程组（2.28）可被进一步写成矩阵的形式 :;如果在方程组中去掉其中已知温度节点的那些方程，由此构成的线性代数方程组中，方程的个数等于温度未知的节点个数，也即方程组的未知数。这样的方程组也是差分格式。也可写成式（2.29）。;综上所述，用有限差分法对式（2.13）~（2.17）所组成的边值问题的数值处理，最终归结成求解线性代数方程组（2.29），方程组的解即各节点的温度。如果整个区域的节点足够多，那么，离散节点的温度分布就近似代替了区域内的连续温度分布。;为便于讨论各种差分格式的优缺点，最好把方程组（2.29）中系数矩阵具体地写出来。但当我们着手书写由式（2.21）、（2.24）~（2.27）组成的代数方程组时，发现它所占的版面太大，造成印刷的困难。;所以，为便于书写，采用图2.3所示的网格 ;得到由全部节点组