随机变量生成算法入门——Wichmann-Hill算法
本文共 682 字，大约阅读时间需要 2 分钟。
Wichmann-Hill算法
1、算法介绍
为了得到长周期的随机数发生器，可将两个周期相差很小的波形序列相加，
如定义3个随机数发生器，它们都是全周期的
则它们合成输出 
等价的乘性发生器中a=16555425264690，m=2.7817*1013。 2、C++代码实现
includeint 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、生成直方图预览
转载地址：http://wpdg.baihongyu.com/
你可能感兴趣的文章
mysql 字段合并问题（group_concat）
查看>>
mysql 字段类型类型
查看>>
MySQL 字符串截取函数,字段截取，字符串截取
查看>>
MySQL 存储引擎
查看>>
mysql 存储过程 注入_mysql 视图 事务 存储过程 SQL注入
查看>>
MySQL 存储过程参数：in、out、inout
查看>>
mysql 存储过程每隔一段时间执行一次
查看>>
mysql 存在update不存在insert
查看>>
Mysql 学习总结（86）—— Mysql 的 JSON 数据类型正确使用姿势
查看>>
Mysql 学习总结（87）—— Mysql 执行计划（Explain）再总结
查看>>
Mysql 学习总结（88）—— Mysql 官方为什么不推荐用雪花 id 和 uuid 做 MySQL 主键
查看>>
Mysql 学习总结（89）—— Mysql 库表容量统计
查看>>
mysql 实现主从复制/主从同步
查看>>
mysql 审核_审核MySQL数据库上的登录
查看>>
mysql 导入 sql 文件时 ERROR 1046 (3D000) no database selected 错误的解决
查看>>
mysql 导入导出大文件
查看>>
mysql 将null转代为0
查看>>
mysql 常用
查看>>
MySQL 常用列类型
查看>>
mysql 常用命令
查看>>