实验一随机数的产生与Monte-Carlo方法(11本系).docVIP

随机数的产生与Monte Carlo方法 姓名 龚智强 学号 5090719077 一、实验目的 1．掌握随机数的产生； 2．应用Monte Carlo方法解决一些简单问题； 3．巩固概率统计知识 二、实验原理 A．随机数的产生 1．定义：设随机变量，来自(的样本，称为来自分布的随机数。 2．定理：设X的分布函数为F（x）， 严格递增， 1）； 2）若, 则的分布函数为F（x）。 3?．产生均匀随机数 常用(伪随机数)具有独立性、均匀性、随机性； * 线性同余发生器 设a,b,M为整数，， 求余运算 *三个较满意的发生器 为任意非负整数。 4．正态分布随机数 若，，令 ，则 ，，且U与V相互独立。 若，则。 5．指数分布随机数 若，则，也可以。 B．计算机随机模拟方法，也称Monte Carlo方法，是计算机仿真技术中的一个有效方法。 步骤： 建立一个概率模型，使它的参数等于问题的解； 应用计算机产生随机数来计算所求参数； 得到所求解的近似解。 用Monte Carlo方法计算常义定积分 三、实验内容 初始值取学号最后两个数字。 初始值为77。 （一）．1．打印10个来自随机数； 0.475789 0.823194 0.796674 0.183235 0.250142 0.0182719 0.204703 0.227195 0.926179 0.094057 2．打印10个来自随机数； 1.11213 0.0945271 -0.604201 -1.76362 1.52042 2.82927 0.204034 1.65234 0.387618 -0.972744 3．打印10个来自随机数如下，取x0=5,a=5/4。要注明分布。 6.39016 5.11816 4.24475 2.79547 6.90053 8.53659 5.25504 7.06542 5.48452 3.78407 4．打印10个来自随机数，取=1/25。要注明分布。 16.1465 43.3175 39.8276 5.0601 7.19675 0.461023 5.72599 6.4432 65.1529 2.46947 （二）．应用Monte Carlo方法解决一、二个简单问题。 F(x)= 在区间（3,4）上的积分实际为12.333333； F(x)=cos（x）在区间（3,4）上的积分实际为-0.897923. 1.取10000个随机数 F(x)= 在区间（3,4）上的积分为12.3101 F(x)=cos(x)在区间（3,4）上的积分为-0.899274 2.取100000个随机数 F(x) = 在区间（3,4）上的积分为12.3076 F(x)=cos(x)在区间（3,4）上的积分为-0.899186 随机数 F(x) = 在区间（3,4）上的积分为12.3084 F(x)=cos(x)在区间（3,4）上的积分为-0.899187 四、实验小结 利用公式求取U（0,1）随机数已经比较接近实际，虽有误差，但是通过Monte Carlo方法解决的问题解和实际解非常接近误差0.14%，因此在误差允许的范围内可以利用随机数产生公式来生成随机数。