c_a.m

来自「该程序是GPS信号C/A码的仿真程序」· M 代码 · 共 47 行

	k1=2;k2=6;delay=5;
    reg=-ones(1,10);%定义寄存器一到十的初始值都为-1
    for j=1:1023;
    MLS(j)=reg(10);%将10号寄存器的输出作为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);
    end
    %将G2与MLS进行延迟检验
    if MLS==g2([delay:1023 1:delay-1])
        disp('OK')
    else
        disp('not matach');
    end
    %在G2序列中找出-1并转换成0,找出1并转换成1
    ind1=find(g2==-1);
    ind2=find(g2==1);
    g2(ind1)=zeros(1,length(ind1));
    g2(ind2)=ones(1,length(ind2));
    temp=g2(1:120);
    x(1)=0;show(1)=temp(1);
    
    
    p=2;
    %将结果显示成方波形式
    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
    %画出仿真波形
    plot(x,show);
    axis([0 length(x)-60 -0.1 1.1]);