偏微分方程数值解 引论1.pptVIP

* * * * * * §2 数值求解微分方程的意义 如果能找到一个（或一族）具有所要求阶连续导数的解析函数，将它代入微分方程（组）中，恰好使得方程（组）的所有条件都得到满足，我们就将它称为这个方程（组）的解析解（也称古典解）。“微分方程的真解”或“微分方程的解”就是指解析解。寻找解析解的过程称为求解微分方程。 微分方程的解在数学意义上的存在性可以在非常一般的条件下得到证明，这已有许多重要的结论。但从实际上讲，人们需要并不是解在数学中的存在性，而是关心某个定义范围内，对应某些特定的自变量的解的取值或是近似值-这样一组数值称为这个微分方程在该范围内的数值解，寻找数值解的过程称为数值求解微分方程。 §2 为什么要研究数值求解方法呢？ 1）在实际问题中我们所能获取的或感兴趣的，往往只是一个特定点上的数据。如空间的温度分布只能一个点一个点地测定，火箭升空传回的控制信息只能以某个确定的时间为间隔，一个个地发送和接受，如此等等。这些离散点上的函数值对于解决实际问题，已经足够了，寻找解析解的一般形式未必必要。 2）在很多情况下，寻找解析解也并无可能。现实问题中归结的微分方程不满足解析解的存在条件的比比皆是，方程中出现的有些函数连续性都无法保证，它们并不存在前述意义的解析解。于是，求数值解便成了在这种情况下解决问题的重要手段了。 3）即使微分方程的解析解存在，以并不意味可以将它表示为初等函数，如多项式、对数函数、指数函数三角函数及它们的不定积分的有限组合形式——显式解。 事实上，有显式解的微分方程只占解析解存在的微分方程中的非常小的一部分。 §3 数值求解方法概述 1）区域剖分 把整个定义域分成若干个小块，以便对每小块上的点或片求出近似值，这样按一定规律对定义域分切的过程称为区域剖分。 2）微分方程的离散 区域剖分完毕后，依据原来的微分方程去形成关于这些离散点或片的函数值的递推公式或方程。这是它们的未知量已不是一个连续函数，而成了若干个离散的未知值的某种组合了，这个步骤称为微分方程离散。 3）初始和边界条件处理 离散后系统是一个递推公式，那它需要若干个初值才能启动。若是一个方程组，那它所含的方程个数一般少于未知量的个数，要想求解还需要补充若干个方程。这些需要补充的初值和方程往往可以通过微分方程的初始条件和边界条件来得到，这就是初始和边界条件处理过程。 4）离散系统的性态研究 我们主要研究：这个系统是否可解，即解的存在性、唯一性问题；它与精确解的差距有多大，这个差距当区域剖分的尺寸趋于零时，是否也会趋于零，趋于零的速度多快，即解的收敛性和收敛速度问题；当外界对数据有所干扰时，所得的解是否会严重背离离散系统的固有的解，即解的稳定性问题。 上述问题说道底是一个误差分析问题，因为如果从实际问题到得出数值解的每一步都没有任何误差的话（当然，这是不可能的），那么数值解就应该是离散点上的精确值，也就不用煞费苦心去讨论上面的问题了。 * 1、模型误差 2、观测误差 3、截断误差 4、舍入误差 误差的来源 1）模型误差 在将实际问题归结为数学模型时需要对问题作一定的简化和假设，由此产生的误差叫模型误差。 2）观测误差 数学模型中需要用到的一些系数，初值等常数来自于测量仪器或统计资料，由于客观条件和仪器精度的限制不可避免有误差，这被称为观测误差。 3）截断误差 将数学模型离散化时由于舍弃一些次要的项而导致的模型问题真解与离散问题真解的误差称为截断误差。 4）舍入误差 在上机实际计算中，由于计算机对所运算的对象按计算机字长四舍五入而产生的最终计算解与离散问题真解的误差是舍入误差。 * 本书主要研究截断误差。一般来说，一个可用方法的截断误差不应该超过模型误差和观测误差 最后，我们可以将离散系统送到计算机上去实际计算了。 * 数值求解微分方程过程示意 微分方程 区域剖分 离散系统的 性态研究 递推计算或解线性代数方程组 微分方程离散 初始和边界条件处理 解的存在性、唯一性 解的收敛性和收敛速度 解的稳定性 得到数值解 求解过程中产生的误差 现实问题 数学模型 离散格式 模型误差 建模 离散 舍入误差 观测模型 截断误差 数值解 计算 * 参考书 1. 《微分方程数值解》，李荣华等编，高等教育出版社，1996 2.《偏微分方程差分方法》，陆金甫等编，高等教育出版社，1988 3.《微分方程数值解》，李立康等编，复旦大学出版社，1999年 4.《微分方程数值解》，余德浩等编，科学出版社，2003年 5. 《Numerical Partial Differential Equations》J.W.Thomas，Springer-Verlag,1997年 6. 《偏微分方程的MAT