mrand2.asv

伪随机数的产生程序 有m序列还有同余序列等等 用与嵌入信息位置的选择或密钥的制作

%%本函数用来生成m序列        不知道哪里错了
%输入格式为  seq=mrand（【1 0 0 1 0】）
%参数说明 connection为连接多项式，1表示相连  0表示断开
%seq为输出的结果m序列
function sep=mrand2(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