star.m

来自「Star plots - MATLAB files for Graphical 」· M 代码 · 共 48 行

function star(v,v0,i,lambda,flagl)
components=  {'Mg'	'Sr'	'Mn'	'Ba'	'Ti'	'Zr'	'Zn'	'V'};

set(0,'DefaultAxesFontName','times');
set(0,'DefaultTextFontName','times');
set(0,'DefaultAxesFontSize',8);
set(0,'DefaultTextFontSize',8);
%set(0,'linewidth',1);


hold off
%data norm

dv=v./v0;
Sv=sum(dv)
dv=dv;
I=sum(v)/sum(v0)*lambda;
%I=sum(v);

%I=(sin(pi/8)*8); 				%in this case r0=1
r0=sqrt(I/(sin(pi/8)*Sv));
a=r0.*dv';


r=[ones(size(a,1),1)*r0 a]';
r=r(:);
x=[0:pi/8:2*pi-pi/8]';
y=sin(x).*r;
x=cos(x).*r;
x(end+1)=x(1);
y(end+1)=y(1);

hold on
for i=1:8
   in=(i-1)*2+1;
   %   fill([0 x(in:in+2)' 0 ],[0 y(in:in+2)' 0 ],i)
       fill([0 x(in:in+2)' 0 ],[0 y(in:in+2)' 0 ],'w','linewidth',1.4)
   %     text(x(in+1)/r(in+1)*r(in)*2,y(in+1)/r(in+1)*r(in)*2,components{i})
   if flagl
      %   text(x(in+1)/r(in+1)*2.5,y(in+1)/r(in+1)*2.5,components{i})
         text(x(in+1)/r(in+1)*2.5,y(in+1)/r(in+1)*2.5,components{i})

   end
end
hold off