f*******i
发帖数: 8492 | 1 【 以下文字转载自 Statistics 讨论区 】
发信人: fantastli (早晨从中午开始), 信区: Statistics
标 题: 双簧包求助，请教该如何处理这批数据
发信站: BBS 未名空间站 (Sat Dec 4 13:44:51 2010, 美东)
我的最终目的是想得到radial distribution function的曲线，类似于附件中的图形。
我设计的程序如下，在一维空间中，20个小球两两互相撞击，每隔一定时间，就记录一
下所有小球距离中间位置的小球的距离。这样，每次可以得到19个数据。
如果程序运行1000次，就有19000个数据点。
另外，如果输入30个小球的话，且程序运行10000次的话，就会有290000个数据点。
理论上，应该模拟的次数越多，就越接近理想的图形（图示中的）。
在图示中，纵坐标是g(r)，定义为g(r)=（dN/dr)/(N/r)；横坐标是r
而我的程序运行完后，会进行类似于excel中的Histogram命令，既以0.1为bin的步长，
把0.1-0.2,0.2-0.3...把所有在此范围内的点计数，最终得到的是以frequency为纵坐标，以
bin的长度为横坐标的图形。
现在，我的问题是这个图形和目标图形的样子差距还是很大的。请问该用哪种数据处理
方法来处理呢？另外，如何消除运行次数，与小球个数的不同所造成的图形差异？谢谢 |
|
