随机变量生成算法入门——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 slave 停了_slave 停止。求解决方法
查看>>
MySQL SQL 优化指南：主键、ORDER BY、GROUP BY 和 UPDATE 优化详解
查看>>
MYSQL sql语句针对数据记录时间范围查询的效率对比
查看>>
mysql sum 没返回,如果没有找到任何值，我如何在MySQL中获得SUM函数以返回'0'？
查看>>
mysql Timestamp时间隔了8小时
查看>>
Mysql tinyint(1)与tinyint(4)的区别
查看>>
mysql union orderby 无效
查看>>
mysql v$session_Oracle 进程查看v$session
查看>>
mysql where中如何判断不为空
查看>>
MySQL Workbench 使用手册：从入门到精通
查看>>
mysql workbench6.3.5_MySQL Workbench
查看>>
MySQL Workbench安装教程以及菜单汉化
查看>>
MySQL Xtrabackup 安装、备份、恢复
查看>>
mysql [Err] 1436 - Thread stack overrun: 129464 bytes used of a 286720 byte stack, and 160000 bytes
查看>>
MySQL _ MySQL常用操作
查看>>
MySQL – 导出数据成csv
查看>>
MySQL —— 在CentOS9下安装MySQL
查看>>
MySQL —— 视图
查看>>
mysql 不区分大小写
查看>>
mysql 两列互转
查看>>