计算化学公社

标题: 求助关于二面角分布概率求法 [打印本页]

作者
Author: 少年爱吃地三鲜 时间: 2021-1-24 21:16
标题: 求助关于二面角分布概率求法
请教一下各位老师，现在我想求一下分子链主链上的二面角分布概率，这个每个时间帧的二面角都可以求出来，现在问题是该如何处理这些结果以得到二面角分布概率？不想通过软件计算，只想知道它内部原理自己写程序算。麻烦老师给点建议。具体求什么就行。
PS:
源程序在这里 http://bbs.keinsci.com/thread-14821-1-1.html    #### ene老师的TCL程序中这一部分没有看懂万望老师指点一二。

set index [expr int(floor([expr ($DIHED + 180) / 2]))] #  为啥/2
set density($index) [expr $density($index) + 1.0]
}
for { set i 1 } { $i < 180 } { incr i } {    #  这个for循环啥意思？
set density($i) [expr $density($i) / $nf * 100]
puts -nonewline $outfile "[expr $i * 2 - 1 - 180 ]"
puts -nonewline $outfile " "
puts $outfile "[expr $density($i)]"
}
close $outfile
puts "All Done!"

作者
Author: lyj714 时间: 2021-1-25 09:59
本帖最后由 lyj714 于 2021-1-25 10:23 编辑

这就是个简单的直方图统计，将360平分为180份，每一段2°，然后统计每个bin出现的频数。图中的分布是将得到的频数除以(360*帧数)得到这种曲线与角度所围面积为1的角度分布图

作者
Author: 少年爱吃地三鲜 时间: 2021-1-25 11:21

lyj714 发表于 2021-1-25 09:59
这就是个简单的直方图统计，将360平分为180份，每一段2°，然后统计每个bin出现的频数。图中的分布是将得到 ...

老师，请问，我统计频数时一共有370个二面角 1001帧，此时我是把这个370个二面角在每一个帧数下都做统计对吗

作者
Author: lyj714 时间: 2021-1-25 11:27

少年爱吃地三鲜发表于 2021-1-25 11:21
老师，请问，我统计频数时一共有370个二面角 1001帧，此时我是把这个370个二面角在每一个帧数下都 ...

看你实际算的是什么东西，如果是370个不同的二面角，那么就分别算，这样当然会得到370个二面角分布图。如果是相同类的二面角，取个平均就行了，这样得到的平均后的一个二面角分布图。

作者
Author: 少年爱吃地三鲜 时间: 2021-1-26 14:29

lyj714 发表于 2021-1-25 11:27
看你实际算的是什么东西，如果是370个不同的二面角，那么就分别算，这样当然会得到370个二面角分布图。如 ...

您说这个取平均，指的是假如分成了 360度每个度数上得统计数量除以370吗？这个平均值是这么求吗

作者
Author: HZW 时间: 2021-4-21 22:42

少年爱吃地三鲜发表于 2021-1-26 14:29
您说这个取平均，指的是假如分成了 360度每个度数上得统计数量除以370吗？这个平均值是这么求吗

你好，请问您这个图片的那篇文献名字能发我一下吗？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)