产生0到1之间均匀分布的一个随机数与随机数序列

1、功能：

产生0到1之间均匀分布的一个随机数
2、方法说明：

设m=216，产生产生0到1之间均匀分布的一个随机数的计算公式如下：

ri=mod(2053ri-1 +13849,m),i=1,2,3……

pi =ri/m

例:连续产生10个0到1之间均匀分布的一个随机数 。r的初值取5.0

1 #include "stdio.h"
2 double rnd1(double *r)
3 { int m;
4 double s,u,v,p;
5 s=65536.0; u=2053.0; v=13849.0;
6 m=(int)(*r/s); *r=*r-m*s;
7 *r=u*(*r)+v; m=(int)(*r/s);
8 *r=*r-m*s; p=*r/s;
9 return(p);
}
main()
{ int i;
double r;
r=5.0;
printf("\n");
for (i=0; i<=9; i++)
printf("%10.7lf\n",rnd1(&r));
printf("\n");
getch();
}

3、功能：
产生0到1之间均匀分布的一个随机数序列

例：产生50个0到1之间均匀分布的一个随机数序列，r的初值取1.0

1
2
3 #include "stdio.h"
4 void rnds(r,p,n)
5 double *r,p[];
6 int n;
7 { int i,m;
8 double s,u,v;
9 s=65536.0; u=2053.0; v=13849.0;
for (i=0; i<=n-1; i++)
{ *r=u*(*r)+v; m=(int)(*r/s);
*r=*r-m*s; p[i]=*r/s;
}
return;
}
main()
{ int i,j;
double p[50],r;
r=1.0;
rnds(&r,p,50);
printf("\n");
for (i=0; i<=9; i++)
{ for (j=0; j<=4; j++)
printf("%10.7lf ",p[5*i+j]);
printf("\n");
}
printf("\n");
getch();
}