24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

q1036210267

木虫 (正式写手)

应助: 26 (小学生)
金币: 2182.2
散金: 806
红花: 4
帖子: 893
在线: 262.1小时
虫号: 1346635
注册: 2011-07-15
性别: GG
专业: 药物化学

[求助] 怎么用Fortran 去做曲线拟合？

曲线拟合的 fortran 程序怎么编？

例如图中：第一列是距离r，第二列是能量E。
已知距离和能量满足公式 E = a（（b/r) **9 - (b/r)**6 ）

怎么用程序拟合得到 a 和 b

谢谢，金币不够我再加[ 来自科研家族寂寞大虫家族 ]

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有41人回复
回收溶剂求助已经有7人回复
硝基苯如何除去已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
三无产品还有机会吗已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

生活中坚持快乐的原则不改变，学习工作中踏实上进！`~`

1楼 2011-11-17 09:45:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★
xzhdty(金币+2): 欢迎常来程序语言讨论 2011-11-17 20:42:44
q1036210267(金币+5): 2011-11-18 14:46:51

用1stOpt吧，比Origin、Matlab等强大、方便好用的多：

CODE:

Variable r,E;

Function E = a*((b/r)^9  - (b/r)^6);

Data;

2.89999999850988        47.5834608565225

   3.09999999850988       -91.0481880549658

   3.29999999850988       -100.890821284316

   3.49999999850988       -100.883832534262

....

均方差(RMSE): 17.3437930194149
残差平方和(SSE): 29479.10131743
相关系数(R): 0.878771772121621
相关系数之平方(R^2): 0.772239827477775

参数最佳估算
---------- -------------
a 906.43777208918
b 2.95729352643941