GPS软件接收机matlab程序学习笔记（一）

clear all;clc;cla;close all;
format longg; %%15位有效数
%% 一、  设定GPS数据文件参数 %%%%%%%%%%%%%%
gps_bd = 0; % 0:gps, 1:bd
n_mean_count = 200; %bias count 
if gps_bd == 0
    % 这个中频数据是 bit8
    filename = '.\data\GPS_sim_5mfs_if1.42m.txt';% 相对路径 当前工程文件夹下的data文件夹下的GPS_sim_5mfs_if1.42m.txt文件
    fs = 5e6; % 采样率★书本260页实际采样在5E6
    fc = 1.42e6;% 采样数据中频★
    sv = 1:32;% 要捕获跟踪的卫星号 若不知道有哪些星可以设置为 sv=1:32 但这样处理速度会比较慢★
    Readtype = 'bit8';%数据读取位宽，要注意 不同的中频信号采样位是否一致，有符号采样时，高位为符号位，低位为幅值位

%     filename = 'D:\9.WorkInTheLab\0.Project\0.GNSS_SignalQualityMonitor\1.data\2771\data\GPSL1_8bit_fs16.369M_fc3.996M.bin';
%     fc = 3.996e6;
%     fs = 16.369e6;
%     sv = 28;
%     Readtype = 'bit8';

%     filename = 'C:\MATLAB7\work\BDSB1-8bit-PwrSpectAbnorm-202006172224.bin';
%     fc = 3.996e6;
%     fs = 16.369e6;
%     sv = 1:32;
%     Readtype = 'bit8';
    
    disp 'gps acq'
else
    filename = 'C:\MATLAB7\work\usbdata.bin';
    fc = 3.999875e6;
    fs = 16.369e6;
    sv = 1:32;
    Readtype = 'bit8';
    disp 'bd acq'
end
%% 二、  任务配置 %%%%%%%%%%%%%%
task_AcqOnly= 0  ;             %捕获标志，为1时只进行捕获，为0时则进行捕获跟踪
task_Acq_Track=1;              %捕获跟踪标志，为1 － 则进行捕获、跟踪；为 0 则调用以前结果进行NAV定位
track_acqresultsName='capture_test.mat';%%捕获跟踪后生成的数据保存名字,task_Acq_Track = 0 时则直接调用此数据进行定位解算。注意每次用不同的中频数据时取不同的名字，以免覆盖。
task_frameFindChk=  1    ;      %找子帧标志， 1则找子帧，其他则不进行                            0
task_nav = 1           ;        %定位解算标志 1则参与定位解算，其他则不进行                      0
plotAcq_AND_Channel=  1  ;      %捕获通道画图标志，1则画图，其他则不进行
plot_Acq_ACF= 1        ;        %捕获三维图和相关累加画图标志，1则画图，其他则不进行
if task_Acq_Track==1            %捕获跟踪标志，为1 － 则进行捕获、跟踪；为 0 则调用以前结果进行NAV
    task_load_results= 0   ;
else
    task_load_results=  0 ; 
e