精选优质文档-倾情为你奉上计算实习报告姓名:蔡欣麟 学号:一、 实习目的1、熟悉偏微分方程数值解的理论知识;2、提高matlab编程能力。3、进一步加深对使用matlab解决相关数学问题的理解。二、 实习内容运用matlab实现对Laplace方程第一边值问题的求解。三、 算法公式考虑如下区域=上求解Poisson方程边值问题Lu=-u=f x,ya,bc,d= (1)un+k(x,y)u=r(x,y)的差分方法。用差分方法求解椭圆边值问题时,对区域离散化统一采用矩形网格剖分。对是矩形区域的情况,可将a,b等分为N等分,记h1=ab/N;将c,d等分为M等分,记h2=cd/M;的离散结果记为h,网格节点记为(i,j),边界结点为i=1和N+1,j=1和M+1的情形,内结点为2iN,2jM的情形。在差分方程中,微分方程的离散化处理都采用差商代替微商的方法。为此,要先假设解函数u(x,y)足够光滑,以便利用Taylor展开式做出差商近似。设(i,j)是内节点,为了