电磁场研学报告

分离变量法计算与电磁场数值计算方法的应用研究——电磁场与电磁波研究型论文姓名:姜梦磊班级：思源1301学号：13274011指导教师：王国栋摘要本文阐述了静态电磁场边值问题的一种求解方法一有限差分法的基本原理。通过该原理进行MATLAB编程计算电磁场电位问题，通过具体实例利用MATLAB编写程序将这种求解方法的结果可视化处理。以分离变量法的计算结果作为准确的解析解，分析有限差分法的近似计算误差，通过误差分析显示有限差分算法的精确程度。AbstractInthispaper,thebasicprincipleofafinitedifferencemethodforsolvingtheboundaryvalueproblemofstaticelectromagneticfieldispresented..Throughthisprinciple,theproblemofelectromagneticfieldpotentialiscalculatedbyMATLABprogramming.Theconvergencefactorofdifferentaccelerationisselectedtocomparetheiterationnumberofsuccessivesuperrelaxationmethod,andthenselecttheoptimalacceleratingconvergencefactor..TheresultofthismethodisvisualizedbyusingMATLABprogram..Theaccuracyoffinitedifferencemethodisanalyzedbyerroranalysis,andtheaccuracyoffinitedifferencemethodisdemonstratedbyerroranalysis..一、引言通过边值条件求解坐标系内的静态电场分布，是电磁场当中的常见问题之一。求解方法一般为泊松方程或拉普拉斯方程，对于满足下列条件之一的边界面S，则电位函数在该区域中除了任意常数外是唯一确定的：1、在全部S上电位已知；2、在全部S上电位的法向导数已知；3、在一部分S上电位函数已知，在其余S上电位的法向导数已知。利用唯一性定理求解边值问题，分为解析法和数值法两大类，常用的解析法有分离变量法、镜像法、复变函数法、保角变换法、格林函数法等，通常可以得到解析函数表示的闭合解，适于求解各种形状规则的边界下的电位函数。解析法的优点是：1、可用一数学公式将所求问题表示出来，从而得到对该问题的精确描述；2、根据所得的公式，可以对该问题进行深入的分析研究，找到该问题中各参数之间的内在联系和变化趋势；3、可作为近似计算和数值解的检验标准。解析法的缺点是具有一定的局限性，不是每一个电磁问题都能找到与其匹配的解答公式。随着大型快速计算机的广泛应用，近似解析计算和数值法的发展迅速。目前，分析电磁工程常用的数值法有矩量法、有限元法和有限差分法。还有解析和近似计算结合的一些近似解析算法，如逐步逼近算法、微扰法、变分法、几何光学法、物理光学法、几何绕射方法和物理绕射方法等。数值法的主要缺点是不能用一数学公式将所研究的电磁问题和各参数之间的内在关系及变化趋势表示出来。二、算法介绍2.1、解析法利用解析法求解静态场边值问题主要是解决已知场量在场域边界上的值，求场域内场的分布。一般是从MaxWell方程出发，结合具体实例根据边界条件求解出标量函数位函数来表示场的分布。其他场量可以通过位函数再进一步求出。其结果是一个解析表达式，而且往往是一个复杂的级数表达式。不能形象的体现场的分布。但是，得出的结果描述了场域中每一点的值是连续，即：能够精确的表示场域中任何一点的场量。所以也叫精确解。为了形象地表达场的分布，我们将位函数的自变量离散化带入位函数表达式，求出位函数的离散值；通过MATLAB绘图画出等位线。2.2、有限差分法有限差分法（FDM），是一种微分方程数值方法，是通过有限差分来近似导数，从而寻求微分方程的近似解。由泰勒展开式可得有限差分法的基本形式。首先假设要近似函数的各级导数都有良好的性质，依照泰勒定理，可以形成以下的泰勒展开式：其中n!表示是n的阶乘，Rn(x)为余数，表示泰勒多项式和原函数之间的差。可以推导函数f一阶导数的近似值：设定x0=a，可得：除以h可得：求解f'(a)：hxRafhafaf)()()()(1'假设相当小，因此可以将f的一阶导数近似为：局部截尾误差为,其中现在，以静电场边值问题为例，说明有限差分法的应用。f(s)为边界点s的点函数，二位场域D和边界L示于图3.2-1中。图3.2-1有限差分的网格分割通常将场域分成足够小的正方形网格,网格线之间的距离为h（步距）,节点4,3,2,1,0上的电位分别用3210,,,和4表示。设函数在0x处可微，则沿x方向在0x处的泰勒公式展开为Fyx2222)(sfL（3.2.1）（3.2.2）0yx02413DLhhnKnKKK000)(!（3.2.3）将1和3分别代入式（3.2.3）,得03330222001)(!31)(!21)(xhxhxh（3.2.4）03330222003)(!31)(!21)(xhxhxh（3.2.5）由（3.2.4）-（3.2.5）得hxxx2)(310（3.2.6）（3.2.4）+（3.2.5）得2301xx22h2x0)(（3.2.7）同理hyyy2)(310（3.2.8）2301222)(0hyyy（3.2.9）将式(3.2.7)、（3.2.9）代入式（3.2.1），得到泊松方程的五点差分格式)(414243210204321FhFh当场域中,0得到拉普拉斯方程的五点差分格式：)(41044321004321（3.2.10）从这个公式我们可以看出，当我们将一个二维无源区场域剖分为一系列正方形网格时，场域内任何一个节点的电位都等于它周围四个节点电位的算术平均值。这就是规则正方形网格内某点的电位所满足的拉普拉斯方程的差分格式，或差分方程[8]。对于场域内的每一个结点，关系式（3.2.10）式都成立，都可以列出。2.3、有限差分法的求解求解差分方程组的解可使用高斯-赛德尔迭代法及改进的超松弛迭代法。2.3.1、高斯-赛德尔迭代法网格节点一般按“自然顺序”排列，即从左到右，再从下到上顺序排列。迭代也按自然顺序进行。xyo12345678图1网格节点排列方法首先对节点(,)iixy取迭代初值(0),ij。再按下式反复迭代(1)(1)(1)()()1,1,,11,,14,1,2,kkkkkijijijijijij其中(),kij上角标（k）表示k次近似值，下脚标i，j表示节点所在位置，即第i行第j列的交点。（在迭代过程中遇到边界点式，需将边界条件,,ijijf带入。）所有内节点满足以下条件时迭代停止(1)(),,kkijijWW是预定的最大允许误差2.3.1、逐次超松弛法高斯-赛德尔迭代法德变形。为加速收敛，相应的迭代公式为(1)()(1)(1)()()(),,1,,11,,1,44kkkkkkkijijijijijijijα称为“加速收敛因子”，且1≤α2。（逐次超松弛法收敛的快慢与α有明显关系。如何选择最佳α，是个复杂问题。）查阅资料了解到最佳值的计算方法，其值为2opt1+sinh（21）三、问题简述如图所示，有一长方形的导体槽，a=20，b=5，设槽的长度为无限长，槽上有两块与槽绝缘的盖板，电位分别为100V和50V，其它板电位为0V，求槽内的电位分布。yx10050ab四、问题求解4.1、思路此问题选用高斯-赛德尔迭代法，首先列出所有边值条件，选择合理的步长，根据高斯-赛德尔迭代法列出循环表达式，进而求出空间内各点的电位值，画出图像，并通过解析法算出理论精确解，并与数值解进行误差分析。判断数值解法的实用性。4.2、数值法—高斯-赛德尔迭代法4.2.1、仿真程序程序流图启动给定边值,12填写域内的初值迭代次数计数=0N1NN按照超松弛迭代法进行一次迭代，求1,nji所有内点相邻二次迭代值的绝对误差是否小于DF？是输出计算结果：迭代次数N数值解ji,停机调用子程序寻求等值点程序代码%应用高斯-赛德尔迭代法求矩形导体槽内的电位分布clc%清除命令行窗口clear%清除工作区closeallh=input('您所选择的步长为：');hx=20/h+1;col_x=hx-2;%计算设置x方向网格节点数hy=5/h+1;line_x=hy-2;%计算设置y方向网格节点数v1=ones(hy,hx);%设置二维数组，且赋初值v1(hy,:)=ones(1,hx)*100;%给定y=a的边界条件值v1(1,:)=ones(1,hx)*50;%给定y=0的边界条件值fori=1:hy;v1(i,1)=0;%给定x=0的边界条件值v1(i,hx)=0;%给定x=a的边界条件值end%m=10;%w=2/(1+sqrt(1-cos(pi/m)*cos(pi/m)));%计算松弛因子maxt=1;t=0;%设置误差和最大误差参量v2=v1;n=0;while(maxt1e-6)%由v1迭代算出v2，迭代精度为10-6n=n+1;%计算迭代次数maxt=0;fori=2:hy-1;%从2到hy-1行循环forj=2:hx-1;%从2到hx-1列循环%v2(i,j)=v1(i,j)+(v1(i,j+1)+v1(i+1,j)+v2(i-1,j)+v2(i,j-1)-4*v1(...i,j))*w/4;%超松弛差分方程v2(i,j)=(v1(i,j+1)+v1(i+1,j)+v2(i-1,j)+v2(i,j-1))/4;%...高斯-赛德尔差分方程t=abs(v2(i,j)-v1(i,j));%收敛精度判据if(tmaxt)maxt=t;endendendv1=v2;endn%打印迭代次数11v2%打印场域内网格点电位的计算结果表figure(1);subplot(1,2,1),mesh(v2);title(['三维曲面图']);%画三维曲面图subplot(1,2,2),contour(v2,25);title(['等电位线图']);%画等电位线图4.2.2仿真结果4.3、解析法—分离变量法4.3.1、分析过程下面应用静电场中的分离变量法进行解析计算，可以得到问题的精确解yx10050ab设()()XxYy，为满足00x和200x的条件，取()sinxXxAkx其中xk需满足sin0xka故有,20xxmmakmka，则1()sin20mmmXxAx又由220xykk，可得()sinhcosh2020mmYyByCy区域内电位通解可表示为1sinsinh+cosh202020mmmmmmxByCy利用边界条件050y1sin=5020mmmCx对上述方程两边同时乘以sin20nx，并对x从0到20积分，得2020001sinsinsin50202020mmnmnxCxdxxdxn为奇数时可得200nCn同理，应用边界条件5100y可得400sinhcosh44nnnnBCn