24小时热门版块排行榜    

查看: 1124  |  回复: 4

fairless

新虫 (初入文坛)

[求助] 1stopt模拟arrhenius equation偏差较大,如何修正比较准确?

昨天第一次采用1stopt 1.5版拟合变温PL谱强度的激活能拟合,
拟合曲线arrhenius equation如下:

采用1stopt拟合后数据如下:
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 1.86186534632594
残差平方和(SSE): 38.1319682463437
相关系数(R): 0.990530299406935
相关系数之平方(R^2): 0.981150274043193
决定系数(DC): 0.980243427395926
卡方系数(Chi-Square): 2.59569288535077
F统计(F-Statistic): 468.460522270874

参数        最佳估算
----------        -------------
a                  36.4492831714766
b                  13.6635397639502
c                  13.6228437473683

====== 结果输出 =====

No        实测值y        计算值y
1        38.88        36.4492832
2        36.68        36.4361053
3        30.73        34.0552933
4        22.65        23.0894411
5        16.85        15.0121520
6        12.03        10.1607705
7        8.77        7.8227848
8        5.91        5.9882358
9        3.55        5.0636281
10        2.31        4.5160420
11        1.64        4.0855435

接近于0K与300K的数据偏差太大。拟合图片如图:

对于曲线拟合不是特别熟悉,所以想求教下列问题:
1.1stopt本身是否可以拟合更准确一些?
2.是否数据本身的问题,不可以实现整个温度段拟合,分开成两个温度范围的拟合?
十分感谢
1stopt模拟arrhenius equation偏差较大,如何修正比较准确?
捕获.PNG
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
代码及数据一起放上来啊,否则很难判断那出错。
2楼2013-12-08 13:23:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fairless

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2013-12-08 13:23:08
代码及数据一起放上来啊,否则很难判断那出错。

不好意思,现在补上,谢谢您
Parameters a,b,c;
Variable x,y;
Function y=a/(1+b*exp(-c/(0.08617*x)));
Data;
5        38.88
15        36.68
30        30.73
50        22.65
70        16.85
95        12.03
120        8.77
160        5.91
200        3.55
240        2.31
290        1.64
3楼2013-12-08 13:48:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
fairless: 金币+5, ★★★★★最佳答案 2013-12-08 19:27:46
就你的数据和公式,1stOpt得出的已经是最优结果了。要想结果更好,只能从公式和数据方面入手。
比如,如果公式加一个参数变为“y=a/(1+b*exp(-c/(0.08617*x^d)))”,效果就要好得多:

均方差(RMSE):0.216105383119834
残差平方和(SSE):0.513716902747072
相关系数(R): 0.999866962140188
相关系数之平方(R^2): 0.999733941979448
确定系数(DC): 0.9997338378858
卡方系数(Chi-Square): 0.0365348597765041
F统计(F-Statistic): 8766.58701162115

参数        最佳估算
----------        -------------
a        -1.87059593829512
b        -1.04724042759144
c        -2.96611774185407E-6
d        -1.74956375131302

====== 结果输出 ======

No        实测值y        计算值y
1        38.88        39.0987620367877
2        36.68        36.4177135677494
3        30.73        30.5781663655639
4        22.65        22.9168257377651
5        16.85        17.0050738036823
6        12.03        11.9466243952217
7        8.77        8.65818947536124
8        5.91        5.48645366678335
9        3.55        3.67508879898935
10        2.31        2.5590085441019
11        1.64        1.68766073913156
4楼2013-12-08 19:23:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fairless

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by dingd at 2013-12-08 19:23:43
就你的数据和公式,1stOpt得出的已经是最优结果了。要想结果更好,只能从公式和数据方面入手。
比如,如果公式加一个参数变为“y=a/(1+b*exp(-c/(0.08617*x^d)))”,效果就要好得多:

均方差(RMSE):0.216105383 ...

好的,我再试试看,谢谢您的回复
5楼2013-12-08 19:27:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fairless 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +7 Lanmanbaby 2026-08-09 14/700 2026-08-09 23:33 by Tide man
[基金申请] 2026国自然放榜时间 +4 布布和一二 2026-08-08 4/200 2026-08-09 22:24 by otj2008
[基金申请] 关于代码变化问题,想知道的进来 +15 且听虎啸 2026-08-07 22/1100 2026-08-09 22:13 by 布布和一二
[基金申请] 这样的filecode谁见过 +10 布布和一二 2026-08-08 21/1050 2026-08-09 17:21 by 天神眷顾
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见