微分方程数值解 西北师范大学.doc

西北师范大学信息与计算科学专业 专业选修课程教学大纲 微分方程数值解 一、说明 （一）课程性质 《微分方程数值解》课是计算科学方向三年级专业限选课, 是一门基本理论与实际应用相结合的课程。主要讲授常微分方程初值问题与椭圆型偏微分边值问题的数值解法。 （二）教学目的 使学生掌握微分方程数值方法中的一些基本概念和基本理论，培养学生具有一定的理论分析能力，并通过上机实习，培养学生具有一定实际解题能力。 （三）教学内容 常微分方程初值问题的数值解法、变分原理、椭圆形方程的有限差分方法、椭圆形方程的有限元方法、离散方程的解法。 （四）教学时数 90学时。理论课周3学时，共54课时。上机实习周2学时，共36学时。 （五）教学方式 课堂讲授与上机实习相结合。 二、本文 理论部分 第一章 常微分方程初值问题 教学要点： Euler折线法、线性多步法、Runge-Kutta法、收敛性与误差估计。 教学时数： 10学时. 教学内容： 第一节 引论（2学时） 近似解析方法，数值解法，适定性，Bellmann不等式。 第二节 Euler折线法（2学时） Euler折线法计算公式，Euler法的局部截断误差与整体截断误差估计，Euler方法的稳定性与收敛性。 第三节 一般单步方法（2学时） Taloy级数方法，一般q阶单步方法，Runge-kutta方法。 第四节 线性多步方法（2学时） Adams外插法，Adams内插法，线性q阶k步方法。 第五节 稳定性、收敛性与误差估计（2学时） 线性差分方程，局部截断误差，相容性，稳定性，收敛性与误差估计。 考核要求： 掌握Euler折线法，Taloy级数方法与Runge-kutta方法, 会对给定的初值问题利用这些方法进行数值计算，并且会对误差进行估计。了解线性q阶k步方法。 第二章 变分原理 教学要点： 二次函数的极值，两点边值问题，广义导数，Sobolev空间，极小位能原理，二阶椭圆边值问题，Ritz-Galerkin方法。 教学时数： 8学时 教学内容： 第一节 二次函数的极值（2学时） 代数方程组与二次函数极值点的关系。 第二节 两点边值问题（2学时） 弦的平衡方程，广义导数，Sobolev空间Hm(I), 极小位能原理，虚功原理。 第三节 二阶椭圆边值问题（2学时） Sobolev空间Hm(G)，Poisson方程第一边值问题的能量泛函，广义解，自然边界条件， Poisson方程混合边值问题的变分方程。 第四节 Ritz-Galerkin方法（2学时） 无限维极值问题的有限维逼近，Ritz方法，Galerkin方法。 考核要求： 掌握两点边值问题与二阶椭圆边值问题的极小位能原理与虚功原理，能写出其能量泛函与变分方程，理解本质边界条件，理解Ritz-Galerkin方法的基本思想。 第三章 椭圆型方程的有限差分法 教学要点： 差分逼近的基本概念，一维边值问题的插分格式，二维边值问题的五点差分格式与三角网的差分格式，差分方程的极值定理与比较定理。 教学时数： 12学时 教学内容： 第一节 差分逼近的基本概念（2学时） 节点，步长，截断误差，差分方程，相容条件，差分方程的解的收敛性与稳定性，差分方程解的先验估计。 第二节 一维差分格式（2学时） 一维边值问题的直接插分格式，积分差分格式，变分-差分方法，及边界条件的处理。 第三节 矩形网的差分格式（2学时） 平面区域上Poisson方程的五点差分格式，极坐标形式的差分格式与边界条件的处理。 第四节 三角网的差分格式（2学时） 平面区域的三角陪剖分，三角陪剖分内点差分方程，三角陪剖分边界点差分方程。 第五节 极值定理（2学时） 差分方程的极值定理，解的唯一性与比较定理，五点格式的收敛速度估计。 第六节 能量不等式（2学时） 先验估计的能量不等式法，差分方程解的存在唯一性及收敛速度估计。 考核要求： 掌握一维边值问题的插分格式构造及对边界条件的处理，掌握平面区域上Poisson方程的五点差分格与三角陪剖分差分格式。掌握差分方程的极值定理。 第四章 椭圆型方程的有限元法 教学要点： 一维问题的线性元法，一维问题的高次元法，二维问题的矩形元法，二维问题的三角形元法，有限元方程的形成与系数矩阵的计算。 教学时数： 12学时 教学内容： 第一节 解一维问题的线性元（2学时） 试探函数空间，从Ritz方法的构造有限元方程，从Galerkin出发