fig702.m

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

% Bar chart for NaCl-NaH2PO4-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];       % common x
yh=[1.332 1.222 1.126 0.951 0.772 0.478 0.291 0.123];   % delta H+
yna=[0 0 0 0 0 0 -16.5 -18.7];    % delta Na+
yfe=[-0.0191 -0.0197 -0.0200 -0.0199 -0.0208 -0.0210 -0.0212 -0.0215]; % delta Fe3+
ycl=[-113.6 -110.77 -107.94 -102.28 -96.62 -90.96 -88.13 -85.3];       % delta Cl-
yp=[82.84 98.86 102.2 53.12 31.9 11.3 6.04 0];           					% dalta H2PO4-

Y=[yh' yna' yfe' ycl' yp'];  % data preparation in a matrix, queue is needed

f2=figure('Color','w','position',[200 100 470 350]);
h1=axes('position',[0 0 1 1],'visible','off');
h2=axes('position',[0.15 0.12 0.8 0.8],'color','w');
w2=bar3(x,Y);
legend(w2,'H^{+}','Na^{+}','Fe^{3+}','Cl^{-}','H2PO4^{-}') % label legend
title('NaH2PO4 \rightarrow H2O sys,  bar3')

set(gca,'ytick',[0 0.5 1 2 3 4 4.5 5])
set(gca,'yticklabel','1|2|4|7|10|15| |21| ')

zlabel('10^{-6} mol')

disp('pls adjust the legend by your mouse then press anykey')
pause

str={'图7.2 NaH2PO4-H2O子体系中各离子的数量变化（方块图）'}; %the note of this picture! 
   set(gcf,'currentAxes',h1)
   text(0.5,0.03,str,'Fontsize',12,'horizontalAlignment','center')
[X,map]=capture(f2);
 imwrite(X,map,'fig702.jpeg') % succeed in 2000.10.10

toc; % end time recording