用乘同余法产生随机数7页.doc

系统辨识与MATLAB仿真附录例2.1 用乘同余法产生随机数 (见光盘 FLch2sjxleg1.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中；k