建模与仿真作业随机数生成.docVIP

建模与仿真作业 随机数生成 本次作业使用MATLAB编写程序。 线性同余法 基本公式： 其中a为乘子，为种子，c为常数，M为模 迭代步骤： （1）设定a, c, M值，并给定初始种子 （2）令 （3） （4） （5），转入（3） 本题采用著名的 Coveyou与 Macpherson 混合同余发生器 及Kobayashi 混合同余发生器 程序代码： 线性同余法函数： function r = MixMOD(x0,n,type) %随机数种子：x0 %产生的随机数个数：n %采用的混合同余法的公式类型：type %产生的随机数序列：r format long; M1 = power(2,31); M2 = power(2,35); a1 = 314159269; a2 = power(5,15); c1 = 453806245; c2 = 1; r = zeros(n,1); x = zeros(n+1,1); x(1) = x0; if type == 1 for i=2:n+1 y = a1*x(i-1)+c1; x(i) = mod(y, M1); r(i-1) = x(i)/M1; end else for i=2:n+1 y =