fig511.m

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

% Na+, Cl- adsorption(desorption) (Vmol/g geotite),  NaCl-HCl-H2O sys at 29C
% Figure 5.11 in 2D plot

% Bar chart for NaCl-HCl-H2O 1,2,4,7,10,15,20,21# 2000.08.13
tic % start time recording
x=[0 0.5 1 2 3 4 4.5 5]; % conmmon x
yna=[-8.62 -9.92 -11.21 -13.81 -16.41 -19.01 -20.31 -21.6]; % delta Na+
ycl=[-53.9 -57.24 -60.58 -67.25 -73.93 -80.6 -83.94 -87.27]; % delta Cl- 

yna=-yna;
ycl=-ycl;  % adsorption into desorption just for an normal figure!


Y=[yna' ycl']; % queue is needed

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]);

w1=bar(x,Y,2,'group'); % figure 1 : linewidth=2, style='group'
legend(w1,'Na^{+}','Cl^{-}') % label legend
ylabel('10^{-6} mol')
title('Na^{+},Cl^{ -} release,  NaCl-HCl-H2O sys, ')
set(gca,'xtick',[-1 0 0.5 1 2 3 4 4.5 5])
set(gca,'xticklabel',' |1#|2#|4#|7#|10#|15#|20#|21#| ')
hold on
 xi=0:0.1:5; % interpolation and plot a curve
   H1=interp1(x,Y(:,1),xi,'linear');
   H2=interp1(x,Y(:,2),xi,'linear'); % Na,Cl lines cause large value
plot(xi,H1,'k') % Na line
   plot(xi,H2,'k') % Cl line
plot([-1 6],[0 0],'g') % the line : y=0
text(2.5,75,'Cl^{-}')
text(2.5,20,'Na^{+}')

r=get(fh,'papersize') % for printing


hold on
plot([0 6],[0 0],'g') % the line : y=0

disp('pls shift the legend to an andiquate position then press anykey')
pause

str={'图5.11 HCl-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,'fig511.jpeg') % succeed in 2000.10.01
toc; % stop time recording