求助关于二面角分布概率求法

少年爱吃地三鲜

请教一下各位老师， 现在我想求一下分子链主链上的二面角分布概率， 这个每个时间帧的二面角都可以求出来，现在问题是该如何处理这些结果以得到二面角分布概率？不想通过软件计算， 只想知道它内部原理自己写程序算。 麻烦老师给点建议。具体求什么就行。
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!"

HZW

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

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

少年爱吃地三鲜

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

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

lyj714

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

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

少年爱吃地三鲜

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

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

lyj714

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