随机变量生成算法入门——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 查询，正数降序排序，负数升序排序
查看>>
MySQL 树形结构 根据指定节点 获取其下属的所有子节点（包含路径上的枝干节点和叶子节点）...
查看>>
mysql 死锁 Deadlock found when trying to get lock； try restarting transaction
查看>>
mysql 死锁(先delete 后insert)日志分析
查看>>
MySQL 死锁了，怎么办？
查看>>
MySQL 深度分页性能急剧下降，该如何优化？
查看>>
MySQL 深度分页性能急剧下降，该如何优化？
查看>>
MySQL 添加列，修改列，删除列
查看>>
mysql 添加索引
查看>>
MySQL 添加索引，删除索引及其用法
查看>>
MySQL 用 limit 为什么会影响性能？
查看>>
MySQL 用 limit 为什么会影响性能？有什么优化方案？
查看>>
MySQL 用户权限管理：授权、撤销、密码更新和用户删除（图文解析）
查看>>
mysql 用户管理和权限设置
查看>>
MySQL 的 varchar 水真的太深了！
查看>>
mysql 的GROUP_CONCAT函数的使用（group_by 如何显示分组之前的数据）
查看>>
MySQL 的instr函数
查看>>
MySQL 的mysql_secure_installation安全脚本执行过程介绍
查看>>
MySQL 的Rename Table语句
查看>>
MySQL 的全局锁、表锁和行锁
查看>>