2021有限差分法实验报告参考

1、有限差分法实验报告参考有限差分法实验报告参考 工 程 电 磁 场 实验报告 有限差分法 用超松弛迭代法求解 接地金属槽内电位的分布V100 迭代顺序可按先行后列，或先列后行进行。 ? 迭代过程遇到边界节点时，代入边界值或边界差分格式，直到所有节点电位满足 SKIPIF 1 0 为止。（2）超松弛迭代法 SKIPIF 1 0 (1-15) 式中： SKIPIF 1 0 加速收敛因子 SKIPIF 1 0 可见：迭代收敛的速度与 SKIPIF 1 0 有明显关系程序源代码#includeiostream.h#includemath.h#includeiomanip.hdouble A55;void

2、 main(void)double BJ55;/数组B用于比较电势int s100;/用于储存迭代次数double d100;/用于记录所有的加速因子d0=1.0;int i,j,N=0,M=0,x;for(i=0;i100;i+)di=0.01*i+d0;/加速因子从1.0到2.0之间的20个数！double w10010;int P,Q; for(P=0;P4;P+) for(Q=0;Q5;Q+)APQ=0;for(P=0;P5;P+)A4P=100;cout数组A的所有元素是:endl;for(i=0;i5;i+)for(j=0;j5;j+)coutAijsetw(6);if(5*i+j

3、+1)%5=0)coutn;int pp=0;for(x=0;x100;x+)dofor(i=0;i5;i+)for(j=0;j5;j+)BJij=Aij;for(i=1;i4;i+) for(j=1;j4;j+) Aij=BJij+(dx/4)*(BJi+1j+BJij+1+Ai-1j+Aij-1-4*BJij);/迭代公式 for(i=1;i4;i+) for(j=1;j4;j+) if(fabs(Aij-BJij)1e-5) pp+; N+; while(pp=9);pp=0;for(i=0;i3;i+)wMi+1=A1i+1;for(i=3;i6;i+)wMi+1=A2i-2;for(

4、i=6;i9;i+)wMi+1=A3i-5;sM=N;M+;N=0;int P,Q; for(P=0;P4;P+) for(Q=0;Q5;Q+)APQ=0;for(P=0;P5;P+)A4P=100;int min=s0;int p,q;cout输出所有的加速因子的迭代次数:n; for(q=1;q100;q+) / coutsqsetw(6);/ if(q%12=0)/ coutn; if(mq) min=sq; p=q; coutendl; if(min=s0)p=0;cout最佳加速因子a=;coutdpn;cout迭代次数为：minn;cout最佳收敛因子对应的各个格内点的电位为：n;for( i=1;i10;i+)coutwpit;if(i%3=0)coutn;coutn;程序框图启动启动赋边界节点已知电位值赋予场域内各节点电位初始值累计迭代次数N=0N=N+1按超松弛法进行一次迭代，求 SKIPIF 1 0所有内点相邻二次迭代值的最大误差是否小于停机打印NY 迭代解程序框图 结果分析 迭代收敛的速度与 SKIPIF 1 0 的关系 收敛因子（ SKIPIF 1 0 ）1.01.71.81.831.851.871.902.0迭代次数（ SKIPIF 1 0 ）1312213