a=0.025,b=0.0125;   %矩形波導的幾何尺寸
t=0,z=0.01;      %時間常數及縱向分量的取值(取不同的縱向分量對應不同橫截面的場圖)
f=1e12;        %工作頻率
e0=8.85e-12;   %介電常數為真空中的
m=1,n=1;       %傳播模式TM11
kc2=(m*pi/a)^2+(n*pi/a)^2;  %截止波數
w=2*pi*f;      %角頻率
k=w*sqrt(e0*4*pi*1e-7);  %波數
y=sqrt(kc2-k^2);        %傳播常數
[X,Y]=meshgrid(0:0.005:0.025,0:0.0015:0.0125); 
Ex=-(y/kc2)*(m*pi/a).*cos(m*pi*X./a).*sin(n*pi.*Y./b).*exp(j*w*t-y*z);     %電場的橫向分量
Ey=-(y/kc2)*(n*pi/b).*sin(m*pi.*X./a).*cos(n*pi.*Y./b).*exp(j*w*t-y*z);     %電場的縱向分量
qe1=sqrt(Ex.^2+Ey.^2);   
Ex=Ex./qe1;
Ey=Ey./qe1;
quiver(X,Y,Ex,Ey,0.5);   %繪制電場分布的矢量圖
hold on;
Hx=(j*w*e0/kc2)*(n*pi/b).*sin(m*pi.*X./a).*cos(n*pi.*Y./b).*exp(j*w*t-y*z);  %磁場的橫向分量
Hy=-(j*w*e0/kc2)*(m*pi/a).*cos(m*pi.*X./a).*sin(n*pi.*Y./b).*exp(j*w*t-y*z); %磁場的縱向分量
qe=sqrt(Hx.^2+Hy.^2);   
Hx=Hx./qe;
Hy=Hy./qe;
quiver(X,Y,Hx,Hy,0.5,'--r');  %繪制磁場分布的矢量圖
legend('電力線分布','磁力線分布')
xlabel('X(單位：米)');
ylabel('Y(單位：米)');
hold off;