抛物型方程的高精度差分格式的构造与理论分析.pdfVIP

河南师范人学硕l。学位论文 第一章引言(差分方法概述) 第一章引言(差分方法概述) 现代的科学研究与工程技术中，大量的数学模型都可以用微分方程来描述。 实际上很多近代自然科学的基本方程本身就是微分方程。从微积分理论形成以 来，人们一直用微分方程来描述解释或预见各种自然现象，不断地取得了碌著的 成效。遗憾的是，绝大多数微分方程(特别是偏微分方程)定解问题的解都不能 以实用的解析形式来表示。这就产生了理论与应用的矛盾：一方面，人们列出了 反映客观现象的各种微分方程，建立了大量实用的数学模型；另一方面，人们又 无法得到这些方程的准确解，以定量地描述客观过程。 随着电子计算机的出现和发展，解决上述矛盾的一门学科——微分方程数 值方法得到了前所未有的发展和应用。虽然常微分方程数值方法的历史可以追溯 到18世纪，一些偏微分方程的数值方法也在本世纪出得到研究，但是，它们发 展成为一门理论上严谨，实用上有效的学科。还是20世纪50年代以来的事f可 参考文献11。特别是近半个世纪以来，随着其他学科和工程技术发展的需要，电 子计算机的诞生和飞速发展，使微分方程的数值方法成为不仅是数学学科，而且 是很多学科领域的一种重要研究手段和方法【可参考文献2】。时至今日，微分方 程的数值方法已有很多，但最基本和常用的就是差分方法和有限元法。 差分方法是求解微分方程定解问题的重要数值方法之一。1928年， 电子计算机的诞生与发展为差分方法提供了强有力的工具，促使这一学科迅速发 展起来，目前，对于线性偏微分方程定解问题，差分方法已经形成了较成熟的算 法格式与软件，并能对其进行较完善的理论分析，对于非线性问题，有效的算法 及理论分析正在迅速发展中【可参考文献4】。 用差分方法求解微分方程定解问题，主要是解决如下几个问题；首先是差 分格式的构造，即解决方程和定解条件的差分逼近问题。其次是差分格式的理论 分析，即讨论差分格式的相容性，稳定性及收敛性问题。三是差分格式的求解问 题，即讨论具有特殊结构系数矩阵的大型线性方程组的求解方法。三个问题既相 互联系义相互制约。只有较好的解决了上述三个问题，才能构造出好的差分格式。 本文将以第一问题为主线展丌论述。 河南师范大学硕十学位论文 第一二章差分格式的构造矗法概述 第二章差分格式的构造方法概述 构造差分格式已有很多方法，常用的有以下几种： 1．直接差化法 即用差商代替微商将微分方程定解问题离散成线性代数 方群组的方法．如古典显、隐格式，加权六点格式，c—N格式等．都是 用这种方法构造出来的。 2．积分插值法对守恒型微分方程进行积分，可以转化成积分方程或微 分积分方程，后者往往保持了前者的守恒性。在此基础上用数值积分 公式来代替积分，就可以得到差分方程。这种构造差分方程的方法叫 积分插值法，使用这种方法，便于将所得格式推广到任意网格和处理 自然边值条件。 3．待定系数法这也是构造差分格式的常用方法。使用这种方法时，先 给出含参数的差分方程或差分算子，然后，适当选择参数，使差分方 程或差分算子逼近所论的微分方程或微分算子的逼近阶尽可能的高， 且有较好的稳定性。 4．算子方法 导数的差商近似眭I可以借助线性算子的运算导出，首先根 据需要引入若干个差分算子，然后建立差分算子和导数算子之问的关 系，利川这些关系就可给出偏导数的差分表达式，用这些差分表达式 代替微分方程中的偏导数就可以得到逼近微分方程的各种差分格式 (方程)。 用差分方法求解微分方程定解问题，需要构造出精度高，计算量小，稳定 性好的差分格式。一般的构造方法，如上述前两种方法很难构造出这样的格式。 下而我们来论述使用后两种方法构造精度高，计算量小，稳定性尚好的差分格式。 下面我们柬论述使用后两种方法构造精度高，计算量小，稳定性尚好的差分格式。 2 河南师范大学硕士学位论文