常微分方程差分解法、入门、多解法

1、毕业论文题 目 抛物型方程的差分解法 学 院 数学科学学院 专 业 信息与计算科学 班 级 计算0802 学 生 王丹丹 学 号 20080901045 指导教师 王宣欣 二一二 年 五 月二十五日- 2 -济南大学毕业论文摘 要偏微分方程的数值解法在数值分析中占有重要的地位，很多科学技术问题的数值计算包括了偏微分方程的数值解问题【1】。近三十多年来，数值解法的理论和方法都有了很大的发展，而且在各个科学技术的领域中应用也愈来愈广泛。本文的研究主要集中在依赖于时间的问题，借助于简单的常系数扩散方程，介绍抛物型方程的差分解法。本文以基本概念和基本方法为主，同时结合算例实现算法。第一部分介绍偏微分方

2、程及差分解法的基本概念，引入本文的研究对象常系数扩散方程：第二部分介绍上述方程的几种差分格式及每种格式的相容性、收敛性与稳定性。第三部分通过算例检验每种差分格式的可行性。关键词：偏微分方程；抛物型；差分格式；收敛性；稳定性；算例ABSTRACTThe numerical solution of partial differential equation holds an important role in numerical analysis .Many problems of compution in the field of science and techology include th

3、e numerical solution of partial differential equation. For more than 30 years, the theory and method of the numerical computation made a great development and its applications in various fields of science and technology are more and more widely. This paper focuses on the problems based on time. I wi

4、ll use object-constant diffusion equation to introduces the finite difference method of parabolic equation. This paper mainly focus on the basic concept ,basic method and simple numerical example.The first part of this paper introduces partial differential equations and basic concepts of finite diff

5、erence method.I will introduce the object-constant diffusion equation for the first time. The second part of this paper introduces several difference schemes of the above equation and their compatibility ,convergence and stability.The third part tests the accuracy of each scheme.Key words:partial di

6、fferential equation；parabolic；difference scheme；convergence；stability；application目 录摘要.IABSTRACT.II目录.III1前言.12基本概念和定理.22.1抛物型方程的基本概念.22.1.1偏微分方程的定义.22.1.2抛物型方程的定义.22.1.3初边值条件的定义.32.2 差分方法的基本思想.32.3网格剖分.42.4截断误差的基本概念.52.5相容性的基本概念.72.6收敛性的基本概念.72.7稳定性的基本概念.82.7.1判断稳定性的直接法.82.7.2判断稳定性的Fourier方法.93常系数扩

7、散方程的差分格式及其相容性、收敛性和稳定性分析.12 3.1向前差分格式.12 3.2向后差分格式.133.3 Crank-Nicolson格式.143.4 Richardson格式.164差分解法的应用.18结论.25参考文献. . .26致谢.27附录.2829- -1前言微积分方程这门学科产生于十八世纪，欧拉在他的著作中最早提出了弦振动的二阶方程，随后不久，法国数学家达朗贝尔也在他的著作论动力学中提出了特殊的偏微分方程2。偏微分方程得到迅速发展是在十九世纪，那时候，数学物理问题的研究繁荣起来了，许多数学家都对数学物理问题的解决做出了贡献。对于偏微分方程的求解，虽然具有明确表达式的解析解是

8、很好的结果，但是能求出解析解的情况却十分有限，即使是最简单的双变量二阶线性常系数偏微分方程，也往往难以得到解析解。这是因为方程的解除了取决于方程本身的复杂度外，还要考虑到边界条件的复杂性3。很复杂的二阶偏微分方程，也许因为边界条件的简单性存在很简单的解析解，但是如果边界条件稍微复杂，就算是二阶常微分方程也没有解析解。从目前的研究现状来看，偏微分方程数值解的理论和方法都日趋成熟，很多学术论文都在力求寻找更为精确且性质良好的求解方法。而且在实际问题中常会遇到多个自变量，非线性的方程或方程组；它们还可能是混合型的偏微分方程（如机翼的跨声速绕流），其解包含着各种间断（如激波间断、按触间断等）4。非线性

