LL用向前差分式计算初边值问题.docVIP

偏 微 分 方 程 课 程 设 计 学号：0683110 姓名：陆莉 指导老师：翟方曼 2010.01 一．题目 用向前差分格式计算如下热传导方程的初边值问题 已知其精确解为 二．理论 作为模型，考虑一维热传导方程： …………（1.1） 其中是正常数，是给定的连续函数。 现在考虑第二类初边值问题的差分逼近： 初始条件：…………（1.2） 边值条件：，，………（1.3） 假设和在相应区域光滑，并且在满足相容条件，使上述问题有惟一充分光滑的解。 1.建立差分格式 （1）.区域剖分 取空间步长和时间步长，其中都是正整数。用两族平行直线和将矩形域分割成矩形网格，网格节点为。以表示网格内点集合，即位于开矩形的网点集合；表示所有位于闭矩形的网点集合；是网格界点集合。 其次，用表示定义在网点的函数， （2）.微分方程的离散，建立相应差分格式 将方程在节点离散化， ， …………（1.4） 对充分光滑的解，由Taylor展式： …………（1.5） …………（1.6） …………（1.7） （1.5）移项得： …………（1.8） （1.6）（1.7）相加得： …………（1.9） 将（1.8）（1.9）代入（1.4）得： …………（1.10） 其中， 舍去，得到逼近（1.1）的向前格式差分方程： ， ……（1.11） 其中