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

计算实习报告 姓名：蔡欣麟学号：2013302736 实习目的 1、熟悉偏微分方程数值解的理论知识; 2、提高Inatlab编程能力。 3、进一步加深对使用matbb解决相关数学问题的理解。 二、 实习内容 运用matlab实现对LaPIaCe方程第一边值问题的求解。 三、 算法公式 考虑如下区域C = ΩU ?Ω上求解POiSSOn方程边值问题 ?u?nLU = —?u = f(xfy) ∈ (a,b) × (Cl ?u ?n =心,y) 的差分方法。 用差分方法求解椭慣I边值问题时，对区域离散化统一采用矩形 网格剖分。对迈是矩形区域的情况，可将G,b］等分为N等分，记 Λ1 =HF∕N;将［c,d］等分为M等分，记h2 = Σd∕M; Q的离散结 果记为Q人，网格节点记为（i, j）,边界结点为i-1和N+l, j=l和 M+1的情形，内结点为2 Si SN, 2≤j≤M的情形。 在差分方程中，微分方程的离散化处理都釆用差商代替微商的方 法。为此，要先假设解函数U （x,y）足够光滑，以便利用TayIor 展开式做出弄商近似。 设（Ij）是内节点，为了表现在该结点处的微商，应用TayIol 展式，有 - 2u(xityi) + U(Xi-I,?)] ?2(xifyi) I h12 ?4u(xifyi) h14 ?6u(xifyi) = ?x2 12 ?x4 360 ?x6