高斯随机数生成算法问题这是网上盛行的一个生成高斯随机数的函数,但它并没有定义均值和方差啊。难道这个程
高斯随机数生成算法问题
这是网上盛行的一个生成高斯随机数的函数,但它并没有定义均值和方差啊。难道这个程序默认为标准正态分布N(0,1)?
#include <stdlib.h>
#include <math.h>
double gaussrand()
{
static double V1, V2, S;
static int phase = 0;
double X;
if(phase == 0) {
do {
double U1 = (double)rand() / RAND_MAX;
double U2 = (double)rand() / RAND_MAX;
V1 = 2 * U1 - 1;
V2 = 2 * U2 - 1;
S = V1 * V1 + V2 * V2;
} while(S >= 1 || S == 0);
X = V1 * sqrt(-2 * log(S) / S);
} else
X = V2 * sqrt(-2 * log(S) / S);
phase = 1 - phase;
return X;
}
[解决办法]
楼主可以根据X = V1 * sqrt(-2 * log(S) / S);得出啊,
要想得到其他正态分布数值,自己传入参数,更改公式
