| 查看: 1718 | 回复: 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人回复
有多少人是今天查系统知道结果的?
已经有19人回复
面上合作单位盖章
已经有9人回复
面上函评意见出来了,像什么等级?
已经有17人回复
为什么国自然不能直接公布
已经有5人回复
学科评审组评审是指会评吗?
已经有4人回复
为什么资助数各大高校都创新高,自己申请怎么就这么难
已经有14人回复
小白求助 投论文要求的highlights应该如何写
已经有4人回复
麻烦专家们看看评委们的意见(F口面上)
已经有12人回复
怎么看青基中了没有啊
已经有6人回复
» 本主题相关价值贴推荐,对您同样有帮助:
基于L-M算法的BP神经网络
已经有3人回复
请教一个贝叶斯推理问题
已经有10人回复
【求助】相关函数归一化后,输出值怎么变成NaN了??
已经有4人回复
【zhou2009个人文集】关于HF、H2O的电荷计算
已经有15人回复
请问origin里面怎样自定义归一化的函数?
已经有3人回复
【求助】场及量子化
已经有13人回复

6楼2012-05-02 21:37:00
老虎大王
木虫 (著名写手)
- 模拟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










回复此楼
投票:
5