24小时热门版块排行榜    

查看: 1591  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

407369800

金虫 (初入文坛)

[求助] 为什么1st opt 每次拟合出的结果都不一样 已有1人参与

小虫正在使用1stopt 拟合一个动力学模型,每次出来的结果都不一样。现在我自己怀疑是不是有过拟合的可能,各位大神给看看呢。感激不尽。
Parameters kx[0,1],kp[0,1],s0,m;
Variables t,y;
Function y=s0-((1/kx+0.1839/kp)*((13.55492*exp(0.26739*t)/(31.14889+0.42925*exp(0.26739*t)))-0.42925))-(((m+0.00252/kp))*(118.09769*ln((31.14889+0.42925*exp(0.26739*t))/31.57814)));
Data;
t,y
0        63.63497
4        59.44704
8        53.82187
12        46.31953
16        34.01546
20        23.00973
24        13.40701
28        8.7705
32      7.1334
36        4.81704
40        4.79449
44        4.45433
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

对于数据和公式模型,自变量前面的待定系数应该是通过拟合是确定的,但是这个待定系数又有几个数进行数学运算,从数学角度讲是有无数解的。
你可以试试把kx或者m取消试试,取消一个后拟合结果唯一。
一团和气
4楼2015-09-15 10:00:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
1/kx+0.1839/kp和m+0.00252/kp可以简化为p1、p2,这样可以得出稳定结果。
迭代数: 16
计算用时(时:分:秒:微秒): 00:00:00:254
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 4.60972692793589
残差平方和(SSE): 254.994988201647
相关系数(R): 0.980557335475958
相关系数之平方(R^2): 0.96149268815571
决定系数(DC): 0.957200935742808
卡方系数(Chi-Square): 15.7773168503866
F统计(F-Statistic): 249.690939747664

参数        最佳估算
----------        -------------
p1                 1
p2                 0.0285957945196896
s0                 55.0667798764375
一团和气
2楼2015-09-14 12:47:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

407369800

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by lipenggg at 2015-09-14 12:47:23
1/kx+0.1839/kp和m+0.00252/kp可以简化为p1、p2,这样可以得出稳定结果。
迭代数: 16
计算用时(时:分:秒:微秒): 00:00:00:254
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收 ...

这样其实就是四个参数变成了三个参数了。
真的是参数设置过多,过拟合的原因吗?
其实我每个参数都有意义的,如果变成三个参数感觉没法用了诶。求指教啊。
3楼2015-09-14 15:04:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

建议你看看s0、kx、kp、m的物理意义,是否可以在其他公式里能确定m或kx或kp。这样公式结果就唯一了。
或者多算几次,将其中的一个认为是合理值确定下来,带入公式再次拟合,这样可以确定其他三个数值的唯一解。
一团和气
5楼2015-09-15 10:11:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研调剂 +4 呼呼?~+123456 2026-03-24 4/200 2026-03-24 23:16 by barlinike
[考研] 281求调剂 +3 Koxui 2026-03-24 4/200 2026-03-24 22:59 by barlinike
[考研] 085601求调剂总分293英一数二 +3 钢铁大炮 2026-03-24 3/150 2026-03-24 22:03 by bingxueer79
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +3 哈哈157349 2026-03-21 3/150 2026-03-24 14:11 by zhyzzh
[考研] 279分求调剂 一志愿211 +18 chaojifeixia 2026-03-19 20/1000 2026-03-24 10:34 by dolphin_ycj
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗? +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +7 vv迷 2026-03-22 7/350 2026-03-23 23:44 by Txy@872106
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +5 niko- 2026-03-22 5/250 2026-03-23 22:01 by fuyu_
[考研] 352求调剂 +3 大米饭! 2026-03-22 3/150 2026-03-22 23:28 by king123!
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 285求调剂 +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by 星空星月
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
信息提示
请填处理意见