matlab_扩频通信中的gold序列特性仿真研究

clear all; byxishu1=[1 1 0 1 1]; byxishu2=[0 1 0 0 1]; %本原多项式系数 mseq1=msequence(byxishu1); mseq2=msequence(byxishu2); balancegoldseq=balancegoldsequence(byxishu1,byxishu2); %******下面计算第一个平衡gold序列与其他由m序列产生的gold序列的互相关值****** ind=find(balancegoldseq ==0); %返回平衡Gold序列中的0元素的位置 balancegoldseq(ind)=-1; %把平衡Gold序列中的0元素修改成－1，使该序列值为双极性{-1,1} r=xiangguan(balancegoldseq (1,:),balancegoldseq (2,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (3,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (4,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (5,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (6,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (7,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (8,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (9,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (10,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (11,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (12,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (13,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (14,:)); %r=xiangguan(balancegoldseq (1,:),balancegoldseq (15,:)); %下面是求第一个平衡gold序列与第1个m序列的互相关值的部分程序 %ind1=find(mseq1==0); %mseq1(ind1)=-1; %r=xiangguan(balancegoldseq (1,:),mseq1(1,:)); %下面是求第一个平衡gold序列与第2个m序列的互相关值的部分程序 %ind2=find(mseq2==0); %mseq2(ind2)=-1; %r=xiangguan(balancegoldseq (1,:),mseq2(1,:)); %******************************************************************** N=length(mseq1);axis=-N+1:N-1; %定义横轴使自相关峰值位于图形的零点 plot(axis,r);xlabel('k');ylabel('R(k)'); grid on;