24小时热门版块排行榜     意得辑论文润色满3千送无线充电宝

查看: 350  |  回复: 7

976246540

新虫 (初入文坛)

[交流] 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
数据较少,可以尝试设置参数的范围进行拟合。

版本较低也可能是一个原因

发自小木虫Android客户端
提供MATLAB、1stOpt数值计算有偿服务,请私信留言。
2楼2020-08-04 21:25:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

铁杆木虫 (著名写手)

★ ★ ★ ★ ★ ★
小木虫: 金币+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
1stopt使用时拟合计算不稳定
123.jpg

QQ:2120156492
3楼2020-08-05 00:00:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

★ ★ ★
小木虫: 金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

976246540

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by dingd at 2020-08-05 09:27:01
应该一致的,只是顺序不同而已(公式中三个部分形式是一样的),另外指数数值计算,有可能出现计算精度差异。

Root of Mean Square Error (RMSE): 0.00167798458106697
Sum of Squared Residual: 2.81563225429 ...

嗯,好像只是顺序变了。但是你的a值和b值比我拟合的要大

发自小木虫Android客户端
5楼2020-08-05 14:26:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

976246540

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by hzlhm at 2020-08-05 00:00:17
如用matlab做的话,结果变化不大,这个跟初值有关系。
第1次计算拟合
a   b   c   g   h   k
-5.48855631156569           4.4979065232128          1.92482409022818         0.469010073268637          0.4 ...

初值是什么不太懂。求赐教,谢谢

发自小木虫Android客户端
6楼2020-08-05 14:33:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
看目标函数残差平方和或均方差,越小越好。

发自手机版WAP
7楼2020-08-05 15:11:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

铁杆木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by 976246540 at 2020-08-05 14:33:48
初值是什么不太懂。求赐教,谢谢
...

初值,就是拟合系数计算最初的值(开始的第一个值)。

发自小木虫Android客户端
QQ:2120156492
8楼2020-08-05 15:14:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 976246540 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 男博士超过35岁毕业参加工作,无缘青基,还可以申请哪些基金? +9 tmgczb 2020-09-30 19/950 2020-09-30 18:12 by spj860iuy
[基金申请] 求各位大神看看上会了吗 +9 fxysnd 2020-09-30 11/550 2020-09-30 17:39 by jyfzhen
[硕博家园] 读博好压抑,想不出好的idea,严重怀疑自己,也没有退路 +5 想要一只喵 2020-09-29 5/250 2020-09-30 17:34 by huanlele
[硕博家园] 深圳还是成都 +10 cxg020928 2020-09-23 10/500 2020-09-30 16:33 by joyair1991
[硕博家园] 高校的教师和科研岗选择 +5 lcn1013 2020-09-29 7/350 2020-09-30 16:24 by 产品可靠
[硕博家园] 国内科研变天了 +8 龚家湾 2020-09-30 9/450 2020-09-30 15:44 by lcn1013
[论文投稿] 关于会议投稿作者漏写 20+3 星缘阿巴查 2020-09-29 7/350 2020-09-30 14:52 by lyfwish
[基金申请] 项目任务书填报常见问题解答 +3 天天向上3035 2020-09-24 7/350 2020-09-30 10:44 by 小明大哥大
[论文投稿] 请帮我看看这个审稿意见是比较正面的还是负面的 +3 yanximing 2020-09-29 6/300 2020-09-30 07:43 by 国际科学编辑
[材料综合] 原位透射的加热速率最快能多少 30+4 子宁访天 2020-09-28 7/350 2020-09-30 06:54 by 迎海蓉城
[有机交流] 求合成方法,Br变OMe 100+3 啄木鸟先森 2020-09-28 5/250 2020-09-29 21:02 by justforjoy
[公派出国] 冒险冲国外高校还是转国内博士? 30+5 ZéphyrTin 2020-09-27 7/350 2020-09-29 20:51 by kwj828
[教师之家] 德国布伦瑞克工业大学, 氢氧PEM电池仿真方向,访问学者机会,德国这边给薪水 +3 gaoxincn 2020-09-24 6/300 2020-09-29 12:57 by gaoxincn
[有机交流] 如何防护氟化氢吡啶? 500+3 lishaomen 2020-09-26 8/400 2020-09-29 09:22 by xiaobenpu
[论文投稿] Reject & Resubmit是否相当于大修 150+3 风之殇k 2020-09-26 4/200 2020-09-28 07:58 by bobvan
[教师之家] 散币销号 +6 lipcaty 2020-09-24 6/300 2020-09-28 07:47 by 昔年。。
[基金申请] 医学部可查了 +3 zhujx81 2020-09-27 6/300 2020-09-27 11:41 by tulv022
[基金申请] 青年基金评审意见分析求助 +11 玄影 2020-09-25 11/550 2020-09-27 10:29 by Doctor_Lei
[基金申请] 青基反馈意见,大家给点意见 +7 wangjunqi 2020-09-24 8/400 2020-09-25 08:14 by 王小艾~
[基金申请] 心电图近半小时飞速增长 +6 DRJM 2020-09-23 13/650 2020-09-24 06:35 by DRJM
信息提示
请填处理意见