GPS信号捕获仿真

产生GPS的C/A码，导航电文，实现捕获仿真，有图，有真相哦

function  CACode = Capture_CACode（IFSignalCACodefreq_doppler）

N = length（IFSignal）;
F0=2e6;                  %中频频率为2MHz
Fs = 5e6;                %采样率
FSN = Fs/1e3;
NA = 5;                  %累积次数
%   生成本地C码
phase_shift=0;           %本地码码偏
index=（1: FSN*NA）/Fs;
index_code=floor（mod（index*1.023e6+phase_shift1023））+1;
index_code（find（index_code == 0）） = 1023;

local_code = CACode;
local_code=2*local_code-1;
local_code=local_code（index_code）;

%   捕获
for findex=1:40
    f = （findex - 20）*500;
    signal_f = IFSignal（1:FSN*NA） .* exp（j*2*pi*（F0-f）*（1:FSN*NA）/5e6）;
    sat_mat = reshape（signal_f（1:FSN*NA） FSN NA）;
    lcode_mat = reshape（local_code（1:FSN*NA） FSN NA）;
    fs = fft（sat_mat）;
    fl = fft（lcode_mat）;
    mu = conj（fs） .* fl;
    cr = （ifft（mu））‘;
    i = 1:NA;
    accu（findex :） = sum（abs（cr（i :）））;
end
%捕获结果
[aabb]=max（accu）;
[ccdd]=max（aa）;
a=cc;
b=bb（dd）;
c=dd;
phase_shift = c - 1            %码相位
freq = -1 * （b - 20） *500      %频偏
CACode = accu;
%figure（1）
subplot（211）;mesh（CACode）;
subplot（212）;plot（CACode（b-1:））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1127  2008-06-30 21:24  程序\Capture_CACode.m

     文件        734  2008-06-30 12:33  程序\Create_CA_Code.m

     文件        732  2008-06-30 14:22  程序\Create_NavigationData.m

     文件        721  2008-06-30 20:57  程序\IF_basebandSignal.m

     文件       1108  2008-06-30 21:23  程序\main.m

     文件        694  2008-06-30 19:47  程序\SpreadSpectrum.m

     目录          0  2009-06-16 23:24  程序

----------- ---------  ---------- -----  ----

                 5116                    7