| 查看: 853 | 回复: 1 | |||
[交流]
【求助】科研需要,请人编两个关于一维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++语言编写类似程序的资料,如果你需要我也可以给你。 |
» 猜你喜欢
售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急
已经有5人回复
有谁可曾问过你过的还好吗?
已经有21人回复
E0414, 我的本子有没有希望?
已经有7人回复
一篇论文同时出现在两个期刊,一模一样,这算不算学术不端,请各位老师斧正。
已经有12人回复
希望面上有个好结果
已经有7人回复
三区计算机方向期刊推荐
已经有5人回复
sci论文二审求助
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 抢金币啦!回帖就可以得到:
多聚组氨酸,大分子量
+1/80
坐标北京,诚征男友,找到会删帖。试试运气!说不定就有了!
+1/80
骨生物材料与侗药调控类器官再生湖南省普通高等学校重点实验室招聘生物医用材料相关博
+2/68
山东省自然科学基金
+1/57
国家青年基金祈福
+1/49
哈工大深圳-材料学院-招收申请审核制2027年春季/秋季入学博士生(2026年9月报名)
+1/36
博士毕业之后申请发明专利
+1/28
湖南师范大学(211)—招收2026年入学的“申请-考核”制博士(有微纳光学经验者优先)
+1/24
上海大学微电子学院杨军教授团队招聘带编专任教师
+1/10
哈尔滨工业大学招收硕士研究生(欢迎环境、市政、生物、化学、农业等专业,长期有效)
+1/10
东北师范大学荒漠与草地生态学方向诚聘博士后青年才俊
+1/10
墨尔本大学(26年QS19)招全奖博士/CSC博士(补齐全奖)/访问学者等-材料/生物医学等
+1/9
密苏里大学生物材料合成生物学博士后招聘
+1/5
密苏里大学生物材料合成生物学博士后招聘
+1/4
近红外二区荧光杂峰
+1/4
哈尔滨工业大学(深圳)赵怡潞课题组诚招博士后、2027学年博士生
+1/3
上海大学微电子学院杨军教授团队招聘带编专任教师
+1/2
招科研助理,提前招27年博士生
+1/1
2026年黄河科技学院纳米功能材料研究所招聘
+1/1
中科院深圳先进技术研究院成会明院士/唐永炳杰青团队博士后招聘
+1/1
2楼2010-12-04 22:57:34











回复此楼