垂直极化波斜入射到分界面、左旋圆极化波传播动态图、矩形波导中TE01模TM11模传播动态图仿真报告

• 垂直极化波斜入射到两种理想介质分界面

图像如下：

• 左旋圆极化波传播动态图

图像如下：

• 矩形波导中TE01模TM11模传播动态图

1)、TE01

2)、TM11

• 实现代码

1. 实验一垂直极化波斜入射到两种理想介质分界面MATLAB仿真-其它文档类资源-CSDN下载

%题目要求电场入射到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



。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。

    

1. 实验二matlab电磁场左旋圆极化波传播动态图-行业报告文档类资源-CSDN下载

%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模TM11模传播动态图MATLAB仿真-其它文档类资源-CSDN下载

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);


。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。