de_sequence.m

来自「一种新的光纤通信调制技术」· M 代码 · 共 18 行

function [a]=de_sequence(fbconnection);
n=length(fbconnection);
N=2^n-1;
register=[1 zeros(1,n-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%获得PN序列
a=zeros(1,length(mseq)+1);
a(1)=0;
for i=1:length(mseq)
    a(i+1)=mseq(i);
end%生成debruijin序列