fig701.m

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

% 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
yh=[61.05 54.95 48.84 36.63 24.42 12.21 6.11 0.005] ;% delta H+
yna=[-8.62 -9.92 -11.21 -13.81 -16.41 -19.01 -20.31 -21.6]; % delta Na+
yfe=[-0.0213 -0.0215 -0.0216 -0.022 -0.0224 -0.0227 -0.0229 -0.0231]; % delta Fe3+
ycl=[-53.9 -57.24 -60.58 -67.25 -73.93 -80.6 -83.94 -87.27]; % delta Cl- 

Y=[yh' yna' yfe' ycl']; % 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]);
h2=bar3(x,Y);
legend(h2,'H^{+}','Na^{+}','Fe^{3+}','Cl^{-}') % label legend
zlabel('10^{-6} mol')
title('NaCl-HCl-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| ')

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

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