有限差分法实验报告.docVIP

工 程 电 磁 场 实验报告 ——有限差分法 老师： XXX 姓名： XXX 学号： XXXXXXX 时间： 2011年3月27日 用超松弛迭代法求解 接地金属槽内电位的分布 一、实验要求 按对称场差分格式求解电位的分布 已知： 给定边值：如图1-7示 图1-7接地金属槽内半场域的网格 给定初值 误范围差: 计算：迭代次数，，将计算结果保存到文件中 实验思想 有限差分法 有限差分法（Finite Differential Method）是基于差分原理的一种数值计算法。其基本思想：将场域离散为许多小网格，应用差分原理，将求解连续函数的泊松方程的问题转换为求解网格节点上的差分方程组的问题。 泊松方程的五点差分格式 当场域中得到拉普拉斯方程的五点差分格式 差分方程组的求解方法 (1) 高斯——赛德尔迭代法 (1-14) 式中： ( 迭代顺序可按先行后列，或先列后行进行。 ( 迭代过程遇到边界节点时，代入边界值或边界差分 格式，直到所有节点电位满足为止。 （2）超松弛迭代法 (1-15) 式中：——加速收敛因子 可见：迭代收敛的速度与有明显关系 程序源代码 #includeiostream.h #includemath.h #includeiomanip.h double A[5][5]; void main(void) { double BJ[5][5];//数组B用于比较电势 int s[100];//用于储存迭代次数 double d[100];//用于记录所有的加速因子 d[0]=1.0; int i,j,N=0,M=0,x; for(i=0;i100;i++) d[i]=0.01*i+d[0];//加速因子从1.0到2.0之间的20个数！ double w[100][10]; int P,Q; for(P=0;P4;P++) for(Q=0;Q5;Q++) A[P][Q]=0; for(P=0;P5;P++) A[4][P]=100; cout数组A的所有元素是:endl; for(i=0;i5;i++) for(j=0;j5;j++) { coutA[i][j]setw(6); if((5*i+j+1)%5==0) cout\n; } int pp=0; for(x=0;x100;x++) { do { for(i=0;i5;i++) for(j=0;j5;j++) BJ[i][j]=A[i][j]; for(i=1;i4;i++) for(j=1;j4;j++) A[i][j]=BJ[i][j]+(d[x]/4)*(BJ[i+1][j]+BJ[i][j+1]+A[i-1][j]+A[i][j-1]-4*BJ[i][j]);//迭代公式 for(i=1;i4;i++) { for(j=1;j4;j++) if(fabs(A[i][j]-BJ[i][j])1e-5) pp++; } N++; }while(pp=9); pp=0; for(i=0;i3;i++) w[M][i+1]=A[1][i+1]; for(i=3;i6;i++) w[M][i+1]=A[2][i-2]; for(i=6;i9;i++) w[M][i+1]=A[3][i-5]; s[M]=N; M++; N=0; int P,Q; for(P=0;P4;P++) for(Q=0;Q5;Q++) A[P][Q]=0; for(P=0;P5;P++) A[4][P]=100; } int min=s[0]; int p,q; cout输出所有的加速因子的迭代次数:\n; for(q=1;q100;q++) { // couts[q]setw(6); // if(q%12==0) // cout\n; if(mins[q]) { min=s[q]; p=q; } } coutendl; if(min==s[0]) p=0; cout最佳加速因子a=; coutd[p]\n; cout迭代次数为：min\n; cout最佳收敛