抛物型微分方程.docVIP

数学与计算科学学院 实 验 报 告 实验项目名称 抛物型方程数值解 所属课程名称 微分方程数值解法 实 验 类 型 验证 实 验 日 期 班 级 信计0902 学 号 姓 名 成 绩 一、实验概述： 【实验目的】 1熟练掌握向前差分格式和Crank-Nicolson格式 2求解系数矩阵为三对角阵的方程组的追赶法 【实验原理】 【实验环境】 Microsoft Visual C++6.0 二、实验内容： 【实验方案】 用向前差分格式和Crank-Nicolson格式求解一维热传导方程初边值问题 0x1, 0tT u(x,0)=sin(pi *x) 0x1 u(0,t)=u(1,t)=0 0tT 【实验过程】（实验步骤、记录、数据、分析） 1，部分算法代码： (1)三角分解 void fenjie(double a,double b,double c,double e[m-1],double d[m-1])//三角分解 { int i; d[0]=a; for(i=2;i=m-1;i++) { e[i-1]=b/d[i-2]; d[i-1]=a-c*e[i-1]; } } void huidai(double g[m-1],double y[m-1],double e[m-1],double d[m-1],double c,double x[m-1])//回代求解 { int i; for(i=2,y[0]=g[0];i=m-1;i++) y[i-1]=g[i-1]-e[i-1]*y[i-2]; //追的过程 for(i=m-2,x[m-2]=y[m-2]/d[m-2];i=1;i--) { x[i-1]=(y[i-1]-x[i]*c)/d[i-1];//赶的过程 } } (2)Crank_Nicolson法 void Crank_Nicolson(double r,double U[m+1][n+1]) { double a,b,c,g[m-1],e[m-1],d[m-1],y[m-1]; double xx[m-1]; int i,j; a=1+r; b=-1.0/2*r; c=b; for(j=0;j=n;j++) { for(i=0;im-1;i++) { if(i==0) g[0]=(1-r)*U[1][j]-b*U[2][j]-b*(U[0][j]+U[0][j+1]); if(i==m-2) g[m-2]=(1-r)*U[m-1][j]-b*U[m-2][j]-b*(U[m][j]+U[m][j+1]); else g[i]=(1-r)*U[i+1][j]-b*(U[i+2][j]+U[i][j]); } fenjie(a,b,c,e,d); huidai(g,y,e,d,c,xx); for(i=0;im-1;i++) U[i+1][j+1]=xx[i]; } } (3)目标函数及解析解 double u(double t,double x)//微分方程的精确解 { double uu; uu=exp(pow(pi,2)*(-t))*sin(pi*x); return uu; } 2，结果分析： 根据上面的实验原理，可以得出：在每一时间层，需要求解的隐式差分方程形成了一个线性代数方程组，它的系数矩阵是三对角形矩阵，即仅在主对角线及其相邻二条对角线上有非零元素，故可以使用三角分解法。使用Crank_Nicolson隐式法，每时间层包含较多的计算工作量，但是其优点为其稳定性要求对步长比的限制大为放宽，而这正是我们所期望的。 【实验结论（结果） 【实验小结 通过本次实验，我重新复习了三角分解法，也对课本上的Crank_Nicolson方法有了更深的了解，同时也提高了自己的上机编程能力，感受到了理论与实践相结合的乐趣。 三、指导教师评语及成绩： 评 语 评语等级 优 良 中 及格 不及格 1.实验报告按时完成,字迹清楚,文字叙述流畅,逻辑性强 2.实验方案设计合理 3.实验过程（实验步骤详细,记录完整,