%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