数学建模社区-数学中国

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

作者: world_xin    时间: 2009-6-8 11:25
标题: 请大家帮忙看看这样取随机数可以吗
因为要同时抽取好多组的随机数,但又不希望用同一个种子文件
$ e7 l& e1 z& b% ^! v因此程序如下,请大家帮忙看看这样做有没有问题,多谢
, `' p0 Z: ]+ o9 Z
1 J4 K& ~# J& |5 q/ f for (i=0;i<bin;i++); V' l$ A! z' x5 v4 P/ h9 G
   {
/ {- [8 S2 b0 \8 }5 _% {    srand((int)time(0)+1000*i);8 J$ V% v$ k0 Q
   for (j=0;j<nn;j++)
3 E) P: R$ l2 ~7 N- `5 y   {  
( c  d% H4 _7 O& z1 g      x = (long double)rand()/RAND_MAX;
% \, H' Q4 j; A2 O0 g5 _$ W9 m      y = (long double)rand()/RAND_MAX*1.574212218771126;
9 y) c$ g) C6 I# s      if(y<profile(x)) 2 J. k1 e% e; q) t0 R$ E1 e
      fprintf(f,"%20.15Lf\n",(x+i)*period);
+ F; r: S# Y+ E% r# R4 L8 T   }
- {/ G- E1 C! d7 W8 m   }
作者: madio    时间: 2009-9-3 23:35
random和 rando多可以产生随机数,包含在stdlib.h里。) \2 |( ]- V0 G0 A
  random函数不是ANSI C标准,不能在gcc,vc等编译器下编译通过。但在C语言中int random(num)可以这样使用,它返回的是0至num-1的一个随机数。 可改用C++下的rand函数来实现。
! N+ S: l# c( h5 B  1、C++标准函数库提供一随机数生成器rand,返回0-RAND_MAX之间均匀分布的伪随机整数。 RAND_MAX必须至少为32767。rand()函数不接受参数,默认以1为种子(即起始值)。随机数生成器总是以相同的种子开始,所以形成的伪随机数列也相同,失去了随机意义。(但这样便于程序调试), t) B4 Y7 J6 N8 c  F
  2、C++中另一函数srand(),可以指定不同的数(无符号整数变元)为种子。但是如果种子相同,伪随机数列也相同。一个办法是让用户输入种子,但是仍然不理想。0 G$ n, H3 h; ^; G  ?- l
  3、 比较理想的是用变化的数,比如时间来作为随机数生成器的种子。 time的值每时每刻都不同。所以种子不同,所以,产生的随机数也不同。
$ F1 K) j( s% _2 O' G5 G  rand运用问题与解答:+ T( X9 E# X9 V. l7 c& ^( c6 G" Q
  问:为什么有时rand()用它都是同一个数?如何产生一个真正随机数?
- f- d! Q( |; O* s* @9 }) P# g  答:之所以rand()每次的随机数都一样是因为rand()函数使用不正确。各种编程语言返回的随机数(确切地说是伪随机数)实际上都是根据递推公式计算的一组数值,当序列足够长,这组数值近似满足均匀分布。如果计算伪随机序列的初始数值(称为种子)相同,则计算出来的伪随机序列就是完全相同的。这个特性被有的软件利用于加密和解密。加密时,可以用某个种子数生成一个伪随机序列并对数据进行处理;解密时,再利用种子数生成一个伪随机序列并对加密数据进行还原。这样,对于不知道种子数的人要想解密就需要多费些事了。当然,这种完全相同的序列对于你来说是非常糟糕的。要解决这个问题,需要在每次产生随机序列前,先指定不同的种子,这样计算出来的随机序列就不会完全相同了。你可以在调用rand()函数之前调用srand( (unsigned)time( NULL ) ),这样以time函数值(即当前时间)作为种子数,因为两次调用rand函数的时间通常是不同的,这样就可以保证随机性了。你也可以使用srand函数来人为指定种子数。Windows 9x/NT的游戏FreeCell就允许用户指定种子数,这样用户如果一次游戏没有成功,下次还可以以同样的发牌结果再玩一次。: F" S9 k% l% O4 G$ E
  例子:% l. Q% y! ]4 J" g: u
  以下为引用的内容:
4 [0 q2 R) B* J" Z  #include <stdlib.h>4 V, {' M0 w/ O; E, C
  #include <iostream.h>
3 J. s: k- l, ]( x, F  #include <conio.h>; p7 z5 F+ v$ D. |6 ]  S
  #include <time.h>
$ P& c2 u7 W9 P) @$ Z$ p% Y& ^$ o  void main()
9 O' u, I3 m. S1 w0 `. v  {
# o. N- ?  K% B! a3 c  srand( (unsigned)time( NULL ) ); //注意不要放在for循环里,否则产生的数是同一个数。
: }  [& v1 K7 s) a9 D  for(int i=0;i<100000;i++)1 ~; z3 _) {1 W9 X4 @6 g, x! x
  {
! D3 l8 |: |: j) r  cout<<rand()<<endl;, [# ?) k/ p0 j3 Z5 A4 P6 D7 {
  }
7 X& A# Q7 [. H: O: W: U  }
作者: 李——建辉    时间: 2012-1-25 20:45
嘿嘿,声明一下:本人看贴和回贴的规则,好贴必看,精华贴必回。5610401069448530913717568348171722629215433605991616398430261146779477050245392




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