电磁场与电磁波大作业
学院:电子工程学院
班级: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 泊松方程和拉普拉斯方程
拉普拉斯方程是一个二阶偏微分方程,可以用解析法、数值分析法、实验模拟和图解法等求解。
电场的位函数是一个标量函数,简称为电位,电位 EMBED Equation.KSEE3 \* MERGEFORMAT 的定义由下式确定
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.1.1)
电位的单位是伏(V),因此电场强度的单位是伏/米(V/m)。
将式1.1.1带入高斯定理的微分形式 EMBED Equation.KSEE3 \* MERGEFORMAT ,得
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.1.2)
(式1.1.2)称为泊松方程,若讨论的区域 EMBED Equation.KSEE3 \* MERGEFORMAT ,则电位微分方程为
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.1.3)
二阶微分方程(式1.1.3)称为拉普拉斯方程。
其中 EMBED Equation.KSEE3 \* MERGEFORMAT 在直角坐标中为
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.1.4)
1.2 边值问题的分类
静电场的计算通常是求场内任一点的电位,一旦电位确定,电场强度和其他的物理量都可由电位求得。在无界空间,如果已知分布电荷的体密度,可以通过积分公式计算任意点的电位。但计算有限区域的电位时,必须使用所讨论区域边界上电位的指定值(称为边值)来确定积分常数;此外,当场域中有不同介质时,还要用到电位在边界上的边界条件。这些用来决定常数的条件,常统称为边界条件。我们把通过微分方程及相关边界条件描述的问题,称为边值问题。
实际上,边界条件(即边值)除了给定电位在边界上的值以为,也可以是电位在边界上的方向导数。根据不同形式的边界条件,边值问题通常分为三类;
第一类边值问题(Dirichlet问题):给定整个边界上的位函数值,即 EMBED Equation.KSEE3 \* MERGEFORMAT ;
第二类边值问题(Neumann问题):给定整个边界上每一点位函数的法向函数,即 EMBED Equation.KSEE3 \* MERGEFORMAT ;
第三类边值问题(混合问题):给定一部分边界上每一点的电位,同时给定另一部分边界上每一点的电位法向函数,即 EMBED Equation.KSEE3 \* MERGEFORMAT , EMBED Equation.KSEE3 \* MERGEFORMAT 。
1.3 直角坐标系中的分离变量法
分离变量法是数学物理方法中应用最广的一种方法,它要求所给的边界面与一个适当的坐标系的坐标面相重合,或分段重合;其次在此坐标系中,待求偏微分方程的解可表示成三个函数的乘积,每一函数仅是一个坐标的函数。这样通过分离变量法就可以把偏微分方程化为常微分方程求解。
在直角坐标系中,拉普拉斯方程为
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.3.1)
设 EMBED Equation.KSEE3 \* MERGEFORMAT 可以表示为三个函数的乘积,即
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.3.2)
其中,X只是x的函数,同时Y只是y的函数,Z只是z的函数。将(式1.3.2)带入(式1.3.1),得
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.3.3)
然后(式1.3.3)各项除以XYZ,得
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.3.4)
以上方程的第一项只是x的函数,第二项只是y的函数,第三项只是z的函数,要这一方程对任一组(x,y,z)成立,这三项必须分别为常数,即
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.3.5a)
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.3.5b)
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.3.5c)
这样,就将偏微分方程化为三个常微分方程, EMBED Equation.KSEE3 \* MERGEFORMAT 是分离常量,都是待定常数,与边界条件有关。它们可以是实数,也可以是虚数,且由(式1.3.4)应有
EMBED Equation.KSEE3 \* MERGEFORMAT (式1.3.6)
以上三个常微分方程(式1.3.5a)、(式1.3.5b)和(式1.3.5c)解的形式,与边界条件有关(即与常数 EMBED Equation.KSEE3 \* MERGEFORMAT 有关),以(式1.3.5a)为例说明X的形式与的关系。
当 EMBED Equation.KSEE3 \* MERGEFORMAT 时,则
EMBED Equation.KSEE3 \* MERGEFORMAT
当 EMBED Equation.KSEE3 \* MERGEFORM