基于MATLAB的均匀平面波仿真

1、课程设计说明书课程设计说明书常用软件课程设计常用软件课程设计题目: 基于基于 MATLABMATLAB 的均匀平面波仿真的均匀平面波仿真 院（部）： 力学与光电物理学院 专业班级： 应用物理 学 号： 学生姓名： 指导教师： 2017 年 7 月 2 日安徽理工大学课程设计（论文）任务书安徽理工大学课程设计（论文）任务书 力学与光电物理 学院 基础与应用物理 教研室学号学号学生姓名学生姓名专业（班级）专业（班级）应物 题目题目基于 MATLAB 的均匀平面波仿真设计设计技术技术参数参数1、平面波知识的复习2、MATLAB 程序的编写3、课程设计说明书的书写设计设计要求要求要求完成设计说明书，内

2、容包括：1、设计背景与意义2、设计原理3、设计总结工工作作量量报告页数不少于 12 页，参考文献不少于六个。工作工作计划计划6 月 19 日到 6 月 21 日进行可行性研究与软件选择；熟悉 MATLAB软件及其运行环境6 月 22 日到 6 月 24 日建立关系表，模块设计6 月 25 日到 6 月 26 日结构设计6 月 27-28 日程序源代码编写、调试程序并编写 Word 文档和实现窗口的截图6 月 29-30 日提交课程设计参考参考资料资料1 康颖等. 大学物理M. 长沙：国防科大出版社，19962 同济大学数学研究室著. 高等数学M. 北京：高等教育出版社，19933 赵凯华. 新

3、概念物理教程M. 北京：高等教育出版社，19944 孙玉发. 电磁场与电磁波M.安徽：合肥工业大学，20065 郝培峰，崔建江，潘峰. 计算机仿真技术M. 北京：机械工业出版社，20096 周建兴，岂兴明，矫津毅. MATLAB 从入门到精通M. 北京：人民邮电出版社，20087 李丽芬，云彩霞，郝鹏伟. 利用 MATLAB 图型技术实现电磁波的可视化J. 信息与电脑，2016：4243.指导教师签字指导教师签字 教研室主任签字教研室主任签字2017 年 6 月 30 日安徽理工大学课程设计（论文）成绩评定表安徽理工大学课程设计（论文）成绩评定表目录目录学生姓名：学生姓名： 学号：学号： 专业

4、班级：专业班级： 应用物理 课程设计题目：课程设计题目： 基于 MATLAB 的均匀平面电磁波仿真 指导教师评语：指导教师评语：成绩：成绩： 指导教师：指导教师： 2017 年 7 月 6 日 摘要摘要.51 绪论绪论.11.1 问题背景.11.2 课题研究意义.12 均匀平面电磁波均匀平面电磁波.32.1 定义与性质.32.2 理想介质中的均匀平面波方程.32.3 平面电磁波的瞬时值形式.63 MATLAB 软件及其基本指令软件及其基本指令.83.1 MATLAB 发展历史 .83.2 MATLAB 的功能与语言特点 .83.3 MATLAB 指令 .94 程序设计与运行程序设计与运行.11

5、4.1 设计思路与框图.114.2 运行结果.125 项目总结项目总结.136 参考文献参考文献.14摘要平面波是指场矢量的等相位面与波传播方向相垂直的无限大平面的一种电磁波12。如果平面波在均匀一致且各向同性的理想介质中将形成均匀平面波。均匀平面波是研究电磁波的基础，研究均匀平面波传输特性有十分重要的实际意义。然而直接观察均匀平面波是很难实现的，所以随着计算机的发展，仿真实验正在不断的发展，仿真软件通过图形化界面联系理论条件与实验过程，同时运用一定的编程达到模拟现实的效果。于是本文用 MATLAB 对均匀平面电磁波在理想介质中的传播进行仿真模拟，从而可以更加形象的学习与理解电磁波的知识。关键

6、词：电磁波; 均匀平面电磁波;理想介质; MATLAB; 仿真 01 绪论1.1 问题背景1.1.1MATLAB1.1.1MATLAB 软件简介软件简介MATLAB 是 matrix&laboratory 两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国 mathworks 公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（

