%%乘同余法计算随机数列
%%公式:Sk+1=A*Sk*mod(M)
M=64;%
A=11;
x0=1;
Sk=x0
%%输出n位随机数列
n=10;
for i=1:n
Sk=mod(A*Sk,M)
end
%%混合同余法计算随机数列
%%迭代公式:Sk+1=(A*Sk+C)mod(M),当C为0时,混合同余法退化为乘同余法
clc;clear
M=64;
A=13;
C=7;
x0=1;
Sk=x0
%%输出n位随机数列
n=10;
for i=1:n
Sk=mod(A*Sk+C,M)
end