| 查看: 1349 | 回复: 8 | |||
astringent铜虫 (著名写手)
|
[交流]
【求助】相关函数曲线异常 已有5人参与
|
我做了相关函数的曲线，可是并不像书上所画的，开始时为1，长时间接近0。按我的数据做的图为蓝色部分，可是一直也没有接近0，横轴的时间是以ps为单位的。粉色的线是我用exp（-a/t）得到的。请大家帮我分析一下，为什么会这样呢？是我的模拟时间不够还是怎么回事？多谢了。 [ Last edited by ghcacj on 2010-12-8 at 16:39 ] |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
求助 求两条曲线相关性 如图
已经有13人回复
- 请问一个有理函数的曲线拟合过程如何做，给出数据，请求拟合，
已经有14人回复
- 如何根据已知函数画函数曲线？谢谢
已经有4人回复
- 【求助】动电位极化曲线异常，请教原因
已经有8人回复
- 【求助】相关函数归一化后，输出值怎么变成NaN了？？
已经有4人回复
- 【求助】如何同时拟合两条曲线或者两个函数！【已解决】
已经有7人回复
2楼2010-12-08 16:23:47
astringent
铜虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 2210.7
- 散金: 1233
- 红花: 4
- 帖子: 1180
- 在线: 354.8小时
- 虫号: 900839
- 注册: 2009-11-12
- 专业: 理论和计算化学
3楼2010-12-08 17:34:06
coolrainbow
木虫 (著名写手)
未来国家冻凉
- 应助: 1 (幼儿园)
- 贵宾: 0.1
- 金币: 3979.5
- 散金: 10
- 红花: 50
- 帖子: 1137
- 在线: 272.3小时
- 虫号: 239560
- 注册: 2006-04-07
- 性别: GG
- 专业: 理论和计算化学
4楼2010-12-08 19:18:30
astringent
铜虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 2210.7
- 散金: 1233
- 红花: 4
- 帖子: 1180
- 在线: 354.8小时
- 虫号: 900839
- 注册: 2009-11-12
- 专业: 理论和计算化学
5楼2010-12-09 03:06:27
coolrainbow
木虫 (著名写手)
未来国家冻凉
- 应助: 1 (幼儿园)
- 贵宾: 0.1
- 金币: 3979.5
- 散金: 10
- 红花: 50
- 帖子: 1137
- 在线: 272.3小时
- 虫号: 239560
- 注册: 2006-04-07
- 性别: GG
- 专业: 理论和计算化学
6楼2010-12-09 11:36:09
astringent
铜虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 2210.7
- 散金: 1233
- 红花: 4
- 帖子: 1180
- 在线: 354.8小时
- 虫号: 900839
- 注册: 2009-11-12
- 专业: 理论和计算化学
|
我的采样确实有点问题。在采样阶段，老板只让我做了100ps，每隔100fs取一次样，所以我只有1000个点，所以采样到到最后时，样本就很少。至于算法，我是这样写的： 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)/nparticle*ntime(ical)/(tcf0(ical) /nparticle*ntime(ical) enddo 有劳你帮我看看，这个算法有没有问题，多谢。 |
7楼2010-12-09 15:55:50
8楼2012-10-27 10:12:20
rusong231
铜虫 (著名写手)
- 应助: 6 (幼儿园)
- 金币: 11731.1
- 散金: 15
- 红花: 4
- 帖子: 2316
- 在线: 263.9小时
- 虫号: 728356
- 注册: 2009-03-22
- 性别: GG
- 专业: 金属材料的合金相、相变及
9楼2012-10-27 10:26:00
