有谁能搞明白这小段代码吗?
//摘自徐士良的算法一书,讲的是高斯白噪声的产生,但是看不懂...,过来请教下高手
double s0=65536.0; //
double w0=2053.0; //
double v0=13849.0; //这3个常量数字不知道怎么来的,郁闷
double r0=0;
double t0=0.0;
int m0;
double noise()
{
t0=0.0;
//for循环里面是什么意思呢?也看不懂什么意思
//不知道作者是怎么得来的
for(int i=0;i <12;i++)
{
r0=w0*r0+v0;
m0=r0/s0;
r0=r0-m0*s0;
t0=t0+r0/s0;
}
return(t0-6.0); //返回白噪声数据
}
[解决办法]
线性同余.
具体参看 <c语言数值酸法程序大全> 2th
(Numerical recipes in c,the art of Scientific Computing)