数学偏微分方程的数值解法

偏微分方程的数值方法

偏微分方程（Partial Differential Equations，简称PDEs）是数学中研究的重要分支，广泛应用于物理学、工程学等领域中。由于一些复杂的PDEs难以找到解析解，因此需要借助数值方法进行求解。本文将介绍偏微分方程的数值解法，包括有限差分法、有限元法和谱方法等。

一、有限差分法（Finite Difference Method）

有限差分法是解偏微分方程最常用的数值方法之一。它将偏微分方程中的导数用差商来近似，将空间离散成若干个小区间和时间离散成若干个小时间步长。通过求解离散化后的代数方程，可以得到原偏微分方程的数值解。

以二维的泊松方程为例，偏微分方程可以表示为：

∂²u/∂x² + ∂²u/∂y² = f(x, y)

其中，u(x, y)为未知函数，f(x, y)为已知函数。我们可以将空间离散成Nx × Ny个小区间，时间离散成Nt个小时间步长。利用中心差分法可以近似表示导数，我们可以得到离散化的代数方程组。

二、有限元法（Finite Element Method）

有限元法是一种重要的数值解PDEs的方法。它将求解区域离散化成一系列的单元，再通过插值函数将每个单元上的未知函数近似表达。然后，利用加权残差方法，将PDEs转化成代数方程组。 在有限元法中，采用形函数来近似未知函数。将偏微分方程转化为弱形式，通过选取适当的形函数和权函数，可以得到离散化后的代数方程组。有限元法适用于求解各种各样的偏微分方程，包括静态和动态、线性和非线性、自由边界和固定边界等问题。

三、谱方法（Spectral Method）

谱方法是一种基于特殊函数（如正交多项式）的数值方法，用于解PDEs。谱方法在求解偏微分方程时，利用高阶连续函数拟合初始条件和边界条件，通过调整特殊函数的系数来近似求解解析解。

谱方法具有高精度和快速收敛的特点，适用于各种偏微分方程求解。与有限差分法和有限元法相比，谱方法的误差减小速度更快。然而，谱方法在处理非线性问题时可能面临困难。

第十章 偏微分方程数值解法

偏微分方程问题，其求解十分困难。除少数特殊情况外，绝

大多数情况均难以求出精确解。因此，近似解法就显得更为重要。本章仅介绍求解各类典型偏微分方程定解问题的差分方法。

§1 差分方法的基本概念

1.1 几类偏微分方程的定解问题

椭圆型方程：其最典型、最简单的形式是泊松（Poisson）方程

),(2222yxfyuxuu

特别地，当0),(yxf时，即为拉普拉斯（Laplace）方程，又称

为调和方程

02222yuxuu

Poisson方程的第一边值问题为

),(),(),(),(),(2222yxyxuyxyxfyuxuyx

其中为以为边界的有界区域，为分段光滑曲线，

称为定解区域，),(yxf，),(yx分别为，上的已知连续函数。

第二类和第三类边界条件可统一表示为

),(),(yxuuyxn

其中n为边界的外法线方向。当0时为第二类边界条件，

0时为第三类边界条件。

抛物型方程：其最简单的形式为一维热传导方程

220(0)uuaatx

方程可以有两种不同类型的定解问题：

初值问题

xxxuxtxuatu)()0,(,0022

初边值问题

221200,0(,0)()0(0,)(),(,)()0uuatTxltxuxxxlutgtultgttT

其中)(x，)(1tg，)(2tg为已知函数，且满足连接条件 )0()(),0()0(21glg

一、 实验目的

1、 了解并掌握不同差分格式的稳定性；

2、 能够掌握比较不同差分格式数值效应的能力。

二、 实验问题

取a=1,2,4, h=0.1, 0.08,用以下几种差分格式求解对流方程

01,0(0,)()0,0txuauxuxfxx

得t=4时数值结果。用图示说明算法的稳定性和间断点附近的计算效果，并进行相应的数值分析。

迎风格式(upwind):

111()nnnnjjjjuuauu

Lax-friedrichs 格式：

1111111()()22nnnnnjjjjjuuuauu

Lax-wendroff 格式：

122111111()(2)22nnnnnnnjjjjjjjuuauuauuu

修正迎风格式:

11(1)nnnjjpjpududu

这里p=[ a],d= a-[ a],/0.8h.为网格比，记号[x]表示不超过x的最大整数。

三、 实验原理

首先取[10,10]x，[0,4]t。按照h=0.1, 0.08划分网格。

再由各种差分格式通过已知的第1层网格点数值可以求出第2层网格点的数值。以此类推，通过逐层的信息最终求得在第4层网格点数值结果。

注意：这里x，t的取值范围应当包含间断点。同时，在所需求的第4层也应当包含间断点。这点要求可以通过初始估算得出。

四、 实验过程

根据要求将网格中x划分为200格，t划分为50格。首先通过matlab构造一个大的零矩阵。然后分别利用差分格式代入并作图。

1、 取a=1

迎风格式(upwind):

Lax-friedrichs 格式：

Lax-wendroff 格式：

修正迎风格式：

2、 取a=2

简析《偏微分方程数值解》教学改革

偏微分方程数值解是本校理学院数学系一门专业课,主要分为有限差分法和有限元法两部分。它为工科专业课程学习和解决实际问题提供必要的数学基础知识及常用的数值解法。随着科技的发展,偏微分方程数值解法在很多领域有着举足轻重的作用,比如在石油、医学以及航天等领域,因而这门课也适用于很多工科专业的学生。但是给数学系的学生讲解这门课,常常使人面临两难的境地:一方面数学专业的学生专业课多数侧重于数学理论,如果把偏微分方程数值解完全上成理论课,倒是迎合了学生一贯的思维模式,但应用方面就会欠缺;如果大量介绍具体应用,由于缺乏相关背景知识,又会使学生迷茫。所以在讲解时,需要把理论和实际合理地结合,即要把理论介绍得深入浅出,又要把某些算法与实际问题联系起来,从而使学生在学习中产生兴趣,是值得思考的问题。另一方面从以往的教学经验来看,偏微分方程数值解的教学存在诸多问题,诸如课时少,内容多;授课一直以板书为主,教学手段比较单一;以往的教学过于偏重理论,学生学过之后不知如何应用;由于课时紧张,而没有安排上机实践等问题。针对上述问题,为了提高教学质量,激发学生学习兴趣,笔者采取了一系列的教改措施,总结如下。

1 课程改革内容

1.1 改革教学内容,选择合适的教材和教辅丛书

关于偏微分方程数值的教材有很多,但一些教材在推导公式时使用的语言晦涩,公式抽象,学生难于理解。在我校的本科生授课时,选择了清华大学陆金甫编写的《偏微分方程数值解》[1]作为教材,以及北京大学出版的李治平编写的《偏微分方程数值解讲义》[2]作为教辅。《偏微分方程数值解》在内容的处理上,体现了由浅入深、循序渐进的原则;在叙述表达上,严谨精练、清晰易读,便于教学与自学。此教材充分反映偏微分方程数值解法中的核心内容。《偏微分方程数值解讲义》的作为教辅的特点是:每章之后配置了相当数量的习题,并在书后附上了大部分习题的答案或提示,更便于学生复习、巩固、理解和拓广所学的知识。