Matlab在电磁场可视化教学和数值计算方面的应用研究..docx

MATLAB课程设计报告题 目Matlab在电磁场与电磁波数值计算和可视化教学方面应用初探学 院电子信息工程学院专 业电子信息工程姓名和学号刘伟 导教师程鸿Matlab在电磁场与电磁波数值计算和可视化教学方面应用初探刘伟电子信息工程 要：本文利用Matlab实现了SOR超松弛法，绘制了单个点电荷的平面电场线与等势线，和等量同号点电荷的电场线和等势线，并利用GUI对平面简谐波进行了仿真，展现了Matlab在电磁场与电磁波数值计算和可视化教学方面应用前景。关键词：Matlab；电磁场与电磁波；数值计算；可视化；GUI；仿真1 引言电磁场与电磁波是电子信息类专业的一门专业核心课。Maxwell用简洁的方程描绘了描述了电磁场所遵循的普遍规律，但其具体问题的求解往往是复杂的，并且很多时候往往不能直接精确求解，需要进行大量的近似计算来获取结果。除此之外，对于数学表达式形式的求解结果，人们又往往难以形象地理解电磁场具体的空间分布。[1]而Matlab由于其强大的功能、简单易学的编程语言和可视化的仿真环境，则为电磁场与电磁波的计算和可视化提供了一定条件。2 Matlab简介Matlab是美国Mathworks公司出品的商业数学软件，它将数值分析、矩阵计算、科学数据可视化、非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。3 Matlab在电磁场数值计算方面的应用实例Matlab在电磁场数值计算方法上的应用主要包含了有限差分法、矩量法、有限元法、时域有限差分法等，通过利用软件计算相同的问题，可以验证结果并加深对算法的理解[2]我们以有限差分法作为例子。实际问题中，一般借助于数值法求电磁场的数值解[3]。然而在求解时，由于节点个数很多，联立差分方程的个数往往可达几百甚至几千个[4]，因此我们需要采用逐次近似的迭代方法解决问题。最常用的迭代法包括Gauss-Seidel迭代法和SOR迭代法，由于Gauss-Sidel迭代法是SOR迭代法的一种特殊形式，我们以后者作为代表。（1-1）其中参数称作松弛因子，显然，当时，解式（1-1）的SOR方法就是Gauss-Seidel迭代法。[5]实现SOR迭代法的自定义函数sor在附件中。4 Matlab在电磁场和电磁波教学可视化实例4.1Matlab与绘图有关的函数4.1.1二维图形绘图语句(1) plot语句：二维图形绘图命令中最基本的指令就是plot。如果输入两个矢量 x、y则plot(x ,y)产生的是y相对于x的图形。(2) fplot语句：当已知函数关系时，采用fplot语句可以更加快速、精确地绘制出指定区间上的函数图像。4.1.2三维曲线绘图语句(1)plot3(x,y,z,s)语句：X、Y、Z是同维数的矢量或矩阵。每组X、Y、Z构成一个点的坐标，各点依次相连，形成一条曲线。如果是矩阵，则它们相应的列构成一条三维曲线的数据点坐标，所以用矩阵可以同时画多条空间曲线。(2)网线和表面图：Matlab中，曲面是用xy平面的各个格点上的z坐标来定义，相邻点用直线连接。因此，为了描绘二元函数形成的空间曲面，首先要建立平面的数据网格。即根据x、y的定义域[],，[]生成平面上的坐标格点, 每个格点对应一组[ x，y] 。生成数据网格的命令是meshgrid，格式为[] =meshgrid(x，y), 使用 meshgrid 命令，将两个独立的矢量 x、y 组成数据网格。mesh是三维网格作图命令，mesh(x,y,z)画出了每一个格点(x, y)上对应的z值。surfl(X ,Y，Z)是画带光照模式的三维曲面图 。 该命令显示一个带阴影的曲面, 融合了周围的、散射的和镜面反射的光照模式。要获得较平滑的颜色过度，应该使用有线性强度变化的色图，如 gray、copper、bone、pink等。[6] 4.2单个点电荷的平面电场线与等势线点电荷是带电体的一种理想模型。如果在研究的问题中，带电体的形状、大小以及电荷分布可以忽略不计，即可将它看作是一个几何点，则这样的带电体就是点电荷。一个实际的带电体能否看作点电荷，不仅和带电体本身有关，还取决于问题的性质和精度的要求。与质点、刚体等概念一样，点电荷是实际带电体的抽象和近似，它是建立具有普遍意义的基本规律的不可或缺的理想模型，又是把复杂多样的实际问题转化或分解为基本问题时必不可少的分析手段。设有一个带有+q电荷的点电荷，则该点电荷在P(x，y)处产生的电势为： （2-1）由电场强度可得,在xO