Efield_xy_plane_cstmatlab_CST超表面_超表面_适用于聚焦超表面_cstmatlab

clear all; em=importdata('orbital_1.txt'); x=em.data(:,1)/1000; % unit:um y=em.data(:,2)/1000; z=em.data(:,3)/1000; % unit:V/m ex_re=em.data(:,4); ey_re=em.data(:,5); ez_re=em.data(:,6); ex_im=em.data(:,7); ey_im=em.data(:,8); ez_im=em.data(:,9); ex=ex_re+1i*ex_im; ey=ey_re+1i*ey_im; % x,y array y1(1)=y(1); x1(1)=x(1); u=1; v=1; ex1(1,1)=ex(1); ey1(1,1)=ey(1); for m=2:1:size(x,1) if x(m)~=x(m-1) v=v+1; u=1; x1(v)=x(m); ex1(u,v)=ex(m); ey1(u,v)=ey(m); elseif x(m)==x(m-1) u=u+1; y1(u)=y(m); ex1(u,v)=ex(m); ey1(u,v)=ey(m); end end % amplitude and phase of circular polarization eR=ex1+1i*ey1; eL=ex1-1i*ey1; AeR=abs(eR); mR=max(max(AeR)); AeR=AeR/mR; AeL=abs(eL); mL=max(max(AeL)); AeL=AeL/mL; PeR=angle(eR)/pi*180+180; PeL=angle(eL)/pi*180+180; AeX=abs(ex1); %PeX=angle(ex1)/pi*180+180; PeX=angle(ex1); fig3=figure(1); %h3=pcolor(x1,y1,AeX); h3=pcolor(x1,y1,AeX.*sin(PeX)); %pcolor绘制伪彩色图 caxis([-20000 20000]); xlabel('\itx\rm (\mum)','FontName','Arial','FontSize',18); ylabel('\ity\rm (\mum)','FontName','Arial','FontSize',18); shading interp; colormap jet; axis([-1.2 1.2 -1.2 1.2]); %坐标轴刻度范围 c=colorbar; c.Label.String = 'Normalized \itE \rm(a.u.)'; c.Label.FontName='Arial'; c.Label.FontSize=18; set(gca,'FontName','Arial','FontSize',18); set(gca,'xtick',-7:2:7,'ytick',-7:2:7); axis square; print(fig3,'OAM RCP to LCP','-dtiff','-r600'); % fig4=figure(2); % h4=pcolor(x1,y1,PeR); % caxis([0 360]); % xlabel('\itx\rm (\mum)','FontName','Arial','FontSize',18); % ylabel('\ity\rm (\mum)','FontName','Arial','FontSize',18); % shading interp; % colormap jet; % colorbar('Ticks',[0,90,180,270,360],'TickLabels',{'0','\pi/2','\pi','3\pi/2','2\pi'},'FontName','Arial','FontSize',16); % set(gca,'FontName','Arial','FontSize',18); % % set(gca,'xtick',-7:2:7,'ytick',-7:2:7); % axis square; % print(fig4,'OAM REP (1,3) to LEP (1,3) Phase','-dtiff','-r600'); % % % coordinate % z1(1)=z(1); m=2; % for n=2:1:size(z,1) % if z(n) ~= z(n-1) % z1(m)=z(n); % m=m+1; % end % end % % pick z position. z>0 transmitted wave % % q=round(size(z1,2)/2)-36; % for q=1 % plot multiple z-plane % for n=1:1:size(z,1) % if z(n)==z1(q) % break % end % end % for m=n+1:1:size(z,1) % if z(m)==z(n) % if y(m)~=y(m-1) % v=1; % u=u+1; % y1(u)=y(m); % ex1(v,u)=ex(m); ey1(v,u)=ey(m); % elseif y(m)==y(m-1) % v=v+1; % x1(v)=x(m); % ex1(v,u)=ex(m); ey1(v,u)=ey(m); % end % else % break % end % end % eR=ex1-1i*ey1; % eL=ex1+1i*ey1; % % amplitude and phase of circular polarization % AeR=abs(eR); % AeL=abs(eL); % PeR=angle(eR)/pi*180; % PeL=angle(eL)/pi*180; % end