数学建模社区-数学中国

标题: 请大家帮忙看看这样取随机数可以吗 [打印本页]

作者: world_xin    时间: 2009-6-8 11:25
标题: 请大家帮忙看看这样取随机数可以吗
因为要同时抽取好多组的随机数,但又不希望用同一个种子文件- a* t6 _$ m' g3 X) S
因此程序如下,请大家帮忙看看这样做有没有问题,多谢, q  @7 U  q  F8 F& r2 Z. z

2 ^9 n- a1 y( ]2 c3 ~, {# Q for (i=0;i<bin;i++)+ w- }. L+ l6 e" C6 ~/ J2 ~
   {
/ D' ~# r9 w' ]6 ~- r! l9 p    srand((int)time(0)+1000*i);! T/ m  g2 h! p& ~5 w! H0 t
   for (j=0;j<nn;j++); c( K* @/ ^2 I$ o8 H$ t- _
   {  
/ Q" V( R0 g& ?; ]& R; v9 I6 @4 b      x = (long double)rand()/RAND_MAX;$ Q! K* B6 G$ b* l* w) S7 }
      y = (long double)rand()/RAND_MAX*1.574212218771126;
' q  v* C. J6 E* F! e4 f# B% N      if(y<profile(x))
6 ~/ S8 t( c$ y- P- c      fprintf(f,"%20.15Lf\n",(x+i)*period);
# h( d3 U: A1 y" C   }
9 `9 P+ e% O% c, q; Z& v/ ^9 N   }
作者: madio    时间: 2009-9-3 23:35
random和 rando多可以产生随机数,包含在stdlib.h里。
- o: r6 C; F$ I/ n  random函数不是ANSI C标准,不能在gcc,vc等编译器下编译通过。但在C语言中int random(num)可以这样使用,它返回的是0至num-1的一个随机数。 可改用C++下的rand函数来实现。! w8 w. r* d) @7 H3 p
  1、C++标准函数库提供一随机数生成器rand,返回0-RAND_MAX之间均匀分布的伪随机整数。 RAND_MAX必须至少为32767。rand()函数不接受参数,默认以1为种子(即起始值)。随机数生成器总是以相同的种子开始,所以形成的伪随机数列也相同,失去了随机意义。(但这样便于程序调试)
1 j# z- d* P( j! w4 u  2、C++中另一函数srand(),可以指定不同的数(无符号整数变元)为种子。但是如果种子相同,伪随机数列也相同。一个办法是让用户输入种子,但是仍然不理想。
" S9 P2 D0 X! T9 ~  3、 比较理想的是用变化的数,比如时间来作为随机数生成器的种子。 time的值每时每刻都不同。所以种子不同,所以,产生的随机数也不同。! Z; f5 p! a* ]* s; _; c5 o) |
  rand运用问题与解答:
. B( K; S2 Q- Q' q! a  问:为什么有时rand()用它都是同一个数?如何产生一个真正随机数?" g7 s! u# n+ Z
  答:之所以rand()每次的随机数都一样是因为rand()函数使用不正确。各种编程语言返回的随机数(确切地说是伪随机数)实际上都是根据递推公式计算的一组数值,当序列足够长,这组数值近似满足均匀分布。如果计算伪随机序列的初始数值(称为种子)相同,则计算出来的伪随机序列就是完全相同的。这个特性被有的软件利用于加密和解密。加密时,可以用某个种子数生成一个伪随机序列并对数据进行处理;解密时,再利用种子数生成一个伪随机序列并对加密数据进行还原。这样,对于不知道种子数的人要想解密就需要多费些事了。当然,这种完全相同的序列对于你来说是非常糟糕的。要解决这个问题,需要在每次产生随机序列前,先指定不同的种子,这样计算出来的随机序列就不会完全相同了。你可以在调用rand()函数之前调用srand( (unsigned)time( NULL ) ),这样以time函数值(即当前时间)作为种子数,因为两次调用rand函数的时间通常是不同的,这样就可以保证随机性了。你也可以使用srand函数来人为指定种子数。Windows 9x/NT的游戏FreeCell就允许用户指定种子数,这样用户如果一次游戏没有成功,下次还可以以同样的发牌结果再玩一次。+ ]: O$ w8 b- w" d; q. d
  例子:4 @* l  i0 E+ [# A
  以下为引用的内容:+ x* G, W" C, U) p0 C2 G: L
  #include <stdlib.h>
' }, g: e! ]: d9 z% c$ g! [  #include <iostream.h>
* t2 W3 a2 ~# N( \1 o  #include <conio.h>; E; W7 x& i( j& B: L1 m7 E7 E! t
  #include <time.h>. V5 w  T" P. B
  void main()
- ~! T0 b: x* z  {! v; y$ I8 c1 T6 ^
  srand( (unsigned)time( NULL ) ); //注意不要放在for循环里,否则产生的数是同一个数。8 b6 C- g$ [! I# q: ]
  for(int i=0;i<100000;i++)
+ N& s/ }: n* u1 x0 ~  {1 a0 F7 F5 @! P9 j# C$ X
  cout<<rand()<<endl;# |/ ]5 E" w$ e
  }" F; i& j3 h" I& g* L
  }
作者: 李——建辉    时间: 2012-1-25 20:45
嘿嘿,声明一下:本人看贴和回贴的规则,好贴必看,精华贴必回。5610401069448530913717568348171722629215433605991616398430261146779477050245392




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5