C언어에서 가우시안 난수 발생시 srand 사용해도 될까요 ?
안녕하세요. 프로그래밍초보 입니다.
대부분 Python기반으로 코딩하다가 속도 Issue때문에 C언어로 코딩중에 있습니다.
<코드>
double gaussianRandom(double AVR, double SD)
{
double V1, V2, S;
int phase = 0;
double X;
if (phase == 0) {
do {
V1 = 2* ((double)rand() / RAND_MAX) -1;
V2 = 2* ((double)rand() / RAND_MAX) -1;
S = V1 * V1 + V2 * V2;
} while (S >= 1 || S == 0);
X = V1 * sqrt(-2* log(S)/S);
}
X = (SD*X) + AVR;
phase = 1 - phase;
return X;
}
코드는 구글링을 통해 얻은 것인데 어떤 원리로 동작하는지를 사실 잘 모르겠습니다..ㅠㅜ
제가 궁금한 것은
rand() 대신에 srand(time(NULL)) 로 수정하여도 가우시안 분포로 난수가 발생할까요 ??
<상황>
for문으로 가우시안 난수발생 후 어떤 연산 A를 거쳐서 B를 뽑는 코드입니다.
가우시안 난수 발생을 1*10^9 / 2*10^9 개 했을 때는 결과가 도출됐는데
3*10^9개 부터는 결과가 도출이 안되더라구요.(아마 2^31을 넘어서 그런듯?)
1*10^9개의 프로그램을 여러개 돌리는 식으로 그 이상의 갯수를 추출하려고 합니다.
rand()함수를 사용하면 프로그램 돌려도 결과가 같을 거라는 것때문..
...
많은 분들이 헷갈리시는데 srand는 난수를 생성하는 함수가 아니라 *난수 생성기를 초기화하는* 함수입니다. 부르려면 난수를 사용하기 전에 처음에 딱 한번 불러야 합니다. "rand 대신에" srand를 부른다는 건 말이 안 됩니다. (일단 srand는 결과값도 없는데요..)
srand를 여러번 부르면 그만큼 훌륭한 난수가 나오는 게 아니라 부를 때마다 다시 초기화가 되어서 아까 나왔던 값이 또 나옵니다. 즉 더 이상 난수라고 할 수 없습니다.
아 최초 1번만 설정하면 되는 것이엿군요.
아 최초 1번만 설정하면 되는 것이엿군요. 이해했습니다. 감사합니다!
댓글 달기