m_sequence.m

这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的

function mseq = m_sequence(fbconnection)
n = length(fbconnection);
N = 2^n-1;
register = [zeros(1,n-1) 1];
mseq(1)=register(n);
for i = 2:N
    newregister(1)=mod(sum(fbconnection.*register),2);
    for j = 2:n
        newregister(j) = register(j-1);
    end
    register = newregister;
    mseq(i) = register(n);
end