一类反应扩散方程组的新的有限差分算法.pdfVIP

学位论文独创性声明 本人所呈交的学位论文是在我导师的指导下进行的研究工作及取得的研究成 果．据我所知，除文中已经引用的内容外，本论文不包含其他个人已经发表或撰写 的研究成果．对本文的研究做出重要贡献的个人和集体，均已在本文中作了明确的 说明并表示谢意． 作者签名 学位论文使用授权声明 本人完全了解华东师范大学有关保留，使用学位论文的规定，学校有权保留学 位论文并向国家主管部门或指定机构送交论文的电子版和纸质版．有权将学位论文 用于非赢利目的的少量复制并允许论文进入学校图书馆被查阅．有权将学位论文的 内容编入有关数据库进行检索．有权将学位论文的标题和摘要出版．保密的学位论 文在解密后适用本规定． 学位论文作者签名 导师签名：-=幽 第一章引言 反应扩散方程(组)是非线性抛物型方程(组)的一类重要的模型，它描述了生 物学中物种数量的迁徙变化，人体或动物等复杂的组织的发育形成过程，人体生理 学中的种种现象以及许多有趣的化学反应，具有深刻的实际背景．扩散是最普遍的 自然现象之一，例如在物理学、化学，生物学等领域中都存在大量的扩散现象．涉 及到扩散现象的实际问题大多可以用反应扩散方程(组)来描述．无论从数学理论还 是应用科学来看，对反应扩散方程的研究都具有重要的意义．特别是近几十年来， 对这类方程的研究吸引了国内外众多的数学工作者，并取得了令人瞩目的进展，大 大丰富了偏微分方程的理论，促进了相关学科的发展． 偏微分方程的数值解法在数值分析中占有重要的地位，很多科学技术问题的数 值计算大都涉及偏微分方程的数值解问题．近三十多年来，偏微分方程数值解的理 论和方法都有了很大的发展，而且在各个科学技术领域中应用也愈来愈广泛．偏微 分方程数值解的方法主要有两类t有限差分方法和有限元方法。这两类方法应用上 有不同的侧重．有限差分法主要集中在依赖于时间的问题(双曲型和抛物型方程)。 而有限元方法则侧重于定态问题(椭圆型问题)．本文讨论用有限差分法求解反应 扩散方程组，论文分为两部分，第一部分是一种基于Crank—Nicolson格式的推 广的方法，这类方法计算较简便；第二部分讨论的是改进的预测．校正格式。它比 传统的预测校正格式有更高的计算精度和稳定性． 对于微分方程初值问题的任何一种计算格式，都必须考虑它的稳定性和收敛 性，这是关系到算法成功与否的极其重要的问题．差分方程的稳定性实际上就是差 分方程的解能否连续依赖初始条件的问题．对一个不稳定的差分格式即便它的精度 再高，也是没有实用价值的．显然，如果微分方程本身没有适定性，即初始条件发生 一些小的变化时，微分方程本身的解就谬以千里，当然不能奢望有一种差分格式能 逼近它，因此，稳定性的讨论总是建立在定解问题的假设下进行的．差分格式的收敛 性，也就是当时间步长和空间步长趋于0时，差分方程的解能否趋于微分方程解的 问题．要使在任一层的差分解都能收敛于其相应的t处的微分方程的解，首先要求在 下一层的差分值精确等于微分方程解的假设下，往上一层计算产生的误差当h一0 华东师范大学硕士学位论文 时也趋于0．换句话说，差分方程收敛的前提是差分格式的截断误差在当h一0时 趋于0，也就是说差分格式应具有相容性．为方便起见，下面给出几个概念． 截断误差； —般地，微分方程的解不会精确地满足差分方程，将差分方程中的各个项。同时用 微分方程的解在相应点的值代入，利用泰勒展开，就会得到一个误差项，这个误差 项就是截断误差． 有限差分格式的相容性： 当h，r一0，截断误差T(％，“)一0，则称差分格式与定解问题是相容的． 有限差分格式的收敛性s 设u(x，t)是偏微分方程的解，t‘々是逼近这个偏微分方程的差分格式的”真解”． 这里所指的真解是在差分格式过程中，忽略了各种类型的误差，比如舍入误差等． 即是说求解差分格式的过程是严格精确的，我们称差分格式是收敛的．如果时间步 长r和空间步长h趋向于0时，e}=u(％，t々)一H}一0，即当时间步长和空间步 长趋于0时，差分格式的解逼近微分方程的解． 有限差分格式的稳定性t 利用有限差分格式进行计算是按时间层逐层推进的．如果考虑二层差分格式，那么 计算第n+l层上的值¨}时，要用到第11层上计算出来的结果值，而计算n层的值 时的