7、如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。1.1.21.1.2 均匀平面电磁波特点：均匀平面电磁波特点：1电场强度 E、磁场强度 H 与传播方向三者之间相互垂直，成右手螺旋关系，传播方向上无电磁分量，和为横电磁波，记为 TEM 波；2电场强度 E 和磁场强度 H 处处同相，两者复振幅之比为媒介的阻抗 n，成为波阻抗；3电场、磁场相位相同，等相位面为平面，等相位面垂直于传播方向； 4顺着传播方向，场矢量的时变状态逐渐滞后。1.2 课题研究意义仿真(Simulation)是对现实系统的某一层次抽象属性的模仿。人们利用不同的模型进行实验，从中得到所需的信息，然后帮助人们

8、对现实世界的某一层次问题做出决策。建立适当的模型，在模型上进行动态实验和研究，从中获取相关信息，然后用一系列有目的、有条件的计算机仿真实验来反映系统的特征，得出数量指标，为决策者提供有关这一过程或系统的定量分析结果，作为决策理论依据。仿真实验是利用计算机编制实验进程的行为。采用仿真技术通过数学建模、设计虚拟仪器、虚拟实验环境观察系统模型各量的变化、过程进展、结果评估的全过程。在仿真实验设计阶段，利用计算机进行数学仿真和数值计算，即修1改、变换模型；在设备研制阶段，用已知实际部件或子系统去代替部分计算机仿真模型进行半实物仿真实验，以提高仿真实验的可信度。在系统开发的阶段，通过半实物仿真实验来修改

9、各部件或子系统的结构和参数。对均匀平面电磁波的仿真，可以更加直观的观察电磁波的传播的过程。由于电磁波比较抽象又不可触摸，同时电磁波又是动态的，每时每刻的位置和状态又在改变。通过 MATLAB 来设计对电磁波仿真，可以比较直观的观察电磁波传播的各个时刻的状态。22 均匀平面电磁波2.1 定义与性质均匀平面波，是指电磁波的场矢量只沿着它的传播方向变化，在与波传播方向垂直的无限大平面内，电场强度 E 和磁感应强度 H 的方向、振幅和相位都保持不变。均匀平面波的波速（又称相速）与媒质特性有关。在自由空间中，波速与光速相等为：m/s1031v800电场和磁场的量值之间有简单的关系：yx0HE称为波阻抗或

