大学毕业设计一席谈之三十四 GPS卫星采样信号的跟踪(2) 程序更新

% 捕获的卫星号:     6    11    19    22

% 卫星信号的频偏值:单位1KHz   1.1000    2.7000    0.9000   -2.2000

% 卫星信号的码相位值:      3020         410        2315         477

% 读取数据

fid = fopen(file_path, 'r');

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

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

fs = 4.092e6*4;

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

%%%%%% 第一通道

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

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

CAcode_1 = generateCAcode(22);  

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

posi_1 = 477;      % posi(22);

% 根据PN序列位置调整到1ms的起始位置