拟合求助：1stopt正版sharedmodel计算求助，我自己写好代码，帮求4组 - 计算模拟 - 理论与计算方法

版块导航: 正在加载中...

客户端APP下载

论文辅导

申博辅导

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

962321702

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 934.9
帖子: 270
在线: 79.3小时
虫号: 1233166
注册: 2011-03-14
专业: 化学反应工程

送红花一朵

引用回帖:

3楼: Originally posted by 独孤神宇 at 2018-10-31 15:33:06
模型公式: a = (p/3251.622+4.36826E-7*p/(3251.622*t)+4.36826E-7*4.63956E-11*p/(3251.622*t^2)+4.36826E-7*3.285*(p/3251.622)/(k12*k13*t+4.36826E-7*p/3251.622+k13*t^2))/3.285
a = (p/3251.622+4 ...

我之前公式弄错了，你能再帮我算一下吗？
参数没有约束
第一组数据：
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

赞一下

回复此楼

11楼2020-02-24 20:38:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 962321702 的主题更新

返回列表