24小时热门版块排行榜

返回列表

tiechong

金虫 (著名写手)

应助: 9 (幼儿园)
金币: 493.7
帖子: 1094
在线: 391.7小时
虫号: 793352

[交流] 【求助】科研需要，请人编两个关于一维radial distribution function matlab小程序

由于研究需要，想请人编两个小程序，用matlab来模拟一维的粒子运动，得到radial distribution function（径向分布函数）。
两个程序分别是在Hard Sphere system下，以及Lennard-Jones system下。
因为，只需要考虑一维空间，既直线运动，所以不会有太大工作量，大概也就几十行代码。

因为我是做纯有机化学的，自己以前从来没有用过matlab，这两天一边看说明，一边编程，编了在hard sphere情况下的程序，结果运行总是出错。我一时半活儿也没找到原因。

希望懂这方面的人能够和我联系，因为我比较着急要。另外，有酬劳。谢谢。

fantastli@hotmail.com， msn联系

1）需要用matlab编译一个程序，来从molecular dynamic上模拟1维的hard sphere system，并且用这个程序来得到radial distribution function，并画图。
2）需要用matlab编译一个程序，来从molecular dynamic上模拟1维的Lennard-Jones system，并且用这个程序来得到radial distribution function，并画图。
要求如下（因为我不是搞分子模拟的，如果你看不懂我写的要求，可以再问我）：

1.我需要程序的code
2.如果可以的话，请写一下简明的构想程序的思路，或者在关键语句前作一下注释
3.这个程序可以生成那种小球之间对撞的动态图形，以及最终得到的方程的曲线图
4.我希望能够尽快拿到，最好在一周以内

对于第一个模型，我已经写好了大部分code，但是运行起来总是提示问题。因为我没学过matlab，所以是边学边做，暂时不知道问题出在哪里。如果你需要，我可以把我的code发给你，你如果能帮我修正好，也可以。

其实这两个程序不是很难。只要构思1维空间就可以。你可以设想有一排小球，只能?左右动。其中第一个和最后一个球是不动的。中间的球的速度是随机的，然后相邻的小球之间会互相撞击，从而改变方向和运动速度。我希望能够得到这些小球互相撞击的演示模型以及生成的动态图），另外需要得到径向分布函数的曲线图。两个模型不同的是，第一个是只有有接触才产生作用力，所以作用力是不连续的（根据文献资料，好像可以把“获得径分布方程”的问题转换成“获得最短碰撞时间”，这个我有些迷惑）；第二个模型小球之间的作用力大小（或势能）是和小球之间的距离有关系，所以作用力是连续的。

另外，我也有一些用C++语言编写类似程序的资料，如果你需要我也可以给你。

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-12-03 12:33:52

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

lijinfeng042

木虫 (小有名气)

仿真EPI: 2
应助: 1 (幼儿园)
金币: 2156.1
帖子: 291
在线: 31.5小时
虫号: 1019062

tiechong(金币+5): 2011-03-05 11:36:34

引用回帖:

Originally posted by tiechong at 2010-12-03 12:33:52:
由于研究需要，想请人编两个小程序，用matlab来模拟一维的粒子运动，得到radial distribution function（径向分布函数）。
两个程序分别是在Hard Sphere system下，以及Lennard-Jones system下。
因为，只需要考 ...

专业性太强一般很难找到合适的

赞一下