工程计算研究论文.docVIP

《工程计算》研究论文 学校： 专业：信息与计算科学 班级：11信息 姓名： 学号： 摘要 科学研究和工程技术中的许多问题均需建立偏微分方程的数学模型。但是除少数特殊情况外，绝大多数偏微分方程问题均难以求出精确解。因此，近似解法就显得尤为重要。而差分法就是一种常用的求偏微分方程近似解的方法。另外，以法国著名数学家和天文学家-拉普拉斯(Laplace)的名字命名的拉普拉斯方程，在科学技术的各个领域有着广泛的应用。本文将给出一些相关的基本概念并通过使用五点差分格式对一个具体的拉普拉斯方程第一边值问题进行求解，以便读者更直观的了解差分法求解偏微分方程的过程。 关键词：五点差分格式；拉普拉斯方程；第一边值问题；Gauss-Seidel迭代法Poisson)方程 但特别地，当时，即为拉普拉斯(Laplace)方程，又称为调和方程：为R2的一个有界区域。 定解条件通常有三类： 第一类边界条件（Dirichlet边界条件）为： 其中为以为边界的有界区域，为分段光滑曲线，称为定解区域，，分别为，上的已知连续函数。 第二类和第三类边界条件可统一表示为： 其中为边界的外法线方向。当时为第二类边界条件，时为第三类边界条件。 二、差分方法 差分方法又称为有限差分方法或网格法，是求偏微分方程定解问题的数值解中应用最广泛的方法之一。 它的基本思想是：先对求解区域作网格剖分，将自变量的连续变化区域用有限离散点(网格点)集代替；将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替；通过用网格点上函数的差商代替导数，将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组(称为差分格式)。如果差分格式有解，且当网格无限变小时其解收敛于原微分方程定解问题的解，则差分格式的解就作为原问题的近似解(数值解)。 因此，用差分方法求偏微分方程定解问题一般需要解决以下问题： 1．选取网格； 2．对微分方程及定解条件选择差分近似，列出差分格式； 3. 求解差分格式； 4. 讨论差分格式解对于微分方程解的收敛性及误差估计。 这里我们先在下面列出几个常用的公式：假设g(x) 在[x0+h , x0-h] 上具有所需的连续导数，则由Taylor公式容易得到（1.1.1）和（1.1.2）式。 （1.1 .1） （1.1.2） 三、以Dirichlet边值问题为例说明差分法求解的一般过程 考虑二维Poisson方程Dirichlet边值问题： 其中 为简单起见，先考虑矩形区域 差分格式的建立 将区间作m等分，记h1 = ( b – a ) / m , xi = a + ih1 , 0≤I ≤m; 将区间[c , d] 作n 等分 ,记 h2 = (d – c ) / n , yj = c + j h2, 0 ≤ j ≤ n. 称h1 为 x 方向的步长h2为 y方向的步长。用两簇平行线 x = xi 0≤ I ≤ m y = yj 0 ≤ j ≤ n 将区域剖分为m*n 个小矩形，称两簇直线的交点 (xi , yj) 为网格结点，如下图所示： 图1 差分示意图 称属于的结点 为内结点，称位于上的结点 为边界结点，显然 为方便起见，记 设 为h上的网格函数。记 称为v 的无穷范数。 在结点处考虑二维Poisson方程Dirichlet边值问题，有 （1.1.3） （1.1.4） 定义h上的网格函数 其中 由（1.1.1），有 将以上两式带入（1.1.3）和（1.1.4），可得 在上式中略去小量项 并用uij代替Uij 得到如下差分格式 （1.1.5） （1.1.6） 称Rij为差分格式（1.1.5）的局部截断误差，它反映了差分格式（1.1.5）对精确解的满足程度，即Rij为差分格式（1.1.5）用精确解代替近似解后等式两边之差 记 则有 2.差分格式的求解 差分格式（1.1.5）和（1.1.6）是以 为未知量的线性方程组。（1.1.5）可改写为 （1.1.7） 记 利用（1.1.6）可将（1.1.7）写为 （1.1.8） 其中 （1.1.8）可进一步写为 （1.1.9） 上述线性方程组的系数矩阵