9、问题的差分法求解是十分困难的。抛物型方程的数值解法目前有傅里叶算法(SSPE)、有限差分法(FDM)、有限元法(FEM)等，每种方法都有自己的适用范围，虽然SSPE 的效率高，但本文将选择使用更容易处理阻抗边界条件的FDM。由于FDM对网格间距要求足够小，计算效率很依赖计算机硬件速度，21 世纪前，大多是FDM的理论推导和误差分析，直到2007 年国际上才出现公开发表的使用FDM求解抛物型方程的实验并得出简单模型的计算结果。随着电子计算机的发展，在解决各种非线性问题中，FDM法得到了很快的发展，并且出现了许多新的思想和方法，如守恒差分格式，时间相关法、分步法等。本文将从基本概念和基本方法入手，

10、通过简单的常系数扩散方程，介绍抛物型方程的差分解法及其简单的实际应用，起到初步介绍偏微分方程数值解法的目的，希望有助于初学者了解相关基本知识，培养进一步学习的兴趣。2基本概念和定理2.1抛物型方程的基本概念2.1.1偏微分方程的定义偏微分方程在科学研究和工程技术中常常会出现，比如核反应和核爆炸过程的数学模型、飞行器设计过程中的空气动力学问题等等。定义2.1 含有未知函数的偏导数的方程称为偏微分方程。如果方程中对于未知函数和它的所有偏导数都是线性的，这样的方程称为线性偏微分方程，否则称它为非线性偏微分方程。方程中出现的偏导数的最高阶是方程的阶。下面举出几个典型的偏微分方程：（1）Laplace方

11、程5其中，该方程称为调和方程。在力学、电学常常遇到的势函数满足这个方程。（2）对流扩散方程其中表示流场中某种物质的浓度，是流速。（3）波动方程其中，而给定。在一些声学、光学和力学的波动问题中常常出现这类方程。2.1.2抛物型方程的定义定义2.2 设，其中可以是时间变量，二阶拟线性方程指，其中，和可以与有关，也可以与和有关。定义2.3 对于二阶拟线性方程，设矩阵是一个的矩阵，如果的特征值至少有一个为零，则该方程称为抛物型方程。考虑常系数扩散方程 （2.1） 其中是扩散过程中某种物质的浓度，是扩散系数。显然它是二阶线性方程，其中，它的矩阵为，的特征值为，所以方程（2.1）是一个抛物型方程。在下文中

12、我们将以该方程为基本模型讨论适用于抛物型方程定解问题的几种差分格式。2.1.3初边值条件的定义定义2.4 对于偏微分方程我们都是在一些特定条件下求方程的解，这样的条件称为定解条件。如果在的某个区域内求解方程，在的边界上给出的条件称为边界条件。在超平面给出的条件称为初始条件。给出了方程和定解条件，就构成了一个定解问题。按照定解条件的不同给法，可将方程（2.1）的定解问题分为两类。第一类，初值问题：空间变量的变化范围是。求具有所需次数偏微分的函数，满足（2.1）和初始条件第二类，初边值问题：空间变量的变化范围是。求具有所需次数偏微分的函数，满足（2.1）和初始条件及边界条件2.2 差分方法的基本思

13、想差分方法又称为有限差分方法或网格法，是求偏微分方程定解问题的数值解中应用最广泛的方法之一。它的基本思想是：先对求解区域作网格剖分，将自变量的连续变化区域用有限离散点（网格点）集代替；将问题中出现的连续变量的函数用定义在网格点上离散变量的函数代替；通过用网格点上函数的差商代替导数，将含连续变量的偏微分方程定解问题化成只含有限个未知数的代数方程组（称为差分格式）。如果差分格式有解，且当网格无限变小时其解收敛于原微分方程定解问题的解，则差分格式的解就作为原问题的近似解（数值解）。衡量一个差分格式是否经济实用，由多方面的因素决定，主要有：1、计算简单。显格式无需解方程组，计算较隐格式简单。但有些隐格

