fgeneratecacode3.m

对C/A码的相关在频域中的matlab实现仿真

function y=fGenerateCAcode3(Svnum)
%产生C/A码的方法二
% Svnum=input('Enter the satellite number =');
g2s=[5;6;7;8;17;18;139;140;141;251;252;254;255;256;257;258;469;470;471;...
    472;473;509;512;513;514;515;516;859;860;861;862];
g2shift=g2s(Svnum,1);
len=1023;%1023;
reg=-1*ones(1,10);
%产生g1
for i=1:len%1023
    g1(i)=reg(10);
    slavel=reg(3)*reg(10);
    reg(1,2:10)=reg(1,1:9);
    reg(1)=slavel;
end
%产生g2
for i=1:len%1023
    g2(i)=reg(10);
    save2=reg(2)*reg(3)*reg(6)*reg(8)*reg(9)*reg(10);
    reg(1,2:10)=reg(1,1:9);
    reg(1)=save2;
end
g2temp(1,1:g2shift)=g2(1,len-g2shift+1:len);
g2temp(1,g2shift+1:len)=g2(1,1:len-g2shift);
g2=g2temp;
%g1和g2卷积后得到C/A码
ca=g1.*g2;
ind1=find(ca==-1);
ind2=find(ca==1);
ca(ind1)=ones(1,length(ind1));
ca(ind2)=zeros(1,length(ind2));
y=ca;