24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 919  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 266求调剂 +4 08电气工程 2026-04-03 4/200 2026-04-08 15:57 by 向9755
[考研] 0702物理学学硕299求调剂 +4 祁柒连 2026-04-06 4/200 2026-04-08 13:56 by wutongshun
[考研] 327求调剂 +11 Xxjc1107. 2026-04-06 11/550 2026-04-08 12:37 by yongzhesheng
[考研] 农学,求调剂,314分 +4 访客记录可爱 2026-04-04 4/200 2026-04-07 21:07 by 等岸
[基金申请] 请问共同通讯和共同一作的认可度问题 10+5 psa1234 2026-04-01 11/550 2026-04-07 16:57 by caiqing
[考研] 生物调剂 +5 橙子橙子橙子啊 2026-04-05 9/450 2026-04-07 15:31 by 上岸快快
[考研] 求调剂 +4 电气小神童 2026-04-04 6/300 2026-04-07 00:14 by guanxin1001
[考研] 求调剂,一志愿郑州大学材料与化工专硕,英二数二342分,求老师收留 +19 v12abo 2026-04-02 21/1050 2026-04-06 09:29 by 蓝云思雨
[考研] 348求调剂 +3 车厘子zzz 2026-04-05 3/150 2026-04-05 20:30 by 啵啵啵0119
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 数一英一274机械调剂 +5 星陨流霞 2026-04-04 6/300 2026-04-05 11:38 by arrow8852
[考研] 283分求调剂 +7 小聂爱学习 2026-04-03 7/350 2026-04-04 21:51 by hemengdong
[考研] 296材料专硕求调剂 +21 202451007219 2026-04-02 22/1100 2026-04-04 21:48 by hemengdong
[考研] 297求调剂 +11 ljy20040718! 2026-04-03 13/650 2026-04-04 09:23 by 来看流星雨10
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 330求调剂 +3 白神呜呼呼 2026-04-02 3/150 2026-04-03 10:15 by 蓝云思雨
[考研] 318求调剂,计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 材料专硕322分 +11 哈哈哈吼吼吼哈 2026-04-01 11/550 2026-04-02 10:52 by lnilvy
[考研] 377求调剂 +3 RASKIN 2026-04-02 3/150 2026-04-02 09:45 by zzchen2000
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
信息提示
请填处理意见