大学毕业设计一席谈之三十四 GPS卫星采样信号的跟踪(5)调节码相位实现跟踪

rxDemod = rxDemod(1:16368000*5)';  % 5s长度

% 4倍抽取

rxDemod = rxDemod(1:4:end);

fs = 4.092e6;

ts = 1/fs;

%%% 环路参数

% Dll环中环路滤波器设计

Bd = 1;                              % 环路噪声带宽

damp = 0.707;                        % 阻尼系数

td = 1/1000;                         % 相关积分时间

Kd = 1;                              % 环路增益

Wd = 2*Bd/(damp+1/(4*damp));         % 自然角频率 

Cd1 = 8*damp*Wd*td/(Kd*(4+4*damp*Wd*td+(Wd*td)^2));  

% 环路滤波器的系数

Cd2 = 4*(Wd*td)^2/(Kd*(4+4*damp*Wd*td+(Wd*td)^2));

% Costas环中环路滤波器设计

Bp = 100;

tp = 1/1000;

Kp = 1;

Wp = 2*Bp/(damp+1/(4*damp));

Cp1 = 8*damp*Wp*tp/(Kp*(4+4*damp*Wp*tp+(Wp*tp)^2));

Cp2 = 4*(Wp*tp)^2/(Kp*(4+4*damp*Wp*tp+(Wp*tp)^2));

%%%%%%%%%%    begin tracking

% 将捕获结果分别送入各个通道！

%%%%%% 第一通道

% 相关最大值（22号卫星）进入第一通道

% 捕获结果显示为第22颗星

CAcode_1 = generateCAcode(22);  

doppler_1 = -2200; % doppler(22)*1000; % 单位 Hz