marand.asv

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

%m序列的产生   又名为：  N级线性最长反馈序列    这个是网上找的  但是结果很奇怪  最好一个如果为0时  输出的数很奇怪
%输入格式seq=mrand([1 0 0 1 0])  1代表相连 0代表断开
%主函数 
function m_rand=marand(fbconnection)
%fbconnection=[0 1 0 0 1];  %特征多项式
m=1;                       %输出周期个数
t=length(fbconnection);    %序列长度
chu=[zeros(1,t-1),1];      %初态
m_out=m_sequence(fbconnection,m,chu);
m_rand=m_out;

function [m_out]=m_sequence(fbconnection,m,chu)
len=length(fbconnection);
n=2^len-1;                %序列长度
register=chu;             %初态
fan=0;
for i=1:n*m               % i循环一次输出一个m_out
    m_out(i)=register(len);
    for k=1:len                  % k循环累加反馈
        fan=fan+fbconnection(k)*register(k);
    end
    fan=mod(fan,2);            % 反馈系数
    newregister(1)=fan;
    for j=2:len               % j循环移位
        newregister(j)=register(j-1);
    end
    register=newregister;
    fan=0; 
end