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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求高版本的1stopt代跑一趟。谢谢
51/1返回列表
查看: 1380  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhaoshazhu

新虫 (小有名气)

[求助] 求高版本的1stopt代跑一趟。谢谢 已有1人参与

我的反应是在固定床反应器中，每一个条件下只知道停留时间和反应进出口的浓度，求九个参数。总共做了12组实验点进行拟合。
Title "Kinetics";
//Parameters k1,k2,k3,k4,k5,k6,k7,k8,k9,;
//Variable t,x1,x2,x3,x4,x5,x6;
//ODEFunction
           r1 = k7*k1*k2*x1*x2 /(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;
           r2 = k8*k1*k4*x1*x4/((1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;
           r3 = k9*k1**k5x1*k5/((1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;
           dx1/dt = -2*r1-2*r2-r3;
           dx2/dt = -r1;
           dx3/dt = r1+r2;
           dx4/dt = r1-r2;
           dx5/dt = r2-r3;
           dx6/dt = r3;
Data
  t    x1in      x2in   x3in    x4in     x5in   x6in     x1out     x2out         x3out      x4out        x5out         x6out
2.05         0.92932        0.00258        0.06810        0        0        0        0.9282         0.0000088         0.0695         0.000031         0.001164         0.001137
1.37         0.93029        0.00155        0.06817        0        0        0        0.9283         0.0000142         0.0697         0.000078         0.001348         0.000573
1.62         0.91758        0.00173        0.08068        0        0        0        0.9150         0.0000168         0.0827         0.000104         0.001727         0.000475
1.98         0.89928        0.00187        0.09884        0        0        0        0.8952         0.0000488         0.1017         0.000289         0.002363         0.000444
1.21         0.91772        0.00158        0.08070        0        0        0        0.9142         0.0000739         0.0831         0.000414         0.001897         0.000280
1.49         0.89895        0.00224        0.09881        0        0        0        0.8952         0.0001940         0.1017         0.000641         0.002019         0.000292
2.70         0.81544        0.00531        0.17925        0        0        0        0.8116         0.0007259         0.1837         0.001304         0.002320         0.000368
1.02         0.92936        0.00254        0.06810        0        0        0        0.9274         0.0000399         0.0703         0.000206         0.001776         0.000361
0.82         0.92918        0.00273        0.06809        0        0        0        0.9272         0.0001846         0.0703         0.000465         0.001608         0.000189
0.97         0.91722        0.00213        0.08065        0        0        0        0.9144         0.0003093         0.0830         0.000790         0.001376         0.000149
1.19         0.89756        0.00379        0.09865        0        0        0        0.8947         0.0007324         0.1017         0.000972         0.001730         0.000210
2.16         0.81436        0.00663        0.17902        0        0        0        0.8103         0.0017511         0.1838         0.001744         0.002119         0.000349
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-01-03 20:31:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
引用回帖:
7楼: Originally posted by zhaoshazhu at 2015-01-05 16:03:57
我第一次是写的r3 = k9*k1*k5*x1*k5/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;，您的代码里也是这样的，分子有两个k5，需要把其中的一个改为x5....

均方差(RMSE):0.00106616879873234
残差平方和(SSE):8.18435453321056E-5
相关系数(R): 0
相关系数之平方(R^2): 0
确定系数(DC): 0
F统计(F-Statistic): 0

参数                  最佳估算
--------------------        -------------
k1        2.65934369026499
k2        58012.3883745
k3        12.3953925501171
k4        910.506489281551
k5        1139098.80593599
k6        310.722381455093
k7        0.630067032131075
k8        77.4677237219543
k9        3.42700408019413E-16
赞一下 回复此楼
8楼2015-01-05 16:52:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhaoshazhu: 金币+10, ★★★很有帮助 2015-01-04 16:59:48
不知道理解的对否：
CODE:
Title "Kinetics";
Parameters k1,k2,k3,k4,k5,k6,k7,k8,k9;
ParameterDomain = [0,];
Variable t,x1,x2,x3,x4,x5,x6;
ConstStr r1 = k7*k1*k2*x1*x2/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2,
         r2 = k8*k1*k4*x1*x4/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2,
         r3 = k9*k1*k5*x1*k5/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;
ODEFunction
           x1' = -2*r1-2*r2-r3;
           x2' = -r1;
           x3' = r1+r2;
           x4' = r1-r2;
           x5' = r2-r3;
           x6' = r3;
Data;
//t        x1in        x2in        x3in        x4in        x5in        x6in
0        0.92932        0.00258        0.06810        0        0        0
2.05        0.9282        0.0000088        0.0695        0.000031        0.001164        0.001137
Data;
0        0.93029        0.00155        0.06817        0        0        0
1.37        0.9283        0.0000142        0.0697        0.000078        0.001348        0.000573
Data;
0        0.91758        0.00173        0.08068        0        0        0
1.62        0.9150        0.0000168        0.0827        0.000104        0.001727        0.000475
Data;
0        0.89928        0.00187        0.09884        0        0        0
1.98        0.8952        0.0000488        0.1017        0.000289        0.002363        0.000444
Data;
0        0.91772        0.00158        0.08070        0        0        0
1.21        0.9142        0.0000739        0.0831        0.000414        0.001897        0.000280
Data;
0        0.89895        0.00224        0.09881        0        0        0
1.49        0.8952        0.0001940        0.1017        0.000641        0.002019        0.000292
Data;
0        0.81544        0.00531        0.17925        0        0        0
2.70        0.8116        0.0007259        0.1837        0.001304        0.002320        0.000368
Data;
0        0.92936        0.00254        0.06810        0        0        0
1.02        0.9274        0.0000399        0.0703        0.000206        0.001776        0.000361
Data;
0        0.92918        0.00273        0.06809        0        0        0
0.82        0.9272        0.0001846        0.0703        0.000465        0.001608        0.000189
Data;
0        0.91722        0.00213        0.08065        0        0        0
0.97        0.9144        0.0003093        0.0830        0.000790        0.001376        0.000149
Data;
0        0.89756        0.00379        0.09865        0        0        0
1.19        0.8947        0.0007324        0.1017        0.000972        0.001730        0.000210
Data;
0        0.81436        0.00663        0.17902        0        0        0
2.16        0.8103        0.0017511        0.1838        0.001744        0.002119        0.000349

参数都要求大于0吧？感觉模型可简化下。

均方差(RMSE):0.00107423565746257
残差平方和(SSE):8.3086721839011E-5

参数                  最佳估算
--------------------        -------------
k1        1.90714390662338
k2        2448.35355734201
k3        2.70694132686204
k4        603360.71453709
k5        0.0435072616240413
k6        129731241697.305
k7        23.7865519895587
k8        0.143945603891662
k9        0.0565240940660268
赞一下 回复此楼
2楼2015-01-04 14:56:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoshazhu

新虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2015-01-04 14:56:16
不知道理解的对否：

Title "Kinetics";
Parameters k1,k2,k3,k4,k5,k6,k7,k8,k9;
ParameterDomain = ;
Variable t,x1,x2,x3,x4,x5,x6;
ConstStr r1 = k7*k1*k2*x1*x2/(1+k1*x1+k2*x2+k3*x3+k4*x4+ ...

实在不好意思。第三个方程是 r3 = k9*k1*k5*x1*x5/(1+k1*x1+k2*x2+k3*x3+k4*x4+k5*x5+k6*x6)^2;可以再帮我算一次吗？参数全部需要大于零的。
赞一下 回复此楼
3楼2015-01-04 16:44:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoshazhu

新虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2015-01-04 14:56:16
不知道理解的对否：

Title "Kinetics";
Parameters k1,k2,k3,k4,k5,k6,k7,k8,k9;
ParameterDomain = ;
Variable t,x1,x2,x3,x4,x5,x6;
ConstStr r1 = k7*k1*k2*x1*x2/(1+k1*x1+k2*x2+k3*x3+k4*x4+ ...

非感谢您帮我。就是这样理解的。我这几个参数都需要计算出来所以模型应该是不能简化。
赞一下 回复此楼
4楼2015-01-04 16:59:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