分布模拟 monter carlo(蒙特卡洛) 标准正态分布 泊松分布.docVIP

老师，我模拟用的是matlab，使用了两种不同的蒙特卡洛方法模拟了标准正太分布，用了不同的第三种蒙特卡洛方法模拟了泊松分布（期望为10),以下是模拟的图。 同时认为暗电流服从泊松分布（取期望为4），模拟了一下。 以下是模拟图： 第一种 方法的标准正态分布： 第二种方法： 泊松分布（期望10） 暗电流取 期望为4 正太分布的两种数据产生方法： 方法一： j=1; k=1; for i=1:200 a=rand; b=rand; w=(2*a-1)^2+(2*b-1)^2; if(w1) z=(-2*log(w)/w)^(1/2); x(j)=a*z; y(k)=b*z; j=j+1; k=k+1; end End 其中产生的x，y向量存的就是符合正太分布。 方法二： j=1; for i=1:2000 a=rand; L=(2*exp(1))/pi; b=-log(a); c=rand; if((b-1)^2=-2*log(c)) norm1(j)=b; j=j+1; end End Norm1中存的就是符合正态分布的。 泊松分布： j=1; a=0;b=25; max=0.3;mean=10; r2=randint(1,20000,[a b]); for i=1:20000 r1=rand; r3=rand; if (r3=max* mean^r2(i)*exp(-mean)/factorial(r2(i))); po_data(j)=r2(i); j=j+1; end end po_data中存的就是符合泊松分布的数据。