24小时热门版块排行榜

返回列表

astringent

铜虫 (著名写手)

应助: 16 (小学生)
金币: 2210.7
散金: 1233
红花: 4
帖子: 1180
在线: 354.8小时
虫号: 900839
注册: 2009-11-12
专业: 理论和计算化学

[交流] 【求助】相关函数归一化后，输出值怎么变成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 ]

回复此楼

» 猜你喜欢

求助:我三月中下旬出站，青基依托单位怎么办？已经有8人回复
Cas 72-43-5需要30g，定制合成，能接单的留言已经有8人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
磺酰氟产物，毕不了业了！已经有5人回复
论文终于录用啦！满足毕业条件了已经有25人回复
2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有3人回复
自荐读博已经有3人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

波函数----并非统计工具而是物理真实【转】已经有40人回复
如何将取向矩阵中归一化的{hkl}转化成互质化的{HKL}，求算法！！！！已经有10人回复
求助如何确定物质溶解时的准确PH值呢等电点怎么测呢已经有14人回复
齐全的药物辅料值得收藏已经有24人回复
港科大 IC Design 授课式的硕士(MSc)，值不值得去读呢？恳求大家的意见已经有59人回复
大家推荐一下实变函数、泛函分析好的教材与视频已经有21人回复
origin进行非线性拟合如何选择函数或者自定义合适的函数已经有19人回复
6-31G基组添加极化和弥散函数已经有14人回复
怎么得到旋轨耦合中不同J值的本征值和波函数？已经有6人回复
为什么画出的对关联函数第三峰刚开始,之后的函数值就全变为零了,怎么解决这个问题已经有5人回复
【求助】相关函数归一化问题已经有6人回复
请问origin里面怎样自定义归一化的函数？已经有3人回复

1楼 2010-12-07 03:59:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

qphll

金虫 (正式写手)

模拟EPI: 10
应助: 18 (小学生)
金币: 2708.7
散金: 2294
红花: 8
帖子: 554
在线: 227.1小时
虫号: 89654
注册: 2005-08-29
性别: GG
专业: 碳素材料与超硬材料

★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+3):thank you 2010-12-07 10:30:54
zh1987hs(金币+3):thank you 2010-12-07 10:31:09

引用回帖:

Originally posted by astringent at 2010-12-07 03:59:01:
我编程计算相关函数时，未归一化时求出的相关函数趋势是对的，可是我加了归一化的语句后（也就是程序中的tcf0(delt)=tcf0(delt)+tcf0(t)*tcf0(t)），输出的相关函数的值就变成NaN了，我查了一下，NaN可能是运算错 ...

一般而言, 如果出现NaN, 那么检查你的输入文件和你程序中的read, 循环语句是不是匹配. 大多数, 这是文件格式的问题.

赞一下(3人)

回复此楼

Life, Love, Laugh.

2楼2010-12-07 04:56:09

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

astringent

铜虫 (著名写手)

应助: 16 (小学生)
金币: 2210.7
散金: 1233
红花: 4
帖子: 1180
在线: 354.8小时
虫号: 900839
注册: 2009-11-12
专业: 理论和计算化学

★
zh1987hs(金币+1):thank you 2010-12-07 10:31:05

引用回帖:

Originally posted by qphll at 2010-12-07 04:56:09:

一般而言, 如果出现NaN, 那么检查你的输入文件和你程序中的read, 循环语句是不是匹配. 大多数, 这是文件格式的问题.

我原来的程序是这样的，此时没出现输出错误的问题，
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) ) ,就出现输出错误了，是怎么回事呢？

赞一下(1人)

回复此楼

3楼2010-12-07 05:01:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