本科学生设计性实验报告(样本).docxVIP

本科学生设计性实验报告 项目组长 华南虎＿ 学号＿ 007（写全） 成 员 华北虎 华东虎 胡戈 专 业 信息与计算科学专业 班级 07 信本 1 班 实验项目名称＿案例分析—— 蒲丰投针问题 指导教师及职称＿曹慧荣＿副教授＿ ＿ 开课学期 2013 至＿ 2014 学年＿一 ＿学期 上课时间 2013 年 11 月 14 日 1 案例分析——蒲丰投针问题 一、摘要 法国科学家蒲丰提出的一种计算圆周率的方法 —— 随机投针法， 即著名 的蒲丰投针问题。本实验给出了蒲丰投针实验 matlab 实现，并且在此基础 上给出圆周率 matlab 的近似计算。 二、实验目的及要求 选择一个典型的问题，进行数学模型的建立及数学理论分析，给出 matlab 程序的实现。 三、实验仪器设备 计算机 4 台 四、实验方案设计 （一）原理描述 1．蒲丰投针问题 1777 年法国科学家蒲丰提出的一种计算圆周率的方法 —— 随机投针法， 即著名的蒲丰投针问题。这一方法的步骤是： step1 ：取一张白纸，在上面画上许多条间距为 d 的平行线； step2 ： 取一根长度为 l （ ld ） 的针，随机地向画有平行直线的纸上 n 次，观察针与直线相交的次数，记为 m； step3 ：计算针与直线相交的频率。 18 世纪，法国数学家布丰和勒可莱尔提出的 “投针问题 ”，记载于布丰 1777 年出版的著作中： “在平面上画有一组间距为 d 的平行线，将一根长度 为 （l ld ）的针任意掷在这个平面上， 球此针与平行线中任一条相交的频率。 ” 布丰本人证明了，这个概率是 p=2l/( πd) （π为圆周率）。利用这个公式可 以用概率的方法得到圆周率的近似值。 2．流程图 根据以上分析，下面给出计算圆周率的近似值的流程图： 2 赋初值（针长、 平行线间距离、 试验总次数等） 循环体（针与线相交时进行累加） 计算针与线相交的频率 计算圆周率的近似值 （二）实验过程设计 step1 ：理解蒲丰投针问题； step2 ：给出圆周率近似计算的公式； step3 ：利用 matlab 模拟蒲丰投针问题。确定近似计算过程中涉及的变 量；两个随机变量的产生方法：给出针的中点到最近平行线的距离 d、针与 线相交的角度的产生方法。 step4 ：给出针与线相交的频率。 step5 ：估计圆周率的值。 （三）实验观测点、观测指标与假设条件 本实验的主要观测点及观测指标是： 针长与平行线相交的次数及频率的 计算。 本实验的假设条件：针长小于平行线的距离。 五、实验内容及步骤 （一）实验调试步骤 1、蒲丰投针问题中圆周率 值的估计 在画有许多间距为 d 的等距平行线的白纸上， 随机投掷一根长为 l ( l d ) 的均匀直针，求针与平行线相交的概率，并计算的 值。 设针与平行线的夹角为 ( 0 ) ，针的中心与最近直线的距离为 3 x(0 x d / 2) ，则针与平行线相交的充要条件为 x l sin 。 2 利用几何概率所以针与平行线相交的概率为 2、蒲丰投针问题的 matlab 实现 根据以上分析给出蒲丰投针问题的 matlab 实现。 在 matlab7.0 中，建立以下脚本文件 n=100000;% 总的试验次数 l=1; % 针的长度 d=2; % 平行线间的距离 m=0; % 针与线相交次数的初始值 for k = l : n % 循环 n 次 x = d/2*rand; % 随机投针，针的中点与离最近平行线间的距离 phi = pi*rand; % 随机产生针与平行线间的交角 0=phi=pi if x 0.5*l*sin(phi) m = m + 1; end end p=m/n; 3、圆周率值的估计 由蒲丰投针问题，给出的相交频率 p, 估计圆周率的值，在上边的 matlab 脚本文件中加入命令 pi_m=2*l/(p*d) 运行，即得结果． 4 当实验总次数越大时，越接近真值。 N=1000:200:100000; chang=length(N); PI=zeros(1,chang); for j=1:chang n=N(j);% 总的试验次数 l=1; % 针的长度 d=2; % 平行线间的距离 m=0; % 针与线相交次数的初始值 for k = 1 : n % 循环 n 次 x = d/2*rand; % 随机投针，针的中点与离最近平行线间的距离 phi = pi*rand; % 随机产生针与平行线间的交角 0=phi=pi if x 0.5*l*sin(phi) m = m + 1; end end p=m/n; pi_m=2*l/(p*d); PI(j)=pi_m; end t=1:length(PI); plot(t,pi*