| 查看: 960 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 常微分方程 |
» 猜你喜欢
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有179人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有23人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
l1003785517: 金币+50, ★★★很有帮助 2017-09-26 13:34:13
感谢参与，应助指数 +1
l1003785517: 金币+50, ★★★很有帮助 2017-09-26 13:34:13
l1003785517
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 11.4
- 散金: 30
- 帖子: 59
- 在线: 22.1小时
- 虫号: 4326948
- 注册: 2016-01-01
- 专业: 能源化工
l1003785517
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 11.4
- 散金: 30
- 帖子: 59
- 在线: 22.1小时
- 虫号: 4326948
- 注册: 2016-01-01
- 专业: 能源化工