我创建了一个使用 rand() 在 C 中生成重复数字的程序。

然而,重复的数字不遵循 Central Limit Theorem

任何人都可以解决这个 rand() 错误问题,或者除了使用 rand() C 库来生成更好的随机数之外,还有其他选择吗?

这是屏幕截图:

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <time.h>
#include <windows.h>


#define TOTAL_THROW 10000000

typedef enum _COINTOSS {
    UNDEFINED = 0,
    HEAD = 1,
    TAIL = 2
} COINTOSS;

COINTOSS toss_coin () {
    int x = rand() % 2;
    if (x == 0) return HEAD;
    else if (x == 1) return TAIL;
}

void main () {
    int x, i, j, v1 = 0, v2 = 200, total = 0;
    int head_range[25] = {0};
    int tail_range[25] = {0};
    int no_range = 0;
    int count = 0;
    int repeated = 0;
    COINTOSS previos_toss = UNDEFINED;
    COINTOSS current_toss;

    srand(time(NULL));

    for (i=0; i<TOTAL_THROW; i++) {
        current_toss = toss_coin();             // current toss
        if (previos_toss == current_toss) {
            count++;
        } else {
            if (current_toss == HEAD) head_range[count] += 1;
            else if (current_toss == TAIL) tail_range[count] += 1;


            previos_toss = current_toss;
            count = 0;
        }

    }

    for (i=24; i>=0; i--) {
        printf("+%d = %d\n", i+1, head_range[i]);
    }

    puts("________________\n");

    for (i=0; i<25; i++) {
        printf("-%d = %d\n", i+1, tail_range[i]);
    }

    printf("\nTOTAL_THROW: %d\n", TOTAL_THROW);


    printf("\nPress [ENTER] to exit. . .");
    getchar();
}

最佳答案

您的问题是使用模数使您的随机数进入所需的范围,它使用较低的位(这是一个经典的问题):

int x = rand() % 2;
rand() 的低位(a linear congruential generator (LCG))不像高位那样随机。这适用于所有 LCG,无论库或语言如何。

对于 [0..N) 的范围,您应该执行以下操作(使用高位):
int r = rand() / ( RAND_MAX / N + 1 );

关于c - rand() 不遵循高斯分布和中心极限定理,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/26008137/

10-13 05:38