陈锋茂[毕业论文]（精）

《陈锋茂[毕业论文]（精）》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、温M大學WENZHOUUNIVERSITY本科毕业论文（2012届）题目：有限差分方法在偏微分方程求解中的应用学院：数学与信息科学学院专业：数学与应用数学班级：08数本一姓名：陈锋茂学号：081093231022012年4月22日引言1有限差分格式11.1网格剖分11.2用Taylor级数展开方法建立差分格式31.3隐式差分格式61.4研究有限差分格式稳定性的Fourier方法81.4.1Fourier方法91.4.2判别准则122偏微分方程数值解152.1显性差分格式152.1.1差分格式的建立152.1.2稳定性分析152.1.3数值例子162.2半显性差分格式172.2.1差分格式的

2、建立172.2.2稳定性分析182.2.3数值算例193时滞偏微分方程组数值解213.1差分格式的建立213.2数值例子22致谢23参考文献24有限差分方法在偏微分方程求解中的应用陈锋茂08数木一摘要：木文主要研究冇限差分方法在偏微分方程求解小的应用，阐述了有限差分方法的基木原理及其求解步骤，并介绍了判断冇限差分格式稳定性的Fourier方法。对于常系数扩散方程，通过引入扩散项的方法，建立了一个3层条件稳定的显性差分格式，通过数值算例验证了该差分格式的可靠性和精确性。对丁•一维抛物型方程,在空间上利用四阶紧致差分逼近公式和加权平均思想，构造了一种高粘度无条件稳定的半显式差分格式，并JIJF

3、ourier分析方法证明了该差分格式是无条件稳定的。数值例子表明，所建立的差分格式是有效的。最后对于生态学上一个常见的时滞偏微分方程组，建立了半显式差分格式，对于具体的数值例子，给出了数值解的三维曲面图。关键字：冇限弟分方法，偏微分方程，数值解，冇限差分格式，时滞偏微分方程组引言计算机的飞速发展，正在日益影响着人们对传统数值分析（即计算方法）的认识。近几十年来，人们越来越认识到计算方法的学习与研究离不开计算机，仅靠数学理论的演绎和推导还不能完整地解决实际中的数值问题，要想研制出好的、实用的算法，必须要集合计算机科学。实践证明，计算方法止在口趋明显地成为数学与计算机科学的交叉性学科。科学计算

4、已经和理论计算、实验并列为三大科学方法⑷。偏微分方程定解问题在科学及工程计算中有着广泛的实际背景，人们可以通过偏微分方程对各种自然现象进行描述、解释和预测⑵。随着计算机科学技术的快速发展，现在已经口J以通过编制高效的程序对各种偏微分方程问题进行数值求解，偏微分方程数值解法在科学和工程计算中将起着越來越大的作用。诸如大型水坝应力分析、数值天气预报、飞行器设计的风洞实验等许多例子都可以说明求解偏微分方程的数值解在科学和各门工程中的应用是多么重要。对偏微分方程进行数值求解已经成为科学与工程计算的核心内容，当今最普遍的求偏微分方程数值解的方法主要有两类，分别为有限元方法和有限差分方法。目前，这两类

5、方法在应用上冇不同的侧重,有限元方法则侧重于定态问题（如椭I员I型方程），而有限差分方法则侧重于依赖时间的问题（如抛物型和双曲型方程）。本文主要研究冇限差分方法在偏微分方程求解中的应用。1有限差分格式1.1网格剖分因为有限差分方法是离散的，而偏微分方程问题是连续的，因此在用有限差分方法对偏微分方程问题进行数值求解时，必须要把连续的问题进行离散化。因此，就耍对所耍求解的区域进行网格剖分，求解区域随着所要求解的问题的不同而不同。下面就用儿个具体的例子来说明怎么对不同的求解区域给出网格剖分，在这个过程屮也我们也引入了一些常用的术语叫例1・1对于抛物型方程和双曲型方程的初边值问题，设其求解区域为D

6、}={（x,r）IOO}.首先我们在XT的上半平面作两族平行于坐标轴的直线，这样就把上半平面划分成了矩形网格。我们称这样的直线为网格线，其交点称为节点或网格点。就一般而言，平行于/轴的直线是可以等距的，设其距离为Ax>0,冇时也可以为力，称为空间步长。而大多数平行于兀轴的直线一般是不等距的，这里为了简便起见，一般地也假设它们是等距的，设其距离为△/〉()，有时也可以记为称其为时间步长。这样，这两族网格线就可以写作x=Xj=7Ax=jh,j=0,1,---,Jt=tn=nt=nr,n=0,1,2,…其中Xj=丿Ar=jh,Ax=/?=—;tn=nAt=nr。网格节点又口J简记

7、为(j,n)。口的网格J剖分见图1.Io例1.2对丁•抛物型方程和双曲型方程的初值问题，求解区域为D2={(x,t)I-oo0}.上述求解区域的网格剖分是曲平行于/轴的直线族x=x.=JAx=jh,j—0,±1,±2,---»与平行于兀轴的直线族t=tn=nt=nr,n=0,1,2,…所构成的，卩的网格剖分见图1・2。例1.3求解椭圆型方程的边值问题⑷。假设其求解区域为x-y平面上的一个有界区域D,其边