文档介绍：工程电磁场实验报告 10 电子袁泳华 2010141062 用超松弛迭代法求解接地金属槽内电位的分布 1、题目要求试用超松弛迭代法求解接地金属槽内电位的分布。已知: cm a4?, mm ah10 4/??给定边值如图所示。给定初值: 0 )0(,? ji?误差范围: 510 ???计算迭代次数, ji,?分布。 1. 编程思路定义一个 5*5 的二维数组, 用于存放场内各点数值。再由已知条件,给区域内各点赋予初始值。迭代 a[i][j]=a[i][j]+n/4*(a[i][j+1]+a[i][j-1]+a[i-1][j]+a[i+ 1][j]-4*a[i][j]) 由此公式给出,通过 for 循环语句实现,利用 M 记录迭代次数。利用 while 语句判断实现迭代结束, while(1) 即循环条件始终为真, 一直进行迭代, 当满足精度要求时迭代结束, 退出循环。输出经过迭代后场点分布的近似值及迭代的总次数。流程图: 0???V 100 ?0?? 0?? NYY 2. 源代码#include<> #include<> void main() { double a[5][5]; int M=0,i,j; double n,t,m,p=2,q; cout<<" 试用超松弛迭代法求解接地金属槽内电位的分布。\n 已知: 赋边界已知电位赋场点初始值累计迭代次数 M=0 迭代次数 M+1 利用公式进行迭代判断是否达到精度要求输出每一个点的点位结束 a=4cm,h=a/4=10mm. 边界条件:除盖板为 100V 外,其余电位均为 0V"<<endl; cout<<" 请输入加速收敛因子( 1≤ n< 2) "<<endl; cin>>n; for(i=1;i<=4;i++) for(j=0;j<=4;j++) { a[i][j]=0; a[0][j]=100; } cout<<" 原电位为:"<<endl; for(i=0;i<=4;i++) { for(int j=0;j<=4;j++) cout<<a[i][j]<<'\t'; cout<<endl; } while(1) { for(i=1;i<4;i++) for(j=1;j<4;j++) { m=a[i][j]; a[i][j]=a[i][j]+n/4*(a[i][j+1]+a[i][j-1]+a[i-1][j]+a[i+1][j]-4*a[i][j]); q=fabs(a[i][j]-m