电磁场大作业

电磁场大作业

一、画出线极化、圆极化、椭圆极化平面波图形

1.线极化

Y

X

Y

clc;clear;close;

t=0:0.01:4*pi;

z=0:0.01:4*pi;

Exm=1;Eym=3;

w=2;k=1;Fy=pi;

%%线极化

subplot(1,2,1)

Ex=Exm*cos(w*t-k*z+Fy); Ey=Eym*cos(w*t-k*z+Fy); plot3(Ex,z,Ey)

title('Ïß¼«»¯')

grid Xlabel('X')

Ylabel('Z')

Zlabel('Y')

subplot(1,2,2)

Ex=Exm*cos(w*t+Fy); Ey=Eym*cos(w*t+Fy); plot(Ex,Ey)

title('Ïß¼«»¯')

grid

Xlabel('X')

Ylabel('Y')

2.圆极化

Y

X

Y

clc;clear;close;

t=0:0.01:4*pi;

z=0:0.01:4*pi;

Em=2;

w=2;k=1;Fy1=pi;Fy2=2*pi; %%圆极化¯

subplot(1,2,1)

Ex=Em*cos(w*t-k*z+Fy1); Ey=Em*sin(w*t-k*z+Fy2); plot3(Ex,z,Ey)

grid Xlabel('X')

Ylabel('Z')

Zlabel('Y')

subplot(1,2,2)

Ex=Em*cos(w*t+Fy1); Ey=Em*sin(w*t+Fy2); plot(Ex,Ey)

grid

Xlabel('X')

Ylabel('Y')

3.椭圆极化

Y

X

Y

clc;clear;close;

t=0:0.01:4*pi;

z=0:0.01:4*pi;

Emx=2;Emy=3;

w=2;k=1;Fy1=3;Fy2=4;

%%椭圆极化

subplot(1,2,1)

Ex=Emx*cos(w*t-k*z+Fy1); Ey=Emy*sin(w*t-k*z+Fy2); plot3(Ex,z,Ey)

grid Xlabel('X')

Ylabel('Z')

Zlabel('Y')

subplot(1,2,2)

Ex=Emx*cos(w*t+Fy1); Ey=Emy*sin(w*t+Fy2); plot(Ex,Ey)

grid

Xlabel('X')

Ylabel('Y')

二、画出线极化平面电磁波垂直入射到理想介质分界面之后

1.图

取n1=3，n2=7,即反射系数为0.4，透射系数为1.4.选取k1=7，k2=3.

Magnetic Field

分解图

合成图

2.代码

clc;clear;close;

z=-2*pi:0.01:0;

Z=0:0.01:2*pi;

x=zeros(1,2*pi*100+1);

y=x;

k1=7;k2=3;

n1=3;n2=7;

T=2*n2/(n1+n2);%%反射系数=T-1

Eio=10;

Ero=Eio*(T-1);

Eto=Eio*T;

figure(1)

%%电场

subplot(2,1,1)

Ei=Eio*exp(-1*j*k1*z);

Er=Ero*exp(j*k1*z);

Et=Eto*exp(-1*j*k2*Z);

%%E=Ei+Er;

E=Eio*(T*exp(-1*j*k1*z)+j*2*(T-1)*sin(k1*z));

plot(z,Ei,z,Er,'*',Z,Et,'--',z,E,'o')

title('Electronic Field')

grid on

%%磁场

subplot(2,1,2)

Hi=Ei/n1;

Hr=-1*Er/n1;

Ht=Et/n2;

%%H=Hi+Hr

H=Eio*(T*exp(-1*j*k1*z)-2*(T-1)*cos(k1*z))/n1;

plot(z,Hi,z,Hr,'*',Z,Ht,'--',z,H,'o')

grid on

title('Magnetic Field')

%%电磁场

figure(2)

plot3(y,z,Ei,y,z,Er,'*',y,Z,Et,'s',Hi,z,x,'--',Hr,z,x,'o',Ht,Z,x,'+') grid

title('Electromagnetic Field')