重调和方程的二阶边值问题的求解.docVIP

微分方程数值解课程设计 ——重调和方程的二阶边值问题的求解 学院、系： 理学院 数学系 专 业： 信息与计算科学 姓 名： 陈剑宇 学 号： 201030470140 任课教师： 黄凤辉2012/12/27 总评成绩： 摘 要 本次课程设计，主要讨论重调和方程二阶边值问题的求解。文章分成四部分：第一部分，介绍如何将重调和方程的二阶边值问题分解，进而用经典的五点差分格式逼近；第二部分，列出对应的MATLAB算法和流程图；第三部分，通过上述方法求解具体的问题，并分析方法的精度和收敛阶；第四部分，将方法推广应用到抛物型方程，然后分析稳定性。 关键词：重调和方程的二阶边值问题；五点差分格式；MATLAB；抛物型方程 目 录 一 引言 4 1.1 背景 4 1.2 问题的提出 5 二 方法介绍 6 2.1 五点差分方法 6 2.2 二阶边值问题分解 8 2.2.1 当a2?4b≥0时 8 2.2.2 当a2?4b0时 9 2.3 推广到抛物型方程 10 三 算法程序 11 3.1 程序列表 11 3.2 算法流程图与介绍 12 3.2.1 重调和方程二阶边值程序 12 3.2.2 抛物型方程程序 14 四 结果与分析 16 4.1 误差对比方法判断重调和方程算法的收敛阶 16 4.1.1 a2?4b≥0的误差对比方法 16 4.1.2 a2?4b0的误差对比方法 18 4.1.3 结论 20 4.2 图像方法判断重调和方程算法的收敛阶 21 4.2.1 a2?4b≥0的图像方法 21 4.2.2 a2?4b0的图像方法 23 4.2.3 结论 24 4.3 特例（大步长，高精度） 25 4.3.1 特殊例子 25 4.3.2 结论 26 4.4 抛物型稳定性条件测试 26 4.4.1 验证稳定性条件 26 4.4.2 结论 28 五 总结 29 六 参考文献 30 引言 1.1 背景 微分方程是构造力学等领域的数学模型的主要。椭圆型抛物型贝努利有限差分、有限元 （1） 其中是Laplace算子，是二维平面上的有限区域，是其光滑边界，a，b非负常数。 数值算列： 1. ，，使问题（1）存在精确解。 2. ，，问题（1）存在精确解u=ysinx+xsiny 分别取不同的a，b值计算，并且包括和两种情形。 ②并将上面的方法推广应用到抛物型方程 数值算列，， ，使问题（11）-（14）存在精确解。 方法介绍 由于区域的不同，会导致算法格式的变化，所以这里为了说明更简单，我们固定区域为矩形区域。其他的一些区域，可由矩形区域变化而来，或者进一步讨论即可得到差分格式，我们这里不再讨论。 2.1 五点差分方法 考虑矩形区域Ω={0xa，0yb}上，二阶线性椭圆型方程 Lu≡=f,其中(x,y) ∈Ω ① 第一边值问题。假设矩形区域Ω网格剖分均匀：h1=a/M，h2=b/N。于是网域包含(M-1)*(N-1)个内点，且均为正则内点。 设（i，j）∈Ω，并且u(x,y)充分光滑，则沿着x和y方向分别用中心差商代替导数有 ② ③ 这里表示u(,)。 而qu和右端项f直接有 qu≡ ④， f= ⑤ 于是方程①在方程的(i,j)点被表示为： 当我们略去截断误差=，就可以得到五点差分格式： ⑥ 其中 由方程⑥我们可以看到，点(i，j)与相邻的四个点（i-1,j）、（i+1,j）、（i,j+1）、（i,j-1）都有关系，于是称之为五点差分格式。 图一 五点差分格式(i,j)关系点图 现在我们再进一步简化问题，令p(x,y)=1,q(x,y)=0时，方程①变为Poisson方程 Lu≡ ⑦ 而我们的五点差分格式⑥也相应地变为 ⑧ 又由于我们题目中的求解域刚好是个正方形，所以我们可以把步长h1=h2=h，用正方形网格剖分区域后，五点差分格式⑧简化为 ⑨ 方程⑨就是我们问题中所运用到的最简化五点差分格式。下面的二阶边值问题和推广到抛物型方程都需要利用到它。 2.2 二阶边值问题分解 重调和方程的二阶边值问题 （1） 其中是Laplace算子，是二维平面上的有限区域，是其光滑边界，a，b非负常数。 虽然直接利用差商逼近导数的方法可以得到