微分方程数值解法2.docVIP

微分方程数值解法 上机作业报告 年级专业：09信息与计算科学 指导老师： 谢树森 刘荣利： 150132009028 赵金鑫： 150132009063 徐弯 ： 150132009053 张润： 150132009061 王冠: 150132009042 赵腾: 150132009064 完成日期：2012-5-9 第三次上机作业 上机任务 用有限体积法求解如下possion方程 边界条件已知 令 边界条件为： 讲求得的解与精确解进行比较。 一次令N=4,8,16,32，取6位小数计算。 用消元法求解，并就处列出差分解与精确解。 就N=32，j=0.25,0.5,0.75及i=0,2,4,…,30,32画出差分解曲线。 算法介绍 采用五点差分格式，取和作矩形剖分，网格节点为 差分格式： 边界条件为： 实验结果 表示精确解，表示数值解，为误差，下同。 节点的排序顺序是按照横纵坐标和从小到大排列的。 N=4时 i j x y u u1 u1-u ...................................................................... 1.0000 1.0000 0.2500 0.2500 0.0253 0.0197 -0.0056 1.0000 2.0000 0.2500 0.5000 0.0000 0.0000 0.0000 1.0000 3.0000 0.2500 0.7500 -0.0253 -0.0197 0.0056 2.0000 1.0000 0.5000 0.2500 0.0358 0.0279 -0.0080 2.0000 2.0000 0.5000 0.5000 0.0000 0.0000 0.0000 2.0000 3.0000 0.5000 0.7500 -0.0358 -0.0279 0.0080 3.0000 1.0000 0.7500 0.2500 0.0253 0.0197 -0.0056 3.0000 2.0000 0.7500 0.5000 0.0000 0.0000 0.0000 3.0000 3.0000 0.7500 0.7500 -0.0253 -0.0197 0.0056 此时矩阵A为 N=8时 i j x y u u1 u1-u ...................................................................... 1.0000 1.0000 0.2500 0.2500 0.0253 0.0227 -0.0026 1.0000 2.0000 0.2500 0.5000 0.0000 -0.0000 -0.0000 1.0000 3.0000 0.2500 0.7500 -0.0253 -0.0227 0.0026 2.0000 1.0000 0.5000 0.2500 0.0358 0.0322 -0.0037 2.0000 2.0000 0.5000 0.5000 0.0000 -0.0000 -0.0000 2.0000 3.0000 0.5000 0.7500 -0.0358 -0.0322 0.0037 3.0000 1.0000 0.7500 0.2500 0.0253 0.0227 -0.0026 3.0000 2.0000 0.7500 0.5000 0.0000 -0.0000 -0.0000