main_corr.m

输入为一个序列或两个序列

%计算m序列和gold序列自相关互相关的主程序
mseq1=msequence(211);%生成序列
mseq2=msequence(217);
goldseq=goldsequence(211,217);
mark1=find(mseq1==0);mseq1(mark1)=-1;%将序列值改为1和-1
mark2=find(mseq2==0);mseq2(mark2)=-1;
mark3=find(goldseq==0);goldseq(mark3)=-1;
r1=correlate(mseq1);%计算相关值
r2=correlate(mseq1,mseq2);
r3=correlate(goldseq(1,:));
r4=correlate(goldseq(1,:),goldseq(2,:));
N=length(mseq1);
subplot(2,1,1);%风别画出相关图行
plot(-N+1:N-1,r1,-N+1:N-1,r2,'r--')
title('m序列');legend('—m序列自相关','--m序列互相关',1);
xlabel('k');ylabel('相关值');
subplot(2,1,2);
plot(-N+1:N-1,r3,-N+1:N-1,r4,'r--')
title('gold序列');legend('—gold序列自相关','--gold序列互相关',1);
xlabel('k');ylabel('相关值');