1、1电 磁 场 与 微 波 技 术 实 验 报 告（二）课程实验：电磁波在介质中传播规律班 级： 姓 名： 指导老师： 实验日期： 2015.11.21 2电磁波在介质中的传播规律1、实验目的：1、用 MATLAB 程序演示了电磁波在无损耗、较小损耗和较大损耗情况下的传播博规律；2、结合图像探讨了电磁波在有耗介质中电场强度和磁场强度的能量变化情况；3、学会使用 Matlab 进行数值计算，并绘出相应的图形，运用 MATLAB 对其进行可视化处理。2、实验原理1、电磁场的波动方程一般情况下，电磁场的基本方程是麦克斯韦方程，而我们讨论的介质是各向同性均匀线性的，即（ ）的情形。麦克斯韦方程组的解既是
2、空间的函数又0,j是时间的函数，而我们只考虑随时间按正弦函数变化的解的形式。对于这种解，其形式可表示成一个与时间无关的复矢量和一个约定时因子 相乘，这里tjexp是角频率。在这种约定下，麦克斯韦方程组便可表示成 1(1) Ej(2) j(3)0(4)对方程（1）两边同取旋度，并将式(2)代入便得（5）2利用如下矢量拉普拉斯算子定义以及方程（3）（6）2方程（5）式变为 3（7） 02k（8）类似地，可得 所满足的方程为（9）02k方程（7）和（9）式称为亥姆霍兹（Helmholtz ）方程，是电磁场的波动方程。2、平面波解一般的电磁波总可用傅里叶分析方法展开成一系列。单色平面波的叠加。所以，对
3、单色平面波的研究具有重要的理论和实际意义。假定波动方程（7）和（8）式的单色平面波的复式量解为 3rktjexp0（10） rtj0（11）式中 , 分别为 , 振幅， 为圆频率， 为波矢量（即电磁波的传播方向） 。0k代表波动的相位因子。tkxjep为了描述均匀平面波的相位在空间的变化快慢，在此引入相速的概念，即平面波等相位的传播速度。很显然等相位面由下面方程决定 1constkrt（12）方程（12）两边对时间 求导可得t(13)kdtrv由式（8）可知1v（14）将（10）和（11）式代入我们上面给出的麦克斯韦方程组可得 34（15）00k0201v（16）（17）0k（18）0由（17
4、）和(18) 可以看出，介质中传播的电磁波是横波，电场与磁场都与传播方向垂直；由（15）和（16）式可知： ， 与 三者相互垂直，且满足右手螺旋关系。0k3、电磁波在线性介质中的传播 1电磁波在线性介质中的传播，即电介质参数和磁导率都为实数的波传播情况。由关系式（8）可知，波数 必为实数。根据平面波解形式（ 10）易知，平面电磁波在线k性介质中传播，只有相位发生变化，无幅值变化。将式（15）写成（19）其中 。而且 的单位是 ，故称为波阻抗。其物理意义是垂直于传播k方向平面上的电场和磁场的比值。在线性介质中，波阻抗 为实数，也就是纯电阻，所以电场和磁场同相。4、电磁波在非线性介质中的传播 1实
5、际中见到的非线性介质是电介质参数为复数的情形，即 ，譬如海水、“j湿地。通常这种介质的损耗是由电导率 引起，故又有 。根据关系式（ 8）有“（20）2/1“jk将复数 写成k（21）jk由式（20）不难推出5（22）2/12“ 12（23）2/12“12由此可知，平面电磁波在非线性介质中传播，除了相位以传播常数 随距离变化外，其幅值也要以衰减常数 随距离指数衰减。此时波阻抗为（24）2/1 “j由此可知，在非线性介质中，一般来说电场和磁场不再同相。下面我们分弱耗和良导体中两种情况进行讨论。在弱耗情况下，即 ，式（22） ， （23） ， （24）可近似210“为（25）（26）“2（27）由此
6、可知，在弱耗情况下，传播常数 与在线性介质中传播下相同，衰减常数 与频 率无关，电场和磁场同相。在良导体下，即 ，式（22） ， （23） ， （24）可近似为2“10（28） 2“（29）（30） 21j6由式（30）可知，在良导体中，电场和磁场不在同相，而是电场始终超前磁场 。由式（29）可知，电磁波在良导体中传播衰减很快，很难深入到4良导体内部。一般电磁场能量集中于良导体表面。为此定义一个趋附深度 ，描述电磁波穿透导体的能力，具体定义式是（31）1即为电磁波幅值减到原来的 0.37 时，所传播的厚度。e3、MATLAB 编程步骤1、设定相关物理量的符号表示；2、编好电场、磁场的表达公式；
7、3、根据点电场强度、磁场强度公式运用 MATLAB 的相关函数plot,mesh 等模拟出电磁波在介质中传播的图像；通过编写 MATLAB 程序，我们可以生成相关的 MATLAB 图像进行可视化，得到图形。四、实验内容1、电磁波在平面上传播程序：t=0:0.2:4*pi;T=meshgrid(t);Z=sin(T);surf(Z);MATLAB 图像：70 2040 6080020406080-1-0.500.512、电磁波在理想介质中传播程序一：grid on;%打开网格x=0:0.2:30;zero=0*ones(size(x);E=ones(size(x)*0;H=ones(sin(x)
8、;t=0;%画动画%for i=1:100 %动画帧数 E=exp(0*x).*cos(20*pi*t-x); %电场表达式 0.05 改为0，就是无损耗了H=exp(0*x).*cos(20*pi*t-x-3*pi/8); %磁场表达式quiver3(x,zero,zero,zero,zero,E,Y); %画电场矢量图hold on; quiver3(x,zero,zero,zero,H,zero,R); %画磁场矢量图ti=title(无损耗介质中电磁波传播, color,k);set(ti,fontsize,20);xlabel(x,fontSize,20);ylabel(y,font
9、Size,20);zlabel(z,fontSize,20); %标注想 x，y，z 轴 axis(0,30,-2.5,2.5,-2.5,2.5); %限定图像范围 view(20+2*i,40); %调整视角pause(0.002) %帧延时t=t+0.004; %时间流逝 hold off; %关闭保持end; %结束循环MATLAB 图像一：80102030-2-1012-2-1012x中中中中中中中中中中yz051015202530-2-1012-2-1012x中中中中中中中中中中yz9051015202530-2 -10 12-2-1012x中中中中中中中中中中yz程序二：clear
10、m=3;x=(0:0.01:1) *m;figure;grid on;hold on;axis(0,m,-1,1,-1,1)data = zeros(size(x);hy = stem(x,data,y.);hz = stem(x,data,r.);n = length(x);i=1;view(3);while 1if indata=data(end),data(1:end-1);elsedata=sin(2*pi*x(i),data(1:end-1);endset(hy,YData,data);set(hz,ZData,data);drawnowpause(0.02)i=i+1;10endM
11、ATLAB 图像二：0123-1-0.500.51-1-0.500.51中中中中中中中中中中3、电磁波在导体（损耗较小的介质）中传播程序与步骤二大题相同，区别在于：E=exp(-0.03*x).*cos(20*pi*t-x); %电场表达式的衰减系数由 0 改为-0。03H=exp(-0.03*x).*cos(20*pi*t-x-3*pi/8); %磁场表达式的衰减系数由 0 改为-0。03MATLAB 图像：110102030 -2-1012-2-1012y中中中中中中中中中中中中xz0102030 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5-2-1012中
12、中中中中中中中中中中中yxz120102030-2 -10 12-2-1012xy中中中中中中中中中中中中z4、电磁波在金属（损耗较大的介质）中传播程序与步骤二大题相同，区别在于：E=exp(-0.1*x).*cos(20*pi*t-x); %电场表达式的衰减系数由 0 改为-0.1H=exp(-0.1*x).*cos(20*pi*t-x-3*pi/8); %磁场表达式的衰减系数由 0 改为-0.1MATLAB 图像：130102030 -2-1012-2-1012y中中中中中中中中中xz0102030-2-1012-2-1012x中中中中中中中中中yz145、0102030-2 -1 01
13、2-2-1012xy中中中中中中中中中z实验总结： 在以前的学习中，我仅只是使用 MATLAB 的数值计算的功能，通过这个实验，对于 MATLAB 强大的仿真功能有了更加深刻的了解，为深层次的学习此软件开了一个很好的头。通过 MATLAB 画出的电磁波在介质中的传播能加深我们对电场、磁场的了解，在画图的过程中，我明白了电磁波在介质中传播是有损耗的，在实际生活中，我们已经离不开电磁波了，电磁波技术革新也与我们的生活息息相关，我们要想达到在理想条件下无损耗的传播电磁波，还需要学习更多的知识，像前辈们一样更加努力。参考文献1 盛新庆. 电磁波述论M. 北京: 科学出版社, 20072 郭硕鸿. 电动力学（第二版）M. 北京: 高等教育出版社, 20063 沙湘月, 伍瑞新. 电磁场理论与微波技术 M. 南京: 南京大学出版社, 2004
