[实验六随机数问题1.docVIP

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

实验六 随机数1.随机数的之间的随机数,通常用表示,即 , 为的概率分布密度,记为。 2.产生区间(0,1)内均匀分布的随机数位数的数字自乘以后,得到一个位数的数字,然后取其中间的是位数,便得到一个新的位数,将新的位数的数字再一次自乘以后,又得到一个位数的数字,然后再取其中间的是位数,便得到第三个新的位数,如此不断进行下去,便得到了一个位字长的随机数序列,而将这个数列中的每个数都除以最大的位整数,便得到了一个分布在之间的伪随机数列。 3.产生区间(0,1)内均匀分布的随机数的公式是:, 。 其中,是乘因子,是模数或。第一式称做以为模数的同余式,即以除后得到的余数记为。,均为之间的整数。从上述构造过程可知(因而对也同样)最多有个相异值, 即这表明()有周期,且。因此,不是真正的随机数序列,但若充分大,在一个周期内能经得起统计检验。周期的长度取决于参数的选择时,。 时, 混合同余法: 递推公式为:, 。 其中,是非负整数。通过适当选取参数可以改善随机数的统计性质。一般取作小于的任意奇数正整数,最好使其与模互素。其他参数的选择与乘同余法相类似,  对随机数的统计检验。  (1)参数检验:子样平均值和理论平均值的差异性检验。  在(0,1)上均匀分布随机变量和的平均值和方差分别为:            如果随机数是X的N个独立观测值,令:            则它们的平均值和方差为:            由中心极限定理知统计量:            渐近服从正态N(0,1)。 给定显著性水平,根据正态分布确定临界值,据此判断与, 与的差异是否显著,从而决定能否把 看成(0,1)上均匀分布的随机变量X的N个独立样本值。   (2)分布均匀性检验: 均匀性检验又称频率检验,它用来检验用随机数(样本值)确定的经验频率和均匀分布频率是否有显著性差异。常用的统计检验方法有检验和累积频率检验(K-S检验)。   (3)独立性检验: 只要是检验随机数中前后各数的统计相关性是否显著。 #include iostream #include ctime #include cstdlib using namespace std; int main() { double random(double,double); srand(unsigned(time(0))); for(int icnt = 0; icnt != 10; ++icnt) cout No. icnt+1 : int(random(0,10)) endl; return 0; } double random(double start, double end) { return start+(end-start)*rand()/(RAND_MAX + 1.0); } /* 运行结果 * No.1: 3 * No.2: 9 * No.3: 0 * No.4: 9 * No.5: 5 * No.6: 6 * No.7: 9 * No.8: 2 * No.9: 9 * No.10: 6 */ 伪随机数生成程序: //rand01.c   #include   static unsigned int RAND_SEED;   unsigned int random(void) 伪随机数   { 伪随机数RAND_SEED=(RAND_SEED*123+59)%65536;   return(RAND_SEED);   }   void random_start(void)   {   int temp[2];   movedata(0x0040,0x006c,FP_SEG(temp),FP_OFF(temp),4);   RAND_SEED=temp[0];   }   main()   {   unsigned int i,n;   random_start();   for(i=0;i10;i++)   printf(%u\t,random());   printf(\n);   }

文档评论(0)

84537592 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档