gold.m

Gold sequence generator

echo on
connections1=[1 0 1 0 0];
connections2=[1 1 1 0 1];
sequence1=s_mlsrs(connections1);
sequence2=s_mlsrs(connections2);
L=2^length(connections1)-1;
for shift_amount=0:L-1,
   temp=[sequence2(shift_amount+1:L) sequence2(1:shift_amount)]; %循环左移
   gold_seq(shift_amount+1,:)=(sequence1+temp)-floor((sequence1+temp)./2).*2; %为了使1+1=0，因为sequence1+temp是十进制运算（1＋1＝2，1＋0＝1）
end;
max_cross_corr=7;
for i=1:L-1,
    for j=i+1:L,
       c1=2*gold_seq(i,:)-1;
       c2=2*gold_seq(j,:)-1;
       for m=0:L-1,
           shifted_c2=[c2(m+1:L) c2(1:m)];
           corr=abs(sum(c1.*shifted_c2));
           if (corr<max_cross_corr),
               max_cross_corr=corr;
               if (max_cross_corr==1),
                   gold1=c1;
                   gold2=shifted_c2;
                   p=i;q=j;
               end;
           end;    
       end;
    end;
end