| 查看: 1146 | 回复: 3 | |||
l1003785517新虫 (小有名气)
|
[求助]
用1stopt计算非线性方程组的未知参数 已有1人参与
|
|
我用1stopt（1.5版本） 计算动力学方程的各参数值，不同温度下对应不同的反应速率，但是我算出来的有的参数为负数，相关系数很低，感觉很不可靠，而且为什么温度也出现在计算的参数里面？第一次用1stopt软件，不是很懂，希望大家给我提供修改意见，感激不尽！附上代码。 Constant R=8.3145; VarConstant T=[563.15,583.15,598.15,608.15]; ConstStr k1=k0*exp(-Ea1/R/T),k2=k*exp(-Ea2/R/T); ShareModel; Variables x(1:3),r1,r2; Parameters a(1:8); Function r1=k1*(x1^a1)/(1+a2*(x2^a3)*(x3^a4)); r2=r1-k2*(x1^a5)/(1+a6*(x2^a7)*(x3^a8)); Data; // x1,x2,x3,r1,r2 6.34514148 0.01212732 0.0255312 0.2864 0.0245 6.10642476 0.02744604 0.2489292 0.2475 0.0214 5.97876876 0.03382884 0.3702024 0.2086 0.0184 5.89260096 0.04340304 0.446796 0.1698 0.0153 5.7796254 0.0542538 0.5489208 0.1077 0.0092 5.67941544 0.06510456 0.63828 0.0926 0.0073 5.59899216 0.06893424 0.7148736 0.0776 0.0054 5.51601576 0.07531704 0.7914672 0.0735 0.0034 5.3711262 0.079785 0.9318888 0.0654 0.0025 5.25751236 0.08489124 1.0403964 0.0574 0.0016 5.14708992 0.08680608 1.148904 0.0493 0.0007 Data; 6.20791128 0.02170152 0.1531872 0.4795 0.0486 6.08536152 0.05489208 0.2425464 0.4473 0.0429 5.90281344 0.06510456 0.414882 0.3829 0.0315 5.71835052 0.07084908 0.5936004 0.3185 0.0201 5.55175944 0.07787016 0.7531704 0.2541 0.0087 5.34495672 0.08680608 0.9510372 0.1962 0.003 5.170068 0.0893592 1.1233728 0.162 0.0029 4.93581924 0.09382716 1.3531536 0.1063 0.0028 4.6945494 0.0989334 1.5893172 0.089 0.0025 4.56753168 0.10467792 1.7105904 0.0717 0.0023 4.4296632 0.1085076 1.8446292 0.0544 0.0021 Data; 6.04706472 0.06127488 0.2744604 0.6009 0.0288 5.76558324 0.08106156 0.5361552 0.5236 0.0242 5.48346348 0.08872092 0.8106156 0.4464 0.0196 5.3264466 0.0925506 0.9638028 0.3691 0.015 5.17517424 0.10340136 1.1042244 0.2919 0.0104 4.8732678 0.111699 1.3978332 0.2265 0.0069 4.70795328 0.11744352 1.5574032 0.1908 0.0053 4.37604768 0.12382632 1.882926 0.1324 0.0028 4.12137396 0.12957084 2.1318552 0.1128 0.0025 3.94265556 0.13595364 2.3041908 0.0932 0.0023 3.75500124 0.13850676 2.489292 0.0735 0.002 Data; 5.64750144 0.08425296 0.6510456 0.5046 0.048 5.3647434 0.1053162 0.9127404 0.4439 0.0389 5.217939 0.1244646 1.0403964 0.3832 0.0298 4.97347776 0.13914504 1.2701772 0.3226 0.0206 4.92624504 0.14169816 1.3148568 0.2619 0.0115 4.59433944 0.14808096 1.6403796 0.2104 0.0063 4.27839084 0.16403796 1.9403712 0.1542 0.0047 4.02818508 0.17169732 2.1829176 0.1227 0.0034 3.8009574 0.175527 2.4063156 0.1162 0.0026 3.4850088 0.1787184 2.7190728 0.1096 0.0017 3.35862936 0.18382464 2.840346 0.103 0.0008 |
回复此楼» 收录本帖的淘帖专辑推荐
 常微分方程 |
» 猜你喜欢
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有673人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
物理学I论文润色/翻译怎么收费?
已经有272人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
l1003785517: 金币+50, ★★★很有帮助 2017-09-26 13:34:13
感谢参与，应助指数 +1
l1003785517: 金币+50, ★★★很有帮助 2017-09-26 13:34:13
|
k0，Ea1，k，Ea2也都是参数？那就有问题了，会出现过拟合现象导致参数组值不唯一，下面参考下： 均方差(RMSE): 0.0229100483751911 残差平方和(SSE): 0.0461885878567165 相关系数(R): 0.972252575246286 相关系数之平方(R^2): 0.945275070073034 修正R平方(Adj. R^2): 0.891494162615898 确定系数(DC): 0.826086223956403 F统计(F-Statistic): -2.82644079611069 参数 最佳估算 -------------------- ------------- k0 3447893.11032413 ea1 91178.73721166 k 4978920.37512857 ea2 91224.8275413718 a1 1.70639628682468 a2 1.15955068548675 a3 7.13855269375096E-21 a4 1.22141382707052 a5 1.45736312486719 a6 1.32314168468273 a7 0.0562010088861979 a8 1.24071362634774 |
2楼2017-09-25 09:33:04
l1003785517
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 11.4
- 散金: 30
- 帖子: 59
- 在线: 22.1小时
- 虫号: 4326948
- 注册: 2016-01-01
- 专业: 能源化工
3楼2017-09-25 21:17:40
l1003785517
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 11.4
- 散金: 30
- 帖子: 59
- 在线: 22.1小时
- 虫号: 4326948
- 注册: 2016-01-01
- 专业: 能源化工
4楼2017-09-26 13:58:55
