24小时热门版块排行榜    

查看: 1350  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工科 320总分 求调剂 +7 梨花珞晚风 2026-03-17 7/350 2026-03-22 10:03 by 2026paper
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 化学工程321分求调剂 +18 大米饭! 2026-03-15 22/1100 2026-03-21 20:20 by HH领袖
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 307求调剂 +10 冷笙123 2026-03-17 10/500 2026-03-21 01:54 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 一志愿中南化学(0703)总分337求调剂 +8 niko- 2026-03-19 9/450 2026-03-20 21:57 by luoyongfeng
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 一志愿华中农业071010,总分320求调剂 +3 困困困困坤坤 2026-03-20 3/150 2026-03-20 20:38 by 学员8dgXkO
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
信息提示
请填处理意见