完整版统计建模与R软件第十讲-(2017).pptVIP

(1) 产生随机数。产生n个相互独立的随机变量抽样序列 (2) 模拟实验，检验不等式x =a/2 * sinα 设上面的实验有k次成功的，则 精选文档 * buffon - function(n, L=1, a=0.8) { alpha = runif(n, 0, pi); x = runif(n, 0, L/2); k -0; for (i in 1:n) { if (x[i] = a/2*sin(alpha[i])); k = k+1; } 2*a*n/(k*L); } 精选文档 * 随机模拟计算的思路： [1] 针对实际问题建立一个简单且便于实现的概率统计模型，使所求的解恰好是所建模型的概率分布或某个数字特征，比如，是某个事件的概率或该模型的期望值； [2] 对模型中的随机变量建立抽样方法，在计算机上进行模拟试验，抽取足够的随机数，并对有关的事件进行统计； [3] 对模拟试验结果加以分析，给出所求解的估计及其精度（方差）的估计； [4] 必要时，还应改进模型以降低估计方差和减少试验费用，提高模拟计算的效率。 精选文档 * 用蒙特卡洛方法进行计算机模拟的步骤： [1] 设计一个逻辑框图，即模拟模型．这个框图要正确反映系统各部分运行时的逻辑关系。 [2] 模拟随