基于MATLAB的均匀平面波仿真.docx

1、基于MATLAB的均匀平面波仿真课程设计说明书常用软件课程设计题目: 基于MTLA得均匀平面波仿真 院（部）： 力学与光电物理学院 专业班级: 应用物理 学 号： 学生姓名： 指导教师： 01年7月2日安徽理工大学课程设计（论文)任务书 力学与光电物理 学院 基础与应用物理 教研室学号学生姓名专业(班级）应物题目基于AAB得均匀平面波仿真设计技术参数1、平面波知识得复习2、ATLAB程序得编写3、课程设计说明书得书写设计要求要求完成设计说明书，内容包括：1、设计背景与意义、设计原理3、设计总结工作量报告页数不少于12页,参考文献不少于六个。工作计划6月19日到6月21日进行可行性研究与软件选择

2、;熟悉MTLA软件及其运行环境月2日到6月24日建立关系表,模块设计6月25日到6月26日结构设计6月2728日程序源代码编写、调试程序并编写Word文档与实现窗口得截图6月2-30日提交课程设计参考资料1 康颖等、 大学物理、 长沙:国防科大出版社，19962 同济大学数学研究室著、 高等数学M、 北京:高等教育出版社，1993 赵凯华、 新概念物理教程、 北京：高等教育出版社,99 孙玉发、 电磁场与电磁波M、安徽：合肥工业大学,2065 郝培峰，崔建江,潘峰、 计算机仿真技术、北京:机械工业出版社，20096 周建兴，岂兴明,矫津毅、 MAL从入门到精通M、 北京：人民邮电出版社,008

3、 李丽芬，云彩霞，郝鹏伟、 利用MATAB图型技术实现电磁波得可视化J、 信息与电脑，201:4243、指导教师签字教研室主任签字207年6月3日安徽理工大学课程设计（论文)成绩评定表学生姓名： 学号： 专业班级：应用物理 课程设计题目: 基于MATLAB得均匀平面电磁波仿真 指导教师评语：成绩： 指导教师： 2017年7月6日 摘要51 绪论1、1 问题背景11、 课题研究意义 1均匀平面电磁波32、 定义与性质32、2 理想介质中得均匀平面波方程32、3 平面电磁波得瞬时值形式 6 MAA软件及其基本指令 83、2 TLA得功能与语言特点83、3 MATAB指令 94 程序设计与运行14、

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

5、波；理想介质; MAB; 仿真 绪论1、 问题背景1、ALAB软件简介MATAB就是matrilartory两个词得组合，意为矩阵工厂（矩阵实验室).就是由美国mthworks公司发布得主要面对科学计算、可视化以及交互式程序设计得高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统得建模与仿真等诸多强大功能集成在一个易于使用得视窗环境中,为科学研究、工程设计以及必须进行有效数值计算得众多科学领域提供了一种全面得解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fotran)得编辑模式,代表了当今国际科学计算软件得先进水平。、1、均匀平面电磁波特点：1电场强度E、

6、磁场强度H与传播方向三者之间相互垂直,成右手螺旋关系，传播方向上无电磁分量,与为横电磁波，记为TEM波；2.电场强度E与磁场强度H处处同相,两者复振幅之比为媒介得阻抗n，成为波阻抗；3.电场、磁场相位相同,等相位面为平面,等相位面垂直于传播方向；4。顺着传播方向，场矢量得时变状态逐渐滞后。1、2 课题研究意义仿真（Simulation）就是对现实系统得某一层次抽象属性得模仿。人们利用不同得模型进行实验，从中得到所需得信息，然后帮助人们对现实世界得某一层次问题做出决策.建立适当得模型,在模型上进行动态实验与研究,从中获取相关信息，然后用一系列有目得、有条件得计算机仿真实验来反映系统得特征,得出数

