| 查看: 1351 | 回复: 9 | ||
[求助]
求高版本的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 |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有255人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有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人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhaoshazhu: 金币+10, ★★★很有帮助 2015-01-04 16:59:48
感谢参与，应助指数 +1
zhaoshazhu: 金币+10, ★★★很有帮助 2015-01-04 16:59:48
|
不知道理解的对否： 参数都要求大于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
3楼2015-01-04 16:44:04
4楼2015-01-04 16:59:24
5楼2015-01-05 10:57:38
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
6楼2015-01-05 14:37:12
7楼2015-01-05 16:03:57
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
|
均方差(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
|
非常感谢，可以帮我算另外一个模型的吗？数据一样的，只是没有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
10楼2015-01-07 21:26:53
