fig601.m

这是一个matlab应用的论文和源代码

% 27 points for NaCl-HCl-H2O system 00.07.18 19:02
x1=[50 45 55 40 50 60 30 50 70 20 30 40 60 80 10 30 50];
x2=[70 90 5 0 10 20 40 60 80 100];
y1=[86.60 77.94 77.94 69.28 69.28 69.28 51.96 51.96 51.96 34.64 34.64 34.64 34.64 34.64];
y2=[17.32 17.32 17.32 17.32 17.32 8.66 0 0 0 0 0 0 0];
x=[x1 x2];
y=[y1 y2];

fh=figure('Color','w','position',[200 100 470 350]);
h1=axes('position',[0 0 1 1],'visible','off');
h2=axes('position',[0.1 0.12 0.8 0.8]);

x0=[0 100 50 0];
y0=[0 0 86.6 0];
c=[0.4 0.8 0.4];  % color

plot(x0,y0); % 3 lines
patch(x0,y0,c); % fill in the triangle with color c
plot(x,y,'o','color','k','Markersize',14) % dotted 27 points

hold on

%    basic lines in blue
lx=100/5;
ly=86.6/5;

for n=1:5
  x0=[50-lx/2*n 50+lx/2*n];
  y0=[86.6-ly*n 86.6-ly*n];
    plot(x0,y0);
     axis equal
  hold on
end

for n=1:5
   x0=[100-10*n 100-20*n];
   y0=[ly*n 0];
   plot(x0,y0);

end

for n=1:5
   x0=[10*n 20*n];
   y0=[ly*n 0];
   plot(x0,y0);
end
axis off % clear the 2 axes
hold on % keep the picture

% label the 27 points
for i=1:27
   a=num2str(i); % transform number into string
   p=0.5;
   if i>9
      p=1.7;
      end
      text(x(i)-p,y(i),a) % label every point
   end
   
   text(-10,0,'H2O')
   text(103,0,'NaCl')
   text(43,91,'NaH2PO4')
   
str={'图6.1 NaCl-NaH2PO4-H2O吸附体系原始溶液各点组成'}; %the note of this picture! 
   set(gcf,'currentAxes',h1)
   text(0.5,0.03,str,'Fontsize',12,'horizontalAlignment','center')
   
   [X,map]=capture(fh);
imwrite(X,map,'fig601.jpeg') % succeed in 2000.10.01