博客
关于我
随机变量生成算法入门——Wichmann-Hill算法
阅读量:375 次
发布时间:2019-03-05

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

Wichmann-Hill算法

1、算法介绍

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

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

在这里插入图片描述
则它们合成输出
在这里插入图片描述
等价的乘性发生器中a=16555425264690,m=2.7817*1013。

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、生成直方图预览

(0,1)均匀分布

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

你可能感兴趣的文章
MySQL(1)的使用 | SQL
查看>>
MySQL(2)DDL详解
查看>>
MySQL(3)DML详解
查看>>
MySQL(4)运算符 | 关联查询详解
查看>>
MySQL(5)条件查询 | 单行函数 | 事务详解
查看>>
Mysql,group by分组查询、order by排序查询、join连接查询、union联合查询
查看>>
Mysql,sql文件导入和导出
查看>>
MYSQL:int类型升级到bigint,对PHP开发语言影响
查看>>
Mysql:mysql 5.X 报错 ERROR 1193 (HY000): Unknown system variable ‘validate_password_length‘
查看>>
MySQL:MySQL执行一条SQL查询语句的执行过程
查看>>
Mysql:SQL性能分析
查看>>
mysql:SQL按时间查询方法总结
查看>>
MySQL:什么样的字段适合加索引?什么样的字段不适合加索引
查看>>
MySQL:判断逗号分隔的字符串中是否包含某个字符串
查看>>
MySQL:某个ip连接mysql失败次数过多,导致ip锁定
查看>>
MySQL:索引失效场景总结
查看>>
Mysql:避免重复的插入数据方法汇总
查看>>
MyS中的IF
查看>>
M_Map工具箱简介及地理图形绘制
查看>>
m_Orchestrate learning system---二十二、html代码如何变的容易
查看>>