任意阶平衡gold码的产生及其互相关性(代码)

clear all;
n=input('请输入移位寄存器的级数n,n是大于等于2的整数，且n不能为4的倍数，为4的倍数不保证结果正确: ');%五级移位%
N=2^n-1;%周期为N%
x=gfprimfd(n,'all');%产生级数为n的本原多项式%
[row,column]=size(x);
if(n==2)%n为2是不能用下面求y的公式，作为特例直接写入%
y=[1,1];
end
if(n>2)
y=x(1:row,2:column);%本来是nXn矩阵，第一列都是1，不要了，现在是nX(n-1)矩阵% end
m=zeros(row,N);%求出所有m序列%
for r=1:row%从1到第row行%
register7=[zeros(1,n-1),1];%寄存器7，表示移位寄存器的初态除了最后一位，其余全0%
m(r,1)=register7(n);%m（r，1）是的第r行第一列的取值%
for i=2:1:N%从2到N，连续取值%
register8(1)=mod(sum(y(r,:).*register7),2);%移位寄存器的向量与反馈系数的向量点乘后，新向量的所有数模二加，得到一个数作为寄存器8的第一位%
for j=2:1:n%从2到N，连续取值%
register8(j)=register7(j-1);%确定后n-1位%
end
register7=register8;%保存移位寄存器现在的值到register1%
m(r,i)=register7(n);%输出的书，从第2位开始%
end %现在得到了第r行的m序列%
end %现在得到了第1到第r行的m序列%
if(n>2)%列出m1最大行数，超过会报错%
fprintf('m序列的行数小于等于%',row)%显示之后要输入m1的要求%
end
if(n==2)%列出m1最大行数，超过会报错%
fprintf('m序列的行数小于等于1')%显示之后要输入m1的要求%
end
fprintf('\n')%换行%
pass=1;%判断m1序列是否有优选对，若有，pass在下面循环中置0%
while(pass==1) %若m1没有优选对，则循环%
nom1=input('请从矩阵y中选出你所要m1序列本愿多项式的行数: ');%输入一个数赋给nom1，number of m1,即m1的行数%
m1=m(nom1,:);%得到了你指定的m1序列%
rm=zeros(1,2^n-1);%两个m序列的互相关性，用来判断m优选对%
s=zeros(1,2^n-1);%m优选对中移位后相同数的个数%
d=zeros(1,2^n-1);%m优选对中移位后相同数的个数%
register2=zeros(1,2^n-1);%寄存器2，用以验证m序列优选对时移位%
cmatch=0;%counter match优选对匹配的数量%
match=zeros(1,N);%匹配优选对所在的行数%
matchs=zeros(1,N);%match sequence匹配优选对的序列%
for k=1:row%对m序列每一行循环%
m2=m(k,:);%把m序列的第k行赋给m2%
register1=m2;%寄存器1，用以验证m序列优选对时移位%
for i=1:2^n-1
A=mod((m1+register1),2);%移位加N次%
for j=1:2^n-1%逐位判断m1与m2是否相同，判断N次%
if A(j)==0%若第j位位相同%
s(i)=s(i)+1;%相同数加1，之后清零%
end
d(i)=N-s(i);%不同数减一，因为是s（i）的因变量，不用清0%
end
rm(i)=s(i)-d(i);%rm是互相关函数，因为是s（i）和d（i）的因变量，不用清0%
s(i)=0;%必须清0，不然k的循环序列就会叠加%
register2(1)=register1(N);%以下三行为循环右移%
register1(2:N)=register1(1:N-1);
register1(1)=register2(1);
end
maxrm=max(abs(rm));%maxrm为互相关值的最大绝对值，因为是rm的因变量，所以不用清0%
if(mod(n,2)==1)%对2除n求余，判断n的奇偶性%
if(maxrm==2^((n+1)/2)+1)%奇数判断互相关最大绝对值是否满足书上条件4—37%
cmatch=cmatch+1;%匹配（即互为优选对）个数加1%
match(1,cmatch)=k;%匹配行数写入match序列中%
matchs(cmatch,:)=m2;%匹配序列写入matchs矩阵中%
end
end
if(mod(n,2)==0)%同上%
if(maxrm==2^((n+2)/2)+1)%偶数判断互相关最大绝对值是否满足书上条件4—37%
cmatch=cmatch+1;
match(1,cmatch)=k;
matchs(cmatch,:)=m2;
end
end
end
if(cmatch~=0)%若m1有优选对，cmatch不等于0%
pass=0;%pass置0，while不在循环%
end
if(cmatch==0)%若m1没有优选对，cmatch等于0%
fprintf('m1没有优选对，请重新选择m1行数')%输出括号中的数%
fprintf('\n')%换行%
end
end
fprintf('下列数字表示与m1互为优选对的m序列的行数')%输出括号中的数%
match(1:cmatch)%因为后面没加分号，执行程序时会显示match序列的第1到第camtch列% choice=input('是否依次显示m1的优选对y/n：','s');%按括号中文字内容从键盘输入一个置%
if(choice=='y')%输入y执行，输入其他任何数不执行%
matchs%matchs，match sequence,表示m1的所有优选对%
end
nom2=input('请从矩阵y中选出你所要m2序列的行数: ');%同中文%
m2=m(nom2,:);%得到m2，用优选对m1,m2，生成下面的gold序列，并选出平衡gold序列% gold=zeros(N+2,N);%一共N+2个gold序列，每个周期为N，设（N+2）*2的矩阵% register3=m2;%给寄存器3赋值m2，下面循环右移要用%
register4=zeros(1,N);%给寄存器4设置成1xN的矩阵，下面循环右移要用%
for i=1:N%循环N次%
gold(i,:)=mod((m1+register3),2);%两m序列逐位模二加%
register4(1)=register3(N);%下面三行为循环右移%
register3(2:N)=register3(1:N-1);
register3(1)=register4(1);
end;
gold(N+1,:)=m1;%m序列包含生成它的2个m序列%
gold(N+2,:)=m2;
c1=zeros(1,N+2);%代表1的个数counter 1%
c0=zeros(1,N+2);%代表0的个数counter 2%
cbg=0;%代表平衡gold序列的个数counter balence gold%
bg=zeros(1,N);%代表平衡gold序列%
for i=1:N+2%开始找，要找N+2行%
for j=1:N%弄出来一个全1的向量，长度为N，命名为judge函数
judge(1,j)=0;
end
s2=mod((gold(i,:)+judge),2);%s2是一个gold序列与全1序列模二加得到的序列% for k=1:N %看有多少个1%
if s2(k)==1
c1(i)=c1(i)+1;
end
end
for k=1:N %有多少个0%
if s2(k)==0
c0(i)=c0(i)+1;
end
end
if c1(i)-c0(i)==1%1的个数比0多1，找到一个平衡序列%
cbg=cbg+1;%平衡gold序列个数加1，可以看是不是约50%或75%是平衡的% bg(cbg,:)=gold(i,:);%这个个平衡序列复制到bg中%
end
end
fprintf('已生成%个平衡gold序列',cbg)
fprintf('\n')。