数学建模社区-数学中国
标题:
请大家帮忙看看这样取随机数可以吗
[打印本页]
作者:
world_xin
时间:
2009-6-8 11:25
标题:
请大家帮忙看看这样取随机数可以吗
因为要同时抽取好多组的随机数，但又不希望用同一个种子文件
7 V# S9 @. H* L+ Y- Z4 t- c
因此程序如下，请大家帮忙看看这样做有没有问题，多谢
- K' \! s t. g! O$ e1 Z8 G
8 D: N# L/ d/ W: v. e
for (i=0;i<bin;i++)
! G5 i/ ]9 ~2 z) P$ H# m
{
! W; x, m6 ]6 w2 J& p
srand((int)time(0)+1000*i);
& t; u r: Y8 T4 d3 L) G
for (j=0;j<nn;j++)
3 I, F( p% }7 I' D0 \1 b
{
; \4 V% ?5 w) K, Q" R: {
x = (long double)rand()/RAND_MAX;
- v/ x9 l) ? s
y = (long double)rand()/RAND_MAX*1.574212218771126;
" D4 f6 g2 f9 u3 F8 P0 b
if(y<profile(x))
6 @1 V$ u% @/ d+ b+ \% s
fprintf(f,"%20.15Lf\n",(x+i)*period);
K" F, J5 ~$ J: L: L" a/ ]! z
}
k: a0 S# s" x
}
作者:
madio
时间:
2009-9-3 23:35
random和 rando多可以产生随机数，包含在stdlib.h里。
9 Z0 m6 K$ x" C# q* f
random函数不是ANSI C标准，不能在gcc,vc等编译器下编译通过。但在C语言中int random(num)可以这样使用，它返回的是0至num-1的一个随机数。 可改用C++下的rand函数来实现。
0 ?) F Z. a a9 d- F
1、C++标准函数库提供一随机数生成器rand，返回0－RAND_MAX之间均匀分布的伪随机整数。 RAND_MAX必须至少为32767。rand()函数不接受参数，默认以1为种子（即起始值）。随机数生成器总是以相同的种子开始，所以形成的伪随机数列也相同，失去了随机意义。（但这样便于程序调试）
# K) Z8 @$ J" f) G$ [2 B5 w
2、C++中另一函数srand（），可以指定不同的数（无符号整数变元）为种子。但是如果种子相同，伪随机数列也相同。一个办法是让用户输入种子，但是仍然不理想。
, c5 k: @3 G9 F9 J% N4 T; z
3、 比较理想的是用变化的数，比如时间来作为随机数生成器的种子。 time的值每时每刻都不同。所以种子不同，所以，产生的随机数也不同。
+ K4 b4 F, ?: @$ {/ V
rand运用问题与解答：
/ J3 n# H+ a$ M3 k! \+ l
问：为什么有时rand()用它都是同一个数？如何产生一个真正随机数？
" {( q4 l+ L& q. o
答：之所以rand()每次的随机数都一样是因为rand()函数使用不正确。各种编程语言返回的随机数（确切地说是伪随机数）实际上都是根据递推公式计算的一组数值，当序列足够长，这组数值近似满足均匀分布。如果计算伪随机序列的初始数值（称为种子）相同，则计算出来的伪随机序列就是完全相同的。这个特性被有的软件利用于加密和解密。加密时，可以用某个种子数生成一个伪随机序列并对数据进行处理；解密时，再利用种子数生成一个伪随机序列并对加密数据进行还原。这样，对于不知道种子数的人要想解密就需要多费些事了。当然，这种完全相同的序列对于你来说是非常糟糕的。要解决这个问题，需要在每次产生随机序列前，先指定不同的种子，这样计算出来的随机序列就不会完全相同了。你可以在调用rand()函数之前调用srand( (unsigned)time( NULL ) )，这样以time函数值（即当前时间）作为种子数，因为两次调用rand函数的时间通常是不同的，这样就可以保证随机性了。你也可以使用srand函数来人为指定种子数。Windows 9x/NT的游戏FreeCell就允许用户指定种子数，这样用户如果一次游戏没有成功，下次还可以以同样的发牌结果再玩一次。
% b) r& L5 @6 b
例子：
0 \, {/ J! v; z- y2 x6 F9 ]7 Z
以下为引用的内容：
4 z0 X2 `- y: u, u( H0 I& C b7 b
#include <stdlib.h>
3 R* B- `! q E
#include <iostream.h>
# K( I; ^6 r) ]) W
#include <conio.h>
* D* H, ]# V0 Q2 }) _" K: M3 L
#include <time.h>
. Z) }' @! u5 D; c- |: z+ w
void main()
6 h7 i/ I2 T' L8 l; |3 v1 @' e
{
8 o" H# y! p( \& u T1 j# C
srand( (unsigned)time( NULL ) ); //注意不要放在for循环里，否则产生的数是同一个数。
. l Q8 G- [7 i4 V
for(int i=0;i<100000;i++)
5 }; L/ b" E4 ~2 l( K& O
{
7 g( b4 t' U% A8 M
cout<<rand()<<endl;
. f2 M; l( I* Z# l P& ^7 Q; X
}
& m0 F, ]) E4 {
}
作者:
李——建辉
时间:
2012-1-25 20:45
嘿嘿，声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。
5610401069448530913717568348171722629215433605991616398430261146779477050245392
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5