4份GPS接收机捕获与跟踪、解调的MATLAB程序

4份GPS接收机捕获与跟踪、解调的MATLAB程序下载，需要的下载

samplesPerCode = round（settings.samplingFreq / ...
                        （settings.codeFreqBasis / settings.codeLength））;

% Create two 1msec vectors of data to correlate with and one with zero DC
signal1 = longSignal（1 : samplesPerCode）;
signal2 = longSignal（samplesPerCode+1 : 2*samplesPerCode）;
% index=ceil（samplesPerCode*（1 : 40920）/40920）;
signal0DC = longSignal - mean（longSignal）; 

% Find sampling period
ts = 1 / settings.samplingFreq;

% Find phase points of the local carrier wave 
phasePoints = （0 : （samplesPerCode-1）） * 2 * pi * ts;

% Number of the frequency bins for the given acquisition band （500Hz steps）
numberOfFrqBins = round（settings.acqSearchBand * 2） + 1;

% Generate all C/A codes and sample them according to the sampling freq.



%--- Initialize arrays to speed up the code -------------------------------
% Search results of all frequency bins and code shifts （for one satellite）
results     = zeros（numberOfFrqBins 2046）;

% Carrier frequencies of the frequency bins
frqBins     = zeros（1 numberOfFrqBins）;


%--- Initialize acqResults ------------------------------------------------
% Carrier frequencies of detected signals
acqResults.carrFreq     = zeros（1 32）;
% C/A code phases of detected signals
acqResults.codePhase    = zeros（1 32）;
% Correlation peak ratios of the detected signals
acqResults.peakMetric   = zeros（1 32）;
IQ1_add=zeros（11023*2）;
IQ2_add=zeros（11023*2）;
fprintf（‘（‘）;

% Perform search for all listed PRN numbers ...
for PRN = settings.acqSatelliteList

%% Correlate signals ======================================================   
    %--- Perform DFT of C/A code ------------------------------------------
    caCodes=generateCAcode（PRN）;
    caCodes1=kron（caCodesones（12））;
%     caCodes2=[caCodes1 0 0];
    caCodeFreqDom = conj（fft（caCodes1））;
    %--- Make the correlation for whole frequency band （for all freq. bins）
    for frqBinIndex = 1:numberOfFrqBins

        %--- Generate carrier wave frequency grid （0.5kHz step） -----------
        frqBins（frqBinIndex） = settings.IF - ...
                               （settings.acqSearchBand/2） * 1000 + ...
                               0.5e3 * （frqBinIndex - 1）;

        %--- Generate local sine and cosine -------------------------------
        sinCarr = sin（frqBins（frqBinIndex） * phasePoints）;
        cosCarr = cos（frqBins（frqBinIndex） * phasePoints）;

        %--- “Remove carrier“ from the signal -----------------------------
        I1      = sinCarr .* signal1;
        Q1      = cosCarr .* signal1;
        I2      = sinCarr .* signal2;
        Q2      = cosCarr .* signal2;
        IQ1=I1 + j*Q1;
        IQ2=I2 + j*Q2;
%         IQ11=IQ1（:index）;
%         IQ22=IQ2（:index）;
        for i=1:1023*2
            IQ1_add（i）=sum（IQ1（（i-1）*8+1:8*i））;
            IQ2_add（i）=sum（IQ2（（i-1）*8+1:8*i））;
        end
        %--- Convert the baseband signal to frequency domain --------------

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

     文件       1301  2006-07-05 17:41  BOC\boc.m

     文件       1150  2006-06-02 16:06  BOC\codegen.m

     文件       1175  2006-07-04 16:53  BOC\cacode.m

     文件       1280  2006-07-04 16:52  BOC\boc.asv

     文件        800  2006-06-23 15:07  BOC\boc_bak.m

     文件       1432  2006-07-04 16:53  BOC\gcode.m

     目录          0  2006-07-20 23:03  BOC

     文件       1432  2006-08-03 22:23  BOC\BOC\boc.asv

     文件       1468  2006-08-20 00:05  BOC\BOC\boc.m

     文件        800  2006-06-23 15:07  BOC\BOC\boc_bak.m

     文件       1175  2006-07-04 16:53  BOC\BOC\cacode.m

     文件       1150  2006-06-02 16:06  BOC\BOC\codegen.m

     文件       1432  2006-07-04 16:53  BOC\BOC\gcode.m

     文件         90  2006-08-20 00:24  BOC\BOC\Untitled2.m

     目录          0  2006-08-20 19:30  BOC\BOC

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

                14903                    16