连续创建 0 到 1 之间的随机数

Question

我正在尝试用 C++ 对股票价格变动进行建模。我需要创建一个介于 0 到 1 之间的随机数。

但似乎随机数生成器的值不断增加，并不是真正随机的。

代码如下所示：

#include<iostream>
#include<ctime>
#include<cstdlib>
#include<time.h>

using namespace std;

int main()
{
    double stockPrice = 25;
    int start = 0, end = 0;
    start = clock();

    srand (time(NULL));
    cout << (double) rand() / (double) (RAND_MAX) << endl;
    system("pause");

    while(stockPrice > 18)
    {
        if(stockPrice == 20)
        {
            double probability = (rand()/(double)RAND_MAX);
            if(probability <= (1/10))
            {
                stockPrice = stockPrice-1;
            }
            else
            {
                stockPrice = stockPrice +1;
            }
        }
        else if (stockPrice < 20)
        {
            double probability = (rand()/(double)RAND_MAX);
            if(probability <= (1/3))
            {
                stockPrice = stockPrice -1;
            }
            else
            {
                stockPrice = stockPrice +1;
            }
        }
        else 
        {
            double probability = (rand()/(double)RAND_MAX);
            if(probability <= (2/3))
            {
                stockPrice = stockPrice -1;
            }
            else
            {
                stockPrice = stockPrice +1;
            }
        }
        cout << stockPrice << endl;
    }

    end = clock();

    double t = (double)(start-end)/CLOCKS_PER_SEC;
    cout << t << endl;
    system("pause");
}

不知道如何解决这个问题..需要一些指导...

Answer 1

需要一些指导...

指导 1：

正确的比较，你应该使用

double probability = (rand()/(double)(RAND_MAX + 1));
                                               ^
                                        for better scaling

因为当前在 if(probability <= (1/10)) 行中，您正在与 0 进行比较，因为将 1/10 1/3 和 2/3 转换为整数

指南 2：

毕竟你可能会使用更好的生成器statistical properties

#include <random>
#include <iostream>

    std::random_device rd;
    std::mt19937 gen(rd());
    std::uniform_real_distribution<> dis(0, 1);
    double uniformOn01 = dis(gen);

Answer 2

if(probability <= (1/3))
    {
    stockPrice = stockPrice -1;
    }
else
    {
    stockPrice = stockPrice +1;
    }

由于probabililty 永远不会是负数，此代码将几乎总是增加stockPrice 的值。唯一不会的情况是probability 为0。那是因为1/3 是整数 除法，其值为0。将所有这些分数更改为1.0/3 之类的东西会好很多。这与随机数生成器的质量无关。有些人在看到rand 时会感到非常紧张，以至于他们什么也没看到。

但是，代码中的缩放存在一个缺陷。而不是

double probability = (rand()/(double)RAND_MAX);

使用

double probability = (rand()/(double)(RAND_MAX + 1));

正如最初所写，每当rand() 产生值RAND_MAX 时，probability 的值将是 1，并且它会更频繁地产生其他值。

Answer 3

在您比较整数除法概率的所有情况下，所有结果都为 0... 而不是您想要的 (2/3) (2/(double)3) 或 2/3.0。