fig609.m

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

% NaCl-NaH2PO4-H2O sys summary 2000 08 04 23:15
% ------------------------------------------------------------------------
% BALANCED soln, X10^6M
% transparency  ! relative unit percentage
z1=[6.8 7.8 14.4 8.6 17.3 40.2 17.3 34.4 88.3 28.5];
z2=[32.6 52.7 90.4 94.0 64.0 91.5 94 95 96 91.7]; 
z3=[100 99.8 99.5 99 98.5 98 97]; 
zmud=[z1 z2 z3];
z=100-zmud;

fh=figure('Color','w','position',[200 100 500 340]);
h1=axes('position',[0 0 1 1],'visible','off');
h2=axes('position',[0.07 0.12 0.32 0.78]);
contour3draw27(z*1.48,'absorptance, %, Equ sol')
% ---------------------------------------------------------------------
% P reading 27

z1=[67.3 71.9 72.4 75.6 76.0 76.3 76.9 81.1 81.5 81.8];
z2=[85.4 89.5 88.8 93.5 87.8 90.1 89.4 95.2 95.0 89.8];
z3=[99.2 99.7 99.6 99.5 99.7 99.5 99.4];
z=[z1 z2 z3];

y=-6.4877+84.352*exp(-z/37.674); % ppm P2O5
y=y*2.5;
y(1)=y(1)*4;
y(2)=y(2)*4;
y(3)=y(3)*4;
y(4)=y(4)*4;
y(5)=y(5)*4;
y(6)=y(6)*4;
y(7)=y(7)*4;
y(8)=y(8)*4;
y(10)=y(10)*4;
y(11)=y(11)*4;
y(12)=y(12)*4;
y(15)=y(15)*4;

% transform ppm(P2O5) to VM (H2PO4-)
for i=1:27
   y(i)=y(i)*1.42E-5*1E6;
end


for i=1:27 % to avoid illegal data
   if y(i)<0
      y(i)=0;
   end
end

h3=axes('position',[0.55 0.12 0.32 0.78]);
contour3draw27(y/8,'[H2PO4^{-}], Equ sol')

str={'图6.9 平衡液的吸光度与平衡液中H2PO4^{ - }浓度的关联'}; %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,'fig609.jpeg') % save as a jpeg file

% 2000.10.01