工程电磁场报告讲解.doc

实 工 验 程 报 电 告 磁 场 自动化学院 2011级 指导老师：陈劲操 1110190216韩濛 1110190219花洁 1110190235杨阳 1110190241朱源 1110190243邹文成 超松弛迭代法求解接地金属槽内电位分布 一、实验任务 用超松弛迭代法求解接地金属槽内的电位分布在完成第一个题目的要求后，分析40*40的情况。 正方形金属槽边长为a,以10mm为h取当中九个点和40*40讨论： 已知： 。 给定边值：除顶边电位为100V外，其余三遍为0V 给定边值： 误差范围： 计算：迭代次数 N =? , 分布。 实验原理与方法 利用有限差分法进行计算 其中φ的迭代公式如下所示： 操作过程如下流程图所示： 其中初值的确定采用迭代法以C++程序先予以求出，再将初值带入进行运算，这样得到的结果往往更加合理、准确。但是本身作为并不准确地计算，初值将未知的值设为0也不会有太多误差。 实验过程与结论 通过C++编程实现对于实验要求的模拟，首先模拟实验要求,模拟所需要的程序代码见所附的五份代码，程序运行时的截图如下图所示： 但是有限差分法在所取的点数在较少的时候其取值往往不够准确，因此，为了观察其误差情况，再多取几组进行模拟，将程序参数略作改动，将h取为4mm进行模拟，程序仍见所附代码，运行的截图如下所示。不难看出，当所取的研究点变多时，相应位置的点所对应的电位的值会有较大的变化，由于所取的h越小，所取的点越多，可以由此研究老师所留的第二个实验问题，研究取40*40个点的问题。 对于40*40个点的问题，仍只需将程序略作修改即可，运行的截图如下图所示，程序的代码见所附的第一份程序。由于所讨论的点较多，所以呈现出来不太美观。但应当能比较清晰的得到每一行的41个值，总计40行。 主要利用第二三组数据，可以较为准确的生成下图所示的等点位图: 所得到的图与上课定性分析得到的图是一致的。 总结与体会 对本实验中的电位进行分析是有限差分法的典型应用，由于计算机具有强大的计算能力，所以，通过编程定量的估计计算各点的电位值也因此成为可能。在今后的学习中，我们将继续努力学习专业知识，并学习利用计算机等分析工具，更好的对知识进行把握与运用。感谢陈老师的指导，感谢在本报告产生过程中提供帮助的同学们，谢谢！ 附录 程序代码： 1110190243 邹文成 第一部分： #includeiostream.h #includemath.h #includeiomanip.h int i,j,n; float M[5][5],N[5][5],a; void L(void) { int m; do { for(i=0;i5;i++) for(j=0;j5;j++) N[i][j]=M[i][j]; for(i=1;i4;i++) for(j=1;j4;j++) M[i][j]=N[i][j]+(a/4)*(N[i+1][j]+N[i][j+1]+M[i-1][j]+M[i][j-1]-4*N[i][j]); for(i=1;i4;i++) { for(j=1;j4;j++) if(fabs(M[i][j]-N[i][j])1.0e-5) m=1; else m=0; } n++; }while(m); } void CreatM() { for(i=0;i5;i++) for(j=0;j5;j++) M[i][j]=0; for(j=0;j5;j++) {M[0][j]=100; } } void main() { CreatM(); cout一开始假设初始电位为：\n; cout100setw(6)100setw(6)100setw(6)100setw(6)100setw(6)\n; cout0setw(6)?setw(6)?setw(6)?setw(6)0\n; cout0setw(6)?setw(6)?setw(6)?setw(6)0\n; cout0setw(6)?setw(6)?setw(6)?setw(6)0\n; cout0setw(6)0setw(6)0setw(6)