第九章 偏微分方程差分方法汇总

偏微分方程数值方法偏微分方程（Partial Differential Equation, PDE）是数学中的一种重要的方程类型，它描述了一个函数的多个变量的变化关系。

解决偏微分方程的数值方法在科学和工程领域有着广泛的应用。

本文将介绍几种常见的偏微分方程数值方法，并对其进行详细阐述。

1. 差分法（Finite Difference Method）：差分法是最早也是最直接的一种数值方法，它基于连续函数在一些点的导数可以用它的前向、后向或中心的差商来近似的思想。

偏微分方程的差分格式包括向前差分法、向后差分法和中心差分法等。

对于二维的偏微分方程，可以采用网格化的方式将空间离散化，然后利用差分法进行近似求解。

2. 有限元法（Finite Element Method）：有限元法是一种基于原始形式或变分形式对偏微分方程进行离散化的方法。

在有限元法中，将求解域分割成许多小的、简单的几何单元，然后在每个单元上构建近似解函数和试验函数。

通过构建弱形式并应用基本的变分原理，可以得到离散化的方程组，并通过求解这个方程组来得到数值解。

3. 有限差分法（Finite Difference Method）：有限差分法是一种将连续的偏微分方程离散化成差分方程的方法。

它与差分法的主要区别在于有限差分法不需要对求解域进行网格化，而是直接在连续的求解域上进行离散化。

将偏微分方程中的导数通过差商来近似，然后通过求解离散化的差分方程来得到数值解。

4. 有限体积法（Finite Volume Method）：有限体积法是一种将偏微分方程离散化为离散体积元的方法。

在有限体积法中，将求解域划分成离散的控制体积，然后通过对控制体积的积分运算，将偏微分方程转化为离散的代数方程组。

然后通过求解得到的代数方程组，可以得到数值解。

以上介绍的只是几种常见的偏微分方程数值方法，实际上还有很多其他的方法，如边界元法（Boundary Element Method）、谱方法（Spectral Method）、逆问题方法（Inverse Problem Method）等。

偏微分方程求解方法总结偏微分方程（Partial Differential Equations，简称PDE）是描述自然界中许多现象的重要数学工具。

求解偏微分方程有许多不同的方法，下面将对其中一些常用的方法进行总结和介绍。

I. 分离变量法（Method of Separation of Variables）分离变量法是求解偏微分方程最常用的方法之一。

它的基本思想是将多个变量的偏微分方程分解成一系列只包含一个变量的常微分方程，再通过求解这些常微分方程来获得原偏微分方程的解。

具体步骤如下：1. 根据问题所给的边界条件和初始条件，确定偏微分方程的类型（椭圆型、双曲型或抛物型）以及边界条件的类型（Dirichlet条件、Neumann条件等）。

2. 假设原方程的解可以表示为一系列只包含一个变量的函数的乘积形式，即 u(x, y) = X(x)Y(y)。

3. 将 u(x, y) 和其各个分量的偏导数代入原偏微分方程，得到关于X(x) 和 Y(y) 的常微分方程。

4. 求解得到 X(x) 和 Y(y) 的表达式，并根据给定的边界条件，确定它们的取值。

5. 最后将 X(x) 和 Y(y) 的表达式代入 u(x, y) 的乘积形式，得到原偏微分方程的解。

分离变量法适用于边界条件分离的情况，并且对于较简单的偏微分方程求解效果较好。

II. 特征线法（Method of Characteristics）特征线法主要用于求解一阶偏微分方程，尤其是双曲型和抛物型偏微分方程。

该方法通过引入新的独立变量和新的变量关系，将原偏微分方程转化为一系列常微分方程来求解。

具体步骤如下：1. 根据偏微分方程的类型，确定要求解的未知函数及其偏导数之间的关系。

2. 引入新的自变量和新的关系式，将偏微分方程化为带有新变量的常微分方程组。

3. 将常微分方程组进行求解，并得到新变量的表达式。

4. 根据新的变量表示原方程的解，进而确定未知函数的表达式。

偏微分方程的解法偏微分方程（Partial Differential Equations，简称PDEs）是数学中的一个重要分支，它描述了多变量函数的偏导数之间的关系。

这些方程在自然科学、工程应用和社会科学等领域都发挥着重要作用。

