24小时热门版块排行榜    

查看: 1599  |  回复: 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的回帖

407369800

金虫 (初入文坛)

引用回帖:
5楼: Originally posted by lipenggg at 2015-09-15 10:11:23
建议你看看s0、kx、kp、m的物理意义,是否可以在其他公式里能确定m或kx或kp。这样公式结果就唯一了。
或者多算几次,将其中的一个认为是合理值确定下来,带入公式再次拟合,这样可以确定其他三个数值的唯一解。

我现在的办法是把后面的m+0.00252/kp简化成p2,然后从前面得到kx和kp的值,最后,再通过p2和kp算出m。
7楼2015-09-15 16:23:23
已阅   回复此楼   关注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

专家顾问 (正式写手)

对于数据和公式模型,自变量前面的待定系数应该是通过拟合是确定的,但是这个待定系数又有几个数进行数学运算,从数学角度讲是有无数解的。
你可以试试把kx或者m取消试试,取消一个后拟合结果唯一。
一团和气
4楼2015-09-15 10:00:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 312求调剂 +5 上岸吧ZJY 2026-03-22 7/350 2026-03-25 22:20 by 544594351
[考研] 材料求调剂 +4 .m.. 2026-03-25 4/200 2026-03-25 21:30 by peike
[考研] 求b区院校调剂 +4 周56 2026-03-24 5/250 2026-03-25 17:12 by yishunmin
[考研] 296求调剂 +4 汪!?! 2026-03-25 7/350 2026-03-25 16:41 by 汪!?!
[考研] 0854电子信息求调剂 324 +4 Promise-jyl 2026-03-23 4/200 2026-03-25 11:36 by Sugarlight
[考研] 085600材料与化工调剂 +9 A-哆啦Z梦 2026-03-23 15/750 2026-03-25 11:18 by Ainin_
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +6 @taotao 2026-03-19 6/300 2026-03-25 08:37 by 木托莫露露
[考研] 291求调剂 +3 HanBeiNingZC 2026-03-24 3/150 2026-03-24 16:34 by barlinike
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
信息提示
请填处理意见