应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 其他 » 【求助】相关函数曲线异常
91/1返回列表
查看: 1325  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

astringent

铜虫 (著名写手)

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

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


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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-12-08 16:00:25
已阅   回复此楼   关注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


可能平衡不够或者时间太短
赞一下(2人) 回复此楼
技术博客：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:

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

你这个图形尾巴的地方应该是采样太少引起的大涨落所致。看看你的算法是不是有问题
赞一下(2人) 回复此楼
技术博客：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 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料化工调剂 +11 今夏不夏 2026-03-01 12/600 2026-03-01 22:33 by 刘兵
[考研] 材料学硕318求调剂 +3 February_Feb 2026-03-01 3/150 2026-03-01 22:21 by gaoxiaoniuma
[考研] 材料类求调剂 +10 wana_kiko 2026-02-28 12/600 2026-03-01 22:10 by 海嵙Y
[考研] 26考研报考西工大材料308分求调剂 +3 weizhong123 2026-03-01 3/150 2026-03-01 21:42 by 公瑾逍遥
[考研] 材料学硕318求调剂 +9 February_Feb 2026-03-01 11/550 2026-03-01 19:47 by 无懈可击111
[考研] 272求调剂 +6 材紫有化 2026-02-28 6/300 2026-03-01 18:58 by 18137688336
[考研] 291分工科求调剂 +9 science饿饿 2026-03-01 10/500 2026-03-01 18:55 by 18137688336
[考研] 材料学调剂 +9 提神豆沙包 2026-02-28 11/550 2026-03-01 18:15 by ms629
[考博] 26申博 +4 想申博！ 2026-02-26 6/300 2026-03-01 17:32 by 想申博！
[考研] 化工专硕348，一志愿985求调剂 +5 弗格个 2026-02-28 8/400 2026-03-01 17:25 by sunny81
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 313求调剂 +3 水流年lc 2026-02-28 3/150 2026-03-01 16:01 by 新能源达人
[考研] 求调剂 +6 repeatt?t 2026-02-28 6/300 2026-03-01 14:37 by Sakura绘
[考研] 材料284求调剂，一志愿郑州大学英一数二专硕 +10 想上岸的土拨鼠 2026-02-28 10/500 2026-03-01 14:12 by yc258
[硕博家园] 2025届双非化工硕士毕业，申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[基金申请] 面上模板改不了页边距吧？ +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
信息提示
关闭
请填处理意见
关闭