14、式左端的系数矩阵式与无关的三对角矩阵，相当于求解三对角矩阵系数的方程组，用消元法也是可取的。2、收敛性和收敛速度。当网格比固定，步长趋于零时，差分解应收敛到真解，并希望有尽可能快的收敛速度。3、稳定性。由于初始数据有误差，并不可避免的有舍入误差，因此必然要关心误差传递时是无限增长还是可以被控制。因此，用差分方法求偏微分方程定解问题一般需要解决以下问题：1、选取网格；2、对微分方程及定解条件选择差分近似，列出差分格式；3、求解差分格式；4、讨论差分格式解对于微分方程解的收敛性及误差估计。下面我们就按照这样的步骤来进行讨论。2.3网格剖分的概念用有限差分方法求解偏微分方程问题必须把连续问题进行离散

15、化，因此首先要对求解区域进行网格剖分6。在求解区域的平面画出两族平行于坐标轴的直线，把平面分成矩形网格。这样的直线称为网格线，其交点称为网格点或节点7。平行于轴的直线可以是等距的，可设距离为，记为，称为空间步长。平行于轴的直线则大多是不等距的，设为，记为，称为时间步长。下面用例子来说明不同区域的剖分。例2.1 双曲型方程和抛物型方程的初值问题，求解区域为：两族网格线可以写为：的网格剖分见图2.1. 图2.1例2.2双曲型方程和抛物型方程的初边值问题，设其求解区域是这个区域的网格由平行于轴的直线族与平行于轴的直线族所构成，其中。的网格剖分见图2.2 图 2.22.4有限差分格式的截断误差为叙述方

16、便，先引入下面的差分记号：向前差分向后差分中心差分二阶中心差分那么对于扩散方程（2.1）的解，关于的向前差分的Taylor级数展开为 （2.2）对变量进行Taylor级数展开为 （2.3）定义2.5 用微分方程的解来替代差分格式中的全部近似解，这样得到的方程两边的差就是截断误差8。如果一个差分格式的截断误差，则称差分格式对是阶精度，对是阶精度。例2.3 考虑下列差分格式的截断误差 ， （2.4）假定是充分光滑的，将带入上式的左侧,利用（2.2）式，（2.3）式有 这样就得到差分格式（2.4）的截断误差为。事实上，对于不在边界上的任何一点，可以定义（2.4）式的截断误差为，由截断误差的定义及以上

17、例子可知，只要网格剖分得很细，即和很小，那么偏微分方程的解近似地满足相应的差分方程。其实，一个有限差分格式的截断误差表示了用（偏微分方程的解）代替（差分方程的解）的差分方程与偏微分方程在点上的差。2.5有限差分格式的相容性从偏微分方程建立差分方程时，总是要求当时差分方程能与微分方程充分“接近”，这就导致了差分方程的一个基本特征，差分格式的相容性。考虑一般性的问题，设为微分算子，初值问题可以叙述为 （2.5）其差分格式可以统一写成 ， （2.6）其中是一个依赖于和的线性算子，称为差分算子，它把定义在第层上的函数变换到定义在第层上的函数。设（2.6）式为（2.5）式的差分格式，则相应的截断误差应是

18、定义2.6 设是定解问题（2.5）的充分光滑解，（2.6）式为求解问题（2.5）所得的差分格式，如果，当时有，则称差分格式（2.6）与定解问题（2.5）是相容的。由此可知，相容性表达了微分方程与差分方程间的关系。2.6有限差分格式的收敛性构造一个差分格式，它是否能在实际中应用，应该考虑当时间步长和空间步长无限缩小时，差分格式的解能否逼近微分方程问题的解，这就是差分格式的收敛性问题。定义2.7 设是偏微分方程的解，是逼近这个微分方程的差分格式的“真解”，所谓真解是指在求解差分格式的过程中忽略了各种类型的误差。如果当时间步长和空间步长趋向于0时，我们称差分格式是收敛的。应该注意的是，收敛性和相容性

