media_perpendicular.m

电磁场圆极化

%电场对理想介质的垂直入射；
%
%copyright 属于研学论坛bbs.matwav.com和plasma

clear all;

gama=.5;
trans=1-gama;

zmin=-10;
zmax=0;

ztmin=0;
ztmax=10;

delz = 0.1;
e0=1;
k1=1;k2=2;
framemax = 1000;
z = [zmin:delz:zmax];
zt= [ztmin:delz:ztmax];
zL=[zmin:delz:ztmax];

yL=0*zL;
y=0*z;
z1=-1:.1:1;
y1=0*z1;
jr=1;jt=1;

for n=1:framemax;
    Ei=e0*cos(2*pi*n/framemax-k1*z);
        Er=-gama*e0*cos(2*pi*n/framemax+k1*z);
    
     for ii=n:1:length(z);
        Ei(ii)=0;
    end
  %     subplot(2,2,2);
        hh1= plot(z,Ei ,                'MarkerFaceColor','r',...
                    'MarkerEdgeColor','k',    'Color','red',...
                        'LineWidth',3,...
                'MarkerSize',50); 
         %    set(hh1,'Marker','o') ;
        %     set(hh1,'MarkerSize',18);

        hold on;
        plot(zL,yL);
        plot(y1,z1,                'MarkerFaceColor','r',...
                    'MarkerEdgeColor','k',    'Color','blue',...
                        'LineWidth',5,...
                'MarkerSize',50); 
        hold off;

        if n>=length(z)+1;
     %   subplot(2,2,2);
     hold on;
        if jr<100;
            jr=jr+1;
            for jj=1:length(z)-jr;
                Er(jj)=0;
            end
        end
       hold on;
        hh2=plot(z,Er,                'MarkerFaceColor','r',...
                    'MarkerEdgeColor','k',    'Color','blue',...
                        'LineWidth',2,...
                'MarkerSize',50); 
        hold on;
        plot(z,Er+Ei,                'MarkerFaceColor','r',...
                    'MarkerEdgeColor','k',    'Color','black',...
                        'LineWidth',5,...
                'MarkerSize',50); 
        hold on;
        Ett=trans*e0*cos(2*pi*n/framemax-k2*zt);
            for ii=n+1-length(zt):1:length(zt);
                Ett(ii)=0;
            end

                plot(zt,Ett,                'MarkerFaceColor','r',...
                    'MarkerEdgeColor','k',    'Color','green',...
                        'LineWidth',4,...
                'MarkerSize',50); 
hold off;
    end
    
 
    axis([zmin ztmax -2*e0 2*e0]);axis on;
    title('红色为正向行波，蓝色为负向行波（反射波），黑色为形成的驻波；绿色为透射波。');
    if n==251;
        grid on;
        pause;
    end
    grid on;
     getframe;

end
%movie(m1,3);