| 查看: 1698 | 回复: 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 ] |
回复此楼» 猜你喜欢
范进中举一文的中心思想
已经有8人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有6人回复
-
我面上完蛋了
已经有13人回复
-
只有每年这种时候来逛逛小木虫
已经有27人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
人气不行了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 基于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
20