goldseq_gen.m

来自「运用群操作进行多用户检测」· M 代码 · 共 24 行

function goldseq = goldseq_gen(offset)
% gold sequences generator
p = offset;
r=6;N=2^r-1; 
s1(1:6)=[1 0 0 0 0 0]; %initial value 1
s2(1:6)=[1 0 0 0 0 0]; %initial value 1
f1=[1 0 0 0 0 1 1]; %特征多项式f＝x^6+x+1;
f2=[1 1 0 0 1 1 1]; %特征多项式f＝x^6+x^5+x^2+x+1;
for n=r+1:N %creat pn 1
s1(n)=mod(sum(s1(n-r:n-1).*f1(1:r)),2);
end
for n=r+1:N %creat pn 2
s2(n)=mod(sum(s2(n-r:n-1).*f2(1:r)),2);
end

if (p == 2^r)
    goldseq = s1;
elseif (p == 2^r+1)
    goldseq = s2;
else
    s=[s1(1+p:N) s1(1:p)]; %creat pn gold
    goldseq = mod(s+s2,2);    
end