用C语言模拟随机数的生成

函数实现出处：http://baike.baidu.com/view/876758.htm

ftime讲解出处：百度百科 ftime

函数定义：int ftime(struct timeb *tp);

函数说明：ftime()将目前日期由tp所指的结构返回。
tp结构定义：
struct timeb{

　　time_t time; /* 为1970-01-01至今的秒数*/

　　unsigned short millitm; /* 千分之一秒即毫秒 */

　　short timezonel; /* 为目前时区和Greenwich相差的时间，单位为分钟 */

　　short dstflag; /* 为日光节约时间的修正状态，如果为非0代表启用日光节约时间修正 */

　　};

　　返回值 ：无论成功或失败都返回0

在函数里，我增加了一些注释，希望大家能看得懂 。

有些人可能会觉得为什么不直接用Initial()函数的返回值就好了，这是因为Initial()函数里是通过 timeb作为种子，因为 timeb 是返回现在与1970.1.1午夜后的时间差，单位是毫秒，而计算机的执行速度很快，执行我们的这个小程序甚至不需要一毫秒，所以给出的数值都会很接近，甚至相同。

该程序的实现也是一个伪随机数。

事实上，在查找资料时，还发现有有另外一种实现方法，该方法较安全，且随机性较大，但较为复杂。有兴趣的可以看一下这个帖子：http://bbs.csdn.net/topics/80515171

#include<stdio.h>
#include<sys/timeb.h>
#include<time.h>
#include<stdlib.h>

//用于返回随机数初值
double Initial()
{
double init;
struct timeb tp;
while(true)
{
ftime(&tp);
init=tp.millitm*0.9876543*0.001;
//tp.millitm的取值范围是	0-999
//因此这边init的取值应该是	[0，0.987654321）
if(init>=0.001)
{
break;
}
}
return init;
}

//返回一个(0 , 1)随机数
double Random()
{
static double rndm=-1.0;
if(rndm==-1.0)
rndm=Initial();
else
rndm=3.99*rndm*(1.0-rndm);
//  0<（rndm*(1-radm)）<=0.25 可以用二次函数证明

return rndm;
}

int main()
{
/* 系统自带生成，产生[0,4]随机数*/
/*
srand((int)time(0));
for(int i=0;i<5;i++)
{
printf("%d  ",rand()%5);
}*/
int t=0;
while(t++<10)
printf("%lf     %lf\n", Random(), Initial());
getchar();
return 0;
}