解决偏微分方程是一个复杂而有挑战性的过程，需要运用多种数学方法和工具来求解。

在本文中，我将为您介绍几种常见的偏微分方程的解法，并提供一些示例以帮助您更好地理解。

以下是本文的主要内容：1. 一阶线性偏微分方程的解法1.1 分离变量法1.2 特征线方法2. 二阶线性偏微分方程的解法2.1 分离变量法2.2 特征值法2.3 Green函数法3. 非线性偏微分方程的解法3.1 平移法3.2 线性叠加法3.3 变换法4. 数值方法解偏微分方程4.1 有限差分法4.2 有限元法4.3 谱方法5. 偏微分方程的应用领域5.1 热传导方程5.2 波动方程5.3 扩散方程在解一阶线性偏微分方程时，我们可以使用分离变量法或特征线方法。

分离变量法的基本思路是将方程中的变量分离，然后通过积分的方式求解每个分离后的常微分方程，最后再将结果合并。

特征线方法则是将方程中的变量替换为新的变量，使得方程中的导数项消失，从而简化求解过程。

对于二阶线性偏微分方程，分离变量法、特征值法和Green函数法是常用的解法。

分离变量法的核心思想与一阶线性偏微分方程相似，将方程中的变量分离并得到常微分方程，然后进行求解。

特征值法则利用特征值和特征函数的性质来求解方程，适用于带有齐次边界条件的问题。

Green函数法则通过引入Green函数来求解方程，其特点是适用于非齐次边界条件的情况。

非线性偏微分方程的解法则更加复杂，常用的方法有平移法、线性叠加法和变换法。

这些方法需要根据具体问题的特点选择合适的变换和求解技巧，具有一定的灵活性和创造性。

除了上述解析解法，数值方法也是解偏微分方程的重要手段。

常用的数值方法包括有限差分法、有限元法和谱方法等。

偏微分方程的有限差分法
有限差分法：是一种数学计算概念，是指在计算过程中，以差分的形势来代替微分，从而使整个计算过程具有有限差分法的出发点，以此达到微分议程和积分微分方式数值解的一种计算过程。

微分方程和积分微分方程数值解的方法。

基本思想是把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分方程组，解此方程组就可以得到原问题在离散点上的近似解。

然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解。

在采用数值计算方法求解偏微分方程时，若将每一处导数由有限差分近似公式替代，从而把求解偏微分方程的问题转换成求解代数方程的问题，即所谓的有限差分法。

有限差分法求解偏微分方程的步骤如下：
1、区域离散化，即把所给偏微分方程的求解区域细分成由有限个格点组成的网格；
2、近似替代，即采用有限差分公式替代每一个格点的导数；
3、逼近求解。

换而言之，这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程（Leon，Lapidus，GeorgeF。

Pinder，1985）。

高等数学中的偏微分方程方法偏微分方程是数学中的一类非常重要的方程。

它们广泛应用于物理、工程和其他领域中，如热传导、电路等等。

因此，研究偏微分方程的方法和技巧尤为重要。

在高等数学中，有许多关于偏微分方程的方法，下面我们来介绍其中的几种。

1. 分离变量法分离变量法是解偏微分方程的一种常用方法。

这种方法的基本思想是假设解可以表示为形式为x、y、z等变量的函数之积的形式，然后通过代入相关偏微分方程中去求解出每个变量的解，最终将这些解组合起来得到总体解。

以拉普拉斯方程为例，其定义如下：$\Delta u=\frac{\partial^2 u}{\partial x^2}+\frac{\partial^2u}{\partial y^2}+\frac{\partial^2 u}{\partial z^2}=0$假设解为$u(x,y,z)=X(x)Y(y)Z(z)$，则可以得到：$\frac{1}{X}\frac{\partial^2 X}{\partialx^2}+\frac{1}{Y}\frac{\partial^2 Y}{\partialy^2}+\frac{1}{Z}\frac{\partial^2 Z}{\partial z^2}=0$由于等式左边是一个只关于x的函数与一个只关于y的函数之和，所以这个等式必须等于常数k。

