在矩形上均匀生成点（x，y）

Question

我想在矩形上均匀地生成点(x,y)。首先你输入 minX 和 maxX 和 minY maxY 然后你统一生成 (x,y)，下面显示的基本代码有没有更好的方法来实现它？ （我需要它以蒙特卡罗方法来制作情节）

#include <iostream>
#include <random>
double drand(const double min = 0., const double max = 1.)
{
    return (max - min) * static_cast<double>(rand()) / static_cast<double>      (RAND_MAX)+min;
}
int main(int argc, char **argv)
{
    for(unsigned short int i=0;i<1000;++i)
    {
        std::cout << "x " << drand(minX, maxX) << std::endl;
        std::cout << "y " << drand(0., maxY) << std::endl;
    }
return 0;
}

Answer 1

您正在使用rand()。因为几乎任何东西都比rand() 好，所以你的障碍并没有那么高。既然您显然拥有 C++11，为什么不使用 C++11 提供的更高级的随机数生成器呢？例如，

#include <random>

double drand(const double min = 0., const double max = 1.)
{
   static std:: mt19937_64 mt_generator(std::random_device());
   static std::uniform_real_distribution<double> u_0_1(0.0, 1.0);

   return min + (max - min) * u_0_1(mt_generator);
}

Answer 2

除非我误读了您的问题，否则我希望看到类似这样的 Java 代码：

int nx = 10;
int ny = 10;
double x = xmin;
double dx = (xmax-xmin)/nx;
double dy = (ymax-ymin)/ny;
int id = 0;
for (int i = 0; i < nx; ++i) {
    double y = ymin;
    for (int j = 0; j < ny; ++j) {
        System.out.println(id, x, y);
        id++;
        y += dy;
    }
    x += dx;
}

我打印了这些值，但您可能希望将它们存储在一个数据结构中，您可以将其提供给您的 Monte Carlo 模拟进行评估。