| 查看: 899 | 回复: 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++语言编写类似程序的资料,如果你需要我也可以给你。 |
» 猜你喜欢
2026年国自然面上资助率
已经有11人回复
微信指数没变化,科研之友没阅读
已经有18人回复
基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊
已经有8人回复
HXDI做水性聚氨酯乳液,是不是特别容易出渣
已经有3人回复
小木虫看见有人已经知道结果了
已经有19人回复
系统今天又提示维护了,估计离放榜不远了
已经有15人回复
申请2027年材料/化学类博士
已经有5人回复
你们的时间戳变了吗
已经有4人回复
今年的WR进展到哪一步了?
已经有9人回复
产物和副产物价值比较
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 抢金币啦!回帖就可以得到:
求助2013年-2011年,武汉哪所学校的实验室合成过单价3万/克(或毫升)的剧毒成品?
+5/2360
今年面上没有消息,想来是没了,散个金攒攒人品,为明年祈福
+1/608
助理研究员——吉林大学-苏州国家实验室光电半导体智能设计与应用团队
+1/479
坐标北京,92男,诚征结婚对象
+1/170
具有超高耐腐蚀性的铝合金阳极氧化后处理工艺
+1/82
化工学院邱惠斌教授联合同仁医院招聘博士后
+1/78
Postdoctoral fellow in Medicinal Chemistry, McGill University
+1/71
祈福面上好运,散金币
+1/48
中国特种设备检测研究院招聘1名高温涂层方向博士后
+2/34
坐标南京
+1/30
深圳大学 招收2027级推免研究生 (金属材料/ 3D打印/ 氢能等方向)
+1/24
西湖大学物理学、光学、电子信息方向博士生招生(长期有效)
+5/20
一株吊兰
+1/12
招生宣讲会倒计时!7月29日20:00,不见不散!
+1/10
四川大学周加境课题组招博士研究生/研究助理(自组装材料/生物医学方向)
+1/6
国家自然科学基金申报者发榜前心理状态调查,麻烦帮忙填一下哇!
+1/5
准现货/原装全新超微整机8x NVIDIA HGX A100 80G
+1/4
国际 EI 会议征稿进行时
+1/2
中科院大连化物所诚聘博士后/科研助理/项目聘用人员
+1/1
浙江大学王宇课题组诚聘博士后/博士研究生/科研助理(招聘长期有效)
+1/1
2楼2010-12-04 22:57:34










回复此楼