24小时热门版块排行榜     石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 1113  |  回复: 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的回帖

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的回帖
相关版块跳转 我要订阅楼主 976246540 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 大学直属学院卸任的副院长退休后还享受副处级的养老待遇吗? +4 苏东坡二世 2024-04-27 4/200 2024-04-27 10:12 by ZHONGWU_U
[论文投稿] 求论文投稿期刊推荐 15+3 mihudawang 2024-04-26 4/200 2024-04-27 09:56 by bobvan
[论文投稿] AC审稿 +3 马winner 2024-04-26 6/300 2024-04-27 09:08 by 马winner
[硕博家园] 聊天 +12 暮色恋伊人 2024-04-22 13/650 2024-04-27 08:40 by WASM
[有机交流] 如何分离原料与产品 5+4 qwerasdf587 2024-04-24 13/650 2024-04-26 23:18 by 119966984
[基金申请] 基金开始函评了吗? +16 wych1103 2024-04-25 16/800 2024-04-26 21:32 by 淀粉搬运工
[考研] 没学上 +6 季向阳 2024-04-26 12/600 2024-04-26 21:06 by 季向阳
[论文投稿] 研二光催化6月底四篇二区什么水平 5+5 wjtab 2024-04-22 15/750 2024-04-26 19:25 by wjtab
[考研] 0854-0855调剂 +8 shangannum1 2024-04-21 12/600 2024-04-26 16:42 by yz仔
[考博] 申博求助 +4 dskabdh 2024-04-24 11/550 2024-04-26 15:54 by dskabdh
[基金申请] 两类问题算是白选了~ +7 jurkat.1640 2024-04-23 12/600 2024-04-26 14:39 by lucky_my2010
[考博] 取博导收留 5+4 zzb777888 2024-04-20 10/500 2024-04-26 08:52 by polymerfriend
[基金申请] 国社科项目,你们学校都限额申报吗? +7 屡战屡败 2024-04-21 10/500 2024-04-25 12:10 by 屡战屡败
[电化学] 耗材发问 +4 Happy C 2024-04-22 4/200 2024-04-25 11:03 by 普通小虫
[基金申请] 化学结清有情发出来了 +3 starboy7286 2024-04-20 3/150 2024-04-24 15:52 by 嘿,黑贝
[考博] 博士招生 +4 zx179 2024-04-24 7/350 2024-04-24 15:01 by H考研成功
[教师之家] 大家访学都是怎么找的啊? +3 luokereng 2024-04-22 3/150 2024-04-24 11:40 by xuechenli
[论文投稿] 期刊推荐 20+4 木颜尘ip 2024-04-22 7/350 2024-04-24 10:06 by bobvan
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +4 我属驴核动力驴 2024-04-22 5/250 2024-04-24 08:47 by 晓目崇
[论文投稿] 研究光催化的,好中的三四区 20+3 sl.0117 2024-04-20 3/150 2024-04-22 09:53 by bobvan
信息提示
请填处理意见