毕业论文
题 目 抛物型方程的差分解法 学 院 数学科学学院 专 业 信息与计算科学 班 级 计算0802 学 生 王丹丹 学 号 20080901045 指导教师 王宣欣
二〇一二 年 五 月二十五日
济南大学毕业论文
摘 要
偏微分方程的数值解法在数值分析中占有重要的地位，很多科学技术问题的数值计算包括了偏微分方程的数值解问题【1】。近三十多年来，数值解法的理论和方法都有了很大的发展，而且在各个科学技术的领域中应用也愈来愈广泛。本文的研究主要集中在依赖于时间的问题，借助于简单的常系数扩散方程，介绍抛物型方程的差分解法。本文以基本概念和基本方法为主，同时结合算例实现算法。
第一部分介绍偏微分方程及差分解法的基本概念，引入本文的研究对象——常系
?u?2u数扩散方程：?a2,x?R,t?0
?t?x第二部分介绍上述方程的几种差分格式及每种格式的相容性、收敛性与稳定性。 第三部分通过算例检验每种差分格式的可行性。
关键词：偏微分方程；抛物型；差分格式；收敛性；稳定性；算例
I - -
济南大学毕业论文
ABSTRACT
The 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 the 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 will 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 difference method.I will introduce the object-constant diffusion equation for the
?u?2ufirst time. ?a2,x?R,t?0
?t?xThe 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 differential equation；parabolic；difference scheme；convergence；stability；
application
II - -
济南大学毕业论文
目 录
摘要..............................................................………………..................I ABSTRACT………………………………………………………………………...…II 目录…………………………………………………………………………………..……III 1前言……….…………………………………………….….……………………….....….1 2基本概念和定理………………………………………………………………………….2
2.1抛物型方程的基本概念.....……………….………….…………………................2
2.1.1偏微分方程的定义……………………………………………………….…2 2.1.2抛物型方程的定义…………………………………………………….……2 2.1.3初边值条件的定义…………………………………………………….……3
2.2 差分方法的基本思想……………………………………………………..………3 2.3网格剖分...................................................................................................................4 2.4截断误差的基本概念...............................................................................................5 2.5相容性的基本概念...................................................................................................7 2.6收敛性的基本概念...................................................................................................7 2.7稳定性的基本概念...................................................................................................8
2.7.1判断稳定性的直接法…………………………………………………...…..8 2.7.2判断稳定性的Fourier方法……………………………………………....….9
3常系数扩散方程的差分格式及其相容性、收敛性和稳定性分析…….........................12 3.1向前差分格式…………………………………………………............................12 3.2向后差分格式........................................................................................................13 3.3 Crank-Nicolson格式..............................................................................................14 3.4 Richardson格式.....................................................................................................16 4差分解法的应用...............................................................................................................18 结论......................................................................................................................................25 参考文献..................................................... .................. .....................................................26 致谢......................………………….……………………..…….…………...…………….27 附录……………………………………………………………………………………...……………..28
III - -
济南大学毕业论文
1前言
微积分方程这门学科产生于十八世纪，欧拉在他的著作中最早提出了弦振动的二
阶方程，随后不久，法国数学家达朗贝尔也在他的著作《论动力学》中提出了特殊的偏微分方程[2]。偏微分方程得到迅速发展是在十九世纪，那时候，数学物理问题的研究繁荣起来了，许多数学家都对数学物理问题的解决做出了贡献。对于偏微分方程的求解，虽然具有明确表达式的解析解是很好的结果，但是能求出解析解的情况却十分有限，即使是最简单的双变量二阶线性常系数偏微分方程，也往往难以得到解析解。这是因为方程的解除了取决于方程本身的复杂度外，还要考虑到边界条件的复杂性[3]。很复杂的二阶偏微分方程，也许因为边界条件的简单性存在很简单的解析解，但是如果边界条件稍微复杂，就算是二阶常微分方程也没有解析解。
从目前的研究现状来看，偏微分方程数值解的理论和方法都日趋成熟，很多学术论文都在力求寻找更为精确且性质良好的求解方法。而且在实际问题中常会遇到多个自变量，非线性的方程或方程组；它们还可能是混合型的偏微分方程（如机翼的跨声速绕流），其解包含着各种间断（如激波间断、按触间断等）[4]。非线性问题的差分法求解是十分困难的。抛物型方程的数值解法目前有傅里叶算法(SSPE)、有限差分法(FDM)、有限元法(FEM)等，每种方法都有自己的适用范围，虽然SSPE 的效率高，但本文将选择使用更容易处理阻抗边界条件的FDM。由于FDM对网格间距要求足够小，计算效率很依赖计算机硬件速度，21 世纪前，大多是FDM的理论推导和误差分析，直到2007 年国际上才出现公开发表的使用FDM求解抛物型方程的实验并得出简单模型的计算结果。随着电子计算机的发展，在解决各种非线性问题中，FDM法得到了很快的发展，并且出现了许多新的思想和方法，如守恒差分格式，时间相关法、分步法等。
本文将从基本概念和基本方法入手，通过简单的常系数扩散方程，介绍抛物型方程的差分解法及其简单的实际应用，起到初步介绍偏微分方程数值解法的目的，希望有助于初学者了解相关基本知识，培养进一步学习的兴趣。
1 - -