modulate.m

主要是关于GPS卫星C/A码信号的产生、捕获以及跟踪的S曲线图

%这是调制信号的程序
%m是伪码的相位偏移(1023以内),doppler是传输中产生的多普勒偏移(在-5000Hz到5000Hz)
function test = modulate(m,doppler)
%产生调制信号的载波
fc =10.23e+6;
fs = 40.92e+6;
t=0:1/fs:(0.001-1/fs);
T = sin((fc+doppler)*2*pi*t);
%产生扩频码
n=1023;
[code,cacode]=generate(n);

cacode_A = cacode(1:m);
cacode_B = cacode(m+1:n);
cacode = [cacode_B,cacode_A];
for j =1:1:n
     ca_aftsample((j-1)*40+1:j*40) = repmat(cacode(j),1,40);
end;
%调制信号并加入白噪声
test=(-1)*ca_aftsample.* T;
%test = awgn(test, -12, 'measured'); %加高斯白噪声，信噪比是-12dB