边值问题的课案.docxVIP

西北农林科技大学理学院应用数学系 《微分方程数值解》结课论文 论文题目 边值问题的研究 2016年 1 月 14 日 一·问题重述 对于下列边值问题： 其中A为学号的倒数第2位，B为学号的倒数第1位。 (1)差分：截断误差、稳定性、收敛半径、递推（隐式）或方程组（显式） (2)有限元：刚度矩阵、算法步骤及代码 二·问题分析 题目明确指出使用差分方法和有限元解法。什么都不管先构造一种差分格式，然后对求解区域做划分将问题离散化，从微分方程的定解问题转化为求线性代数方程的解，以便于能够使用计算机进行计算。 在这里选用的是中心差分法，同时将边界进行处理，同时用Ritz有限元法和Galerkin法有限元法尝试去得到结果，最后再去比较两种解法所得到结果的精确性，分析相容性和截断误差等等。 三·解题过程 1·首先建立差分格式，考虑两点的边值问题，由题目知道 建立中心差分格式如下 对求解区间做网格划分，在a到b之间取N+1个节点，定义为xi（i取1到N） 即将区间I=[a,b]分为N个小区间 由此得到区间的一个网格剖分。记。用表示网格内点，，…，的集合，表示内点和界点的集合。 取相邻节点的中点，称为半整数点。由节点又构成的一个对偶剖分。 用差商代替微商，将方程（1.1）在内点离散化. 逼近边值问题（1.1）（1.2）的差分方程为： 当网格均匀，即时差分方程简化为 这相当于用一阶中心