| 查看: 1406 | 回复: 6 | |||
astringent铜虫 (著名写手)
|
[交流]
【求助】相关函数归一化问题 已有4人参与
|
|
最近一直在看自相关函数的程序，有个问题想不通。 自相关函数的定义是A=(a(t)*a(0))/(a(0)*a(0)),当对体系中所有粒子求和时，是先算出每个例子的(a(t)*a(0)) /(a(0)*a(0))然后求和？还是先分别求a(t)*a(0)以及a(0)*a(0)和然后再相除呢？两者有什么区别？请大家指点。 [ Last edited by ghcacj on 2010-12-8 at 09:27 ] |
回复此楼» 猜你喜欢
本人考085602 化学工程 专硕
已经有9人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有7人回复
-
326求调剂
已经有3人回复
-
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 基于L-M算法的BP神经网络
已经有3人回复
- 请教一个贝叶斯推理问题
已经有10人回复
- 【求助】相关函数归一化后，输出值怎么变成NaN了？？
已经有4人回复
- 【zhou2009个人文集】关于HF、H2O的电荷计算
已经有15人回复
- 请问origin里面怎样自定义归一化的函数？
已经有3人回复
- 【求助】场及量子化
已经有13人回复
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2010-12-07 10:59:29
astringent
铜虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 2210.7
- 散金: 1233
- 红花: 4
- 帖子: 1180
- 在线: 354.8小时
- 虫号: 900839
- 注册: 2009-11-12
- 专业: 理论和计算化学
|
请问老虎大王，我在写程序的时候直接把t=0的相关函数做分母了，也就是 tcf(delt)=tcf(delt)+tcf(ntel)*tcf0(t)/ tcf0(t)*tcf0(t)，可是输出错误，请问，我这样写为什么不可以呢？应该怎样写？ 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(t)*tcf0(t) enddo endif enddo enddo do ical=1,npoints tcf(ical)=tcf(ical)/(npoints*ntime(t)) write(*,*) ical, tcf(ical) enddo |
3楼2010-12-07 15:57:21
astringent
铜虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 2210.7
- 散金: 1233
- 红花: 4
- 帖子: 1180
- 在线: 354.8小时
- 虫号: 900839
- 注册: 2009-11-12
- 专业: 理论和计算化学
4楼2010-12-08 06:34:00
astringent
铜虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 2210.7
- 散金: 1233
- 红花: 4
- 帖子: 1180
- 在线: 354.8小时
- 虫号: 900839
- 注册: 2009-11-12
- 专业: 理论和计算化学
|
为什么我用两种方法得到的结果有很大差别呢，劳烦指教。 下面是我的程序 这是我的程序： 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 if(delt.lt.npoints) tcf(delt)=tcf(delt)+tcf(ntel)*tcf0(t) tcf0(delt)=tcf0(delt)+tcf0(t)*tcf0(t) acf(delt)=tcf(delt)+tcf(ntel)*tcf0(t)/ tcf0(t)*tcf0(t) endif enddo endif enddo enddo do ical=1,npoints tcf(ical)=(tcf(ical)/(npoints*ntime(t))/(tcf0(ical)/(npoints*ntime(t)) acf(ical)=acf(ical)/(npoints*ntime(t)) write(*,*) ical, tcf(ical) write(*,*) ical, acf(ical) enddo 下面是我的输出结果（我只截取了一部分） 0.0000000000000000 1.0000000000000000 0.0000000000000000 1.14207400639561443E-004 1.0000000149011612 0.67655786350148372 1.0000000149011612 7.72679149727596715E-005 2.0000000298023224 0.62686567164179108 2.0000000298023224 7.15926989083817998E-005 3.0000000447034836 0.60360360360360366 3.0000000447034836 6.89359985842397961E-005 4.0000000596046448 0.58610271903323263 4.0000000596046448 6.69372680485647080E-005 5.0000000745058060 0.57012195121951215 5.0000000745058060 6.51121460963353292E-005 6.0000000894069672 0.55828220858895705 6.0000000894069672 6.37599598662582241E-005 7.0000001043081284 0.54629629629629639 7.0000001043081284 6.23910799790196834E-005 8.0000001192092896 0.51552795031055898 8.0000001192092896 5.88771071620099369E-005 9.0000001341104507 0.52187499999999998 9.0000001341104507 5.96019872087711250E-005 10.000000149011612 0.52515723270440251 10.000000149011612 5.99768424742350986E-005 |
5楼2010-12-10 00:17:07
6楼2012-05-02 21:37:00
manutdhk
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 5 (幼儿园)
- 金币: 5939.4
- 散金: 70
- 红花: 5
- 帖子: 94
- 在线: 97.3小时
- 虫号: 983544
- 注册: 2010-03-27
- 专业: 核技术及其应用
7楼2012-05-02 22:38:17
