|
|
参考程序如下，会用MATLAB的话，很简单
clear;clc;
inputK=[9.3 3.8 1.6];
inputK=10.^inputK;
ph=[0:0.01:14];
h=exp(-ph*log(10));
z=cumprod(inputK'*h);
z0=1./(1+sum(z));
z1=[];
for i=1:size(inputK,2)
z1(i, =z(i,.*z0;
end;
plot(ph,z0,ph,z1);hold on;xlabel('pH');grid on; |
|
回复此楼
