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

Wichmann-Hill算法

1、算法介绍

为了得到长周期的随机数发生器，可将两个周期相差很小的波形序列相加，

如定义3个随机数发生器，它们都是全周期的

2、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;}
   

3、使用matlab画直方图

u=[…(从C++程序生成文件中导入数据)]histogram(u)

4、生成直方图预览