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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » 用1stopt计算非线性方程组的未知参数
41/1返回列表
查看: 1001  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 收录本帖的淘帖专辑推荐

常微分方程

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2017-09-24 20:48:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

l1003785517

新虫 (小有名气)
恩，是的，都是未知参数。能把源代码给我看看吗？我的1.5版本运算出来参数都有负数，不知道是不是我的代码写错了。

发自小木虫IOS客户端
赞一下 回复此楼
3楼2017-09-25 21:17:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

l1003785517

新虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2017-09-25 09:33:04
k0，Ea1，k，Ea2也都是参数？那就有问题了，会出现过拟合现象导致参数组值不唯一，下面参考下：
均方差(RMSE): 0.0229100483751911
残差平方和(SSE): 0.0461885878567165
相关系数(R): 0.972252575246286
相关系 ...

恩，是的，都是未知参数。能把源代码给我看看吗？我的1.5版本运算出来参数都有负数，不知道是不是我的代码写错了。
赞一下 回复此楼
4楼2017-09-26 13:58:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 l1003785517 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +7 yangfz 2026-03-17 7/350 2026-03-17 23:57 by 星空星月
[考研] 能源材料化学课题组招收硕士研究生8-10名 +3 脱颖而出 2026-03-16 6/300 2026-03-17 21:19 by z1z2z3879
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 085600材料与化工求调剂 +5 绪幸与子 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[考研] 求调剂，总分315，考的生物医药，一志愿湖南师范大学。调剂到任何专业都可以 +4 小丁想进步 2026-03-11 5/250 2026-03-17 16:05 by 外星文明
[考研] 工科材料085601 279求调剂 +3 困于星晨 2026-03-17 3/150 2026-03-17 14:08 by ms629
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] （081700）化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 283求调剂，材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
信息提示
关闭
请填处理意见
关闭