conductor_perpendicular.m

电磁场圆极化

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

%

clear all;
zmin=-10;
zmax=0;
delz = 0.1;
e0=1;k1=1;
framemax = 1000;
z = [zmin:delz:zmax];
y=0*z;
z1=-1:.1:1;
y1=0*z1;
jjj=1;

for n=1:framemax;
    Ei=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(z,y);
        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;
        Er=-e0*cos(2*pi*n/framemax+k1*z);
        if jjj<100;
            jjj=jjj+1;
            for jj=1:101-jjj;
                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 off;
    end
    
    if n>208;
S=e0^2*sin(2*k1*z)*sin(4*pi*n/framemax);
hold on;
plot(z,abs(S),                'MarkerFaceColor','r',...
                    'MarkerEdgeColor','k',    'Color','green',...
                        'LineWidth',5,...
                'MarkerSize',50); 
            hold off;
        end
        
    axis([zmin 0 -2*e0 2*e0]);axis on;
    title('红色为正向行波，蓝色为负向行波（反射波），黑色为形成的驻波。；绿色为驻波能量密度');
    if n==251;
        grid on;
        pause;
    end
    grid on;
     getframe;

end
%movie(m1,3);