fig606607.m

来自「这是一个matlab应用的论文和源代码」· M 代码 · 共 67 行

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

% ------------------------------------------------------------------------
% 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;

% ------------------------------------------------------------------------
% BALANCED 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;

fh1=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]);
contourdraw27(zp1,'Ini [H2PO4^{-}], 10^{-6}M')

h3=axes('position',[0.55 0.12 0.32 0.8]);
contourdraw27(zp2,'Equ [H2PO4^{-}], 10^{-6}M')

str={'图6.6 吸附前后H2PO4^{ - }浓度变化（等高线图）'}; %the note of this picture! 
   set(gcf,'currentAxes',h1)
   text(0.5,0.03,str,'Fontsize',12,'horizontalAlignment','center')

[X,map]=capture(fh1);
 imwrite(X,map,'fig606.jpeg') % succeed in 2000.7.31


disp('the 3D picture will begin, pls press anykey')
pause

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

disp('n=-30 is good')

draw27(zp1,'Ini [H2PO4^{-}], 10^{-6}M')

h3=axes('position',[0.55 0.12 0.32 0.8]);
draw27(zp2,'Equ [H2PO4^{-}], 10^{-6}M')
set(h3,'Zlim',[0 1200])

str={'图6.7 吸附前后H2PO4^{ - }浓度变化（三维网格图）'}; %the note of this picture! 
   set(gcf,'currentAxes',h1)
   text(0.5,0.03,str,'Fontsize',12,'horizontalAlignment','center')

[X,map]=capture(fh2);
 imwrite(X,map,'fig607.jpeg') % succeed in 2000.10.01

% n=-30 best