Variable p,c;
Function P=(c-P/H)*(10^(A*(c-P/H)+B));
Data;
//p        c
0        0
0.001        0.307598804
0.002        0.538975997
0.003        0.777628765
0.004        0.959350561
0.005        1.069744162
0.006        1.190592446
0.008        1.340828541
0.011        1.512368753
0.014        1.643259589
0.019        1.769609554
0.023        1.903576465
0.03        2.027165658
0.037        2.151392907
0.048        2.294367083
0.06        2.425558327
0.08        2.598588163
0.1        2.723808495
0.136        2.912380082
0.17        3.049909679
0.217        3.206076964
0.279        3.386393547
0.388        3.644346885
0.444        3.758352826

clear,clc
p=[0  0.001  0.002  0.003  0.004  0.005  0.006  0.008  0.011  0.014  0.019  0.023  0.03  0.037  0.048  0.06  0.08  0.1  0.136  0.17  0.217  0.279  0.388  0.444];
c=[0  0.307598804  0.538975997  0.777628765  0.959350561  1.069744162  1.190592446  1.340828541  1.512368753  1.643259589  1.769609554  1.903576465  2.027165658  2.151392907  2.294367083  2.425558327 2.598588163  2.723808495  2.912380082  3.049909679  3.206076964 3.386393547  3.644346885  3.758352826];
cc=c';
pp=p';
F=@(q,x)(x(:,1)-x(:,2)/q(1)).*(10.^(q(2)*(x(:,1)-x(:,2)/q(1)))+q(3))-x(:,2);   
x=[cc pp];

q0=[1    1     1];
warning off
q=nlinfit(x,zeros(size(x,1),1),F,q0);
disp('H A B分别为')
disp(num2str(q));
plot(x(:,1),x(:,2),'r^');hold on;
ezplot(@(x,y)F(q,[x,y]),[0,3.8,-0.1,0.5]);