随机变量生成算法入门——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 InnoDB存储引擎中缓冲池Buffer Pool、Redo Log、Bin Log、Undo Log、Channge Buffer
查看>>
MySQL InnoDB引擎的锁机制详解
查看>>
Mysql INNODB引擎行锁的3种算法 Record Lock Next-Key Lock Grap Lock
查看>>
mysql InnoDB数据存储引擎 的B+树索引原理
查看>>
mysql innodb通过使用mvcc来实现可重复读
查看>>
mysql insert update 同时执行_MySQL进阶三板斧（三）看清“触发器 (Trigger)”的真实面目...
查看>>
mysql interval显示条件值_MySQL INTERVAL关键字可以使用哪些不同的单位值？
查看>>
Mysql join原理
查看>>
MySQL Join算法与调优白皮书（二）
查看>>
Mysql order by与limit混用陷阱
查看>>
Mysql order by与limit混用陷阱
查看>>
mysql order by多个字段排序
查看>>
MySQL Order By实现原理分析和Filesort优化
查看>>
mysql problems
查看>>
mysql replace first,MySQL中处理各种重复的一些方法
查看>>
MySQL replace函数替换字符串语句的用法（mysql字符串替换）
查看>>
mysql replace用法
查看>>
Mysql Row_Format 参数讲解
查看>>
mysql select， from ，join ，on ，where groupby,having ,order by limit的执行顺序和书写顺序
查看>>
MySQL Server 5.5安装记录
查看>>