mrand.m

来自「伪随机数发生器: “利用随机序列进行间隔控制」· M 代码 · 共 14 行

function [seq]=mrand(connection)
m=length(connection);
L=2^m-1;
registers=[zeros(1,m-1) 1];
seq(1)=registers(m);
for i=2:L
    new_reg_cont(1)=connection(1)*seq(i-1);
    for j=2:m
        new_reg_cont(j)=bitxor(registers(j-1),connection(j)*seq(i-1));
    end
    registers=new_reg_cont;
    seq(i)=registers(m);
end