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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (5722)
>考研 (2572)
>导师招生 (1877)
>虫友互识 (547)
>文献求助 (241)
>硕博家园 (171)
>休闲灌水 (145)
>招聘信息布告栏 (124)
>考博 (119)
>论文投稿 (48)
>博后之家 (44)
>基金申请 (38)
>教师之家 (36)
>公派出国 (27)
>海外博后 (23)
>绿色求助(高悬赏) (22)

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

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

回复此楼

» 猜你喜欢

［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有81人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有2人回复
070200求调剂，一志愿某211，288分已经有18人回复
法国博士后职位已经有0人回复
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。已经有0人回复
广州大学光电信息工程专业调剂，招收物理学专业学生已经有1人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有0人回复

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

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

1028942887

顶一下

发自小木虫Android客户端

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

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

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31016.3
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 856.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
金币: 31016.3
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 856.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的回帖

独孤神宇

版主 (知名作家)

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

【答案】应助回帖

模型公式: a = (p/5344.079+5.23156E-7*p/(5344.079*t)+5.23156E-7*7.23603E-11*p/(5344.079*t^2)+5.23156E-7*3.285*(p/5344.079)/(k12*k13*t+5.23156E-7*p/5344.079+k13*t^2))/3.285
      a = (p/5344.079+5.23156E-7*p/(5344.079*t)+5.23156E-7*7.23603E-11*p/(5344.079*t^2)+5.23156E-7*3.285*(p/5344.079)/((-30.8828030721688)*(-4289195.3669547)*t+5.23156E-7*p/5344.079+(-4289195.3669547)*t^2))/3.285
模型公式: t = (1+k13*3.285*t/(k12*k13*t+5.23156E-7*p/5344.079+k13*t^2))*(5.23156E-7*p/(5344.079*t)+5.23156E-7*3.285*(p/5344.079)/(k12*k13*t+5.23156E-7*p/5344.079+k13*t^2)+2*5.23156E-7*7.23603E-11*p/(5344.079*t^2)+9.38162E-14/t)
      t = (1+(-4289195.3669547)*3.285*t/((-30.8828030721688)*(-4289195.3669547)*t+5.23156E-7*p/5344.079+(-4289195.3669547)*t^2))*(5.23156E-7*p/(5344.079*t)+5.23156E-7*3.285*(p/5344.079)/((-30.8828030721688)*(-4289195.3669547)*t+5.23156E-7*p/5344.079+(-4289195.3669547)*t^2)+2*5.23156E-7*7.23603E-11*p/(5344.079*t^2)+9.38162E-14/t)

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

参数       最佳估算
-------------------- -------------
k12 -30.8828030721688
k13 -4289195.3669547
t0 1.98782097449534E-9
t1 2.54773452408826E-9
t2 3.15576671056898E-9
t3 3.987183582274E-9

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

文件: 数据文件-1
No 目标 a 计算 a
1 0.538 0.538
2 0.562 0.562
3 0.583 0.583
4 0.598 0.598

赞一下

回复此楼

数值计算

5楼2018-10-31 15:34:15

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

独孤神宇

版主 (知名作家)

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

【答案】应助回帖

模型公式: a = (p/6791.551+4.79116E-7*p/(6791.551*t)+4.79116E-7*7.57401E-11*p/(6791.551*t^2)+4.79116E-7*3.285*(p/6791.551)/(k12*k13*t+4.79116E-7*p/6791.551+k13*t^2))/3.285
      a = (p/6791.551+4.79116E-7*p/(6791.551*t)+4.79116E-7*7.57401E-11*p/(6791.551*t^2)+4.79116E-7*3.285*(p/6791.551)/(0*5009.21978072262*t+4.79116E-7*p/6791.551+5009.21978072262*t^2))/3.285
模型公式: t = (1+k13*3.285*t/(k12*k13*t+4.79116E-7*p/6791.551+k13*t^2))*(4.79116E-7*p/(6791.551*t)+4.79116E-7*3.285*(p/6791.551)/(k12*k13*t+4.79116E-7*p/6791.551+k13*t^2)+2*4.79116E-7*7.57401E-11*p/(6791.551*t^2)+2.49148E-13/t)
      t = (1+5009.21978072262*3.285*t/(0*5009.21978072262*t+4.79116E-7*p/6791.551+5009.21978072262*t^2))*(4.79116E-7*p/(6791.551*t)+4.79116E-7*3.285*(p/6791.551)/(0*5009.21978072262*t+4.79116E-7*p/6791.551+5009.21978072262*t^2)+2*4.79116E-7*7.57401E-11*p/(6791.551*t^2)+2.49148E-13/t)

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

参数       最佳估算
-------------------- -------------
k12 0
k13 5009.21978072262
t0 7.51487920320897E-7
t1 8.22783417437169E-7
t2 8.95866448673935E-7
t3 9.913208023056E-7

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

文件: 数据文件-1
No 目标 a 计算 a
1 0.465 0.465
2 0.495 0.495
3 0.517 0.517
4 0.533 0.533