19、是两个完全不同的概念。对于一个相容的差分格式，从定义判断其是否收敛是不容易的，通常需要寻求一些判别差分格式的收敛准则。不加证明地，我们给出如下定理。定理2.1（Lax等价定理） 给定一个适定的线性初值问题以及与其相容的差分格式，则差分格式的稳定性是差分格式收敛性的充分必要条件在应用中，差分格式的相容性是容易验证的，只要使截断误差趋于0就可以了。有了Lax等价定理，我们可以着重于差分格式的稳定性的讨论，一般不再讨论收敛性问题。差分格式一旦具有稳定性，就可以用差分格式计算出偏微分方程的近似解来。2.7有限差分格式的稳定性2.7.1判断稳定性的直接法利用有限差分格式进行计算时是按时间逐层推进的，计算

20、第层上的值要用到第层上的计算结果，而前一步的舍入误差必然会影响到的值。所以要分析这种误差传播的情况，我们希望误差的影响不会越来越大以至于掩盖差分格式的解的面貌，这就是所谓稳定性问题。差分格式可以统一写成（2.6）式，重复应用该式，有 为了度量误差及其他应用，引入范数 现在给出差分格式（2.6）的稳定性描述。定义2.8 设 有一个误差，则就有一个误差。如果存在一个正常数，使得当时，一致地有 （2.7）则称差分格式（2.6）是稳定的。 2.7.2判断稳定性的Fourier方法按差分格式稳定性的定义来直接验证其稳定性往往比较复杂，对于现行常系数偏微分方程初值问题可以用Fourier变换来进行求解和研

21、究。Fourier积分和Fourier变换是很多数学分支用到的有利工具，本文用于差分方法的分析。1、Fourier变换的定义：从Fourier级数出发进行推导可以得到Fourier积分，所谓Fourier积分公式，是指定义在上的函数的一个关系式，设，有关系式。令，则称为的Fourier变换9。2、线性常系数方程的Fourier方法由于解及初值只在网格点上有意义，为了应用Fourier方法讨论方程（2.1）的第一类初值问题，需要扩充这些函数的定义域，使它们在整个实轴R上都有定义。令这样，对任意都有定义。例2.4 以差分格式 （2.9）为例讨论Fourier方法。（2.9）式可以写为 （2.10）

22、对（2.10）式两边用Fourier积分，可以得到消去可得以上方法推广到一般形式的差分格式可以得到式中为增长因子，显然差分格式（2.9）的增长因子为。不加证明地，我们给出以下结论：差分格式稳定的充分必要条件是存在常数使得当时，有上式也称为von Neumann条件10。3、线性常系数方程组的Fourier方法稳定性的概念及相关Fourier方法的推导可以推广到线性常系数方程组。这在下面我们讨论Richardson差分格式时需要用到。一般的差分方程组可以写为其中。由于，即和满足一定关系，故在中仅标出。令则有 （2.11）其中称为差分算子，上式为一个差分格式。由（2.11）式可以得出如果不依赖于，

23、即为常系数差分方程组，则可利用Fourier积分得到，其中，称为增长矩阵。类似于差分方程的情况，我们有如下结论：差分格式（2.11）稳定的充要条件是存在常数使得当时，有定理2.2 差分格式（2.11）稳定的必要条件是当时，有，其中表示的特征值，为常数。定理2.311 当只有一个元素时，von Neumann条件是差分格式稳定的充要条件。对于Fourier方法，还有很多定理和推论可以帮助我们更加容易的判断不同性质的差分格式的稳定性，本文仅给出了需要用到的基本定理。3常系数扩散方程的差分格式及其相容性、收敛性和稳定性分析用Taylor级数展开方法求解偏微分方程实际上等价于用差商来近似微商得到相应的

