建模和仿真作业-随机数生成.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 座机电话号码9; a2 power 5,15 ; c1 座机电话号码5; 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 a2*x i-1 +c2; x i mod y, M2 ; r i-1 x i /M2; end end format short; 对（0，