10、媒质的本质阻抗。0在任意时刻，任意点上的电能密度与磁能密度相等，各占电磁总能量的一半。对于均匀平面波，电磁能量沿着波的传播方向流动，其流动速度即为波速v。2.2 理想介质中的均匀平面波方程均匀平面电磁波满足麦克斯韦方程组为：00DBtBEtDH将和代入得，EDHB300HEtHEtEH）（）（）（）（4-23-22-21-2 方程前两个式子都含有 E 和 H，不便于求解。我们可将这两个方程经过一系列数学变换综合成每个方程只含有一个变换的方程式，对第二个式子取旋度得(2-5)()()(HDttHE将方程第一个式子代入（2-5）得：(2-6)(22HtEE）（上式可变化为(2-7)222tEE同理

11、，有，(2-8)222tHH将带入上两式得，1(2-9)22221tEE(2-10)22221tHH式（2-9）和（2-10）就是理想介质中的电磁场方程。 也是 E 和 H 应满足的波动方程。在平面电磁波中，电磁波沿着与等相位平面垂直的方向传播。假设电磁波沿 x 轴方向传播，则各场最只是空间坐标 x 和时间坐标 t 的函数，所以(2-9)和(2-10)式可化简为4(2-11)222221tExE(2-12)222221tHxH方程涉及矢量问题，给求解带来了麻烦，为得到不含矢量的平面波方程，将方程在直角坐标系中展开得，tEHyyyyxxzxyyzxxyzetEetEetEeyHxHexHzHez

12、HyH)()(）（由于电磁波沿 x 轴传播，则 E 和 H 沿 y 轴和 z 轴方向没有变化，上式分解后得，(2-13)(2-14)tEx0tExHyz(2-15)tExHzy同理，对于方程则有，tHE(2-16)(2-17)tHx0tHxEyz(2-18)tHxEzy由式（2-13）和（2-16）知，Ex和 Hx是与时间无关的常量，可令其值为零。因此对于均匀平面波，E 和 H 都只有与波的传播方向垂直的分量，这种电磁称为横电磁波。将（2-15）和（2-17）对 x 求导并分别将式（2-17）和（2-15）代入，得：222221tHxHyy222221tExEyy5同理，又有：222221tH

13、xHzz2222z21tExEz以上两组方程就是理想介质中的均匀平面波方程。若取 y 轴与 H 方向一致，则只需前一组方程;若取 y 轴与 E 方向一致，则只需后一组方程。2.3 平面电磁波的瞬时值形式假设一个均匀平面电磁波沿 x 轴传播，电场 E 方向与 z 轴平行，由此，来推导平面电磁波的瞬时值形式。由亥姆霍兹方程：022kE由于电场方向与 z 轴平行，可设，)(xEeEzz得0222zzEkxE其解为kzjkzjzeAeA21E由边界条件决定，21, AA111jmeEA 222jmeEA 于是，瞬时值形式为：（2-19）)cos(),(zkxtEetxEmzz与电场 E 相伴的磁场 H

14、 可由求得，HjE6EeHx1其瞬时值形式：)cos(1ykxtEeHmy其中是媒质的本征阻抗。于是，均匀平面波的瞬时值的表达式的一般形式为：)cos(t)(x,H)cos(t)(x,EzzymyzmzkxtEekxtEe73 MATLAB 软件及其基本指令MATLAB 是 MATrix LABoratory 的缩写，是一款由美国 MathWorks 公司出品的商业数学软件。MATLAB 是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。MATLAB 作为高性能、交互式的科学计算工具，具有非常友好的图形界面，这使得 MATLAB 的应用非常广泛；同时MATLA

15、B 也提供了强大的绘图功能。3.1 MATLAB 发展历史70 年代中期，Cleve Moler 和他的同事开发了 LINPACK 和 EISPACK 的Fortran 子程序库；70 年代末期，Cleve Moler 在新墨西哥大学给学生开线性代数，为学生编写了接口程序，这程序取名为 MATLAB，即 MATrix LABoratory；1983 年春天，工程师 John Little 与 Moler、Steve Bangert 一起开发了第二代专业版 MATLAB；1984 年，MathWorks 公司成立，MATLAB 正是推向市场；1992 年，学生版 MATLAB 推出；1993 年

16、，Microsoft Windows 版MATLAB 面世；1995 年，推出 Linux 版。3.2 MATLAB 的功能与语言特点3.2.13.2.1MATLABMATLAB 主要功能主要功能1、基本的数据处理2、优化和解方程3、动态过程仿真：实时的和非实时的4、数据来源：Excel、数据库、A/D 等等5、嵌入式的控制：Pc/104 和 DSP86、神经元网络、小波分析、GA 等等7、虚拟现实仿真3.2.23.2.2MATLABMATLAB 语言特点语言特点MATLAB 语言是一种交互性的数学脚本语言，其语法与 C/C+类似。它支持包括逻辑（boolen）、数值（numeric）、文本（

17、text）、函数柄（function handle）和异质数据容器（heterogeneous container）在内的 15 种数据类型，每一种类型都定义为矩阵或阵列的形式（0 维至任意高维）执行 MATLAB 代码的最简单方式是在 MATLAB 程序的命令窗口（Command Window）的提示符处（ ）输入代码，MATLAB 会即时返回操作结果（如果有的话）。此时, MATLAB 可以看作是一个交互式的数学终端，简单来说，一个功能强大的“计算器”。MATLAB 代码同样可以保存在一个以.m 为后缀名的文本文件中，然后在命令窗口或其它函数中直接调用。MATLAB 语言具有下述显着特点：

18、1、具有强大的矩阵运算能力：Matrix Laboratory（矩阵实验室），使得矩阵运算非常简单。2、是一种演算式语言3、MATLAB 的基本数据单元是既不需要指定维数，也不需要说明数据类型的矩阵（向量和标量为矩阵的特例），而且数学表达式和运算规则与通常的习惯相同。4、MATLAB 语言编程简单，使用方便。3.3 MATLAB 指令3.3.1基本指令1、 clc ：擦去一页命令窗口光标回屏幕左上角2、clear ：从工作空间清除所有变量3、Plot（plot3）：绘制二（三）维图形Plot 函数可以接一些参数，来改变所画图像的属性（颜色，图像元素9等）。下面是一些属性的说明bblue（蓝色)

19、. point（点） - solid（实线）ggreen（绿色） o circle（圆圈） : dotted（点线)rred（红色）x x-mark（叉号） -. dashdot (点画线）ccyan（墨绿色） + plus(加号） - dashed(虚线）mmagenta（紫红色） * star（星号）yyellow（黄色） s square（正方形）kblack（黑色） d diamond（菱形）3.3.23.3.2程序可能用到的的指令程序可能用到的的指令grid on在画图的时候添加网格线。hold on在当前图的轴（坐标系）中画了一幅图，再画另一幅图时，原来的图还在，与新图共存，都看得

20、到axis(xmin xmax ymin ymax) 用来设置 axes 的样式，包括坐标轴范围，可读比例zeros(size(x)生成与 X 相同大小的全零矩阵stem(stem3)绘制二维（三维）函数针状图pause(t)暂停 t 秒后继续执行程序104 程序设计与运行 4.1 设计思路与框图 图 4-1 程序设计框图开始初始化参量设置时间 t 的初始值计算 E(z,t),H(z,t)的值刷新屏幕绘制坐标系绘制图形并着色设置质点形状位置图形保持tT?t=t+1结束114.2 程序编写打开 MATLAB 界面如图 4-2，新建脚本图 4-2 MATLAB 打开界面首先，对瞬时值表达式中的常量

21、进行定义并赋值，代码如下：u0=4*pi*1e-7; %自由空间中的磁导率e0=1e-9/(36*pi); %自由空间中的电介质常数Z0=(u0/e0)0.5; %自由空间中的波阻抗f=1e8; %电磁波的频率w=2*pi*f; k=w*(u0*e0)0.5; %波数phi_E=0; %初始相位设为 0phi_H=0;EE=20; %电场振幅HH=EE/Z0; %磁场振幅代码界面如图 4-312图 4-3 常量赋值界面利用电磁波瞬时值表达式画出电磁波图，并且写一个循环，使电磁波传播能狗像动画一样动起来。代码如下：for t=0:1:300 Ez=EE*cos(k*x-w*t*1e-9+phi_

22、E); Hy=HH*cos(k*x-w*t*1e-9+phi_H); plot3(x,m0,Ez,b,LineWidth,2); hold on;grid on; plot3(x,Hy,m0,r,LineWidth,2); hold off xlabel(传播方向) ylabel(磁场 Hy) zlabel(电场 Ez) title(平面电磁波传播示意图,t=,num2str(t),ns,fontsize,14) drawnow pause(0.01)end代码截图界面如图 4-4.134-4 完整代码4.3 运行结果运行结果如下图4-2 程序运行结果14如图，蓝色的是电场传播的波，红色的是磁

23、场传播的波。由图可得的结论有：1、电场和磁场在空间相互垂直与传播方向，E、H、x 轴满足右手螺旋关系（TEM 波）。2、电场、磁场相位变化相同3、电场、磁场的振幅不随传播距离增加而改变155 项目总结这次课程设计是对我们学习常用软件课程结果的一次大检验。通过这次基于 MATLAB 的均匀平面波仿真课程设计让我们更熟悉地掌握该软件的功能，灵活运用 MATLAB 软件，加强对 MATLAB 软件强大的图形处理功能，掌握利用MATLAB 绘图功能做出相关函数曲线，从而运用 MATLAB 分析。在熟悉掌握编写 MATLAB 程序和操作的同时培养了我们的独立思考能力，钻研精神，解决问题能力和动手能力。本次课程设计中通过查阅资料，阅读程序并读写程序对 MATLAB 均匀平面电磁波的仿真有了更深的了解。我同时也认识到了 MATLAB 功能并不只是图形的绘制及波形的处理，有着很多方面的运用，如绘制函数，处理音频、图像数据，创建用户界面等功能，实为一个功能强大的软件。每一次课给我一种新的体验与感受，这次课程设计一样，平时都是啃课本，很单调很乏味，但是，这次课设给了我们动手和自主学习的机会，当然，在这过程中不免会遇到困难，当时也会很困惑、沮丧，但问题总要解决的，所以自己还是会去想办法，我想这也是一种收获。还有，经过了课程设计，我也深刻体会到 MATLAB 功能的强大以及学好 MATLAB