| 查看: 1620 | 回复: 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 ] |
回复此楼» 猜你喜欢
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 波函数----并非统计工具而是物理真实 【转】
已经有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
20