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

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

962321702

银虫 (小有名气)

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

[求助] 拟合求助：1stopt正版sharedmodel计算求助，我自己写好代码，帮求4组已有1人参与

代码：
parameters  k12,k13;
constant k9=9.92471*10^(-15), k10=4.36826*10^(-7), k11=4.63956*10^(-11), h=3251.622, m=3.285;
variable a,p;
parvariable t=[0.000001,0.000000000001];
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

对应其他三组数据：
k9       2.86026*10^(-14)    9.38162*10^(-14)    2.49148*10^(-13)
k10       5.01604*10^(-7)       5.23156*10^(-7)       4.79116*10^(-7)
k11       5.99747*10^(-11)    7.23603*10^(-11)    7.57401*10^(-11)
h          4078.309                   5344.079                   6791.551
m       3.285                         3.285                         3.285
a,p
//0.598 34.5
0.612 46.5
0.636 60
0.66 78

//0.538 34.5
0.562 46.5
0.583 60
0.598 78

//0.465 34.5
0.495 46.5
0.517 60
0.533 78@月只蓝@beefly

回复此楼

» 猜你喜欢

【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有283人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

1楼2018-10-31 14:02:16

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

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

赞一下

回复此楼

10楼2020-02-24 20:28:39

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

查看全部 11 个回答

1028942887

顶一下

发自小木虫Android客户端

2楼2018-10-31 15:25:27

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

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
962321702: 金币+200, ★★★★★最佳答案 2018-11-01 22:34:27

模型公式: 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.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)/(238.258176720815*3.69008634786073e-5*t+4.36826E-7*p/3251.622+3.69008634786073e-5*t^2))/3.285
模型公式: t = (1+k13*3.285*t/(k12*k13*t+4.36826E-7*p/3251.622+k13*t^2))*(4.36826E-7*p/(3251.622*t)+4.36826E-7*3.285*(p/3251.622)/(k12*k13*t+4.36826E-7*p/3251.622+k13*t^2)+2*4.36826E-7*4.63956E-11*p/(3251.622*t^2)+9.92471E-15/t)
      t = (1+3.69008634786073e-5*3.285*t/(238.258176720815*3.69008634786073e-5*t+4.36826E-7*p/3251.622+3.69008634786073e-5*t^2))*(4.36826E-7*p/(3251.622*t)+4.36826E-7*3.285*(p/3251.622)/(238.258176720815*3.69008634786073e-5*t+4.36826E-7*p/3251.622+3.69008634786073e-5*t^2)+2*4.36826E-7*4.63956E-11*p/(3251.622*t^2)+9.92471E-15/t)

迭代数: 24
计算用时(时:分:秒:微秒): 00:00:01:647
优化算法: 通用全局优化算法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0
残差平方和(SSR): 0
相关系数(R): 1
相关系数之平方(R^2): 1
修正R平方(Adj. R^2): -100
确定系数(DC): 1

参数       最佳估算
-------------------- -------------
k12 238.258176720815
k13 3.69008634786073E-5
t0 3.09325866287863E-7
t1 3.70106253183998E-7
t2 4.45708163439644E-7
t3 5.1208433542262E-7

====== 结果输出 ======

文件: 数据文件-1
No 目标 a 计算 a
1 0.638 0.638
2 0.667 0.667
3 0.684 0.684
4 0.713 0.713

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

962321702

数值计算

3楼2018-10-31 15:33:06

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

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

模型公式: a = (p/4078.309+5.01604E-7*p/(4078.309*t)+5.01604E-7*5.99747E-11*p/(4078.309*t^2)+5.01604E-7*3.285*(p/4078.309)/(k12*k13*t+5.01604E-7*p/4078.309+k13*t^2))/3.285
      a = (p/4078.309+5.01604E-7*p/(4078.309*t)+5.01604E-7*5.99747E-11*p/(4078.309*t^2)+5.01604E-7*3.285*(p/4078.309)/(3.55859929088445e-10*6027.61803316927*t+5.01604E-7*p/4078.309+6027.61803316927*t^2))/3.285
模型公式: t = (1+k13*3.285*t/(k12*k13*t+5.01604E-7*p/4078.309+k13*t^2))*(5.01604E-7*p/(4078.309*t)+5.01604E-7*3.285*(p/4078.309)/(k12*k13*t+5.01604E-7*p/4078.309+k13*t^2)+2*5.01604E-7*5.99747E-11*p/(4078.309*t^2)+2.86026E-14/t)
      t = (1+6027.61803316927*3.285*t/(3.55859929088445e-10*6027.61803316927*t+5.01604E-7*p/4078.309+6027.61803316927*t^2))*(5.01604E-7*p/(4078.309*t)+5.01604E-7*3.285*(p/4078.309)/(3.55859929088445e-10*6027.61803316927*t+5.01604E-7*p/4078.309+6027.61803316927*t^2)+2*5.01604E-7*5.99747E-11*p/(4078.309*t^2)+2.86026E-14/t)

迭代数: 24
计算用时(时:分:秒:微秒): 00:00:01:616
优化算法: 通用全局优化算法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0
残差平方和(SSR): 0
相关系数(R): 1
相关系数之平方(R^2): 1
修正R平方(Adj. R^2): -100
确定系数(DC): 1

参数       最佳估算
-------------------- -------------
k12 3.55859929088445E-10
k13 6027.61803316927
t0 6.94110112812778E-7
t1 7.847386865608E-7
t2 8.49821236885364E-7
t3 9.23503826465456E-7

====== 结果输出 ======

文件: 数据文件-1
No 目标 a 计算 a
1 0.598 0.598
2 0.612 0.612
3 0.636 0.636
4 0.66 0.66

赞一下

回复此楼

数值计算

4楼2018-10-31 15:33:36

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

查看全部 11 个回答