7、量指标，为决策者提供有关这一过程或系统得定量分析结果,作为决策理论依据。仿真实验就是利用计算机编制实验进程得行为.采用仿真技术通过数学建模、设计虚拟仪器、虚拟实验环境观察系统模型各量得变化、过程进展、结果评估得全过程。在仿真实验设计阶段，利用计算机进行数学仿真与数值计算,即修改、变换模型;在设备研制阶段，用已知实际部件或子系统去代替部分计算机仿真模型进行半实物仿真实验,以提高仿真实验得可信度.在系统开发得阶段，通过半实物仿真实验来修改各部件或子系统得结构与参数。对均匀平面电磁波得仿真，可以更加直观得观察电磁波得传播得过程。由于电磁波比较抽象又不可触摸，同时电磁波又就是动态得,每时每刻得位置与状

8、态又在改变。通过MATLB来设计对电磁波仿真，可以比较直观得观察电磁波传播得各个时刻得状态. 均匀平面电磁波2、1 定义与性质均匀平面波，就是指电磁波得场矢量只沿着它得传播方向变化,在与波传播方向垂直得无限大平面内，电场强度与磁感应强度H得方向、振幅与相位都保持不变。均匀平面波得波速（又称相速）与媒质特性有关.在自由空间中,波速与光速相等为:电场与磁场得量值之间有简单得关系:称为波阻抗或媒质得本质阻抗.在任意时刻,任意点上得电能密度与磁能密度相等，各占电磁总能量得一半。对于均匀平面波,电磁能量沿着波得传播方向流动，其流动速度即为波速。、 理想介质中得均匀平面波方程均匀平面电磁波满足麦克斯韦方程

