微分方程数值解上机报告第二章报告.doc

第二章实验报告 实验六 一、实验目的 通过用Jacobi迭代、Guass-Seidel迭代和SOR迭代等算法求解Dirichlet问题 的数值解，深刻领会微分方程数值解的实质，熟练掌握各算法的计算机实现过程，并能从理论及实验结果分析各种算法的优缺点。 二、实验内容 利用五点差分格式近似Dirichlet问题 取步长1，试用Jacobi迭代、Guass-Seidel迭代和SOR迭代求解。 三、程序 // test.cpp : Defines the entry point for the console application. // #include"cmath" #include"iostream" #include"iomanip" const int N=7; const int M=9; const double pi=3.1415926; double UpFy(double x,double y); double DoFy(double x,double y); double LeFx(double x,double y); double RiFx(double x,double y); using namespace std; int main() { static double Matrix_1[N][M],Matrix_2[N][M]; double x1,x2,y1,y2;//区域边界 double dx,dy;//步长 double Max,e=1.0;//判断标准(误差) int i,j; cout<<"请输入矩形区域区间:\n"; cout<<"<x1,x2>: "; cin>>x1>>x2; cout<<"<y1,y2>: "; cin>>y1>>y2; dx=(x2-x1)/(M-1); dy=(y2-y1)/(N-1); for(i=0;i<N;i++) { Matrix_1[i][0]=LeFx(x1,i*dy+y1); Matrix_1[i][M-1]=RiFx(x2,i*dy+y1); } for(j=0;j<M;j++) { Matrix_1[0][j]=DoFy(x1+dx*j,y1); Matrix_1[N-1][j]=UpFy(x1+dx*j,y2); } while(e>0.00001) { for(i=1;i<N-1;i++) for(j=1;j<M-1;j++) { Matrix_2[i][j]=(Matrix_1[i-1][j]+Matrix_1[i+1][j]+ Matrix_1[i][j+1]+Matrix_1[i][j-1])/4.0; } Max=fabs(Matrix_1[1][1]-Matrix_2[1][1]);//求出误差项最大的元素作为判断标准 for(i=1;i<N-1;i++) { for(j=1;j<M-1;j++) if(fabs(Matrix_2[i][j]-Matrix_1[i][j])>Max) Max=fabs(Matrix_2[i][j]-Matrix_1[i][j]); } e=Max; for(i=1;i<N-1;i++)//为下一次循环做准备 for(j=1;j<M-1;j++) Matrix_1[i][j]=Matrix_2[i][j]; } for(i=0;i<N;i++) { for(j=0;j<M;j++) cout<<setw(8)<<setiosflags(ios::left)<<Matrix_1[i][j]<<" "; cout<<endl; } return 0; } double UpFy(double x,double y) { return 0; } double DoFy(double x,double y) { if(((int)(x)%4==0)&&(x>1))//当x为4的整数倍时为了提高精度，令x=0 x=0; return (sin(x/4.0*pi)); } double LeFx(double x,double y) { return (y*(3-y)); } double RiFx(double x,double y) { return 0; } 四、数值结果及分析 请输入矩形区域区间: <x1,x2>: 0 4 <