主要是关于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