| 查看: 1204 | 回复: 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.是否数据本身的问题，不可以实现整个温度段拟合，分开成两个温度范围的拟合？ 十分感谢  捕获.PNG |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有124人回复
-
物理学I论文润色/翻译怎么收费?
已经有214人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
fairless
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1.5
- 散金: 10
- 帖子: 10
- 在线: 26.6小时
- 虫号: 2309468
- 注册: 2013-03-01
- 专业: 凝聚态物性 II ：电子结构
5楼2013-12-08 19:27:31
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2013-12-08 13:23:08
fairless
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1.5
- 散金: 10
- 帖子: 10
- 在线: 26.6小时
- 虫号: 2309468
- 注册: 2013-03-01
- 专业: 凝聚态物性 II ：电子结构
3楼2013-12-08 13:48:28
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★
fairless: 金币+5, ★★★★★最佳答案 2013-12-08 19:27:46
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
200
