24小时热门版块排行榜    

查看: 1806  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 06:55 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:46 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 06:22 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 02:40 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:25 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:52 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见