试验二电磁波在介质中的传播规律

1、电磁场与微波技术实验报告课程实验：电磁波在介质中传播规律班 级姓 名2015.11.21指导老师:实验日期：电磁波在介质中的传播规律一、实验目的:1、用MATLAB程序演示了电磁波在无损耗、较小损耗和较大损耗情况下的传播博规律；2、结合图像探讨了电磁波在有耗介质中电场强度和磁场强度的能量变化情况；3、 学会使用Matlab进行数值计算，并绘出相应的图形，运用MATLAB对其进行可视 化处理。二、实验原理1、电磁场的波动方程一般情况下，电磁场的基本方程是麦克斯韦方程，而我们讨论的介质是各向同性均匀线性的，即（0, j 0）的情形。麦克斯韦方程组的解既是空间的函数又是时间的函数，而我们只考虑随时间

2、按正弦函数变化的解的形式。对于这种解，其 形式可表示成一个与时间无关的复矢量和一个约定时因子exo j t相乘，这里 是角频率。在这种约定下，麦克斯韦方程组便可表示成1E j H(1)H j EE 0H0对方程（1）两边同取旋度，并将式（2）代入便得(5)利用如下矢量拉普拉斯算子定义以及方程（3）(6)(7)(8)(9)2 EEE方程（5）式变为2 E k2 E类似地，可得B所满足的方程为k2 B方程（7）和（9）式称为亥姆霍兹（Helmholtz ）方程，是电磁场的波动方程。2、平面波解一般的电磁波总可用傅里叶分析方法展开成一系列。单色平面波的叠加。所以，对单色平面波的研究具有重要的理论和实

3、际意义。假定波动方程（7 ）和（8）式的单色平面波的复式量解为3(10)Eg exp jB Bcexo j t k r（11）式中E0, B0分别为E, B振幅， 为圆频率，k为波矢量（即电磁波的传播方向）。exp j kx t代表波动的相位因子。为了描述均匀平面波的相位在空间的变化快慢，在此引入相速的概念，即平面波等相位的传播速度。很显然等相位面由下面方程决定1方程（12 ）两边对时间t求导可得(13)drvdt k由式（8）可知(14)其中而且 的单位是 ，故称为波阻抗。其物理意义是垂直于传播将（10 ）和（11 ）式代入我们上面给出的麦克斯韦方程组可得E°kB0(15)1B0k

4、2E0(16)vkE00(17)k B00(18)由（17）和（18）可以看出，介质中传播的电磁波是横波，电场与磁场都与传播方向垂直;由（15 ）和（16 ）式可知：E0 , B0与k三者相互垂直，且满足右手螺旋关系。3、电磁波在线性介质中的传播1电磁波在线性介质中的传播，即电介质参数和磁导率都为实数的波传播情况。由关系式（8）可知，波数k必为实数。根据平面波解形式（10）易知，平面电磁波在线性介质中传播，只有相位发生变化，无幅值变化。将式（15）写成(19)k EH方向平面上的电场和磁场的比值。在线性介质中，波阻抗为实数，也就是纯电阻，所 以电场和磁场同相。4、电磁波在非线性介质中的传播1、

5、IH实际中见到的非线性介质是电介质参数为复数的情形，即j ，譬如海水、湿地。通常这种介质的损耗是由电导率引起，故又有根据关系式（8 ）有” 1/2kJ '1 L将复数k写成kj由式（20 ）不难推出1/211 " 21 - 11/2':"22 V1 1(20)(21)(22)(23)随距离变化外，其由此可知，平面电磁波在非线性介质中传播，除了相位以传播常数幅值也要以衰减常数 随距离指数衰减。此时波阻抗为(24)由此可知，在非线性介质中，一般来说电场和磁场不再同相。下面我们分弱耗和良导体 中两种情况进行讨论。在弱耗情况下，即10 2，式（22），（ 23），（

6、 24 ）可近似为(26)(27)由此可知，在弱耗情况下，传播常数 与在线性介质中传播下相同，衰减常数与频率无关，电场和磁场同相。在良导体下，即r 102，式（22），（ 23），（ 24 ）可近似为(28)(29)(30)由式（30）可知，在良导体中，电场和磁场不在同相，而是电场始终超 前磁场一。由式（29）可知，电磁波在良导体中传播衰减很快，很难深入到4良导体内部。一般电磁场能量集中于良导体表面。为此定义一个趋附深度 描述电磁波穿透导体的能力，具体定义式是(31 )即为电磁波幅值减到原来的e 10.37时，所传播的厚度三、MATLAB编程步骤1、设定相关物理量的符号表示;2、编好电场、磁场

7、的表达公式;3、根据点电场强度、磁场强度公式运用MATLAB的相关函数plot,mesh 等模拟出电磁波在介质中传播的图像；通过编写MATLAB程序，我们可以生成相关的 MATLAB图像进行可视 化，得到图形。四、实验内容1、电磁波在平面上传播程序：t=0:0.2:4*pi;T=meshgrid(t);Z=si n( T);surf(Z);MATLAB 图像：10.50-0.5-180iiUHil r!liIII6060404020202、电磁波在理想介质中传播程序pause(0.002)%帧延时grid on ;%打开网格 x=0:0.2:30;zero=0* on es(size(x);E

8、=on es(size(x)*0;H=on es(si n(x);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=ti tle('无损耗介质中电磁波传播','col

9、or' ,'k');set(ti, 'fontsize',20);xlabel( 'x' ,'fontSize' ,20);ylabel( 'y' ,'fontSize' ,20);zlabel( 'z' ,'fontSize' ,20);%标注想x，y，z轴axis(0,30,-2.5,2.5,-2.5,2.5);%限定图像范围view(20+2*i,40);%调整视角t=t+0.004;%时间流逝30hold off ;end ;MATLAB%关闭保持%结

10、束循环图像一:21z 0-1-2程序二:理想介质中电磁波传播21-1-2-21015202530理想介质中电磁波传播clear 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 = len gth(x);i=1;view(3);while 1if i> ndata=data(e nd),data(1:e nd-1); elsedata=

11、si n(2*pi*x(i),data(1:e nd-1);endset(hy,'YData',data);set(h z,'ZData',data);draw nowpause(0.02)i=i+1;endMATLAB 图像二:3、电磁波在导体(损耗较小的介质)中传播 程序与步骤二大题相同，区别在于：E=exp(-0.03*x).*cos(20*pi*t-x);-0。03H=exp(-0.03*x).*cos(20*pi*t-x-3*pi/8);-0。03%电场表达式的衰减系数由%磁场表达式的衰减系数由0改为0改为MATLAB 图像:损耗较小介质中电磁波传播-

12、212001210损耗较小介质中电磁波传播y损耗较小介质中电磁波传播y4、电磁波在金属(损耗较大的介质)中传播程序与步骤二大题相同，区别在于：E=exp(-0.1*x).*cos(20*pi*t-x);-0.1H=exp(-0.1*x).*cos(20*pi*t-x-3*pi/8);%电场表达式的衰减系数由0改为-0.1MATLAB 图像:%磁场表达式的衰减系数由0改为21z 0-1-23021z 0-1-2电磁波在金属中传播0x电磁波在金属中传播210-1-20五、实验总结在以前的学习中，我仅只是使用MATLAB的数值计算的功能，通过这个实验, 对于MATLAB强大的仿真功能有了更加深刻的了解，为深层次的学习此软件 开了一个很好的头。通过MATLAB画出的电磁波在介质中的传播能加深我们 对电场、磁场的了解，在画图的过程中，我明白了电磁波