赞一下

回复此楼

数值计算

6楼2018-10-31 15:34:33

已阅回复此楼关注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[9.25236E-11,1.13084E-10],k13[2.21332E+08,2.70517E+08];  //拟合有机胺质子化常数，氨基甲酸盐稳定常数
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[2.37085E-10,2.89770E-10],k13[1.41991E+08,1.73545E+08];  //拟合有机胺质子化常数，氨基甲酸盐稳定常数
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[7.28625E-10,8.90542E-10],k13[8.36054E+07,1.02184E+08];  //拟合有机胺质子化常数，氨基甲酸盐稳定常数
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[1.97175E-09,2.40992E-09],k13[5.22724E+07,6.38885E+07];  //拟合有机胺质子化常数，氨基甲酸盐稳定常数
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

赞一下

回复此楼

7楼2018-11-01 20:50:30

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

独孤神宇

版主 (知名作家)

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

引用回帖:

7楼: Originally posted by 962321702 at 2018-11-01 20:50:30
我之前参数没有给范围，所有拟合的结果不好，能重新帮我算一下吗？
第一组数据：
Parameters k12,k13; //拟合有机胺质子化常数，氨基甲酸盐稳定常数
Constant k9=9.92471E-15, k10=4.36826E-07, k11=4.63956E ...

已经发给你了。

回复此楼

数值计算

8楼2018-11-01 20:56:57

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

新星之光1

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 165.4
帖子: 67
在线: 6.8小时
虫号: 10731884
注册: 2018-10-17
性别: GG
专业: 遥感机理与方法

大佬，你的1stopt安装包发我一份，网上好多没有注册

发自小木虫Android客户端

赞一下

回复此楼

9楼2018-11-07 19:08:30

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

962321702

银虫 (小有名气)

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

送红花一朵

引用回帖:

之前公式弄错了，能帮我再算一下吗？
第一组数据：
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的回帖

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

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 0703化学 +20	妮妮ninicgb 2026-04-04	23/1150	2026-04-06 23:06 by chenzhimin
[考研] 求调剂 +5	熊二想上岸 2026-04-06	5/250	2026-04-06 22:27 by chenzhimin
[考研] 285求调剂 +5	mapmath 2026-04-06	6/300	2026-04-06 17:18 by 蓝云思雨
[考研] 一志愿南航，数一英一学硕317求调剂！！ +6	Acaciad 2026-04-04	6/300	2026-04-06 12:13 by 考研学校招点人
[考研] 319分085702安全工程求调剂 +4	rious 2026-04-05	4/200	2026-04-06 11:08 by 静静静静静静静�
[考研] 0817化学工程与技术求调剂，一志愿中海洋319 +14	lv945 2026-04-04	14/700	2026-04-06 10:20 by 蓝云思雨
[考研] 272分求调剂 +4	wangyile2233 2026-04-02	4/200	2026-04-05 22:21 by 286640313
[考研] 262求调剂 +7	天下第一文 2026-04-04	8/400	2026-04-05 21:31 by 激流勇渡
[考研] 求调剂求调剂 +8	121. 2026-04-02	8/400	2026-04-05 20:15 by lys0704
[考研] 求调剂 +10	Hll胡 2026-04-04	10/500	2026-04-05 20:09 by nepu_uu
[考研] 材料专硕(0856) 339分求调剂 +10	哈哈哈鹅哈哈哈 2026-04-04	10/500	2026-04-05 18:51 by 蓝云思雨
[考研] 307分材料专业求调剂 +7	Hll胡 2026-04-05	7/350	2026-04-05 18:47 by 无际的草原
[考研] 313求调剂 +3	海日海日 2026-04-04	3/150	2026-04-05 07:48 by 544594351
[考研] 一志愿华北电力大学（北京），材料科学与工程学硕265，求调剂 +11	yelck 2026-04-03	12/600	2026-04-04 19:52 by dongzh2009
[考研] 309求调剂 +4	快乐的小白鸽 2026-04-04	5/250	2026-04-04 15:55 by cql1109
[考研] 387求调剂 +4	爱吃片豆土 2026-04-03	5/250	2026-04-04 08:10 by 岸上的一条鱼
[考研] 兽医调剂 +3	wh119216 2026-04-02	3/150	2026-04-03 19:34 by zrongyan
[硕博家园] 求老师收留 +9	lllq123 2026-04-03	9/450	2026-04-03 13:48 by 呼吸都是减肥
[考研] 一志愿华南师范大学-22408计算机-292分-求华南师范大学调剂 +4	爱读书的小鳄鱼 2026-04-02	4/200	2026-04-02 18:35 by 求调剂zz
[考研] 江苏苏北高校诚邀调剂同学 +3	zzll406 2026-03-31	3/150	2026-03-31 16:54 by 及时行乐fan

24小时热门版块排行榜

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

» 猜你喜欢

【答案】应助回帖

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

【答案】应助回帖

【答案】应助回帖

【答案】应助回帖

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