24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1389  |  回复: 8

astringent

铜虫 (著名写手)

[交流] 【求助】相关函数曲线异常 已有5人参与

我做了相关函数的曲线,可是并不像书上所画的,开始时为1,长时间接近0。按我的数据做的图为蓝色部分,可是一直也没有接近0,横轴的时间是以ps为单位的。粉色的线是我用exp(-a/t)得到的。请大家帮我分析一下,为什么会这样呢?是我的模拟时间不够还是怎么回事?多谢了。


[ Last edited by ghcacj on 2010-12-8 at 16:39 ]
回复此楼

» 猜你喜欢

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

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

huai

木虫 (职业作家)

什么量的相关函数?
2楼2010-12-08 16:23:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

astringent

铜虫 (著名写手)

引用回帖:
Originally posted by huai at 2010-12-08 16:23:47:
什么量的相关函数?

hydrogen bond
3楼2010-12-08 17:34:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+3):谢谢 2010-12-09 14:31:03
如果你算得是纯水的相关函数的话,那么你的算法有可能是错的

这是我算得一个相关函数




这是个文献上的nature,379-4


可能平衡不够或者时间太短
技术博客:http://hi.baidu.com/coolrainbow/blog
4楼2010-12-08 19:18:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

astringent

铜虫 (著名写手)

引用回帖:
Originally posted by coolrainbow at 2010-12-08 19:18:30:
如果你算得是纯水的相关函数的话,那么你的算法有可能是错的

这是我算得一个相关函数


[img]http://pic.muchong.com/201012/8/239560_191517.jp ...

我算不是纯水,是氨基酸周围的水氢键的相关函数。
你说的平衡不够或者时间太短,是说我抽样阶段的模拟时间短吗?
从我的图来看,一点接近零的趋势都没有啊。
5楼2010-12-09 03:06:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+2): 谢谢 2010-12-09 14:30:54
引用回帖:
Originally posted by astringent at 2010-12-09 03:06:27:

我算不是纯水,是氨基酸周围的水氢键的相关函数。
你说的平衡不够或者时间太短,是说我抽样阶段的模拟时间短吗?
从我的图来看,一点接近零的趋势都没有啊。

你这个图形尾巴的地方应该是采样太少引起的大涨落所致。看看你的算法是不是有问题
技术博客:http://hi.baidu.com/coolrainbow/blog
6楼2010-12-09 11:36:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

astringent

铜虫 (著名写手)

引用回帖:
Originally posted by coolrainbow at 2010-12-09 11:36:09:



你这个图形尾巴的地方应该是采样太少引起的大涨落所致。看看你的算法是不是有问题

我的采样确实有点问题。在采样阶段,老板只让我做了100ps,每隔100fs取一次样,所以我只有1000个点,所以采样到到最后时,样本就很少。至于算法,我是这样写的:
do i=1,10000 (粒子数)
     do ntel=1,npoints (共有npoints个数据)
          if (确定timeorigin的条件) then
             t0=t0+1
             time0(t0)=ntel
             tcf0(t0)=tcf(ntel)         
             do t=1,t0
                  delt=ntel-time0(t)+1
                  tcf(delt)=tcf(delt)+tcf(ntel)*tcf0(t)  
                  tcf0(delt)=tcf0(delt)+tcf0(t)*tcf0(t)
             enddo
           endif
    enddo
enddo
do ical=1,npoints
tcf(ical)=(tcf(ical)/nparticle*ntime(ical)/(tcf0(ical) /nparticle*ntime(ical)
enddo
有劳你帮我看看,这个算法有没有问题,多谢。
7楼2010-12-09 15:55:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

葵sunny

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
不知道楼主可否告诉一下相关函数怎么算呢,怎么出来上面类似的图呢,谢谢啦
8楼2012-10-27 10:12:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rusong231

铜虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
最简单的方法直接用me软件计算就行了,然后将时间步长取短一些,周期性晶胞取大些时候计算时间取长点,不一定

[ 发自手机版 http://muchong.com/3g ]
9楼2012-10-27 10:26:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 astringent 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 272分材料子求调剂 +37 Loy0361 2026-04-10 48/2400 2026-04-12 10:28 by lhj2009
[考研] 电子信息270求调剂 +15 terminal469 2026-04-07 15/750 2026-04-12 09:44 by 逆水乘风
[考研] 347求调剂 +4 mhyqyy 2026-04-06 4/200 2026-04-12 02:27 by 秋豆菜芽
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 339求调剂 +7 hanwudada 2026-04-11 7/350 2026-04-11 23:06 by 852137818
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 296求调剂 +14 汪!?! 2026-04-08 15/750 2026-04-11 20:28 by dongdian1
[考研] 本人女孩 +7 吼吼, 2026-04-10 9/450 2026-04-11 14:45 by ACS Nano——
[考研] 297求调剂 +9 Kwgyz 2026-04-09 9/450 2026-04-11 10:09 by zhq0425
[考研] 346,工科0854求调剂,专硕 +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by 猪会飞
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 一志愿华南理工大学331分材料求调剂 +9 天下ww 2026-04-09 9/450 2026-04-10 22:58 by Ftglcn90
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[硕博家园] 0856材料化工求调剂,一志愿211,初试成绩349 +5 江淮北月 2026-04-05 5/250 2026-04-10 16:26 by 高维春
[考研] 生物与医药273求调剂 +18 荔题南墙 2026-04-05 19/950 2026-04-10 08:14 by kangsm
[考研] 一志愿中国科学院上海有机所,有机化学356分找调剂 +11 Nadiums 2026-04-09 11/550 2026-04-09 18:04 by lijunpoly
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +3 慕绝cc 2026-04-09 3/150 2026-04-09 09:57 by liuhuiying09
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +11 Mr. Z 2026-04-05 12/600 2026-04-08 16:15 by luoyongfeng
[考研] 336求调剂,一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
信息提示
请填处理意见