因此，我们可以得到：$\frac{1}{X}\frac{\partial^2 X}{\partial x^2}=k_1$，$\frac{1}{Y}\frac{\partial^2 Y}{\partial y^2}=k_2$，$\frac{1}{Z}\frac{\partial^2 Z}{\partial z^2}=k_3$然后我们可以对每一个方程分别求解得到：$X(x)=Ae^{\sqrt{k_1}x}+Be^{-\sqrt{k_1}x}$，$Y(y)=Ce^{\sqrt{k_2}y}+De^{-\sqrt{k_2}y}$，$Z(z)=Ee^{\sqrt{k_3}z}+Fe^{-\sqrt{k_3}z}$最终得到的总体解形式为：$u=\sum_{n=1}^{\infty} C_ne^{(-\sqrt{k_1^2+k_2^2+k_3^2})r}sin(n_1x)sin(n_2y)sin(n_3z)$2. 特征线法特征线法是一种常用于解决一阶偏微分方程的方法。

差分方程模型的理论和方法青岛建筑工程学院胡京爽引言1、差分方程：差分方程反映的是关于离散变量的取值与变化规律。

通过建立一个或几个离散变量取值所满足的平衡关系，从而建立差分方程。

差分方程就是针对要解决的目标，引入系统或过程中的离散变量，根据实际背景的规律、性质、平衡关系，建立离散变量所满足的平衡关系等式，从而建立差分方程。

通过求出和分析方程的解，或者分析得到方程解的特别性质（平衡性、稳定性、渐近性、振动性、周期性等），从而把握这个离散变量的变化过程的规律，进一步再结合其他分析，得到原问题的解。

2、应用：差分方程模型有着广泛的应用。

实际上，连续变量可以用离散变量来近似和逼近，从而微分方程模型就可以近似于某个差分方程模型。

差分方程模型有着非常广泛的实际背景。

在经济金融保险领域、生物种群的数量结构规律分析、疾病和病虫害的控制与防治、遗传规律的研究等许许多多的方面都有着非常重要的作用。

可以这样讲，只要牵涉到关于变量的规律、性质，就可以适当地用差分方程模型来表现与分析求解。

3、差分方程建模：在实际建立差分方程模型时，往往要将变化过程进行划分，划分成若干时段，根据要解决问题的目标，对每个时段引入相应的变量或向量，然后通过适当假设，根据事物系统的实际变化规律和数量相互关系，建立每两个相邻时段或几个相邻时段或者相隔某几个时段的量之间的变化规律和运算关系（即用相应设定的变量进行四则运算或基本初等函数运算或取最运算等）等式（可以多个并且应当充分全面反映所有可能的关系），从而建立起差分方程。

或者对事物系统进行划分，划分成若干子系统，在每个子系统中引入恰当的变量或向量，然后分析建立起子过程间的这种量的关系等式，从而建立起差分方程。

在这里，过程时段或子系统的划分方式是非常非常重要的，应当结合已有的信息和分析条件，从多种可选方式中挑选易于分析、针对性强的划分，同时，对划分后的时段或子过程，引入哪些变量或向量都是至关重要的，要仔细分析、选择，尽量扩大对过程或系统的数量感知范围，包括对已有的、已知的若干量进行结合运算、取最运算等处理方式，目的是建立起简洁、深刻、易于求解分析的差分方程。

偏微分方程的有限差分法及地球物理应用有限差分法是一种常用的数值求解偏微分方程的方法。

它将连续的偏微分方程转化为离散的差分方程，通过近似求解差分方程，得到偏微分方程的数值解。

这种方法在地球物理学中有着广泛的应用，如地震波传播模拟、电磁场分布计算等领域。

首先，假设我们要研究地震波在地下介质中的传播，可以采用波动方程来描述地震波的传播过程。

波动方程可以写成：∂^2u/∂t^2 = c^2∇^2u其中，u是地震波场，c是地下介质中的波速。

为了用有限差分法求解波动方程，我们需要将连续的空间和时间离散化。

假设我们将空间离散化为网格点（i,j,k），其中i,j,k分别代表空间的x,y,z方向，将时间离散化为时间步长Δt。

对波动方程进行近似，我们可以得到：(u(i,j,k,t+Δt) - 2u(i,j,k,t) + u(i,j,k,t-Δt))/Δt^2 = c^2(u(i+1,j,k,t) + u(i-1,j,k,t) + u(i,j+1,k,t) + u(i,j-1,k,t) +u(i,j,k+1,t) + u(i,j,k-1,t) - 6u(i,j,k,t))/Δx^2将此差分方程应用于地震波传播模拟，我们可以得到地震波场在空间和时间上的离散解。

