boccode.m

gps/galileo信号源频谱分析

function BOC_Code = BOCcode(phcode,fd,Loop_time)

global time_unit;
global fs;
global nn;
global F_if;
global CA_freq;
global tc;
global CA;
global F_Carrier;
global CA_length;

ts = (1+(fd/F_Carrier))/fs;

if(phcode < 0)
    offset = (1023+phcode)/(CA_freq);   
else
    offset = phcode/(CA_freq);
end

CA_S=[];
for i = 0:nn-1
    time_dyn = mod((ts*(i+nn*(Loop_time-1))+offset),CA_length); 
    n = ceil(time_dyn/tc)+1;
    if (n > length(CA))
        n = n-1;
    else
        n = n;
    end
    CA_S(i+1) = CA(n);
end
CA_S;
BOC_Code = CA_S;