- 垂直极化波斜入射到两种理想介质分界面
图像如下:
- 左旋圆极化波传播动态图
图像如下:
- 矩形波导中TE01模TM11模传播动态图
1)、TE01
2)、TM11
- 实现代码
%题目要求电场入射到xo平面
clear;
clc;
u0=4*pi*1e-7; %自由空间中的磁导率
e0=1e-9/(36*pi); %自由空间中的电介质常数
f=1e8; %电磁波的频率
w=2*pi*f;
Ei=8*1.41; %入射波幅度
R=0.268; %反射系数
T=1.268; %透射系数
Er=Ei*R; %反射波幅值
Et=Ei*T; %透射波幅值
k1=w*(4*u0*e0)^0.5; %真空中(介质1)波数
k2=w*(2*u0*e0)^0.5; %介质2波数
gifname='E_view_40,0,0.gif'; %取图用
figure(1)
for t=0:2:100 %为了消除波数与频率之间的数量级带来的影响,时间单位为ns
%入射波部分
%===入射电场部分
xei=-40:0.1:-0.5; %入射电场采样范围及采样步长
m1=zeros(size(xei));
[yei,zei]=meshgrid(-40:0.1:-0.5,-40:0.1:-0.5); %生成y,z取样数组网络
Eyi=1*(1.73*yei+zei)+Ei*0.5*cos(w*t*1e-9-k1*0.5*(1.73*yei+zei)); %电场y方向上分量,来源教材公式6.8-22
%时域分解来源于5.7-14,后面反射,折射亦同的相同
%其中1.73*(1.73*yei+zei)相当于根号3y尖,为方向系数,是为了与画图函数plot3匹配而设置
Ezi=-1.73*(1.73*yei+zei)+Ei*0.5*cos(w*t*1e-9-k1*0.5*(1.73*yei+zei)); %电场z方向上分量
plot3(m1,Eyi,Ezi,'LineWidth',1.5); %绘制入射电场
hold on
%反射波部分
%===反射电场部分
xer=0:0.1:40;
m2=zeros(size(xer));
[yer,zer]=meshgrid(0:0.1:40,0:0.1:40);
Eyr=1*(1.73*yer+zer)+Er*0.5*cos(k1*0.5*(1.73*yer+zer)-w*t*1e-9);
Ezr=1.73*(1.73*yer+zer)-Er*0.5*cos(k1*0.5*(1.73*yer+zer)-w*t*1e-9);
plot3(m2,Eyr,Ezr,'LineWidth',1);
hold on
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
%Electromagnetic wave simulation
%假定电磁波沿x轴方向传播,电磁波在y轴方向和z轴方向有分量
Ey_m=8; %电磁波y分量幅值
Ex_m=8; %电磁波z分量幅值
E=(Ey_m^2+Ex_m^2)^0.5; %求出电磁波的总幅值
w=10; %电磁波的角速度为10
a_1=pi; %y分量电磁波的初相位
a_2=pi/2; %z分量电磁波的初相位
%a_1和a_2的大小差值,决定了相位差
%位相差与分量的幅值进一步决定了电磁波为何种程度的极化电磁波(线极化,圆极化,椭圆极化)
z=0:0.01:4; %设置沿x轴方向传播的点位
m_0=zeros(size(z)); %m_0为与x矢量同尺寸的0坐标参考点位
figure
for t=0:1000
Ey=Ey_m*cos(w*z+w*t*1e-2+a_1) %沿x轴方向传播的波,在y-x平面上的波函数
Bx=Ex_m*cos(w*z+w*t*1e-2+a_2) %沿x轴方向传播的波,在z-x平面上的波函数
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
1)、TE01
clear
a=4 ;b=3 ;
ma=linspace(0,a,15);
mb=linspace(0,b,15);
[x,y]=meshgrid(ma,mb);
m=0 ;n=1;
Ex=(n/b)*cos(m*pi*x/a).*sin(n*pi*y/b);
Ey=-(m/a)*sin(m*pi*x/a).*cos(n*pi*y/b);
Ez=0;
%微波技术原理
scale=.9;
handle=quiver(x,y,Ex,Ey,scale,'r')
set(handle,'LineWidth',1)
%plot(Hx,Hy,'r')
title('矩形波导TE01波的电磁场结构')
%title('linspace(0,a,15)')
%title('LineWidth=3')
hold on
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
2)、TM11
% 程序2 TM波
%矩形波导的电磁场线
clear
a=4 ;b=3 ;
ma=linspace(0,a,15);
mb=linspace(0,b,15);
[x,y]=meshgrid(ma,mb);
m=1; n=1;
Ex=-(m/a)*cos(m*pi*x/a).*sin(n*pi*y/b);
Ey=-(n/b)*sin(m*pi*x/a).*cos(n*pi*y/b);
scale=1;
quiver(x,y,Ex,Ey,scale,'r')
%plot(Hx,Hy,'r')
title('矩形波导TM11波的电磁场线')
hold on
Hx1=(n/b)*sin(m*pi*x/a).*cos(n*pi*y/b);
Hy1=-(m/a)*cos(m*pi*x/a).*sin(n*pi*y/b);
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。