粒子物理与核物理中的数.pdfVIP

粒子物理与核物理实验中的数 据分析 杨振伟 清华大学 第四讲:蒙特卡罗方法 24/03/2010 1 上一讲回顾 概率的基本概念 随机变量与概率密度函数 随机变量的平均值与方差 能不通过实验对随机变量进行研究吗？ 24/03/2010 2 本讲要点  蒙特卡罗方法  随机数产生子  任意分布抽样之函数变换法与舍选法  蒙特卡罗方法中的精度问题  在粒子物理与核物理中的应用 24/03/2010 3 蒙特卡罗方法简介 蒙特卡罗方法就是利用一系列随机数来计算各种概率大小和 随机变量均值等等的数值分析技术。通常的步骤为： 1) 产生一系列在[0，1]之间均匀分布的随机数 r , r ,..., r。 1 2 m 2) 利用这些随机数按某些概率密度函数f (x) 抽样生成我 们感兴趣的另一随机序列x x ,...,x 。 1, 2 n 3) 利用这些 值来估计f (x) 的一些特性，例如：通过找 x 到在区间 的 比例，给出积分值 b f (x)dx 。 a ≤x ≤b x ∫ a 第一层面上的应用： 蒙特卡罗计算 = 积分 第二层面上的应用： 蒙特卡罗变量 = “模拟的数据” 24/03/2010 4 随机数的产生 用物理方法产生 用数学方法产生 真正的随机数 伪随机数 •不可重复 •可以重复 •产生速度慢 •产生的速度快 24/03/2010 5 真随机数与伪随机数 美国兰德（RAND）公司在1950年代，利用真空管中 产生的噪音制作了一个含十万个真正的随机数表，并运用于 其开展的所有模拟研究中。 真正的随机数与伪随机数之间的区别在于：数据串是否 具有可压缩性，即能否用更短的形式来表示。 真正的随机数是不可压缩的，非常不规则，以至于无法 用更短的形式来表示它。 在粒子物理与核物理研究中，随机数的可重复性经常也 是非常有用的，尤其是程序的调试(debugging)。 24/03/2010