离散系统仿真—随机数生成及检验(韩宝安).docxVIP

下载本文档

20
0
约5.78万字
约 14页
2016-08-19 发布于湖北
举报
版权申诉

离散系统仿真—随机数生成及检验(韩宝安).docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

离散系统仿真—随机数生成及检验(韩宝安)

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(拒绝原假设，即所生成的随机数不符合要求！);System.out.println(-----------------------------);System.out.println(均匀性分析:);int Segment=200;int N[]=newint[Segment];double ka=0;double X2=232.912;//卡方分布的分位点值int Num;for(int i=1;i=Segment;i++){Num=0;for(int j=1;j=n;j++){if(Un[j]=(1.0/Segment*(i-1)) Un[j]=1.0/Segment*i){Num++;}}N[i-1]=Num;ka=ka+Math.pow(N[i-1]-n/Segment,2)/(n/Segment);}System.out.print(各区间随机数个数：);for(int i=0;iSeg