欢迎光临散文网 会员登陆 & 注册

布丰投针试验——几何概型

2021-12-05 10:04 作者:匆匆-cc  | 我要投稿

        18世纪,布丰提出以下问题:

        设一块含有平行且等距的木纹的地板,随意抛一枚针,求针和木纹相交的概率。

        这就是著名的布丰投针试验

        我们将这个问题抽象为:

        一组等距平行线上任意画定长线段,求该线段与平行线相交的概率。

布丰投针试验——模拟图

        为解题需要,我们需要设以下参数: 

        x:线段中心距离与之最近的直线的距离;

        θ:线段与平行线的夹角;

         l:线段长;

        d:平行线的间距。

        这里,由于投针试验为随机试验,显然有:

x%20%5Cin%20%5B0%2C%5Cfrac%7Bd%7D%7B2%7D%20%5D,且 x 为线均匀分布

%5Ctheta%20%5Cin%20%5B0%2C%5Cfrac%7B%5Cpi%20%7D%7B2%7D%20%5D,且 θ 为角均匀分布

        为方便起见,我们仅拿其中一根针来研究。

        由图可知,相交条件为:

x%20%5Cleq%20%5Cfrac%7Bl%7D%7B2%7D%20%5Csin%20%5Ctheta%20%20

布丰投针试验——参数关系示意图

        根据一定的非线性规划知识,这表示为下图中ABE区域的面积,而样本空间区域面积为下图中ABCD的面积。

布丰投针试验——面积区域示意图

        为了求解此题概率,我们引入几何概型

几何概型:

        若试验只可能出现无限多种基本事件,且各种基本事件发生概率仅与构成该事件区域的长度(面积、体积)有关,则称此概率模型为几何概型。

P_%7B(A)%7D%20%20%3D%20%5Cfrac%7B%E6%9E%84%E6%88%90%E4%BA%8B%E4%BB%B6A%E7%9A%84%E5%8C%BA%E5%9F%9F%E9%95%BF%E5%BA%A6%EF%BC%88%E9%9D%A2%E7%A7%AF%E3%80%81%E4%BD%93%E7%A7%AF%EF%BC%89%7D%7B%E8%AF%95%E9%AA%8C%E6%A0%B7%E6%9C%AC%E7%A9%BA%E9%97%B4%E7%9A%84%E5%8C%BA%E5%9F%9F%E9%95%BF%E5%BA%A6%EF%BC%88%E9%9D%A2%E7%A7%AF%E3%80%81%E4%BD%93%E7%A7%AF%EF%BC%89%7D%20

        举一个简单的例子,用蒙特卡罗方法估算 π 的值,python 程序实现如下:

        程序运行结果如下。

蒙特卡罗方法求 π——程序试验

        可以看到,在试验频数增加的情况下,数值逐渐向 π 逼近。

        蒙特卡罗方法求 π 的本质即为几何概型。

        如图,根据几何概型概率的定义,有:

P%20%20%3D%20%5Cfrac%7BS%E6%89%87%E5%BD%A2%7D%7BS%E6%AD%A3%E6%96%B9%E5%BD%A2%7D%20%3D%20%5Cfrac%7B%5Cfrac%7B1%7D%7B4%7D%5Cpi%20a%5E2%20%7D%7Ba%5E2%20%7D%20%3D%5Cfrac%7B%5Cpi%20%7D%7B4%7D%20


蒙特卡罗方法求 π——几何解释

        是不是很巧妙?

        我们回到布丰投针试验。

        由非线性规划图,我们无法比较%5Cfrac%7Bd%7D%7B2%7D%20%5Cfrac%7Bl%7D%7B2%7D%20的大小关系,因此需要分类讨论。

    ①d%20%5Cgeq%20l

        由图可知,

P%20%3D%20%5Cfrac%7B%5Cint_%7B0%7D%5E%7B%5Cfrac%7B%5Cpi%7D%7B2%7D%7Dxd%5Ctheta%7D%7B%5Cfrac%7B%5Cpi%7D%7B2%7D%5Ccdot%5Cfrac%7Bd%7D%7B2%7D%7D%20%3D%5Cfrac%7B2l%5Cint_%7B0%7D%5E%7B%5Cfrac%7B%5Cpi%7D%7B2%7D%7D%5Csin%20%5Ctheta%20d%5Ctheta%20%7D%7B%5Cpi%20d%7D%20%3D%5Cfrac%7B2l%7D%7B%5Cpi%20d%7D%20

注:%5Cint_%7B%7D%5E%7B%7D%20%5Csin%20x%20dx%20%3D%20-%5Ccos%20x%20%2BC

布丰投针试验——情况1

        若取d%20%3D%202l,即可得到:

P%20%3D%20%5Cfrac%7B1%7D%7B%5Cpi%20%7D%20

        令人震惊!一个看起来仿佛与 π 没有什么关系的问题,最终结果出现了 π !

    ②d%20%3C%20l

        由图可知,直线与曲线有一交点,设交点横坐标为%5Ctheta%20_%7B0%7D%20,两面积区域分别为区域ABCF与区域ABCD。

        我们有:

%5Cfrac%7Bl%7D%7B2%7D%5Csin%20%5Ctheta%20_%7B0%7D%20%20%3D%20%5Cfrac%7Bd%7D%7B2%7D%20

%5Csin%20%5Ctheta%20_%7B0%7D%20%20%3D%20%5Cfrac%7Bd%7D%7Bl%7D%20

P%20%3D%20%5Cfrac%7B%5Cint_%7B0%7D%5E%7B%5Ctheta%20_%7B0%7D%20%7D%20xd%5Ctheta%20%2B(%5Cfrac%7B%5Cpi%20%7D%7B2%7D%20-%5Ctheta%20_%7B0%7D%20)%5Ccdot%20%5Cfrac%7Bd%7D%7B2%7D%20%7D%7B%5Cfrac%7B%5Cpi%20%7D%7B2%7D%5Ccdot%20%5Cfrac%7Bd%7D%7B2%7D%20%20%7D%20%3D%5Cfrac%7B2l%20%5Cint_%7B0%7D%5E%7B%5Ctheta%20_%7B0%7D%7D%5Csin%20%5Ctheta%20d%5Ctheta%20%2B%5Cpi%20d%20-%202d%20%5Ccdot%20%5Ctheta_%7B0%7D%20%7D%7B%5Cpi%20d%7D%20%3D%201-%5Cfrac%7B2%7D%7B%5Cpi%7D%5Carcsin%5Cfrac%7Bd%7D%7Bl%7D%20%2B%5Cfrac%7B2l%7D%7B%5Cpi%20d%7D(1-%5Csqrt(1-%5Cfrac%7Bd%5E2%20%7D%7Bl%5E2%7D))

布丰投针试验——情况2

        注:%5Clim%5Cnolimits_%7B%5Cfrac%7Bd%7D%7Bl%7D%20%5Cto0%7D%20P%20%3D%201%20,即线段非常长时线段与平行线必然有交点。

        以下是关于布丰投针试验的python程序模拟。

        测试结果较好。

布丰投针试验——程序试验

后记:

        创作缘由就是数学高考紫五三下的拓展介绍。

        初次尝试,耗时较长;有所纰漏,还请原谅。

布丰投针试验——几何概型的评论 (共 条)

分享到微博请遵守国家法律