布丰的投针试验
公元1777年的一天,法国科学家布丰(D.Buffon1707-1788)的家里宾客满堂,原来他们是应主人的邀请前来观看一次奇特试验的。
试验开始,但见年已古稀的布丰先生兴致勃勃地拿出一张纸来,纸上预先画好了一条条等距离的平行线。接着他又抓出一大把原先准备好的小针,这些小针的长度都是平行线间距离的一半。然后布丰先生宣布:“请诸位把这些小针一根一根往纸上扔吧!不过,请大家务必把扔下的针是否与纸上的平行线相交告诉我。”
客人们不知布丰先生要干什么,只好客随主意,一个个加入了试验的行列。一把小针扔完了,把它捡起来又扔。而布丰先生本人则不停地在一旁数着、记着,如此这般地忙碌了将近一个钟头。最后,布丰先生高声宣布:“先生们,我这里记录了诸位刚才的投针结果,共投针2212次,其中与平行线相交的有704次。总数2212与相交数704的比值为3.142。”说到这里,布丰先生故意停了停,并对大家报以神秘的一笑,接着有意提高声调说:“先生们,这就是圆周率π的近似值!”
众宾哗然,一时议论纷纷,个个感到莫名其妙。“圆周率π?这可是与圆半点也不沾边的呀!”
布丰先生似乎猜透了大家的心思,得意洋洋地解释道:“诸位,这里用的是概率的原理,如果大家有耐心的话,再增加投针的次数,还能得到π的更精确的近似值。不过,要想弄清其间的道理,只好请大家去看敝人的新作了。”说着布丰先生扬了扬自己手上的一本《或然算术试验》的书。
π在这种纷纭杂乱的场合出现,实在是出乎人们的意料,然而它却是千真万确的事实。由于投针试验的问题,是布丰先生最先提出的,所以数学史上就称它为布丰问题。布丰得出的一般结果是:如果纸上两平行线间相距为d,小针
长为l,投针的次数为n,所投的针当中与平行线相交的次数是m,那么当n相当大时有:
在上面故事中,针长l等于平行线距离d的一半,所以代入上面公式简化
我想,喜欢思考的读者,一定想知道布丰先生投针试验的原理,下面就是一个简单而巧妙的证明。
找一根铁丝弯成一个圆圈,使其直径恰好等于平行线间的距离d。可以想象,对于这样的圆圈来说,不管怎么扔下,都将和平行线有两个交点。因此,如果圆圈扔下的次数为n次,那么相交的交点总数必为2n。
现在设想把圆圈拉直,变成一条长为πd的铁丝。显然,这样的铁丝扔下时与平行线相交的情形要比圆圈复杂些,可能有4个交点、3个交点、2个交点、1个交点,甚至于都不相交。
由于圆圈和直线的长度同为πd,根据机会均等的原理,当它们投掷次数较多,且相等时,两者与平行线组交点的总数可望是一样的。这就是说,当长为πd的铁丝扔下n次时,与平行线相交的交点总数应大致为2n。
现在再来讨论铁丝长为l的情形。当投掷次数n增大的时候,这种铁丝跟平行线相交的交点总数m应当与长度l成正比,因而有:m=k l式中k是比例系数。
为了求出k来,只需注意到,对于l=πd的特殊情形,有m=2n。
于
2nl m
dπ
=
这便是著名的布丰公式。
亲爱的读者,你不妨一试。
课本P187-2
随便说出3个正数,以这3个数为边长一定能围成一个三角形吗?一定能围成一个钝角三角形(其中最大边的平方
大于另两边的平方和)吗?估计能围成一个钝角三角形的概率.
本题仍是利用试验的方法估计随机事件发生的概率,随便说出三个正数,以这三个正数为边不一定能组成一个三角形,如不能以1,3,5三个数为边长组成三角形;当然也不一定能组成一个钝角三角形;能围成一个钝角三角形的概率的估计值因人而异,因试验次数而异.事实上,不妨设所取三数为(a ,b ,c(0 应满足:a+b>c ,a2+b2 c a +>1,2 22 2c b c a + <1. 设, , a b x y c c ==则有x+y>1,x2+y2<1,其中0<x ≤1,0<y ≤1 在直角坐标系中,任意a ,b ,c 所对应的图形为正方形OABC 及其内部,而能构成钝角三角形的a ,b ,c 对应的图形是图中的阴影部分。 因此,所求概率为11 4 21 π- =42 -π 抓阄中的概率问题 日常生活中,人们经常通过抓阄对一些事情作出决策。例如,现在有一张去科学宫的参观券,小明、小华、小彬3 个同学都想去,为了公平,可以做3个阄,其中一个阄做上标记,谁抓中做了标记的阄即可得到去科学宫的参观券。可他们3人认为最后抓阄的人没有任何选择的余地,认为抓阄对后抓的人不利,都不愿意最后抓阄。他们的想法正确吗? 这个问题挺复杂。对于复杂问题,不妨先动手试一试,亲身感受一下,或许能得出问题的结论呢! 可以做3个阄(其中一个阄做上标记),3个同学为一组,安排好抓阄顺序,具体地抓抓看。多抓几次,统计一下各人抓中有标记的阄的次数,看看3人抓中有标记的阄的概率如何?当然,时间长了,同学们可能“无意”中会记得各个阄的特征,为了保证抓阄的随机性,可以通过摸球或计算器出示的随机数等进行模拟试验。 通过模拟试验,同学们也许已经得到了问题的结论。但这毕竟是一种感性的认识,能否对此进行理性的分析呢? “可这个问题挺复杂的。” “是的,但我们可以先考虑一个简单的情况呀!”试想只有2人抓阄,同学们不难明白,不管谁先抓,2人抓中有 标记的阄的概率应都是21 。 “那3个人抓阄呢?”不妨依小华、小彬、小明的顺序抓阄,显然原来有3个阄,其中一个做了标记,小华抓中 有标记阄的可能性是31 ,抓不中的可能性是32 ,可用图1表示。只有在小华抓不中的情况下,小彬才有可能抓得有标 记的阄,而且这时他抓得有标记的阄的可能是50%,因此,小彬抓得有标记的阄的可能性是这32 可能性下的21 ,即 31 2132=?。同样,小明得到有标记的阄的可能性也是31 ,如图2。 31 32 31 31 31 图1 图2 可见,3个人抓阄,抓中有标记的阄的可能性与抓阄的顺序并无关系,那n 个人抓阄呢?有兴趣的同学可以通过模拟试验感受一下,也可以仿照上面的思路分析分析。 生活与概率 公元1052年4月,侬智高起兵反宋。当朝皇帝宋仁宗决定派遣大将狄青去平定叛乱。当时路途艰险,军心不稳,狄青取胜的把握不大。为了鼓舞士气,狄青便设坛拜神,说:“这次出兵讨伐叛军,胜败没有把握,是吉是凶,只好由神明决定了。是吉的话,那我随便掷100个铜钱,神明保佑,正面定然会全部朝上;只要有一个背面朝上,那我们就难以制敌,只好回朝了。” 左右官员诚惶诚恐,劝道:“大将军,运气再好,100个铜钱,总不会个个正面朝上,如果有背面朝上,岂不动摇军心?如果不战而回朝,那更是违抗圣旨。请大将军三思而行!”此时的狄青已是胸有成竹,叫心腹拿来一袋铜钱,在千万人的注视下,举手一挥,把铜钱全部抛向空中,100个铜钱居然鬼使神差地全部朝上。顿时,全军欢呼,声音响彻山野。由于士兵个个认定神灵护佑,战斗中奋勇争先,仅一次战役,就收回了失地,大功告成。 那么,那100个铜钱究竟是怎么回事呢?原来,狄青那100个铜钱正反两面都是正面的图案,使得正面朝上的机会为100%,从而鼓舞了士气,大军获胜。 以上只是古人利用简单的概率知识获利。其实,从古到今,概率就与人们的生活息息相关。如今,还有许多不法分子利用人们对概率的不了解牟取暴利。下面,我们就以“机会型”赌博,简要地讲一下如何计算概率以及概率的重要性。 “机会型”赌博规则如下:每个参加者每次先付赌金1元,然后将3枚骰子一起掷出。他可以赌某一个点数,譬如赌“1”点。如果三枚骰子中出现一个“1”点,庄家除把赌金发还外,再奖一元;如果出现两个“1”点,发还赌金外,再奖两元;如果全是“1”,那么发还赌金,再奖三元。 看起来,一枚骰子赌“1”点,取胜的可能性是61 ;那么两枚骰子就是31 的可能性,三枚就是21 。即使是一元对一元的奖励,机会也是均等的,何况还可能是2倍、3倍奖励的可能性,自然对参加者有利。其实,这只是一个假象。 我们来计算一下,三枚骰子一起掷,会出现怎样的情况。见表1。 表1 1”点的情况:出现“1”点的骰子可能是第一枚,也可能是第二枚或第三枚,共有三种可能;而其余两枚不出现“1”点的可能性有5×5=25种,所以共有3×25=75种可能。这75种可能出现时,它可获2元,那么总共可获75×2=150元。再来看出现两枚“ 1”点的可能性:可以出现在第一枚和第二枚,也可以是第一枚和第三枚,还可以是第二枚和第三枚,也是三种可能;而另一枚骰子不出现“1”点只有5种可能,所以共有15种可能。这时,每次他可获3元,共45元。最后,三枚都出现“1”点的只有一种可能,这时,它可获4元。 这样,216次,他共获150+45+4=199元。但每次先付一元,他一共付了216元。所以,一般来说,他会输216-199=17元,占总金额的7.9%。 我们再来看看庄家的情况。根据前面的分析过程,假使有6人参加赌博,每人分别赌“1”、“2”、……“6”点,并假定每人进行了216次,则庄家共收了6×216=1296元,一共付出了720+450+24=1194元,净赚1296-1194=102元,占总金额的7.9%。 通过概率的计算,我们看到赢的一定是庄家。看清了赌博的真面目,我们就应该抵制赌博。 同样我们可以利用概率计算动物的寿命,以乌龟的寿命为例,如表2: 根据表2内容,再计算出,活满20岁的乌龟有0.87÷0.92×100%=95%的概率可活到80岁,活满120岁的乌龟有0.39÷0.87×100%=50%的概率可活到200岁。 同理,通过大量调查数据获得人类的寿命表,保险公司便可算出保险费率。 以上两个例子说明,概率与人们的生活息息相关,只要你熟练地掌握了概率的知识,并应用到日常生活中去,我想你就能做到较好地把握机会,将胜算牢牢地掌握在自己的手中。 ——选自首届初中生数学学习《“时代之星”实践与创新论文大赛》江苏教育出版社(有改动) 象棋比赛阵容 少年宫请来了一位象棋大师,他对少年象棋队的队员们做了一些辅导之后,决定与少年棋手来几盘棋赛。大师的棋艺高出少年棋手好多好多,怎么能比呢?不要紧,大师下的是盲棋——不看棋盘,由别人将对手的走着告诉大师,大师再把自己的走着告诉这个人,由他代走。 比赛作了这样的约定:由少年象棋队挑出两名队员,轮流与大师赛棋,共赛三盘。如果能连胜大师两盘,就算少年棋队胜。注意:是连胜两盘,不是共胜两盘。 假定少年棋手甲能胜大师的概率是0.75,乙能胜大师的概率是0.5,那么少年棋队应该用“甲—乙—甲”,还是用“乙—甲—乙”的阵容来对付大师呢? “当然用‘甲—乙—甲’阵容啦!甲是我队最好的队员嘛!”少年棋队的队员们一致这样看。 其实,“甲—乙—甲”阵容战胜大师(连胜两盘)的概率比“乙—甲—乙”阵容战胜大师的概率要小一些。 为什么呢?我们在这里只做一些直观的解释。 用“甲—乙—甲”阵容参战,最佳的棋手可以上场两次,看来好像是有利的。但是,我们现在的规则是:连胜两盘才能算少年队赢。用这个阵容,即使甲胜了两盘,也没用,因为不是“连胜”两盘。 要连胜两盘,必须在第二盘比赛中取胜,因此第二盘比赛是关键。而“乙—甲—乙”阵容,就是把最佳选手安排在最关键的场合,所以是较好的方案。 系统建模与仿真题目:Buffon实验的仿真 院系: 电子工程学院 专业:信息对抗技术 班级:021231 姓名:余颖智 学号:02123021 指导老师:刘洋 完成时间:2015年4月 西安电子科技大学 基于MATLAB的投针实验仿真 摘要 在求证圆周率的过程中经过割圆术后,出现的投针试验以求出圆周率,目前利用MATLAB数学建模的仿真实验,运用到计算机中,简化其随机实验的操作量大,运算慢等特点。不同针距相同实验量运算后得出不同的π,其针距与线间距离相等,所得值接近于π。 目录 摘要 (2) 二、实验内容 (4) 三、建模流程图 (5) 四、程序主要代码 (6) 五、运行结果 (6) 六、结论 (7) 一、实验原理 1777年法国科学家布丰提出的一种计算圆周率的方法——随机投针法,即著名的布丰投针问题。该投针实验主要有如下三个步骤:(一)取一张白纸,在上面画许多条间距为a的平行线;(二)取一根长度为l(l 三、建模流程图 四、程序主要代码 str(handles.edit1,'string'); %取得变量,定义变量,变量初始化 n = str2double(str); str = get(handles.edit2,'string'); l = str2double(str); str = get(handles.edit3,'string'); a = str2double(str); counter = 0; %变量初始化 phi = 0; frequency = 0; Pi = 0; x = unifrnd(0,a/2,1,n);%产生n个(0,a/2)之间均匀分布的随机数,这里a/2是投针的中点到最近的平行线的距离 phi = unifrnd(0,pi,1,n);% 产生n个(0,pi)之间均匀分布的随机数,这里pi是投针与最近平行线的角度 for i=1:n if x(i) 系统建模与仿真 基于MATLAB的布丰实验模拟 姓名:石星宇 学号: 02123010 指导教师:刘洋 2015年4月9日 目录 基于MATLAB的布丰实验模拟 .................................................................... - 1 - 一、实验原理......................................................................................... - 1 - 二、编程模拟......................................................................................... - 1 - 1、程序流程图............................................................................... - 1 - 2、程序代码................................................................................... - 2 - 三、实验结果......................................................................................... - 2 - 基于MATLAB 的布丰实验模拟 一、实验原理 找一根铁丝弯成一个圆圈,使其直径恰恰等于平行线间的距离a 。可以想象得到,对于这样的圆圈来说,不管怎么扔下,都将和平行线有两个交点。因此,如果圆圈扔下的次数为n 次,那么相交的交点总数必为n 2。现在设想把圆圈拉直,变成一条长为a π的铁丝。显然,这样的铁丝扔下时与平行线相交的情形要比圆圈复杂些,可能有4个交点,3个交点,2个交点,1个交点,甚至于都不相交。由于圆圈和直线的长度同为a π,根据机会均等的原理(即等概率事件),当它们投掷次数较多,且相等时,两者与平行线组交点的总数期望也是一样的。这就是说,当长为a π的铁丝扔下n 次时,与平行线相交的交点总数应大致为n 2。现在转而讨论铁丝长为l 的情形。当投掷次数n 增大的时候,这种铁丝跟平行线相交的交点总数k 应当与长度l 成正比,因而有:l k λ=,式中λ是比例系数。为了求出λ来,只需注意到,对于a l π=的特殊情形,有n k 2=。于是求得a n πλ2=。代入前式就有:a m πln 2≈从而ak nl 2≈π。 二、编程模拟 1、程序流程图 参数初始化 产生位置随机数; 产生角度随机数 判断相交 1+=k k 1+=n n 是 否 判断结束 公元1777年的一天,法国科学家D·布丰(D·buffon,1707~1788)的家里宾客满堂,原来他们是应主人的邀请前来观看一次奇特试验的. 试验开始,但见年已古稀的布丰先生兴致勃勃地拿出一张纸来,纸上预先画好了一条条等距离的平行线.接着他又抓出一大把原先准备好的小针,这些小针的长度都是平行线间距离的一半.然后布丰先生宣布:“请诸位把这些小针一根一根往纸上扔吧!不过,请大家务必把扔下的针是否与纸上的平行线相交告诉我.” 客人们不知布丰先生要干什么,只好客随主意,一个个加入了试验的行列.一把小针扔完了,把它捡起来又扔.而布丰先生本人则不停地在一旁数着、记着,如此这般地忙碌了将近一个钟头.最后,布丰先生高声宣布:“先生们,我这里记录了诸位刚才的投针结果,共投针2212次,其中与平行线相交的有704次.总数2212与相交数704的比值为3.142.”说到这里,布丰先生故意停了停,并对大家报以神秘的一笑,接着有意提高声调说:“先生们,这就是圆周率π的近似值!” 众宾哗然,一时议论纷纷,个个感到莫名其妙;“圆周率π?这可是与圆半点也不沾边的呀!” 布丰先生似乎猜透了大家的心思,得意洋洋地解释道:“诸位,这里用的是概率的原理,如果大家有耐心的话,再增加投针的次数,还能得到π的更精确的近似值.不过,要想弄清其间的道理,只好请大家去看敝人的新作了.”随着布丰先生扬了扬自己手上的一本《或然算术试验》的书. π在这种纷纭杂乱的场合出现,实在是出乎人们的意料,然而它却是千真万确的事实.由于投针试验的问题,是布丰先生最先提出的,所以数学史上就称它为布丰问题.布丰得出的一般结果是:如果纸上两平行线间相距为d,小针长为l,投针的次数为n,所投的 针当中与平行线相交的次数是m,那么当n相当大时有:π≈2ln dm .在上面故事中,针长l 等于平行线距离d的一半,可以代入上面公式简化.我想,喜欢思考的读者一定想知道布丰先生投针试验的原理,下面就是一个简单而巧妙的证明. 找一根铁丝弯成一个圆圈,使其直径恰恰等于平行线间的距离d.可以想象得到,对于这样的圆圈来说,不管怎么扔下,都将和平行线有两个交点.因此,如果圆圈扔下的次数为n次,那么相交的交点总数必为2n. 现在设想把圆圈拉直,变成一条长为πd的铁丝.显然,这样的铁丝扔下时与平行线相交的情形要比圆圈复杂些,可能有4个交点,3个交点,2个交点,1个交点,甚至于都不相交. 由于圆圈和直线的长度同为πd,根据机会均等的原理,当它们投掷次数较多且相等时,两者与平行线组交点的总数可望是一样的.这就是说,当长为πd的铁丝扔下n次时,与平 概率论与数理统计实验 蒲丰投针与蒙特卡罗法 班级应数12级01班 学号2012444086 姓名张旭东 蒲丰投针与蒙特卡罗法 张旭东2012444086 (重庆科技学院数学与应用数学,重庆沙坪坝) 【摘要】通过设计一个投针实验使这个事件的概率和未知量π有关,然后通过重复实验,以频率估计概率,即可求得未知参数π的近似解。这种方法称为随机模拟法,也称为蒙特卡罗法。一般来说,实验次数越多所得的近似值就越接近真值。可以利用MATLAB来大量重复地模拟所设计的随机实验。 【关键词】随机模拟;投针实验;重复实验 1 引言 蒲丰投针问题是由法国科学家蒲丰(Buffon)在1777年提出的,它是概率中非常有代表性的问题,它是第一个用几何形式表达概率问题的例子,其结论具有很强的理论与实际意义。蒲丰针问题的解决不仅较典型的反应了集合概率的特征及处理方法,而且还可以由此领略到从“概率土壤”上开出的一朵瑰丽的鲜花——蒙特卡洛(Monte-Carlo)方法。 蒙特卡罗(Monte Carlo)方法,也称计算机模拟方法,是一种基于“随机数”的计算方法,大数定律为近年来发展迅速的随机计算机和随机模拟方法提供了理论基础。 MATLAB是一个适合多学科,具有多种工作平台的功能强大的大型软件。MATLAB已经成为线性代数、自动控制理论、数理统计、数字信号处理、时间序列分析、动态系统仿真等高级课程的进本教学工具,Matlab随机数发生器的种类丰富且用法简便。 本文介绍了利用随机模拟方法和大数定律的相关理论解决蒲丰投针问题计算π的近似值。 2 有关数学实验的有关基础 定理(贝努力大数定律) 设n μ是n 重贝努力实验中事件A 出现的次数,P 是事件A 每次实验中出现的概率,即P(A)=p,则对任意的 ε>0,有 3 实验 蒲丰投针问题 在平面上画有等距离的一些平行线,平行线间的距离为a(a>0),向平面上随机投一长为l(l基于MATLAB的布丰投针实验仿真
布丰投针实验模拟
苏科版-数学-九年级上册-知识拓展 布丰的投针试验
蒲丰投针实验模拟
蒲丰氏投针问题的模拟过程