计算机投针实验程序,蒲丰投针问题
18世纪,蒲丰提出以下问题:设我们有一个以平行且等距木纹铺成的地板(如概述图),随意抛一支长度比木纹之间距离小的针,求针和其中一条木纹相交的概率。并以此概率,布丰提出的一种计算圆周率的方法——随机投针法。这就是蒲丰投针问题(又译“布丰投针问题”)。
中文名
蒲丰投针问题
外文名
Buffon's needle problem
别 名
布丰投针问题提出者
布丰
提出时间
1777年
适用领域
应用学科
计算圆周率,概率学
蒲丰投针问题投针步骤
编辑
语音
法国数学家布丰(1707-1788)最早设计了投针试验。
这一方法的步骤是:
1) 取一张白纸,在上面画上许多条间距为a的平行线。
2) 取一根长度为l(l≤a) 的针,随机地向画有平行直线的纸上掷n次,观察针与直线相交的次数,记为m。
3)计算针与直线相交的概率.
18世纪,法国数学家布丰提出的“投针问题”,记载于布丰1777年出版的著作中:“在平面上画有一组间距为a的平行线,将一根长度为l(l≤a)的针任意掷在这个平面上,求此针与平行线中任一条相交的概率。”
布丰本人证明了,这个概率是:
(其中π为圆周率)
由于它与π有关,于是人们想到利用投针试验来估计圆周率的值。
布丰惊奇地发现:有利的扔出与不利的扔出两者次数的比,是一个包含π的表示式.如果针的长度等于a/2,那么扔出的概率为1/π.扔的次数越多,由此能求出越为精确的π的值。
蒲丰投针问题实验数据
下面是利用这个公式,用概率的方法得到圆周率的近似值的一些资料。试验者时间投掷次数相交次数圆周率估计值
Wolf1850年500025323.1596
1855年32041218.53.1554
C.De Morgan1860年600382.53.137
1884年10304893.1595
Lazzerini1901年340818083.1415929
1925年25208593.1795
公元1901年,意大利数学家拉兹瑞尼宣称进行了多次的投针试验,每次投针数为3408次,平均相交数为1808次,给出π的值为3.1415929——准确到小数后6位。不过,不管拉兹瑞尼是否实际上投过针,他的实验还是受到了美国犹他州奥格登的国立韦伯大学的L·巴杰的质疑.通过几何、微积分、概率等广泛的范围和渠道发现π,这是着实令人惊讶的!
布丰投针实验是第一个用几何形式表达概率问题的例子,他首次使用随机实验处理确定性数学问题,为概率论的发展起到一定的推动作用。[1]
蒲丰投针问题证明
证明一:找一根铁丝弯成一个圆圈,使其直径恰恰等于平行线间的距离d。可以想象得到,对于这样的圆圈来说,不管怎么扔下,都将和平行线有两个交点。因此,如果圆圈扔下的次数为n次,那么相交的交点总数必为2n。设想把圆圈拉直,变成一条长为πd的铁丝。显然,这样的铁丝扔下时与平行线相交的情形要比圆圈复杂些,可能2个交点,1个交点,甚至于都不相交。由于圆圈和直线的长度同为πd,根据机会均等的原理,当它们投掷次数较多,且相等时,两者与平行线组交点的总数期望也是一样的。这就是说,当长为πd的铁丝扔下n次时,与平行线相交的交点总数应大致为2n。
转而讨论铁丝长为l的情形。当投掷次数n增大的时候,这种铁丝跟平行线相交的最大的交点总数m应当与长度l成正比,因而有:m=kl,式中k是比例系数。
为了求出k来,注意到l=πd时的特殊情形,有m=2n。于是求得。
代入前式就有:,将此结论推广到l=a/2,那么最多也只有一个交点,m与n的比值是针与直线相交的概率。但此证明较不严谨,例如圆和直线期望相等,铁丝与平行线的交点成正比。接下来用概率论和微积分提供严谨的证明。
证明二:由于向桌面投针是随机的,所以用二维随机变量(X,Y)来确定它在桌上的具体位置。设X表示针的中点到平行线的距离,Y表示针与平行线的夹角,如果时,针与直线相交。
并且X在服从均匀分布,Y在服从均匀分布,XY相互独立,由此可以写出(X,Y)的概率密度函数
因此所求概率
蒲丰投针问题蒙特卡罗方法
编辑
语音
像投针实验一样,用通过概率实验所求的概率来估计我们感兴趣的一个量,这样的方法称为蒙特卡罗方法(Monte Carlo method)。当由于这类模型含有不确定的随机因素,分析起来通常比确定性的模型困难.有的模型难以作定量分析,得不到解析的结果,或者是虽有解析结果,但计算代价太大以至不能使用.在这种情况下,可以考虑采用 Monte Carlo 方法,蒙特卡罗方法是在第二次世界大战期间随着计算机的诞生而兴起和发展起来的。这种方法在应用物理、原子能、固体物理、化学、生态学、社会学以及经济行为等领域中得到广泛利用。
利用钝角三角形的边长计算圆周率
此外,随便说出3个正数,以这3个正数为边长可以围成一个钝角三角形的概率P也与π有关,这个概率为 (π-2)/4,证明如下:
设这三个正数为x,y,z,不妨设x≤y≤z,对于每一个确定的z,则必须满足x+y>z,x^2+y^2﹤z^2,容易证明这两个式子即为以这3个正数为边长可以围成一个钝角三角形的充要条件,用线性规划可知满足题设的可行域为直线x+y=z与圆x^2+y^2=z^2;围成的弓形,总的可行域为一个边长为z的正方形,则可以围成一个钝角三角形的概率P=S弓形/S正方形=(πz^2/4-z^2/2)/z^2=(π-2)/4.因为对于每一个z,这个概率都为(π-2)/4,因此对于任意的正数x,y,z,有P=(π-2)/4,命题得证。
为了估算π的值,我们需要通过实验来估计它的概率,这一过程可交由计算机编程来实现,事实上x+y>z,x^2+y^2﹤z^2等价于(x+y-z)(x^2+y^2-z^2)﹤0,因此只需检验这一个式子是否成立即可。若进行了m次随机试验,有n次满足该式,当m足够大时,n/m趋近于(π-2)/4,令n/m=(π-2)/4,解得π=4n/m+2,即可估计出π值。
值得注意的是这里采用的方法:设计一个适当的试验,它的概率与我们感兴趣的一个量(如π)有关,然后利用试验结果来估计这个量,随着计算机等现代技术的发展,这一方法已经发展为具有广泛应用性的蒙特卡罗方法。
计算机模拟
Monte Carlo方法是计算机模拟的基础,它的名字来源于世界著名的赌城——摩纳哥的蒙特卡洛, 其历史起源于 1777 年法国科学家蒲丰提出的一种计算圆周π 的方法——随机投针法,即著名的蒲丰投针问题。
Monte Carlo方法的基本思想是首先建立一个概率模型,使所求问题的解正好是该模型的参数或其他有关的特征量. 然后通过模拟一统计试验, 即多次随机抽样试验 (确定 m和 n) ,统计出某事件发生的百分比。只要试验次数很大,该百分比便近似于事件发生的概率.这实际上就是概率的统计定义。利用建立的概率模型,求出要估计的参数。蒙特卡洛方法属于试验数学的一个分支。
MATLAB语言编程实现:
l=1;
n=1000;
d=2;
m=0;
for k=l:n
x=unifrnd(0,d/2);
p=unifrnd(0,pi);
if x<0.5*sin(p)
m=m+1;
else
end
end
p=m/n
pi_m=1/p
运行,即得结果。
c++语言编程实现:
#include
#include
#include
#include
using namespace std;
int main()
{
longi,in,N=1000000;
doublex,y,pi;
srand(time(NULL));
for(i=0,in=0;i
{
x=2.0*rand()/RAND_MAX-1;
y=2.0*rand()/RAND_MAX-1;
if((x*x+y*y)<=1)
in++;
}
pi=4.0*in/N;
cout<
getch();
return 1;
}
蒙特卡洛方法适用范围很广泛,它既能求解确定性的问题,也能求解随机性的问题以及科学研究中的理论问题.例如利用蒙特卡洛方法可以近似地计算定积分,即产生数值积分问 题。
任意曲边梯形面积的近似计水塘的面积.应该怎样做呢?测量方法如下:假定水塘位于一块面积已知的矩形农田之中。如图1所示:随机地向这块农田扔石头使得它们都落在农田内。被扔到农田中的石头可能溅上了水,也可能没有溅上水,估计被“溅上水的”石头量占总的石头量的百分比。试想如何利用这估计的百分比去近似计算该水塘面积?
图1 蒙特卡洛方法(8.2)
词条图册
更多图册
参考资料
1.
马复.义务教育课程标准实验教科书数学九年级上册:北京师范大学出版社,2012.6:185-187
计算机投针实验程序,蒲丰投针问题相关推荐
- 计算机投针实验程序,蒲丰投针实验SAS程序
投针实验 data a; do i=1 to 5000; /*实验5000次*/ pi=constant('pi');/*得出pi真实值*/ x=0.5*ranuni(0);/*产生随机数,a=1*/ ...
- 蒲丰投针计算机模拟ppt,蒲丰投针实验模课件.doc
概率论与数理统计实验 蒲丰投针与蒙特卡罗法 班级 应数12级01班 学号 2012444086 姓名 张旭东 蒲丰投针与蒙特卡罗法 张旭东 2012444086 (重庆科技学院 数学与应用数学 ,重庆 ...
- java实现蒲丰投针求_蒲丰投针问题
问题描述 平面上画有间隔为d 的等距平行线,向平面任意投掷一枚长为 l 的针,求针与平行线相交的概率. 解: 以 \(x\) 表示针的中点与最近一条平行线的距离.又以 \(\varphi\) 表示针与 ...
- Matlab模拟蒲丰投针实验计算Π值
文章目录 1.前言 2.实验目的 3.实验原理 (1)问题描述 (2)问题求解 4.实验过程 5.实验结果 6.实验结论 7.Matlab代码 1.前言 学习了概率论,决定自己动手编程实现蒲丰投针实验 ...
- c语言模拟随机投点计算概率,我写的蒲丰投针的程序(用概率学计算π),有问题求教...
谢谢各位的回答. 此题我已经解决了. 问题是我提议没理解清楚,估计还误导大家了. 蒲丰投针的平行线是一组的,我以为是2条.呵呵. 那么这样的话相交条件改一下,就行了. 下面是我的程序,有兴趣的朋友可以 ...
- 强化学习 蒲丰投针实验 蒙特卡洛算法
目录 一.蒲丰投针实验 1.1背景故事 1.2原理介绍 二.蒙特卡洛方法 2.1方法介绍 2.2经验轨迹 2.3在线学习与离线学习 2.4数学原理 一.蒲丰投针实验 1.1背景故事 1777年,蒲丰请 ...
- 【游戏仿真实验】Unity仿真蒲丰投针实验,丢针计算圆周率,丢了一百万次针得出的结果是...
文章目录 一.前言 二.制作场景 1.制作针模型 2.桌面制作 3.平行线制作 三.物理仿真 1.桌面无反弹 2.针掉落 3.针滚动问题 4.针架到平行线上的问题 5.针与针相互影响的问题 四.UI界 ...
- 一文读懂蒲丰投针实验
这个博客源于概率论复习期间,蒙特卡罗方法的思想起源,这种求π\piπ的思路非常的巧妙 附:历史上用蒲丰投针实验估计圆周率的实验记录,来源 蒲丰投针 蒲丰投针实验是法国数学家.自然科学家"乔治 ...
- 蒲丰投针实验的一些理解
由于要做ppt这里把自己的代码先暂时放上来,不当之处请多指正! 实验过程: 1)取一张白纸,在上面画上许多条间距为a的平行线. 2) 取一根长度为l(l≤a) 的针,随机地向画有平行直线的纸上掷n次, ...
- 蒲丰投针问题的matlab仿真
吃饭的时候突然想到蒲丰投针问题. 在许多地方都遇到过此针,高中概统,大学概统,随机过程的概率部分,今天读Shannon的传记,他女儿不小心把针洒地上了. 问题描述: 在平面中有一组间隔为a的平行线,随 ...
最新文章
- java:final,finally,finalize区别
- 栈应用_检测成对符号是否正确使用(代码、分析、汇编)
- 单片机小精灵t2_搭建S5P4418 ARM环境下 GPU OPENGL ES开发环境(适用 NANOPI2,3,M2,M3,T2,T3)...
- WebApp的前端所遇问题
- FLEX Array和ArrayCollection的区别
- kafka调试工具kafkacat的使用
- python权威指南 pdf_Python高级编程(第2版) 中文完整pdf扫描版[76MB]
- 传媒大学媒体中心资源批量获取工具的制作
- 机器学习基础 | 有监督学习篇
- java之多态(一)
- java银行存款案例
- 计算机三级考点2:管理和运营宽带城域网的关键技术
- OpenGL学习之旅(三)——Bezier曲线理论与OpenGL实现
- 三天搞定射频识别技术(一)1.3
- 内部存储器——③主存储器
- Failed to check the status of the service. No provider available for the service【已解决】
- 阿里云学生成长计划资格考试分享
- Mac上php环境配置
- 比越南妹划算 1元团购美女
- 完美解决Win10“无法登陆到你的账户”问题,无法登录账户的全方面解决方案!
热门文章
- 程序员多次下载的浏览器竟然是它!
- 报错Found existing installation: tensorflow 1.2.1
- HCI 基础知识之一
- 我们项目组的Team Building
- IT公司内推 互联网内推 bat内推
- 【牛客】CPU的运算速度与许多因素有关,下面______是提高速度的有效措施?
- 【古代文学论文】沈德潜诗学思想的调和格调与文化意蕴(节选)
- 如何区分精确率(precision)、准确率(accuracy)和召回率(recall)
- scikit-learn功能大全
- 向云再出发:如数据般飞驰的内蒙古