24、差分格式，除此之外也可以用积分的方法。下面用Taylor级数展开的方法来讨论扩散方程初值问题 （3.1） （3.2）的差分格式。3.1向前差分格式假定偏微分方程初值问题的解是充分光滑的。由Taylor级数展开有 （3.3） 1、差分格式的推导过程利用（3.3）中第一式和第六式有：如果是（3.1）式的光滑解，即是满足的光滑函数，那么，方程（3.1）可以用如下差分方程来近似： ， （3.4） 这就是微分方程（3.1）的向前差分格式。它所联系的网点分布为：2、截断误差在前面的分析中，我们曾经利用这种差分格式来讨论截断误差的概念，得知其截断误差为。3、稳定性分析用Fourier方法来讨论其稳定性，容易

25、求出（3.1）式的增长因子是，其中，如果，那么有，即满足von Neumann条件，所以向前差分格式的稳定性条件是。这种情况的差分格式是条件稳定的。4、收敛性分析由Lax定理可知，当差分格式稳定时一定是收敛的，所以当时，向前差分格式收敛。3.2向后差分格式1、差分格式的推导过程上面构造的差分格式是显式的，即在时间层上的每一个可以独立地根据在时间层上的值得出。但是如果采用，和（3.3）中的六式，就可以得到扩散方程（3.1）的另一个差分格式 （3.5）这是微分方程（3.1）的向后差分格式。它所联系的网点分布为:向后差分格式与向前差分格式不同，在新时间层上包含了3个未知量，因此不能由直接计算出，这种

26、差分格式称为隐式格式。2、截断误差考虑其截断误差，其中。因此其截断误差为。3、稳定性分析先把差分格式变形为其中，令，并把它带入上面方程并消去公因子，容易求出（3.5）式的增长因子为由于，所以对任何网格比都有。由定理2.3知，差分格式（3.5）是稳定的。4、收敛性分析由Lax定理可知，当差分格式稳定时一定是收敛的，所以差分格式（3.5）收敛。3.3 Crank-Nicloson格式1、差分格式的推导过程将向前差分格式与向后差分格式作算术平均，即得Crank-Nicolson格式： （3.6）Crank-Nicolson格式又叫六点对称格式12，因为这种格式所联系的网点分布为:Crank-Nico

27、lson格式也是隐式格式。2、截断误差令将截断误差在处展开，则得因此其截断误差为。3、稳定性分析仍然用Fourier方法。先把差分格式变形为令，带入上面方程并消去公因子，容易求出（3.6）式的增长因子为由于，显然对任何网格比都有，由定理2.3知，差分格式（3.6）是稳定的。4、收敛性分析由Lax定理可知，当差分格式稳定时一定是收敛的，所以差分格式（3.6）收敛。3.4 Richardson格式1、差分格式的推导过程利用（3.3）式中的第二式和第六式可以得到Richardson格式： （3.7）从（3.7）式可以看出，这种格式联系到三个时间层，称为三层格式。一般地，一个多于二层的差分格式成为多层

28、差分格式。它所联系的网点分布为：2、截断误差这个格式是为了提高差分方程的截断误差所设计的，利用上述方法我们可以很容易的求出其截断误差为。单从截断误差这一角度来考虑，这个格式是比较好的，但是以后的分析可以看到，它是没有实用价值的。3、稳定性分析首先将（3.7）式变形为讨论多层差分格式，一般先化成与其等价的二层差分方程组。Richardson差分方程的等价二层差分方程组为如果令，则上面的方程组可以写成设，将其带入上式并消去公因子，可以得因此增长矩阵为其特征值为取则有由定理2.2可知Richardson格式是不稳定的。这种情况的差分格式称为绝对不稳定的。4、Richardson格式的改进1953年D

