function y=junyun(n) 
%0-1的均勻分布，n代表數據量，
%
y=ones(1,n); %設單位序列x
x=ones(1,n); %設單位序列y
m=100000;    %取模值
x0=mod(ceil(m*rand(1,1)),m); %對隨機數進行取模
x0=floor(x0/2); %除2運算
x0=2*x0+1;    
u=21;%混合同余法 長周期相關，高度不均勻性
x(1)=x0; 
for i=1:n-1 
x(i+1)=u*x(i)+3; 
x(i+1)=mod(x(i+1),m); 
x(i)=x(i)/m; %x(n)單位化 
end 
y=x;