数学建模社区-数学中国
标题:
请大家帮忙看看这样取随机数可以吗
[打印本页]
作者:
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