29、u Fort和Frankle对Richardson格式进行了修改，提出了如下的DuFort-Frankle格式13这是一个无条件稳定的差分格式，但是其相容性是有条件的。在这里对DuFort-Frankle格式我们不做过多讨论。事实证明，我们无法构造出无条件相容和无条件稳定的显示格式。4差分解法的应用下面我们用简单的例子来考察上述差分格式的数值结果。例4.1 考虑扩散方程14 （4.1）已知其解析解为取，为时间步长，为网格比。解：由于向前差分格式稳定性的条件为，所以当时，必须取，在计算时我们取。下面我们分别采用向前差分格式、向后差分格式、Crank-Nicloson格式和Richardson格式

30、进行求解，并将结果在与处与解析解进行比较。（1）向前差分格式为写成便于计算的格式为 （4.2）初边值条件的格式为由于向前差分格式为显格式，所以可以逐步求解取，利用初值和边值可以算出第一层的值。取，利用和边值，通过式（4.2）可以算出。如此下去，可以逐层计算出所有。计算结果如表4.1：=0.005=0.01=0.02=0.03=0.04=0.05=0.06=0.07=000000000=0.10.2938930.2795080.2528180.2286760.2068390.1870880.1692230.153063=0.20.5590170.5316570.4808880.4349670.3

31、934320.3558620.321880.291144=0.30.7694210.7317630.6618860.5986810.5415120.4898020.443030.400725=0.40.9045080.8602390.7780930.7037920.6365860.5757970.5208140.47108=0.50.9510560.9045080.8181360.7400110.6693460.6054290.5476160.495323=0.60.9045090.8602390.7780930.7037920.6365860.5757970.5208140.47108=0.

32、70.7694210.7317630.6618860.5986810.5415120.4898020.443030.400725=0.80.5590170.5316570.4808880.4349670.3934320.3558620.3218810.291144=0.90.2938930.2795090.2528180.2286760.2068390.1870880.1692230.153063 表4.1 当时，数值解与解析解比较如表4.2：时刻数值解解析解误差0.0050.5590170.559483-0.0004660.010.5316570.532544-0.0008870.020.4

33、808880.482495-0.0016070.030.4349670.437149-0.0021820.040.3934320.396065-0.0026330.050.3558620.358842-0.002980.060.321880.325117-0.0032370.070.2911440.294562-0.003418 表4.2当时，数值解与解析解比较如图4.1：图4.1从结果来看，当时，向前差分格式是稳定的，而且误差符号是相同的。由于求解是逐层进行的，当取不同的值时，误差都是随着时间层的增大而增大的，并且通过分析可知，每个网点上的误差对于以后各层网点上的影响都不差过其绝对值。所以取

34、进行比较是比较有代表性的，这在以下各个方法的分析中也是适用的。这个格式的优点是简单、易算，但由于截断误差为，精度不算高，又它仅当时才收敛和稳定，收敛的阶为的一阶，所以想要算得略为精确些，就要缩小和。注意到最大为，于是缩小一半，得缩小成四分之一，这就大大增加了工作量。（2）向后差分格式为写成便于计算的格式 （4.3）初边值条件的格式为由于向后差分格式为隐格式，不能逐步求解。当时，由式（4.3）可以得到若干个九元的线性方程组，该方程组的系数矩阵严格对角占优，方程可解，本文采用追赶法求解。所求结果如表4.3：=0.005=0.01=0.02=0.03=0.04=0.05=0.06=0.07=0000

35、00000=0.10.2896160.2692070.2332180.2184150.1938210.1737980.1567280.14176=0.20.5404280.4975950.4325950.4072250.3644780.3284680.2969910.269007=0.30.6965280.6403170.572880.5452940.4945210.4484520.406750.369057=0.40.627650.6706160.6524630.6281910.5750390.5234960.4758580.432366=0.50.8630160.7868470.69666

36、40.6625460.6021180.5478570.4983710.453203=0.60.8244160.7507410.6618920.6286630.5706890.5192420.4725210.42991=0.70.5325320.5672860.5493910.5283240.4833750.4404870.4010340.364982=0.80.4966960.4533390.4048870.3858490.3511470.3196570.2909690.264824=0.90.2786830.2526760.2176710.2052980.1851170.1680730.15

