clear all;
em=importdata('my-326-0.05a-rushe-3.1.txt');
x=em.data(:,1); % unit:um
y=em.data(:,2);
z=em.data(:,3);
% 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;
eR=ex+1i*ey;
eL=ex-1i*ey;
% x,y array
x1(1)=x(1); z1(1)=z(1); u=1; v=1;
eR1(1,1)=eR(1); eL1(1,1)=eL(1);
for m=2:1:size(x,1)
if x(m)~=x(m-1)
v=v+1;
u=1;
x1(v)=x(m);
eR1(u,v)=eR(m); eL1(u,v)=eL(m);
elseif x(m)==x(m-1)
u=u+1;
z1(u)=z(m);
eR1(u,v)=eR(m); eL1(u,v)=eL(m);
end
end
% amplitude and phase of circular polarization
AeR=abs(eR1);
AeL=abs(eL1);
PeR=angle(eR1)/pi*180+180;
PeR1=angle(eR1);
% mR=max(max(AeR));
% AeR=AeR/mR;
A=AeR.^2;
z2=z1';
x2=x1';
B=A';
fig1=figure(1); %figure(1)第一个图
%h3=pcolor(x1,z1,A);
%caxis([0 5]);
h3=pcolor(x1,z1,PeR);
caxis([0 360]);
%h3=pcolor(x1,z1,PeR1);
%caxis([0 pi]);
%h3=pcolor(x1,z1,AeR.*sin(PeR1));
%caxis([-3.5 3.5]);
xlabel('\itx\rm (um)','FontName','Times New Roman','FontSize',16); %x轴标识
ylabel('\itz\rm (um)','FontName','Times New Roman','FontSize',16); %y轴标识
shading interp; %对图形对象的颜色着色进行色彩的插值处理,使色彩平滑过渡 ;
colormap jet; %大值为红,小值为蓝色
c=colorbar; %定义色列参数
%c.Label.String = 'Amplitude';
c.Label.FontName='Times New Roman';
c.Label.FontSize=20;
set(gca,'FontName','Times New Roman','FontSize',20); %调整图形的坐标属性
%set(gca,'YDir','reverse')%对Y方向反转
set(gca,'tickdir','out')%刻度外侧
axis equal;
axis([-7.65 7.65 0 10]); %坐标轴刻度范围
set(gca, 'Layer', 'top')
print(fig1,'y Amp','-dpng','-r600'); %输出图像
评论0