第四章偏微分方程的有限差分法PPT课件.pptxVIP

Harbin Institute of Technology Yangkun kyang@hit.edu.cn 1/75 第四章 偏微分方程的有限差分法 4.1 有限差分法原理 4.2 热传导方程的差分解法 4.3 波动方程的差分解法 Harbin Institute of Technology Yangkun kyang@hit.edu.cn 2/75 4.1 有限差分法原理 抛物线形 双曲型 椭圆形 不可逆过程 可逆过程 平衡过程 热传导方程 波动方程 位势方程 物理学中许多物理规律都用偏微分方程描述，偏微分方程主要分为以下三类： 上式中a,c,f以及未知函数u为定义在求解区域上的实（复）函数 Harbin Institute of Technology Yangkun kyang@hit.edu.cn 3/75 4.1 有限差分法原理 有限差分解法 差分近似代替微分，差商近似代替微商 这样就把求解区域内连续分布函数离散化成求网络节点上的分立函数值，从而把所需求解的微分方程变为一组相应的差分方程，进一步可以求解离散节点上的函数值。 数学基础 泰勒（Taylor）展开 Harbin Institute of Technology Yangkun kyang@hit.edu.cn 4/75 4.1 有限差分法原理 差商公式的构造 利用泰勒级数展开定义差商 Harbin Institute of Technology Yangkun kyang@hit.edu.cn 5/75 4.1 有限差分法原理 误差为O(h)差商公式： 一阶向前差商： 一阶向后差商： Harbin Institute of Technology Yangkun kyang@hit.edu.cn 6/75 4.1 有限差分法原理 二阶向前差商：式(2)-式(1)X2 误差为O(h2)差商公式： Harbin Institute of Technology Yangkun kyang@hit.edu.cn 7/75 4.1 有限差分法原理 二阶向后差商: 式(2)-式(1)X2 Harbin Institute of Technology Yangkun kyang@hit.edu.cn 8/75 4.1 有限差分法原理 一阶向前差商： Harbin Institute of Technology Yangkun kyang@hit.edu.cn 9/75 4.1 有限差分法原理 一阶向后差商： Harbin Institute of Technology Yangkun kyang@hit.edu.cn 10/75 4.1 有限差分法原理 一阶中心差商： Harbin Institute of Technology Yangkun kyang@hit.edu.cn 11/75 4.1 有限差分法原理 二阶中心差商： Harbin Institute of Technology Yangkun kyang@hit.edu.cn 12/75 4.1 有限差分法原理 差分格式的收敛性和稳定性 收敛性： 稳定性： 当步长h−→0时，差分方程的解趋向于微分方程的解。 误差在运算过程中不会失控，即累计误差不会无限增加。 Harbin Institute of Technology Yangkun kyang@hit.edu.cn 13/75 4.1 有限差分法原理 从数学上讲，没有限制的微分方程会有无穷多个解，不能构成一个定解问题。 从物理上讲，描述物理问题的微分方程仅适用于描述在一个连续体或物理场的内部发生的物理过程，仅靠这些微分方程不足以确定物理过程的具体特征。 因此，要