24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1385  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 335求调剂 +5 想上岸呀!! 2026-04-12 5/250 2026-04-12 07:41 by laoshidan
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 085400 328分 求调剂 +10 喂你一个大橙子 2026-04-09 14/700 2026-04-11 19:53 by lqspecial
[考研] 复试调剂 +9 积极向上; 2026-04-10 11/550 2026-04-11 09:25 by 猪会飞
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 22408 366分,本科211,一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 085404 298分求调剂 +10 呼啦呼啦呼呼呼 2026-04-10 11/550 2026-04-10 16:44 by wangy0907
[考研] 材料复试求调剂 +20 xhhdjdjsjks 2026-04-09 20/1000 2026-04-10 10:25 by 孙小小12457
[考研] 生物学调剂,一志愿西南大学348,Top期刊一区二作、二区三作,三等奖学金三次 +4 candyyyi 2026-04-09 4/200 2026-04-09 18:39 by l_paradox
[考研] 367求调剂 +10 hffQAQ 2026-04-09 10/500 2026-04-09 18:06 by lijunpoly
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 一志愿郑州大学085600求调剂 +21 吃的不少 2026-04-05 24/1200 2026-04-08 16:47 by sunhuadong
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +10 努力奋斗112 2026-04-07 10/500 2026-04-08 15:01 by screening
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 277求调剂 +4 考研调剂lxh 2026-04-06 6/300 2026-04-08 10:40 by 逆水乘风
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 277求调剂 +5 考研调剂lxh 2026-04-05 5/250 2026-04-05 19:03 by chy09050039
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
信息提示
请填处理意见