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

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

2、的函数，而我们只考虑随时间按正弦函数变化的解的形式。对于这种解，其 形式可表示成一个与时间无关的复矢量和一个约定时因子 exp j t 相乘，这里 是 角频率。在这种约定下，麦克斯韦方程组便可表示成 1E j （1） j （2） 0（3） 0（4）对方程（ 1）两边同取旋度，并将式 （2）代入便得 2 （ 5）利用如下矢量拉普拉斯算子定义以及方程（ 3）2 （6）方程（ 5）式变为 2222 k 2 0 （7）k （ 8）类似地，可得 所满足的方程为222 k 2 0（9）方程（ 7）和（ 9）式称为亥姆霍兹（ Helmholtz）方程，是电磁场的波动方程。2、平面波解一般的电磁波总可用傅里叶

3、分析方法展开成一系列。单色平面波的叠加。所以，对 单色平面波的研究具有重要的理论和实际意义。假定波动方程（ 7）和（ 8）式的单色平 面波的复式量解为 3 0 ex pj t k r（ 10） 0 ex pj t k r（ 11）式中 0 , 0分别为 , 振幅， 为圆频率， k 为波矢量（即电磁波的传播方向）exp j kx t 代表波动的相位因子。为了描述均匀平面波的相位在空间的变化快慢，在此引入相速的概念，即平面波等12)(13)14)相位的传播速度。很显然等相位面由下面方程决定 1t kr const方程（ 12）两边对时间 t 求导可得dr dt k由式（ 8）可知将（10）和（ 1

4、1）式代入我们上面给出的麦克斯韦方程组可得30 k 2 0 v2（16）k 0 0（17）k 0 0（18）0 k015)由（17）和（18）可以看出，介质中传播的电磁波是横波，电场与磁场都与传播方向垂直；由（ 15）和（ 16）式可知： 0，0与k 三者相互垂直，且满足右手螺旋关系。3、电磁波在线性介质中的传播 1电磁波在线性介质中的传播，即电介质参数和磁导率都为实数的波传播情况。由 关系式（ 8）可知，波数 k 必为实数。根据平面波解形式（ 10）易知，平面电磁波在线 性介质中传播，只有相位发生变化，无幅值变化。将式（ 15）写成k （ 19）其中 。而且 的单位是 ，故称为波阻抗。其物理

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

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

7、定相关物理量的符号表示；2、编好电场、磁场的表达公式；3、根据点电场强度、磁场强度公式运用 MATLAB 的相关函数 plot,mesh 等模拟出电磁波在介质中传播的图像；通过编写 MATLAB 程序，我们可以生成相关的 MATLAB 图像进行可视 化，得到图形。四、实验内容1、电磁波在平面上传播 程序： t=0:0.2:4*pi;T=meshgrid(t);Z=sin(T);surf(Z);MATLAB 图像：60604040202010.50-0.5-180002、电磁波在理想介质中传播程序一：grid on ; %打开网格 x=0:0.2:30;zero=0*ones(size(x);E

8、=ones(size(x)*0;H=ones(sin(x);t=0;%画动画 %for i=1:100E=exp(0*x).*cos(20*pi*t-x); 就是无损耗了H=exp(0*x).*cos(20*pi*t-x-3*pi/8); quiver3(x,zero,zero,zero,zero,E, hold on ;quiver3(x,zero,zero,zero,H,zero, ti=title( 无损耗介质中电磁波传播 set(ti, fontsize ,20);xlabel(x, fontSize,20);ylabel(y, fontSize,20);zlabel(z, fontS

9、ize,20);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 ;%动画帧数%电场表达式 0.05 改为 0 ，Y );%磁场表达式%画电场矢量图R );color, k );%画磁场矢量图%标注想 x， y，z 轴%限定图像范围%调整视角 %帧延时 %时间流逝 %关闭保持%结束循环MATLAB 图像一：理想介质中电磁波传播理想介质中电磁波传播10-1-2-2101510152202530理想介质中电磁波传播y 程序二：clear m=3; x=(0:0.01:1) *m;

10、 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 in data=data(end),data(1:end-1);else data=sin(2*pi*x(i),data(1:end-1);end set(hy,YData,data); set(hz,ZData,data); drawnow pause(0.02) i=i+1;endMATLAB

11、图像二：%电场表达式的衰减系数%磁场表达式的衰减系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 图像：10损耗较小介质中电磁波传播损耗较小介质中电磁波传播zy11损耗较小介质中电磁波传播y4、电磁波在金属(损耗较大的介质)中传播程序与步骤二大题相同，区别在于 ：%电场表达式的衰减系数%磁场表达式的衰减系数E=exp(-0.1*x).*cos(20*pi*t

12、-x);由 0 改为 -0.1H=exp(-0.1*x).*cos(20*pi*t-x-3*pi/8);由 0 改为 -0.1MATLAB 图像：12电磁波在金属中传播-212001102电磁波在金属中传播0-210102012301321z0-1-23021z0-1-2210101-220-21030120y电磁波在金属中传播五、实验总结在以前的学习中，我仅只是使用 MATLAB的数值计算的功能，通过这个实验， 对于 MATLAB强大的仿真功能有了更加深刻的了解， 为深层次的学习此软件开 了一个很好的头。 通过 MATLAB画出的电磁波在介质中的传播能加深我们对电 场、磁场的了解，在画图的过程中，我明白了电磁波在介质中传播是有损耗 的，在实际生活中，我们已经离不