大学毕业设计一席谈之三十四 GPS卫星采样信号的跟踪(3)降采样跟踪

N = 16.368e6;

% 读取数据

fid = fopen(file_path, 'r');

rxDemod = fread(fid, N*15, 'int8');  % 读取N*5个int8数据，保留符号

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

% 4倍抽取

rxDemod = rxDemod(1:4:end); % 四选一方式！！！

fs = 4.092e6;

ts = 1/fs;

%%% 锁相环参数

damp = 0.707;

Bp = 100;     % 噪声带宽

tp = 1/1000;  % 时间间隔 1ms调整一次

Kp = 1;

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

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

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

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

%%%%%% 第一通道