电磁场与电磁波-矩形波导中电力线的绘制

矩形波导中的模式以及场的分布情况

矩形波导管的横截面是矩形的空心金属波导管。在其横截面上，波导的长度为a,高度为b。假设波导中填充材料的介电常数为$\epsilon$，磁导率为$\mu$，波导可以看成是理想导体。
在矩形波导里面的传输的电磁波可以分为两种模式，即TE模式和TM 模式。当波导中传播TM电磁波的时候，磁场的分量位语电磁波传播方向的横截面内，即$H_z=0$，也叫E波；当波导中传播TE电磁波的时候，电场的分量垂直于电磁波传播的横截面内，即$E_z=0$，也叫H 波。对于波导中两种模式的公式，这里不在进行赘述，我们直接开始绘制矩形波导端面 的电力线和磁力线。

首先我们可以利用各个模式中电场和磁场的分布函数，利用quiver函数就可以进行电力线和磁力线的绘制。所谓的断面，也就是我们$z=0$的地方的力线分布。此时我们就不需要再去考虑纵向场的分布，而只要考虑横向电磁场的分布情况就好了。

1. TM模式矩形波导的电力线
   这里我们首先假设一下参数
   $$t=0,z=0,\gamma /k_c^2=1$$
   此时再利用分布函数进行进行绘图设置，代码如下：

a=2;
b=1;
E0=1;
m=2;n=2;                             %电磁波的模式是TM22
[X,Y]=meshgrid(0:0.1:2,0:0.1:1);     %这里设置波导的横截面网格
Ex=-m*pi*a*E0*cos(m*pi/a*X).*sin(n*pi/b*Y);
Ey=-n*pi/b*E0*sin(m*pi/a*X).*cos(n*pi/b*Y);
quiver(X,Y,Ex,Ey); %这里开始绘制电力线
hold on 
plot(0:0.1:2,zeros(1,21),'r','lineWidth',4);%对矩形波导进行绘制
plot(0:0.1:2,b*ones(1,21),'r','linewidth',4);
plot(zeros(1,11),0:0.1:1,'r','linewidth',4);
plot(a*ones(1,11),0:0.1:1,'r','linewidth',4);
axis image

• TM模式矩形波导的磁力线
  这里我们首先假设一下参数：
  这里我们首先假设一下参数
  $$t=0,z=0,j\omega \gamma /k_c^2=1$$

a=2;
b=1;
E0=1;
m=2;n=2;                             %电磁波的模式是TM22
[X,Y]=meshgrid(0:0.1:2,0:0.1:1);     %这里设置波导的横截面网格
Hx=n*pi/b*E0*sin(m*pi/a*X).*cos(n*pi/b*Y);
Hy=-m*pi/a*E0*cos(m*pi/a*X).*sin(n*pi/b*Y);
quiver(X,Y,Hx,Hy); %这里开始绘制磁力线
hold on 
plot(0:0.1:2,zeros(1,21),'r','lineWidth',4);%对矩形波导进行绘制
plot(0:0.1:2,b*ones(1,21),'r','linewidth',4);
plot(zeros(1,11),0:0.1:1,'r','linewidth',4);
plot(a*ones(1,11),0:0.1:1,'r','linewidth',4);
axis image

• TE模式矩形波导的电力线
  采用数值计算方法去得到横向的电场分布，并绘制电力线，部分参数的设置和上面的相同。

a=2;
b=1;
H0=1;
m=2;n=2;                             %电磁波的模式是TE22
[X,Y]=meshgrid(0:0.1:2,0:0.1:1);     %这里设置波导的横截面网格
Hz=H0*cos(m*pi/a*X).*cos(n*pi/b*Y);
[Hx,Hy]=gradient(Hz,0.1,0.1);
Hx=-Hx;Hy=-Hy;
Ex=Hy;
Ey=-Hx;
quiver(X,Y,Ex,Ey);
hold on 
plot(0:0.1:2,zeros(1,21),'r','lineWidth',4);%对矩形波导进行绘制
plot(0:0.1:2,b*ones(1,21),'r','linewidth',4);
plot(zeros(1,11),0:0.1:1,'r','linewidth',4);
plot(a*ones(1,11),0:0.1:1,'r','linewidth',4);
axis image

• TE模式矩形波导的磁力线

a=2;
b=1;
H0=1;
m=2;n=2;                             %电磁波的模式是TE22
[X,Y]=meshgrid(0:0.1:2,0:0.1:1);     %这里设置波导的横截面网格
Hz=H0*cos(m*pi/a*X).*cos(n*pi/b*Y);
[Hx,Hy]=gradient(Hz,0.1,0.1);
Hx=-Hx;Hy=-Hy;
quiver(X,Y,Hx,Hy);
hold on 
plot(0:0.1:2,zeros(1,21),'r','lineWidth',4);%对矩形波导进行绘制
plot(0:0.1:2,b*ones(1,21),'r','linewidth',4);
plot(zeros(1,11),0:0.1:1,'r','linewidth',4);
plot(a*ones(1,11),0:0.1:1,'r','linewidth',4);
axis image