clc;
clear all;
%產(chǎn)生C/A碼的方法一
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%,%%%%%%%
    k1=2;k2=6;delay=5;    %定義參數(shù)k1、k2以及延遲
    Reg=-ones(1,10);      %定義寄存器1～10的初始值都為-1
                        %通過循環(huán)產(chǎn)生長度為1024的MLS序列以及G2序列
    for j=1:1023
      MLS(j)=Reg(10);    %將士0號寄存器的輸出作為MLS輸出
    modulo=Reg(2)*Reg(3)*Reg(6)*Reg(8)*Reg(9)*Reg(10);
    Reg(2:10)=Reg(1:9);
    Reg(1)=modulo;
    g2(j)=Reg(k1)*Reg(k2);    %將參數(shù)k1和k2所代表的寄存器模2相加后作為G2輸出
                              %g2即為輸出
    end
    %將G2與MLS進(jìn)行延遲檢驗(yàn)
     if MLS == g2([delay:1023 1:delay-1])
     disp('OK') ;
    else  
    disp('not match')
    end
    %在G2序列中找出-1并轉(zhuǎn)換為0,找出1并轉(zhuǎn)換為1
    ind1=find(g2==1);
    ind2=find(g2==1);
    g2(ind1)=ones(1,length(ind1));
    g2(ind2)=zeros(1,length(ind2));
    temp=g2(1:120);
x(1)=0;
Show(1)=temp(1);
    P=2;
    %下面的循環(huán)是為了將結(jié)果顯示成方波形式
    for i=2:length(temp)
       if temp(i)==temp(i-1)
       x(P)=i-1;
   Show(P)=temp(i-1);
       x(P+1)=i-1+0.01;
       Show(P+1)=temp(i);
       P=P+2;
else
    Show(P)=temp(i);
    x(P)=i;
    P=P+1;
    end
    end
    %畫出仿真結(jié)果圖
    plot(x,Show+1);
    axis([0  length(x)-60 -0.1 1.1]);
    grid on;