sjs.m

随机序列产生程序(用乘同余法产生随机数)

%用乘同余法产生随机数
A=6;  N=100;      %初始化；
x0=1;  M=255;
for k=1:N          %乘同余法递推100次；
   x2=A*x0;        %x2和x0分别表示xi和xi-1；
  x1=mod (x2,M);    %将x2存储器的数除以M，取余数放x1（xi）中；
   v1=x1/256;       %将x1存储器的数除以256得到小于1的随机数放v1中；
   v(:,k)=v1;        % 将v1中的数（ ）存放在矩阵存储器v的第k列中，v(:,k)
%表示行不变、列随递推循环次数变化；     
   x0=x1;           %xi-1= xi；
   v0=v1;
end                 %递推100次结束；
v2=v                %该语句末无‘；’，实现矩阵存储器v中随机数放在v2中，%且可直接显示在MATLAB的window中；
k1=k;
%grapher             %以下是绘图程序；
k=1:k1;
plot(k,v,k,v,'r');
xlabel('k'), ylabel('v');tktle('(0-1)均匀分布的随机序列')