|
|
 送红花一朵
我之前公式弄错了，你能再帮我算一下吗？
参数没有约束
第一组数据：
Parameters k12,k13; //拟合有机胺质子化常数，氨基甲酸盐稳定常数
Constant k9=9.92471E-15, k10=4.36826E-07, k11=4.63956E-11, H=3251.622, m=3.285; //离子积，CO2水解，HCO3水解，亨利系数，MEA初始浓度
Variable a,p; //负载，分压为变量
ParVariable t[0.000000000001,0.000001]; //氢离子浓度为参数变量
SharedModel;
Function a=(p/H+k10*p/(H*t)+k10*k11*p/(H*t^2)+k10*m*(p/H)/(k12*k13*t+k10*p/H+k13*t^2))/m;
t*(1+k13*m*t/(k12*k13*t+k10*p/H+k13*t^2))=(k10*p/(H*t)+k10*m*(p/H)/(k12*k13*t+k10*p/H+k13*t^2)+2*k10*k11*p/(H*t^2)+k9/t);
Data;
0.638 34.5
0.667 46.5
0.684 60
0.713 78
第二组数据：
Parameters k12,k13; //拟合有机胺质子化常数，氨基甲酸盐稳定常数
Constant k9=2.86026E-14, k10=5.01604E-07, k11=5.99747E-11, H=4078.309, m=3.285; //离子积，CO2水解，HCO3水解，亨利系数，MEA初始浓度
Variable a,p; //负载，分压为变量
ParVariable t[0.000000000001,0.000001]; //氢离子浓度为参数变量
SharedModel;
Function a=(p/H+k10*p/(H*t)+k10*k11*p/(H*t^2)+k10*m*(p/H)/(k12*k13*t+k10*p/H+k13*t^2))/m;
t*(1+k13*m*t/(k12*k13*t+k10*p/H+k13*t^2))=(k10*p/(H*t)+k10*m*(p/H)/(k12*k13*t+k10*p/H+k13*t^2)+2*k10*k11*p/(H*t^2)+k9/t);
Data;
0.598 34.5
0.612 46.5
0.636 60
0.66 78
第三组数据：
Parameters k12,k13; //拟合有机胺质子化常数，氨基甲酸盐稳定常数
Constant k9=9.38162E-14, k10=5.23156E-07, k11=7.23603E-11, H=5344.079, m=3.285; //离子积，CO2水解，HCO3水解，亨利系数，MEA初始浓度
Variable a,p; //负载，分压为变量
ParVariable t[0.000000000001,0.000001]; //氢离子浓度为参数变量
SharedModel;
Function a=(p/H+k10*p/(H*t)+k10*k11*p/(H*t^2)+k10*m*(p/H)/(k12*k13*t+k10*p/H+k13*t^2))/m;
t*(1+k13*m*t/(k12*k13*t+k10*p/H+k13*t^2))=(k10*p/(H*t)+k10*m*(p/H)/(k12*k13*t+k10*p/H+k13*t^2)+2*k10*k11*p/(H*t^2)+k9/t);
Data;
0.538 34.5
0.562 46.5
0.583 60
0.598 78
第四组数据：
Parameters k12,k13; //拟合有机胺质子化常数，氨基甲酸盐稳定常数
Constant k9=2.49148E-13, k10=4.79116E-07, k11=7.57401E-11, H=6791.551, m=3.285; //离子积，CO2水解，HCO3水解，亨利系数，MEA初始浓度
Variable a,p; //负载，分压为变量
ParVariable t[0.000000000001,0.000001]; //氢离子浓度为参数变量
SharedModel;
Function a=(p/H+k10*p/(H*t)+k10*k11*p/(H*t^2)+k10*m*(p/H)/(k12*k13*t+k10*p/H+k13*t^2))/m;
t*(1+k13*m*t/(k12*k13*t+k10*p/H+k13*t^2))=(k10*p/(H*t)+k10*m*(p/H)/(k12*k13*t+k10*p/H+k13*t^2)+2*k10*k11*p/(H*t^2)+k9/t);
Data;
0.465 34.5
0.495 46.5
0.517 60
0.533 78 |
|
回复此楼