数学实验-实验报告-概率与频率

数学实验报告实验序号：8日期：6/5班级信科姓名学号实验名称概率与频率问题背景描述：概率，又称为几率、或然率，是反映某种事件发生的可能性大小的一种数量指标。它介于0和1之间。这里的事件是指随机现象中出现的某个可能结果。实验目的：概率论是研究随机现象统计规律的一门数学分支学科，它有着悠久的历史。通过本实验的学习，加深对频率和概率等概念的理解和认识，并帮助掌握一些概率统计的原理。实验原理与数学模型：相关函数（命令）简介1．rand(m,n):生成m×n的随机矩阵,每个元素都在(0,1)间,生成方式为均匀分布2．randn(m,n):生成m×n的随机矩阵,每个元素都在(0,1)间,生成方式为正态分布3.randperm(m)：生成一个1~n的随机整数排列。4．perms(1:n)：生成1到n的全排列，共n！个。5．一系列取整的函数：(1)fix(x)：截尾法取整；(2)floor(x)：退一法取整(不超过x的最大整数)；(3)ceil(x)：进一法取整(=floor(x)+1)；(4)round(x)：四舍五入法取整。6．unique(a)：合并a中相同的项。7．switch表达式case情况1命令系列1case情况2命令系列2……otherwise命令系列end8．prod(x)：向量x的所有分量元素的积。9.unidrnd(n)生成一个1到n的随机整数。实验所用软件及版本：Matlab2009主要内容（要点）：1．通过实验，填写完成表格2~6的数据3．用MonteCarlo方法求两平面曲线y=x2x≥0与y=1-x2及y4．分析附录中的[程序丙]和[程序丁]的设计本意。请问他们为什么都是错误的？5．设计一个三维投点的蒙特卡罗法计算π。并比较运行结果与二维投点的蒙特卡罗法的运行结果，哪个更准确些。提示：随机投点落在单位正方体的内切球体内部。实验过程记录（含基本步骤、主要程序清单及异常情况记录等）：1．通过实验，填写完成表格2~6的数据实验1：随机投掷均匀骰子，验证各点数出现的概率是否为1/6表2试验次数/n100001000010000100001000010000国徽朝上频率0.49680.50780.49360.49990.50070.5004国徽朝下频率0.50310.49780.49910.49430.50170.5019从实验结果我们可以看出[程序乙]的误差要小很多，所以我们有理由认为[程序乙]正确，另一方面，分析[程序甲]和[程序乙]的不同之处：（1）[程序甲]没有分别用变量x和y事先定义rand(1)*a和rand(1)*b（2）[程序甲]的if条件句：rand(1)*b>=(rand(1)*a)^2&rand(1)*b<1-(rand(1)*a)^2[程序乙]的if条件句：y<=1-x^2&y>x^2即：rand(1)*b<=1-(rand(1)*a)^2&rand(1)*b>(rand(1)*a)^2

可以看出[程序甲]和[程序乙]的取等情况及不等式的顺序不同，不过很显然，这两种逻辑并不影响实验结果。经过分析[程序甲]和[程序乙]的不同之处我们可以认为，由于[程序甲]没有用变量x和y事先定义rand(1)*a和rand(1)*b而引起甲乙两结果不同，所以Monte

Carlo投点法在使用过程中应事先定义，再进行if语句的运行。4．分析附录中的[程序丙]和[程序丁]的设计本意。请问他们为什么都是错误的？[程序丙]结果[程序丁]结果通过分析对比[程序丙]和[程序丁]与[程序乙]的区别，我们可以看出：[程序丙]的a的赋值是错误的，曲线y=x2x≥0与y=1-x2的交点横坐标为[程序丁]不仅没有事先定义rand(1)*a和rand(1)*b，而且[程序丁]的if条件句rand(1)<1-rand(1)^2&rand(1)>=rand(1)^2也是错误的，rand(1)没有乘以a或b，使得结果偏小很多。5．设计一个三维投点的蒙特卡罗法计算π。并比较运行结果与二维投点的蒙特卡罗法的运行结果，哪个更准确些。提示：随机投点落在单位正方体的内切球体内部。试验次数n100000100000100000100000100000100000（二维）所得π的近似值3.13843.14523.13823.13853.14223.1321（三维）所得π的近似值3.14833.12673.13983.12953.14523.1216通过对比二维与三维投点的蒙特卡罗法的运行结果可以发现，二维投点的蒙特卡罗法的运行结果更加准确。实验结果报告与实验总结：通过本实验加深了我们对频率和概率等概念的理解和认识，而且我们可以体会到运用经典的蒙特卡罗投点法可以近似求解无理数π或是不规则曲面面积等，从频率与概率的角度来解决数学问题也是一个很好的思路。思考与深入：本次实