java语言产生泊松分布

1、用 java 语言产生符合泊松分布的随机数一、实验目的利用随机数函数rand来产生的随机数具有一定局限性，rand函数 运用线性同余法产生的是伪随机数，即每次都一样。必须加上时间作 为种子，因此不直接利用 rand 函数二、线性同余法的选取同余法是Lehmer于1951年提出来的，此方法是利用数论中的同 余运算原理来产生随机数，有线性同余法、非线性同余法等，其中线 性同余法法又分为加同余法、乘同余法以及混合同余法。同余法是现 在发展迅速且使用普遍的方法之一。首先介绍下线同余法。线性同余法的递推式为：x = ax + c (mod m)(n = 0,1,2, )n+1n式中参数a，c和m分别称为

2、乘、增量和模，X0为种子。如果这些 参数和种子(初值) 都指定序列也就确定下来了。通常取x=訥=0丄2,)作为区间(0, 1)上均匀分布U (0,1)的随机数。从上 可以看岀当c = 0，为乘同余法；当a =1，c丰0时，为加同余法，否则称 为混合同余法。按惯例，当强调使用某方法产生随机数时，常使用某 方法(随机数)发生器的称呼。线性同余法有如下特点:0 x = c) n+; while (p = c);return u;/线性同余法生成U【0,1)随机数 public static double randomU() seed = (long) (32719 * seed) % (m - 1)

3、;/ System.out.println(seed:+seed);return seed / (m - 1);public static void main(String args) / TODO Auto-generated method stub/ 随机数种子生成与初始化seed = System.currentTimeMillis();m = Math.pow(2, 31); double uu = new double1000;int N=1000;/ 生成10个泊松分布分布随机数for (int i = 0; i N; i+) Sys tem.out.prin tln(“ 第+ i

4、 + 个)；uui = P_rand(0.2); /计算均值与方差 double sum = 0, sum2 = 0, ave = 0, fangca = 0； for (int i = 0； i N； i+) sum = sum + uui；sum2 = sum2 + uui * uui；System.out.println(u: + i + : + uui)； ave = (double) sum / N；Sys tem.o ut .pri ntln(均值:+ ave);for (int i = 0； i N； i+) / sum = sum + uui;sum2 = sum2 + (uu

5、i - ave) * (uui - ave); fangca = sum2 / N;System. out.println(方差:+ fangca);通过随机性检验，得到均值和方差在预想值附近，得到以下实验结果慎Package E導帶口宠慎Package E導帶口宠宙口 .JRE SysSrfYi Librarv ri.aS L.71Hi PoEmmJa怕瓏 phscke 回gna；public c1as& Pat&ion T禽障an afis蛋 Outline”-_ - .： .： ；. - -, L ：_ -ra. in O Con wk- S3呻terrriFmtedx Pion App

6、lhtio町 住许。巾时秆1讯7阴曲常刑binjgvaw思毎 丘01丰324 F牛8点9逞3|u=.eaill09S8&9J3L?&4b ?=0.6811090973174 n=6.0 c=&. fll730753e?7981fl gU-4. B341?5S1 LB4&-49E p*B bS34135EL1E414AB n-B.B C-fl. B-LB7 3B77BS1Su-Q.BOSE5723S1763a9B p-B.H3717G5339939152 3 n-1.3 c-B. BlE73075Sfl77331B S2+u-e. 14611390332976&7-1 p-6.1+5113903

7、329774 rr-0.0 C-&, S18730753S779813u.7MSa：!6i441 酹酩 pP J&6QP轴43】航佻 n-0.fi c-.&tfl7?a7t207?9915 *41U-i3,574E7E盟5韶巧84+ pT刃轉丁胡皓;闘了駅4 “网花 31077530773313 J5+u=0.45394-6119714&4.n=0.fi 口飢吐B75S75307JSSHSSfitu-4.E3B73B9fl3235952 p-a .fl3E7B5i3K23S352B n-B.S c-fl. BK.B7 =7&36775313U-R,072559231141W5 p-0.0563

8、735595175 n-1-0 C-fr. BlS730755fl?73310 J7ru=0.39463fiSi72625673 p-0.594363472625673 n=e.S c=0a313757530775313 JsSi-1.垃33朮阳抽酣1昭 尸葩空抽彷昭阳旳钿n-6.a c-H.Blfl73fl?S3fi779ei2 泗亍u-fi,6197&W2H2533&3 p-?.6197WS2012S3303 n-0-B -6,3107753718u;l ?fi-&445&5729B176jS9au ;1.146U39&3329766U4U!3 :_7eeSBiB4&iJ.Lfl7AEU ；4也巧；朝托T0WR沾妙u：5 ：i, 3043315L1371Wu：6 *.7256929114114925u;7;-e.3946360i?2fi25673u：3:fl.l2S3e-S5a5939S10HU!9 ：0.51979201293303SI：&. 3141771