离散系统仿真—随机数生成.docx

下载文档 降价啦

3
0
约5.61万字
约 12页
2016-08-19 发布于湖北
举报
版权申诉
保障服务

离散系统仿真—随机数生成.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

离散系统仿真—随机数生成

import java.util.*;publicclass RanddomNum {/***@paramargs*/publicstaticvoid main(String[] args) {// TODO Auto-generated method stub//随机数发生器LCGint n=1000;//循环次数int b=31;//模数指数double Xn[]=newdouble[n+1];//定义XnXn[0]=10.0;//初始值double Un[]=newdouble[n+1];//定义Unlong a=214159269;//乘子long c=453806245;//增量double m=Math.pow(2, b);//模数double Sum=0.0;//样本和值double u,S2=0.0;//样本均值和方差System.out.println(随机数发生器生成的随机数如下：);for(int i=0;in;i++){Xn[i+1]=(a*Xn[i]+c)%m;Un[i+1]=Xn[i]/m;System.out.println(第+(i+1)+个随机数：++Un[i+1]);Sum+=Un[i+1];}u=Sum/n;//计算样本均值System.out.println(样本均值u=+u);for(int i=0;in;i++){Un[i+1]=Xn[i]/m;S2+=Math.pow(Un[i+1]-u,2)/(n-1);//计算样本方差}System.out.println(样本方差S2=+S2);//统计分析double E=1.0/2;//总体均值，亦即样本均值期望double Var=1.0/12;//总体方差double Varu=1.0/(12*n);//样本方差double ES2=Var;//样本均值方差的期望double VarS2=1.0/(180*n);//样本均值方差的方差double confidence_level=0.95;//置信度double v1=(u-E)/Math.sqrt(Varu);//样本均值统计量double v2=(S2-ES2)/Math.sqrt(VarS2);//样本方差统计量double Z=1.96;//Z1-confidence_level/2的值if(Math.abs(v1)=Z Math.abs(v2)=Z)System.out.println(不拒绝原假设);elseSystem.out.println(拒绝原假设，即所生成的随机数不符合要求！);System.out.println(-----------------------------);Sum=0.0;S2=0.0;System.out.println(Java中自带随机函数生成的随机数如下：);for(int i=0;in;i++){Un[i+1]=Math.random();System.out.println(第+(i+1)+个随机数：++Un[i+1]);Sum+=Un[i+1];}u=Sum/n;//样本均值System.out.println(样本均值u=+u);for(int i=0;in;i++){Un[i+1]=Xn[i]/m;S2+=Math.pow(Un[i+1]-u,2)/(n-1);}System.out.println(样本方差S2=+S2);v1=(u-E)/Math.sqrt(Varu);//样本均值统计量v2=(S2-ES2)/Math.sqrt(VarS2);//样本方差统计量if(Math.abs(v1)=Z Math.abs(v2)=Z)System.out.println(不拒绝原假设);elseSystem.out.println(拒绝原假设，即所生成的随机数不符合要求！);}}结果:随机数发生器生成的随机数如下：第1个随机数：4.6566128730773926E-9第2个随机数：0.2085768091492355第3个随机数：0.18907278031110764第4个随机数：0.630544401705265第5个随机数：0.35256192088127136第6个随机数：0.4644889086484909第7个随机数：0.34608858823776245第8个随机数：0.2775612026453018第9个随机数：0.47259870171546936第10个随机数：0.7010539621114731第11个随机数：0.26666679978370667第12个随机数：0.11956797540187836第13个随机数：0.4191962517797947第14个随机数：0.060020774602890