| 查看: 3442 | 回复: 11 | ||||
tandz木虫 (正式写手)
|
[求助]
如何拟合morse potential的参数?
|
|||
| 结合能E与距离r 已知. 请教各位，如何写个script拟合morse potential E(r) = D*(1-exp(a*(r-rm))^2的参数 D, a, rm呢? 本人几乎没有编程水平，请大家帮忙，非常感谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 精华 |
» 猜你喜欢
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有130人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有3人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68137.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
tandz: 金币+50, ★★★很有帮助, 太感谢了！ 2013-04-04 17:03:15
tandz: 金币+20, ★★★★★最佳答案 2013-04-06 17:18:56
csgt0: 金币+2, 谢谢 2013-04-07 15:11:49
感谢参与，应助指数 +1
tandz: 金币+50, ★★★很有帮助, 太感谢了！ 2013-04-04 17:03:15
tandz: 金币+20, ★★★★★最佳答案 2013-04-06 17:18:56
csgt0: 金币+2, 谢谢 2013-04-07 15:11:49
|
MATLAB代码可以吗？ 假如已经获得了一组r，E的数据： r：1 2 3 9 E：1 5 8 18 MATLAB程序如下，结果见附图1。 function feixianxingnihe clear all;clc format long tspan=[1 2 3 9]; %r的数据，在此输入 xexp=[1 5 8 18]; %E的数据，在此输入 k0=[1 1 1]; lb=-[1 1 1]*1e5; ub=[1 1 1]*1e5; %------------------------------------------------------------------------- % 使用函数lsqnonlin()进行参数估计 OPTIONS=optimset('MaxFunEvals',1000); [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,OPTIONS,tspan,xexp); ci = nlparci(k,residual,jacobian); %residual; fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n') fprintf('\n\t参数 D = %.16f',k(1)) fprintf('\n\t参数 a = %.16f',k(2)) fprintf('\n\t参数rm = %.16f',k(3)) y=KineticsEqs(tspan,k); R2=1-sum((xexp-y).^2)./sum((xexp-mean(y)).^2); fprintf('\n\t相关系数之平方R^2 = %.16f',R2); figure plot(tspan,KineticsEqs(tspan,k),'b',tspan,xexp,'or'),legend('计算值','实验值','Location','Best') %------------------------------------------------------------------------- function f = ObjFunc(k,tspan,xexp) f=KineticsEqs(tspan,k)-xexp; %------------------------------------------------------------------------ function xt = KineticsEqs(t,k) xt=k(1)*(1-exp(k(2).*(t-k(3)))).^2; 做方程拟合，1stopt代码简单易懂，1stopt代码如下： Parameters D,a,rm; Variable r,E; Function E=D*(1-exp(a*(r-rm)))^2; Data; 1 1 2 5 3 8 9 18 计算结果见附图2。 可以看到两种方法计算结果相当接近。  附图1.jpg  附图2.jpg |
5楼2013-04-04 16:27:04
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2013-04-03 10:15:16
tandz
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1780.2
- 散金: 38
- 红花: 1
- 帖子: 449
- 在线: 703.8小时
- 虫号: 636814
- 注册: 2008-10-26
- 专业: 理论和计算化学
3楼2013-04-03 10:25:28
tandz
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1780.2
- 散金: 38
- 红花: 1
- 帖子: 449
- 在线: 703.8小时
- 虫号: 636814
- 注册: 2008-10-26
- 专业: 理论和计算化学
4楼2013-04-03 19:21:49
tandz
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1780.2
- 散金: 38
- 红花: 1
- 帖子: 449
- 在线: 703.8小时
- 虫号: 636814
- 注册: 2008-10-26
- 专业: 理论和计算化学
6楼2013-04-04 17:06:25
tandz
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1780.2
- 散金: 38
- 红花: 1
- 帖子: 449
- 在线: 703.8小时
- 虫号: 636814
- 注册: 2008-10-26
- 专业: 理论和计算化学
7楼2013-04-06 16:06:21
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68137.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
tandz8楼2013-04-06 18:08:39
tandz
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1780.2
- 散金: 38
- 红花: 1
- 帖子: 449
- 在线: 703.8小时
- 虫号: 636814
- 注册: 2008-10-26
- 专业: 理论和计算化学
9楼2013-04-07 09:58:22
10楼2013-10-11 14:37:54
5