电磁场与电磁波大作业
学院：电子工程学院
班级：021231
指导老师：侯建强
组长：
组员：
基于MATLAB的电磁场数值分析
摘要使用计算机进行电磁场数值分析已成为电磁场的工程开发、科研和教学的重要手段。本文介绍了电磁场数值分析的基本理论，并且基于MATLAB PDE工具箱实现了的静态场的边值型问题的求解。实验结果表明，MATLAB使电磁场问题的求解迅速、简单、方便。
关键词：MATLAB 数值分析法边值型问题
Electromagnetic Field Numerical Analysis Based on
MATLAB
Abstract:Using computers to analyze electromagnetic field has been an important method of the development of projects, research and teaching. The essay introduces some basic theories of electromagnetic field numerical analysis. And basing on MATLAB PDE tool, the electromagnetic field boundary value problem has been solved. Furthermore, the results show that it is easier, more prompt and more convenient to figure it out with the software, MATALAB.
Keywords: MATLAB, Electromagnetic Field Numerical Analysis, boundary value problem
目录
0 引言 (2)
1 数值分析法基本原理 (2)
1.1 泊松方程和拉普拉斯方程 (2)
1.2 边值问题的分类 (3)
1.3 直角坐标系中的分离变量法 (3)
2 例题分析 (5)
3 MATLAB实现 (6)
4 结论 (7)
5 结束语 (8)
6 参考文献 (8)
0 引 言
MATLAB 是一种用于数值计算、可视化及编程的高级语言和交互式环境。使用 MATLAB ，可以分析数据，开发算法，创建模型和应用程序。借助其语言、工具和内置数学函数，您可以探求多种方法，比电子表格或传统编程语言（如 C/C++ 或 Java ™）更快地求取结果。MATLAB 应用广泛，其中包括信号处理和通信、图像和视频处理、控制系统、测试和测量、计算金融学及计算生物学等众多应用领域。在各行业和学术机构中，有一百多万工程师和科学家使用 MATLAB 这一技术计算语言。MATLAB 偏微分方程工具箱(PDE toolbox)可以实现对二维问题高速、准确的求解过程，通过使用用户界面或者M 文件，画出所需要的任何区域，输入方程类型和有关系数，可以显示解的图形和数值解。
1 数值分析法基本原理
当电荷或者电流分布已知时，可以通过积分来计算电场或磁场。但实际上我们通常要处理两种类型的静电场问题。一种是已知场源（电荷、电流分布）直接计算空间个点的场强或位函数，这类问题叫做分布型问题。另一种是已知空间某给定区域内的长远分布和该区域边界上的位函数（或其法向导数），求区域内位函数的分布，这类问题叫作边值型问题。求解这些边值型问题空间电场、磁场的分布可以化为求解给定边界条件下的位函数的拉普拉斯方程或泊松方程，即求解边值问题。
1.1 泊松方程和拉普拉斯方程
拉普拉斯方程是一个二阶偏微分方程，可以用解析法、数值分析法、实验模拟和图解法等求解。
电场的位函数是一个标量函数，简称为电位，电位ϕ的定义由下式确定
ϕ-∇=E （式1.1.1）
电位的单位是伏（V ），因此电场强度的单位是伏/米（V/m ）。
将式1.1.1带入高斯定理的微分形式0
ερ=⋅∇E ，得 0
2ερϕϕ-=∇=∇⋅∇ （式1.1.2）
（式1.1.2）称为泊松方程，若讨论的区域0=ρ，则电位微分方程为
20ϕ∇= （式1.1.3）
二阶微分方程（式1.1.3）称为拉普拉斯方程。
其中2
∇在直角坐标中为 22
22222
z y x ∂∂+∂∂+∂∂=∇ （式1.1.4）
1.2 边值问题的分类
静电场的计算通常是求场内任一点的电位，一旦电位确定，电场强度和其他的物理量都可由电位求得。在无界空间，如果已知分布电荷的体密度，可以通过积分公式计算任意点的电位。但计算有限区域的电位时，必须使用所讨论区域边界上电位的指定值（称为边值）来确定积分常数；此外，当场域中有不同介质时，还要用到电位在边界上的边界条件。这些用来决定常数的条件，常统称为边界条件。我们把通过微分方程及相关边界条件描述的问题，称为边值问题。
实际上，边界条件（即边值）除了给定电位在边界上的值以为，也可以是电位在边界上的方向导数。根据不同形式的边界条件，边值问题通常分为三类；
第一类边值问题（Dirichlet 问题）：给定整个边界上的位函数值，即0()r ϕϕ=； 第二类边值问题（Neumann 问题）：给定整个边界上每一点位函数的法向函数，即0n
ϕϕ∂=∂； 第三类边值问题（混合问题）：给定一部分边界上每一点的电位，同时给定另一部分边界上每一点的电位法向函数，即10r ϕ
ϕ=，2
0r n ϕψ∂=∂。
1.3 直角坐标系中的分离变量法
分离变量法是数学物理方法中应用最广的一种方法，它要求所给的边界面与一个适当的坐标系的坐标面相重合，或分段重合；其次在此坐标系中，待求偏微分方程的解可表示成三个函数的乘积，每一函数仅是一个坐标的函数。这样通过分离变量法就可以把偏微分方程化为常微分方程求解。
