Matlab 电磁波垂直入射、透射和反射

版本

Matlab2016a

内容

code

e1r=1;
u1r=1;
o1=0;
e2r=4;
u2r=1;
o2=0;

e0=10^-9/(36*pi);
u0=4*pi*10^-7;

f=1e9;
T=1/f;
c=3*1e8;
L1=c/f;
w=2*pi*f;

k1=w*((e1r*e0)*(u1r*u0))^(1/2);
k2=w*((e2r*e0)*(u2r*u0))^(1/2);
L2=2*pi/k2;
Eim=100;
n1=120*pi*(u1r/e1r)^(1/2);
n2=120*pi*(u2r/e2r)^(1/2);
R=(n2-n1)/(n2+n1);
Tra=1+R;


%1

% t=0:0.01*T:5*T;
% z1=linspace(-5*L1,0,length(t));
% z2=linspace(0,5*L2,length(t));
% for t=0:0.01*T:5*T;
%     Ei=Eim*cos(w*t-k1*z1);
%     Er=Eim*R*cos(w*t+k1*z1);
%     Et=Eim*Tra*cos(w*t-k2*z2);
%     plot(z1,Ei);
%     grid on;
%     hold on;
%     plot(z1,Er);
%     hold on;
%     grid on;
%     plot(z2,Et);
%     grid on;
%     hold on;
%     hold off;
%     drawnow;
% end

%2
t=0:0.01*T:5*T;
z1=linspace(-5*L1,0,length(t));
z2=linspace(0,5*L2,length(t));
for t=0:0.005*T:5*T;
    Ei=Eim*cos(w*t-k1*z1);
    Er=Eim*R*cos(w*t+k1*z1);
    Et=Eim*Tra*cos(w*t-k2*z2);
    
    Etra=Eim*Tra*cos(w*t-k1*z1);
    Est=Eim*2*R*sin(k1*z1)*cos(w*t+1/2*pi);
    
    Esum=Etra+Est;
    plot(z1,Etra);
     axis([-1.5 0 -150 150]);
    grid on;
    hold on;
    plot(z1,Est);
     axis([-1.5 0 -150 150]);
    hold on;
    grid on;
    plot(z1,Esum);
     axis([-1.5 0 -150 150]);
    hold on;
    grid on;
    hold off;
    drawnow;
end


%3
% t=0:0.01*T:5*T;
% z1=linspace(-5*L1,0,length(t));
% z2=linspace(0,5*L2,length(t));
% for t=0:0.005*T:5*T;
%     Ei=Eim*cos(w*t-k1*z1);
%     Er=Eim*R*cos(w*t+k1*z1);
%     Et=Eim*Tra*cos(w*t-k2*z2);
%     
%     E1ts=Ei+Er;
%     E1z=Eim*(1+R*R+2*R*cos(2*k1*z1)).^(1/2);  
%     
%     plot(z1,E1z);
%     axis([-1.5 0 -150 150]);
%     grid on;
%     hold on;
%     plot(z1,E1ts);
%     axis([-1.5 0 -150 150]);
%     hold on;
%     grid on;
%     hold off;
%     drawnow;
% end