蒙特卡罗方法B.ppt

第三章 由已知分布的随机抽样 f(x) →Xf §3.1 随机抽样的基本方法 直接抽样法 离散分布的直接抽样 连续分布的直接抽样 例3-3 (a) f(x)=2x clear all; N=10000; for i=1:N x=rand; cfx=x; while randcfx x=rand; cfx=x; end X(i)=x; end hist(X) §3.2 复合分布的抽样 加抽样方法 乘抽样方法 减抽样方法 乘加抽样方法 乘减抽样方法 §3.3 变量变换法抽样 变量变换方法 变量变换的分布密度函数 §3.4 随机抽样的其他方法 近似抽样方法 偏倚抽样方法 多维分布抽样方法 §3.5 粒子位置分布的随机抽样 圆均匀分布 环均匀分布 球均匀分布 球壳均匀分布 指数分布（自由程分布） §3.6 粒子方向分布的随机抽样 平行束入射角余弦分布 各向同性散射方位角余弦分布 各向同性散射散射角余弦分布 质心系各向同性的实验室系散射角余弦分布 作业 作分布的抽样: 2. φ为[0,2π]均匀分布，求f(cosφ)。 求 ξ1/ξ2 的分布。 石墨的密度为2.3g/cm3，中子与 C 原子核的作用截面为4.5b（1b=10-24cm2），求中子在石墨中的自由程的抽样。 第四章 计数的统计分布 §4.1 二项式分布 公式的推导 计数的期望值和方差 §4.2 泊松分布 N01, p1 §4.3 高斯分布 m = N0*p 1 * 1 0.999924 0.912 0.899 F(xi) 0.999924~1 0.0076% e+ e– 4 0.912~0.999924 8.79% π0 γ 3 0.899~0.912 1.3% π+ π– 2 0~0.899 89.9% π+ π– π0 1 F(xi) 分支比 f(xi) 衰变方式 i 例3-1：ω介子衰变 clear all; N=10000; F1=0.899; F2=0.912; F3=0.999924; F4=1.0; for i=1:N y=rand; if yF1 n(i)=1; elseif yF2 n(i)=2; elseif yF3 n(i)=3; else n(i)=4; end end m=hist(n,1:4) clear all; N=10000; for i=1:N z(i)=rand^(1/2); end hist(z) 例3-2：连续分布的直接抽样 clear all; N=10000; for i=1:N z(i)=rand^(1/3); end hist(z) clear all; N=10000; lambda=10; for i=1:N z(i)=-1/lambda*log(rand); end hist(z,40) 二. 舍选抽样法 clear all; N=10000; for i=1:N x=10*rand-5; cfx=exp(-x^2/2); while randcfx x=10*rand-5; cfx=exp(-x^2/2); end X(i)=x; end hist(X,40) (b) 正态分布舍选抽样 clear all; N=20000; for i=1:N x=rand; y=rand; if x=0.25 x1=y^(1/2); else x1=y^(1/3); end z(i)=x1; end hist(z,30) 例3-4 (a) 加分布 clear all; N=10000; for i=1:N x1=rand; x2=rand; y=-log(1-x1); while x2exp(-(y-1)^2/2) x1=rand; x2=rand; y=-log(1-x1); end x3=rand; if x3=0.5 x=y; else x=-y; end z(i)=x; end hist(z,30) (b)乘分布 clear all; N=10000; for i=1:N x=rand; x1=rand; y=sqrt(x); while