가우시안 랜덤변수 발생기

double generateGaussian(float mean, float sigma) {
double v1, v2, s;

do {
v1 = 2 * ((double)rand() / RAND_MAX) - 1;      // -1.0 ~ 1.0 까지의 값
v2 = 2 * ((double)rand() / RAND_MAX) - 1;      // -1.0 ~ 1.0 까지의 값
s = v1 * v1 + v2 * v2;
} while (s >= 1 || s == 0);

s = sqrt((-2 * log(s)) / s);

double value = v1* s;
return value*sigma + mean ;
}

 

참조 : http://mwultong.blogspot.com/2006/11/c-vc-gaussian-random-number.html

C언어 VC++] 가우시안 랜덤 난수 발생; 가우스 분포 Gaussian Random Number