| 查看: 995 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 常微分方程 |
» 猜你喜欢
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有268人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有0人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
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
10