www.pudn.com > Arnold.rar > wei.m, change:2014-07-12,size:505b


function b=wei() 
% x=input('请输入2个大素数p和q,且两者同余:'); 
% s=input('请输入一个随机种子数S,要求S与n互素(n=p*q):'); 
% m=input('请输入随机数的位数(二进制):'); 
x=[383,503]; 
s=101355; 
m=256*256; 
c=zeros(1,m); 
b=c; 
n=x(1)*x(2); 
x=mod(s*s,n); 
c(1)=mod(x^2,n); 
b(1)=mod(c(1),2); 
for i=2:m 
    c(i)=mod(c(i-1)^2,n); 
    b(i)=mod(c(i),2); 
end 
% disp('输出伪随机数:'); 
% disp(b); 
% disp('输出他的十进制数:'); 
% a=0; 
% for i=1:m 
%     a=a+b(i)*(2^(m-i)); 
% end 
% disp(a);