Re: 求高斯白噪声的C程序

sigma = 方差;
x = 均值
pointU=(double)rand()/(double)(RAND_MAX + eps); //eps为一很小的数，保证计算p
ointR时不会出现除零错误
pointv=(double)rand()/(double)(RAND_MAX + eps);
pointR = sigma * sqrt(2*log(1.0/((1.0 - pointU))));
pointGauss = pointR * cos(2* PI * pointv) ;
return  x + pointGauss;

.4

那样会有log(0)错误吧，不过改改上面两行程序也可以
具体原因我自己也给忘了
-__________________-!!

嗯，要是为了避免 log0 ，倒是的确可以用 1 来减一下...
——不过为啥要取倒数，这个比较费解... hehe..
.4

我自己现在看着也挺费解...忘了当时怎么想的了