本文共 868 字，大约阅读时间需要 2 分钟。

Wichmann-Hill算法

算法介绍

Wichmann-Hill算法是一种生成长周期随机数的方法。其基本思想是将两个周期相差很小的波形序列叠加，通过这种方式可以生成一个更长周期的随机数发生器。具体来说，定义三个随机数发生器，它们都是全周期的。然后将它们合成输出，形成一个等价的乘性发生器。

C++代码实现

#include 
   
    int main() {    float x[10001], y[10001], z[10001], u[10001];    x[0] = 1;    y[0] = 2;    z[0] = 3;    // 种子可取任意整数    for (int i = 0; i < 10000; ++i) {        x[i+1] = fmod(171 * x[i], 30269);        y[i+1] = fmod(170 * y[i], 30307);        z[i+1] = fmod(172 * z[i], 30323);        u[i+1] = fmod(x[i+1] / 30269 + y[i+1] / 30307 + z[i+1] / 30323, 1);    }    FILE *fp;    fp = fopen("/*输入自定义路径*/", "w");    for (int i = 1; i <= 10000; ++i) {        fprintf(fp, "%f ", u[i]);    }    fclose(fp);    return 0;}
   

使用MATLAB画直方图

将C++生成的文件导入MATLAB后，使用histogram函数绘制直方图即可。以下是具体步骤：

通过这些步骤，可以直观地查看随机数分布的情况。

生成直方图预览

此外，以下直方图展示了生成的随机数在[0,1)区间内的均匀分布特性。

转载地址：http://wpdg.baihongyu.com/