| 查看: 5437 | 回复: 95 | |||||||
| 【奖励】 本帖被评价85次，作者freewain增加金币 66.8 个 | |||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||
[资源]
提供一个 Birch-Murnaghan 拟合的算法（更新1）
|
|||||||
|
提供一个计算 Birch-Murnaghan EOS 的思路和 MatLab 算法，请大家指正。更新： 最近发现，提供的这个算法并不太好，虽然原理上没有问题，但是在数值算法中有时会产生误差。造成这个的原因在于中间涉及的一些数太小，或者说需要的精度太高。当我试图用 C 或者 Fortran 实现这种算法时，发现连 double 精度都无法保证正确结果。 现在看来，最好的方法还是用 least squares 方法直接拟合 Birch-Murnaghan 方程。我最近会试着用 Fortran 写一下这个功能，如果成功，会更新这个帖子放到这里。 很抱歉之前的不成熟的算法耽误大家的时间了。 另外，这里有一个拟合 Murnaghan 方程的 Fortran 代码，作者是 SIESTA 的维护者之一 Postnikov。http://www.home.uni-osnabrueck.de/apostnik/Software/mur_fit.f [ Last edited by freewain on 2013-4-22 at 21:21 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Birch-Murnaghan Regression.pdf
2012-03-06 05:36:43, 94.58 K
» 收录本帖的淘帖专辑推荐
 @个人收集@ | 第一性原理计算-tg | Photochemistry | 分子动力学 |
| 表界面计算 |
» 猜你喜欢
职称评审没过，求安慰
已经有49人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
三无产品还有机会吗
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- XPS谱中，用284.8校正峰位是不是拟合的C峰中一定有一个是284.8的？我迷茫了……
已经有6人回复
- 请问一个有理函数的曲线拟合过程如何做，给出数据，请求拟合，
已经有14人回复
- 晶格常数优化
已经有7人回复
- 请教一个数据非线性拟合的问题
已经有11人回复
- 【公式拟合问题】如何拟合下面这样的一个公式，多谢大牛。
已经有6人回复
- 用origin如何在一个图中对多组连续的点进行线性拟合呢？
已经有5人回复
- 【求助】E-V曲线晶体状态方程拟合问题
已经有8人回复
- 【求助】帮忙拟合一下
已经有7人回复
- 【一个数据拟合的问题，比较麻烦】
已经有17人回复
- 【求助】拟合晶格常数
已经有14人回复
- 【求助】我过著名的药物动力学分析软件 DAS 2.0 数据拟合部分用的是什么算法？
已经有7人回复
- 【求助】如何使得matlab拟合的曲线强制经过一个点
已经有8人回复
- 【原创】elk中的状态方程拟合工具eos
已经有128人回复
16楼2012-03-13 09:23:00
12楼2012-03-07 21:48:00
13楼2012-03-08 03:21:48
简单回复
2012-03-06 09:19
回复
五星好评 顶一下，感谢分享！
2012-03-06 09:43
回复
五星好评 顶一下，感谢分享！
cheerian5楼
2012-03-06 11:00
回复
五星好评 
乐石||6楼
2012-03-06 16:45
回复
五星好评 顶一下，感谢分享！
nkhjh7楼
2012-03-07 09:25
回复
五星好评 顶一下，感谢分享！
zhxd10288楼
2012-03-07 14:08
回复
五星好评 顶一下，感谢分享！
thinkllh9楼
2012-03-07 15:32
回复
五星好评 顶一下，感谢分享！
luxia021010楼
2012-03-07 15:50
回复
五星好评 顶一下，感谢分享！
tonyngch11楼
2012-03-07 16:14
回复
三星好评 顶一下，感谢分享！
2012-03-12 16:53
回复
五星好评 顶一下，感谢分享！