五点菱形格式-拉普拉斯一边值问题-计算实习报告.docVIP

计算实习报告 姓名：蔡欣麟 学号：2013302736 实习目的 1、熟悉偏微分方程数值解的理论知识； 2、提高matlab编程能力。 3、进一步加深对使用matlab解决相关数学问题的理解。 实习内容 运用matlab实现对Laplace方程第一边值问题的求解。 算法公式 考虑如下区域Ω=Ω∪? Lu=-?u=f x,y∈ ?u 的差分方法。 用差分方法求解椭圆边值问题时，对区域离散化统一采用矩形网格剖分。对Ω是矩形区域的情况，可将[a,b]等分为N等分，记h1=ab/N；将[c,d]等分为M等分，记h2=cd/M；Ω的离散结果记为 在差分方程中，微分方程的离散化处理都采用差商代替微商的方法。为此，要先假设解函数u（x,y）足够光滑，以便利用Taylor展开式做出差商近似。 设（i,j）是内节点，为了表现在该结点处的微商，应用Taylor展式，有 1 1 用以上两式中的沿x和y方向的二阶中心差商直接代替方程（1）中的uxx -1 由于式（2）中只出现u在点（i,j）及其4个邻点共五个节点处的值；另外，对一切内节点（2）都适用；称（2）为五点差分格式。 五点差分格式的截断误差为 R 程序设计 计算实例 在区域求解Laplace方程第一边值问题：(要求取x方向与y方向取相同的步长)： (结果：取h＝0.125，利用五点菱形格式，用Gauss-Seidel进行迭代计算得数值结果按