2017春季计算物理第二部分第1讲-北京大学物理学院.PDF

2017春季计算物理第二部分 第5讲 李强 北京大学物理学院中楼411 qliphy0@ /personnel/member/LiQ.xml 5. 随机数及简明概率论 真、伪随机数 产生方法：(1) 取中法, (2) 线性同余产生器, (3) 反馈位移寄存器 C/C++ Random generator CERN Root三种Random generator 随机性的检验: 均匀性，独立性 简明概率论 平均值，方差 中心极限定理 常见分布函数: Gauss, Landau, Chi2分布 置信度 随机性的检验: correlation function 抽样 直接抽样, 变换抽样，取舍取样 加、减、乘抽样 引言: 布冯实验 (ℎ, dℎd ≤ 随机数 伪随机数 伪随机数的产生 取中法 线性同余法 线性同余法 4, 1, 3, 0, 2, 4, 1, 3, …. 11, 38, 11, …. 线性同余法 C++的线性同余法 C/C++的线性同余法 #include cstdlib C/C++的线性同余法 if(fx0.1) { f[0]+=1.0; } #include iostream else if(fx0.2) { f[1]+=1.0; #include cmath } #include cstdlib …… #include ctime else{ using namespace std; f[9]+=1.0; int main() } { mean+=fx; double mean=0.0; variance+=fx*fx; } double variance=0.0; mean=mean/double(NN); double chi2=0.0; variance=sqrt(variance/double(NN)-mean*mean); double f[10];