| 查看: 1112 | 回复: 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人回复
-
vasp极化强度
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有247人回复
-
vasp计算的一些问题
已经有0人回复
-
国科温州研究院 项目聘用有意“二战”的优秀本科毕业生
已经有31人回复
-
SCAPS-1D学习交流
已经有2人回复
-
深圳大学跟哈工大联合博士
已经有8人回复
-
Journal of alloys and compounds期刊
已经有1人回复
5楼2020-08-05 14:26:02
独孤神宇
版主 (知名作家)
- 应助: 488 (硕士)
- 贵宾: 0.008
- 金币: 29542.8
- 散金: 802
- 红花: 120
- 沙发: 1
- 帖子: 5339
- 在线: 850.2小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
2楼2020-08-04 21:25:35
hzlhm
至尊木虫 (著名写手)
- 应助: 373 (硕士)
- 金币: 16178.8
- 红花: 52
- 帖子: 2570
- 在线: 572.1小时
- 虫号: 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
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15356.8
- 散金: 101
- 红花: 234
- 帖子: 3409
- 在线: 1222.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
20