| 查看: 2004 | 回复: 4 | |||
astringent铜虫 (著名写手)
|
[交流]
【求助】相关函数归一化后，输出值怎么变成NaN了？？ 已有3人参与
|
|
我编程计算相关函数时，未归一化时求出的相关函数趋势是对的，可是我加了归一化的语句后（也就是程序中的tcf0(delt)=tcf0(delt)+tcf0(t)*tcf0(t)），输出的相关函数的值就变成NaN了，我查了一下，NaN可能是运算错误，请大家帮忙看看，是我的归一化条件加的不对吗？？先谢谢了 下面是我的程序： 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)/tcf0(ical) [ Last edited by ghcacj on 2010-12-8 at 09:28 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 波函数----并非统计工具而是物理真实 【转】
已经有40人回复
- 如何将取向矩阵中归一化的{hkl}转化成互质化的{HKL}，求算法！！！！
已经有10人回复
- 求助 如何确定物质溶解时的准确PH值呢 等电点怎么测呢
已经有14人回复
- 齐全的药物辅料 值得收藏
已经有24人回复
- 港科大 IC Design 授课式的硕士(MSc)，值不值得去读呢？恳求大家的意见
已经有59人回复
- 大家推荐一下实变函数、泛函分析好的教材与视频
已经有21人回复
- origin进行非线性拟合如何选择函数或者自定义合适的函数
已经有19人回复
- 6-31G基组添加极化和弥散函数
已经有14人回复
- 怎么得到旋轨耦合中不同J值的本征值和波函数？
已经有6人回复
- 为什么画出的对关联函数第三峰刚开始,之后的函数值就全变为零了,怎么解决这个问题
已经有5人回复
- 【求助】相关函数归一化问题
已经有6人回复
- 请问origin里面怎样自定义归一化的函数？
已经有3人回复
qphll
金虫 (正式写手)
- 模拟EPI: 10
- 应助: 18 (小学生)
- 金币: 2708.7
- 散金: 2294
- 红花: 8
- 帖子: 554
- 在线: 227.1小时
- 虫号: 89654
- 注册: 2005-08-29
- 性别: GG
- 专业: 碳素材料与超硬材料
2楼2010-12-07 04:56:09
astringent
铜虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 2210.7
- 散金: 1233
- 红花: 4
- 帖子: 1180
- 在线: 354.8小时
- 虫号: 900839
- 注册: 2009-11-12
- 专业: 理论和计算化学
★
zh1987hs(金币+1):thank you 2010-12-07 10:31:05
zh1987hs(金币+1):thank you 2010-12-07 10:31:05
|
我原来的程序是这样的，此时没出现输出错误的问题， 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) enddo endif enddo enddo do ical=1,npoints tcf(ical)=tcf(ical)/nparticle*ntime(t) enddo 但是我需要算归一化的相关函数，于是我就把tcf(delt)=tcf(delt)+tcf(ntel)*tcf0(t) 改成了tcf(delt)=tcf(delt)+（tcf(ntel)*tcf0(t)）/(tcf0(t)*tcf0(t) ) ,就出现输出错误了，是怎么回事呢？ |
3楼2010-12-07 05:01:33
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
4楼2010-12-07 10:44:19
astringent
铜虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 2210.7
- 散金: 1233
- 红花: 4
- 帖子: 1180
- 在线: 354.8小时
- 虫号: 900839
- 注册: 2009-11-12
- 专业: 理论和计算化学
5楼2010-12-08 06:31:37
40