fig801.m

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

% NaCl-NaH2PO4-H2O sys summary 2000 08 04 23:15

% ------------------------------------------------------------------------
% ORIGINAL soln, X10^6M
% H+
z1=[6.729 6.165 6.067 5.675 5.584 5.470 4.786 4.632 5.000 3.881];
z2=[3.837 3.784 3.772 4.045 2.421 2.328 2.244 2.517 2.437 1.496];
z3=[0.679 0.495 0.448 0.448 0.495 0.476 0.411];
z=[z1 z2 z3];
zh1=z*1E-6;


% ------------------------------------------------------------------------
% ORIGINAL soln, X10^6M

% H2PO4-
z1=[1500 1350 1350 1200 1200 1200 900 900 900 600];
z2=[600 600 600 600 300 300 300 300 300 150];
z3=[0 0 0 0 0 0 0];
z=[z1 z2 z3];
zp1=z;
% ------------------------------------------------------------------------
% ************************************************************************
% ------------------------------------------------------------------------
% equ. soln, X10^6M

% H2PO4-
z1=[1085.8 855.2 831.7 689 672 659.3 634.4 470.3 113.9 440.5];
z2=[320.2 192.2 53.3 20 243.5 43.6 48.8 8.9 10.2 119.8];
z3=[0 0 0 0 0 0 0];
z=[z1 z2 z3];
zp2=z;

% ------------------------------------------------------------------------
%  Ion adsorption mol per g of geothite :  Vmol/g

zp=(zp1-zp2)*0.2;

fh=figure('Color','w','position',[200 100 500 350]);
h1=axes('position',[0 0 1 1],'visible','off');
h2=axes('position',[0.07 0.12 0.32 0.8]);
draw27801(zh1,'H^{+} Ini, M')

h3=axes('position',[0.55 0.12 0.32 0.8]);
draw27801(zp,'P adsptn, 10^{-6}mol/g')
set(h3,'Zlim',[0 200])

str={'图8.1  原始溶液中H^{ + }浓度与磷吸附量的关系'}; %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,'fig801.jpeg') % succeed in 2000.11.09