| 查看: 1948 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
1stopt使用时拟合计算不稳定 已有3人参与
|
|||
|
我在用1stopt对干燥方程进行拟合时，每次计算之后所得常数都会发生变化，数值不稳定。所用方程为经验方程可能出现过拟合情况。请大佬帮忙看一下，是我软件版本低的情况吗？马上快毕业了，比较急求大佬指点万分感谢。 Variable t,m; Parameters a,k,b,g,c,h; Function m=a*exp(-k*t)+b*exp(-g*t)+c*exp(-h*t); Data; 1 0.761468058 2 0.572266754 3 0.429168622 4 0.306584911 5 0.221415875 6 0.153590152 7 0.104075754 8 0.066698653 9 0.038105337 9.5 0.028177319 第一次计算拟合 "Type your title here" 迭代数: 144 计算用时(时:分:秒:微秒): 00:00:05:264 优化算法: 准牛顿法(BFGS) + 通用全局优化法 计算结束原因: 达到收敛判断标准 均方差(RMSE): 0.00167811805641684 残差平方和(SSE): 2.81608021127223E-5 相关系数(R): 0.999974701263212 相关系数之平方(R^2): 0.999949403166451 决定系数(DC): 0.99994940316645 卡方系数(Chi-Square): 4.48676514663589E-5 F统计(F-Statistic): 158104.6611854 参数 最佳估算 ---------- ------------- a -3.91637356524864 k 0.409125168268213 b 4.93177468581318 g 0.374949664548045 c -0.0267488857371107 h -0.0180087835359877 ====== 结果输出 ===== No 实测值m 计算值m 1 0.761468058 0.7611145 2 0.572266754 0.5741946 3 0.429168622 0.4253856 4 0.306584911 0.3095410 5 0.221415875 0.2208491 6 0.153590152 0.1538050 7 0.104075754 0.1036216 8 0.066698653 0.0663422 9 0.038105337 0.0388043 10 0.028177319 0.0278935 第二次计算拟合 "Type your title here" 迭代数: 196 计算用时(时:分:秒:微秒): 00:00:05:617 优化算法: 准牛顿法(BFGS) + 通用全局优化法 计算结束原因: 达到收敛判断标准 均方差(RMSE): 0.00167811998989461 残差平方和(SSE): 2.81608670048389E-5 相关系数(R): 0.999974701204915 相关系数之平方(R^2): 0.999949403049858 决定系数(DC): 0.999949403049858 卡方系数(Chi-Square): 4.48675487738193E-5 F统计(F-Statistic): 158104.296839769 参数 最佳估算 ---------- ------------- a -0.0267533848791293 k -0.0179991625982107 b -3.88507793264088 g 0.409264735572213 c 4.90048287680716 h 0.374841529171177 ====== 结果输出 ===== No 实测值m 计算值m 1 0.761468058 0.7611144 2 0.572266754 0.5741946 3 0.429168622 0.4253856 4 0.306584911 0.3095410 5 0.221415875 0.2208491 6 0.153590152 0.1538051 7 0.104075754 0.1036217 8 0.066698653 0.0663422 9 0.038105337 0.0388043 10 0.028177319 0.0278935 第三次计算拟合 "Type your title here" 迭代数: 201 计算用时(时:分:秒:微秒): 00:00:05:610 优化算法: 准牛顿法(BFGS) + 通用全局优化法 计算结束原因: 达到收敛判断标准 均方差(RMSE): 0.00167810849402698 残差平方和(SSE): 2.81604811772551E-5 相关系数(R): 0.999974701551534 相关系数之平方(R^2): 0.99994940374308 决定系数(DC): 0.999949403743079 卡方系数(Chi-Square): 4.48681005484139E-5 F统计(F-Statistic): 158106.463142428 参数 最佳估算 ---------- ------------- a -0.0267168568795642 k -0.0180845713175149 b -4.08185733520843 g 0.408430070251824 c 5.09722825774971 h 0.37550527538769 ====== 结果输出 ===== No 实测值m 计算值m 1 0.761468058 0.7611144 2 0.572266754 0.5741945 3 0.429168622 0.4253857 4 0.306584911 0.3095411 5 0.221415875 0.2208491 6 0.153590152 0.1538050 7 0.104075754 0.1036216 8 0.066698653 0.0663422 9 0.038105337 0.0388044 10 0.028177319 0.0278935 |
回复此楼» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有204人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有3人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17808.5
- 红花: 53
- 帖子: 2879
- 在线: 606.3小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
独孤神宇: 金币+5, 鼓励交流 2020-08-05 08:29:42
小木虫: 金币+0.5, 给个红包，谢谢回帖
独孤神宇: 金币+5, 鼓励交流 2020-08-05 08:29:42
|
如用matlab做的话，结果变化不大，这个跟初值有关系。 第1次计算拟合 a b c g h k -5.48855631156569 4.4979065232128 1.92482409022818 0.469010073268637 0.44404365198402 0.511719978797606 第2次计算拟合 -5.48855631156569 4.4979065232128 1.92482409022818 0.469010073268637 0.44404365198402 0.511719978797606 第3次计算拟合 -5.48855631156569 4.4979065232128 1.92482409022818 0.469010073268637 0.44404365198402 0.511719978797606  123.jpg |
3楼2020-08-05 00:00:17
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
2楼2020-08-04 21:25:35
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
独孤神宇: 金币+2, 鼓励交流 2020-08-05 09:31:02
小木虫: 金币+0.5, 给个红包，谢谢回帖
独孤神宇: 金币+2, 鼓励交流 2020-08-05 09:31:02
|
应该一致的，只是顺序不同而已（公式中三个部分形式是一样的），另外指数数值计算，有可能出现计算精度差异。 Root of Mean Square Error (RMSE): 0.00167798458106697 Sum of Squared Residual: 2.81563225429849E-5 Correlation Coef. (R): 0.999974705287561 R-Square: 0.999949411214945 Adjusted R-Square: 0.999934957276358 Determination Coef. (DC): 0.999949411214944 Parameter Best Estimate ---------- ------------- a 84.9030954163188 k 0.390313878805315 b -83.8883211096539 g 0.392092581899161 c -0.0261161738829758 h -0.0196386740358328 |
4楼2020-08-05 09:27:01
5楼2020-08-05 14:26:02
