【毕业设计(论文)】二维热传导方程有限差分法的matlab实现.pdfVIP

第1章 前 言 1.1问题背景 在史策教授的《一维热传导方程有限差分法的MATLAB实现》和曹刚教授的《一维 偏微分方程的基本解》中，对偏微分方程的解得MATLAB实现问题进行过研究，但只停 留在一维中，而实际中二维和三维的应用更加广泛。诸如粒子扩散或神经细胞的动作 电位。也可 以作为某些金融现象的模型，诸如布莱克-斯科尔斯模型与 Ornstein-uhlenbeck过程。热方程及其非线性的推广形式也被应用与影响分析。 在科学和技术发展过程中，科学的理论和科学的实验一直是两种重要的科学方法 和手段。虽然这两种科学方法都有十分重要的作用，但是一些研究对象往往由于他们 的特性（例如太大或太小，太快或太慢）不能精确的用理论描述或用实验手段来实现。 自从计算机出现和发展以来，模拟那些不容易观察到的现象，得到实际应用所需要的 数值结果，解释各种现象的规律和基本性质。 科学计算在各门自然科学和技术科学与工程科学中其越来越大的作用，在很多重 要领域中成为不可缺少的重要工具。而科学与工程计算中最重要的内容就是求解科学 研究和工程技术中出现的各种各样的偏微分方程或方程组。 解偏微分方程已经成为科学与工程计算的核心内容，包括一些大型的计算和很多 已经成为常规的计算。为什么它在当代能发挥这样大的作用呢？第一是计算机本身有 了很大的发展；第二是数值求解方程的计算法有了很大的发展，这两者对人们计算能 力的发展都是十分重要的。 1.2问题现状 近三十年来，解偏微分方程的理论和方法有了很大的发展，而且在各个学科技术 的领域中应用也愈来愈广泛，在我国，偏微分方程数值解法作为一门课程，不但在计 算数学专业，而且也在其他理工科专业的研究生的大学生中开设。同时，求解热传导 方程的数值算法也取得巨大进展，特别是有限差分法方面，此算法的特点是在内边界 处设计不同于整体的格式，将全局的隐式计算化为局部的分段隐式计算。而且精度上 更好。 目前，在欧美各国MATLAB 的使用十分普及。在大学的数学、工程和科学系科，MATLAB 1 苏佳园：二维热传导方程有限差分法的MATLAB实现 被用作许多课程的辅助教学手段，MATLAB 也成为大学生们必不可少的计算工具，甚至 是一项必须掌握的基本技能。 在我国，MATLAB在各大专院校的应用日益普遍，许多专业已把MATLAB作为基本计 算工具。在科研机构和工业界，MATLAB正得到越来越广泛的应用。 MATLAB 具有强大的图形绘制功能，为科学计算和图形处理提供了很大的方便。我 们只需制定的绘图方式，再提供绘图数据，有程序指令就可以得到形象、直观的图形 结果。因此，近些年越来越多的人开始使用MATLAB来求解数值计算和图形处理技术， 我们也可以绘制出热传导方程数值解的二维、三维图形，从而可以更好的理解热传导 方程的意义。 1.3 问题解决 目前，对于求解偏微分方程有很多方法，但差分法和有限元离散法式主要解决问 题的两种方法。一般来说，用差分法来接偏微分方程，解得得结果就是方程的准确解 函数再借点上的近似值。而用变分近似的方法求解，是将近似解表示成有限维子空间 中基函数的线性组合。有限元法也是基于变分原理，由于选择了特殊的基函数，使它 能适用于一般的区域。这种基函数是与区域的剖分有关的，近似解 表示为基函数的u 线性组合，二线性组合中的系数，又是剖分节点上 或其导数的近似值。u 有关一维热传导方程的有限差分法求解的MATLAB实现，西安建筑科技大学的史策 教授已经解决，本文借鉴史老师的求解思想，对二维热传导方程进行转换，再对解法 编程实现，从而进一步对热传导方程进行探讨。二维热传导方程求解在现实生活中的 应用也更加广泛，所以有很好的现实意义。 2 第2章 预备知识 [8] 定义2.1 含有未知函数u(x ,x ,,x ,t)的偏导数的方程称为偏微分方程。 1 2 n [8] 定义2.2 方程 u  u  u