uinform
%产生一组[0,1]分布的随机数,采用余数法
%从文献中得知下列参数组较为适用:
%y的初值为1,w=7,模M=10^10(伪随机数周期为5*10^7)
%y的初值为1,w=5^13,模M=2^36(伪随机数周期为2^34,约2*10^10)
%y的初值为1,w=5^17,模M=2^42(伪随机数周期为2^40,约10^12)
%--------------------------------------------------------
%function x=uniform(y,w,M,n)
%n为要产生的随机数个数
riddling
%产生一个N(0,1)正态分布随机数
%采用筛选法,精度较高
%其他方法可参阅《现代应用数学手册--概率统计与随机过程卷》清华大学出版社马振华主编
%function y=riddling()
normal
%产生n个N(a,b)正态分布随机数
%其中a为均值,b为方差
%function x=normal(a,b,n)
function x=normal(a,b,n)