二维静态电磁场边值问题的数值计算汇编.docVIP

二维静态电磁场边值问题的数值计算 作者：高汉超 (陕西理工学院 物理与电信工程学院 物理1103班 汉中723001) 指导教师：潘峰 [摘要]基于静态电磁场边值的求解问题的实例，分别通过解析法，数值法，有限差分法以及有限元法进行了求解，并得到了电位分布示意图，求出了数值解，并与解析法得到的精确解进行比较，得出了用数值法求解电磁场问题基本满足工程需要的结论。采用数值计算时，对如何减小计算机内存，取消冗余数据，也做了进一步的讨论。最后阐述了三种方法的各自的特征。 [关键词]电磁场;边值;数值解;数值分析;解析法 引言 ???场分布不随着时间的变化而变化的场被称为静态场，静态场的求解对电磁场的分布是至关重要的。事实上知区域内的电荷分，求解静态电磁场就是对已布以及已知区域边界电位及电荷分布进行求解，对已知区域内的场量的分布进行求解，这类问题的求解也被叫做静电场边值问题的求解。静态场的分布求解是基于已知区域边界条件及电荷分布基础上，对满足边界分布的拉普拉斯方程或者泊松方程的求解。通常情况下，求解的方法是数值法和解析法。本文通过对静态电磁场边值的求解，帮助人们加深对电磁场的理解。 1静态电磁场边值求解的解析法 ? 通过解析法对静态电磁场的边值问题进行求解，1864年Maxwell在前人理论和实验的基础上建立了统一的电磁场理论，并用数学模型揭示了自然界一切宏观电磁现象所遵循的普遍规律，这就是Maxwell方程组。笼统而言，所有的宏观电磁问题都可以归结为Maxwell方程组在各种边界条件下的求解问题。解析法包括建立和求解偏微分方程或积分方程。严格求解偏微分方程的经典方法是分离变量法，即在可分离变量的坐标系中求解Maxwell方程组或其退化形式，最后得到解析解。实际上主要针对在场域边界上场量的值已知，而对场域内的场的分布进行求解的情况。求解的结果通常是解析表达式的形式，并且一般较为复杂。解析法虽然不能对体现场的分布进行形象的描述，然而，由于其结果能够对场域内的每一个点的分量都能够进行精确的表述 图1 金属槽 探究一无限长直接地金属槽，其三壁电势为零，顶盖与三壁绝缘且电势为，其中V0=100V，截面长宽分别为a=10cm和b=5cm，如图1所示。求金属槽内的电势分布[]。 分析　金属槽无限长，故槽内电势与坐标z无关。由于槽内各点上电荷密度ρ=0，槽内电势满足二维直角坐标系中的拉普拉斯方程及其边界条件: 应用分离变量法，得到满足方程(1)和边界条件式(2)—式(4)的解的形式为 代入边界条件（5）得 比较系数得： 槽内电势的解析解为 2.数值法 电磁场数值计算是求解电磁场问题重要的方法之一。它将电磁场原本连续的场域问题转换成离散系统，并对其求解数值解。通过在场域离散化的模型上求得各个点上的数值解，近似逼近连续场域的真实解[2]。数值法的出现，使电磁场问题的分析研究从经典方法进入到离散系统的数值分析方法，从而使许多解析法很难解决的复杂的电磁场问题，有可能通过电磁场的计算机辅助分析获得高精度的离散解，同时可极大地促进各种电磁场数值计算方法的发展。有限差分法、有限元法是电磁场数值计算中最常用的两种方法。　上面的例子属于规则形状的第一类边值问题，通过分离变量法已得到精确的解析解，为了与解析解作比较，以验证数值计算的精度，还以此为例对槽内电势进行数值分析。 2.1静态电磁场边值求解的有限差分法有限差分法是将偏微分方程中的偏导函数用差商形式来表示，将所求电磁场的区域中计算无限多个点的函数值变为计算有限多个点上的函数(这一过程称之为离散化)，求出数值解的方法[3]。 2.1.1有限差分法的计算 本文选用MATLAB来编写程序，MATLAB是近年来十分流行的通用性很强的优秀软件，它的程序简单明了，容易看懂[]。而且MATLAB还具有一些更方便的特殊功能，如有专门实现偏微分方程数值求解的工具箱PDEToolbox等，使用这些工具箱能够直观、快速、准确、形象地描述数值计算的结果。 为简单起见，取步长h=1，x、y方向的网格数为m=10，n=5，共有10×5=50个网孔，11×6=66个节点，其中槽内节点(电势代求点)有9×4=36个，边界节点(电势已知点)66-36=30个。采用1/4有限差分形式： 设迭代精度为10-6，利用MATLAB编制的主要计算程序如下: hx=11;hy=6; %设置网络节点 v1=ones(hy,hx); %设置行列二维数组 for j=1:hx %上下两行的 Dinichlet 边界条件 vl(hy,j)=100*sin(pi*(j-1)/(hx-1); v1(1,j)=0; end for i=1:hy %左右两列的 Dirichlet 边界条件