fgeneratecacode3.m

gps c/a码程序以及模拟卫星定位的源代码、地球模型

function y=fGenerateCAcode3(svnum)
g2s=[5;6;7;8;17;18;139;140;141;251;252;254;255;256;257;258;469;470;471;...
    472;473;474;509;512;513;514;515;516;859;860;861;862];
g2shift=g2s(svnum,1);
reg=-1*ones(1,10);
for i=1:100
    g1(i)=reg(10);
    slave1=reg(3)*reg(10);
    reg(1,2:10)=reg(1:1:9);
    reg(1)=slave1;
end
reg=-1*ones(1,10);
Table=[2 6;
       3 7;
       4 8;
       5 9;
       1 9;
       2 10;
       1 8;
       2 9;
       3 10;
       2 3;
       3 4;
       5 6;
       6 7;
       7 8;
       8 9;
       9 10;
       1 4;
       2 5;
       3 6;
       4 7;
       5 8;
       6 9;
       1 3;
       4 6;];
   for i=1:100
       save2=reg(2)*reg(3)*reg(6)*reg(8)*reg(9)*reg(10);
       g2(i)=reg(Table(svnum,1))*reg(Table(svnum,2));
       reg(1,2:10)=reg(1:1:9);
       reg(1)=save2;
   end
ss_ca=g1.*g2;
ca=ss_ca;
ind1=find(ca==-1);
ind2=find(ca==1);
ca(ind1)=ones(1,length(ind1));
ca(ind2)=zeros(1,length(ind2));
y=ca;