R语言对布丰投针(蒲丰投针)实验进行模拟和动态可视化生成GIF动画

布丰投针是几何概率领域中最古老的问题之一。它最早是在1777年提出的。

它涉及将针头放到衬有衬纸的纸上,并确定针头越过页面上一行的可能性。引人注目的结果是概率与pi的值直接相关。R程序将根据上段所述的情况估算pi的值并使用gganimate进行动态可视化。

第1部分

对于A部分,我们创建一个数据帧,该数据帧将在3个不同的间隔上生成随机值,这些间隔将代表x,y的范围以及每个落针点的角度。这是一个易于实现的随机数情况,需要使用runif函数。此功能要求输入数量,后跟一个间隔。生成数字后,我们会将值保存到数据框中。


课程

R语言数据分析挖掘必知必会

面对扑面而来的数据浪潮,包含Google、Facebook等国际企业,都已采用R语言进行数据分析

探索课程 ➜

第2部分

我们绘制第一部分中的针。重要的是不要在这个问题上出现超过2条水平线。它使我们可以进行较小的检查以了解此处描绘的几何特性的一般概念。

话虽如此,让我们注意我们决定在每个方向上将图形扩展1个单位。原因是想象一个针的尾巴从y = 1开始,其角度为pi / 2。我们需要假设该方向的范围最大为2。

第3部分

在下面,将基于阅读布冯针和基本几何原理的知识,查看pi的估算值。

第4部分

运行代码后,我们收到以下答案。> buffon(X)[1] 3.846154

第5部分

如前几节所述,当我们放下更多的针头时,我们期望以最小的可变性获得更准确的答案。从Approxpi函数运行代码后,我们收到了平均值= 3.172314和方差0.04751391的值。对于这样一个简单的实验,它对pi进行了很高的估计。

接下来对模拟次数从500~600的预测进行动态可视化,红色表示针投放到了直线上:

参考资料

Schroeder,L.(1974年)。布冯针问题:许多数学概念的激动人心的应用。


可下载资源

​非常感谢您阅读本文,如需帮助请联系我们!


关于作者

Kaizong Ye拓端研究室(TRL)的研究员。

本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。


随时关注您喜欢的主题

在wechat上关注我们

最新洞察

技术干货

Leave A Reply

电子邮件地址不会被公开。

 
QQ在线咨询
售前咨询热线
15121130882
售后咨询热线
0571-63341498
error: Content is protected
%d 博主赞过: