电磁场研学说课.doc

分离变量法计算与电磁场数值计算方法的应用研究 ——电磁场与电磁波研究论文 ：思源：：王国栋 bstract In this paper, the basic principle of a finite difference method for solving the boundary value problem of static electromagnetic field is presented.. Through this principle, the problem of electromagnetic field potential is calculated by MATLAB programming. The convergence factor of different acceleration is selected to compare the iteration number of successive super relaxation method, and then select the optimal accelerating convergence factor.. The result of this method is visualized by using MATLAB program.. The accuracy of finite difference method is analyzed by error analysis, and the accuracy of finite difference method is demonstrated by error analysis.. 一、引言 边值条件求解坐标系的电场分布是电磁场当中的常见问题之一。 泊松方程或拉普拉斯方程对于满足下列条件之一的面则在该区域中除了任意常数外是唯一确定的： 1在全部上电位已知 2、在全部S上电位的法向导数已知 3、在一部分S上，在其余S上电位的法向导数已知。 定理求解问题，分为解析法和数值法两大类，常用的解析法有分离变量法、镜像法、复变函数法、保变换法格林函数法等，通常可以得到解析函数表示的闭合解，适于求解各种形状规则的边界下的电位函数 3、可作为近似计算和数值解的检验标准。解析法的缺点是具有一定的局限性，不是每一个电磁问题都能找到与其匹配的解答公式。 随着大型快速计算机的广泛应用，近似解析计算和数值法的发展迅速。目前，分析电磁工程常用的数值法有矩量法、有限元法和有限差分法。还有解析和近似计算结合的一些近似解析算法，如逐步逼近算法、微扰法、变分法、几何光学法、物理光学法、几何绕射方法和物理绕射方法等。数值法的主要缺点是不能用一数学公式将所研究的电磁问题和各参数之间的内在关系及变化趋势表示出来。 算法介绍 2.1、解析法 利用解析法求解静态场边值问题主要是解决已知场量在场域边界上的值，求场域内场的分布。一般是从MaxWell方程出发，结合具体实例根据边界条件求解出标量函数位函数来表示场的分布。其他场量可以通过位函数再进一步求出。其结果是一个解析表达式，而且往往是一个复杂的级数表达式。不能形象的体现场的分布。但是，得出的结果描述了场域中每一点的值是连续，即：能够精确的表示场域中任何一点的场量。所以也叫精确解。为了形象地表达场的分布，我们将位函数的自变量离散化带入位函数表达式，求出位函数的离散值；通过MATLAB绘图画出等位线。 2.2、有限差分法 由泰勒展开式可得有限差分法的基本形式首先假设要近似函数的各级导数都有良好的性质，依照泰勒定理，可以形成以下的泰勒展开式： 其中n!表示是n的阶乘，Rn(x)为余数，表示泰勒多项式和原函数之间的差。可以推导函数f一阶导数的近似值： 设定x0=a，可得： 除以h可得： 求解f(a)： 假设相当小，因此可以将f的一阶导数近似为： 为 , 其中 为例，说明有限差分法的应用。f(s)为边界点s的点函数，二位场域D和边界L示于图3.2-1中。 图3.2-1 有限差分的网格分割 通常将场域分成足够小的正方形网格,网格线之间的距离为（步距）,节点上的电位分别用和表示。 设函数在处可微，则沿方向在处的泰勒公式展开为 （3.2.3） 将和分别代入式（3.2.3）,得 （3.2.4） （3.2.5） 由（3.2.4）-（3.2.5）得 （3.2.6） （3.2.4）+（3.2.5）得 （3.2.7） 同理