9、组为: 将与代入得， 方程前两个式子都含有E与H，不便于求解。我们可将这两个方程经过一系列数学变换综合成每个方程只含有一个变换得方程式,对第二个式子取旋度得 （25）将方程第一个式子代入（25）得: (26）上式可变化为 （2）同理,有, （28）将带入上两式得， (2) (2-10）式（2)与（210）就就是理想介质中得电磁场方程。也就是E与H应满足得波动方程。在平面电磁波中，电磁波沿着与等相位平面垂直得方向传播。假设电磁波沿x轴方向传播,则各场最只就是空间坐标与时间坐标t得函数，所以（2）与(10)式可化简为 (2-1） (22）方程涉及矢量问题，给求解带来了麻烦,为得到不含矢量得平面波方

10、程,将方程在直角坐标系中展开得， 由于电磁波沿x轴传播，则与H沿y轴与z轴方向没有变化，上式分解后得, （1） (2-14) （1）同理，对于方程 则有， （2-16） （27） （18）由式（23）与(216)知，x与x就是与时间无关得常量，可令其值为零。因此对于均匀平面波,E与H都只有与波得传播方向垂直得分量,这种电磁称为横电磁波.将(2-15）与（21)对x求导并分别将式(2-17)与(2-15）代入，得： 同理,又有： 以上两组方程就就是理想介质中得均匀平面波方程。若取轴与H方向一致，则只需前一组方程；若取轴与方向一致,则只需后一组方程。2、 平面电磁波得瞬时值形式假设一个均匀平面电磁

11、波沿轴传播,电场E方向与z轴平行,由此,来推导平面电磁波得瞬时值形式.由亥姆霍兹方程： 由于电场方向与轴平行，可设，得 其解为由边界条件决定， 于就是,瞬时值形式为： (19）与电场相伴得磁场H可由求得，其瞬时值形式:其中就是媒质得本征阻抗。于就是，均匀平面波得瞬时值得表达式得一般形式为： ALAB软件及其基本指令MATAB就是Ari LABoatr得缩写，就是一款由美国athWrks公司出品得商业数学软件。MATL就是一种用于算法开发、数据可视化、数据分析以及数值计算得高级技术计算语言与交互式环境.MATLAB作为高性能、交互式得科学计算工具,具有非常友好得图形界面，这使得MATAB得应用非

12、常广泛;同时MATLAB也提供了强大得绘图功能.3、1 MALA发展历史70年代中期,Cve oler与她得同事开发了INPACK与EISPCK得otr子程序库； 70年代末期,Clve Mole在新墨西哥大学给学生开线性代数，为学生编写了接口程序,这程序取名为MALAB，即MAr ABoraory; 183年春天,工程师Johnittle与Mlr、Seve Bangert一起开发了第二代专业版MATLAB； 14年，Mathoks公司成立,MATLAB正就是推向市场; 199年，学生版MATLA推出;99年，Mcroft Wndows版MATLA面世；195年,推出Linux版。3、2 MA

13、TLA得功能与语言特点3、2、ATLAB主要功能1、基本得数据处理2、优化与解方程3、动态过程仿真：实时得与非实时得4、数据来源:Excel、数据库、A/D等等5、嵌入式得控制：c104与DS6、神经元网络、小波分析、GA等等7、虚拟现实仿真3、2、2 ATLAB语言特点MTLAB语言就是一种交互性得数学脚本语言，其语法与CC+类似。它支持包括逻辑（boolen）、数值（nueric）、文本(ex)、函数柄（fuction ade）与异质数据容器(hetrogeous cntner）在内得种数据类型,每一种类型都定义为矩阵或阵列得形式(维至任意高维） 执行MATLAB代码得最简单方式就是在MA

14、TAB程序得命令窗口(and Widow）得提示符处()输入代码,MATLAB会即时返回操作结果(如果有得话)。此时,AAB可以瞧作就是一个交互式得数学终端，简单来说，一个功能强大得“计算器”。MATLAB代码同样可以保存在一个以、m为后缀名得文本文件中,然后在命令窗口或其它函数中直接调用。 MATLAB语言具有下述显着特点：1、具有强大得矩阵运算能力：MtriLaoraoy（矩阵实验室），使得矩阵运算非常简单.2、就是一种演算式语言3、MA得基本数据单元就是既不需要指定维数，也不需要说明数据类型得矩阵(向量与标量为矩阵得特例），而且数学表达式与运算规则与通常得习惯相同。 4、MATLAB语言

15、编程简单,使用方便.3、3 MAAB指令3、3、1 基本指令、clc ：擦去一页命令窗口光标回屏幕左上角、cer ：从工作空间清除所有变量3、lot(t3)：绘制二(三)维图形Po 函数可以接一些参数，来改变所画图像得属性(颜色,图像元素等）。下面就是一些属性得说明b lue(蓝色) 、point(点） ld（实线)g reen（绿色） oircl（圆圈) ： oted（点线）r re(红色） x x-mark(叉号）-、 dshdot(点画线） cya(墨绿色）+pls(加号）- dashd(虚线)mmaen(紫红色） * tr（星号）yyellow(黄色) ssqure(正方形)k lak

16、（黑色） diand(菱形)3、3、2程序可能用到得得指令gridon在画图得时候添加网格线.hd on在当前图得轴(坐标系)中画了一幅图，再画另一幅图时，原来得图还在，与新图共存，都瞧得到axis(in a yin yma） 用来设置aes得样式,包括坐标轴范围，可读比例zeros(ize(x)）生成与X相同大小得全零矩阵st(ste3） 绘制二维（三维）函数针状图paue(t）暂停t秒后继续执行程序4 程序设计与运行 、1 设计思路与框图 图1程序设计框图4、 程序编写打开MTAB界面如图4-2，新建脚本图4 MATLAB打开界面首先，对瞬时值表达式中得常量进行定义并赋值,代码如下:u0=

17、4pi1e7; %自由空间中得磁导率e01e-/(36); %自由空间中得电介质常数0=(u0/0)、； 自由空间中得波阻抗f=1e8; 电磁波得频率w2*pif； =w*(u00)0、5; %波数hiE=0； 初始相位设为0phi_=0；E20； 电场振幅HE/Z0； 磁场振幅代码界面如图43图43 常量赋值界面利用电磁波瞬时值表达式画出电磁波图,并且写一个循环，使电磁波传播能狗像动画一样动起来。代码如下：ort=：1:30 z=Ec(k*te9phiE）； Hy=HH*c(kx*t*1e+phi_H）； plot3（x，m0，Ez,LneWith，2)； olon;ri ; plot（x，

18、Hy，0，r,LneWidth，2）； old off xel（传播方向) labl（磁场Hy） abel(电场） title(平面电磁波传播示意图，t=，um2sr（t）,fotize,14） dawnow paue(、1）nd代码截图界面如图44、4完整代码4、3 运行结果运行结果如下图4-2程序运行结果如图，蓝色得就是电场传播得波，红色得就是磁场传播得波.由图可得得结论有:1、电场与磁场在空间相互垂直与传播方向,E、H、轴满足右手螺旋关系（TM波）.2、电场、磁场相位变化相同3、电场、磁场得振幅不随传播距离增加而改变5项目总结这次课程设计就是对我们学习常用软件课程结果得一次大检验。通过这

19、次基于MTA得均匀平面波仿真课程设计让我们更熟悉地掌握该软件得功能，灵活运用MATLAB软件,加强对ATLB软件强大得图形处理功能，掌握利用MATL绘图功能做出相关函数曲线,从而运用MATLAB分析。在熟悉掌握编写MLA程序与操作得同时培养了我们得独立思考能力，钻研精神，解决问题能力与动手能力.本次课程设计中通过查阅资料，阅读程序并读写程序对MATA均匀平面电磁波得仿真有了更深得了解。我同时也认识到了MATLAB功能并不只就是图形得绘制及波形得处理,有着很多方面得运用，如绘制函数，处理音频、图像数据,创建用户界面等功能，实为一个功能强大得软件.每一次课给我一种新得体验与感受，这次课程设计一样,

20、平时都就是啃课本，很单调很乏味，但就是，这次课设给了我们动手与自主学习得机会，当然，在这过程中不免会遇到困难，当时也会很困惑、沮丧,但问题总要解决得，所以自己还就是会去想办法，我想这也就是一种收获。还有，经过了课程设计，我也深刻体会到ATLAB功能得强大以及学好MATLAB得重要性，在今后得学习生活中，我会继续学习使用它。6 参考文献1 康颖等、 大学物理M、 长沙：国防科大出版社,1962 同济大学数学研究室著、 高等数学M、北京：高等教育出版社,193 赵凯华、 新概念物理教程、 北京:高等教育出版社，1944 孙玉发、 电磁场与电磁波M、安徽：合肥工业大学,20065 郝培峰,崔建江，潘

21、峰、 计算机仿真技术M、 北京：机械工业出版社，2009 周建兴，岂兴明，矫津毅、 MLB从入门到精通、 北京:人民邮电出版社,207 李丽芬,云彩霞,郝鹏伟、 利用MATLA图型技术实现电磁波得可视化J、 信息与电脑，216:24、M文件：clea；clcu0=i7; %自由空间中得磁导率e0=1e9/(36i）; 自由空间中得电介质常数0=(u0/e0)0、5; 自由空间中得波阻抗f=e8; %电磁波得频率w2pif; =w（u0*0）0、5; 波数h_E=0； 初始相位设为0p_H=0；=0； %电场振幅H=EEZ0； 磁场振幅x=0：、：0；m0=zeros(siz(x）；freort=0:：30 z=EE*cos(kx-te-+piE); HyHH*c(x-*1e9+phi_); plo（,m0，Ez，b,ieWidth，2）； hold on;grido; plot(，H，m,r，Lnth，2); ho of xael(传播方向) ylabel(磁场Hy） labe(电场Ez) title(平面电磁波传播示意图,t=，u2sr（t），ns,fonsize，14) dwow use(0、01)e