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

登录注册

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

24小时热门版块排行榜

返回列表

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

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

[求助] 求高版本的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

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有292人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1stopt 语句修改已经有23人回复
求助 1stopt低版本方程的拟合（含若干未知参数）已经有10人回复
关于用1stpro拟合数值的问题。已经有4人回复
1stopt自变量超多时是否有简单语法？AND 自变量多数据量大的话1stopt都会算的特慢么已经有7人回复
1stopt数据拟合求助已经有8人回复
求大神跑下1stOpt，一个sharemodel，谢谢！已经有7人回复
跪求高手帮忙看看这组数据如何进行拟合已经有7人回复
求高版本1stopt跑一个程序已经有7人回复
1stopt解方程组已经有12人回复
用1stOpt进行多模型多变量非线性拟合得不到结果已经有4人回复
求用高版本的1stOpt拟合方程已经有8人回复
1stopt 拟合计算已经有6人回复
求高版本1stopt拟合，已经有11人回复
1stOpt程序的一个小问题已经有5人回复
跪求1stOpt2.5以上版本，或者请高人帮忙拟合包含复数的函数的参数已经有8人回复

1楼 2015-01-03 20:31:42

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

zhaoshazhu

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

引用回帖:

8楼: Originally posted by dingd at 2015-01-05 16:52:09
均方差(RMSE):0.00106616879873234
残差平方和(SSE):8.18435453321056E-5
相关系数(R): 0
相关系数之平方(R^2): 0
确定系数(DC): 0
F统计(F-Statistic): 0

参数最佳估算
-------------------- ...

非常感谢，可以帮我算另外一个模型的吗？数据一样的，只是没有x6组分，模型不一样
Title "Kinetics";
Parameters k1,k2,k3,k4,k5,k6,k7,k8,k9;
ParameterDomain = [0,];
Variable t,x1,x2,x3,x4,x5;
ConstStr
r1 = k(6).*((k(1).*x(1).*k(2).*x(2)).^0.5)./...
(1+2.*(k(2).*x(2)).^0.5+2.*(k(4).*x(4)).^0.5+(k(1).*x(1)).^0.5+k(3).*x(3)+k(5).*x(5)).^2;
r2 = k(7).*((k(1).*x(1).*k(4).*x(4)).^0.5)./...
(1+2.*(k(2).*x(2)).^0.5+2.*(k(4).*x(4)).^0.5+(k(1).*x(1)).^0.5+k(3).*x(3)+k(5).*x(5)).^2;

r3 = k(8).*k(5).*k(1).*x(5).*x(1)./...
(1+2.*(k(1).*x(1)).^0.5+2.*(k(4).*x(4)).^0.5+(k(2).*x(2)).^0.5+k(3).*x(3)+k(5).*x(5)).^2;
ODEFunction
         x1' = -2*r1-2*r2-r3;
         x2' = -r1;
         x3' = r1+r2;
         x4' = r1-r2;
         x5' = r2-r3;
Data;
//t       x1in       x2in       x3in       x4in       x5in
0       0.92932       0.00258       0.06810       0       0
2.05       0.9282       0.0000088       0.0695       0.000031       0.001164
Data;
0       0.93029       0.00155       0.06817       0       0
1.37       0.9283       0.0000142       0.0697       0.000078       0.001348
Data;
0       0.91758       0.00173       0.08068       0       0
1.62       0.9150       0.0000168       0.0827       0.000104       0.001727
Data;
0       0.89928       0.00187       0.09884       0       0
1.98       0.8952       0.0000488       0.1017       0.000289       0.002363
Data;
0       0.91772       0.00158       0.08070       0       0
1.21       0.9142       0.0000739       0.0831       0.000414       0.001897
Data;
0       0.89895       0.00224       0.09881       0       0
1.49       0.8952       0.0001940       0.1017       0.000641       0.002019
Data;
0       0.81544       0.00531       0.17925       0       0
2.70       0.8116       0.0007259       0.1837       0.001304       0.002320
Data;
0       0.92936       0.00254       0.06810       0       0
1.02       0.9274       0.0000399       0.0703       0.000206       0.001776
Data;
0       0.92918       0.00273       0.06809       0       0
0.82       0.9272       0.0001846       0.0703       0.000465       0.001608
Data;
0       0.91722       0.00213       0.08065       0       0
0.97       0.9144       0.0003093       0.0830       0.000790       0.001376
Data;
0       0.89756       0.00379       0.09865       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
2.16       0.8103       0.0017511       0.1838       0.001744       0.002119

赞一下

回复此楼

9楼2015-01-05 19:16:57

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

查看全部 10 个回答

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 125.7
散金: 10
帖子: 51
在线: 38小时
虫号: 1812564
注册: 2012-05-12
专业: 环境化工

引用回帖:

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的回帖