有限差分法在地球物理中有着广泛的应用。

例如，它可以用于模拟地震波在地下介质中的传播，帮助研究地震灾害的发生机制和地下构造的特征。

通过调整网格的大小和时间步长，可以模拟不同频率的地震波传播过程，从而了解地震波在不同介质中的传播规律。

此外，有限差分法还可以应用于电磁场的计算。

例如，在电磁勘探中，可以利用有限差分法求解麦克斯韦方程，计算电磁场在地下介质中的传播和散射过程。

通过模拟电磁场的分布情况，可以帮助研究地下矿产资源的寻找和勘探。

需要注意的是，有限差分法在应用过程中还需要考虑边界条件的处理。

通常情况下，边界条件是已知的，例如地震波在地表的边界条件可以假设为自由表面，电磁场计算中的边界条件可以假设为电场和磁场的边界条件等。

有限差分法解偏微分方程综述绪论有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。

在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。

有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。

该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。

有限差分法以Taylor 级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。

对于有限差分格式，从格式的精度来划分，有一阶格式、二阶格式和高阶格式。

从差分的空间形式来考虑，可分为中心格式和逆风格式。

考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。

目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。

差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。

其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。

通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。

有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。

数值模拟偏微分方程的三种方法介绍（有限差分方法、有限元方法、有限体积方法）I．三者简介有限差分方法（Finite Difference Methods）是数值模拟偏微分方程最早采用的方法，至今仍被广泛使用。

该方法包括区域剖分和差商代替导数两个步骤。

首先将求解区域划分为差分网格，用有限个网格节点代替连续的求解区域。

其次，利用Taylor级数展开等方法将偏微分方程中的导数项在网格节点上用函数值的差商代替进行离散，从而建立以网格节点上的值为未知量的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且十分成熟的数值方法。

差商代替导数后的格式称为有限差分格式，从格式的精度来考虑，有一阶格式、二阶格式和高阶格式。

从差分的空间离散形式来考虑，有中心格式和迎风格式。

对于瞬态方程，考虑时间方向的离散，有显格式、隐格式、交替显隐格式等。

目前常见的差分格式，主要是以上几种格式的组合，不同的组合构成不同的差分格式。

差分方法主要适用于结构网格，网格的大小一般根据问题模型和Courant 稳定条件来决定。

有限元方法（Finite Element Methods）的基础是虚位移原理和分片多项式插值。

该方法的构造过程包括以下三个步骤。

首先，利用虚位移原理得到偏微分方程的弱形式，将计算区域划分为有限个互不重叠的单元（三角形、四边形、四面体、六面体等），在每个单元上选择合适的节点作为求解函数的插值点，将偏微分方程中的变量改写成由各变量或其导数的节点值与所选用的分片插值基函数组成的线性表达式，得到微分方程的离散形式。

利用插值函数的局部支集性质及数值积分可以得到未知量的代数方程组。

有限元方法有较完善的理论基础，具有求解区域灵活（复杂区域）、单元类型灵活（适于结构网格和非结构网格）、程序代码通用（数值模拟软件多数基于有限元方法）等特点。

有限元方法最早应用于结构力学，随着计算机的发展已经渗透到计算物理、流体力学与电磁学等各个数值模拟领域。

偏微分方程离散差分格式差分方法等偏微分方程（Partial Differential Equations，简称PDE）是一类涉及多个独立变量的微分方程，其中至少一个是时间变量。

这类方程广泛应用于物理、工程、金融等领域，解析解往往难以获得，因此需要使用数值方法进行求解。

差分方法是其中一种常用的数值方法，将连续的变量和算子替换为离散的差分近似，从而将偏微分方程转化为代数方程组求解。

差分方法的基本思想是将连续的自变量和函数替换为离散的自变量和函数。

设自变量x的取值范围是[a,b]，将其等分为N个点，即x_i=a+i·△x，其中△x=(b-a)/N。

常见的差分格式有前向差分、后向差分和中心差分。

下面以一维热传导方程为例，介绍差分方法的基本思想和常用格式。