37、28830.139126 表4.3当时，数值解与解析解比较如表4.4时刻数值解解析解误差0.0050.5404280.559483-0.0190550.010.4975950.532544-0.0349490.020.4325950.482495-0.04990.030.4072250.437149-0.0299240.040.3644780.396065-0.0315870.050.3284680.358842-0.0303740.060.2969910.325117-0.0281260.070.2690070.294562-0.025555 表4.4当时，数值解与解析解比较如图4.2：图4

38、.2从结果来看，向后差分格式是稳定的，它的最大优点在于稳定性、收敛性对于一切网格比都成立，网格比增大意味着工作量减少，所以隐式格式当时的工作量比显示格式要少，然而从截断误差的角度来看，也不能太大，否则太大，误差也就太大，所以更高精度的格式也是十分需要的。（3）Crank-Nicloson格式为 写成便于计算的格式为 （4.4）初边值条件的格式为从式（4.4）可以看出，Crank-Nicloson格式也是一个隐格式，不能逐步求解。与向后差分方法类似，由式（4.4）可以得到若干个九元的线性方程组，该方程组的系数矩阵严格对角占优，方程可解，同样采用追赶法求解。所求结果如表4.5：=0.005=0.0

39、1=0.02=0.03=0.04=0.05=0.06=0.07=000000000=0.10.2942540.2801960.2540630.2303680.2088820.1894010.1717360.155719=0.20.5597040.5329650.4832570.4381860.3973180.3602610.3266610.296195=0.30.7703670.7335630.6651470.6031110.5468610.4958570.4496110.407677=0.40.905620.8623550.7819270.7090.6428740.5829150.52854

40、90.479253=0.50.9522260.9067340.8221660.7454860.6759580.6129130.5557490.503917=0.60.9056210.8623550.7819270.7090.6428740.5829150.5285490.479253=0.70.7703670.7335630.6651470.6031110.5468610.4958570.4496110.407677=0.80.5597040.5329650.4832570.4381860.3973180.3602610.3266610.296195=0.90.2942540.2801960.

41、2540630.2303680.2088820.1894010.1717360.155719 表4.5当时，数值解与解析解比较如表4.6：时刻数值解解析解误差0.0050.5597040.5594830.0002210.010.5329650.5325440.0004210.020.4832570.4824950.0007620.030.4381860.4371490.0010370.040.3973180.3960650.0012530.050.3602610.3588420.0014190.060.3266610.3251170.0015440.070.2961950.2945620.00

42、1633 表4.6当时，数值解与解析解比较如图4.3：图4.3从结果来看，Crank-Nicloson格式是稳定的，误差也比较小，收敛阶是的二阶，且随着时间层的增长误差逐渐增大，这是因为在求解方程组的过程中，各方程之间的表达式是相互关联的，前面方程的解对后面的解产生影响，而本身数值解法是存在误差的，所以随着求解过程的深入，误差会越来越大。但是总体来看，绝对误差逐渐增大，相对误差的变化幅度还是比较小的。（4）Richardson格式为写成便于计算的格式为 （4.5）初边值条件的格式为分析（4.5）式可知，Richardson格式是一个三层显示格式，能够逐层进行计算。令，我们可以看到的计算不仅需要

43、初值，还需要，在这里，我们引用向前差分格式中计算所得的。计算结果如表4.7：=0.005=0.01=0.02=0.03=0.04=0.05=0.06=0.07=000000000=0.10.2938930.2802480.2541590.2304330.2079480.174759-0.02023-2.27598=0.20.5590170.5330650.483460.4385880.3992320.3804670.5747993.317385=0.30.7694210.73370.6654080.6034280.5464840.4862750.339571-0.64845=0.40.9045080.8625170.7822410.7094320.6429270.5744310.376832-2.07622=0.50.9510560.9069050.8225090.7461150.6785190.6406880.9480076.360319=0.60.9045090.8625160.7822290.7092960.6412270.5516960.059067-6.64484=0.70.7694210.7337010.665420.6035920.548