tl rass 1 - polynome.m

来自「3种伪随机数列发生器的matlab编码以及其各种特性的分析」· M 代码 · 共 13 行

function res = polynome(n)
res = zer(n,1);
masque = [ 0 0 0 1 0 0 0 0 0 0 0 0 1 0 1 1 ];
masque = masque';
memoire = [ 0 0 0 1 0 0 0 1 0 0 0 0 1 0 1 1 ];
memoire = memoire';
for k = 1:n
       res(k)=rem(double(masque&memoire)'*double(masque&memoire),2);
       memoire(2:16)=memoire'(1:15);
       memoire(1)=res(k);
end