一维热传导方程是描述材料温度分布随时间变化的方程，其数学表达式为∂u/∂t=a·∂²u/∂x²，其中u(x,t)表示温度分布，a是热传导系数。

为了使用差分方法求解该方程，我们需要将偏导数用近似的差分形式替代。

常用的差分格式是中心差分格式，我们将二阶导数的中心差分表示为(∂²u/∂x²)_i=(u_(i+1)-2u_i+u_(i-1))/△x²。

将此近似代入热传导方程，则可以得到u_i^(n+1)=u_i^n+a·△t/△x²·(u_(i+1)^n-2u_i^n+u_(i-1)^n)，其中u_i^n表示在x_i处、t_n时刻的温度，△t表示时间步长。

上述离散方程是一个差分方程，可以通过迭代计算求解。

首先，我们需要给定初始条件u(x,0)，即温度在初始时刻的分布。

然后，使用上述离散方程迭代计算下一个时间步的温度分布，直到达到所需的时间范围。

差分方法的稳定性和精度主要取决于离散精度和时间步长。

差分格式的离散精度决定了近似解和精确解之间的误差大小，一般而言，中心差分格式具有二阶精度。

偏微分方程的数值解法差分法
偏微分方程是描述自然现象和工程问题的重要数学工具。

它们出现在许多领域，如物理学、化学、工程学等。

由于解析求解偏微分方程的方法往往非常困难，因此需要数值方法来求解。

差分法是偏微分方程数值解法中的一种常用方法。

它的基本思想是通过将区域离散化为网格，将偏微分方程转化为离散化方程组。

然后使用迭代算法求解方程组，得到数值解。

差分法的主要优点是易于理解和实现。

通过选取不同的差分格式和网格划分方法，可以得到不同精度和稳定性的数值解。

此外，差分法还可以方便地处理不规则区域和非线性问题。

在使用差分法求解偏微分方程时，需要注意选择合适的网格划分和差分格式。

同时，还需要考虑数值解的稳定性和精度，以及计算效率等问题。

总之，差分法是求解偏微分方程的常用数值方法，对于解决实际问题具有重要的应用价值。

- 1 -。

差分知识点总结一、差分的概念差分是一种数学运算方法，用来计算函数在两个相近的点之间的变化量。

差分的基本思想是利用两个相近点之间的函数值的差来近似表示函数在这一区间的变化率。

差分主要应用在数值计算、微分方程数值解法、离散化微分方程和差分方程等领域。

二、差分的方法1. 前向差分前向差分是指用函数在点x和x+h处的函数值之差来近似表示函数在点x处的导数。

前向差分的公式为：f'(x) ≈ (f(x+h) - f(x)) / h2. 后向差分后向差分是指用函数在点x和x-h处的函数值之差来近似表示函数在点x处的导数。

后向差分的公式为：f'(x) ≈ (f(x) - f(x-h)) / h3. 中心差分中心差分是指用函数在点x+h和x-h处的函数值之差来近似表示函数在点x处的导数。

中心差分的公式为：f'(x) ≈ (f(x+h) - f(x-h)) / 2h4. 二阶中心差分二阶中心差分是指用函数在点x+h、x和x-h处的函数值之差来近似表示函数在点x处的二阶导数。

二阶中心差分的公式为：f''(x) ≈ (f(x+h) - 2f(x) + f(x-h)) / h^25. 前向差分法和后向差分法的优缺点前向差分法和后向差分法都是利用简单的迭代方式得到节点之间的差值。

前向差分法计算简单，但是会使误差更大；后向差分法计算较为繁琐，但是误差相对较小。

6. 应用差分方法广泛用于微分方程和差分方程的数值解法，离散化微分方程，数值积分等方面，其基本思想是用差分概念近似表示数学模型的微分和积分运算。

三、差分方法的误差分析1. 截断误差在差分近似计算中，由于只取有限个点的函数值，使得近似结果与真实结果之间存在一定的误差，这种误差称为截断误差。

2. 离散化误差差分方法中最主要的误差来源是离散化误差。

因为使用差分方法时，通常需要将连续的问题离散化为一个离散的问题，这个离散化的过程会使得结果与真实结果之间存在误差。

