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

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

2、。麦克斯韦方程组的解既是空间的函数又 是时间的函数，而我们只考虑随时间按正弦函数变化的解的形式。寸于这种解，其 形式可表示成一个与时间无关的复矢量和一个约定时因子exp j t相乘，这里 是 角频率。在这种约定下，麦克斯韦方程组便可表示成1(1)(1)HjE（2 2）E0 0（3 3）H0 0（4 4）寸方程（ 1 1 ）两边同取旋度，并将式 （2 2）代入便得E2E（ 5 5）利用如下矢量拉普拉斯算子定义以及方程（ 3 3）6 6）方程（ 5 5）式变为2Ek2E0(7(7)3k k类似地，可得B所满足的方程为2Bk2B0方程（7 7）和（9 9）式称为亥姆霍兹（HelmholtzHelmh

3、oltz）方程，是电磁场的波动方程。2 2、平面波解一般的电磁波总可用傅里叶分析方法展开成一系列。单色平面波的叠加。所以，对单色平面波的研究具有重要的理论和实际意义。假定波动方程（7 7）和（8 8）式的单色平面波的复式量解为3Egexp jB Bgexo j t k r（1111）式中E0, ,B0分别为E, ,B振幅，为圆频率，k k 为波矢量（即电磁波的传播方向）。exp j kx t代表波动的相位因子。为了描述均匀平面波的相位在空间的变化快慢，在此引入相速的概念，即平面波等相位的传播速度。很显然等相位面由下面方程决定1t t krkr constconst方程（1212）两边对时间 t

4、 t 求导可得drvdt k由式（8 8）可知1 1(8)(8)(9(9)(10(10)( (13)(13)(14(14)EokBo(15(15)4V厂将（1010）和（1111）式代入我们上面给出的麦克斯韦方程组可得3EokBo(15(15)5Bok kEokBoO（1818）由（1717）和（1818 ）可以看出，介质中传播的电磁波是横波，电场与磁场都与传播方向垂直；由（1515）和（1616）式可知：Eo，Bo与 k k 三者相互垂直，且满足右手螺旋关系。3 3、电磁波在线性介质中的传播1电磁波在线性介质中的传播，即电介质参数和磁导率都为实数的波传播情况。由 关系式（8 8）可知，波数

5、k k 必为实数。根据平面波解形式（1O1O）易知，平面电磁波在线 性介质中传播，只有相位发生变化，无幅值变化。将式（1515）写成kEH（1919）其中。而且 的单位是 ，故称为波阻抗。其物理意义是垂直于传播k方向平面上的电场和磁场的比值。在线性介质中，波阻抗为实数，也就是纯电阻，所以电场和磁场同相。4 4、电磁波在非线性介质中的传播1实际中见到的非线性介质是电介质参数为复数的情形，即j，譬如海水、湿地。通常这种介质的损耗是由电导率引起，故又有一。根据关系式（8 8）有Eo(16(16)(17(17)EokBo(15(15)6将复数 k k 写成k j由式（2O2O）不难推出(21(21)7

6、,_ 1/2, 2 1 1 1 12、由此可知，平面电磁波在非线性介质中传播，除了相位以传播常数幅值也要以衰减常数 随距离指数衰减。此时波阻抗为(24(24)由此可知，在非线性介质中，一般来说电场和磁场不再同相。下面我们分弱耗和良导体 中两种情况进行讨论。在弱耗情况下，即 一102，式（2222），（2323），（ 2424）可近似为I I(25(25)(26(26)由此可知，在弱耗情况下，传播常数 与在线性介质中传播下相同，衰减常数与频率102，式（2222），（ 2323），（2424）可近似为（2828）1/2(22(22)(23(23)随距离变化外，其无关，电场和磁场同相。在良导体下，

7、即8（2929）（3030）9由式(3030)可知， 在良导体中， 电场和磁场不在同相， 而是电场始终超 前磁场_。 由式(2929)可知，电磁波在良导体中传播衰减很快，很难深入到4良导体内部。一般电磁场能量集中于良导体表面。为此定义一个趋附深度，描述电磁波穿透导体的能力，具体定义式是1(3131)三、MATLA编程步骤设定相关物理量的符号表示;2 2、编好电场、磁场的表达公式;3 3、根据点电场强度、磁场强度公式运用 MATLABMATLAB 的相关函数 plot,meshplot,mesh 等模拟出电磁波在介质中传播的图像；通过编写 MATLABMATLAB 程序，我们可以生成相关的 MA

8、TLABMATLAB 图像进行可视 化，得到图形。四、实验内容电磁波在平面上传播 程序：t=0:0.2:4*pi;T=meshgrid(t);Z=si n( T);surf(Z);MATLABMATLAB 图像：即为电磁0.370.37 时，所传播的厚度1102 2、电磁波在理想介质中传播程序一：grid on ; %打开网格x=0:0.2:30;zero=0* on es(size(x);E=on es(size(x)*0;H=on es(si n(x);t=0;%画动画%for i=1:100E=exp(0*x).*cos(20*pi*t-x);就是无损耗了H=exp(0*x).*cos(

9、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, fon tsize ,20);xlabel(x, fontSize,20);ylabel(y, fontSize,20);zlabel(z, fontSize,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 ;%标注想 x， y，z 轴%

10、限定图像范围%调整视角 %帧延时 %时间流逝%关闭保持%结束循环MATLABMATLAB 图像一:6060404020200 00.50-0.5-18080%动画帧数%电场表达式0.05改为 0 ,Y);%磁场表达式%画电场矢量图R);color,k);%画磁场矢量图理想介质中电磁波传播11-2理想介质中电磁波传播2010-2-23025201510理想介质中电磁波传播12程序二:clearm=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

11、.); 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=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;2 1 0-1-2 M1015-22025123013endMATLABMATLAB 图像二:3 3、电磁波在导体(损耗较小的介质)中传播 程序与步骤二大题相同，区别在于：E=exp(-

12、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图像：理想介质中电磁波传播3-1 0%电场表达式的衰减系数%磁场表达式的衰减系损耗较小介质中电磁波传播ii损耗较小介质中电磁波传播z-22001210损耗较小介质中电磁波传播iiy损耗较小介质中电磁波传播16y4 4、电磁波在金属(损耗较大的介质)中传播程序与步骤二大题相同，区别在于：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图像：%电场表达式的衰减系数%磁场表达式的衰减系数17电磁波在金属中传播电磁波在金属中传播21z0-1-23021z0-1-2218电磁波在金属中传播五、实验总结在以前的学习中，我仅只是使用 MATLABMATLAB 勺数值计算的功能，通过这个实验， 对于 MATLAMATLA 强大的仿真功能有了更加深刻的了解， 为深层次的学习此软件开 了一个很好的头。通过 MATLAMATLA 画出的电磁波在介质中的传播能加深我们对电 场、磁场的了解，在画图的过程中，我明白了电磁波在介质中传播是有损耗 的，在实际生活中，我们已经离不开电磁波了，电磁波技