24小时热门版块排行榜    

查看: 1807  |  回复: 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

金虫 (初入文坛)

引用回帖:
4楼: Originally posted by lipenggg at 2015-09-15 10:00:49
对于数据和公式模型,自变量前面的待定系数应该是通过拟合是确定的,但是这个待定系数又有几个数进行数学运算,从数学角度讲是有无数解的。
你可以试试把kx或者m取消试试,取消一个后拟合结果唯一。

确实把参数减少一个结果就唯一了,应该就是所谓的过拟合了
6楼2015-09-15 15:14:26
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 07:05 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 06:11 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:25 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 面上意见出来了 +8 黄鸟于飞Chao 2026-08-29 16/800 2026-08-29 20:45 by 黄鸟于飞Chao
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
信息提示
请填处理意见