偏微分方程初值问题差分方法finite difference method of initial value problem of partial differential equation 一种求解偏微分方程初值问题的主要数值方法。

许多连续介质的运动过程都可表示成含时间的偏微分方程。

最简单的有双曲型的对流方程(1) 和抛物型的扩散方程(2)式中α和σ是常数。

当u的初始状态(设为t =0时的状态)给定后，常要研究这些过程在t>0后的演化,在数学上就是给定初值(3)后, 求微分方程在｜x｜<∞,t>0时的解。

这种问题叫作初值问题初值问题(1)、(3)的解为(4)初值问题的差分方法包含下列步骤和问题，先把问题的求解区域进行网格剖分，再在格子点上按适当的数值微分公式把问题中的微商换成差商，从而把微分问题离散化,得到差分格式,最后求出差分格式的数值解。

差分格式的解的存在性和惟一性,有时并不显然,需要论证。

解的求法和解法的数值稳定性也需要研究。

此外，还要估计差分问题的解与微分问题的解的差别，研究在网格步长趋于零时前者对后者的收敛性以及差分问题的解是否连续地依赖于初值，即稳定性的问题。

微分方程定解问题的离散化(差分格式的建立)以初值问题(1)、(2)为例,在-平面上作两族平行于坐标轴的直线对初值问题(1)、(2)的求解区域进行网格剖分,如图[初值问题(1)、(2)求解区域的网格剖分]所示（网格线之间的距离，也可以是不相等的）。

网格线的交点叫作格点，在格点(x+Δx，t+Δt)上，由数值微分公式及方程(1)，可得(1)的另一种表达式∆+∆) ]略去其中(2x t 项,即得一个差分方程。

原问题的解(,)自然不能满足这个差分方程用表示这差分方程的解，则适合(6)把微分方程(1)的充分光滑的解代∆+∆)，叫作差分方程的截断误差, 它对空间是二阶的,对时间是入式(6)，其误差为(2x t一阶的。

式(6)就是一个把微分方程(1)离散化而得的差分方程。

第9章 偏微分方程的差分方法含有偏导数的微分方程称为偏微分方程。

由于变量的增多和区域的复杂性，求偏微分方程的精确解一般是不可能的，经常采用数值方法求方程的近似解。

偏微分方程的数值方法种类较多，最常用的方法是差分方法。

差分方法具有格式简单，程序易于实现，计算量小等优点，特别适合于规则区域上偏微分方程的近似求解。

本章将以一些典型的偏微分方程为例，介绍差分方法的基本原理和具体实现方法。

9.1椭圆型方程边值问题的差分方法9.1.1 差分方程的建立最典型的椭圆型方程是Poisson （泊松）方程G y x y x f yux u u ∈=∂∂+∂∂-≡∆-),(),,()(2222 （9.1）G 是x ,y 平面上的有界区域，其边界Γ为分段光滑的闭曲线。

当f (x ,y )≡0时，方程（9.1）称为Laplace(拉普拉斯）方程。

椭圆型方程的定解条件主要有如下三种边界条件第一边值条件 ),(y x u α=Γ （9.2） 第二边值条件),(y x nuβ=∂∂Γ （9.3） 第三边值条件 ),()(y x ku nuγ=+∂∂Γ （9.4） 这里，n 表示Γ上单位外法向，α(x,y ),β(x,y ),γ(x,y )和k (x,y )都是已知的函数，k (x,y )≥0。

满足方程（9.1）和上述三种边值条件之一的光滑函数u (x ,y )称为椭圆型方程边值问题的解。

用差分方法求解偏微分方程，就是要求出精确解u (x ,y )在区域G 的一些离散节点（x i ，y i ）上的近似值u i ,j ≈(x i ，y i ）。

差分方法的基本思想是，对求解区域G 做网格剖分，将偏微分方程在网格节点上离散化，导出精确解在网格节点上近似值所满足的差分方程，最终通过求解差分方程，通常为一个线性方程组，得到精确解在离散节点上的近似值。

设G ={0<x <a , 0<y <b }为矩形区域，在x ,y 平面上用两组平行直线x =ih 1, i =0,1,…,N 1, h 1=a /N 1 y =jh 2, j =0,1,…,N 2, h 2=b /N 2将G 剖分为网格区域，见图9-1。