椭圆型微分方程.docxVIP

纟汐理工只联 数学与计算科学学院 实验报告 实验项目名称 所属课程名称 实验项目名称 所属课程名称 实验类型 实验日期 椭圆型方程数值解 微分方程数值解法 验证 信计0902 一、实验概述： 【实验目的】 1掌握椭圆型方程的五点差分方法 2掌握求解线性方程组的Jacobi迭代法 ［实验原理】 【实验环境】 Microsoft Visual C++6.0 二、实验内容： [实验方案】 用五点差分格式的Jacobi迭代法求解单位正方形区域Q上的Laplace方程Dirichlet边值问 题 务茅Ogg) 1ST} u(儿y) h = /(匕刃=兀2-尸 【实验过程】(实验步骤、记录、数据、分析) 部分主要算法代码： 仃)Dirichlet边值、系数矩阵A以及向量序列K的初步赋值 for(i=0;i=m;i++)//Dirichlet 边界值 { U[i][O]=f(x[i],y[O]); U[i][m]=f(x[i],y[m]); U[0][i]二 f(x[O],y[i]); } for(i=O;in;i++) for(j=0;jn;j++) { a[i][j]=O; KLi]=O; for(i=0;in;i++)// 系数矩阵 for(j=0;jn;j++) { if(i==j) a[i][j]=4; else if(j-i==l j%(m-l )!=0) { a[i]UJ=-l