- 1、本文档共14页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
32钟伦利科技创新概论报告
湖南工程学院
科技创新概论报告
课程名称 电磁场模拟
专 业 自动化
班 级 1191班
学 号 201101029132
姓 名 钟伦利
2014年 6月 16 日
目 录
一、数值模拟问题 1
二、理论算法 2
三、程序及模拟结果 6
四、对电磁场模拟的理解 8
附录:实验Matlab程序 10
一、数值模拟问题
第二题:一个二维正方形静电场域电位函数为φ(x,y),边界条件如图1.1所示,试用超松弛法确定二维静电场域的电位分布。
二、理论算法
本课题用MATLAB实现了有限差分法解静电场边值问题的算法,将偏微分方程的问题化为线性方程组问题,并使用了迭代法进行线性方程组的数值解。经典场的边值问题在数学上表达为泊松方程和拉普拉斯方程,但解偏微分方程往往是困难的。幸而很多时候对于具体问题我们需要的不是解析解,而是数值解,所以可以考虑用连续变量离散化的方法求出数值解,在足够的精度上进行逼进,这就引出了有限差分法。
有限差分法(Finite Differential Method)是基于差分原理的一种数值计算法。其基本思想:将场域离散为许多小网格,应用差分原理,将求解连续函数?的泊松方程的问题换为求解网格节点上?的差分方程组的问题。
三、程序及模拟结果
依据有限差分法的理论知识,以及程序框图,在MATLAB上编写程序(见附录),并且编译,运行,得到二维静电场域的电位分布,如图3.1和图3.2所示。
图3.1 三维静电场域的电位分布
图3.2 二维静电场域的电位分布
四、对电磁场模拟的理解
在以往的探究过程中我们发现电磁场有:理论性强、理论模型抽象、无法直观到电磁波的动态变化过程的这些难点。自从麦克斯韦预言了电磁场的基本方程以及电磁波的存在以来,电磁波就成为了人们认识自然的一种重要工具。我们通过将肉眼不可见的电磁波进行可视化并进行各方面的应用。以往我们在研究电磁场和电磁波的空间传播的时候不仅要面对大量繁琐的数学推导过程,而且还有难以绘制的三维图形,往往描绘出来的仅仅是简单的示意图。而Mat lab软件的引入,合理化地解决了这些。
电磁波是时变电磁场相互转换而来的,它在空间中传播的状态时不被我们肉眼所看到,但是可以依靠麦克斯韦方程以及边界调节来模拟出来,有一点必须要提及的是,电磁波是出于被动状态的,它每时每刻的位置都在发生着变化。这点可以由均匀平面电磁波的电场E的表达式得知:电场E是沿着正轴Z传播的,它是时间T和变量Z的函数。行波的动态模拟均匀平面电磁波在无线大理想介质中传播时,能量始终都沿着一个确定的方向传播着,称之为行波。首先依据EH的参数建立三维坐标系。因为电场、磁场的方向始终不变,所以使电场的图形与X轴平行,磁场的图形与Y轴平行。然后,分别以T和Z的自变量,EH的大小为因变量设置双重循环,并根据它们的数值,利用Mat lab提供的实现动态动画效果的指令和函数作出动态图形。但这样得到的效果,只能是两个图形随着时间而移动。为更形象地表现行波随着传播举例的变化,特设置两个不同质点于电场与磁场的波形上,使它们的移动规律也满足电磁场方程。这样,就能够通过图形及质点的瞬间变化的样子,再现动态行波了。
由于在入射波的作用下,理想截止I垂直入射到理想导体II表面上时,由于在入射波的作用下,理想导体表面产生了感应电流和感应电荷,这些电流和电荷激发的层次波,在介质I中形成反射波。入射波与反射波在介质I中叠加,形成驻波。行波与驻波的一个显著区别在于行波的能量是沿着一个明确的方向传播的,而驻波不能形成电磁能量的传播。动态模拟驻波形成过程的Mat lab程序要计算的是Ei(z,t)、Er(z,t)以及合成波E(z,t)在Mat lab命令运行窗口中,我们除了能够看到动态变化的行波和驻波,还能够看到不同的波色彩的变化。
附录:实验Matlab程序
clear
hx=20;
hy=20;%确定网格
v1=ones(hy,hx);
m=19;
n=19;
for i=1:hx
v1(1,i)=50;
v1(20,i)=100;
end
for j=1:hy
v(j,1)=0;
v(j,20)=100;
end
t1=(cos(pi/m)+cos(pi/n))/2;
w=2/(1+sqrt(1-t1^2));
v2=v1;
maxt=1;
1亿VIP精品文档
相关文档
最近下载
- 2023年四年级下语文期中复习·心田花开.pdf VIP
- 2023年全国统一高考英语试卷(全国甲卷)(含答案与解析).pdf
- 【课件】歌唱的情感与表达+当那一天来临+课件高中音乐人音版(2019)+必修2+歌唱.pptx VIP
- 03s702钢筋混凝土化粪池图集(高清).pdf
- 在线网课学习课堂《生命科学基础(北理 )》单元测试考核答案.docx
- 电视画面编辑(本科)王晓红.ppt
- 2021-2022学年威海市乳山市八年级下学期期中数学试卷(含答案解析).docx
- 特种设备生产(含安装、改造、维修)单位质量安全风险管控清单 .docx
- 人工智能导论 人工智能理论基础 人工智能理论基础.pptx
- 搅拌设备设计.ppt
文档评论(0)