基于matlab的gold码的生成及相关性分析

clc;N=7;%以7级寄存器为例，并组其中的一组优选对：211，,217connections=gfprimfd(N,'all');f1=connections(4,:);%取一组本原多项式序列，211f2=connections(16,:);%取另一组本原多项式序列,217registers1=[0000001];%给定寄存器的初始状态registers2=[0000001];%取相同的初始状态L=2^N-1;%周期长度sum2=0;sum1=0;fork=1:Lseq1(k)=registers1(N);%第一组m序列seq2(k)=registers2(N);%第二组序列forj=1:N%进行模2加sum1=sum1+f1(j+1)*registers1(j);%各级寄存器送参与模2加的值sum1=mod(sum1,2);sum2=sum2+f2(j+1)*registers2(j);%各级寄存器送参与模2加的值sum2=mod(sum2,2);endfort=N:-1:2%寄存器移位registers1(t)=registers1(t-1);registers2(t)=registers2(t-1);endregisters1(1)=sum1;registers2(1)=sum2;sum2=0;sum1=0;enddisp(f1);disp(f2);z=seq1+seq2;%m序列的相加gold=mod(z,2);%模2运算gold=1-2*gold;%转换为2值电平disp(gold);R=xcorr(gold,'unbiased');%自相关R=R/max(R);%归一化figure;plot(R);title('gold序列的自相关函数');s=fftshift(abs(fft(gold,2*L)).^2);%求功率谱s=s/max(s);figure;plot(s);title('gold序列的功率谱');