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

1、精选文档电磁场与微波技术实验报告课程实验：电磁波在介质中传播规律级:姓 名: 指导老师: 实验日期： 2015.11.212015.11.21精选文档2 EEE(6(6)方程（5 5）式变为2MATLABMATLAB寸其进行可视化电磁波在介质中的传播规律、实验目的:1 1、用MATLAB?MATLAB?序演示了电磁波在无损耗、较小损耗和较大损耗情况下的传播博规律;2 2、结合图像探讨了电磁波在有耗介质中电场强度和磁场强度的能量变化情况;3 3、学会使用MatlabMatlab进行数值计算，并绘出相应的图形，运用 处理。二、实验原理1 1、电磁场的波动方程一般情况下，电磁场的基本方程是麦克斯韦方

2、程，而我们讨论的介质是各向同性均匀线性的，即（0,j0）的情形。麦克斯韦方程组的解既是空间的函数又是时间的函数，而我们只考虑随时间按正弦函数变化的解的形式。对于这种解，其形式可表示成一个与时间无关的复矢量和一个约定时因子ex） j t相乘，这里 是角频率。在这种约定下，麦克斯韦方程组便可表示成 1对方程（1 1）两边同取旋度，并将式（2 2）代入便得精选文档(14(14)v v将（1010）和（1111）式代入我们上面给出的麦克斯韦方程组可得3利用如下矢量拉普拉斯算子定义以及方程（3 3）精选文档(14(14)v v将（1010）和（1111）式代入我们上面给出的麦克斯韦方程组可得3(9(9)

3、(10(10)式中Eo, , Bo分别为E, , B振幅,为圆频率，k k为波矢量（即电磁波的传播方向）。(12(12)(13(13) )2 E k2 E 0类似地，可得B所满足的方程为2 B k2 B 0方程（7 7）和（9 9）式称为亥姆霍兹（HelmholtzHelmholtz）方程，是电磁场的波动方程。2 2、平面波解一般的电磁波总可用傅里叶分析方法展开成一系列。单色平面波的叠加。所以，对单色平面波的研究具有重要的理论和实际意义。假定波动方程（7 7）和（8 8）式的单色平面波的复式量解为3EE0 exp j(11(11)exp j kx t代表波动的相位因子。为了描述均匀平面波的相位

4、在空间的变化快慢，在此引入相速的概念，即平面波等相位的传播速度。很显然等相位面由下面方程决定1t t krkr constconst方程（1212）两边对时间t t求导可得dr _ dt k由式（8 8）可知精选文档由式（2020）不难推出（1717）和（1818）可以看出，（1515）和（1616）式可知:介质中传播的电磁波是横波，电场与磁场都与传播方向垂直;Eo， Bo与k k三者相互垂直，且满足右手螺旋关系。EokBo(15(15)Bok(16(16)k Eo(17(17)k Bo(18(18)3 3、电磁波在线性介质中的传播1电磁波在线性介质中的传播，即电介质参数和磁导率都为实数的波传

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

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

7、介质中传播下相同，衰减常数与频率(28(28)1/2(22(22)1/2精选文档具体定义式是由式（3030）可知，在良导体中，电场和磁场不在同相，而是电场始终超 前磁场。由式（2929）可知，电磁波在良导体中传播衰减很快，很难深入到4良导体内部。一般电磁场能量集中于良导体表面。为此定义一个趋附深度 描述电磁波穿透导体的能力，即为电磁波幅值减到原来的e 10.370.37时，所传播的厚度。三、MATLAB程步骤1 1、设定相关物理量的符号表示;2 2、编好电场、磁场的表达公式;3 3、根据点电场强度、磁场强度公式运用 MATLABMATLAB的相关函数plot,meshplot,mesh 等模拟

8、出电磁波在介质中传播的图像；通过编写MATLABMATLAB程序，我们可以生成相关的 MATLABMATLAB图像进行可视 化，得到图形。四、实验内容1 1、电磁波在平面上传播程序:t=0:0.2:4* pi; T=meshgnd(t);Z=s in (T); surf(Z);MATLABMATLAB 图像:(31(31)精选文档-1 ,8010.50-0.580400 060程序2 2、电磁波在理想介质中传播grid on ; %打开网格 x=0:0.2:30;zero=0* on es(size(x); E=on es(size(x)*0; H=on es(s in( x);t=0;%画动

9、画%for i=1:100E=ex p(0*x).*cos(20* pi*t-x);就是无损耗了%动画帧数%电场表达式0.05 改为0 ,H=ex p(0*x).*cos(20* pi*t-x-3* pi/8); quiver3(x,zero,zero,zero,zero,E, hold on ;Y）；%磁场表达式%画电场矢量图quiver3(x,zero,zero,zero,H,zero, ti=title( set(ti,Rxlabel(ylabel(无损耗介质中电磁波传播fon tsize,20);，fon tSize,fontSizey,20);,20);，color);,k）；%画磁

10、场矢量图zlabel(axis(0,30,-2.5,2.5,-2.5,2.5);view(20+2*i,40);pause(0.002)，fon tSize,20)；%标注想x,y，z轴%限定图像范围%调整视角%帧延时t=t+0.004;hold off ;end ;%时间流逝%关闭保持%结束循环MATLABMATLAB图像精选文档精选文档-201015202530理想介质中电磁波传播理想介质中电磁波传播30-1一 一也一一一卞一一一二一一S . - - _ -工;工样 .I应精选文档-201015202530精选文档理想介质中电磁波传播程序二:clear m=3;x=(0:0.01:1) *

11、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 nz 0-r-2X1 -1 -03.-202530，-z/-2X-财 5y 10 15F 20X2精选文档data=data(e nd),data(1:e nd-1); elsedata=si n(2* pi*x(i),data(1:e nd-1);endset(hy,YData,data);

12、set(h z, ZData,data); draw nowp ause(0.02)i=i+1;精选文档理想介质中电磁波传播3 3、电磁波在导体（损耗较小的介质）中传播 程序与步骤二大题相同，区别在于endMATLABMATLAB 图像二:-0.5 .-0.511-1 00.50-11精选文档%电场表达式的衰减系数%磁场表达式的衰减系E=ex p(-0.03*x).*cos(20* pi*t-x); 由0改为-0。03H=ex p(-0.03*x).*cos(20* pi*t-x-3* pi/8); 数由0改为-0。03MATLAB 图像:精选文档损耗较小介质中电磁波传播精选文档-1 -2%电

13、场表达式的衰减系数%磁场表达式的衰减系数损耗较小介质中电磁波传播4 4、电磁波在金属（损耗较大的介质）中传播 程序与步骤二大题相同，区别在于E=ex p(-0.1*x).*cos(20* pi*t-x);由0改为-0.1H=ex p(-0.1*x).*cos(20* pi*t-x-3* pi/8);由0改为-0.1MATLAB图像:/1020-11302精选文档精选文档0yxzy精选文档2233z 0-电磁波在金属中传播30 x11 五、实验总结:在以前的学习中，我仅只是使用 MATLABMATLAB勺数值计算的功能，通过这个实验， 对于MATLAMATLA强大的仿真功能有了更加深刻的了解，为深层次的学习此软件开 了一个很好的头。通过MATLAMATLA画出的电磁波在介质中的传播能加深我们对电 场、磁场的了解，在画图的过程中，我明白了电磁波在介质中传播是有损耗 的，在实际生活中，我们已经