《电磁场论》实验课（二）.docVIP

实验课（二） 实验课（二）为关于自由空间和有界空间电磁波问题的计算和图形演示的综合实验，涵盖了自由和有界空间电磁波问题的基本理论和计算公式，加深对电磁波在自由和有界空间传播的电磁特性的理解和掌握。借助专业数学软件Mathematica对电磁学问题进行编程求解和形象化图形展示, 以综合培养学生应用数学软件进行科学计算和形象化演示的能力和意识。 实验要求： 1. 熟悉自由和有界空间电磁波场的基本理论和计算公式； 2. 掌握电磁波传播特性和波导结构电磁场分布的计算方法； 3. 进一步了解专业数学软件(Mathematica)的编程和图形演示功能。 实验内容： 实验1、平面电磁波的叠加 设有两沿z轴方向极化的单色平面电磁波沿正x轴传播，电场幅度分别为E01=2*E02=1vm, 频率分别f1和f2,变化范围为0.1～2，初相位分别为?1和?2, 变化范围为0～2π，动态演示两电磁波及其叠加波在空间原点处的时域波形随频率和初相位的变化。 解：设该两平面波在原点处的电场分别为： ， 则根据波的叠加原理，叠加波的表达式为： 程序代码： Clear[Global`*]; E01=1; E02=E01/2; E1[t_,f1_,theta1_]:=E01*Sin[2*Pi*f1*t+theta1]; E2[t_,f2_,theta2_]:=E02*Sin[2*Pi*f2*t+theta2]; Animate[Column[{Plot[{E1[t,f1,theta1],E2[t,f2,theta2]},{t,-2*Pi,2*Pi},PlotRange-2,AxesLabel-{t,E1,E2}],Plot[E1[t,f1,theta1]+E2[t,f2,theta2],{t,-2*Pi,2*Pi},PlotRange-2,AxesLabel-{t,E}]}],{f1,0.1,2},{f2,0.1,2},{theta1,0,2*Pi},{theta2,0,2*Pi},AnimationRunning-False] 输出结果： 注意：可以通过鼠标调节电磁波的频率和初相位。 练习题：上述为电磁波在空间原点处的时域波形变化，尝试修改程序演示在t=0时刻电磁波的空间波形的变化情况。其中E1和E2的表达式变成：，, 波矢 k1和 k2的变化范围为1～10，初相位分别为θ1和θ2, 变化范围为0～2π，空间区域范围z为0～10。 实验2、电磁波的极化 已知真空中极化电磁波的电场为，试求此电磁波的磁场，瞬时坡印廷矢量和时间平均坡印廷矢量。当, , , ,,时，绘图演示电场矢量随时间的变化情况。 解：该电磁波电场的复数形式为： 其对应的磁场计算公式为： 瞬时坡印廷矢量为： 时间平均坡印廷矢量： 程序代码： Clear[Global`*]; Ee={Ex*Exp[I*(k*z-w*t+phix)],Ey*Exp[I*(k*z-w*t+phiy)],0}; Needs[VectorAnalysis`]; SetCoordinates[Cartesian[x,y,z]]; (*磁场矢量*) H=-1/(I*w*mu0)*Curl[Ee]; Print[H=,H] Ereal=Re[Ee]//ComplexExpand; Hreal=Re[H]//ComplexExpand; (*瞬时坡印廷矢量*) S=Cross[Ereal,Hreal]//Simplify; Print[S=,S] (*时间平均坡印廷矢量*) Sav=1/2*Re[Cross[Ee,H/.Complex[0,n_]?-Complex[0,n]]]//ComplexExpand//Simplify; Print[Sav=,Sav] (*电场作图*) Ex=1; Ey=2; phix=0;phiy=Pi/2; z=0; w=2*Pi; Exreal[t_]=Ereal[[1]]; Eyreal[t_]=Ereal[[2]]; Animate[Graphics[{Red,Arrow[{{0,0},{Exreal[t],Eyreal[t]}}]},PlotRange?{{-2,2},{-2,2}},AspectRatio?Automatic],{t,0,1},AnimationRunning?False] 输出结果： 练习题： 上述电磁波为何种极化波？尝试修改程序中的参数，使得该电磁波变为左旋圆极化波，并用动画图形演示之。 实验3、电磁波在介质界面处的反射和透射 如下图所示有一将两个理想非磁性介质(μ＝μ0)分开的平界面：左边介质1的折射率为n1=2; 右边介质2的折射率为n2。试根据垂直极化入射电磁波下的菲涅尔公式，分别画出该入射