线性同余法产生（0，1）均匀分布的随机数

原理

maltab代码

 clear all
close all
clc

a=16807;
c=0;
M=2^32;
X(1)=10^4;

N=4096;
for n=1:N-1
   X(n+1)=mod(a*X(n),M);
end
X=X/(M-1);
figure(1)
plot(X)

figure(2)%%
hist(X,32)

[counts,centers]=hist(X,32);
prob=counts/N;
figure(3)
plot(centers,prob)
ylim([0,1]);